版塊導(dǎo)航: 正在加載中...

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實(shí)名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號(hào)能夠正常使用，請(qǐng)盡快對(duì)帳號(hào)進(jìn)行手機(jī)號(hào)驗(yàn)證，感謝您的理解與支持！

24小時(shí)熱門版塊排行榜

返回列表

歐陽少恭

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 49
帖子: 4
在線: 2.8小時(shí)
蟲號(hào): 1361266
注冊(cè): 2011-08-05
專業(yè): 無機(jī)材料化學(xué)

[求助] Lammps 模擬NPT系統(tǒng)出錯(cuò)

這是我的input file，想要在NVT的基礎(chǔ)上接著運(yùn)行NPT，但是使用如下fix命令后無論是溫度還是壓強(qiáng)都沒有被fix，如果刪除fix rigid 那么分子之間的bond就消失了。。。使用fix rigid/npt 又說 invalid style。。。只有這種形式才能運(yùn)行，但是溫度和壓強(qiáng)又不fix。。。請(qǐng)各位大神幫幫忙~~~~~

input file：
————————————————————————————————————————————————————————————————————————
dimension 2
boundary p p p

neighbor 0.3 bin
neigh_modify delay 1
timestep 0.002

units lj
atom_style full

read_restart restart.dimer_final

group dimer molecule <> 1 3200

thermo_style custom step temp press vol enthalpy pe ke etotal lx ly lz atoms

pair_style    hybrid/overlay lj/cut 6.0 yukawa 1.0 6.0
pair_coeff    * * lj/cut 1.0 1.0
pair_coeff    * * yukawa 2.250

fix  10 dimer rigid molecule
fix 1 all press/berendsen iso 0.0018 0.0018 1.0
fix 2 all npt temp 0.045 0.045 1.0 iso 0.0018 0.0018 1.0

neigh_modify exclude molecule dimer

fix    5 all enforce2d
——————————————————————————————————————————————————————————————————————————

output file 就這樣：
————————————————————————————————————————————————————————————————————————
LAMMPS (9 Nov 2011)
Reading restart file ...
  orthogonal box = (0 0 0) to (357.77 357.77 2.2)
  1 by 5 by 1 processor grid
  6400 atoms
  3200 bonds
Finding 1-2 1-3 1-4 neighbors ...
  1 = max # of 1-2 neighbors
  0 = max # of 1-3 neighbors
  0 = max # of 1-4 neighbors
  1 = max # of special neighbors
6400 atoms in group dimer
3200 rigid bodies with 6400 atoms
WARNING: One or more atoms are time integrated more than once (modify.cpp:240)
Setting up run ...
Memory usage per processor = 4.26651 Mbytes
Step Temp Press Volume Enthalpy PotEng KinEng TotEng Lx Ly Lz Atoms
11000000  0.045404677 -0.0048495 127999.37  -0.16623931 -0.1032962  0.034046413 -0.069249787    357.77    357.77       2.2    6400
11001000  0.035714946 0.0025690215    149326.3 -0.030062248  -0.11678389  0.026780629 -0.090003259 386.42761 386.42761       2.2    6400
11002000  0.030264793 0.0023758957 176161.74 -0.033412274  -0.12150331  0.022693866 -0.098809449 419.71626 419.71626       2.2    6400
11003000  0.029004975 0.0022623105 192742.28 -0.033463396  -0.12334429  0.021749199 -0.1015951 439.02424 439.02424       2.2    6400
11004000  0.028354983 0.0022907102 208286.37 -0.028235237  -0.12404763  0.021261807  -0.10278582 456.38402 456.38402       2.2    6400
11005000 0.02802591 0.0023240801 230934.88 -0.019681906  -0.12455808  0.021015054  -0.10354302 480.55684 480.55684       2.2    6400
11006000  0.027851188 0.0022193808    248575.1 -0.017691289  -0.12477576  0.020884039  -0.10389173 498.57306 498.57306       2.2    6400
11007000  0.027480945 0.0022516638 249899.94 -0.015969955  -0.12449678  0.020606415  -0.10389037 499.89993 499.89993       2.2    6400
11008000  0.027463389 0.0022537957 268442.03 -0.0096305365  -0.12475714  0.020593251  -0.10416389 518.11391 518.11391       2.2    6400
————————————————————————————————————————————————————————————————————————————————
多謝各位大神。。。

回復(fù)此樓

» 猜你喜歡

求調(diào)劑院校信息已經(jīng)有4人回復(fù)
085600材料與化工306 已經(jīng)有4人回復(fù)
286求調(diào)劑已經(jīng)有10人回復(fù)
328求調(diào)劑，英語六級(jí)551，有科研經(jīng)歷已經(jīng)有9人回復(fù)
一志愿北京化工大學(xué)070300 學(xué)碩336求調(diào)劑已經(jīng)有4人回復(fù)
286分人工智能專業(yè)請(qǐng)求調(diào)劑愿意跨考！已經(jīng)有8人回復(fù)
資源與環(huán)境調(diào)劑申請(qǐng)(333分) 已經(jīng)有5人回復(fù)
280求調(diào)劑已經(jīng)有12人回復(fù)
269專碩求調(diào)劑已經(jīng)有5人回復(fù)
材料學(xué)碩301分求調(diào)劑已經(jīng)有7人回復(fù)

» 本主題相關(guān)商家推薦: (我也要在這里推廣)

1樓 2015-06-10 04:32:02

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

智能機(jī)器人

Robot (super robot)

我們都愛小木蟲

找到一些相關(guān)的精華帖子，希望有用哦~

lammps中LJ勢(shì)能下單位問題已經(jīng)有5人回復(fù)
Lammps計(jì)算油水界面，油相密度、平衡等問題求指導(dǎo) 已經(jīng)有6人回復(fù)
關(guān)于npt系綜設(shè)定？已經(jīng)有5人回復(fù)
lammps并行安裝最后一步make serial出錯(cuò) 已經(jīng)有17人回復(fù)
lammps做拉伸試驗(yàn)，求大神進(jìn)來指導(dǎo)下（烈日炎炎，有勞了`(*∩_∩*)′）已經(jīng)有14人回復(fù)
lammps中的npt系綜下逐步降溫問題已經(jīng)有3人回復(fù)
lammps中界面的建立已經(jīng)有9人回復(fù)
Forcite 弛豫已經(jīng)有6人回復(fù)
分子模擬中MSD的問題討論已經(jīng)有25人回復(fù)
急！lammps安裝修改makefile.g++文件后運(yùn)行出錯(cuò)。已經(jīng)有12人回復(fù)
lammps實(shí)驗(yàn)拉伸Al納米線的in文件的錯(cuò)誤提示好久找不出？已經(jīng)有8人回復(fù)
關(guān)于lammps fix npt/deform做拉伸時(shí)溫度的變化已經(jīng)有6人回復(fù)
LAMMPS溫度－步長曲線震蕩已經(jīng)有11人回復(fù)
lammps fix npt中控溫控壓方法已經(jīng)有5人回復(fù)
NPT系綜動(dòng)力學(xué) 已經(jīng)有19人回復(fù)
LAMMPS運(yùn)行結(jié)果出現(xiàn)nan 已經(jīng)有20人回復(fù)
第一原理動(dòng)力學(xué)的系綜選擇，和各種系綜的輸入文件Example 已經(jīng)有13人回復(fù)
NPT系綜的平衡問題已經(jīng)有10人回復(fù)
【討論】npt壓力控制問題已經(jīng)有10人回復(fù)
【討論】系綜問題已經(jīng)有32人回復(fù)
【求助】NVT弛豫問題已經(jīng)有3人回復(fù)
【討論】LAMMPS的replicate command與creat box區(qū)別已經(jīng)有3人回復(fù)
【討論】lammps運(yùn)行eam/fs出錯(cuò) 已經(jīng)有38人回復(fù)

點(diǎn)擊這里搜索更多相關(guān)資源

科研從小木蟲開始，人人為我，我為人人

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主歐陽少恭的主題更新