 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 查看: 6761 | 回復: 29 | |||
mgqqlwq銅蟲 (著名寫手)
|
[求助]
求教我的二氧化硅體系的弛豫應該采用nvt系綜還是npt系綜 已有3人參與
|
||
|
我用VMD建立了SiO2的盒子，LAMMPS里面三個邊界都選為周期性的，初始化和定義力場之后，利用minimize和npt（或者nvt）系綜對體系進行弛豫，in文件如下，區(qū)別是一種采用nvt系綜而另一種采用npt系綜，根據(jù)結果來看兩種弛豫的主要區(qū)別在于100000步之后的壓力區(qū)別比較大。本人小白有兩個問題求教一下各位模擬的大神和壇友。 1. minimize應該是在0K下的結構最優(yōu)化嗎？velocity create之后體系原子的溫度已經(jīng)是300K了，這時候做最小化是否合理？ 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 從結果看主要的區(qū)別是壓力值差別非常大。請問對體系弛豫的時候應該采用哪一種系綜呢？ 我對minimize命令和系綜的使用概念還不是很清楚，還請各位模擬高手多多指導和幫助。 我的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 |
回復此樓» 猜你喜歡
346求調劑[0856]
已經(jīng)有7人回復
-
一志愿山大07化學 332分 四六級已過 本科山東雙非 求調劑！
已經(jīng)有3人回復
-
310求調劑
已經(jīng)有3人回復
-
機械專碩299求調劑至材料
已經(jīng)有4人回復
-
070300化學319求調劑
已經(jīng)有7人回復
-
08工科 320總分 求調劑
已經(jīng)有6人回復
-
一志愿天津大學化學工藝專業(yè)（081702）315分求調劑
已經(jīng)有12人回復
-
307求調劑
已經(jīng)有3人回復
-
265求調劑
已經(jīng)有3人回復
-
301求調劑
已經(jīng)有10人回復
brucefan
專家顧問 (著名寫手)
-
專家經(jīng)驗: +218 - 模擬EPI: 4
- 應助: 422 (碩士)
- 金幣: 2294.8
- 散金: 1720
- 紅花: 88
- 帖子: 1423
- 在線: 670.2小時
- 蟲號: 3048965
- 注冊: 2014-03-14
- 性別: GG
- 專業(yè): 凝聚態(tài)物性 II ：電子結構
- 管轄: 分子模擬
7樓2018-01-06 15:43:24
19樓2018-01-08 03:35:27
hongyingwu
至尊木蟲 (文壇精英)
- 應助: 148 (高中生)
- 金幣: 17761
- 紅花: 28
- 帖子: 29761
- 在線: 452.5小時
- 蟲號: 4388826
- 注冊: 2016-01-30
- 專業(yè): 動物生理及行為學
3樓2018-01-05 09:51:24
hongyingwu
至尊木蟲 (文壇精英)
- 應助: 148 (高中生)
- 金幣: 17761
- 紅花: 28
- 帖子: 29761
- 在線: 452.5小時
- 蟲號: 4388826
- 注冊: 2016-01-30
- 專業(yè): 動物生理及行為學
【答案】應助回帖
★ ★ ★ ★ ★
感謝參與，應助指數(shù) +1
mgqqlwq: 金幣+5 2018-01-06 16:20:45
感謝參與，應助指數(shù) +1
mgqqlwq: 金幣+5 2018-01-06 16:20:45
2樓2018-01-05 09:49:44
mgqqlwq
銅蟲 (著名寫手)
- 應助: 12 (小學生)
- 金幣: 12831.1
- 散金: 108
- 紅花: 14
- 帖子: 1282
- 在線: 207.6小時
- 蟲號: 2436148
- 注冊: 2013-04-25
- 專業(yè): 凝聚態(tài)物性I:結構、力學和
4樓2018-01-05 11:00:47
mgqqlwq
銅蟲 (著名寫手)
- 應助: 12 (小學生)
- 金幣: 12831.1
- 散金: 108
- 紅花: 14
- 帖子: 1282
- 在線: 207.6小時
- 蟲號: 2436148
- 注冊: 2013-04-25
- 專業(yè): 凝聚態(tài)物性I:結構、力學和
5樓2018-01-05 23:00:07
6樓2018-01-06 05:47:46
mgqqlwq
銅蟲 (著名寫手)
- 應助: 12 (小學生)
- 金幣: 12831.1
- 散金: 108
- 紅花: 14
- 帖子: 1282
- 在線: 207.6小時
- 蟲號: 2436148
- 注冊: 2013-04-25
- 專業(yè): 凝聚態(tài)物性I:結構、力學和
8樓2018-01-06 16:26:44
mgqqlwq
銅蟲 (著名寫手)
- 應助: 12 (小學生)
- 金幣: 12831.1
- 散金: 108
- 紅花: 14
- 帖子: 1282
- 在線: 207.6小時
- 蟲號: 2436148
- 注冊: 2013-04-25
- 專業(yè): 凝聚態(tài)物性I:結構、力學和
9樓2018-01-06 16:54:26
brucefan
專家顧問 (著名寫手)
-
專家經(jīng)驗: +218 - 模擬EPI: 4
- 應助: 422 (碩士)
- 金幣: 2294.8
- 散金: 1720
- 紅花: 88
- 帖子: 1423
- 在線: 670.2小時
- 蟲號: 3048965
- 注冊: 2014-03-14
- 性別: GG
- 專業(yè): 凝聚態(tài)物性 II ：電子結構
- 管轄: 分子模擬
10樓2018-01-06 17:49:10
相關版塊跳轉
我要訂閱樓主 mgqqlwq 的主題更新
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 310求調劑 +3 | baibai1314 2026-03-16 | 3/150 |
|
|---|---|---|---|---|
|
[考研]
 271材料工程求調劑
+8
|
.6lL 2026-03-18 | 8/400 |
|
|
|
[考研] 一志愿西南交大，求調劑 +5 | 材化逐夢人 2026-03-18 | 5/250 |
|
|
|
[考研] 294求調劑材料與化工專碩 +15 | 陌の森林 2026-03-18 | 15/750 |
|
|
|
[考研] 350求調劑 +5 | weudhdk 2026-03-19 | 5/250 |
|
|
|
[考研] A區(qū)線材料學調劑 +5 | 周周無極 2026-03-20 | 5/250 |
|
|
|
[考研] 081700化工學碩調劑 +3 | 【1】 2026-03-16 | 3/150 |
|
|
|
[考研] 0703化學 305求調劑 +4 | FY_yy 2026-03-14 | 4/200 |
|
|
|
[考研] 材料，紡織，生物（0856、0710），化學招生啦 +3 | Eember. 2026-03-17 | 9/450 |
|
|
|
[考研] 0703化學336分求調劑 +6 | zbzihdhd 2026-03-15 | 7/350 |
|
|
|
[考博] 26博士申請 +3 | 1042136743 2026-03-17 | 3/150 |
|
|
|
[考研] 326求調劑 +5 | 上岸的小葡 2026-03-15 | 6/300 |
|
|
|
[考研] 材料專碩326求調劑 +6 | 墨煜姒莘 2026-03-15 | 7/350 |
|
|
|
[考研] 有沒有道鐵/土木的想調劑南林，給自己招師弟中～ +3 | TqlXswl 2026-03-16 | 7/350 |
|
|
|
[考研] 302求調劑 +4 | 小賈同學123 2026-03-15 | 8/400 |
|
|
|
[考研] 304求調劑 +4 | ahbd 2026-03-14 | 4/200 |
|
|
|
[考研] 070305求調劑 +3 | mlpqaz03 2026-03-14 | 4/200 |
|
|
|
[考研] 085601材料工程315分求調劑 +3 | yang_0104 2026-03-15 | 3/150 |
|
|
|
[考研] 080500，材料學碩302分求調劑學校 +4 | 初識可樂 2026-03-14 | 5/250 |
|
|
|
[考研] 本科南京大學一志愿川大藥學327 +3 | 麥田耕者 2026-03-14 | 3/150 |
|