| 1 | 1/1 | 返回列表 |
| 查看: 1852 | 回復(fù): 0 | |||
gn02530640銀蟲 (小有名氣)
|
[求助]
NAMD的conf文件,關(guān)于 NVE 和 NVT ensemble 的設(shè)定差別?
|
|
以下2個(gè)conf 文件,一個(gè)是教程里本身包含的 NVE conf 文件,一個(gè)是在網(wǎng)絡(luò)上找到的 NVT conf文件,自己我比對后實(shí)在看不出有什么差異,有了解NAMD的蟲友能說明一下嗎? 非常感謝 ubq-nve.conf ############################################################# ## JOB DESCRIPTION ## ############################################################# # Equilibration of Ubiquitin in a Water Sphere # Use for calculating temperature fluctuations # WARNING: DOES NOT SAVE DCD (and normally you need it) ############################################################# ## ADJUSTABLE PARAMETERS ## ############################################################# structure ../common/ubq_ws.psf coordinates ../common/ubq_ws.pdb set temperature 310 set outputname ubq-nve # Continuing a job from the restart files if {1} { set inputname ../1-2-sphere/ubq_ws_eq binCoordinates $inputname.restart.coor binVelocities $inputname.restart.vel ;# remove the "temperature" entry if you use this! extendedSystem $inputname.restart.xsc } firsttimestep 0 ############################################################# ## SIMULATION PARAMETERS ## ############################################################# # Input paraTypeCharmm on parameters ../common/par_all27_prot_lipid.inp #temperature $temperature # Force-Field Parameters exclude scaled1-4 1-4scaling 1.0 cutoff 12.0 switching on switchdist 10.0 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 if { 0 } { langevin on ;# do langevin dynamics langevinDamping 5 ;# damping coefficient (gamma) of 5/ps langevinTemp $temperature langevinHydrogen off ;# don't couple langevin bath to hydrogens } # Output outputName $outputname #restartfreq 50 ;# 500steps = every 1ps #dcdfreq 50 outputEnergies 50 outputPressure 50 ############################################################# ## EXTRA PARAMETERS ## ############################################################# # Spherical boundary conditions sphericalBC on sphericalBCcenter 30.3081743413, 28.8049907121, 15.353994423 sphericalBCr1 26.0 sphericalBCk1 10 sphericalBCexp1 2 ############################################################# ## EXECUTION SCRIPT ## ############################################################# # Minimization # minimize 100 # reinitvels $temperature run 500000 ;# 1000ps NVT.conf !if you set up a simulation with langevin dynamics but without pressure !controls, it will be in the NVT ensemble. For example, the following is !a modification of the ubiquitin equilibration file from the namd !tutorial to be in the NVT ensemble: ############################################################# ## JOB DESCRIPTION ## ############################################################# # Minimization and Equilibration of # Ubiquitin in a Water Box ############################################################# ## ADJUSTABLE PARAMETERS ## ############################################################# structure ../common/ubq_wb.psf coordinates ../common/ubq_wb.pdb set temperature 310 set outputname ubq_wb_eq firsttimestep 0 ############################################################# ## SIMULATION PARAMETERS ## ############################################################# # Input paraTypeCharmm on parameters ../common/par_all27_prot_lipid.inp temperature $temperature # 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 $temperature langevinHydrogen off ;# don't couple langevin bath to hydrogens # Periodic Boundary Conditions cellBasisVector1 42. 0. 0. cellBasisVector2 0. 44. 0. cellBasisVector3 0. 0 47. cellOrigin 31. 29. 17.5 wrapAll on # PME (for full-system periodic electrostatics) PME yes PMEGridSizeX 32 PMEGridSizeY 32 PMEGridSizeZ 32 useGroupPressure yes ;# needed for rigidBonds useFlexibleCell no useConstantArea no #Note the lack of langevinpiston, which is the pressure control # Output outputName $outputname restartfreq 500 ;# 500steps = every 1ps dcdfreq 250 xstFreq 250 outputEnergies 100 outputPressure 100 ############################################################# ## EXTRA PARAMETERS ## ############################################################# ############################################################# ## EXECUTION SCRIPT ## ############################################################# # Minimization minimize 100 reinitvels $temperature run 2500 ;# 5ps |
| 1 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 286分調(diào)劑 +20 | Faune 2026-03-30 | 22/1100 |
|
|---|---|---|---|---|
|
[考研] 283分材料與化工求調(diào)劑 +7 | 羅KAKA 2026-04-02 | 7/350 |
|
|
[考研] 282求調(diào)劑 +11 | 呼吸都是減肥 2026-04-01 | 11/550 |
|
|
[考研] 找調(diào)劑 +6 | 楚喬喬 2026-04-01 | 6/300 |
|
|
[考研] 289求調(diào)劑 +23 | 新時(shí)代材料 2026-03-27 | 26/1300 |
|
|
[考研] 070305高分子化學(xué)與物理 304分求調(diào)劑 +14 | c297914 2026-03-28 | 14/700 |
|
|
[考研] 生物學(xué)327,求調(diào)劑 +5 | 書上的梅子 2026-04-01 | 6/300 |
|
|
[考研] 08工科275分求調(diào)劑 +13 | AaAa7420 2026-03-31 | 13/650 |
|
|
[考研] 一志愿南昌大學(xué)324求調(diào)劑 +12 | hanamiko 2026-03-27 | 12/600 |
|
|
[考研] 298求調(diào)劑 +4 | 什么是胖頭魚 2026-03-30 | 6/300 |
|
|
[考研] 353求調(diào)劑 +4 | 拉鉤不許變 2026-04-01 | 4/200 |
|
|
[考研] 086502化學(xué)工程342求調(diào)劑 +7 | 阿姨復(fù)古不過 2026-03-27 | 7/350 |
|
|
[考研] 08工科275求調(diào)劑,可跨考。 +5 | AaAa7420 2026-03-31 | 5/250 |
|
|
[考研] 一志愿南昌大學(xué)324求調(diào)劑 +7 | hanamiko 2026-03-30 | 7/350 |
|
|
[考研] 張芳銘-中國農(nóng)業(yè)大學(xué)-環(huán)境工程專碩-298 +9 | 手機(jī)用戶 2026-03-26 | 9/450 |
|
|
[考研] 江蘇蘇北高校誠邀調(diào)劑同學(xué) +3 | zzll406 2026-03-31 | 3/150 |
|
|
[考研] 考研調(diào)劑求助 +7 | 13287130938 2026-03-31 | 7/350 |
|
|
[考研] 一志愿大連理工大學(xué),機(jī)械工程學(xué)碩,341 +3 | 西瓜田的守望者 2026-03-30 | 3/150 |
|
|
[有機(jī)交流]
10+3
|
kaobao456 2026-03-29 | 4/200 |
|
|
[考研] 調(diào)劑 +3 | 李嘉圖·S·路 2026-03-27 | 3/150 |
|