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, ReaxFFDPD, 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的输出文件是一种简单的文本格式,我们也可以通过其它的工具来换专这些格式。
.我们可以用以下几个软件来完成高质量的可视工作:VMDAtomEyePymolRaster3dRasMol的,可能,它们大多数是并行最后要说一下的是,以下这些也是自由分子动力学包LAMMPS联合起来使用以完成模拟工作。
也适合来完成你的研究工作,当然也可以与CHARMMAMBERNAMDNWCHEMDL_POLYTinkerTinker是专们用于模拟生物分子的。
,NWCHEM,CHARMM,AMBER,NAMD二、开始本部分主要描述如何创建和运行LAMMPS。
1.在LAMMPS发行包理含有:READMELICENSEBench:测式任务Doc:文本Examples:简单的测试任务Potentials:嵌入原子方法与力场文件Src:源代码Tools: 前处理与后处理工具假如你下载的是windows可执行文件的话,你里面只有一个文件(并行与非并行两种)Lmp-windows.exe2. 编译LAMMPS之前的工作:编译LAMMPS不是一个繁琐的工作。
首先你可能要写一个makefile文件,里面要选择编译器,附加的一些将要用到的库等。
事先装上MPI或FFT等库。
编译出一个可执行LAMMPS:在SRC目录里头含有C++源文件和头文件。
当然也包括一个高水平的Makefile,在MAKE目录里头有几个低水平的Makefile.*files分别适有不同的平台。
进入SRC目录,输入make或gmake,你将会看到一列的可选项。
假如其中有一种符合你的机器,你可以输入像下面一样的命令:Make linuxGmake mac注意,在一个多处理器或多核处理器的平台上你可以进行平行编译,在make命令中使用“-j”选项就可以,这样编译起来会更快一些。
.lmp-linux的可执行文件。
在此过程中不发生错误的话,你可以得到一个类似于在编译过程当中将会发生的常见错误:说明你*为名的文件的话,)如果编译过程当中发生错误,并提示不能找到一个含有通配符(1。
如还不行的话,中使用通配符。
那就偿式使用gmakemake器允许makefile 机器上的make对像。
如:选项,用Makefile.list作为就试试加入-fMake makelistf makefile.list linux Make –f makefile.lst macGmake –时,可能由于对机器的设置不正确,会导致一些错误。
假如你的当你使用低水平的makefile(2)。
使用任何一个与你机器相Makefile.foo的话,你将要在MAKE目录中创建一个,平台叫“foo”近的文件作为开始总是一个不错的选择。
)如你在链接的时候出现库丢失或少了依赖关系的话,可能是由于:(3package libaray. 你编译的包需要一个附加的库,但却没有事先编译需要的你要链接的库在你的系统中不存在。
没有连接到必要的系统库makefile.foo. 后两种问题出现,你就需要修改你的低水平makefile.foo:编辑一个新的低水平,不论你写成什么,这一行将会出现在屏幕上,如果你只输入foo在#后的句子中,替换(1)命令的话。
make包括优编译器列出编译器与链接器的设置,部分为你的C++在“complier/linker settings”(2)如果你,当然你也可以用MPICC你可以在任何UNIX系统中使用G++编译器。