<div>Dear all,</div><div><br></div><div>I'm doing an "slow-growth" alchemical free energy perturbation calculation of the formation of a disulfide bridge between two Cysteines with Gromacs.</div><div><br></div>
<div>I've had tried different ways to combine the topology of both state A and state B, and finally settled with the most direct way -- to "mutate" the atoms that have different partial charges in the two states, and transform the two hydrogen atoms into "dummy" atoms. As suggested by Gromacs manual and some messages from the internet, I put explicitly the OPLS parameters for the bonds, pairs, angles and dihedrals changed from state A to state B, and the grompp didn't give a warning. But when I was testing the production simulation on two processors, there was a warning of fatal error in the log file:</div>
<div><br></div><div>=========================================================================================</div><div>Initializing Domain Decomposition on 2 nodes</div><div>Dynamic load balancing: auto</div><div>Will sort the charge groups at every domain (re)decomposition</div>
<div>Initial maximum inter charge-group distances:</div><div> two-body bonded interactions: 3.774 nm, LJC Pairs NB, atoms 12 205</div><div> multi-body bonded interactions: 0.640 nm, Ryckaert-Bell., atoms 1013 417</div>
<div>Minimum cell size due to bonded interactions: 4.151 nm</div><div>Using 0 separate PME nodes</div><div>Scaling the initial minimum size with 1/0.8 (option -dds) = 1.25</div><div>Optimizing the DD grid for 2 cells with a minimum initial size of 5.189 nm</div>
<div>The maximum allowed number of cells is: X 0 Y 0 Z 0</div><div><br></div><div>-------------------------------------------------------</div><div>Program mdrun, VERSION 4.5.4</div><div>Source code file: domdec.c, line: 6436</div>
<div><br></div><div>Fatal error:</div><div>There is no domain decomposition for 2 nodes that is compatible with the given box and a minimum cell size of 5.18876 nm</div><div>Change the number of nodes or mdrun option -rdd or -dds</div>
<div>Look in the log file for details on the domain decomposition</div><div>For more information and tips for troubleshooting, please check the GROMACS</div><div>website at <a href="http://www.gromacs.org/Documentation/Errors">http://www.gromacs.org/Documentation/Errors</a></div>
<div><br></div><div>=============================================================================================</div><div><br></div><div>The atoms indexed 12 and 205 are a CD1 in Leu1 and another CB in Ala12 respectively. I checked the topology and found no entries containing both these atoms. What's weird is that there's no way for these two atoms to have bonded interactions, but it says in the log that they have a two-body bonded interaction with a distance of 3.774 nm, which I cannot understand. </div>
<div><br></div><div>Does anyone have an explanation what this could mean? Any suggestion is appreciated!</div><div><br></div><div>Thanks in advance!</div><div><br></div><div><br></div><div>Xiaoxiao He</div><div>Oct. 15, 2011</div>
<div><br></div><div><br></div><div><br></div><div>Attached are the input parameters for the production md:</div><div>=====================================================================</div><div>-Input Parameters:</div>
<div> integrator = sd</div><div> nsteps = 2000000</div><div> init_step = 0</div><div> ns_type = Grid</div><div> nstlist = 10</div><div> ndelta = 2</div>
<div> nstcomm = 10</div><div> comm_mode = Linear</div><div> nstlog = 100</div><div> nstxout = 100000</div><div> nstvout = 100000</div><div> nstfout = 100000</div>
<div> nstcalcenergy = 10</div><div> nstenergy = 100</div><div> nstxtcout = 1000</div><div> init_t = 0</div><div> delta_t = 0.0005</div><div> xtcprec = 1000</div>
<div> nkx = 36</div><div> nky = 36</div><div> nkz = 36</div><div> pme_order = 4</div><div> ewald_rtol = 1e-05</div><div> ewald_geometry = 0</div>
<div> epsilon_surface = 0</div><div> optimize_fft = TRUE</div><div> ePBC = xyz</div><div> bPeriodicMols = FALSE</div><div> bContinuation = FALSE</div><div> bShakeSOR = FALSE</div>
<div> etc = No</div><div> nsttcouple = -1</div><div> epc = Berendsen</div><div> epctype = Isotropic</div><div> nstpcouple = 10</div><div> tau_p = 1</div>
<div> ref_p (3x3):</div><div> ref_p[ 0]={ 1.00000e+00, 0.00000e+00, 0.00000e+00}</div><div> ref_p[ 1]={ 0.00000e+00, 1.00000e+00, 0.00000e+00}</div><div> ref_p[ 2]={ 0.00000e+00, 0.00000e+00, 1.00000e+00}</div>
<div> compress (3x3):</div><div> compress[ 0]={ 4.50000e-05, 0.00000e+00, 0.00000e+00}</div><div> compress[ 1]={ 0.00000e+00, 4.50000e-05, 0.00000e+00}</div><div> compress[ 2]={ 0.00000e+00, 0.00000e+00, 4.50000e-05}</div>
<div> refcoord_scaling = No</div><div> posres_com (3):</div><div> posres_com[0]= 0.00000e+00</div><div> posres_com[1]= 0.00000e+00</div><div> posres_com[2]= 0.00000e+00</div><div> posres_comB (3):</div>
<div> posres_comB[0]= 0.00000e+00</div><div> posres_comB[1]= 0.00000e+00</div><div> posres_comB[2]= 0.00000e+00</div><div> andersen_seed = 815131</div><div> rlist = 1.3</div><div>
rlistlong = 1.3</div>
<div> rtpi = 0.05</div><div> coulombtype = PME</div><div> rcoulomb_switch = 0</div><div> rcoulomb = 1.3</div><div> vdwtype = Switch</div><div> rvdw_switch = 0.8</div>
<div> rvdw = 0.9</div><div> epsilon_r = 1</div><div> epsilon_rf = 1</div><div> tabext = 1</div><div> implicit_solvent = No</div><div> gb_algorithm = Still</div>
<div> gb_epsilon_solvent = 80</div><div> nstgbradii = 1</div><div> rgbradii = 1</div><div> gb_saltconc = 0</div><div> gb_obc_alpha = 1</div><div> gb_obc_beta = 0.8</div>
<div> gb_obc_gamma = 4.85</div><div> gb_dielectric_offset = 0.009</div><div> sa_algorithm = Ace-approximation</div><div> sa_surface_tension = 2.05016</div><div> DispCorr = EnerPres</div>
<div> free_energy = yes</div><div> init_lambda = 0</div><div> delta_lambda = 5e-07</div><div> n_foreign_lambda = 0</div><div> sc_alpha = 0.5</div><div> sc_power = 1</div>
<div> sc_sigma = 0.3</div><div> sc_sigma_min = 0.3</div><div> nstdhdl = 10</div><div> separate_dhdl_file = no</div><div> dhdl_derivatives = yes</div><div> dh_hist_size = 0</div>
<div> dh_hist_spacing = 0.1</div><div> nwall = 0</div><div> wall_type = 9-3</div><div> wall_atomtype[0] = -1</div><div> wall_atomtype[1] = -1</div><div> wall_density[0] = 0</div>
<div> wall_density[1] = 0</div><div> wall_ewald_zfac = 3</div><div> pull = no</div><div> disre = No</div><div> disre_weighting = Conservative</div><div> disre_mixed = FALSE</div>
<div> dr_fc = 1000</div><div> dr_tau = 0</div><div> nstdisreout = 100</div><div> orires_fc = 0</div><div> orires_tau = 0</div><div> nstorireout = 100</div>
<div> dihre-fc = 1000</div><div> em_stepsize = 0.01</div><div> em_tol = 10</div><div> niter = 20</div><div> fc_stepsize = 0</div><div> nstcgsteep = 1000</div>
<div> nbfgscorr = 10</div><div> ConstAlg = Lincs</div><div> shake_tol = 0.0001</div><div> lincs_order = 4</div><div> lincs_warnangle = 30</div><div> lincs_iter = 1</div>
<div> bd_fric = 0</div><div> ld_seed = 1993</div><div> cos_accel = 0</div><div> deform (3x3):</div><div> deform[ 0]={ 0.00000e+00, 0.00000e+00, 0.00000e+00}</div><div>
deform[ 1]={ 0.00000e+00, 0.00000e+00, 0.00000e+00}</div><div> deform[ 2]={ 0.00000e+00, 0.00000e+00, 0.00000e+00}</div><div> userint1 = 0</div><div> userint2 = 0</div><div>
userint3 = 0</div><div> userint4 = 0</div><div> userreal1 = 0</div><div> userreal2 = 0</div><div> userreal3 = 0</div><div> userreal4 = 0</div>
<div>grpopts:</div><div> nrdf: 3903 22866</div><div> ref_t: 300 300</div><div> tau_t: 0.1 0.1</div><div>anneal: No No</div><div>ann_npoints: 0 0</div>
<div> acc:<span class="Apple-tab-span" style="white-space:pre">        </span> 0 0 0</div><div> nfreeze: N N N</div><div> energygrp_flags[ 0]: 0</div><div> efield-x:</div>
<div> n = 0</div><div> efield-xt:</div><div> n = 0</div><div> efield-y:</div><div> n = 0</div><div> efield-yt:</div><div> n = 0</div><div> efield-z:</div><div> n = 0</div><div> efield-zt:</div>
<div> n = 0</div><div> bQMMM = FALSE</div><div> QMconstraints = 0</div><div> QMMMscheme = 0</div><div> scalefactor = 1</div><div>qm_opts:</div><div> ngQM = 0</div>
<div><br></div><div>======================================================================================================</div><br clear="all"><div><br></div>-- <br><span style="border-collapse:collapse;font-family:arial, sans-serif;font-size:13px"><div>
<span style="border-collapse:collapse;font-family:arial, sans-serif;font-size:13px">Xiaoxiao HE</span></div>Department of Chemistry<br>The Hong Kong University of Science and Technology<br>Clear Water Bay, Kowloon<br>Hong Kong</span><div>
<span style="border-collapse:collapse;font-family:arial, sans-serif;font-size:13px"><br></span></div><div><span style="border-collapse:collapse;font-family:arial, sans-serif;font-size:13px">Email: <a href="mailto:xxiaohe@ust.hk" target="_blank">xxiaohe@ust.hk</a> Tel: 23585896</span></div>
<br>