The Gapsys *et al.* soft-core is formulated to act on the level of van der Waals and electrostatic forces:
the non-bonded interactions are linearized at a point defined as, :math:`r_{scLJ}` or :math:`r_{scQ}`, respectively.
The linearization point depends on the state of the system as controlled by the :math:`\lambda` parameter and
-two parameters :math:`\alpha_Q` (set with ``sc-gapsys-scale-linpoint-Q``) and :math:`\alpha_{LJ}` (set with ``sc-gapsys-scale-linpoint-LJ``).
+two parameters :math:`\alpha_Q` (set with ``sc-gapsys-scale-linpoint-q``) and :math:`\alpha_{LJ}` (set with ``sc-gapsys-scale-linpoint-lj``).
The dependence on :math:`\lambda` guarantees that the end-states are properly represented by their hard-core potentials.
:numref:`Fig. %s <fig-gapsyssc>` illustrates the behaviour of the linearization point, forces and integrated potential energies with respect
-to the parameters :math:`\alpha_Q` and :math:`\alpha_{LJ}`. The optimal choices of the parameter values have been systematically explored in :ref:`185 <refGapsys2012>`. These recommended values are set by default when ``sc-function=gapsys`` is selected: ``sc-gapsys-scale-linpoint-Q=0.3`` and ``sc-gapsys-scale-linpoint-LJ=0.85``.
+to the parameters :math:`\alpha_Q` and :math:`\alpha_{LJ}`. The optimal choices of the parameter values have been systematically explored in :ref:`185 <refGapsys2012>`. These recommended values are set by default when ``sc-function=gapsys`` is selected: ``sc-gapsys-scale-linpoint-q=0.3`` and ``sc-gapsys-scale-linpoint-lj=0.85``.
.. _fig-gapsyssc:
than :mdp:`sc-sigma`.
Used only with `sc-function=beutler`
-.. mdp:: sc-gapsys-scale-linpoint-LJ
+.. mdp:: sc-gapsys-scale-linpoint-lj
(0.85)
for `sc-function=gapsys` it is the unitless alphaLJ parameter.
van der Waals interactions.
Used only with `sc-function=gapsys`
-.. mdp:: sc-gapsys-scale-linpoint-Q
+.. mdp:: sc-gapsys-scale-linpoint-q
(0.3) [nm/e^2]
For `sc-function=gapsys` the alphaQ parameter
result in the standard hard-core Coulombic interactions.
Used only with `sc-function=gapsys`
-.. mdp:: sc-gapsys-sigma-LJ
+.. mdp:: sc-gapsys-sigma-lj
(0.3) [nm]
for `sc-function=gapsys` the soft-core sigma for particles
fep->sc_r_power = get_ereal(&inp, "sc-r-power", 6.0, wi);
fep->sc_sigma = get_ereal(&inp, "sc-sigma", 0.3, wi);
fep->bScCoul = (getEnum<Boolean>(&inp, "sc-coul", wi) != Boolean::No);
- fep->scGapsysScaleLinpointLJ = get_ereal(&inp, "sc-gapsys-scale-linpoint-LJ", 0.85, wi);
- fep->scGapsysScaleLinpointQ = get_ereal(&inp, "sc-gapsys-scale-linpoint-Q", 0.3, wi);
- fep->scGapsysSigmaLJ = get_ereal(&inp, "sc-gapsys-sigma-LJ", 0.3, wi);
+ fep->scGapsysScaleLinpointLJ = get_ereal(&inp, "sc-gapsys-scale-linpoint-lj", 0.85, wi);
+ fep->scGapsysScaleLinpointQ = get_ereal(&inp, "sc-gapsys-scale-linpoint-q", 0.3, wi);
+ fep->scGapsysSigmaLJ = get_ereal(&inp, "sc-gapsys-sigma-lj", 0.3, wi);
fep->dh_hist_size = get_eint(&inp, "dh_hist_size", 0, wi);
fep->dh_hist_spacing = get_ereal(&inp, "dh_hist_spacing", 0.1, wi);
fep->separate_dhdl_file = getEnum<SeparateDhdlFile>(&inp, "separate-dhdl-file", wi);
sc-r-power = 6
sc-sigma = 0.3
sc-coul = no
-sc-gapsys-scale-linpoint-LJ = 0.85
-sc-gapsys-scale-linpoint-Q = 0.3
-sc-gapsys-sigma-LJ = 0.3
+sc-gapsys-scale-linpoint-lj = 0.85
+sc-gapsys-scale-linpoint-q = 0.3
+sc-gapsys-sigma-lj = 0.3
separate-dhdl-file = yes
dhdl-derivatives = yes
dh_hist_size = 0
sc-r-power = 6
sc-sigma = 0.3
sc-coul = no
-sc-gapsys-scale-linpoint-LJ = 0.85
-sc-gapsys-scale-linpoint-Q = 0.3
-sc-gapsys-sigma-LJ = 0.3
+sc-gapsys-scale-linpoint-lj = 0.85
+sc-gapsys-scale-linpoint-q = 0.3
+sc-gapsys-sigma-lj = 0.3
separate-dhdl-file = yes
dhdl-derivatives = yes
dh_hist_size = 0
sc-r-power = 6
sc-sigma = 0.3
sc-coul = no
-sc-gapsys-scale-linpoint-LJ = 0.85
-sc-gapsys-scale-linpoint-Q = 0.3
-sc-gapsys-sigma-LJ = 0.3
+sc-gapsys-scale-linpoint-lj = 0.85
+sc-gapsys-scale-linpoint-q = 0.3
+sc-gapsys-sigma-lj = 0.3
separate-dhdl-file = yes
dhdl-derivatives = yes
dh_hist_size = 0
sc-r-power = 6
sc-sigma = 0.3
sc-coul = no
-sc-gapsys-scale-linpoint-LJ = 0.85
-sc-gapsys-scale-linpoint-Q = 0.3
-sc-gapsys-sigma-LJ = 0.3
+sc-gapsys-scale-linpoint-lj = 0.85
+sc-gapsys-scale-linpoint-q = 0.3
+sc-gapsys-sigma-lj = 0.3
separate-dhdl-file = yes
dhdl-derivatives = yes
dh_hist_size = 0
sc-r-power = 6
sc-sigma = 0.3
sc-coul = no
-sc-gapsys-scale-linpoint-LJ = 0.85
-sc-gapsys-scale-linpoint-Q = 0.3
-sc-gapsys-sigma-LJ = 0.3
+sc-gapsys-scale-linpoint-lj = 0.85
+sc-gapsys-scale-linpoint-q = 0.3
+sc-gapsys-sigma-lj = 0.3
separate-dhdl-file = yes
dhdl-derivatives = yes
dh_hist_size = 0
sc-r-power = 6
sc-sigma = 0.3
sc-coul = no
-sc-gapsys-scale-linpoint-LJ = 0.85
-sc-gapsys-scale-linpoint-Q = 0.3
-sc-gapsys-sigma-LJ = 0.3
+sc-gapsys-scale-linpoint-lj = 0.85
+sc-gapsys-scale-linpoint-q = 0.3
+sc-gapsys-sigma-lj = 0.3
separate-dhdl-file = yes
dhdl-derivatives = yes
dh_hist_size = 0
sc-r-power = 6
sc-sigma = 0.3
sc-coul = no
-sc-gapsys-scale-linpoint-LJ = 0.85
-sc-gapsys-scale-linpoint-Q = 0.3
-sc-gapsys-sigma-LJ = 0.3
+sc-gapsys-scale-linpoint-lj = 0.85
+sc-gapsys-scale-linpoint-q = 0.3
+sc-gapsys-sigma-lj = 0.3
separate-dhdl-file = yes
dhdl-derivatives = yes
dh_hist_size = 0
sc-r-power = 6
sc-sigma = 0.3
sc-coul = no
-sc-gapsys-scale-linpoint-LJ = 0.85
-sc-gapsys-scale-linpoint-Q = 0.3
-sc-gapsys-sigma-LJ = 0.3
+sc-gapsys-scale-linpoint-lj = 0.85
+sc-gapsys-scale-linpoint-q = 0.3
+sc-gapsys-sigma-lj = 0.3
separate-dhdl-file = yes
dhdl-derivatives = yes
dh_hist_size = 0
sc-r-power = 6
sc-sigma = 0.3
sc-coul = no
-sc-gapsys-scale-linpoint-LJ = 0.85
-sc-gapsys-scale-linpoint-Q = 0.3
-sc-gapsys-sigma-LJ = 0.3
+sc-gapsys-scale-linpoint-lj = 0.85
+sc-gapsys-scale-linpoint-q = 0.3
+sc-gapsys-sigma-lj = 0.3
separate-dhdl-file = yes
dhdl-derivatives = yes
dh_hist_size = 0
PS("separate-dhdl-file", enumValueToString(fep->separate_dhdl_file));
PS("dhdl-derivatives", enumValueToString(fep->dhdl_derivatives));
PS("sc-function", enumValueToString(fep->softcoreFunction));
- PR("sc-gapsys-scale-linpoint-LJ", fep->scGapsysScaleLinpointLJ);
- PR("sc-gapsys-scale-linpoint-Q", fep->scGapsysScaleLinpointQ);
- PR("sc-gapsys-sigma-LJ", fep->scGapsysSigmaLJ);
+ PR("sc-gapsys-scale-linpoint-lj", fep->scGapsysScaleLinpointLJ);
+ PR("sc-gapsys-scale-linpoint-q", fep->scGapsysScaleLinpointQ);
+ PR("sc-gapsys-sigma-lj", fep->scGapsysSigmaLJ);
};
static void pr_pull(FILE* fp, int indent, const pull_params_t& pull)