 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 1 | 1/1 | 返回列表 |
| 查看: 1791 | 回復(fù): 0 | |||
cgzhang_gg銅蟲 (著名寫手)
|
[交流]
關(guān)于 run 0 command 重要性
|
|
lammps 里面的run 0 command 具有重要的作用：不妨看下面的例子： ********************************************************************************************************************* # TiO2 Rutile variable a loop 4 log vacan.formation.tect.O.$a units metal atom_style charge boundary p p p lattice custom 4.493 a1 1.0 0.0 0.0 a2 0.0 1.0 0.0 a3 0.0 0.0 0.669& basis 0.0 0.0 0.0 basis 0.5 0.5 0.5 & basis 0.303 0.303 0.0 basis 0.697 0.697 0.0 & basis 0.803 0.197 0.5 basis 0.197 0.803 0.5 region mybox block -4 4 -4 4 -4 4 create_box 2 mybox create_atoms 2 box basis 1 1 basis 2 1 & basis 3 2 basis 4 2 basis 5 2 basis 6 2 # mass after simulation box is defined, before a simulation is run. mass 1 47.867 # Ti mass 2 15.9994 # O set type 1 charge 2.196 # Ti set type 2 charge -1.098 # O # potential pair_style hybrid/overlay table linear 150000 coul/long 11.0 pair_coeff 1 1 table TiO2.pair.table Ti_Ti 11.0 pair_coeff 1 2 table TiO2.pair.table Ti_O 11.0 pair_coeff 2 2 table TiO2.pair.table O_O 11.0 pair_coeff * * coul/long # right kspace_style pppm 1.0e-5  Sample Text variable E equal pe #variable cE equal $E variable N equal atoms #print "@@@@before minimum! cE=${cE}" neighbor 0.3 bin neigh_modify delay 2 every 1 check yes thermo 1 min_style cg minimize 1.0e-20 1.0e-20 10000 100000 #etol ,ftol, max iterations, max number print "@@@@@ 1. atoms and total potential: $N $E" # output pe after a minimum clear next a #jump in.TiO2 muti_point ******************************************************************************************************************************* in.TiO2 例子是進行了四次相同的結(jié)構(gòu)優(yōu)化，對同一個構(gòu)型，得到的正確的結(jié)果應(yīng)該相同的： @@@@@ 1. atoms and total potential: 3072 -40730.760848533624085 @@@@@ 1. atoms and total potential: 3072 -40730.70939078804804 @@@@@ 1. atoms and total potential: 3072 -40730.70939078804804 @@@@@ 1. atoms and total potential: 3072 -40730.70939078804804 但經(jīng)過測試，只有當(dāng)在in.TiO2文件中，在重新設(shè)置構(gòu)型后，進行一次 run 0 command 之后，才可以得到相同的結(jié)果，如上。如果沒有書寫 run 0 command，則得到的結(jié)果不相同，根據(jù)不同的機器，結(jié)果不同。如下： @@@@@ 1. atoms and total potential: 3072 -40730.76084853617067 @@@@@ 1. atoms and total potential: 3072 -45348.092323018747265 @@@@@ 1. atoms and total potential: 3072 -42739.32519361760933 @@@@@ 1. atoms and total potential: 3072 -40900.023453094305296 自己的理解：run 0 command 有一些初始化的作用。當(dāng)不用run 0 command 時pe 可能是不確定的。不知道大家是怎么理解的？ |
回復(fù)此樓» 猜你喜歡
0703化學(xué)調(diào)劑
已經(jīng)有12人回復(fù)
-
0703化學(xué) 305求調(diào)劑
已經(jīng)有4人回復(fù)
-
0703化學(xué)調(diào)劑，求各位老師收留
已經(jīng)有10人回復(fù)
-
271材料工程求調(diào)劑
已經(jīng)有5人回復(fù)
-
281求調(diào)劑（0805）
已經(jīng)有16人回復(fù)
-
304求調(diào)劑
已經(jīng)有6人回復(fù)
-
材料工程專碩調(diào)劑
已經(jīng)有6人回復(fù)
-
一志愿天大材料與化工（085600）總分338
已經(jīng)有4人回復(fù)
-
085700資源與環(huán)境308求調(diào)劑
已經(jīng)有3人回復(fù)
-
求材料調(diào)劑
已經(jīng)有8人回復(fù)
找到一些相關(guān)的精華帖子，希望有用哦~
【求助】vasp安裝編譯錯誤，請大家?guī)兔聪履抢锍鲥e了。
已經(jīng)有6人回復(fù)
- 【原創(chuàng)】64位linux5.4下安裝MS5.0全過程及MS linux原文件lic文件（已搜索，無重復(fù)）
已經(jīng)有43人回復(fù)
點擊這里搜索更多相關(guān)資源科研從小木蟲開始，人人為我，我為人人
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 cgzhang_gg 的主題更新
| 1 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 281求調(diào)劑（0805） +6 | 煙汐憶海 2026-03-16 | 16/800 |
|
|---|---|---|---|---|
|
|
[考研] 一志愿天津大學(xué)化學(xué)工藝專業(yè)（081702）315分求調(diào)劑 +10 | yangfz 2026-03-17 | 10/500 |
|
|
[考研] 一志愿中國海洋大學(xué)，生物學(xué)，301分，求調(diào)劑 +4 | 1孫悟空 2026-03-17 | 4/200 |
|
|
|
[考研] 266求調(diào)劑 +5 | 陽陽哇塞 2026-03-14 | 9/450 |
|
|
|
[考研] 311求調(diào)劑 +6 | 26研0 2026-03-15 | 6/300 |
|
|
|
[考研] 收復(fù)試調(diào)劑生 +4 | 雨后秋荷 2026-03-18 | 4/200 |
|
|
|
[考研] 材料，紡織，生物（0856、0710），化學(xué)招生啦 +3 | Eember. 2026-03-17 | 9/450 |
|
|
|
[考研] 307求調(diào)劑 +3 | 冷笙123 2026-03-17 | 3/150 |
|
|
|
[考研] 268求調(diào)劑 +6 | 簡單點0 2026-03-17 | 6/300 |
|
|
|
[考研] 有沒有道鐵/土木的想調(diào)劑南林，給自己招師弟中～ +3 | TqlXswl 2026-03-16 | 7/350 |
|
|
|
[考研] 材料工程專碩274一志愿211求調(diào)劑 +6 | 薛云鵬 2026-03-15 | 6/300 |
|
|
|
[考研] 285化工學(xué)碩求調(diào)劑（081700） +9 | 柴郡貓_ 2026-03-12 | 9/450 |
|
|
|
[考研] 東南大學(xué)364求調(diào)劑 +5 | JasonYuiui 2026-03-15 | 5/250 |
|
|
|
[考研] 326求調(diào)劑 +3 | mlpqaz03 2026-03-15 | 3/150 |
|
|
|
[考研] 22408總分284求調(diào)劑 +3 | InAspic 2026-03-13 | 3/150 |
|
|
|
[考研] 070305求調(diào)劑 +3 | mlpqaz03 2026-03-14 | 4/200 |
|
|
|
[考研] 288求調(diào)劑 +4 | 奇點0314 2026-03-14 | 4/200 |
|
|
|
[考研] 一志愿哈工大材料324分求調(diào)劑 +5 | 閆旭東 2026-03-14 | 5/250 |
|
|
|
[考研] 求材料調(diào)劑 085600英一數(shù)二總分302 前三科235 精通機器學(xué)習(xí) 一志愿哈工大 +4 | 林yaxin 2026-03-12 | 4/200 |
|
|
|
[考研] 311求調(diào)劑 +3 | 冬十三 2026-03-13 | 3/150 |
|