 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 查看: 3043 | 回復(fù): 16 | ||
[求助]
lammps 簡(jiǎn)單模型弛豫時(shí)周期邊界出問題。 已有2人參與
|
|
問題：利用lammps建立一個(gè)矩形，x和y方向?yàn)橹芷�，z方向?yàn)榉侵芷冢?p p m) ,利用NVT弛豫（300K），發(fā)現(xiàn)周期邊界處的原子混亂; 如果把(p p m) 改為(p m m)，發(fā)現(xiàn)還是周期邊界方向出現(xiàn)混亂，如下圖。不知什么原因？謝謝！ lammps腳本 #initialize units metal boundary p p m atom_style atomic #============================================================================ # difinition of variable #============================================================================ variable T0 equal 300 #K variable dis equal 15 variable edge equal 20 variable h equal $r+3+${dis} variable subup equal 0 variable subdw equal -20 variable fixedline equal ${subdw}*0.75 variable bxl equal -${r}-${edge} variable bxr equal ${r}+${edge} variable byl equal -${r}-${edge} variable byr equal ${r}+${edge} variable bzl equal ${subdw}-${edge} variable bzr equal $h+$r+${edge} #============================================================================ # build group #============================================================================ region box block ${bxl} ${bxr} ${byl} ${byr} ${bzl} ${bzr} units box create_box 2 box #box 中原子類型分類有6種 lattice fcc 3.61505 #orient x 1 0 0 orient y 0 1 1 orient z 0 -1 1 region rsubtrate1 block INF INF INF INF ${subdw} ${subup} units box region rsubtrate intersect 2 box rsubtrate1 create_atoms 1 region rsubtrate group gsubtrate region rsubtrate region rsubactive1 block INF INF INF INF ${fixedline} ${subup} units box region rsubactive intersect 2 box rsubactive1 group gsubactive region rsubactive group gsubfixed subtract gsubtrate gsubactive set group gsubfixed type 2 pair_style eam/alloy pair_coeff * * CuAg.eam.alloy Cu Cu #============================================================================ # setting #============================================================================ timestep 0.001 #fs thermo_style custom step temp etotal ke pe press vol neighbor 0.5 bin neigh_modify every 5 delay 0 check yes #============================================================================ # # initial velocities #============================================================================ thermo 100 velocity all create ${T0} 482748 dist gaussian fix 1 all nvt temp ${T0} ${T0} 0.1 drag 0.2 #保持初始溫度，在NVT下弛豫 dump 2 all atom 100 relax_ball.lammpstrj run 800 周期為 p p m 周期為 p m m |
回復(fù)此樓» 猜你喜歡
招收生物學(xué)/細(xì)胞生物學(xué)調(diào)劑
已經(jīng)有5人回復(fù)
-
300求調(diào)劑，材料科學(xué)英一數(shù)二
已經(jīng)有7人回復(fù)
-
070305高分子化學(xué)與物理 304分求調(diào)劑
已經(jīng)有7人回復(fù)
-
289求調(diào)劑
已經(jīng)有13人回復(fù)
-
299求調(diào)劑
已經(jīng)有7人回復(fù)
-
一志愿哈爾濱工業(yè)大學(xué)材料與化工方向336分
已經(jīng)有9人回復(fù)
-
081200-11408-276學(xué)碩求調(diào)劑
已經(jīng)有6人回復(fù)
-
調(diào)劑求院校招收
已經(jīng)有5人回復(fù)
-
調(diào)劑310
已經(jīng)有8人回復(fù)
-
一志愿北京理工大學(xué)本科211材料工程294求調(diào)劑
已經(jīng)有5人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- 更易上手的帶圖形界面的VASP軟件MedeA
已經(jīng)有14人回復(fù)
- python之蟒開啟理論計(jì)算化學(xué)的新時(shí)代
已經(jīng)有22人回復(fù)
2樓2016-08-24 07:53:35
|
謝謝回復(fù)！按理說(shuō)300k馳譽(yù)模型是很工整的。如果是非周期就不會(huì)亂。如圖2，y軸。感覺少了什么約束。 發(fā)自小木蟲Android客戶端 |
3樓2016-08-24 08:21:15
4樓2016-08-24 08:22:44
追夢(mèng)2016
銅蟲 (小有名氣)
- 應(yīng)助: 18 (小學(xué)生)
- 金幣: 81.2
- 紅花: 4
- 帖子: 122
- 在線: 25.2小時(shí)
- 蟲號(hào): 4326873
- 注冊(cè): 2016-01-01
- 性別: GG
- 專業(yè): 微/納機(jī)械系統(tǒng)
5樓2016-08-24 16:06:01
6樓2016-08-25 02:34:28
yangdongjiao
新蟲 (小有名氣)
- 應(yīng)助: 2 (幼兒園)
- 金幣: 1797.8
- 紅花: 2
- 帖子: 79
- 在線: 357.4小時(shí)
- 蟲號(hào): 3779820
- 注冊(cè): 2015-04-01
- 性別: MM
- 專業(yè): 微/納機(jī)械系統(tǒng)
7樓2016-08-25 08:51:40
8樓2016-08-25 09:52:06
yangdongjiao
新蟲 (小有名氣)
- 應(yīng)助: 2 (幼兒園)
- 金幣: 1797.8
- 紅花: 2
- 帖子: 79
- 在線: 357.4小時(shí)
- 蟲號(hào): 3779820
- 注冊(cè): 2015-04-01
- 性別: MM
- 專業(yè): 微/納機(jī)械系統(tǒng)
9樓2016-08-26 08:35:53
enrity
鐵桿木蟲 (正式寫手)
- 應(yīng)助: 2 (幼兒園)
- 金幣: 6668.1
- 散金: 6
- 紅花: 4
- 帖子: 693
- 在線: 129.1小時(shí)
- 蟲號(hào): 1982653
- 注冊(cè): 2012-09-07
- 專業(yè): 固體力學(xué)
|
這個(gè)是你的region和create-atoms命令的問題，邊界是周期性的時(shí)候必須保證尺寸是晶格常數(shù)的整數(shù)倍，并且原子跨過(guò)邊界后不會(huì)發(fā)生重疊，你的結(jié)果是不滿足這些導(dǎo)致的，如果想省事，create-atom用參數(shù)box代替region 發(fā)自小木蟲IOS客戶端 |
tspillar10樓2016-08-28 00:43:29
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 tspillar 的主題更新
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[碩博家園] 求調(diào)劑 330分 085600材料與化工 +3 | gqhhh 2026-03-22 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 070300化學(xué)354求調(diào)劑 +3 | 101次希望 2026-03-28 | 3/150 |
|
|
|
[考研] 本科雙非材料，跨考一志愿華電085801電氣，283求調(diào)劑，任何專業(yè)都可以 +6 | 芝士雪baoo 2026-03-28 | 7/350 |
|
|
|
[考研] 311（085601）求調(diào)劑 +4 | liziyeyeye 2026-03-28 | 4/200 |
|
|
|
[考研] 085701環(huán)境工程，267求調(diào)劑 +16 | minht 2026-03-26 | 16/800 |
|
|
|
[考研]
 材料與化工考研調(diào)劑
+17
|
孅華 2026-03-22 | 17/850 |
|
|
|
[考研] 340求調(diào)劑 +5 | jhx777 2026-03-27 | 5/250 |
|
|
|
[考研] 352分 化工與材料 +5 | 海納百川Ly 2026-03-27 | 5/250 |
|
|
|
[考研] 279 分 求調(diào)劑 +4 | 睡個(gè)好覺_16 2026-03-24 | 4/200 |
|
|
|
[考研] 348求調(diào)劑 +4 | 小懶蟲不懶了 2026-03-27 | 5/250 |
|
|
|
[考研] 315調(diào)劑 +4 | 0860求調(diào)劑 2026-03-26 | 5/250 |
|
|
|
[考研] 一志愿北化求調(diào)劑 +3 | Jsman 2026-03-22 | 3/150 |
|
|
|
[考研] 329求調(diào)劑 +5 | 1() 2026-03-22 | 5/250 |
|
|
|
[考研] 086000生物與醫(yī)藥292求調(diào)劑 +6 | 小小陳小小 2026-03-22 | 9/450 |
|
|
|
[考研] 081700 調(diào)劑 267分 +11 | 迷人的哈哈 2026-03-23 | 11/550 |
|
|
|
[考研] 一志愿河工大 081700 276求調(diào)劑 +4 | 地球繞著太陽(yáng)轉(zhuǎn) 2026-03-23 | 4/200 |
|
|
|
[考研] 303求調(diào)劑 +6 | 藍(lán)山月 2026-03-25 | 6/300 |
|
|
|
[考研] 求調(diào)劑 +3 | 李李不服輸 2026-03-25 | 3/150 |
|
|
|
[考研]
材料考研調(diào)劑生
+3
|
黃粱一夢(mèng)千年 2026-03-24 | 3/150 |
|
|
|
[考研] 341求調(diào)劑(一志愿湖南大學(xué)070300) +5 | 番茄頭--- 2026-03-22 | 6/300 |
|