Hello,
I did REMD according to the tutorial. After the simulation is over, an error is reported when using ˇ° extract_300k_traj.xˇ± script to extract the temperature trajectory:
remdtraj requested in trajin but this NETCDF trajectory does not contain temperatures! Trajectory will be processed as replica traj.
The ˇ°extract_300k_traj.xˇ± script is as follows:
#!/bin/csh
ptraj 5EK4_1.prmtop << EOF
trajin remd.mdcrd.001 remdtraj remdtrajtemp 300
trajout remd.300K.mdcrd nobox
go
EOF
The following are the key input files when I perform REMD:
&cntrl
irest=0, ntx=1,
nstlim=500, dt=0.002,
irest=0, ntt=3, gamma_ln=1.0,
temp0=XXXXX, ig=RANDOM_NUMBER,
ntc=2, ntf=2, nscm=1000,
ntb=0, igb=5,
cut=999.0, rgbmax=999.0,
ntpr=100, ntwx=1000, ntwr=10000,
nmropt=1,
numexchg=50000,
ntr=1,
restraint_wt=50.0,
restraintmask=':1-347,374-390.CA',
/
&wt TYPE='END'
/
DISANG=5EK4_chir.dat
REMD simulations are run for a temperature range between 290 - 472K. The starting structure of the 30 replicas of REMD is the 5EKC protein with 30 different structures. 50 ns for each replica.
remd.groupfile:
-O -rem 1 -remlog rem.log -i remd.mdin.001 -o remd.mdout.001 -c 5EK4_1.inpcrd -r remd.rst.001 -x remd.mdcrd.001 -inf remd.mdinfo.001 -p 5EK4_1.prmtop -ref 5EK4_1.inpcrd
-O -rem 1 -remlog rem.log -i remd.mdin.002 -o remd.mdout.002 -c 5EK4_2.inpcrd -r remd.rst.002 -x remd.mdcrd.002 -inf remd.mdinfo.002 -p 5EK4_1.prmtop -ref 5EK4_2.inpcrd
ˇˇ
-O -rem 1 -remlog rem.log -i remd.mdin.028 -o remd.mdout.028 -c 5EK4_28.inpcrd -r remd.rst.028 -x remd.mdcrd.028 -inf remd.mdinfo.028 -p 5EK4_1.prmtop -ref 5EK4_28.inpcrd
-O -rem 1 -remlog rem.log -i remd.mdin.029 -o remd.mdout.029 -c 5EK4_29.inpcrd -r remd.rst.029 -x remd.mdcrd.029 -inf remd.mdinfo.029 -p 5EK4_1.prmtop -ref 5EK4_29.inpcrd
-O -rem 1 -remlog rem.log -i remd.mdin.030 -o remd.mdout.030 -c 5EK4_30.inpcrd -r remd.rst.030 -x remd.mdcrd.030 -inf remd.mdinfo.030 -p 5EK4_1.prmtop -ref 5EK4_30.inpcrd
The command to run REMD is:
srun -n 30 pmemd.cuda.MPI -ng 30 -groupfile remd.groupfile
The output file is as follows:
rem.log:
......
# exchange 50000
1 -1.00 483.93 -4919.86 472.00 472.00 0.00 -1
2 1.01 425.03 -5615.09 421.90 428.80 0.28 -1
3 -1.00 357.92 -6665.98 363.30 363.30 0.26 -1
4 -1.00 431.49 -5213.03 443.00 443.00 0.29 -1
5 0.99 310.04 -7672.94 316.80 311.30 0.25 -1
6 -1.00 381.14 -6215.04 395.00 395.00 0.27 -1
7 0.99 457.53 -5069.34 457.60 450.20 0.29 -1
8 0.99 308.95 -7850.23 305.80 300.50 0.24 -1
9 1.01 310.87 -7623.60 311.30 316.80 0.25 -1
10 -1.00 435.48 -5500.29 435.90 435.90 0.29 -1
11 0.99 422.83 -5597.49 428.80 421.90 0.28 -1
12 -1.00 320.61 -7401.61 322.30 322.30 0.25 -1
13 1.01 448.75 -5117.91 450.20 457.60 0.28 -1
14 1.01 300.06 -7902.66 300.50 305.80 0.25 -1
15 -1.00 407.83 -6045.13 401.60 401.60 0.27 -1
16 -1.00 411.97 -5993.03 408.30 408.30 0.28 -1
17 -1.00 342.54 -7093.49 345.30 345.30 0.26 -1
18 -1.00 374.51 -6618.57 369.50 369.50 0.27 -1
19 -1.00 406.04 -5767.74 415.00 415.00 0.28 -1
20 -1.00 347.37 -7103.56 339.40 339.40 0.26 -1
21 -1.00 328.74 -7117.47 333.60 333.60 0.25 -1
22 -1.00 354.91 -6798.13 357.20 357.20 0.26 -1
23 -1.00 377.97 -6340.68 382.10 382.10 0.28 -1
24 0.99 294.19 -7967.01 295.20 290.00 0.25 -1
25 1.01 291.39 -7976.23 290.00 295.20 0.25 -1
26 -1.00 326.46 -7318.99 327.90 327.90 0.25 -1
27 -1.00 380.62 -6394.33 375.70 375.70 0.26 -1
28 -1.00 384.89 -6096.66 388.50 388.50 0.28 -1
29 -1.00 349.34 -6909.41 351.20 351.20 0.26 -1
30 -1.00 468.11 -5049.65 465.00 465.00 0.32 -1
remd.mdout.010 :
......
NSTEP = 25000000 TIME(PS) = 50000.000 TEMP(K) = 433.81 PRESS = 0.0
Etot = 1429.7036 EKtot = 6800.4443 EPtot = -5370.7407
BOND = 1641.3302 ANGLE = 4534.2074 DIHED = 5377.9737
1-4 NB = 1573.3763 1-4 EEL = 17498.4438 VDWAALS = -3220.7094
EELEC = -28666.4125 EGB = -4444.1999 RESTRAINT = 335.2497
EAMBER (non-restraint) = -5705.9904
TEMP0 = 435.9000 REPNUM = 10 EXCHANGE# = 50000
------------------------------------------------------------------------------
NMR restraints: Bond = 0.000 Angle = 0.000 Torsion = 18.340
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Fri Apr 23 2021 - 18:30:02 PDT
