Section containing the all coarse-graining options The following table contains the input options for CSG,

Property Name

Default Value

Description

fmatch

Force matching options

fmatch.constrainedLS

boolean variable: false - simple least squares, true -
constrained least squares. For details see the VOTCA paper.
Practically, both algorithms give the same results, but
simple least squares is faster. If you are a mathematician
and you think that a spline can only then be called a
spline if it has continuous first and second derivatives,
use constrained least squares.

fmatch.dist

1e-5

Accuracy for evaluating the difference in bead positions.
Default is 1e-5

fmatch.frames_per_block

number of frames, being used for block averaging. Atomistic
trajectory, specified with –trj option, is divided into
blocks and the force matching equations are solved separately
for each block. Coarse-grained force-field, which one
gets on the output is averaged over those blocks.

nbsearch

grid

Grid search algorithm, simple (N square search) or grid

bonded

Interaction specific option for bonded interactions, see
the cg.non-bonded section for all options

bonded.dlpoly

bonded.dlpoly.header

Header of the interaction in dlpoly TABBND or TABANG file.
The header should be a unique set of the interaction-site
names, and these should match the corresponding names
specified in the mapping file.

bonded.name

Name of the bonded interaction. The name can be arbitrary
but should be unique. For bonded interactions, this should
match the name specified in the mapping file.

bonded.periodic

0

set to 1 when calculating bond dihedral potentials with
csg_fmatch -> enforces periodicity of potential. (default
is 0)

non-bonded

Interaction specific option for non-bonded interactions

non-bonded.dlpoly

non-bonded.dlpoly.header

Header of the interaction in dlpoly TABLE file. The header
should be a unique pair of the interaction-site names,
and these should match the corresponding names specified
in the mapping file.

non-bonded.name

Name of the interaction. The name can be arbitrary but
should be unique. For bonded interactions, this should
match the name specified in the mapping file.

non-bonded.type1

Bead type 1 of non-bonded interaction.

non-bonded.type2

Bead type 2 of non-bonded interaction.

non-bonded.bondtype

Internal alias for “non-bonded” or “bonded”, set automatically

non-bonded.min

Lower bound of interval for potential table in which calculations
are performed. Should be set based on reference distributions.

non-bonded.max

Upper bound of interval for potential table in which calculations
are performed. Should be set based on reference distributions.

non-bonded.step

Step size of interval for potential table in which calculations
are performed. If step site is too small, lots of statistics
is needed ( long runs ). If it’s too big, features in
the distribtuion/potentials might get lost.

non-bonded.fmatch

Force matching options

non-bonded.fmatch.min

Minimum value of interval for distribution sampled in
atomistic MD simulation. One can get this number by looking
at the distribution function for this interaction. For
non-bonded interactions it’s the distance to the rdf start.
For CG bonds and angles the variable has the similar meaning
( note, that for angles it is specified in radians ).

non-bonded.fmatch.max

Maximum value of interval for distribution sampled in
atomistic MD simulation. One can get this number by looking
at the distribution function for this interaction. For
non-bonded interactions it’s the cut-off of the interaction.

non-bonded.fmatch.step

grid spacing for the spline, which represents the interaction.
This parameter should not be too big, otherwise you might
lose some features of the interaction potential, and not
too small either, otherwise you will have unsampled bins
which result in an ill-defined equation system and NaNs
in the output.

non-bonded.fmatch.out_step

Grid spacing for the output grid. Normally, one wants
to have this parameter smaller than fmatch.step, to have
a smooth curve, without additional spline interpolation.
As a rule of thumb we normally use fmatch.out_step which
is approximately 5 times smaller than fmatch.step.

non-bonded.re

Relative entropy options

non-bonded.re.function

Functional form for the potential. Available functional
forms: lj126 (Lennard-Jones 12-6), ljg (Lennard-Jones
12-6 plus Gaussian), and cbspl (uniform cubic B-splines).

non-bonded.re.cbspl

options specific to cbspl function form

non-bonded.re.cbspl.nknots

Number of knot values to be used for the cbspl functional
form. Uniform grid size of the CBSPL depends on this parameter;
for fixed potential range more the nknots smaller the
grid spacing. Make sure grid spacing is sufficiently large
and enough CG simulation steps are performed such that
the bins at distance greater than the minimum distance
are sampled sufficiently otherwise ill-defined system
of equation would give NaNs in the output.