當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖

kevin8871

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 28.5
帖子: 8
在線: 4.1小時(shí)
蟲號(hào): 1075839

[交流] 【求助】用NAMD 做NPT模擬遇到下面的問題

本人用NAMD做NPT模擬，遇到了下面的問題
FATAL ERROR: Periodic cell has become too small for original patch grid!
Possible solutions are to restart from a recent checkpoint,
increase margin, or disable useFlexibleCell for liquid simulation.
FATAL ERROR: Periodic cell has become too small for original patch grid!
Possible solutions are to restart from a recent checkpoint,
increase margin, or disable useFlexibleCell for liquid simulation.
這是為什么？下面是我的conf file
structure       3HT_20.psf
coordinates    3HT_20_centered.pdb

set outputname    NPT file for 3 HT

#continuing a run
set inputname    nimi
binVelocities    3HT_20_mini.restart.vel
bincoordinates    3HT_20_mini.restart.coor
extendedsystem    3HT_20_mini.restart.xsc
firsttimestep    6000
numsteps          2506000

#############################################################
## SIMULATION PARAMETERS                                  ##
#############################################################

# Input

paraTypeCharmm       on
parameters          ht_monomer.par
# Force-Field Parameters
exclude          scaled1-4
1-4scaling       1.0
cutoff             12.
switching          on
switchdist       10.
pairlistdist       13.5

# Integrator Parameters
timestep          2.0  ;# 2fs/step
rigidBonds       all  ;# needed for 2fs steps
nonbondedFreq    1
fullElectFrequency  2
stepspercycle    10

# Constant Temperature Control
langevin          on ;# do langevin dynamics
langevinDamping    5    ;# damping coefficient (gamma) of 5/ps
langevinTemp       300
langevinHydrogen off ;# don't couple langevin bath to hydrogens

#PBC
wrapall    on
cellBasisVector1 31.312999725341797 0.0 0.0
cellBasisVector2    0.0 32.558000564575195 0.0
cellBasisVector3    0.0 0 45.066001892089844
cellOrigin       -3.0093722216406604e-6 1.4230939996195957e-5 2.136936927854549e-6

wrapAll          on

# PME (for full-system periodic electrostatics)
PME                yes
PMEGridSizeX       42
PMEGridSizeY       43
PMEGridSizeZ       55

# Constant Pressure Control (variable volume)
useGroupPressure    yes ;# needed for rigidBonds
useFlexibleCell    no  ;# no for water box yes for membrane
useConstantArea    no  ;# no for waterbox, maybe for membrane

langevinPiston       on
langevinPistonTarget  1.01325 ;#  in bar -> 1 atm
langevinPistonPeriod  100.
langevinPistonDecay 50.
langevinPistonTemp 300

# Output
outputName       $outputname

restartfreq       2500    ;# 2500steps = every 5ps
dcdfreq          1000    ;# 1000steps = every 2ps
xstFreq          1000
outputEnergies    500
outputPressure    500

#############################################################
## EXTRA PARAMETERS                                     ##
#############################################################

#############################################################
## EXECUTION SCRIPT                                     ##
#############################################################

# Minimization

run 2500000 ;# 5ns

回復(fù)此樓

» 猜你喜歡

材料求調(diào)劑一志愿哈工大總分298分，前三科223分已經(jīng)有4人回復(fù)
086502化學(xué)工程342求調(diào)劑已經(jīng)有3人回復(fù)
求調(diào)劑已經(jīng)有9人回復(fù)
352分化工與材料已經(jīng)有5人回復(fù)
085602 化工專碩 338分求調(diào)劑已經(jīng)有10人回復(fù)
0856材料化工調(diào)劑總分330 已經(jīng)有10人回復(fù)
330一志愿中國海洋大學(xué) 化學(xué)工程 085602 有讀博意愿求調(diào)劑已經(jīng)有4人回復(fù)
一志愿哈爾濱工業(yè)大學(xué)材料與化工方向336分已經(jīng)有6人回復(fù)
材料求調(diào)劑一志愿哈工大324 已經(jīng)有6人回復(fù)
085600 286分材料求調(diào)劑已經(jīng)有4人回復(fù)

» 本主題相關(guān)價(jià)值貼推薦，對您同樣有幫助:

關(guān)于超級計(jì)算機(jī)的安裝權(quán)限問題：不能用sudo命令，該如何是好？已經(jīng)有18人回復(fù)
word參考文獻(xiàn)-尾注引用，出了點(diǎn)問題，求高手解答已經(jīng)有3人回復(fù)
簽證時(shí)用身份證的問題已經(jīng)有12人回復(fù)
1,3-二羰基化合物單碘代的問題請做過類似反應(yīng)的指教下已經(jīng)有5人回復(fù)
【求助】在NAMD運(yùn)行中使用CHARMM格式的力場文件，但力場參數(shù)來自AMBER力場，會(huì)出問題?/a> 已經(jīng)有4人回復(fù)

【求助】做NAMD的1K4C.pdb時(shí)候出現(xiàn)的錯(cuò)誤已經(jīng)有30人回復(fù)
【求助】Gromacs和NAMD比Amber快很多嗎？如何選擇？已經(jīng)有7人回復(fù)
【專家會(huì)診】bay 討論namd,gmx處理水溶液,蛋白等MD 模擬;討論水及水通道已經(jīng)有50人回復(fù)
【求助】NAMD和Gromacs哪個(gè)簡單？已經(jīng)有5人回復(fù)
【求助】namd tutorial遇到問題已經(jīng)有10人回復(fù)
【求助】NAMD安裝問題已經(jīng)有19人回復(fù)
【在線答疑】關(guān)于NAMD/VMD 已經(jīng)有210人回復(fù)

» 搶金幣啦！回帖就可以得到:

查看全部散金貼

1樓 2010-11-28 15:14:36

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

wally8962

木蟲 (著名寫手)

應(yīng)助: 9 (幼兒園)
金幣: 5366.7
帖子: 1062
在線: 471.1小時(shí)
蟲號(hào): 600024

引用回帖:

Originally posted by kevin8871 at 2010-12-05 13:36:13:
按照wally的建議改了下,多跑了很多,但是到最后同一個(gè)問題又出現(xiàn)了,請問應(yīng)該怎么解決?

繼續(xù)改大。
有的時(shí)候需要改到1000、500

贊一下

回復(fù)此樓

8樓2010-12-05 14:02:16

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

查看全部 11 個(gè)回答

leewz922

金蟲 (小有名氣)

應(yīng)助: 6 (幼兒園)
金幣: 1141.8
帖子: 101
在線: 89.8小時(shí)
蟲號(hào): 609327

★ ★ ★
zh1987hs(金幣+3):謝謝 2010-12-02 22:11:55

把下面這幾行去掉，試試
cellBasisVector1 31.312999725341797 0.0 0.0
cellBasisVector2    0.0 32.558000564575195 0.0
cellBasisVector3    0.0 0 45.066001892089844
cellOrigin       -3.0093722216406604e-6 1.4230939996195957e-5 2.136936927854549e-6

贊一下(1人)

回復(fù)此樓

2樓2010-12-02 20:55:35

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

kevin8871

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 28.5
帖子: 8
在線: 4.1小時(shí)
蟲號(hào): 1075839

我已經(jīng)去掉了，我先做了NVT，然后用NVT的做NPT，但是我做到后面的時(shí)候出錯(cuò)了，我查看了下，我的box的邊長變成了16 16 25，但是我的cutoff 是12，我覺得邊長應(yīng)該大于cutoff的兩倍，所以最少應(yīng)該是24，請問你有什么辦法修改這個(gè)嘛，我現(xiàn)在是不知道怎么弄了。

贊一下

回復(fù)此樓

3樓2010-12-03 14:30:48

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

wally8962

木蟲 (著名寫手)

應(yīng)助: 9 (幼兒園)
金幣: 5366.7
帖子: 1062
在線: 471.1小時(shí)
蟲號(hào): 600024

★ ★ ★
zh1987hs(金幣+3):thank you 2010-12-03 17:01:47
kevin8871(金幣+2): 2010-12-05 09:27:13

langevinPistonPeriod 100.
langevinPistonDecay 50.

太小。改大。

ug建議的是200、100.實(shí)際上可能需要更大。

贊一下(1人)

回復(fù)此樓

4樓2010-12-03 16:52:48

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

最具人氣熱帖推薦 [查看全部]	作者	回/看	最后發(fā)表

[考研] 352分化工與材料 +5	海納百川Ly 2026-03-27	5/250	2026-03-28 03:39 by fmesaito
[考研] 322求調(diào)劑 +4	舊吢 2026-03-24	4/200	2026-03-27 15:38 by 不吃魚的貓
[考研] 269專碩求調(diào)劑 +10	金恩貝 2026-03-21	10/500	2026-03-27 15:10 by caszguilin
[考研] 22408 359分調(diào)劑 +3	Qshers 2026-03-27	3/150	2026-03-27 12:22 by wxiongid
[考研] 334求調(diào)劑 +3	雨清天晴 2026-03-21	3/150	2026-03-27 11:30 by 不吃魚的貓
[考研] 324求調(diào)劑 +5	hanamiko 2026-03-26	5/250	2026-03-27 10:33 by wangjy2002
[考研] 324求調(diào)劑 +4	wysyjs25 2026-03-21	4/200	2026-03-26 20:38 by fmesaito
[考研] 生物學(xué)學(xué)碩，一志愿湖南大學(xué)，初試成績338 +4	YYYYYNNNNN 2026-03-26	4/200	2026-03-26 19:00 by macy2011
[考研] 一志愿211 初試270分求調(diào)劑 +6	谷雨上岸 2026-03-23	7/350	2026-03-26 18:55 by 不吃魚的貓
[考研] 總分293求調(diào)劑 +6	加一一九 2026-03-25	8/400	2026-03-26 13:30 by yujianx
[考研] 一志愿上海交大生物與醫(yī)藥專碩324分，求調(diào)劑 +6	jiajunX 2026-03-22	6/300	2026-03-25 23:05 by licg0208
[考研] 303求調(diào)劑 +6	藍(lán)山月 2026-03-25	6/300	2026-03-25 22:47 by 418490947
[考研] 334分一志愿武理-080500 材料求調(diào)劑 +4	李李不服輸 2026-03-25	4/200	2026-03-25 21:26 by 星空星月
[考研] 考研一志愿蘇州大學(xué)初始315（英一）求調(diào)劑 +3	sbdksD 2026-03-24	4/200	2026-03-25 18:16 by xcjcqu
[考研] 一志愿武理085500機(jī)械專業(yè)總分300求調(diào)劑 +3	an10101 2026-03-24	7/350	2026-03-25 00:00 by 山鬼0-
[考研] 291求調(diào)劑 +3	HanBeiNingZC 2026-03-24	3/150	2026-03-24 16:34 by barlinike
[考博] 26申博自薦 +3	whh869393 2026-03-24	3/150	2026-03-24 09:55 by 21018060
[考研] 284求調(diào)劑 +3	yanzhixue111 2026-03-23	6/300	2026-03-23 22:58 by pswait
[考研] 336化工調(diào)劑 +4	王大坦1 2026-03-23	5/250	2026-03-23 18:32 by allen-yin
[考研] 293求調(diào)劑 +3	濤濤Wjt 2026-03-22	5/250	2026-03-22 22:21 by jiangpengfei