第五讲分子动力学模拟的Lammps实现
LAMMPS软件与分子模拟的实现
LAMMPS软件与分子模拟的实现LAMMPS (Large-scale Atomic/Molecular Massively Parallel Simulator) 是一个基于粒子动力学原理的分子模拟软件。
它使用分子动力学模型来模拟原子、分子或其他粒子在不同温度、压力和相互作用条件下的行为。
它是一个高效、可扩展和灵活的软件，可以模拟从数百到数百万个粒子的多种物理和化学现象。
1. 引入粒子和相互作用模型: LAMMPS实现了多种粒子和相互作用模型。
用户可以指定模拟系统中的粒子类型，包括原子、分子和其他粒子类型。
LAMMPS支持多种相互作用力场模型，如Lennard-Jones和Coulomb 相互作用，以及更复杂的模型如多体相互作用。
2. 粒子动力学模拟: LAMMPS使用经典的牛顿力学原理来模拟粒子在时间和空间上的演化。
它迭代破解了每个粒子所受到的力，并计算粒子的速度和位置。
它使用了一些高效的算法和数据结构来提高模拟效率，如Verlet积分算法和空间分解技术。
3. 温度和压力控制: LAMMPS可以在模拟过程中控制系统的温度和压力。
它采用了多种算法来模拟温度和压力，如Nose-Hoover算法、Berendsen热浴、Langevin动力学和Parrinello-Rahman方法。
这些算法可以在模拟过程中维持系统的平衡状态。
4.边界条件和周期性边界条件:LAMMPS支持各种不同的边界条件。
它可以模拟有限尺寸系统，也可以模拟无限尺寸系统。
对于无限尺寸系统，LAMMPS采用了周期性边界条件，以模拟系统中的无限复制。
5.输入和输出:LAMMPS提供了灵活的输入和输出功能。
用户可以通过输入文件来设置模拟系统的参数，如初始位置、速度、力场模型和模拟时间。
LAMMPS会将模拟结果输出到文件中，用户可以对结果进行分析和后处理。
6.并行计算:LAMMPS是一个并行化的软件，可以在多个计算节点上并行计算，以提高计算效率。
全原子分子动力学模型 lammps
全原子分子动力学模型 lammps全原子分子动力学模型LAMMPS，是一款非常优秀的分子模拟软件。
它是一款免费的并依托开源社区共同开发的分子模拟软件，在学术界和工业界都具有广泛的应用。
LAMMPS包含许多强大的功能和工具，能够模拟分子、多体相互作用、材料能量和温度等方面，是材料科学、化学、生物学等领域研究的重要工具之一。
下面我们来具体了解一下如何使用LAMMPS进行分子模拟。
第一步：软件安装与配置首先，我们需要前往LAMMPS的官方网站进行下载和安装。
下载的版本可以根据自己的需要选择，一般来说最新的版本越稳定也越实用。
安装之后，我们需要配置环境变量，以便在终端或命令行中可以直接使用LAMMPS。
第二步：建立分子模型在使用LAMMPS进行分子模拟之前，我们需要首先建立分子模型。
这可以通过算法或者数据实验等方式实现。
具体来说，我们需要确定分子的数目、类型、位置等信息。
对于这些信息，可通过多种科学方法获取。
我们建立好分子模型之后，需要将其写入到LAMMPS的输入文件中。
输入文件包含了我们的模型、模拟参数、计算方式和输出等信息，是LAMMPS模拟的核心。
第三步：设置模拟参数LAMMPS除了支持模型参数输入外，还提供了一个非常强大的用户交互机制，以便更灵活地控制模型。
在这里，我们可以设置温度、压力、能量、力场、约束等不同的模拟参数。
不同的模型需要根据具体应用需求进行不同参数的调整，比如需要考虑不同的温度、压力等等。
第四步：运行模拟当我们设置好了LAMMPS的输入文件和模拟参数之后，就可以开始利用LAMMPS进行模拟了。
一般来说，我们可以采用命令行操作，以便更精确地控制模拟进程。
模拟完成之后，我们可以根据之前设置的输出选项进行相应的结果分析。
LAMMPS支持多种输出格式，方便进行分析和后续处理。
总结：通过以上步骤，我们可以看到使用LAMMPS进行分子模拟的过程非常清晰和简单。
LAMMPS强大的功能和灵活性，可以帮助我们快速、准确地获取分子的性质和行为，是当今分子模拟研究领域的重要工具之一。
份分子动力学模拟资源lammps+MS适合初学者
完成安装
配置环境变量
下载LAMMPS和MS安装包 运行安装程序
选择安装选项
测试安装是否成功 验证安装是否成功
LAMMPS和MS配置参数
• LAMMPS配置参数： - 内存需求：根据模拟系统大小和精度要求进行配置 - 输入文件：包含系统参数、初始 构型等信息的文件 - 输出文件：模拟过程中的轨迹、能量等信息的输出文件 - 命令行参数：用于控制模拟过 程的各种参数
MS：Materials Studio，一款材料科学模拟软件，用于模拟材料的物理和化学性质
LAMMPS和MS的共同点：都可以进行分子动力学模拟，但LAMMPS更侧重于模拟分子和 原子的运动，而MS更侧重于模拟材料的物理和化学性质
LAMMPS和MS的区别：LAMMPS更适用于初学者，因为它的界面更简洁，操作更简 单，而MS则更适合有一定基础的用户，因为它的功能更丰富，可以模拟更复杂的材料 性质。
LAMMPS和MS应用领域
材料科学：模拟材 料微观结构，研究 材料性能
生物物理：模拟生 物大分子，研究蛋 白质折叠、DNA 复制等
化学物理：模拟化 学反应，研究化学 反应机理
纳米科学：模拟纳 米材料，研究纳米 材料的性质和应用
LAMMPS和MS安装与配置
LAMMPS和MS安装步骤ຫໍສະໝຸດ 解压安装包输入安装路径
• MS常用命令： - create：创建分子 - add：添加分子 - delete：删除分子 - move：移动分子
• - create：创建分子 • - add：添加分子 • - delete：删除分子 • - move：移动分子
• MS常用参数： - atom_style：原子类型 - bond_style：键类型 - angle_style：角类型 - dihedral_style：二面角类型 - improper_style：非平面角类型
用Lammps实现模拟
Lammps 计算输入文件
• # create geometry 创建初始几何构形 • • • • • • • • • • • • • • Lattice hex 0.93 # 指定晶格类型（二维 hex ）和晶格常数 Region box block 0 100 0 40 -0.25 0.25 # 定义一个区域 create_box 5 box # 在指定区域建立一个 simulation box,5 表示原子类型的种类数 create_atoms 1 box # 在 simulation box 中创建类型为 1 的原子（原子位置初始化） Mass 1 1.0 Mass 2 1.0 Mass 3 1.0 Mass 4 1.0 Mass 5 1.0 # 指定单个原子的质量
A. 常用分子动力学模拟程序
• LAMMPS 免费 / 一般性分子模拟软件。 兼容当前大多数的势能模型，编程水平高，计算效率高。可以 模拟软材料和固体物理系统。 • Materials Explorer 立足于 Windows 平台的多功能分子动力学软件。拥有强大 的分子动力学计算及 Monte Carlo 软件包，是结合应用领域来 研究材料工程的有力工具。 Materials Explorer 可以用来研究 有机物、高聚物、生物大分子、金属、陶瓷材料、半导体等晶 体、非晶体、溶液，流体，液体和气体相变、膨胀、压缩系数 、 抗 张 强 度 、 缺 陷 等 。 Materials Explorer 软 件 中 包 含 2Body ， 3Body ， EAM ， AMBER 等 63 个力场可供用户选 择。 Materials Explorer 软件拥有完美的图形界面，方便使用 者 操 作 。
Lammps 计算输入文件
• • • • • • • • • • • • • • • • # initial velocities 初始化速度 compute new mobile temp # 定义温度的计算（可动区域内统计平均） compute new2 mobile stress/atom # 定义原子应力的计算（整个区域） Velocity mobile create 0.01 887723 temp new # 按指定的温度（ 0.01 ）计算方法，初始化原子的速度 Velocity upper set 0.0 0.3 0.0 #upper 原子组 y 方向的速度为 0.3 Velocity mobile ramp vy 0.0 0.3 y 1.25 38.75 sum yes #mobile 原子的速初始度从 0 到 0.3 线性变化 # fixe boundary setforce NULL 0.0 0.0 # 边界 boundary 上力条件，钢化原子，便于加载！！
分子动力学 lammps 数据处理
分子动力学 lammps 数据处理标题：分子动力学模拟在材料科学中的应用引言：分子动力学模拟是一种常用的计算方法，可用于研究材料科学中的各种物理现象和化学反应。
本文将介绍分子动力学模拟的基本原理和在材料科学中的应用，旨在为读者提供一个全面的了解和深入的思考。
1. 分子动力学模拟的基本原理分子动力学模拟是通过数值计算来模拟原子或分子在给定条件下的运动和相互作用。
它基于牛顿运动定律和经典力学原理，将系统的势能函数和初始状态输入计算机程序中，通过数值积分方法计算出系统在不同时间点的状态。
2. 分子动力学模拟在材料科学中的应用2.1 材料的力学性能研究分子动力学模拟可以模拟材料的强度、韧性、断裂行为等力学性能。
通过模拟应力-应变曲线和材料的断裂过程，可以预测材料的力学行为，为材料设计和工程应用提供指导。
2.2 材料的热力学性质研究分子动力学模拟可以模拟材料在不同温度下的热膨胀、热传导等热力学性质。
通过模拟材料的原子振动和相互作用，可以计算材料的热导率、热膨胀系数等参数，为材料的热管理和热稳定性设计提供理论基础。
2.3 材料的界面和表面性质研究分子动力学模拟可以模拟材料的界面和表面性质，如界面结合能、表面能等。
通过模拟原子的迁移和重新排列，可以计算界面和表面的结构稳定性和能量变化，为材料的涂层、薄膜和纳米结构设计提供指导。
2.4 材料的化学反应研究分子动力学模拟可以模拟材料的化学反应，如催化反应、电化学反应等。
通过模拟反应物的结构和能量变化，可以计算反应的速率常数和反应机理，为材料的催化和能源转化应用提供理论支持。
3. 结论分子动力学模拟在材料科学中具有广泛的应用前景，可以预测材料的力学性能、热力学性质、界面和表面性质以及化学反应。
通过模拟计算，可以加深对材料微观结构和相互作用的理解，为材料设计和应用提供重要的指导和优化方案。
随着计算机计算能力的不断提升，分子动力学模拟将在材料科学研究中发挥越来越重要的作用。
lammps 蒙特卡洛化学反应
lammps 蒙特卡洛化学反应LAMMPS（Large-scale Atomic/Molecular Massively Parallel Simulator）是一种用于分子动力学模拟的开源软件，它可以模拟多种化学反应，其中包括蒙特卡洛化学反应。
本文将介绍LAMMPS 中的蒙特卡洛化学反应模拟方法及其应用。
蒙特卡洛方法是一种基于随机抽样的数值计算方法，可以用来模拟和分析复杂的化学反应过程。
在LAMMPS中，蒙特卡洛化学反应模拟是通过在分子动力学模拟中引入蒙特卡洛步骤来实现的。
在蒙特卡洛化学反应模拟中，分子系统中的每个分子都被赋予一个随机的能量状态，并通过随机抽样的方式来选择反应类型和反应路径。
通过在模拟过程中不断更新分子的能量状态和反应路径，可以模拟出分子间的化学反应。
在LAMMPS中，蒙特卡洛化学反应模拟的基本步骤如下：1. 初始化系统：包括定义分子的初始位置和能量状态。
2. 选择反应类型：根据反应的类型和反应路径，通过随机抽样的方式选择要进行的反应。
3. 计算反应概率：根据选择的反应类型和反应路径，计算反应的概率。
4. 进行反应：根据计算得到的反应概率，决定是否进行反应。
如果进行反应，则更新分子的能量状态和反应路径。
5. 更新分子状态：根据反应的结果，更新分子的位置和能量状态。
6. 重复步骤2至步骤5，直到达到模拟的时间或步数。
蒙特卡洛化学反应模拟在化学研究中具有广泛的应用。
它可以用来研究化学反应的速率常数、平衡常数以及反应路径等。
通过模拟不同的反应条件和反应路径，可以得到一系列的化学反应数据，进而对实验中观测到的化学反应进行解释和预测。
蒙特卡洛化学反应模拟还可以用来研究分子间的相互作用和聚集行为。
通过模拟不同的分子结构和反应条件，可以了解分子在不同环境中的行为，并且可以用于设计新的分子材料和催化剂。
蒙特卡洛化学反应模拟是一种强大的工具，可以用来模拟和分析复杂的化学反应过程。
通过LAMMPS软件的应用，可以实现高效的蒙特卡洛化学反应模拟，并为化学研究提供重要的理论指导和预测能力。
lammps 算例
LAMMPS 算例简介LAMMPS（Large-scale Atomic/Molecular Massively Parallel Simulator）是一个用于模拟原子、分子和离子的分子动力学程序。
它可以模拟各种材料的性质和行为，包括固体、液体和气体。
LAMMPS 使用分子动力学方法，通过模拟原子之间的相互作用和运动来研究材料的宏观性质。
本文将介绍一个 LAMMPS 的算例，以帮助读者理解如何使用 LAMMPS 进行分子动力学模拟，并展示一些常见的应用场景。
算例背景在这个算例中，我们将模拟一个固体材料的拉伸变形过程。
我们将使用 LAMMPS 来模拟原子之间的相互作用和运动，并观察材料在不同应变下的力学响应。
算例步骤1. 准备输入文件首先，我们需要准备一个输入文件，该文件包含了模拟所需的参数、原子坐标和相互作用势函数。
下面是一个示例输入文件的内容：# 输入文件示例# 设置模拟的尺寸和周期性边界条件dimension 3boundary p p punits metalatom_style atomic# 设置原子类型和质量read_data datafile.dat# 设置相互作用势函数pair_style lj/cut 2.5pair_coeff 1 1 1.0 1.0 2.5# 设置模拟参数timestep 0.001thermo 100run 10002. 运行模拟接下来，我们需要运行模拟。
在命令行中输入以下命令来运行 LAMMPS：lammps -in input_file.in其中，input_file.in是我们准备的输入文件。
3. 分析结果模拟运行完成后，我们可以通过分析输出文件来获取模拟结果。
LAMMPS 会生成一个包含模拟过程中各个时间步的能量、力和原子坐标等信息的输出文件。
我们可以使用 Python 或其他数据处理工具来分析输出文件，并绘制出力学响应曲线、原子位移等结果。
算例结果下图是一个示例结果的力学响应曲线：通过模拟，我们可以观察到随着应变的增加，材料的应力也随之增加，直到达到材料的极限。
2.LAMMPS分子动力学模拟-in文件编写
二、LAMMPS分子动力学模拟-in文件编写教程1.说明：in文件是LAMMPS软件的运行程序文件，该文件程序描述了模拟需求指令。
所有模拟指令需根据LAMMPS用户手册，即LAMMPS Users Manual （）文件进行编写。
2.以下将根据一个简单案例进行in文件基本结构说明，该案例中的结构并不固定，可根据需要进行调整。
3.in文件案例：-------------------------------------------------模型基本指令设置------------------------------ # Lennard-Jones crystal （#符号表示不执行该条指令）units real （此命令用于设置模拟的单位类型，有lj or real or metal or si orcgs or electron or micro or nano多种类型，每种类型有各自的单位设定，在后续程序编写中要注意所有数据的单位）boundary p p p （该指令用于设置模型每个维度的边界类型，p为periodic边界，三个p代表x,y,z三个方向都是周期边界）atom_style full （定义在模拟中使用的原子类型，样式的选择决定了data文件中分子结构数据所包含的要素）-------------------------------------------------分子结构模型设置------------------------------#read_data X.data （读入包含lammps运行模拟所需信息的数据文件，data文件中包含了原子坐标、种类、键、角和所带电荷等信息；分子结构也可以通过set，box等指令在in文件中进行设定）read_restart poly.restart.100000 （读入前次模拟保存的运行结果文件，从中断的模拟位置重新启动模拟）---------------------------------原子间作用势类型和参数设定------------------------------ pair_style lj/cut/coul/cut 12 12 （设置用来计算原子对相互作用的势能公式）pair_coeff 1 1 10 10 （根据指定的原子对势能函数，设置势能参数）pair_coeff 2 2 100 10pair_coeff 1 2 10 10------------------------------------键、角类型和参数设定------------------------------------- bond_style harmonicbond_coeff 1 450 1.0bond_coeff 2 500 1.45angle_style harmonicangle_coeff 1 55 109.0angle_coeff 2 55 109.28dihedral_style harmonicdihedral_coeff 1 0.062 1 3dihedral_coeff 2 0.062 1 3---------------------------系统能量最小化方法和参数设定-------------------------------- #min_style sd （选择执行最小化命令时要使用的最小化算法）#minimize 1.0e-5 1.0e-5 100 100 （通过迭代调整原子坐标，实现系统的能量最小化，该指令设置迭代终止条件）--------------------------------------其它模拟相关指令设定----------------------------------- #velocity all create 300.0 200000 （设置或改变一组原子的速度）#velocity all scale 300.0fix 1 all nvt temp 873.0 873.0 1 （设置NVT系综）#fix 2 all temp/rescale 50 673 673 10 1.0 （通过重新调整原子群的速度来重置原子群的温度）compute KE all ke/atom （对一组原子执行计算）variable temp atom c_KE/0.0001292355 （此命令将数值或公式计算结果指配给变量名，以便稍后在输入脚本或模拟过程中使用该变量进行计算）fix 6 all ave/time 10 10000 100000 v_temp file tem.profile （输出时间平均计算结果，写入一个命名为tem.profile的文件）-----------------------------------模拟结果输出相关指令设定------------------------------- timestep 1 （设置分子模拟的时间步长大小）thermo_style custom time temp press density pe （设置将热力学数据打印到屏幕和日志文件的样式和内容）thermo_modify lost ignore flush yesthermo 50restart 100000 poly.restart （每隔这么多个时间步写出一个包含当前模拟数据的重新启动文件）dump 1 all atom 100000 mmpstrj （每100000个时间步将Atom数据转储到mmpstrj文件）run 50000000 （指定运行的步数）。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
Lammps计算输入文件
• pair_style eam • #定义原子之间相互作用势 • pair_coeff * * Ni_u3.eam • #定义作用势的参数，通过文件输入 • Neighbor 2.0 bin • #定义近邻表参数(r+rc),rc=2.0埃 • neigh_modify delay 5 • #定义更新近邻表的加载步间隔
B、使用L-J势模拟裂纹的扩展
leftupper 裂纹
upper
y
leftlower
lower x
Lammps计算输入文件
• # 2d LJ crack simulation(问题的基本初始化)
• dimension 2 • #2维分子动力学模拟 • boundary s s p • #x,y方向不存在周期边界条件，z方向周期边界条件。
Lammps计算输入文件
• # run运行计算 • timestep 0.003 • #时间间隔步 • Thermo 200 • #每200步输出热动力学统计量 • thermo_modify temp new • #计算温度通过new指示的方法计算
• neigh_modify exclude type 2 3 • #原子2，3之间作用取消（也就是通过不使他们在近邻表中出
• Mass 1 1.0 • Mass 2 1.0 • Mass 3 1.0 • Mass 4 1.0 • Mass 5 1.0 • #指定单个原子的质量
Lammps计算输入文件
• # LJ potentials(指定原子作用势) • pair_style lj/cut 2.5 • #指定lj势，截断半径为2.5 • pair_coeff * * 1.0 1.0 2.5 • #指定lj势参数
现实现）
• Dump 1 all atom 500 dump.crack • #每隔500步将原子信息写入文件dump.crack
• Dump 2 mobile custom 500 dump2.crack tag x y z c_new2[2]
• run
5000
• #进行5000步的模拟
C.使用EAM势模拟Ni的循环行为
• # define groups（便于加载） • Region 1 block INF INF INF 1.25 INF INF • Group lower region 1 • #定义lower组（便于施加外加速度） • Region 2 block INF INF 38.75 INF INF INF • Group upper region 2 • #定义upper组（便于施加外加速度） • Group boundary union lower upper • #定义总边界组 • Group mobile subtract all boundary • #定义可动原子组（便于统计温度）
Lammps计算输入文件
• #弛豫和循环加载 • run 2000 • #不加任何载荷计算2000步 • velocity upper set 0 0 0.4 • velocity lower set 0 0 -0.4 • run 10000
• velocity upper set 0 0 -0.4 • velocity lower set 0 0 0.4 • run 20000
Lammps计算输入文件
• Lattice fcc 3.52 orient x 1 0 0 orient y 0 1 1 orient z & 0 -1 1 origin 0.5 0 0
• #指定晶格类型、晶格常数、x,y,z坐标的晶向 • # &续行符、晶胞的起始点（0.5表示半个晶格常数） • create_atoms 1 box • #在box区域内建立类型为1的原子 • pair_style eam • #定义原子之间作用势为嵌入原子势 • pair_coeff * * Ni_u3.eam • #指定嵌入势参数文件
z
u
u
y
x
u
常温30K条件下
金属材料模拟中Lammps的单位
Lammps计算输入文件
• # 3d metal 拉伸模拟 • #模拟条件的初始化 • Units metal • #指定模拟中的单位类型 • boundary s s s • #指定模拟的边界条件 • atom_style atomic • #指定原子类型，原子的属性 • Lattice fcc 3.52 • #指定材料的晶格类型和晶格常数 • region box block -5 5 -5 5 -15 15 • #xlo,xhi,ylo,yhi,zlo,zhi =区域box的上下限 • create_box 1 box • #建立只有1中原子类型的simulation box • create_atoms 1 box • #在simulation box中创建类型为1的原子（原子位置初始化）
Lammps计算输入文件
• region lower block INF INF INF INF INF -13 • #生成lower区域，便于加载 • region upper block INF INF INF INF 13 INF • #生成upper区域，便于加载 • group lower region lower • #定义lower组 • group upper region upper • #定义upper组 • group boundary union lower upper • #定义boundary组 • group mobile subtract all boundary • #定义mobile组，便于计算统计温度（调温）
• # initial velocities初始化速度 • compute new mobile temp • #定义温度的计算（可动区域内统计平均） • compute new2 mobile stress/atom • #定义原子应力的计算（整个区域） • Velocity mobile create 0.01 887723 temp new • #按指定的温度（0.01）计算方法，初始化原子的速度 • Velocity upper set 0.0 0.3 0.0 • #upper原子组y方向的速度为0.3 • Velocity mobile ramp vy 0.0 0.3 y 1.25 38.75 sum yes • #mobile原子的速初始度从0到0.3线性变化 • # fixes施加约束 • fix 1 all nve • #nve系综的积分算法 • fix 2 boundary setforce NULL 0.0 0.0 • #边界boundary上力条件，钢化原子，便于加载！！
• #定义左上、左下原子组（便于指定裂纹的存在）
• set
group leftupper type 2
• set
group leftlower type 3
• set
group lower type 4
• set
group upper type 5
• #指定原子类型（便于指定裂纹的存在）
Lammps计算输入文件
z <011>
Z‘ <001>
<011> y
<010> Y’
<100> x
Lammps计算输入文件
• # 3d metal shear simulation 3维金属剪切模拟
• Units metal • #采用金属材料单位 • boundary s s p • #施加z方向的周期边界条件
• atom_style atomic • #定义原子之间的相互作用（没有键键作用） • Lattice fcc 3.52 • #指定晶格为fcc，晶格常数3.52埃 • Region box block 0 16.0 0 10.0 0 2.828427 • #定义区域box （x上下限，y上下限，z上下限） • create_box 3 box • #在box区域内创建包含3种原子类型的simulation box
为随机数
Lammps计算输入文件
• dump 2 all custom 200 d_t.txt tag x y z c_new2 • #定义输出，每两百步输出一次 • #tag为原子标号 • #x，y，z为原子坐标 • #c_new2为new2所计算的原子适配度 • #这一参数可以判断原子是否是理想晶格 • #或者在位错和堆垛层错上 • timestep 0.001 • #时间步长间隔（ps=10e-12秒） • #时间步长在飞秒量级10e-15秒 • #小于原子自振周期的十分之一（0.1ps/10）
• atom_style atomic • #原子类型（金属） • Neighbor 0.3 bin • #建近邻表参数（rc+0.3）bin表示为近邻表类型。 • neigh_modify delay 5 • #间隔多少载荷步重新形成近邻表
Lammps计算输入文件
• # create geometry创建初始几何构形
• Lattice hex 0.93 • #指定晶格类型（二维hex）和晶格常数 • Region box block 0 100 0 40 -0.25 0.25 • #定义一个区域 • create_box 5 box • #在指定区域建立一个simulation box,5表示原子类型的种类数 • create_atoms 1 box • #在simulation box中创建类型为1的原子（原子位置初始化）
Lammps计算输入文件
• #定义控温和原子速度初始化 • fix 4 mobile temp/rescale 10 30.0 30.0 5.0 1.0 • #使用直接调温法调整mobile原子的速度，每10步一调 • #使得目标温度为30+-1K，温度在30+-5rescale之外调 • fix_modify 4 temp new3 • #调温时，计算温度采用news计算 • Velocity mobile create 30 887723 temp new3 • #初始化mobile group中原子的初速度 • #通过温度30K初始化，温度通过new3计算，887723