| 5 | 1/1 | 返回列表 |
| 查看: 871 | 回復: 5 | ||
| 當前只顯示滿足指定條件的回帖,點擊這里查看本話題的所有回帖 | ||
[求助]
分子動力學開始瞬間,整個體系突然縮小 已有1人參與
|
||
|
我將所研究的系統(tǒng)放在一個立方體水箱中,做了20000步優(yōu)化之后,從0K逐步升溫,每次升溫30K, 跑5000步,但是在優(yōu)化完成后開始跑分子動力學的瞬間,整個水箱(包括里面的研究體系)突然縮小,從VMD的動畫中可以明顯看出。 不知道這是否是正,F(xiàn)象,如果不是,究竟是何原因?qū)е碌,如何解決,不知道有沒有人遇到過類似的問題? 謝謝大家! |
木蟲 (正式寫手)
木蟲 (正式寫手)
|
您好!感謝您的回復,我的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ā)表 | |
|---|---|---|---|---|
|
[考研] 299求調(diào)劑 +3 | 霧與海 2026-04-02 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 330分求調(diào)劑 +13 | qzenlc 2026-03-29 | 13/650 |
|
|
[考研] 一志愿085600中科院寧波所276分求調(diào)劑 +26 | 材料學257求調(diào)劑 2026-03-28 | 27/1350 |
|
|
[考研] 材料求調(diào)劑一志愿哈工大324 +10 | 閆旭東 2026-03-28 | 12/600 |
|
|
[考研] 309求調(diào)劑 +8 | 呆菇不是戴夫 2026-04-02 | 8/400 |
|
|
[考研] 348求調(diào)劑 +6 | 吳彥祖24k 2026-04-02 | 6/300 |
|
|
[考研] 初試301,代碼085701環(huán)境工程,本碩一致,四六級已過,有二區(qū)一作,共發(fā)表5篇論文 +6 | axibli 2026-04-01 | 6/300 |
|
|
[考研] 一志愿北京理工大學本科211材料工程294求調(diào)劑 +8 | mikasa的圍巾 2026-03-28 | 9/450 |
|
|
[考研] 292分,材料與化工,申請調(diào)劑 +19 | 程晴之 2026-04-01 | 23/1150 |
|
|
[考研] 265求調(diào)劑11408 +4 | 劉小鹿lu 2026-03-27 | 4/200 |
|
|
[考研] 材料科學與工程339求調(diào)劑 +11 | hyz0119 2026-03-31 | 12/600 |
|
|
[考研] 353求調(diào)劑 +4 | 拉鉤不許變 2026-04-01 | 4/200 |
|
|
[考研] 08工科275求調(diào)劑,可跨考。 +5 | AaAa7420 2026-03-31 | 5/250 |
|
|
[考研] 一志愿中國科學院大學265求調(diào)劑 +8 | 恬淡ye 2026-03-31 | 9/450 |
|
|
[考研] 08工科,295,接受跨專業(yè)調(diào)劑 +6 | lmnlzy 2026-03-31 | 6/300 |
|
|
[考研] 340求調(diào)劑 +4 | 希望如此i 2026-03-31 | 4/200 |
|
|
[考研] 081200-11408-276學碩求調(diào)劑 +4 | 崔wj 2026-03-31 | 4/200 |
|
|
[考研] 福建理工大學材料學院先進合金團隊招收考研調(diào)劑學生 +3 | 大華金商都 2026-03-30 | 4/200 |
|
|
[考研] 340求調(diào)劑 +5 | jhx777 2026-03-27 | 5/250 |
|
|
[考研] 調(diào)劑 +3 | 李嘉圖·S·路 2026-03-27 | 3/150 |
|