LAMMPS手册-中文版讲解

LAMMPS手册-中文解析一、简介

本部分大至介绍了LAMMPS的一些功能和缺陷。

1．什么是LAMMPS?

LAMMPS是一个经典的分子动力学代码，他可以模拟液体中的粒子，固体和汽体的系综。他可以采用不同的力场和边界条件来模拟全原子，聚合物，生物，金属，粒状和粗料化体系。LAMMPS可以计算的体系小至几个粒子，大到上百万甚至是上亿个粒子。

LAMMPS可以在单个处理器的台式机和笔记本本上运行且有较高的计算效率，但是它是专门为并行计算机设计的。他可以在任何一个按装了C++编译器和MPI 的平台上运算，这其中当然包括分布式和共享式并行机和Beowulf型的集群机。LAMMPS是一可以修改和扩展的计算程序，比如，可以加上一些新的力场，原子模型，边界条件和诊断功能等。

通常意义上来讲，LAMMPS是根据不同的边界条件和初始条件对通过短程和长程力相互作用的分子，原子和宏观粒子集合对它们的牛顿运动方程进行积分。高效率计算的LAMMPS通过采用相邻清单来跟踪他们邻近的粒子。这些清单是根据粒子间的短程互拆力的大小进行优化过的，目的是防止局部粒子密度过高。在并行机上，LAMMPS采用的是空间分解技术来分配模拟的区域，把整个模拟空间分成较小的三维小空间，其中每一个小空间可以分配在一个处理器上。各个处理器之间相互通信并且存储每一个小空间边界上的”ghost”原子的信息。LAMMPS(并行情况)在模拟3维矩行盒子并且具有近均一密度的体系时效率最高。

2．LAMMPS的功能

总体功能：

可以串行和并行计算

分布式MPI策略

模拟空间的分解并行机制

开源

高移植性C++语言编写

MPI和单处理器串行FFT的可选性（自定义）

可以方便的为之扩展上新特征和功能

只需一个输入脚本就可运行

有定义和使用变量和方程完备语法规则

在运行过程中循环的控制都有严格的规则

只要一个输入脚本试就可以同时实现一个或多个模拟任务

粒子和模拟的类型：

（atom style命令）原子粗粒化粒子DNA 全原子聚合物，有机分子，蛋白质，联合原子聚合物或有机分子金属粒子材料粗粒化介观模型延伸球形与椭圆形粒子点偶极粒子刚性粒子所有上面的杂化类型力场：）（命令：pair style, bond style, angle style, dihedral style, improper style, kspace style, tabulated.

class2(COMPASS)L-J, Buckingham, Morse, Yukawa, soft, 对相互作用势：Coulombic, point-dipole.

带点对相互作用势：Stillinger-Weber, modified EAM(MEAM), EAM, 多体作用势：Finnis/Sinclair EAM,

Tersoff, AIREBO, ReaxFF

DPD, GayBerne, Resquared, Colloidal, DLVO 粗粒化作用势：granular, Peridynamics

介观作用势：harmonic, FENE, Morse, nonlinear, class2, quartic

键势能：harmonic, CHARMM, cosine, cosine/squared, class2(COMPASS) 键角势能：harmonic, CHARMM, multi-harmonic, helix, OPLS, class2(COMPASS) 二面角势能：harmonic, CVFF, class2(COMPASS) 不合理势能：all-atom, united-atom,

bead-spring, breakable 聚合物势能：水势能：TIP3P，TIP4P，SPC

隐式溶剂势能：hydrodynamic lubrication, Debye

长程库伦与分散：Ewald, PPPM, Ewald/N(针对长程L-J作用)

可以有与普适化力场如CHARMM，AMBER，OPLS，GROMACS相兼容的力场可以采用GPU加速的成对类型

杂化势能函数：multiple pair, bond, angle, dihedral, improper potentials(多对势能处于更高的优先级)

原子创建：

（命令：read_data, lattice, create-atoms, delete-atoms, displace-atoms, replicate）

从文件中读入各个原子的坐标

在一个或多个晶格中创建原子．

删除几何或逻辑原子基团

复制已存在的原子多次

替换原子

系综，约束条件，边界条件：

（命令：fix）

二维和三维体系

正角或非正角模拟空间

常NVE，NVT，NPT，NPH积分器

原子基团与几何区域可选择不同的温度控制器

有Nose/Hoover和Berendsen压力控制器来控制体系的压力（任一维度上）

模拟合子的变形（扭曲与剪切）

简谐（unbrella）束缚力

刚体约束

摇摆键与键角约束

各种边界环境

非平行太分子动力学NEMD

各种附加边界条件和约束

积分器：

Velocity-verlet积分器

Brown积分器

rRESPA继承时间延化积分器

刚体积分器

共轭梯度或最束下降算法能量最小化器

输出：

（命令：dump, restart）

热力学信息日志

原子坐标，速度和其它原子量信息的文本dump文件

二进制重启文件

各原子量包括：能量，压力，中心对称参数，CAN等

用户自定义系统宽度或各原子的计算信息

每个原子的时间与空间平均

系统宽量的时间平均

原子图像，XYZ，XTC，DCD，CFG格式

数据的前处理与后处理：

包里提供了一系列的前处理与后处理工具

模拟的设置，分LAMMPS它可以进行pizza.py, 另外，可以使用独立发行的工具组．

析，作图和可视化工作。特别功能：MD模拟实时的可视化与交互式-连续体模拟与有限元方法结合进行原子POEMS库中提供了刚体积分工具在并行裉火并行复制动力学模拟MC对低密度液体直接使用介观建模Peridynamic 目

标型与无目标型分子动力学双温度电子模型 LAMMPS不具备的功能：是对牛顿

运动方程积分的工具，所以很多必要的数据前处理与后处由于LAMMPS

：理功能是LAMMPS核心不具备的。其原因为的小巧性保证LAMMPS 后处理不能进行并行运算前处理与这些功能可以有其它工具来完成原代码开发的局限性

不能：特别地，LAMMPS通过图形用户界面来工作创建分子体系自动的加上力场系数MD模拟提供智能化的数据分析为MD的可视化为输出数据作图

输入一系列的原子类型，原子坐标，分子拓朴信息和所有原LAMMPS我们需要为。子与键的力场参数。LAMMPS不会自动的为我们创建分子体系与力场参数命令来为固态晶格加上原子。可以能creat-atomsLAMMPS提供了对与原子体系，等命来加上小数目的力场参数。对于分子体系或更复pair coeff,bond coeff, angle coeff过输出文件来做到这杂的模拟体系，我们通常会用其它工具来创建或者是转换LAMMPS 有的还会写一些自已的代码来完成这项任务。些事情。我们需要为之提供上面个拓朴信息与力场，对于一个复杂的分子体系（如，蛋白质）或其它的分子建模器来完成这些任务，并或CHARMMAMBER参数。所以我们建议用所允许的输入格式。把之输到一个文件中去。然后，改变其格式以达到LAMMPS同样，LAMMPS的输出文件是一种简单的文本格式，我们也可以通过其它的工具来换专这些格式。．

我们可以用以下几个软件来完成高质量的可视工作：VMD

AtomEye

Pymol

Raster3d

RasMol

的，可能，它们大多数是并行最后要说一下的是，以下这些也是自由分子动力学包LAMMPS联合起来使用以完成模拟工作。也适合来完成你的研究工作，当然