lammps案例 團簇融化過程模擬案例學習

2021-10-25 17:07:53 字數 1715 閱讀 7843

分享乙個團簇融化的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 團簇...