 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 4 | 1/1 | 返回列表 |
| 查看: 1711 | 回復(fù): 3 | |||
/name-?=p/銅蟲 (初入文壇)
|
[交流]
lammps在npt平衡后換成nve溫度一直升高是怎么回事？ 已有2人參與
|
|
各位大神請教一個問題 用lammps模擬計算3C-SiC熱導(dǎo)率,時間步長1fs 在1000K下，先在NVT下模擬10萬步，使系統(tǒng)達(dá)到平衡然后換上NVE系綜后在運行125萬步。但是當(dāng)換上NVE系綜不久后，溫度就開始升高，到結(jié)束時溫度都達(dá)到幾千甚至上萬K，這么高溫度碳化硅固體早就沒了。125萬步不是設(shè)定步數(shù)，125萬步時原子丟失了，在低溫時 溫度也會升高 想問這是哪里出問題了 這是in文件 晶格常數(shù)設(shè)置的讀坐標(biāo)文件 盒子邊長大小是十個晶格常數(shù)也就是一共8000個原子 variable T equal 1000 # 溫度為 400 K #variable A equal 4.348 # 晶格常數(shù)大概為 6.35A variable DT equal 0.001 # 積分步長為 0.3 fs # 初始化坐標(biāo)、勢函 數(shù)、速度等數(shù)據(jù) boundary p p p read_data SiC-beta.txt pair_style tersoff pair_coeff * * SiC.tersoff C Si min_style cg #minimize 1.0e-9 1.0e-9 10000 30000 velocity all create ${T} 123456 dist gaussian # 根據(jù)溫度初始化速度 # 在NVT系綜下平衡體系 fix NVT all nvt temp ${T} ${T} 0.1 # Nose-Hoover 熱浴 timestep ${DT} # 設(shè)置積分步長 thermo_style custom step temp press # 觀察溫度和壓強的變化 thermo 1000 # 100 步輸出一次 # 根據(jù)以上設(shè)置跑若干步，目的是使系統(tǒng)達(dá)到熱力學(xué)平衡態(tài) run 100000 # 跑 10000 步 # 在產(chǎn)出階段將用NVE系綜 unfix NVT # 撤掉之前的 NVT 系綜 fix NVE all nve # 換上 NVE 系綜 # 在產(chǎn)出階段計算熱流 compute KE all ke/atom # 計算單原子動能 compute PE all pe/atom # 計算單原子勢能 compute V all stress/atom NULL virial # 計算單原子位力 compute J all heat/flux KE PE V # 計算單原子熱流 # 在產(chǎn)出階段計算熱流自關(guān)聯(lián)函數(shù) <J(0)J(t)> variable Ns equal 100 # 取樣間隔 variable Nc equal 100 # 關(guān)聯(lián)數(shù)據(jù)量（最大關(guān)聯(lián)時間為 Ns*Nc*DT） variable Np equal ${Ns}*${Nc}*10 # 產(chǎn)出步數(shù)（這是一個好的選擇） fix HAC all ave/correlate ${Ns} ${Nc} ${Np} c_J[1] c_J[2] c_J[3] & type auto file hac.txt #�。。。。。。。。。。。。。�！ variable M equal 14 # 熱流自關(guān)聯(lián)函數(shù)的個數(shù) variable N equal ${M}*${Np} # 產(chǎn)出步數(shù) # 跑產(chǎn)出階段 |
回復(fù)此樓» 猜你喜歡
A區(qū)一本交叉課題組，低分調(diào)劑，招收機械電子信息通信等交叉方向
已經(jīng)有15人回復(fù)
-
材料學(xué)碩318求調(diào)劑
已經(jīng)有18人回復(fù)
-
江西師范大學(xué)王昆-盧章輝課題組招化學(xué)方向碩士研究生兩名
已經(jīng)有3人回復(fù)
-
沒上岸的看過來
已經(jīng)有5人回復(fù)
-
290求調(diào)劑
已經(jīng)有5人回復(fù)
-
材料與化工，291，求調(diào)劑
已經(jīng)有4人回復(fù)
-
301求調(diào)劑
已經(jīng)有4人回復(fù)
-
一志愿東北大學(xué)080500，英一數(shù)二343分求調(diào)劑
已經(jīng)有5人回復(fù)
-
274環(huán)境工程求調(diào)劑
已經(jīng)有6人回復(fù)
-
304求調(diào)劑
已經(jīng)有3人回復(fù)
/name-?=p/
銅蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 186.9
- 紅花: 1
- 帖子: 49
- 在線: 9.4小時
- 蟲號: 7570339
- 注冊: 2017-12-05
- 專業(yè): 物理學(xué)I
2樓2020-05-12 18:06:31
lqpzh 
新蟲 (小有名氣)
- 應(yīng)助: 4 (幼兒園)
- 金幣: 111.6
- 散金: 5
- 帖子: 90
- 在線: 5.8小時
- 蟲號: 9150633
- 注冊: 2018-06-21
- 專業(yè): 金屬材料的合金相、相變及
3樓2020-05-13 07:15:26
4樓2020-05-18 08:52:22
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 /name-?=p/ 的主題更新
| 4 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 304求調(diào)劑 +3 | 曼殊2266 2026-03-04 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 085701環(huán)境工程 求調(diào)劑 +7 | xiiiia 2026-03-04 | 7/350 |
|
|
|
[考研] 267調(diào)劑求助 +5 | 聰少OZ 2026-03-04 | 5/250 |
|
|
|
[考研] 化學(xué) 0703求調(diào)劑 總分293 一志愿211 +9 | 土土小蟲 2026-03-03 | 11/550 |
|
|
|
[基金申請] 沒有青基直接申請面上，感覺自己瘋了 +6 | kevin63t 2026-03-02 | 12/600 |
|
|
|
[考研] 成績276，專業(yè)代碼0856求調(diào)劑 +10 | 小陳朵 2026-03-03 | 10/500 |
|
|
|
[考研] 學(xué)碩材料275調(diào)劑 +9 | 路三三 2026-03-03 | 9/450 |
|
|
|
[考研] 292求調(diào)劑 +9 | yhk_819 2026-02-28 | 9/450 |
|
|
|
[考博] 26申博-目前4篇SCI一作 +3 | chen_2024 2026-03-02 | 3/150 |
|
|
|
[考研] 324求調(diào)劑 +4 | wxz2 2026-03-03 | 5/250 |
|
|
|
[考研]
 高分子化學(xué)與物理調(diào)劑
+6
|
好好好1233 2026-02-28 | 16/800 |
|
|
|
[考研] 材料284求調(diào)劑，一志愿鄭州大學(xué)英一數(shù)二專碩 +15 | 想上岸的土撥鼠 2026-02-28 | 15/750 |
|
|
|
[考研] 一志愿山東大學(xué)材料與化工325求調(diào)劑 +5 | 半截的詩0927 2026-03-02 | 5/250 |
|
|
|
[基金申請] 此成果不能導(dǎo)入原因：元數(shù)據(jù)必填信息不完整，可 進(jìn)行補充。 +4 | Kittylucky 2026-03-02 | 5/250 |
|
|
|
[考研] 281求調(diào)劑 +5 | 2026計算機_誠心 2026-03-01 | 8/400 |
|
|
|
[考研] 284求調(diào)劑 +10 | 天下熯 2026-02-28 | 11/550 |
|
|
|
[考研] 274求調(diào)劑 +3 | cgyzqwn 2026-03-01 | 7/350 |
|
|
|
[考研] 材料學(xué)調(diào)劑 +10 | 提神豆沙包 2026-02-28 | 12/600 |
|
|
|
[考研] 290求調(diào)劑 +9 | 材料專碩調(diào)劑； 2026-02-28 | 11/550 |
|
|
|
[考研] 0856材料求調(diào)劑 +4 | 麻辣魷魚 2026-02-28 | 4/200 |
|