pair_style oxdna2/excv command

pair_style oxdna2/stk command

pair_style oxdna2/hbond command

pair_style oxdna2/xstk command

pair_style oxdna2/coaxstk command

pair_style oxdna2/dh command

Syntax

pair_style style1

pair_coeff * * style2 args
  • style1 = hybrid/overlay oxdna2/excv oxdna2/stk oxdna2/hbond oxdna2/xstk oxdna2/coaxstk oxdna2/dh
  • style2 = oxdna2/stk or oxdna2/dh
  • args = list of arguments for these two particular styles
oxdna2/stk args = T 6.0 0.4 0.9 0.32 0.6 1.3 0 0.8 0.9 0 0.95 0.9 0 0.95 2.0 0.65 2.0 0.65
  T = temperature (oxDNA units, 0.1 = 300 K)
oxdna2/dh args = T rhos qeff
  T = temperature (oxDNA units, 0.1 = 300 K)
  rhos = salt concentration (mole per litre)
  qeff = effective charge (elementary charges)

Examples

pair_style hybrid/overlay oxdna2/excv oxdna2/stk oxdna2/hbond oxdna2/xstk oxdna2/coaxstk oxdna2/dh
pair_coeff * * oxdna2/excv    2.0 0.7 0.675 2.0 0.515 0.5 2.0 0.33 0.32
pair_coeff * * oxdna2/stk     0.1 6.0 0.4 0.9 0.32 0.6 1.3 0 0.8 0.9 0 0.95 0.9 0 0.95 2.0 0.65 2.0 0.65
pair_coeff * * oxdna2/hbond   0.0 8.0 0.4 0.75 0.34 0.7 1.5 0 0.7 1.5 0 0.7 1.5 0 0.7 0.46 3.141592653589793 0.7 4.0 1.5707963267948966 0.45 4.0 1.5707963267948966 0.45
pair_coeff 1 4 oxdna2/hbond   1.0678 8.0 0.4 0.75 0.34 0.7 1.5 0 0.7 1.5 0 0.7 1.5 0 0.7 0.46 3.141592653589793 0.7 4.0 1.5707963267948966 0.45 4.0 1.5707963267948966 0.45
pair_coeff 2 3 oxdna2/hbond   1.0678 8.0 0.4 0.75 0.34 0.7 1.5 0 0.7 1.5 0 0.7 1.5 0 0.7 0.46 3.141592653589793 0.7 4.0 1.5707963267948966 0.45 4.0 1.5707963267948966 0.45
pair_coeff * * oxdna2/xstk    47.5 0.575 0.675 0.495 0.655 2.25 0.791592653589793 0.58 1.7 1.0 0.68 1.7 1.0 0.68 1.5 0 0.65 1.7 0.875 0.68 1.7 0.875 0.68
pair_coeff * * oxdna2/coaxstk 58.5 0.4 0.6 0.22 0.58 2.0 2.891592653589793 0.65 1.3 0 0.8 0.9 0 0.95 0.9 0 0.95 40.0 3.116592653589793
pair_coeff * * oxdna2/dh      0.1 1.0 0.815

Description

The oxdna2 pair styles compute the pairwise-additive parts of the oxDNA force field for coarse-grained modelling of DNA. The effective interaction between the nucleotides consists of potentials for the excluded volume interaction oxdna2/excv, the stacking oxdna2/stk, cross-stacking oxdna2/xstk and coaxial stacking interaction oxdna2/coaxstk, electrostatic Debye-Hueckel interaction oxdna2/dh as well as the hydrogen-bonding interaction oxdna2/hbond between complementary pairs of nucleotides on opposite strands.

The exact functional form of the pair styles is rather complex. The individual potentials consist of products of modulation factors, which themselves are constructed from a number of more basic potentials (Morse, Lennard-Jones, harmonic angle and distance) as well as quadratic smoothing and modulation terms. We refer to (Snodin) and the original oxDNA publications (Ouldridge-DPhil) and (Ouldridge) for a detailed description of the oxDNA2 force field.

Note

These pair styles have to be used together with the related oxDNA2 bond style oxdna2/fene for the connectivity of the phosphate backbone (see also documentation of bond_style oxdna2/fene). Almost all coefficients in the above example have to be kept fixed and cannot be changed without reparametrizing the entire model. Exceptions are the first coefficient after oxdna2/stk (T=0.1 in the above example) and the coefficients after oxdna2/dh (T=0.1, rhos=1.0, qeff=0.815 in the above example). When using a Langevin thermostat e.g. through fix langevin or fix nve/dotc/langevin the temperature coefficients have to be matched to the one used in the fix.

Example input and data files for DNA duplexes can be found in examples/USER/cgdna/examples/oxDNA/ and /oxDNA2/. A simple python setup tool which creates single straight or helical DNA strands, DNA duplexes or arrays of DNA duplexes can be found in examples/USER/cgdna/util/. A technical report with more information on the model, the structure of the input file, the setup tool and the performance of the LAMMPS-implementation of oxDNA can be found here.


Restrictions

These pair styles can only be used if LAMMPS was built with the USER-CGDNA package and the MOLECULE and ASPHERE package. See the Making LAMMPS section for more info on packages.