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联合起来使用以完成模拟工作。

也适合来完成你的研究工作，当然
也可以与CHARMM
AMBER
NAMD
NWCHEM
DL_POLY
Tinker
Tinker是专们用于模拟生物分子的。

，NWCHEM，CHARMM，AMBER，NAMD
二、开始
本部分主要描述如何创建和运行LAMMPS。

1．在LAMMPS发行包理含有：
README
LICENSE
Bench：测式任务
Doc：文本
Examples：简单的测试任务
Potentials：嵌入原子方法与力场文件
Src：源代码
Tools: 前处理与后处理工具
假如你下载的是windows可执行文件的话，你里面只有一个文件（并行与非并行两种）
Lmp-windows.exe
2. 编译LAMMPS
之前的工作：
编译LAMMPS不是一个繁琐的工作。

首先你可能要写一个makefile文件，里面要选择编译器，附加的一些将要用到的库等。

事先装上MPI或FFT等库。

编译出一个可执行LAMMPS：
在SRC目录里头含有C++源文件和头文件。

当然也包括一个高水平的Makefile，在MAKE目录里头有几个低水平的Makefile.*files分别适有不同的平台。

进入SRC目录，输入make或gmake，你将会看到一列的可选项。

假如其中有一种符合你的机器，你可以输入像下面一样的命令：Make linux
Gmake mac
注意，在一个多处理器或多核处理器的平台上你可以进行平行编译，在make命令中使用“-j”选项就可以，这样编译起来会更快一些。

．
lmp-linux的可执行文件。

在此过程中不发生错误的话，你可以得到一个类似于在编译过程当中将会发生的常见错误：说明你*为名的文件的话，）如果编译过程当中发生错误，并提示不能找到一个含有通配符（1。

如还不行的话，中使用通配符。

那就偿式使用gmakemake器允许makefile 机器上的make对像。

如：选项，用Makefile.list作为就试试加入-fMake makelist
f makefile.list linux Make –f makefile.lst mac
Gmake –时，可能由于对机器的设置不正确，会导致一些错误。

假如你的当你使用低水平的makefile(2)。

使用任何一个与你机器相Makefile.foo的话，你将要在MAKE目录中创建一个,平台叫“foo”近的文件作为开始总是一个不错的选择。

）如你在链接的时候出现库丢失或少了依赖关系的话，可能是由于：（3package libaray. 你编译的包需要一个附加的库，但却没有事先编译需要的你要链接的库在你的系统中不存在。

没有连接到必要的系统库makefile.foo. 后两种问题出现，你就需要修改你的低水平makefile.foo:
编辑一个新的低水平，不论你写成什么，这一行将会出现在屏幕上，如果你只输入foo在#后的
句子中，替换（1）
命令的话。

make包括优编译器列出编译器与链接器的设置，部分为你的C++在“complier/linker settings”（2）如果你，当然你也可以用MPICC你可以在任何UNIX系统中使用G++编译器。

化符号。

的话。

的系统中安装了MPI 变量中列出来。

如过在编译过程当中需要符加的库的话，你必需在LIB编译器创建一个源文件的依赖关系列表，当源文件或头文件改变的设置可以让
C++DEPFLAGS来实现。

时候可以加快编译速度。

有些编译器不能创建依赖关系列表，或者你
可以用选项—D。

如果你的编译器不能创建依赖关系文件的话，那么你就需要创建一个可以使用-DG++ Makefile.storm（它用一系列的不需要依赖文件的规划）相对应。

Makefile.foo来与system-specific settings”部分有四个小部分：3）“（INC变量，包括一些与系统相关的条件选项。

A LMP—库。

MPI个MPI变量用于指定
B 3如你要进行并行计算的话，那么你必须在你的平台上安上MPI库。

如你想用MPI内置C++编译器的话，你可以让这三个变量空着，如你不用MPIC
C 的话，那么，你要指定MPI.h(MPI_INC)文件在哪，MPI(MPI_PA TH)库在哪，还有库名(MPI_LIB). 如果你想自已安装MPI的话，我们建议用MPICH1.2或2.0。

LAM MPI也可以。

如果我的是大
平始的话，你的供应商已经为你装上了MPI，其可能比MPICH或LAM更快，你可以把找出来
并与之链接。

如你用LAM或MPICH，你必需要设置他并编译他使之适合你的平台。

如果你想在单处理器的机器上运行的话，你可以用STIBS库，这样你就可以不用在你的系统中
安装MPI库。

防照makefile.serial，看是如果设置这三个变量的。

当然你在编译LAMMPS之前
你必需创建STUBS库。

在STUBS目录中，输入make，不出错的话你将会得到一个libmpi.a文件可供链接到LAMMPS。

当出错，你则要修改STUBS下的MAKEFILE。

STUBS/MPI.CPP有一个CPU计时器MPI_Wtime()可以调用gettimeofday(). 如你的系统不支持gettimeofday()，则你就要插入一句代码来调用另一个计时器，要注意的是，clock()函数在一个小时之后会归0，所以对于一个长时间的LAMMPS模拟来说这是不够用的。

选项PPPM命令来计算长程库伦作用时使用kspace-style当要用到库，FFT变量用于指定C FTT．时要用到。

XXX—FFT库。

可以能过开关—DFFT要使用此选项，你必须要在你的机器上安装
一个一维的。

没有办法的情况下可以用供应商提SCSL，或FFTWDEC来指定，其中XXX=INTEL，，SGI，本。

FFTW是一个快速的，可移植性的库，它可以在任何一个平台上运行。

最好2.1.X供
的库。

makefile.fooFFT库，你都要在编译FFTW库时只要用./configure;make就可以。

不任你是
用哪一种中正确的设置咱们的FFT_INC,FFT_PATH,FFT_LIB。

设FFT库。

这种情况下，你可以把FFT_INC当然，你如果不用PPPM的话，你将没有必要安装KSPACE包剔除。

成-DFFT-NONE 并让其它几个变量空着。

你也可以在编译LAMMPS时把时其中有一个或几个包除非你在编译LAMMPSSYSLIB和SYSPATH变量你可以忽略，D 几个变量的前辍。

和SYSPA TH要用到附加的系统库。

所有这些包都的的名称都将会是SYSLIB SYSPA TH则是路径，只有当这些库为非默认路径时才有设定。

SYSLIB变量将列出系统库。

之后，你只)和预编译好了所有的其它库最后，当你正确的写好了makefile.foo(MPI,FFT,包库等目录下输入下面其中一个命令就可以了要在SRCMake foo
Gmake foo
的可执行文件。

不出意外，你将会得到lmp_foo 附加建义LAMMPS
为多平台编译（1）。

每一个目标都有他自已的目标路径，目录下为多平台编译LAMMPS你可
以在同一个SRC Obj_name 用于存贮指定系统的目标文件。

清理2）（在编译时创建的目标
文件。

make clean-foo将会清理LAMMPS输入make clear-all或LAMMPS
编译为windows）（3
可执行文件。

下载页面上可以下载已经编译好了的windows在LAMMPS平台上从windowswindows可执行文件不能满足你的要求的时候，你可以在如已经编译好了
的./src/MAKE/Windows. LAMMPS。

但是不建义这么做。

见源文件编译LAMMPS 3 如何有选择性的编译/剔除包加入LAMMPS编译之前可以加入或剔除所有或部分包。

在如果不是这样，这此包中相互关联的包必须包括在一起。

和OPT包。

只有两种况下是列外，GPU包须asphereGPU 包里的全部文件，GPU那么在和OPT子目录下的所有包都不能安装。

要安装manybody包须首先安装。

OPT下的所有文件，那么kspace与安装。

要安装当然，由于某些模拟用不到其中的某些包，那么你想把这些包剔除这是可以的。

这样你就可以不用编译一些额外的库，当然你的可执行文件也将会更小，运行起来也会更快。

默认情况下，LAMMPS只包含
kspace,manybody,molecule这三个包。

可以通过如make yes-name, make no-name来加入和剔除一些包，当然你也可以用make
yes-standard, make no-standard,make yes –user, make no-user, make yes-all, make no-all加入或剔除可以看到多个选项。

Make package各种包。

．。