| 查看: 2222 | 回復(fù): 18 | ||||
| 本帖產(chǎn)生 1 個 模擬EPI ,點(diǎn)擊這里進(jìn)行查看 | ||||
xiaowu787木蟲 (正式寫手)
|
[交流]
【討論】高溫高壓的體系需要注意的環(huán)節(jié) 已有4人參與
|
|||
|
希望LAMMPS模擬的模型中心產(chǎn)生一個高溫高壓(有時間壽命)的局部區(qū),然后觀測溫度和壓力的向周圍擴(kuò)展的情況。可是溫度高壓強(qiáng)高時會出現(xiàn)lost atoms,現(xiàn)在 第一個疑問是:周期邊界條件下的image atom怎么不起作用了--周期性的邊界條件不是可以保證粒子數(shù)不變? 第二: 局部升溫的選擇: 1. fix all nve+fix 局部 temp/rescale,可是壓強(qiáng)無法控制 2. fix all nve+fix 局部 nvt or npt? 符合實際情況嗎? 3. fix all npt+fix 局部 nvt or npt? 符合實際情況嗎? 4. fix all nvt+fix 局部 nvt or npt? 符合實際情況嗎? 那種選擇更好一些? |
動力學(xué)或熱力學(xué)計算 |
木蟲 (著名寫手)
木蟲 (正式寫手)
木蟲 (小有名氣)
|
記得先前有位版友也提到過 高溫的時候 lost atoms。 我不是用lammps的,所以具體不清楚什么原因。但是據(jù)經(jīng)驗分析,lost atoms,是不是因為粒子速度過大,在步長之內(nèi)的位移過大,以至于越出了邊界? 溫度的擴(kuò)散容易模擬和控制,但是壓強(qiáng)的擴(kuò)散就麻煩一點(diǎn)了。一般來說,壓強(qiáng)可以通過外加應(yīng)力,或者改變晶格參數(shù)來實現(xiàn)。如果不是恒壓,應(yīng)該選取NPT的系綜,而且體積大小、形狀都允許發(fā)生變化。 建議,先做溫度擴(kuò)散,再做壓力擴(kuò)散,然后再綜合考慮。 |
木蟲 (正式寫手)
木蟲 (小有名氣)
木蟲 (小有名氣)
木蟲 (正式寫手)
木蟲 (小有名氣)
木蟲 (正式寫手)
|
時間步長設(shè)置 0.01fs 300K初始弛豫,時間步長改為 1fs 后開始升溫2000K(整體NVE+部分temp/rescale(速度標(biāo)度))也可以進(jìn)行。但是有個問題--由于整體NVE,導(dǎo)致在部分進(jìn)行速度標(biāo)度是,溫度向周圍傳遞的很快(形成連續(xù)的溫度梯度),希望的溫度梯度最好不是很連續(xù)--中間快速升溫時,兩端溫度上升的越慢越好,當(dāng)隨著時間的推移,溫度真正傳遞到這里時才出現(xiàn)真正的升溫過程。 [ Last edited by xiaowu787 on 2010-9-30 at 09:54 ] |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|