 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 查看: 6763 | 回復(fù): 29 | |||
mgqqlwq銅蟲 (著名寫手)
|
[求助]
求教我的二氧化硅體系的弛豫應(yīng)該采用nvt系綜還是npt系綜 已有3人參與
|
||
|
我用VMD建立了SiO2的盒子，LAMMPS里面三個(gè)邊界都選為周期性的，初始化和定義力場(chǎng)之后，利用minimize和npt（或者nvt）系綜對(duì)體系進(jìn)行弛豫，in文件如下，區(qū)別是一種采用nvt系綜而另一種采用npt系綜，根據(jù)結(jié)果來看兩種弛豫的主要區(qū)別在于100000步之后的壓力區(qū)別比較大。本人小白有兩個(gè)問題求教一下各位模擬的大神和壇友。 1. minimize應(yīng)該是在0K下的結(jié)構(gòu)最優(yōu)化嗎？velocity create之后體系原子的溫度已經(jīng)是300K了，這時(shí)候做最小化是否合理？ 2. 采用nvt系綜100000步之后的thermo信息是： Step Temp E_pair E_mol TotEng Press 100000 299.99624 -1786635.5 0 -1776212.1 38516.104 而采用npt系綜的thermo信息是： Step Temp E_pair E_mol TotEng Press Volume 100000 300.35095 -1789725.2 0 -1779289.5 -126.32433 4191230.7 從結(jié)果看主要的區(qū)別是壓力值差別非常大。請(qǐng)問對(duì)體系弛豫的時(shí)候應(yīng)該采用哪一種系綜呢？ 我對(duì)minimize命令和系綜的使用概念還不是很清楚，還請(qǐng)各位模擬高手多多指導(dǎo)和幫助。 我的in文件如下： echo both variable T equal 300 units metal dimension 3 boundary p p p atom_style full read_data silicondioxide.data pair_style tersoff pair_coeff * * SiO.tersoff O Si timestep 0.001 thermo 10 velocity all create $T 38168 rot yes mom yes dist gaussian dump 1 all custom 1000 dump.silicondioxide.minimize id type x y z min_style cg minimize 0 0 100000 100000 reset_timestep 0 undump 1 dump 2 all custom 5000 dump.silicondioxide.nvt id type x y z fix nvt_relax all nvt temp $T $T 0.1 drag 0.2 #fix npt_relax all npt temp $T $T 0.1 iso 0 0 1 drag 0.2 #兩種方法的區(qū)別是另一種采用NPT系綜 run 100000 unfix nvt_relax |
回復(fù)此樓» 猜你喜歡
【考研調(diào)劑】化學(xué)專業(yè) 281分，一志愿四川大學(xué)，誠心求調(diào)劑
已經(jīng)有11人回復(fù)
-
一志愿深大，0703化學(xué)，總分302，求調(diào)劑
已經(jīng)有4人回復(fù)
-
材料學(xué)碩333求調(diào)劑
已經(jīng)有3人回復(fù)
-
材料求調(diào)劑
已經(jīng)有3人回復(fù)
-
333求調(diào)劑
已經(jīng)有5人回復(fù)
-
求助
已經(jīng)有6人回復(fù)
-
278求調(diào)劑
已經(jīng)有9人回復(fù)
-
297求調(diào)劑
已經(jīng)有12人回復(fù)
-
0703化學(xué)297求調(diào)劑
已經(jīng)有3人回復(fù)
-
材料學(xué)碩301分求調(diào)劑
已經(jīng)有6人回復(fù)
brucefan
專家顧問 (著名寫手)
-
專家經(jīng)驗(yàn): +218 - 模擬EPI: 4
- 應(yīng)助: 422 (碩士)
- 金幣: 2294.8
- 散金: 1720
- 紅花: 88
- 帖子: 1423
- 在線: 670.2小時(shí)
- 蟲號(hào): 3048965
- 注冊(cè): 2014-03-14
- 性別: GG
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
- 管轄: 分子模擬
7樓2018-01-06 15:43:24
19樓2018-01-08 03:35:27
hongyingwu
至尊木蟲 (文壇精英)
- 應(yīng)助: 148 (高中生)
- 金幣: 17761
- 紅花: 28
- 帖子: 29761
- 在線: 452.5小時(shí)
- 蟲號(hào): 4388826
- 注冊(cè): 2016-01-30
- 專業(yè): 動(dòng)物生理及行為學(xué)
3樓2018-01-05 09:51:24
hongyingwu
至尊木蟲 (文壇精英)
- 應(yīng)助: 148 (高中生)
- 金幣: 17761
- 紅花: 28
- 帖子: 29761
- 在線: 452.5小時(shí)
- 蟲號(hào): 4388826
- 注冊(cè): 2016-01-30
- 專業(yè): 動(dòng)物生理及行為學(xué)
【答案】應(yīng)助回帖
★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
mgqqlwq: 金幣+5 2018-01-06 16:20:45
感謝參與，應(yīng)助指數(shù) +1
mgqqlwq: 金幣+5 2018-01-06 16:20:45
2樓2018-01-05 09:49:44
mgqqlwq
銅蟲 (著名寫手)
- 應(yīng)助: 12 (小學(xué)生)
- 金幣: 12831.1
- 散金: 108
- 紅花: 14
- 帖子: 1282
- 在線: 207.6小時(shí)
- 蟲號(hào): 2436148
- 注冊(cè): 2013-04-25
- 專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和
4樓2018-01-05 11:00:47
mgqqlwq
銅蟲 (著名寫手)
- 應(yīng)助: 12 (小學(xué)生)
- 金幣: 12831.1
- 散金: 108
- 紅花: 14
- 帖子: 1282
- 在線: 207.6小時(shí)
- 蟲號(hào): 2436148
- 注冊(cè): 2013-04-25
- 專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和
5樓2018-01-05 23:00:07
6樓2018-01-06 05:47:46
mgqqlwq
銅蟲 (著名寫手)
- 應(yīng)助: 12 (小學(xué)生)
- 金幣: 12831.1
- 散金: 108
- 紅花: 14
- 帖子: 1282
- 在線: 207.6小時(shí)
- 蟲號(hào): 2436148
- 注冊(cè): 2013-04-25
- 專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和
8樓2018-01-06 16:26:44
mgqqlwq
銅蟲 (著名寫手)
- 應(yīng)助: 12 (小學(xué)生)
- 金幣: 12831.1
- 散金: 108
- 紅花: 14
- 帖子: 1282
- 在線: 207.6小時(shí)
- 蟲號(hào): 2436148
- 注冊(cè): 2013-04-25
- 專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和
9樓2018-01-06 16:54:26
brucefan
專家顧問 (著名寫手)
-
專家經(jīng)驗(yàn): +218 - 模擬EPI: 4
- 應(yīng)助: 422 (碩士)
- 金幣: 2294.8
- 散金: 1720
- 紅花: 88
- 帖子: 1423
- 在線: 670.2小時(shí)
- 蟲號(hào): 3048965
- 注冊(cè): 2014-03-14
- 性別: GG
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
- 管轄: 分子模擬
10樓2018-01-06 17:49:10
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 mgqqlwq 的主題更新
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 求調(diào)劑 +3 | 要好好無聊 2026-03-21 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 299求調(diào)劑 +4 | 某某某某位 2026-03-21 | 4/200 |
|
|
|
[考研] 346求調(diào)劑[0856] +4 | WayneLim327 2026-03-16 | 7/350 |
|
|
|
[考研] 一志愿山大07化學(xué) 332分 四六級(jí)已過 本科山東雙非 求調(diào)劑！ +3 | 不想理你 2026-03-16 | 3/150 |
|
|
|
[考研] 310求調(diào)劑 +3 | baibai1314 2026-03-16 | 3/150 |
|
|
|
[考研] 一志愿華中科技大學(xué)，080502，354分求調(diào)劑 +5 | 守候夕陽CF 2026-03-18 | 5/250 |
|
|
|
[考研] 一志愿重慶大學(xué)085700資源與環(huán)境專碩，總分308求調(diào)劑 +3 | 墨墨漠 2026-03-18 | 3/150 |
|
|
|
[考研] 295求調(diào)劑 +4 | 一志愿京區(qū)211 2026-03-18 | 6/300 |
|
|
|
[考研] 324求調(diào)劑 +5 | lucky呀呀呀鴨 2026-03-20 | 5/250 |
|
|
|
[考研] 一志愿 西北大學(xué) ，070300化學(xué)學(xué)碩，總分287，雙非一本，求調(diào)劑。 +4 | 晨昏線與星海 2026-03-19 | 4/200 |
|
|
|
[考研] 求調(diào)劑，一志愿:南京航空航天大學(xué)大學(xué) ，080500材料科學(xué)與工程學(xué)碩，總分289分 +4 | @taotao 2026-03-19 | 4/200 |
|
|
|
[考研] 0817 化學(xué)工程 299分求調(diào)劑 有科研經(jīng)歷 有二區(qū)文章 +22 | rare12345 2026-03-18 | 22/1100 |
|
|
|
[考研] 298-一志愿中國農(nóng)業(yè)大學(xué)-求調(diào)劑 +9 | 手機(jī)用戶 2026-03-17 | 9/450 |
|
|
|
[考研] 材料工程專碩274一志愿211求調(diào)劑 +6 | 薛云鵬 2026-03-15 | 6/300 |
|
|
|
[考研] 考研調(diào)劑 +3 | 淇ya_~ 2026-03-17 | 5/250 |
|
|
|
[考研] 機(jī)械專碩325，尋找調(diào)劑院校 +3 | y9999 2026-03-15 | 5/250 |
|
|
|
[考研] 333求調(diào)劑 +3 | 文思客 2026-03-16 | 7/350 |
|
|
|
[考研] 中科院材料273求調(diào)劑 +4 | yzydy 2026-03-15 | 4/200 |
|
|
|
[考研] 070303 總分349求調(diào)劑 +3 | LJY9966 2026-03-15 | 5/250 |
|
|
|
[考研] 本科南京大學(xué)一志愿川大藥學(xué)327 +3 | 麥田耕者 2026-03-14 | 3/150 |
|