[gmx-users] Using Charmm 36 forcefield on Gromacs 4.6 or higher
akk5r
akk5r at virginia.edu
Thu Aug 1 22:10:08 CEST 2013
Hi All,
I am trying to run a membrane protein simulation on Gromacs 4.6 using the
Charmm 36 force field parameters.
I found the following parameters on the gromacs mailing list for Charmm 36:
*; nblist cut-off
rlist = 1.2
; long-range cut-off for switched potentials
rlistlong = 1.4
; OPTIONS FOR ELECTROSTATICS AND VDW
; Method for doing electrostatics
coulombtype = PME
rcoulomb-switch = 0
rcoulomb = 1.2
; Relative dielectric constant for the medium and the reaction field
epsilon_r = 1
epsilon_rf = 1
; Method for doing Van der Waals
vdw-type = switch
; cut-off lengths
rvdw-switch = 0.8
rvdw = 1.2
; Apply long range dispersion corrections for Energy and Pressure
DispCorr = No *
Putting these setting into my .mdp file, I then get this error on grompp:
*ERROR 1 [file nvt.mdp]:
With Verlet lists only cut-off LJ interactions are supported*
I then turned off the Verlet cut-off scheme and was able to create a .tpr
file.
Then I ran my .tpr on mdrun and got the following error:
*Program mdrun, VERSION 4.6.1
Source code file: /home/akk5r/Downloads/gromacs-4.6.1/src/kernel/runner.c,
line: 824
Fatal error:
OpenMP threads have been requested with cut-off scheme Group, but these are
only supported with cut-off scheme Verlet
For more information and tips for troubleshooting, please check the GROMACS
website at http://www.gromacs.org/Documentation/Errors*
At this point, I don't know what to do. I have to run this job on a cluster
system due to the size of the system. Does any one know how to run Charmm 36
simulations on Gromacs 4.6 or higher using Verlet cut-off scheme?
Best Wishes,
Ali
--
View this message in context: http://gromacs.5086.x6.nabble.com/Using-Charmm-36-forcefield-on-Gromacs-4-6-or-higher-tp5010290.html
Sent from the GROMACS Users Forum mailing list archive at Nabble.com.
More information about the gromacs.org_gmx-users
mailing list