Lammps大学物理开放性实验PPT.pptx

大学物理开放性实验铝的熔化转变1.Lampps的安装和使用查看电脑的属性，安装64位Lampps。安装路径C:/Program files/Lammps 64-bit打开我的电脑，输入cmd，打开管理员窗口(C：\Windows\System32\cmd.exe)。 输入cd..进入前一个文件，cd为当前文件，dir查看当前文件里都有哪些子文件。 运行一个输入文件时，例如melt文件，melt文件在Examples将melt文件复制到bin文件里，管理员控制窗口运行到C:/Program files/Lammps 64-bitbin；之后输入lmp_serial -in in.melt之后回车（注意输入的时候空格键），文件输出运行结果。2.lammps基本知识的了解Fix为定义对部分原子或者对系统的操作。Fix-npt：恒定对N（原子数量），P（压强），T（温度）对时间积分Fix-nvt：恒定对N，V（体积），T对时间积分?drag = style name of this fix command?x,y,z = coord to drag atoms towards3.使用范例In.melt输入文件# 3d Lennard-Jones meltunits lj # 单位，指定为lammps 里的非金属类的单位，长度为?，能量为eV。atom_style atomic # 原子模式lattice fcc 0.8442 # 晶格常数0.8442region box block 0 10 0 10 0 10 # x,y,z 各方向上的晶胞重复单元数，也即区域大小create_box 1 box # 将上述区域指定为模拟的盒子create_atoms 1 box # 将原子按晶格填满盒子mass 1 1.0velocity all create 3.0 87287pair_style lj/cut 2.5pair_coeff 1 1 1.0 1.0 2.5neighbor 0.3 binneigh_modify every 20 delay 0 check nofix 1 all nve#dump id all atom 50 dump.melt#dump 2 all image 25 image.*.jpg type type # axes yes 0.8 0.02 view 60 -30#dump_modify 2 pad 3#dump 3 all movie 25 movie.mpg type type # axes yes 0.8 0.02 view 60 -30#dump_modify 3 pad 3thermo 50 #每隔50 步输出热力学结果run 250 #运行250 步in.melt文件输出结果数据对应的含义数据第一列是步数（step）第二列是温度（temperature)第三列是相互作用能量（pairwise energy）第五列是总能量（total energy）第六列是压力（press）。用数据绘制表格铝的熔化转变: 对于铝，LAMMPS 建立 8×8×5的 FCC 晶格体系;充分弛豫后利用 N持压强为零，使体系从 T=2.5K开始加热，直至发生熔化转变。下面是铝熔化的输入文件： 在运行过程中遇到的问题，与旧版不同，在新版中，运行基本命令fix时，需要加入关键字。 例如：fix 1 all nvt temp $x $x 1.0 drag 0.2 这个语句中nvt后面需要加入temp，temp是keyword； 而在 fix 1 all npt temp couple iso $x 1500 4.00 xyz 0.0 中，在npt后面加入temp couple iso； 在运行这个文件时需要EAM势文件，在potentials中寻找合适的势文件例如本实验用的是 Al_mm.eam.fs。数据导入的方法在管理员窗口全选，复制到记事本里，另存为数据*txt文件打开EXCEl，选择数据窗口，导入数据，选择数据源，选择数据*txt文件。选定后，选择需要的数据进行数据处理。对于铝，其熔点为 934.477K, 在我们的模拟中，其在 1110.0K附近发生一级相变，单位体积发生突变。其相对于平衡熔点增大了了 18.9%。范例上面为范例的图片下面是自己拟合