 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 3/1 | 返回列表 |
| 查看: 6753 | 回復(fù): 29 | |||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | |||
mgqqlwq銅蟲 (著名寫手)
|
[求助]
求教我的二氧化硅體系的弛豫應(yīng)該采用nvt系綜還是npt系綜 已有3人參與
|
||
|
我用VMD建立了SiO2的盒子，LAMMPS里面三個(gè)邊界都選為周期性的，初始化和定義力場之后，利用minimize和npt（或者nvt）系綜對體系進(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ū)別是壓力值差別非常大。請問對體系弛豫的時(shí)候應(yīng)該采用哪一種系綜呢？ 我對minimize命令和系綜的使用概念還不是很清楚，還請各位模擬高手多多指導(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ù)此樓» 猜你喜歡
311求調(diào)劑
已經(jīng)有9人回復(fù)
-
302求調(diào)劑
已經(jīng)有4人回復(fù)
-
被我言中：新模板不強(qiáng)調(diào)格式了，假專家開始管格式了
已經(jīng)有4人回復(fù)
-
303求調(diào)劑
已經(jīng)有5人回復(fù)
-
考研求調(diào)劑
已經(jīng)有4人回復(fù)
-
085601專碩，總分342求調(diào)劑，地區(qū)不限
已經(jīng)有4人回復(fù)
-
085601材料工程專碩求調(diào)劑
已經(jīng)有5人回復(fù)
-
一志愿天津大學(xué)化學(xué)工藝專業(yè)（081702）315分求調(diào)劑
已經(jīng)有6人回復(fù)
-
能源材料化學(xué)課題組招收碩士研究生8-10名
已經(jīng)有6人回復(fù)
-
296求調(diào)劑
已經(jīng)有5人回復(fù)
brucefan
專家顧問 (著名寫手)
-
專家經(jīng)驗(yàn): +218 - 模擬EPI: 4
- 應(yīng)助: 422 (碩士)
- 金幣: 2294.8
- 散金: 1720
- 紅花: 88
- 帖子: 1423
- 在線: 670.2小時(shí)
- 蟲號: 3048965
- 注冊: 2014-03-14
- 性別: GG
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
- 管轄: 分子模擬
mgqqlwq22樓2018-01-08 14:36:58
hongyingwu
至尊木蟲 (文壇精英)
- 應(yīng)助: 148 (高中生)
- 金幣: 17761
- 紅花: 28
- 帖子: 29761
- 在線: 452.5小時(shí)
- 蟲號: 4388826
- 注冊: 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
hongyingwu
至尊木蟲 (文壇精英)
- 應(yīng)助: 148 (高中生)
- 金幣: 17761
- 紅花: 28
- 帖子: 29761
- 在線: 452.5小時(shí)
- 蟲號: 4388826
- 注冊: 2016-01-30
- 專業(yè): 動(dòng)物生理及行為學(xué)
3樓2018-01-05 09:51:24
mgqqlwq
銅蟲 (著名寫手)
- 應(yīng)助: 12 (小學(xué)生)
- 金幣: 12831.1
- 散金: 108
- 紅花: 14
- 帖子: 1282
- 在線: 207.6小時(shí)
- 蟲號: 2436148
- 注冊: 2013-04-25
- 專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和
4樓2018-01-05 11:00:47
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 303求調(diào)劑 +3 | 睿08 2026-03-17 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 能源材料化學(xué)課題組招收碩士研究生8-10名 +3 | 脫穎而出 2026-03-16 | 6/300 |
|
|
|
[考研] 268求調(diào)劑 +7 | 好運(yùn)連綿不絕 2026-03-12 | 8/400 |
|
|
|
[考研] 材料專碩306英一數(shù)二 +7 | z1z2z3879 2026-03-16 | 9/450 |
|
|
|
[考研] 326求調(diào)劑 +5 | 上岸的小葡 2026-03-15 | 6/300 |
|
|
|
[考研] 考研化學(xué)學(xué)碩調(diào)劑，一志愿985 +4 | 張vvvv 2026-03-15 | 6/300 |
|
|
|
[考研] 化學(xué)工程321分求調(diào)劑 +11 | 大米飯！ 2026-03-15 | 14/700 |
|
|
|
[考研] 274求調(diào)劑0856材料化工 +13 | z2839474511 2026-03-11 | 14/700 |
|
|
|
[考研] 材料與化工專碩調(diào)劑 +5 | heming3743 2026-03-16 | 5/250 |
|
|
|
[教師之家] 焦慮 +7 | 水冰月月野兔 2026-03-13 | 9/450 |
|
|
|
[考研] 材料與化工 323 英一+數(shù)二+物化，一志愿：哈工大 本人本科雙一流 +4 | 自由的_飛翔 2026-03-13 | 5/250 |
|
|
|
[考研] 285 求調(diào)劑 資源與環(huán)境 一志愿北京化工大學(xué) +3 | 未名考生 2026-03-10 | 3/150 |
|
|
|
[考研] 304求調(diào)劑 +6 | Mochaaaa 2026-03-12 | 7/350 |
|
|
|
[考研] 308求調(diào)劑 +5 | 是Lupa啊 2026-03-11 | 5/250 |
|
|
|
[考研] 304求調(diào)劑 +7 | 7712b 2026-03-13 | 7/350 |
|
|
|
[考研] 329求調(diào)劑 +3 | miaodesi 2026-03-12 | 4/200 |
|
|
|
[考研] 290求調(diào)劑 +7 | ADT 2026-03-12 | 7/350 |
|
|
|
[考研] 材料301分求調(diào)劑 +5 | Liyouyumairs 2026-03-12 | 5/250 |
|
|
|
[考研] 一志愿山大07化學(xué) 332分 四六級已過 本科山東雙非 求調(diào)劑！ +3 | 不想理你 2026-03-12 | 3/150 |
|
|
|
[考研] 求調(diào)劑 資源與環(huán)境 285 +3 | 未名考生 2026-03-10 | 3/150 |
|