Lammps大学物理开放性实验PPT

大学物理开放性实验铝的熔化转变1.Lampps的安装和使用查看电脑的属性，安装64位Lampps。安装路径C:/Programfiles/Lammps64-bit20151116.打开我的电脑，输入cmd，打开管理员窗口(C：\Windows\System32\cmd.exe)。输入cd..进入前一个文件，cd为当前文件，dir查看当前文件里都有哪些子文件。运行一个输入文件时，例如melt文件，melt文件在Examples将melt文件复制到bin文件里，管理员控制窗口运行到C:/Programfiles/Lammps64-bit20151116/bin；之后输入lmp_serial-inin.melt之后回车（注意输入的时候空格键），文件输出运行结果。2.lammps基本知识的了解Fix为定义对部分原子或者对系统的操作。Fix-npt：恒定对N（原子数量），P（压强），T（温度）对时间积分Fix-nvt：恒定对N，V（体积），T对时间积分drag=stylenameofthisfixcommandx,y,z=coordtodragatomstowards3.使用范例In.melt输入文件#3dLennard-Jonesmeltunits lj#单位，指定为lammps里的非金属类的单位，长度为Å，能量为eV。atom_style atomic#原子模式lattice fcc0.8442#晶格常数0.8442region boxblock010010010#x,y,z各方向上的晶胞重复单元数，也即区域大小create_box 1box#将上述区域指定为模拟的盒子create_atoms 1box#将原子按晶格填满盒子mass 11.0velocity allcreate3.087287pair_style lj/cut2.5pair_coeff 111.01.02.5neighbor 0.3binneigh_modify every20delay0checknofix 1allnve#dump idallatom50dump.melt#dump 2allimage25image.*.jpgtypetype&# axesyes0.80.02view60-30#dump_modify 2pad3#dump 3allmovie25movie.mpgtypetype&# axesyes0.80.02view60-30#dump_modify 3pad3thermo 50#每隔50步输出热力学结果run 250#运行250步in.melt文件输出结果数据对应的含义数据第一列是步数（step）第二列是温度（temperature)第三列是相互作用能量（pairwiseenergy）第五列是总能量（totalenergy）第六列是压力（press）。用数据绘制表格铝的熔化转变:对于铝，LAMMPS建立8×8×5的FCC晶格体系;充分弛豫后利用N持压强为零，使体系从T=2.5K开始加热，直至发生熔化转变。下面是铝熔化的输入文件：

在运行过程中遇到的问题，与旧版不同，在新版中，运行基本命令fix时，需要加入关键字。例如：fix1allnvttemp$x$x1.0drag0.2这个语句中nvt后面需要加入temp，temp是keyword；而在

fix1allnpttemp

coupleiso$x15004.00xyz0.0中，在npt后面加入tempcoupleiso；在运行这个文件时需要EAM势文件，在potentials中寻找合适的势文件例如本实验用的是Al_mm.eam.fs。数据导入的方法在管理员窗口全选，复制到记事本里，另存为数据*txt文件打开EXCEl，选择数据窗口，导入数据，选择数据源，选择数据*txt文件。选定后，选择需要的数据进行数据处理。对于铝，其熔点为934.477K,在我们的模拟中，其在1110.0K附近发生一级相变，单位体积发生突变。其相对于平衡熔点增大了了18.9%。范例上面为范例的图片下面是自己拟合的图片发生一级相变的温度1099.95