 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 3044 | 回復: 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 |
回復此樓» 猜你喜歡
299求調劑
已經(jīng)有8人回復
-
一志愿北京理工大學本科211材料工程294求調劑
已經(jīng)有6人回復
-
300求調劑，材料科學英一數(shù)二
已經(jīng)有8人回復
-
招收生物學/細胞生物學調劑
已經(jīng)有5人回復
-
070305高分子化學與物理 304分求調劑
已經(jīng)有7人回復
-
289求調劑
已經(jīng)有13人回復
-
一志愿哈爾濱工業(yè)大學材料與化工方向336分
已經(jīng)有9人回復
-
081200-11408-276學碩求調劑
已經(jīng)有6人回復
-
調劑求院校招收
已經(jīng)有5人回復
-
調劑310
已經(jīng)有8人回復
» 本主題相關價值貼推薦，對您同樣有幫助:
- 更易上手的帶圖形界面的VASP軟件MedeA
已經(jīng)有14人回復
- python之蟒開啟理論計算化學的新時代
已經(jīng)有22人回復
yangdongjiao
新蟲 (小有名氣)
- 應助: 2 (幼兒園)
- 金幣: 1797.8
- 紅花: 2
- 帖子: 79
- 在線: 357.4小時
- 蟲號: 3779820
- 注冊: 2015-04-01
- 性別: MM
- 專業(yè): 微/納機械系統(tǒng)
9樓2016-08-26 08:35:53
2樓2016-08-24 07:53:35
3樓2016-08-24 08:21:15
4樓2016-08-24 08:22:44
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 295材料工程專碩求調劑 +7 | 1428151015 2026-03-27 | 7/350 |
|
|---|---|---|---|---|
|
[考研] 壓國家一區(qū)線，求導師收留，有恩必謝！ +7 | 迷人的哈哈 2026-03-28 | 7/350 |
|
|
|
[考研] 275求調劑 +10 | Micky11223 2026-03-25 | 14/700 |
|
|
|
[考研] 材料與化工272求調劑 +9 | 阿斯蒂芬2004 2026-03-28 | 9/450 |
|
|
|
[考研] 311求調劑 +3 | 希望上岸阿小楊 2026-03-23 | 3/150 |
|
|
|
[考研] 291求調劑 +7 | 孅華 2026-03-22 | 7/350 |
|
|
|
[考研] 266分求材料化工冶金礦業(yè)等專業(yè)的調劑 +4 | 哇呼哼呼哼 2026-03-26 | 4/200 |
|
|
|
[考研] 307求調劑 +8 | 超級伊昂大王 2026-03-24 | 9/450 |
|
|
|
[考研] 274求調劑 +17 | 顧九笙要謙虛 2026-03-24 | 23/1150 |
|
|
|
[考研] 調劑 +3 | 李嘉圖·S·路 2026-03-27 | 3/150 |
|
|
|
[考研] 調劑推薦 +5 | 清酒714 2026-03-26 | 6/300 |
|
|
|
[考研] 324求調劑 +8 | hanamiko 2026-03-26 | 10/500 |
|
|
|
[考研] 材料科學與工程 317求調劑 +4 | JKSOIID 2026-03-26 | 4/200 |
|
|
|
[考研] 化學調劑一志愿上海交通大學336分-本科上海211 +4 | 小魚愛有機 2026-03-25 | 4/200 |
|
|
|
[考研] 生物技術與工程 +3 | 1294608413 2026-03-25 | 4/200 |
|
|
|
[考研] 0854AI CV方向招收調劑 +4 | 章小魚567 2026-03-23 | 4/200 |
|
|
|
[考研] 347求調劑 +4 | L when 2026-03-25 | 4/200 |
|
|
|
[考研] 284求調劑 +15 | Zhao anqi 2026-03-22 | 15/750 |
|
|
|
[考研] 調劑 +4 | 13853210211 2026-03-24 | 4/200 |
|
|
|
[考研] 275求調劑 +6 | shansx 2026-03-22 | 8/400 |
|