Lammps大學(xué)物理開放性實(shí)驗(yàn)PPT

大學(xué)物理開放性實(shí)驗(yàn)鋁的熔化轉(zhuǎn)變1.Lampps的安裝和使用查看電腦的屬性，安裝64位Lampps。安裝路徑C:/Programfiles/Lammps64-bit20151116.打開我的電腦，輸入cmd，打開管理員窗口(C：\Windows\System32\cmd.exe)。輸入cd..進(jìn)入前一個(gè)文件，cd為當(dāng)前文件，dir查看當(dāng)前文件里都有哪些子文件。運(yùn)行一個(gè)輸入文件時(shí)，例如melt文件，melt文件在Examples將melt文件復(fù)制到bin文件里，管理員控制窗口運(yùn)行到C:/Programfiles/Lammps64-bit20151116/bin；之后輸入lmp_serial-inin.melt之后回車（注意輸入的時(shí)候空格鍵），文件輸出運(yùn)行結(jié)果。2.lammps基本知識(shí)的了解Fix為定義對(duì)部分原子或者對(duì)系統(tǒng)的操作。Fix-npt：恒定對(duì)N（原子數(shù)量），P（壓強(qiáng)），T（溫度）對(duì)時(shí)間積分Fix-nvt：恒定對(duì)N，V（體積），T對(duì)時(shí)間積分drag=stylenameofthisfixcommandx,y,z=coordtodragatomstowards3.使用范例In.melt輸入文件#3dLennard-Jonesmeltunits lj#單位，指定為lammps里的非金屬類的單位，長(zhǎng)度為?，能量為eV。atom_style atomic#原子模式lattice fcc0.8442#晶格常數(shù)0.8442region boxblock010010010#x,y,z各方向上的晶胞重復(fù)單元數(shù)，也即區(qū)域大小create_box 1box#將上述區(qū)域指定為模擬的盒子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步輸出熱力學(xué)結(jié)果run 250#運(yùn)行250步in.melt文件輸出結(jié)果數(shù)據(jù)對(duì)應(yīng)的含義數(shù)據(jù)第一列是步數(shù)（step）第二列是溫度（temperature)第三列是相互作用能量（pairwiseenergy）第五列是總能量（totalenergy）第六列是壓力（press）。用數(shù)據(jù)繪制表格鋁的熔化轉(zhuǎn)變:對(duì)于鋁，LAMMPS建立8×8×5的FCC晶格體系;充分弛豫后利用N持壓強(qiáng)為零，使體系從T=2.5K開始加熱，直至發(fā)生熔化轉(zhuǎn)變。下面是鋁熔化的輸入文件：

在運(yùn)行過程中遇到的問題，與舊版不同，在新版中，運(yùn)行基本命令fix時(shí)，需要加入關(guān)鍵字。例如：fix1allnvttemp$x$x1.0drag0.2這個(gè)語(yǔ)句中nvt后面需要加入temp，temp是keyword；而在

fix1allnpttemp

coupleiso$x15004.00xyz0.0中，在npt后面加入tempcoupleiso；在運(yùn)行這個(gè)文件時(shí)需要EAM勢(shì)文件，在potentials中尋找合適的勢(shì)文件例如本實(shí)驗(yàn)用的是Al_mm.eam.fs。數(shù)據(jù)導(dǎo)入的方法在管理員窗口全選，復(fù)制到記事本里，另存為數(shù)據(jù)*txt文件打開EXCEl，選擇數(shù)據(jù)窗口，導(dǎo)入數(shù)據(jù)，選擇數(shù)據(jù)源，選擇數(shù)據(jù)*txt文件。選定后，選擇需要的數(shù)據(jù)進(jìn)行數(shù)據(jù)處理。對(duì)于鋁，其熔點(diǎn)為934.477K,在我們的模擬中，其在1110.0K附近發(fā)生一級(jí)相變，單位體積發(fā)生突變。其相對(duì)于平衡熔點(diǎn)增大了了18.9%。范例上面為范例的圖片下面是自己擬合的圖片發(fā)生一級(jí)相變的溫度1099.95