| 5 | 1/1 | 返回列表 |
| 查看: 868 | 回復: 5 | ||
| 當前只顯示滿足指定條件的回帖,點擊這里查看本話題的所有回帖 | ||
[求助]
分子動力學開始瞬間,整個體系突然縮小 已有1人參與
|
||
|
我將所研究的系統(tǒng)放在一個立方體水箱中,做了20000步優(yōu)化之后,從0K逐步升溫,每次升溫30K, 跑5000步,但是在優(yōu)化完成后開始跑分子動力學的瞬間,整個水箱(包括里面的研究體系)突然縮小,從VMD的動畫中可以明顯看出。 不知道這是否是正,F(xiàn)象,如果不是,究竟是何原因導致的,如何解決,不知道有沒有人遇到過類似的問題? 謝謝大家! |
木蟲 (正式寫手)
木蟲 (正式寫手)
|
您好!感謝您的回復,我的conf文件如下: ############################################################# ## JOB DESCRIPTION ## ############################################################# # Minimization and Equilibration of # the whole complex in a Water Box ############################################################# ## ADJUSTABLE PARAMETERS ## ############################################################# structure ionized.psf coordinates ionized.pdb set temperature 300 set outputname whole_out firsttimestep 0 ############################################################# ## SIMULATION PARAMETERS ## ############################################################# # Input paraTypeCharmm on parameters par_all27_prot_lipid_na.inp temperature $temperature # Force-Field Parameters exclude scaled1-4 1-4scaling 1.0 cutoff 14.0 switching on switchdist 10.0 pairlistdist 16.0 # Integrator Parameters timestep 2 ;# 2fs/step rigidBonds all ;# needed for 2fs steps nonbondedFreq 1 fullElectFrequency 2 stepspercycle 10 # Constant Temperature Control langevin on ;# do langevin dynamics langevinDamping 1 ;# damping coefficient (gamma) of 1/ps langevinTemp $temperature langevinHydrogen off ;# don't couple langevin bath to hydrogens # Periodic Boundary Conditions cellBasisVector1 152.1 0.0 0.0 cellBasisVector2 0.0 158.9 0.0 cellBasisVector3 0.0 0.0 152.6 cellOrigin -52.6 -16.9 -29.8 wrapAll on # PME (for full-system periodic electrostatics) PME yes PMEGridSpacing 1.0 #manual grid definition #PMEGridSizeX 45 #PMEGridSizeY 45 #PMEGridSizeZ 48 # Constant Pressure Control (variable volume) useGroupPressure yes ;# needed for rigidBonds useFlexibleCell yes useConstantArea no langevinPiston on langevinPistonTarget 1.01325 ;# in bar -> 1 atm langevinPistonPeriod 100.0 langevinPistonDecay 50.0 langevinPistonTemp $temperature # Output outputName $outputname restartfreq 500 ;# 500steps = every 1ps dcdfreq 250 xstFreq 250 outputEnergies 100 outputPressure 100 ############################################################# ## EXTRA PARAMETERS ## ############################################################# ############################################################# ## EXECUTION SCRIPT ## ############################################################# # Minimization minimize 20000 for { set TEMP 0 } { $TEMP < 300 } { incr TEMP 30 } { reinitvels $TEMP langevinTemp $TEMP run 5000 } reinitvels $temperature langevinTemp 300 langevinPistonTemp 300 run 1000000 ;# 2ns |
木蟲 (正式寫手)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 312求調劑 +3 | 上岸吧ZJY 2026-03-22 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 各位老師您好:本人初試372分 +5 | jj涌77 2026-03-25 | 6/300 |
|
|
[考研] 求調劑 +3 | 李李不服輸 2026-03-25 | 3/150 |
|
|
[考研] 293求調劑 +7 | 加一一九 2026-03-24 | 7/350 |
|
|
[考研] 292求調劑 +4 | 鵝鵝鵝額額額額?/a> 2026-03-25 | 4/200 |
|
|
[考研] 一志愿中南大學化學學碩0703總分337求調劑 +6 | niko- 2026-03-22 | 6/300 |
|
|
[考研] B區(qū)考研調劑 +4 | yqdszhdap- 2026-03-22 | 5/250 |
|
|
[考研] 食品專碩 一志愿雙一流 328 +3 | xiaom99 2026-03-21 | 4/200 |
|
|
[有機交流]
20+3
|
FENGSHUJEI 2026-03-23 | 5/250 |
|
|
[考研] 335求調劑 +4 | yuyu宇 2026-03-23 | 5/250 |
|
|
[考研] 341求調劑(一志愿湖南大學070300) +5 | 番茄頭--- 2026-03-22 | 6/300 |
|
|
[考研] 284求調劑 +3 | yanzhixue111 2026-03-23 | 6/300 |
|
|
[論文投稿] 急發(fā)核心期刊論文 +3 | 賢達問津 2026-03-23 | 5/250 |
|
|
[考研] 352求調劑 +3 | 大米飯! 2026-03-22 | 3/150 |
|
|
[考研] 311求調劑 +3 | 26研0 2026-03-20 | 3/150 |
|
|
[考研] 一志愿南大,0703化學,分數(shù)336,求調劑 +3 | 收到VS 2026-03-21 | 3/150 |
|
|
[考研] 296求調劑 +4 | www_q 2026-03-20 | 4/200 |
|
|
[考研] 22 350 本科985求調劑,求老登收留 +3 | 李軼男003 2026-03-20 | 3/150 |
|
|
[考研] 南京大學化學376求調劑 +3 | hisfailed 2026-03-19 | 6/300 |
|
|
[考研] 中南大學化學學碩337求調劑 +3 | niko- 2026-03-19 | 6/300 |
|