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