Merge branch 'master' into fix-compiler-warnings

doc/src/Eqs/fix_bond_react.tex

@@ -0,0 +1,9 @@
+\documentstyle[12pt]{article}
+\pagestyle{empty}
+\begin{document}
+
+\begin{eqnarray*}
+   k = AT^{n}e^{\frac{-E_{a}}{k_{B}T}}
+\end{eqnarray*}                           
+
+\end{document}

doc/src/PDF/pair_gayberne_extra.tex

@@ -48,8 +48,8 @@ $$ \mathbf{G}_{12} = \mathbf{A}_1^T \mathbf{S}_1^2 \mathbf{A}_1 +
 \mathbf{G}_2. $$
 
 Let the relative energy matrices $\mathbf{E}_i = \mbox{diag}
-(\epsilon_{ia}, \epsilon_{ib}, \epsilon_{ic})$ be given by 
-the relative well depths (dimensionless energy scales 
+(\epsilon_{ia}^{-1/\mu}, \epsilon_{ib}^{-1/\mu}, \epsilon_{ic}^{-1/\mu})$
+be given by the relative well depths (dimensionless energy scales 
 inversely proportional to the well-depths of the respective 
 orthogonal configurations of the interacting molecules). The 
 $\chi$ orientation-dependent energy based on the user-specified 
@@ -62,8 +62,8 @@ $$ \hat{\mathbf{r}}_{12} = { \mathbf{r}_{12} } / |\mathbf{r}_{12}|, $$
 
 and
 
-$$ \mathbf{B}_{12} = \mathbf{A}_1^T \mathbf{E}_1^2 \mathbf{A}_1 +
-\mathbf{A}_2^T \mathbf{E}_2^2 \mathbf{A}_2 = \mathbf{B}_1 +
+$$ \mathbf{B}_{12} = \mathbf{A}_1^T \mathbf{E}_1 \mathbf{A}_1 +
+\mathbf{A}_2^T \mathbf{E}_2 \mathbf{A}_2 = \mathbf{B}_1 +
 \mathbf{B}_2. $$
 
 Here, we use the distance of closest approach approximation given by the

doc/src/fix_bond_react.txt

@@ -266,7 +266,7 @@ either 'none' or 'charges.' Further details are provided in the
 discussion of the 'update_edges' keyword. The fourth optional section
 begins with the keyword 'Constraints' and lists additional criteria
 that must be satisfied in order for the reaction to occur. Currently,
-there are two types of constraints available, as discussed below.
+there are three types of constraints available, as discussed below.
 
 A sample map file is given below:
 
@@ -320,6 +320,27 @@ the central atom). Angles must be specified in degrees. This
 constraint can be used to enforce a certain orientation between
 reacting molecules.
 
+The constraint of type 'arrhenius' imposes an additional reaction
+probability according to the temperature-dependent Arrhenius equation:
+
+:c,image(Eqs/fix_bond_react.jpg)
+
+The Arrhenius constraint has the following syntax:
+
+arrhenius {A} {n} {E_a} {seed} :pre
+
+where 'arrhenius' is the required keyword, {A} is the pre-exponential
+factor, {n} is the exponent of the temperature dependence, {E_a} is
+the activation energy ("units"_units.html of energy), and {seed} is a
+random number seed. The temperature is defined as the instantaneous
+temperature averaged over all atoms in the reaction site, and is
+calculated in the same manner as for example
+"compute_temp_chunk"_compute_temp_chunk.html. Currently, there are no
+options for additional temperature averaging or velocity-biased
+temperature calculations. A uniform random number between 0 and 1 is
+generated using {seed}; if this number is less than the result of the
+Arrhenius equation above, the reaction is permitted to occur.
+
 Once a reaction site has been successfully identified, data structures
 within LAMMPS that store bond topology are updated to reflect the
 post-reacted molecule template. All force fields with fixed bonds,
@@ -395,7 +416,7 @@ will apply to all reactions.
 Computationally, each timestep this fix operates, it loops over
 neighbor lists (for bond-forming reactions) and computes distances
 between pairs of atoms in the list. It also communicates between
-neighboring processors to coordinate which bonds  are created and/or
+neighboring processors to coordinate which bonds are created and/or
 removed. All of these operations increase the cost of a timestep. Thus
 you should be cautious about invoking this fix too frequently.
 

doc/src/pair_oxdna.txt

@@ -23,7 +23,7 @@ style1 = {hybrid/overlay oxdna/excv oxdna/stk oxdna/hbond oxdna/xstk oxdna/coaxs
 style2 = {oxdna/excv} or {oxdna/stk} or {oxdna/hbond} or {oxdna/xstk} or {oxdna/coaxstk}
 args = list of arguments for these particular styles :ul
 
-  {oxdna/stk} args = seq T xi kappa 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
+  {oxdna/stk} args = seq T xi kappa 6.0 0.4 0.9 0.32 0.75 1.3 0 0.8 0.9 0 0.95 0.9 0 0.95 2.0 0.65 2.0 0.65
     seq = seqav (for average sequence stacking strength) or seqdep (for sequence-dependent stacking strength)
     T = temperature (oxDNA units, 0.1 = 300 K)
     xi = temperature-independent coefficient in stacking strength
@@ -36,7 +36,7 @@ args = list of arguments for these particular styles :ul
 
 pair_style hybrid/overlay oxdna/excv oxdna/stk oxdna/hbond oxdna/xstk oxdna/coaxstk
 pair_coeff * * oxdna/excv    2.0 0.7 0.675 2.0 0.515 0.5 2.0 0.33 0.32
-pair_coeff * * oxdna/stk     seqdep 0.1 1.3448 2.6568 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 * * oxdna/stk     seqdep 0.1 1.3448 2.6568 6.0 0.4 0.9 0.32 0.75 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 * * oxdna/hbond   seqdep 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 oxdna/hbond   seqdep 1.077 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 oxdna/hbond   seqdep 1.077 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
@@ -51,7 +51,10 @@ excluded volume interaction {oxdna/excv}, the stacking {oxdna/stk}, cross-stacki
 and coaxial stacking interaction {oxdna/coaxstk} as well
 as the hydrogen-bonding interaction {oxdna/hbond} between complementary pairs of nucleotides on
 opposite strands. Average sequence or sequence-dependent stacking and base-pairing strengths
-are supported "(Sulc)"_#Sulc1.
+are supported "(Sulc)"_#Sulc1. Quasi-unique base-pairing between nucleotides can be achieved by using 
+more complementary pairs of atom types like 5-8 and 6-7, 9-12 and 10-11, 13-16 and 14-15, etc. 
+This prevents the hybridization of in principle complementary bases within Ntypes/4 bases 
+up and down along the backbone.
 
 The exact functional form of the pair styles is rather complex.
 The individual potentials consist of products of modulation factors,

doc/src/pair_oxdna2.txt

@@ -24,7 +24,7 @@ style1 = {hybrid/overlay oxdna2/excv oxdna2/stk oxdna2/hbond oxdna2/xstk oxdna2/
 style2 = {oxdna2/excv} or {oxdna2/stk} or {oxdna2/hbond} or {oxdna2/xstk} or {oxdna2/coaxstk} or {oxdna2/dh}
 args = list of arguments for these particular styles :ul
 
-  {oxdna2/stk} args = seq T xi kappa 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
+  {oxdna2/stk} args = seq T xi kappa 6.0 0.4 0.9 0.32 0.75 1.3 0 0.8 0.9 0 0.95 0.9 0 0.95 2.0 0.65 2.0 0.65
     seq = seqav (for average sequence stacking strength) or seqdep (for sequence-dependent stacking strength)
     T = temperature (oxDNA units, 0.1 = 300 K)
     xi = temperature-independent coefficient in stacking strength
@@ -41,7 +41,7 @@ args = list of arguments for these particular styles :ul
 
 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     seqdep 0.1 1.3523 2.6717 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/stk     seqdep 0.1 1.3523 2.6717 6.0 0.4 0.9 0.32 0.75 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   seqdep 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   seqdep 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   seqdep 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
@@ -57,7 +57,10 @@ excluded volume interaction {oxdna2/excv}, the stacking {oxdna2/stk}, cross-stac
 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. Average sequence or sequence-dependent stacking and base-pairing strengths
-are supported "(Sulc)"_#Sulc2.
+are supported "(Sulc)"_#Sulc2. Quasi-unique base-pairing between nucleotides can be achieved by using 
+more complementary pairs of atom types like 5-8 and 6-7, 9-12 and 10-11, 13-16 and 14-15, etc. 
+This prevents the hybridization of in principle complementary bases within Ntypes/4 bases 
+up and down along the backbone.
 
 The exact functional form of the pair styles is rather complex.
 The individual potentials consist of products of modulation factors,