 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 2/1 | 返回列表 |
| 查看: 6755 | 回復(fù): 29 | ||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | ||
mgqqlwq銅蟲 (著名寫手)
|
[求助]
求教我的二氧化硅體系的弛豫應(yīng)該采用nvt系綜還是npt系綜 已有3人參與
|
|
|
我用VMD建立了SiO2的盒子，LAMMPS里面三個(gè)邊界都選為周期性的，初始化和定義力場之后，利用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)劑
已經(jīng)有7人回復(fù)
-
一志愿天大材料與化工（085600）總分338
已經(jīng)有3人回復(fù)
-
085600材料與化工
已經(jīng)有5人回復(fù)
-
085600材料與化工調(diào)劑 324分
已經(jīng)有6人回復(fù)
-
286求調(diào)劑
已經(jīng)有9人回復(fù)
-
焦慮
已經(jīng)有12人回復(fù)
-
344求調(diào)劑
已經(jīng)有6人回復(fù)
-
266求調(diào)劑
已經(jīng)有9人回復(fù)
-
化學(xué)工程321分求調(diào)劑
已經(jīng)有18人回復(fù)
-
314求調(diào)劑
已經(jīng)有8人回復(fù)
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é)和
20樓2018-01-08 11:24:07
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
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
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
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 266求調(diào)劑 +5 | 陽陽哇塞 2026-03-14 | 9/450 |
|
|---|---|---|---|---|
|
|
[考研] 085601專碩，總分342求調(diào)劑，地區(qū)不限 +5 | share_joy 2026-03-16 | 5/250 |
|
|
|
[考研] 311求調(diào)劑 +6 | 26研0 2026-03-15 | 6/300 |
|
|
|
[考研] 0817 化學(xué)工程 299分求調(diào)劑 有科研經(jīng)歷 有二區(qū)文章 +7 | rare12345 2026-03-18 | 7/350 |
|
|
|
[考研] 一志愿天津大學(xué)化學(xué)工藝專業(yè)（081702）315分求調(diào)劑 +9 | yangfz 2026-03-17 | 9/450 |
|
|
[考研] 0703化學(xué)調(diào)劑 +4 | pupcoco 2026-03-17 | 7/350 |
|
|
|
[考研] 一志愿中國海洋大學(xué)，生物學(xué)，301分，求調(diào)劑 +3 | 1孫悟空 2026-03-17 | 3/150 |
|
|
|
[考研] 生物學(xué)071000 329分求調(diào)劑 +3 | 我愛生物生物愛?/a> 2026-03-17 | 3/150 |
|
|
|
[基金申請(qǐng)] 被我言中：新模板不強(qiáng)調(diào)格式了，假專家開始管格式了 +4 | beefly 2026-03-14 | 4/200 |
|
|
|
[考研] 283求調(diào)劑 +3 | 聽風(fēng)就是雨； 2026-03-16 | 3/150 |
|
|
|
[考研]
 11408 一志愿西電，277分求調(diào)劑
+3
|
zhouzhen654 2026-03-16 | 3/150 |
|
|
|
[考研] 0703一志愿211 285分求調(diào)劑 +5 | ly3471z 2026-03-13 | 5/250 |
|
|
|
[考研] 一志愿華中師范071000，325求調(diào)劑 +6 | RuitingC 2026-03-12 | 6/300 |
|
|
|
[考研] 0856求調(diào)劑 +3 | 劉夢(mèng)微 2026-03-15 | 3/150 |
|
|
|
[考博] 東華理工大學(xué)化材專業(yè)26屆碩士博士申請(qǐng) +6 | zlingli 2026-03-13 | 6/300 |
|
|
|
[考研] 070305求調(diào)劑 +3 | mlpqaz03 2026-03-14 | 4/200 |
|
|
|
[考研] 297求調(diào)劑 +4 | 學(xué)海漂泊 2026-03-13 | 4/200 |
|
|
|
[考研] ［0860］321分求調(diào)劑，ab區(qū)皆可 +4 | 寶貴熱 2026-03-13 | 4/200 |
|
|
|
[考研] 315求調(diào)劑 +9 | 小羊小羊_ 2026-03-11 | 10/500 |
|
|
|
[考研] 295求調(diào)劑 +3 | 小匕仔汁 2026-03-12 | 3/150 |
|