LAMMPS手冊-中文版講解

LAMMPS手冊-中文解析之五兆芳芳創(chuàng)作一、簡介本部分大至介紹了LAMMPS的一些功效和缺陷.什么是LAMMPS?LAMMPS是一個經(jīng)典的份子動力學代碼，他可以模擬液體中的粒子，固體和汽體的系綜.他可以采取不合的力場和鴻溝條件來模擬全原子，聚合物，生物，金屬，粒狀和粗料化體系.LAMMPS可以計較的體系小至幾個粒子，大到上百萬甚至是上億個粒子.LAMMPS可以在單個處理器的臺式機和筆記本本上運行且有較高的計較效率，但是它是專門為并行計較機設計的.他可以在任何一個按裝了C++編譯器和MPI的平臺上運算，這其中當然包含散布式和同享式并行機和Beowulf型的集群機.LAMMPS是一可以修改和擴展的計較程序，比方，可以加上一些新的力場，原子模型，鴻溝條件和診斷功效等.通常意義上來講，LAMMPS是按照不合的鴻溝條件和初始條件對通太短程和長程力相互作用的份子，原子和宏不雅粒子薈萃對它們的牛頓運動方程進行積分.高效率計較的LAMMPS通過采取相鄰清單來跟蹤他們鄰近的粒子.這些清單是按照粒子間的短程互拆力的大小進行優(yōu)化過的，目的是避免局部粒子密度太高.在并行機上，LAMMPS采取的是空間分化技巧來分派模擬的區(qū)域，把整個模擬空間分紅較小的三維小空間，其中每一個小空間可以分派在一個處理器上.各個處理器之間相互通信并且存儲每一個小空間鴻溝上的《卜^二原子的信息.LAMMPS（并行情況）在模擬3維矩行盒子并且具有近均一密度的體系時效率最高.LAMMPS的功效總體功效：可以串行和并行計較散布式MPI戰(zhàn)略模擬空間的分化并行機制開源高移植性C++語言編寫MPI和單處理器串行FFT的可選性（自定義）可以便利的為之擴展上新特征和功效只需一個輸入腳本就可運行有定義和使用變量和方程完備語律例則在運行進程中循環(huán)的控制都有嚴格的法則只要一個輸入腳本試就可以同時實現(xiàn)一個或多個模擬任務粒子和模擬的類型：（atomstyle命令）原子粗?；Ｗ尤泳酆衔铮袡C份子，蛋白質(zhì)，DNA聯(lián)合原子聚合物或有機份子金屬粒子資料粗?；椴谎拍Ｐ脱由烨蛐闻c橢圓形粒子點偶極粒子剛性粒子所有上面的雜化類型力場：(命令：pairstyle,bondstyle,anglestyle,dihedralstyle,improperstyle,kspacestyle)對相互作用勢：L-J,Buckingham,Morse,Yukawa,soft,class2(COMPASS),tabulated.帶點對相互作用勢：Coulombic,point-dipole.多體作用勢：EAM,Finnis/SinclairEAM,modifiedEAM(MEAM),Stillinger-Weber,Tersoff,AIREBO,ReaxFF粗?；饔脛荩篋PD,GayBerne,Resquared,Colloidal,DLVO介不雅作用勢：granular,Peridynamics鍵勢能：harmonic,FENE,Morse,nonlinear,class2,quartic鍵角勢能：harmonic,CHARMM,cosine,cosine/squared,class2(COMPASS)二面角勢能：harmonic,CHARMM,multi-harmonic,helix,OPLS,class2(COMPASS)不合理勢能：harmonic,CVFF,class2(COMPASS)聚合物勢能：all-atom,united-atom,bead-spring,breakable水勢能：TIP3P,TIP4P,SPC隱式溶劑勢能：hydrodynamiclubrication,Debye長程庫倫與分離：Ewald,PPPM,Ewald/N（針對長程L-J作用）可以有與普適化力場如CHARMM，AMBER，OPLS，GROMACS相兼容的力場可以采取GPU加快的成對類型雜化勢能函數(shù)：multiplepair,bond,angle,dihedral,improperpotentials^對勢能處于更高的優(yōu)先級）原子創(chuàng)建：（命令：read_data,lattice,create-atoms,delete-atoms,displace-atoms,replicate）從文件中讀入各個原子的坐標在一個或多個晶格中創(chuàng)建原子刪除幾何或邏輯原子基團復制已存在的原子多次替換原子系綜，約束條件，鴻溝條件：（命令：fix）二維和三維體系正角或非正角模擬空間常NVE，NVT，NPT，NPH積分器原子基團與幾何區(qū)域可選擇不合的溫度控制器有Nose/Hoover和Berendsen壓力控制器來控制體系的壓力（任一維度上）模擬合子的變形（扭曲與剪切）簡諧（unbrella）束縛力剛體約束搖擺鍵與鍵角約束各類鴻溝情況非平行太份子動力學NEMD各類附加鴻溝條件和約束積分器：Velocity-verlet積分器Brown積分器rRESPA承繼時間延化積分器剛體積分器共軛梯度或最束下降算法能量最小化器輸出：（命令：dump,restart）熱力學信息日志原子坐標，速度和其它原子量信息的文本dump文件二進制重啟文件各原子量包含：能量，壓力，中心對稱參數(shù)，CAN等用戶自定義系統(tǒng)寬度或各原子的計較信息每個原子的時間與空間平均系統(tǒng)寬量的時間平均原子圖像，XYZ,XTC,DCD,CFG格局數(shù)據(jù)的前處理與后處理：包里提供了一系列的前處理與后處理東西另外，可以使用獨立刊行的東西組pizza.py它可以進行LAMMPS模擬的設置，闡發(fā)，作圖和可視化任務.特別功效：實時的可視化與交互式MD模擬與有限元辦法結(jié)合進行原子-連續(xù)體模擬在POEMS庫中提供了剛體積分東西并行裉火并行復制動力學對低密度液體直接使用MC模擬Peridynamic介不雅建模目標型與無目標型份子動力學雙溫度電子模型LAMMPS不具備的功效：由于LAMMPS是對牛頓運動方程積分的東西，所以良多需要的數(shù)據(jù)前處理與后處理功效是LAMMPS焦點不具備的.其原因為：包管LAMMPS的小巧性前處理與后處理不克不及進行并交運算這些功效可以有其它東西來完成原代碼開發(fā)的局限性特別地，LAMMPS不克不及：通過圖形用戶界面來任務創(chuàng)建份子體系自動的加上力場系數(shù)為MD模擬提供智能化的數(shù)據(jù)闡發(fā)MD的可視化為輸出數(shù)據(jù)作圖我們需要為LAMMPS輸入一系列的原子類型，原子坐標，份子拓樸信息和所有原子與鍵的力場參數(shù).LAMMPS不會自動的為我們創(chuàng)建份子體系與力場參數(shù).對與原子體系，LAMMPS提供了creat-atoms命令來為固態(tài)晶格加上原子.可以能過paircoeff,bondcoeff,anglecoeff等命來加上小數(shù)目的力場參數(shù).對于份子體系或更龐雜的模擬體系，我們通常會用其它東西來創(chuàng)建或是轉(zhuǎn)換LAMMPS輸出文件來做到這些事情.有的還會寫一些自已的代碼來完成這項任務.對于一個龐雜的份子體系（如，蛋白質(zhì)），我們需要為之提供上面?zhèn)€拓樸信息與力場參數(shù).所以我們建議用CHARMM或AMBER或其它的份子建模器來完成這些任務，并把之輸?shù)揭粋€文件中去.然后，改動其格局以達到LAMMPS所允許的輸入格局.同樣，LAMMPS的輸出文件是一種復雜的文本格局，我們也可以通過其它的東西來換專這些格局.我們可以用以下幾個軟件來完成高質(zhì)量的可視任務：VMDAtomEyePymolRaster3dRasMol最后要說一下的是，以下這些也是自由份子動力學包，它們大多數(shù)是并行的，可能也適合來完成你的研究任務，當然也可以與LAMMPS聯(lián)合起來使用以完成模擬任務.CHARMMAMBERNAMDNWCHEMDL_POLYTinkerCHARMM，AMBER，NAMD，NWCHEM，Tinker是專們用于模擬生物份子的.二開始、本部分主要描述如何創(chuàng)建和運行LAMMPS.1在LAMMPS刊行包理含有：．READMELICENSEBench:測式任務口05文本Examples:復雜的測試任務Potentials:嵌入原子辦法與力場文件$^：源代碼Tools:前處理與后處理東西假設你下載的是windows可執(zhí)行文件的話，你里面只有一個文件（并行與非并行兩種）L2.編譯LAMMPS之前的任務:編譯LAMMPS不是一個繁瑣的任務.首先你可能要寫一個makefile文件，里面要選擇編譯器，附加的一些將要用到的庫等.事先裝上MPI或FFT等庫.編譯出一個可執(zhí)行LAMMPS：在SRC目錄里頭含有C++源文件和頭文件.當然也包含一個高水平的Makefile，在MAKE目錄里頭有幾個低水平的Makefile.*files辨別適有不合的平臺.進入SRC目錄，輸入make或gmake，你將會看到一列的可選項.假設其中有一種合適你的機械，你可以輸入像下面一樣的命令：MakelinuxGmakemac注意，在一個多處理器或多核處理器的平臺上你可以進行平行編譯，在make命令中使用“-j”選項就可以，這樣編譯起來會更快一些.在此進程中不產(chǎn)生錯誤的話，你可以得到一個類似于lmp-linux的可執(zhí)行文件.在編譯進程當中將會產(chǎn)生的罕有錯誤：（1）如果編譯進程當中產(chǎn)生錯誤，并提示不克不及找到一個含有通配符*為名的文件的話，說明你機械上的make器允許makefile中使用通配符.那就償式使用gmake.如還不成的話，就試試參加-f選項，用Makefile.list作為make對像.如：MakemakelistMake-fmakefile.listlinuxGmake-fmakefile.lstmac（2）當你使用低水平的makefile時，可能由于對機械的設置不正確，會導致一些錯誤.假設你的平臺叫江^”，的話，你將要在MAKE目錄中創(chuàng)建一個Makefile.foo.使用任何一個與你機械相近的文件作為開始總是一個不錯的選擇.（3）如你在鏈接的時候出現(xiàn)庫丟失或少了依賴關系的話，可能是由于：你編譯的包需要一個附加的庫，但卻沒有事先編譯需要的packagelibaray.你要鏈接的庫在你的系統(tǒng)中不存在.沒有連接到需要的系統(tǒng)庫后兩種問題出現(xiàn)，你就需要修改你的低水平makefile.foo.編輯一個新的低水平makefile.foo:（1）在#后的句子中，替換foo,不管你寫成什么，這一即將會出現(xiàn)在屏幕上，如果你只輸入make命令的話.（2）在“complier/linkersettings”部分為你的C++編譯器列出編譯器與鏈接器的設置，包含優(yōu)化符號.你可以在任何UNIX系統(tǒng)中使用G++編譯器.當然你也可以用MPICC，如果你的系統(tǒng)中裝置了MPI的話.如過在編譯進程當中需要符加的庫的話，你必須在LIB變量中列出來.DEPFLAGS設置可以讓C++編譯器創(chuàng)建一個源文件的依賴關系列表，當源文件或頭文件改動的時候可以放慢編譯速度.有些編譯器不克不及創(chuàng)建依賴關系列表，或你可以用選項—D來實現(xiàn).G++可以使用-D.如果你的編譯器不克不及創(chuàng)建依賴關系文件的話，那么你就需要創(chuàng)建一個Makefile.foo來與Makefile.storm（它用一系列的不需要依賴文件的計劃）相對應.（3）“system-specificsettings”部分有四個小部分：ALMP—INC變量，包含一些與系統(tǒng)相關的條件選項.B3個MPI變量用于指定MPI庫.如你要進行并行計較的話，那么你必須在你的平臺上安上MPI庫.如你想用MPI內(nèi)置C++編譯器的話，你可以讓這三個變量空著，如你不必MPICC的話，那么，你要指定MPI.h(MPI_INC)文件在哪，MPI(MPI_PATH庫在哪，還有庫名(MPI_LIB).如果你想自己裝置MPI的話，我們建議用MPICH1.2或2.0.LAMMPI也可以.如果我的是大平始的話，你的供給商已經(jīng)為你裝上了MPI,其可能比MPICH或LAM更快，你可以把找出來并與之鏈接.如你用LAM或MPICH，你必須要設置他并編譯他使之適合你的平臺.如果你想在單處理器的機械上運行的話，你可以用STIBS庫，這樣你就可以不必在你的系統(tǒng)中裝置MPI庫.防照makefile.serial，看是如果設置這三個變量的.當然你在編譯LAMMPS之前你必須創(chuàng)建STUBS庫在STUBS目錄中，輸入make，不出錯的話你將會得到一個libmpi.a文件可供鏈接到LAMMPS.當出錯，你則要修改STUBS下的MAKEFILE.STUBS/MPI.CPP有一個CPU計時器MPI_Wtime(可以調(diào)用gettimeofday()如你的系統(tǒng)不支持gettimeofday(),則你就要拔出一句代碼來調(diào)用另一個計時器，要注意的是，。1。。的函數(shù)在一個小時之后會歸0，所以對于一個長時間的LAMMPS模擬來說這是不敷用的.CFTT變量用于指定FFT庫，當要用到kspace-style命令來計較長程庫倫作用時使用PPPM選項時要用到.要使用此選項，你必須要在你的機械上裝置一個一維的FFT庫.可以能過開關一DFFT—XXX來指定，其中XXX=INTELDEC，SGI，SCSL，或尸尸丁亞.沒有辦法的情況下可以用供給商提供的庫下尸丁亞是一個快速的，可移植性的庫，它可以在任何一個平臺上運行.最好2.1.X本.編譯FFTW庫時只要ffi./configure;make就可以.不任你是用哪一種FFT庫，你都要在makefile.foo中正確的設置咱們的FFT_INC,FFT_PATH,FFT_LIB.當然，你如果不必PPPM的話，你將沒有需要裝置FFT庫.這種情況下，你可以把FFT_INC設成-DFFT-NONE并讓其它幾個變量空著.你也可以在編譯LAMMPS時把KSPACE包剔除.D幾個SYSLIB和SYSPATH變量你可以疏忽，除非你在編譯LAMMPS時其中有一個或幾個包要用到附加的系統(tǒng)庫.所有這些包都的的名稱都將會是SYSLIB和SYSPATH變量的前輟.SYSLIB變量將列出系統(tǒng)庫.SYSPATH則是路徑，只有當這些庫為非默認路徑時才有設定.最后，當你正確的寫好了makefile.foo和預編譯好了所有的其它庫（MPI,FFT,包庫等）之后，你只要在SRC目錄下輸入下面其中一個命令就可以了MakefooGmakefoo不出意外，你將會得到lmp_foo的可執(zhí)行文件.附加建義（1）為多平臺編譯LAMMPS你可以在同一個SRC目錄下為多平臺編譯LAMMPS.每一個目標都有他自己的目標路徑，Obj_name用于存貯指定系統(tǒng)的目標文件.2）清理輸入makeclear-a