分享乙個團簇融化的lammps模擬案例,由carsten svaneborg博士編寫。
模擬單位為lj,2d模擬體系,週期性邊界條件,力場為最常見的lj力場。
模擬思路:
(1)設定乙個4040的模擬盒子,在盒子**設定乙個2020的區域,在中間區域內隨機填充原子。
(2)能量最小化消除隨機原子的重疊。
(3)設定lj力場
(4)設定nve系綜,使用langevin控溫方式進行控溫。這一步也可以直接用nvt系綜代替。
(5)執行並輸出軌跡檔案。
案例比較簡單,需要注意的是:採用隨機方法生成原子,原子之間容易重疊,一定要在正式模擬之前消除原子重疊,否則原子距離過近導致排斥力過大,容易產生「lost atoms」錯誤。
消除重疊可使用「delete_atoms verlap」和「minimize」命令。
案例**已注釋,僅供參考。
#原子數量
variable npart equal 500
#lj單位溫度
variable tempstart equal 0.5
variable tempstop equal 0.5
#模型基本設定
units lj
dimension 2
atom_style atomic
boundary p p p
neighbor 0.5 bin
neigh_modify every 1 delay 0 check yes
#模擬盒子設定
region box block -
2020-20
20-0.10.1
create_box 1 box
#2維計算設定
fix 2d all enforce2d
#在box**設定box2
region box2 block -
1010-10
10-0.10.1
#在box2中隨機生成原子
create_atoms 1 random $
324523 box2
#摩爾質量
mass *
1#lj力場設定
pair_style lj/cut 2.5
pair_coeff *
*1.0
1.0#能量最小化
minimize 1e-4
1e-4
1000
1000
reset_timestep 0
#儲存原子軌跡
dump img all atom 1000 dump.xyz
#nve系綜,langevin控溫
fix integrator all nve
fix dynamics all langevin $ $
1.02352325
#設定螢幕輸出熱力學資訊
thermo_style custom step temp ke pe
thermo 100
#模擬步長設定
timestep 0.01
#執行20000步
run 20000
融雲專案案例開發教程
先來理下思路,對比前面的幾個教程,我們時這樣實現的,在訊息列表頁我們連線一次融雲,開啟會話頁面後又載入了一次融雲,這樣開啟會話頁面後就會導致訊息列表頁的監聽事件取消,同時還會存在一些位置的bug,導致了開發難度。我們這次改變下思路,把所有的融雲類的處理全部放在訊息列表頁,然後通過apicloud的a...
lammps案例 石墨烯捲曲到Fe奈米線
石墨烯是一種二維材料,力學效能優異,自被發現以來一直是研究的熱點。石墨烯在范德華力的作用下,可自動的包覆到其他材料,實現自動組裝。這篇推文復現 diverse nanowires activated self scrolling of graphene nanoribbons 中模擬方法 對作者表示...
分形 DLA團簇模型的摸擬
分形 dla 團簇模型的摸擬 emilmatthew emilmatthew 126.com 1前言 dla 團簇模型是分形成長中的乙個很易於實現的模型 當然啦,背後的理論工作是非常高有難度的,涉及到代數,泛函分析,非線性偏微分方程等知識。我們就拿著別人的研究出的成果,來感受一下其中的美妙吧j 團簇...