| 查看: 2223 | 回復(fù): 18 | ||||
| 本帖產(chǎn)生 1 個 模擬EPI ,點擊這里進行查看 | ||||
xiaowu787木蟲 (正式寫手)
|
[交流]
【討論】高溫高壓的體系需要注意的環(huán)節(jié) 已有4人參與
|
|||
|
希望LAMMPS模擬的模型中心產(chǎn)生一個高溫高壓(有時間壽命)的局部區(qū),然后觀測溫度和壓力的向周圍擴展的情況。可是溫度高壓強高時會出現(xiàn)lost atoms,現(xiàn)在 第一個疑問是:周期邊界條件下的image atom怎么不起作用了--周期性的邊界條件不是可以保證粒子數(shù)不變? 第二: 局部升溫的選擇: 1. fix all nve+fix 局部 temp/rescale,可是壓強無法控制 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)的位移過大,以至于越出了邊界? 溫度的擴散容易模擬和控制,但是壓強的擴散就麻煩一點了。一般來說,壓強可以通過外加應(yīng)力,或者改變晶格參數(shù)來實現(xiàn)。如果不是恒壓,應(yīng)該選取NPT的系綜,而且體積大小、形狀都允許發(fā)生變化。 建議,先做溫度擴散,再做壓力擴散,然后再綜合考慮。 |
木蟲 (正式寫手)
木蟲 (小有名氣)
木蟲 (小有名氣)
木蟲 (正式寫手)
木蟲 (小有名氣)
木蟲 (正式寫手)
|
時間步長設(shè)置 0.01fs 300K初始弛豫,時間步長改為 1fs 后開始升溫2000K(整體NVE+部分temp/rescale(速度標(biāo)度))也可以進行。但是有個問題--由于整體NVE,導(dǎo)致在部分進行速度標(biāo)度是,溫度向周圍傳遞的很快(形成連續(xù)的溫度梯度),希望的溫度梯度最好不是很連續(xù)--中間快速升溫時,兩端溫度上升的越慢越好,當(dāng)隨著時間的推移,溫度真正傳遞到這里時才出現(xiàn)真正的升溫過程。 [ Last edited by xiaowu787 on 2010-9-30 at 09:54 ] |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 311求調(diào)劑 +13 | 冬十三 2026-03-15 | 14/700 |
|
|---|---|---|---|---|
|
[考研] 0805 316求調(diào)劑 +3 | 大雪深藏 2026-03-18 | 3/150 |
|
|
[考研] 一志愿南大,0703化學(xué),分?jǐn)?shù)336,求調(diào)劑 +3 | 收到VS 2026-03-21 | 3/150 |
|
|
[考研] 278求調(diào)劑 +9 | 煙火先于春 2026-03-17 | 9/450 |
|
|
[考研] 336求調(diào)劑 +5 | rmc8866 2026-03-21 | 5/250 |
|
|
[考研] 306求調(diào)劑 +4 | chuanzhu川燭 2026-03-18 | 4/200 |
|
|
[考研] 265求調(diào)劑 +9 | 梁梁校校 2026-03-17 | 9/450 |
|
|
[考研] 求調(diào)劑 +3 | Ma_xt 2026-03-17 | 3/150 |
|
|
[考研] 296求調(diào)劑 +6 | www_q 2026-03-18 | 10/500 |
|
|
[考研] 308求調(diào)劑 +3 | 阿姐阿姐家啊 2026-03-18 | 3/150 |
|
|
[考研] 317求調(diào)劑 +5 | 申子申申 2026-03-19 | 9/450 |
|
|
[考研] 287求調(diào)劑 +7 | 晨昏線與星海 2026-03-19 | 8/400 |
|
|
[考研] 藥學(xué)383 求調(diào)劑 +3 | 藥學(xué)chy 2026-03-15 | 5/250 |
|
|
[考研] 材料學(xué)求調(diào)劑 +4 | Stella_Yao 2026-03-20 | 4/200 |
|
|
[考研] 求調(diào)劑 +3 | eation27 2026-03-20 | 3/150 |
|
|
[考研] 招收調(diào)劑碩士 +4 | lidianxing 2026-03-19 | 12/600 |
|
|
[考研] 一志愿中國海洋大學(xué),生物學(xué),301分,求調(diào)劑 +5 | 1孫悟空 2026-03-17 | 6/300 |
|
|
[考研] 考研求調(diào)劑 +3 | 橘頌. 2026-03-17 | 4/200 |
|
|
[考研] 材料專碩326求調(diào)劑 +6 | 墨煜姒莘 2026-03-15 | 7/350 |
|
|
[考研] 機械專碩325,尋找調(diào)劑院校 +3 | y9999 2026-03-15 | 5/250 |
|