 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 3047 | 回復: 16 | ||
| 當前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖 | ||
[求助]
lammps 簡單模型弛豫時周期邊界出問題。 已有2人參與
|
||
|
問題：利用lammps建立一個矩形，x和y方向為周期，z方向為非周期，(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 |
回復此樓» 猜你喜歡
337求調(diào)劑
已經(jīng)有3人回復
-
329求調(diào)劑，一志愿西北工業(yè)大學，材料工程（085601）
已經(jīng)有10人回復
-
一志愿北京化工大學材料與化工（085600）296求調(diào)劑
已經(jīng)有16人回復
-
291求調(diào)劑
已經(jīng)有6人回復
-
299求調(diào)劑
已經(jīng)有10人回復
-
面上5B能上會嗎？
已經(jīng)有4人回復
-
322求調(diào)劑：一志愿湖南大學 材料與化工（085600），已過六級。
已經(jīng)有4人回復
-
085600 材料與化工 329分求調(diào)劑
已經(jīng)有14人回復
-
材料科學與工程求調(diào)劑
已經(jīng)有3人回復
-
289求調(diào)劑
已經(jīng)有5人回復
» 本主題相關價值貼推薦，對您同樣有幫助:
- 更易上手的帶圖形界面的VASP軟件MedeA
已經(jīng)有14人回復
- python之蟒開啟理論計算化學的新時代
已經(jīng)有22人回復
3樓2016-08-24 08:21:15
2樓2016-08-24 07:53:35
4樓2016-08-24 08:22:44
5樓2016-08-24 16:06:01
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿北京化工大學材料與化工（085600）296求調(diào)劑 +16 | 稻妻小編 2026-03-26 | 16/800 |
|
|---|---|---|---|---|
|
|
[考研] 085600 材料與化工 329分求調(diào)劑 +14 | Mr. Z 2026-03-25 | 14/700 |
|
|
[考研] 一志愿211，335分，0856，求調(diào)劑院校和導師 +5 | 傾____蕭 2026-03-27 | 6/300 |
|
|
|
[考研] 070305高分子化學與物理 304分求調(diào)劑 +12 | c297914 2026-03-28 | 12/600 |
|
|
|
[碩博家園] 招收生物學/細胞生物學調(diào)劑 +4 | IceGuo 2026-03-26 | 5/250 |
|
|
|
[考研] 一志愿211 初試270分 求調(diào)劑 +8 | 谷雨上岸 2026-03-23 | 9/450 |
|
|
|
[考研] 0703 化學 求調(diào)劑，一志愿山東大學 342 分 +4 | Shern—- 2026-03-28 | 4/200 |
|
|
|
[考研] 數(shù)一英一271專碩（085401）求調(diào)劑，可跨 +7 | 前行必有光 2026-03-28 | 8/400 |
|
|
|
[考研] 071000生物學求調(diào)劑，初試成績343 +7 | 小小甜面團 2026-03-25 | 7/350 |
|
|
|
[考研] 328求調(diào)劑 +7 | 嗯滴的基本都 2026-03-27 | 7/350 |
|
|
|
[考研] 求調(diào)劑推薦 材料 304 +15 | 荷包蛋hyj 2026-03-26 | 15/750 |
|
|
|
[考研] 085600材料與化工調(diào)劑 +10 | A-哆啦Z夢 2026-03-23 | 16/800 |
|
|
|
[考研] 化學308分求調(diào)劑 +8 | 你好明天你好 2026-03-23 | 9/450 |
|
|
|
[考研] 考研調(diào)劑 +10 | 呼呼？~+123456 2026-03-24 | 10/500 |
|
|
|
[考研] 求調(diào)劑，一志愿 南京航空航天大學大學 ，080500材料科學與工程學碩 +4 | @taotao 2026-03-26 | 5/250 |
|
|
|
[考研] 085601求調(diào)劑總分293英一數(shù)二 +4 | 鋼鐵大炮 2026-03-24 | 4/200 |
|
|
|
[考研] 一志愿河工大 081700 276求調(diào)劑 +4 | 地球繞著太陽轉(zhuǎn) 2026-03-23 | 4/200 |
|
|
|
[考研] 一志愿天津大學339材料與化工求調(diào)劑 +3 | 江往賣魚 2026-03-26 | 3/150 |
|
|
|
[考研] 機械學碩總分317求調(diào)劑！�。。� +4 | Acaciad 2026-03-25 | 4/200 |
|
|
|
[考研] 化工專碩求調(diào)劑 +3 | question挽風 2026-03-24 | 3/150 |
|