Merge remote-tracking branch 'lammps-ro/master' into lammps-icms
|
Before Width: | Height: | Size: 599 KiB After Width: | Height: | Size: 9.0 KiB |
9
doc/Eqs/polymorphic1.tex
Normal file
@ -0,0 +1,9 @@
|
||||
\documentclass[12pt]{article}
|
||||
|
||||
\begin{document}
|
||||
|
||||
$$
|
||||
E=\frac{1}{2}\sum_{i=1}^{i=N}\sum_{j=1}^{j=N}\left[\left(1-\delta_{ij}\right)\cdot U_{IJ}\left(r_{ij}\right)-\left(1-\eta_{ij}\right)\cdot F_{IJ}\left(r_{ij}\right)\cdot V_{IJ}\left(r_{ij}\right)\right]
|
||||
$$
|
||||
|
||||
\end{document}
|
||||
|
Before Width: | Height: | Size: 610 KiB After Width: | Height: | Size: 7.0 KiB |
10
doc/Eqs/polymorphic2.tex
Normal file
@ -0,0 +1,10 @@
|
||||
\documentclass[12pt]{article}
|
||||
|
||||
\begin{document}
|
||||
|
||||
$$
|
||||
X_{ij}=\sum_{k=i_1,k\neq i,j}^{i_N}W_{IK}\left(r_{ik}\right)\cdot G_{JIK}\left(\theta_{jik}\right)\cdot P_{IK}\left(\Delta r_{jik}\right)
|
||||
\label{X_eq2}
|
||||
$$
|
||||
|
||||
\end{document}
|
||||
|
Before Width: | Height: | Size: 572 KiB After Width: | Height: | Size: 2.4 KiB |
10
doc/Eqs/polymorphic3.tex
Normal file
@ -0,0 +1,10 @@
|
||||
\documentclass[12pt]{article}
|
||||
|
||||
\begin{document}
|
||||
|
||||
$$
|
||||
\Delta r_{jik}=r_{ij}-\xi_{IJ}\cdot r_{ik}
|
||||
\label{Dr_eq3}
|
||||
$$
|
||||
|
||||
\end{document}
|
||||
|
Before Width: | Height: | Size: 638 KiB After Width: | Height: | Size: 26 KiB |
17
doc/Eqs/polymorphic4.tex
Normal file
@ -0,0 +1,17 @@
|
||||
\documentclass[12pt]{article}
|
||||
|
||||
\begin{document}
|
||||
|
||||
\begin{eqnarray*}
|
||||
\left\{\begin{array}{l}
|
||||
\eta_{ij}=\delta_{ij},\xi_{IJ}=0 \\
|
||||
U_{IJ}\left(r\right)=A_{IJ}\cdot\epsilon_{IJ}\cdot \left(\frac{\sigma_{IJ}}{r}\right)^q\cdot \left[B_{IJ}\cdot \left(\frac{\sigma_{IJ}}{r}\right)^{p-q}-1\right]\cdot exp\left(\frac{\sigma_{IJ}}{r-a_{IJ}\cdot \sigma_{IJ}}\right) \\
|
||||
V_{IJ}\left(r\right)=\sqrt{\lambda_{IJ}\cdot \epsilon_{IJ}}\cdot exp\left(\frac{\gamma_{IJ}\cdot \sigma_{IJ}}{r-a_{IJ}\cdot \sigma_{IJ}}\right) \\
|
||||
F_{IJ}\left(X\right)=-X \\
|
||||
P_{IJ}\left(\Delta r\right)=1 \\
|
||||
W_{IJ}\left(r\right)=\sqrt{\lambda_{IJ}\cdot \epsilon_{IJ}}\cdot exp\left(\frac{\gamma_{IJ}\cdot \sigma_{IJ}}{r-a_{IJ}\cdot \sigma_{IJ}}\right) \\
|
||||
G_{JIK}\left(\theta\right)=\left(cos\theta+\frac{1}{3}\right)^2
|
||||
\end{array}\right.
|
||||
\end{eqnarray*}
|
||||
|
||||
\end{document}
|
||||
|
Before Width: | Height: | Size: 675 KiB After Width: | Height: | Size: 29 KiB |
17
doc/Eqs/polymorphic5.tex
Normal file
@ -0,0 +1,17 @@
|
||||
\documentclass[12pt]{article}
|
||||
|
||||
\begin{document}
|
||||
|
||||
\begin{eqnarray*}
|
||||
\left\{\begin{array}{l}
|
||||
\eta_{ij}=\delta_{ij},\xi_{IJ}=1 \\
|
||||
U_{IJ}\left(r\right)=\frac{D_{e,IJ}}{S_{IJ}-1}\cdot exp\left[-\beta_{IJ}\sqrt{2S_{IJ}\left(r-r_{e,IJ}\right)}\right]\cdot f_{c,IJ}\left(r\right) \\
|
||||
V_{IJ}\left(r\right)=\frac{S_{IJ}\cdot D_{e,IJ}}{S_{IJ}-1}\cdot exp\left[-\beta_{IJ}\sqrt{\frac{2}{S_{IJ}}\left(r-r_{e,IJ}\right)}\right]\cdot f_{c,IJ}\left(r\right) \\
|
||||
F_{IJ}\left(X\right)=\left(1+X\right)^{-\frac{1}{2}} \\
|
||||
P_{IJ}\left(\Delta r\right)=exp\left(2\mu_{IK}\cdot \Delta r\right) \\
|
||||
W_{IJ}\left(r\right)=f_{c,IK}\left(r\right) \\
|
||||
G_{JIK}\left(\theta\right)=\gamma_{IK}\left[1+\frac{c_{IK}^2}{d_{IK}^2}-\frac{c_{IK}^2}{d_{IK}^2+\left(h_{IK}+cos\theta\right)^2}\right]
|
||||
\end{array}\right.
|
||||
\end{eqnarray*}
|
||||
|
||||
\end{document}
|
||||
|
Before Width: | Height: | Size: 604 KiB After Width: | Height: | Size: 8.9 KiB |
13
doc/Eqs/polymorphic6.tex
Normal file
@ -0,0 +1,13 @@
|
||||
\documentclass[12pt]{article}
|
||||
|
||||
\begin{document}
|
||||
|
||||
\begin{eqnarray*}
|
||||
f_{c,IJ}=\left\{\begin{array}{lr}
|
||||
1, & r\leq r_{s,IJ} \\
|
||||
\frac{1}{2}+\frac{1}{2} cos \left[\frac{\pi \left(r-r_{s,IJ}\right)}{r_{c,IJ}-r_{s,IJ}}\right], & r_{s,IJ}<r<r_{c,IJ} \\
|
||||
0, & r \geq r_{c,IJ} \\
|
||||
\end{array}\right.
|
||||
\end{eqnarray*}
|
||||
|
||||
\end{document}
|
||||
|
Before Width: | Height: | Size: 651 KiB After Width: | Height: | Size: 57 KiB |
25
doc/Eqs/polymorphic7.tex
Normal file
@ -0,0 +1,25 @@
|
||||
\documentclass[12pt]{article}
|
||||
|
||||
\begin{document}
|
||||
|
||||
\begin{eqnarray*}
|
||||
\left\{\begin{array}{l}
|
||||
\eta_{ij}=\delta_{ij},\xi_{IJ}=1 \\
|
||||
U_{IJ}\left(r\right)=\left\{\begin{array}{lr}
|
||||
A_{IJ}\cdot exp\left(-\lambda_{1,IJ}\cdot r\right)\cdot f_{c,IJ}\left(r\right), & r\leq r_{s,1,IJ} \\
|
||||
A_{IJ}\cdot exp\left(-\lambda_{1,IJ}\cdot r\right)\cdot f_{c,IJ}\left(r\right)\cdot f_{c,1,IJ}\left(r\right), & r_{s,1,IJ}<r<r_{c,1,IJ} \\
|
||||
0, & r\ge r_{c,1,IJ}
|
||||
\end{array}\right. \\
|
||||
V_{IJ}\left(r\right)=\left\{\begin{array}{lr}
|
||||
B_{IJ} \cdot exp\left(-\lambda_{2,IJ}\cdot r\right)\cdot f_{c,IJ}\left(r\right), & r\le r_{s,1,IJ} \\
|
||||
B_{IJ} \cdot exp\left(-\lambda_{2,IJ}\cdot r\right)\cdot f_{c,IJ}\left(r\right)+A_{IJ}\cdot exp\left(-\lambda_{1,IJ}\cdot r\right)\cdot & \\ ~~~~~~ f_{c,IJ}\left(r\right)\cdot \left[1-f_{c,1,IJ}\left(r\right)\right], & r_{s,1,IJ}<r<r_{c,1,IJ} \\
|
||||
B_{IJ} \cdot exp\left(-\lambda_{2,IJ}\cdot r\right)\cdot f_{c,IJ}\left(r\right)+A_{IJ}\cdot exp\left(-\lambda_{1,IJ}\cdot r\right)\cdot & \\ ~~~~~~ f_{c,IJ}\left(r\right) & r \ge r_{c,1,IJ}
|
||||
\end{array}\right. \\
|
||||
F_{IJ}\left(X\right)=\left[1+\left(\beta_{IJ}\cdot X\right)^{n_{IJ}}\right]^{-\frac{1}{2n_{IJ}}} \\
|
||||
P_{IJ}\left(\Delta r\right)=exp\left(\lambda_{3,IK}\cdot \Delta r^3\right) \\
|
||||
W_{IJ}\left(r\right)=f_{c,IK}\left(r\right) \\
|
||||
G_{JIK}\left(\theta\right)=1+\frac{c_{IK}^2}{d_{IK}^2}-\frac{c_{IK}^2}{d_{IK}^2+\left(h_{IK}+cos\theta\right)^2}
|
||||
\end{array}\right.
|
||||
\end{eqnarray*}
|
||||
|
||||
\end{document}
|
||||
|
Before Width: | Height: | Size: 594 KiB After Width: | Height: | Size: 10 KiB |
13
doc/Eqs/polymorphic8.tex
Normal file
@ -0,0 +1,13 @@
|
||||
\documentclass[12pt]{article}
|
||||
|
||||
\begin{document}
|
||||
|
||||
\begin{eqnarray*}
|
||||
f_{c,1,IJ}=\left\{\begin{array}{lr}
|
||||
1, & r\leq r_{s,1,IJ} \\
|
||||
\frac{1}{2}+\frac{1}{2} cos \left[\frac{\pi \left(r-r_{s,1,IJ}\right)}{r_{c,1,IJ}-r_{s,1,IJ}}\right], & r_{s,1,IJ}<r<r_{c,1,IJ} \\
|
||||
0, & r \geq r_{c,1,IJ} \\
|
||||
\end{array}\right.
|
||||
\end{eqnarray*}
|
||||
|
||||
\end{document}
|
||||
|
Before Width: | Height: | Size: 661 KiB After Width: | Height: | Size: 13 KiB |
17
doc/Eqs/polymorphic9.tex
Normal file
@ -0,0 +1,17 @@
|
||||
\documentclass[12pt]{article}
|
||||
|
||||
\begin{document}
|
||||
|
||||
\begin{eqnarray*}
|
||||
\left\{\begin{array}{l}
|
||||
\eta_{ij}=1-\delta_{ij},\xi_{IJ}=0 \\
|
||||
U_{IJ}\left(r\right)=\phi_{IJ}\left(r\right) \\
|
||||
V_{IJ}\left(r\right)=1 \\
|
||||
F_{II}\left(X\right)=-2F_I\left(X\right) \\
|
||||
P_{IJ}\left(\Delta r\right)=1 \\
|
||||
W_{IJ}\left(r\right)=f_{K}\left(r\right) \\
|
||||
G_{JIK}\left(\theta\right)=1
|
||||
\end{array}\right.
|
||||
\end{eqnarray*}
|
||||
|
||||
\end{document}
|
||||
@ -860,7 +860,7 @@ variables that have to match the VMD installation on your system.</p>
|
||||
source code for individual scripts for comments about what they do.</p>
|
||||
<p>Here are screenshots of the vizplotgui_tool.py script in action for
|
||||
different visualization package options. Click to see larger images:</p>
|
||||
<a data-lightbox="group-c530f122-6019-45af-8a35-e7cdc70ddb55"
|
||||
<a data-lightbox="group-default"
|
||||
href="_images/screenshot_gl.jpg"
|
||||
class=""
|
||||
title=""
|
||||
@ -870,7 +870,7 @@ different visualization package options. Click to see larger images:</p>
|
||||
width="25%"
|
||||
height="auto"
|
||||
alt=""/>
|
||||
</a><a data-lightbox="group-f712591c-ea31-483b-86da-3905db788920"
|
||||
</a><a data-lightbox="group-default"
|
||||
href="_images/screenshot_atomeye.jpg"
|
||||
class=""
|
||||
title=""
|
||||
@ -880,7 +880,7 @@ different visualization package options. Click to see larger images:</p>
|
||||
width="25%"
|
||||
height="auto"
|
||||
alt=""/>
|
||||
</a><a data-lightbox="group-4a2be907-5b47-4889-89de-f686acc6e915"
|
||||
</a><a data-lightbox="group-default"
|
||||
href="_images/screenshot_pymol.jpg"
|
||||
class=""
|
||||
title=""
|
||||
@ -890,7 +890,7 @@ different visualization package options. Click to see larger images:</p>
|
||||
width="25%"
|
||||
height="auto"
|
||||
alt=""/>
|
||||
</a><a data-lightbox="group-e62d1b66-7f66-44c5-ba32-c4396e2529ed"
|
||||
</a><a data-lightbox="group-default"
|
||||
href="_images/screenshot_vmd.jpg"
|
||||
class=""
|
||||
title=""
|
||||
|
||||
|
Before Width: | Height: | Size: 599 KiB After Width: | Height: | Size: 9.0 KiB |
|
Before Width: | Height: | Size: 610 KiB After Width: | Height: | Size: 7.0 KiB |
|
Before Width: | Height: | Size: 638 KiB After Width: | Height: | Size: 26 KiB |
|
Before Width: | Height: | Size: 675 KiB After Width: | Height: | Size: 29 KiB |
|
Before Width: | Height: | Size: 604 KiB After Width: | Height: | Size: 8.9 KiB |
|
Before Width: | Height: | Size: 651 KiB After Width: | Height: | Size: 57 KiB |
|
Before Width: | Height: | Size: 594 KiB After Width: | Height: | Size: 10 KiB |
|
Before Width: | Height: | Size: 661 KiB After Width: | Height: | Size: 13 KiB |
2
doc/_static/css/theme.css
vendored
@ -3743,7 +3743,7 @@ hr {
|
||||
display: block;
|
||||
height: 1px;
|
||||
border: 0;
|
||||
border-top: 1px solid #e1e4e5;
|
||||
border-top: 2px solid #e1e4e5;
|
||||
margin: 24px 0;
|
||||
padding: 0;
|
||||
}
|
||||
|
||||
@ -137,11 +137,11 @@
|
||||
</ul>
|
||||
<pre class="literal-block">
|
||||
<em>id</em> value = <em>yes</em> or <em>no</em>
|
||||
<em>map</em> value = <em>array</em> or <em>hash</em>
|
||||
<em>first</em> value = group-ID = group whose atoms will appear first in internal atom lists
|
||||
<em>sort</em> values = Nfreq binsize
|
||||
Nfreq = sort atoms spatially every this many time steps
|
||||
binsize = bin size for spatial sorting (distance units)
|
||||
<em>map</em> value = <em>array</em> or <em>hash</em>
|
||||
<em>first</em> value = group-ID = group whose atoms will appear first in internal atom lists
|
||||
<em>sort</em> values = Nfreq binsize
|
||||
Nfreq = sort atoms spatially every this many time steps
|
||||
binsize = bin size for spatial sorting (distance units)
|
||||
</pre>
|
||||
</div>
|
||||
<div class="section" id="examples">
|
||||
|
||||
@ -136,13 +136,13 @@
|
||||
</ul>
|
||||
<pre class="literal-block">
|
||||
args = none for any style except <em>body</em> and <em>hybrid</em>
|
||||
<em>body</em> args = bstyle bstyle-args
|
||||
bstyle = style of body particles
|
||||
bstyle-args = additional arguments specific to the bstyle
|
||||
see the <a class="reference internal" href="body.html"><em>body</em></a> doc page for details
|
||||
<em>template</em> args = template-ID
|
||||
template-ID = ID of molecule template specified in a separate <a class="reference internal" href="molecule.html"><em>molecule</em></a> command
|
||||
<em>hybrid</em> args = list of one or more sub-styles, each with their args
|
||||
<em>body</em> args = bstyle bstyle-args
|
||||
bstyle = style of body particles
|
||||
bstyle-args = additional arguments specific to the bstyle
|
||||
see the <a class="reference internal" href="body.html"><em>body</em></a> doc page for details
|
||||
<em>template</em> args = template-ID
|
||||
template-ID = ID of molecule template specified in a separate <a class="reference internal" href="molecule.html"><em>molecule</em></a> command
|
||||
<em>hybrid</em> args = list of one or more sub-styles, each with their args
|
||||
</pre>
|
||||
<ul class="simple">
|
||||
<li>accelerated styles (with same args) = <em>angle/cuda</em> or <em>angle/kk</em> or <em>atomic/cuda</em> or <em>atomic/kk</em> or <em>bond/kk</em> or <em>charge/cuda</em> or <em>charge/kk</em> or <em>full/cuda</em> or <em>full/kk</em> or <em>molecular/kk</em></li>
|
||||
|
||||
@ -138,22 +138,22 @@
|
||||
</ul>
|
||||
<pre class="literal-block">
|
||||
<em>x</em> args = <em>uniform</em> or Px-1 numbers between 0 and 1
|
||||
<em>uniform</em> = evenly spaced cuts between processors in x dimension
|
||||
numbers = Px-1 ascending values between 0 and 1, Px - # of processors in x dimension
|
||||
<em>x</em> can be specified together with <em>y</em> or <em>z</em>
|
||||
<em>y</em> args = <em>uniform</em> or Py-1 numbers between 0 and 1
|
||||
<em>uniform</em> = evenly spaced cuts between processors in y dimension
|
||||
numbers = Py-1 ascending values between 0 and 1, Py - # of processors in y dimension
|
||||
<em>y</em> can be specified together with <em>x</em> or <em>z</em>
|
||||
<em>z</em> args = <em>uniform</em> or Pz-1 numbers between 0 and 1
|
||||
<em>uniform</em> = evenly spaced cuts between processors in z dimension
|
||||
numbers = Pz-1 ascending values between 0 and 1, Pz - # of processors in z dimension
|
||||
<em>z</em> can be specified together with <em>x</em> or <em>y</em>
|
||||
<em>shift</em> args = dimstr Niter stopthresh
|
||||
dimstr = sequence of letters containing "x" or "y" or "z", each not more than once
|
||||
Niter = # of times to iterate within each dimension of dimstr sequence
|
||||
stopthresh = stop balancing when this imbalance threshhold is reached
|
||||
<em>rcb</em> args = none
|
||||
<em>uniform</em> = evenly spaced cuts between processors in x dimension
|
||||
numbers = Px-1 ascending values between 0 and 1, Px - # of processors in x dimension
|
||||
<em>x</em> can be specified together with <em>y</em> or <em>z</em>
|
||||
<em>y</em> args = <em>uniform</em> or Py-1 numbers between 0 and 1
|
||||
<em>uniform</em> = evenly spaced cuts between processors in y dimension
|
||||
numbers = Py-1 ascending values between 0 and 1, Py - # of processors in y dimension
|
||||
<em>y</em> can be specified together with <em>x</em> or <em>z</em>
|
||||
<em>z</em> args = <em>uniform</em> or Pz-1 numbers between 0 and 1
|
||||
<em>uniform</em> = evenly spaced cuts between processors in z dimension
|
||||
numbers = Pz-1 ascending values between 0 and 1, Pz - # of processors in z dimension
|
||||
<em>z</em> can be specified together with <em>x</em> or <em>y</em>
|
||||
<em>shift</em> args = dimstr Niter stopthresh
|
||||
dimstr = sequence of letters containing "x" or "y" or "z", each not more than once
|
||||
Niter = # of times to iterate within each dimension of dimstr sequence
|
||||
stopthresh = stop balancing when this imbalance threshhold is reached
|
||||
<em>rcb</em> args = none
|
||||
</pre>
|
||||
<ul class="simple">
|
||||
<li>zero or more keyword/value pairs may be appended</li>
|
||||
@ -161,7 +161,7 @@
|
||||
</ul>
|
||||
<pre class="literal-block">
|
||||
<em>out</em> value = filename
|
||||
filename = write each processor's sub-domain to a file
|
||||
filename = write each processor's sub-domain to a file
|
||||
</pre>
|
||||
</div>
|
||||
<div class="section" id="examples">
|
||||
@ -265,7 +265,7 @@ the processor that owns them. The leftmost diagram is the default
|
||||
partitioning of the simulation box across processors (one sub-box for
|
||||
each of 16 processors); the middle diagram is after a “grid” method
|
||||
has been applied.</p>
|
||||
<a data-lightbox="group-1ea1982f-2af8-42b4-be30-2e7feb8e77f5"
|
||||
<a data-lightbox="group-default"
|
||||
href="_images/balance_uniform.jpg"
|
||||
class=""
|
||||
title=""
|
||||
@ -275,7 +275,7 @@ has been applied.</p>
|
||||
width="25%"
|
||||
height="auto"
|
||||
alt=""/>
|
||||
</a><a data-lightbox="group-a4b3cf56-edc0-41c1-a22e-57acb6d303e2"
|
||||
</a><a data-lightbox="group-default"
|
||||
href="_images/balance_nonuniform.jpg"
|
||||
class=""
|
||||
title=""
|
||||
@ -285,7 +285,7 @@ has been applied.</p>
|
||||
width="25%"
|
||||
height="auto"
|
||||
alt=""/>
|
||||
</a><a data-lightbox="group-c0259437-bded-4589-b518-4425480db921"
|
||||
</a><a data-lightbox="group-default"
|
||||
href="_images/balance_rcb.jpg"
|
||||
class=""
|
||||
title=""
|
||||
|
||||
@ -136,7 +136,7 @@
|
||||
</ul>
|
||||
<pre class="literal-block">
|
||||
args = none for any style except <em>hybrid</em>
|
||||
<em>hybrid</em> args = list of one or more styles
|
||||
<em>hybrid</em> args = list of one or more styles
|
||||
</pre>
|
||||
</div>
|
||||
<div class="section" id="examples">
|
||||
|
||||
@ -136,9 +136,9 @@
|
||||
</ul>
|
||||
<pre class="literal-block">
|
||||
<em>p</em> is periodic
|
||||
<em>f</em> is non-periodic and fixed
|
||||
<em>s</em> is non-periodic and shrink-wrapped
|
||||
<em>m</em> is non-periodic and shrink-wrapped with a minimum value
|
||||
<em>f</em> is non-periodic and fixed
|
||||
<em>s</em> is non-periodic and shrink-wrapped
|
||||
<em>m</em> is non-periodic and shrink-wrapped with a minimum value
|
||||
</pre>
|
||||
</div>
|
||||
<div class="section" id="examples">
|
||||
|
||||
@ -169,8 +169,8 @@ parameter = <em>x</em> or <em>y</em> or <em>z</em> or <em>xy</em> or <em>xz</em>
|
||||
</ul>
|
||||
<pre class="literal-block">
|
||||
<em>units</em> value = <em>lattice</em> or <em>box</em>
|
||||
lattice = distances are defined in lattice units
|
||||
box = distances are defined in simulation box units
|
||||
lattice = distances are defined in lattice units
|
||||
box = distances are defined in simulation box units
|
||||
</pre>
|
||||
</div>
|
||||
<div class="section" id="examples">
|
||||
|
||||
@ -137,9 +137,9 @@
|
||||
</ul>
|
||||
<pre class="literal-block">
|
||||
<em>mode</em> value = <em>single</em> or <em>multi</em> = communicate atoms within a single or multiple distances
|
||||
<em>cutoff</em> value = Rcut (distance units) = communicate atoms from this far away
|
||||
<em>group</em> value = group-ID = only communicate atoms in the group
|
||||
<em>vel</em> value = <em>yes</em> or <em>no</em> = do or do not communicate velocity info with ghost atoms
|
||||
<em>cutoff</em> value = Rcut (distance units) = communicate atoms from this far away
|
||||
<em>group</em> value = group-ID = only communicate atoms in the group
|
||||
<em>vel</em> value = <em>yes</em> or <em>no</em> = do or do not communicate velocity info with ghost atoms
|
||||
</pre>
|
||||
</div>
|
||||
<div class="section" id="examples">
|
||||
|
||||
@ -139,7 +139,7 @@
|
||||
</ul>
|
||||
<pre class="literal-block">
|
||||
<em>theta</em> = tabulate angles
|
||||
<em>eng</em> = tabulate angle energies
|
||||
<em>eng</em> = tabulate angle energies
|
||||
</pre>
|
||||
</div>
|
||||
<div class="section" id="examples">
|
||||
|
||||
@ -139,7 +139,7 @@
|
||||
</ul>
|
||||
<pre class="literal-block">
|
||||
<em>type</em> = atom type of the body particle
|
||||
<em>integer</em> = 1,2,3,etc = index of fields defined by body style
|
||||
<em>integer</em> = 1,2,3,etc = index of fields defined by body style
|
||||
</pre>
|
||||
</div>
|
||||
<div class="section" id="examples">
|
||||
|
||||
@ -139,8 +139,8 @@
|
||||
</ul>
|
||||
<pre class="literal-block">
|
||||
<em>dist</em> = bond distance
|
||||
<em>eng</em> = bond energy
|
||||
<em>force</em> = bond force
|
||||
<em>eng</em> = bond energy
|
||||
<em>force</em> = bond force
|
||||
</pre>
|
||||
</div>
|
||||
<div class="section" id="examples">
|
||||
|
||||
@ -164,30 +164,30 @@ style = <em>bin/1d</em> or <em>bin/2d</em> or <em>bin/3d</em> or <em>type</em> o
|
||||
</ul>
|
||||
<pre class="literal-block">
|
||||
<em>region</em> value = region-ID
|
||||
region-ID = ID of region atoms must be in to be part of a chunk
|
||||
<em>nchunk</em> value = <em>once</em> or <em>every</em>
|
||||
once = only compute the number of chunks once
|
||||
every = re-compute the number of chunks whenever invoked
|
||||
<em>limit</em> values = 0 or Nc max or Nc exact
|
||||
0 = no limit on the number of chunks
|
||||
Nc max = limit number of chunks to be <= Nc
|
||||
Nc exact = set number of chunks to exactly Nc
|
||||
<em>ids</em> value = <em>once</em> or <em>nfreq</em> or <em>every</em>
|
||||
once = assign chunk IDs to atoms only once, they persist thereafter
|
||||
nfreq = assign chunk IDs to atoms only once every Nfreq steps (if invoked by <a class="reference internal" href="fix_ave_chunk.html"><em>fix ave/chunk</em></a> which sets Nfreq)
|
||||
every = assign chunk IDs to atoms whenever invoked
|
||||
<em>compress</em> value = <em>yes</em> or <em>no</em>
|
||||
yes = compress chunk IDs to eliminate IDs with no atoms
|
||||
no = do not compress chunk IDs even if some IDs have no atoms
|
||||
<em>discard</em> value = <em>yes</em> or <em>no</em> or <em>mixed</em>
|
||||
yes = discard atoms with out-of-range chunk IDs by assigning a chunk ID = 0
|
||||
no = keep atoms with out-of-range chunk IDs by assigning a valid chunk ID
|
||||
mixed = keep or discard such atoms according to spatial binning rule
|
||||
<em>bound</em> values = x/y/z lo hi
|
||||
x/y/z = <em>x</em> or <em>y</em> or <em>z</em> to bound sptial bins in this dimension
|
||||
lo = <em>lower</em> or coordinate value (distance units)
|
||||
hi = <em>upper</em> or coordinate value (distance units)
|
||||
<em>units</em> value = <em>box</em> or <em>lattice</em> or <em>reduced</em>
|
||||
region-ID = ID of region atoms must be in to be part of a chunk
|
||||
<em>nchunk</em> value = <em>once</em> or <em>every</em>
|
||||
once = only compute the number of chunks once
|
||||
every = re-compute the number of chunks whenever invoked
|
||||
<em>limit</em> values = 0 or Nc max or Nc exact
|
||||
0 = no limit on the number of chunks
|
||||
Nc max = limit number of chunks to be <= Nc
|
||||
Nc exact = set number of chunks to exactly Nc
|
||||
<em>ids</em> value = <em>once</em> or <em>nfreq</em> or <em>every</em>
|
||||
once = assign chunk IDs to atoms only once, they persist thereafter
|
||||
nfreq = assign chunk IDs to atoms only once every Nfreq steps (if invoked by <a class="reference internal" href="fix_ave_chunk.html"><em>fix ave/chunk</em></a> which sets Nfreq)
|
||||
every = assign chunk IDs to atoms whenever invoked
|
||||
<em>compress</em> value = <em>yes</em> or <em>no</em>
|
||||
yes = compress chunk IDs to eliminate IDs with no atoms
|
||||
no = do not compress chunk IDs even if some IDs have no atoms
|
||||
<em>discard</em> value = <em>yes</em> or <em>no</em> or <em>mixed</em>
|
||||
yes = discard atoms with out-of-range chunk IDs by assigning a chunk ID = 0
|
||||
no = keep atoms with out-of-range chunk IDs by assigning a valid chunk ID
|
||||
mixed = keep or discard such atoms according to spatial binning rule
|
||||
<em>bound</em> values = x/y/z lo hi
|
||||
x/y/z = <em>x</em> or <em>y</em> or <em>z</em> to bound sptial bins in this dimension
|
||||
lo = <em>lower</em> or coordinate value (distance units)
|
||||
hi = <em>upper</em> or coordinate value (distance units)
|
||||
<em>units</em> value = <em>box</em> or <em>lattice</em> or <em>reduced</em>
|
||||
</pre>
|
||||
</div>
|
||||
<div class="section" id="examples">
|
||||
|
||||
@ -140,14 +140,14 @@
|
||||
</ul>
|
||||
<pre class="literal-block">
|
||||
<em>pair</em> args = pstyle pparam I J v_delta
|
||||
pstyle = pair style name, e.g. lj/cut
|
||||
pparam = parameter to perturb
|
||||
I,J = type pair(s) to set parameter for
|
||||
v_delta = variable with perturbation to apply (in the units of the parameter)
|
||||
<em>atom</em> args = aparam I v_delta
|
||||
aparam = parameter to perturb
|
||||
I = type to set parameter for
|
||||
v_delta = variable with perturbation to apply (in the units of the parameter)
|
||||
pstyle = pair style name, e.g. lj/cut
|
||||
pparam = parameter to perturb
|
||||
I,J = type pair(s) to set parameter for
|
||||
v_delta = variable with perturbation to apply (in the units of the parameter)
|
||||
<em>atom</em> args = aparam I v_delta
|
||||
aparam = parameter to perturb
|
||||
I = type to set parameter for
|
||||
v_delta = variable with perturbation to apply (in the units of the parameter)
|
||||
</pre>
|
||||
<ul class="simple">
|
||||
<li>zero or more keyword/value pairs may be appended</li>
|
||||
@ -155,11 +155,11 @@
|
||||
</ul>
|
||||
<pre class="literal-block">
|
||||
<em>tail</em> value = <em>no</em> or <em>yes</em>
|
||||
<em>no</em> = ignore tail correction to pair energies (usually small in fep)
|
||||
<em>yes</em> = include tail correction to pair energies
|
||||
<em>volume</em> value = <em>no</em> or <em>yes</em>
|
||||
<em>no</em> = ignore volume changes (e.g. in <em>NVE</em> or <em>NVT</em> trajectories)
|
||||
<em>yes</em> = include volume changes (e.g. in <em>NpT</em> trajectories)
|
||||
<em>no</em> = ignore tail correction to pair energies (usually small in fep)
|
||||
<em>yes</em> = include tail correction to pair energies
|
||||
<em>volume</em> value = <em>no</em> or <em>yes</em>
|
||||
<em>no</em> = ignore volume changes (e.g. in <em>NVE</em> or <em>NVT</em> trajectories)
|
||||
<em>yes</em> = include volume changes (e.g. in <em>NpT</em> trajectories)
|
||||
</pre>
|
||||
</div>
|
||||
<div class="section" id="examples">
|
||||
|
||||
@ -140,8 +140,8 @@
|
||||
</ul>
|
||||
<pre class="literal-block">
|
||||
<em>pair</em> value = <em>yes</em> or <em>no</em>
|
||||
<em>kspace</em> value = <em>yes</em> or <em>no</em>
|
||||
<em>boundary</em> value = <em>yes</em> or <em>no</em>
|
||||
<em>kspace</em> value = <em>yes</em> or <em>no</em>
|
||||
<em>boundary</em> value = <em>yes</em> or <em>no</em>
|
||||
</pre>
|
||||
</div>
|
||||
<div class="section" id="examples">
|
||||
|
||||
@ -138,11 +138,11 @@
|
||||
</ul>
|
||||
<pre class="literal-block">
|
||||
<em>extra</em> value = N
|
||||
N = # of extra degrees of freedom to subtract
|
||||
<em>dynamic</em> value = <em>yes</em> or <em>no</em>
|
||||
yes/no = do or do not recompute the number of atoms contributing to the temperature
|
||||
<em>thermo</em> value = <em>yes</em> or <em>no</em>
|
||||
yes/no = do or do not add contributions from fixes to the potential energy
|
||||
N = # of extra degrees of freedom to subtract
|
||||
<em>dynamic</em> value = <em>yes</em> or <em>no</em>
|
||||
yes/no = do or do not recompute the number of atoms contributing to the temperature
|
||||
<em>thermo</em> value = <em>yes</em> or <em>no</em>
|
||||
yes/no = do or do not add contributions from fixes to the potential energy
|
||||
</pre>
|
||||
</div>
|
||||
<div class="section" id="examples">
|
||||
|
||||
@ -139,10 +139,10 @@
|
||||
</ul>
|
||||
<pre class="literal-block">
|
||||
<em>dist</em> = pairwise distance
|
||||
<em>eng</em> = pairwise energy
|
||||
<em>force</em> = pairwise force
|
||||
<em>fx</em>,*fy*,*fz* = components of pairwise force
|
||||
<em>pN</em> = pair style specific quantities for allowed N values
|
||||
<em>eng</em> = pairwise energy
|
||||
<em>force</em> = pairwise force
|
||||
<em>fx</em>,*fy*,*fz* = components of pairwise force
|
||||
<em>pN</em> = pair style specific quantities for allowed N values
|
||||
</pre>
|
||||
</div>
|
||||
<div class="section" id="examples">
|
||||
|
||||
@ -137,73 +137,73 @@
|
||||
<li>input = one or more atom attributes</li>
|
||||
</ul>
|
||||
<div class="highlight-python"><div class="highlight"><pre>possible attributes = id, mol, proc, type, mass,
|
||||
x, y, z, xs, ys, zs, xu, yu, zu, ix, iy, iz,
|
||||
vx, vy, vz, fx, fy, fz,
|
||||
q, mux, muy, muz, mu,
|
||||
radius, diameter, omegax, omegay, omegaz,
|
||||
angmomx, angmomy, angmomz,
|
||||
shapex,shapey, shapez,
|
||||
quatw, quati, quatj, quatk, tqx, tqy, tqz,
|
||||
end1x, end1y, end1z, end2x, end2y, end2z,
|
||||
corner1x, corner1y, corner1z,
|
||||
corner2x, corner2y, corner2z,
|
||||
corner3x, corner3y, corner3z,
|
||||
nbonds,
|
||||
vfrac, s0,
|
||||
spin, eradius, ervel, erforce,
|
||||
rho, drho, e, de, cv,
|
||||
i_name, d_name
|
||||
x, y, z, xs, ys, zs, xu, yu, zu, ix, iy, iz,
|
||||
vx, vy, vz, fx, fy, fz,
|
||||
q, mux, muy, muz, mu,
|
||||
radius, diameter, omegax, omegay, omegaz,
|
||||
angmomx, angmomy, angmomz,
|
||||
shapex,shapey, shapez,
|
||||
quatw, quati, quatj, quatk, tqx, tqy, tqz,
|
||||
end1x, end1y, end1z, end2x, end2y, end2z,
|
||||
corner1x, corner1y, corner1z,
|
||||
corner2x, corner2y, corner2z,
|
||||
corner3x, corner3y, corner3z,
|
||||
nbonds,
|
||||
vfrac, s0,
|
||||
spin, eradius, ervel, erforce,
|
||||
rho, drho, e, de, cv,
|
||||
i_name, d_name
|
||||
</pre></div>
|
||||
</div>
|
||||
<div class="highlight-python"><div class="highlight"><pre>id = atom ID
|
||||
mol = molecule ID
|
||||
proc = ID of processor that owns atom
|
||||
type = atom type
|
||||
mass = atom mass
|
||||
x,y,z = unscaled atom coordinates
|
||||
xs,ys,zs = scaled atom coordinates
|
||||
xu,yu,zu = unwrapped atom coordinates
|
||||
ix,iy,iz = box image that the atom is in
|
||||
vx,vy,vz = atom velocities
|
||||
fx,fy,fz = forces on atoms
|
||||
q = atom charge
|
||||
mux,muy,muz = orientation of dipole moment of atom
|
||||
mu = magnitude of dipole moment of atom
|
||||
radius,diameter = radius,diameter of spherical particle
|
||||
omegax,omegay,omegaz = angular velocity of spherical particle
|
||||
angmomx,angmomy,angmomz = angular momentum of aspherical particle
|
||||
shapex,shapey,shapez = 3 diameters of aspherical particle
|
||||
quatw,quati,quatj,quatk = quaternion components for aspherical or body particles
|
||||
tqx,tqy,tqz = torque on finite-size particles
|
||||
end12x, end12y, end12z = end points of line segment
|
||||
corner123x, corner123y, corner123z = corner points of triangle
|
||||
nbonds = number of bonds assigned to an atom
|
||||
mol = molecule ID
|
||||
proc = ID of processor that owns atom
|
||||
type = atom type
|
||||
mass = atom mass
|
||||
x,y,z = unscaled atom coordinates
|
||||
xs,ys,zs = scaled atom coordinates
|
||||
xu,yu,zu = unwrapped atom coordinates
|
||||
ix,iy,iz = box image that the atom is in
|
||||
vx,vy,vz = atom velocities
|
||||
fx,fy,fz = forces on atoms
|
||||
q = atom charge
|
||||
mux,muy,muz = orientation of dipole moment of atom
|
||||
mu = magnitude of dipole moment of atom
|
||||
radius,diameter = radius,diameter of spherical particle
|
||||
omegax,omegay,omegaz = angular velocity of spherical particle
|
||||
angmomx,angmomy,angmomz = angular momentum of aspherical particle
|
||||
shapex,shapey,shapez = 3 diameters of aspherical particle
|
||||
quatw,quati,quatj,quatk = quaternion components for aspherical or body particles
|
||||
tqx,tqy,tqz = torque on finite-size particles
|
||||
end12x, end12y, end12z = end points of line segment
|
||||
corner123x, corner123y, corner123z = corner points of triangle
|
||||
nbonds = number of bonds assigned to an atom
|
||||
</pre></div>
|
||||
</div>
|
||||
<div class="highlight-python"><div class="highlight"><pre>PERI package per-atom properties:
|
||||
vfrac = ???
|
||||
s0 = ???
|
||||
vfrac = ???
|
||||
s0 = ???
|
||||
</pre></div>
|
||||
</div>
|
||||
<div class="highlight-python"><div class="highlight"><pre>USER-EFF and USER-AWPMD package per-atom properties:
|
||||
spin = electron spin
|
||||
eradius = electron radius
|
||||
ervel = electron radial velocity
|
||||
erforce = electron radial force
|
||||
spin = electron spin
|
||||
eradius = electron radius
|
||||
ervel = electron radial velocity
|
||||
erforce = electron radial force
|
||||
</pre></div>
|
||||
</div>
|
||||
<div class="highlight-python"><div class="highlight"><pre>USER-SPH package per-atom properties:
|
||||
rho = ???
|
||||
drho = ???
|
||||
e = ???
|
||||
de = ???
|
||||
cv = ???
|
||||
rho = ???
|
||||
drho = ???
|
||||
e = ???
|
||||
de = ???
|
||||
cv = ???
|
||||
</pre></div>
|
||||
</div>
|
||||
<pre class="literal-block">
|
||||
<a class="reference internal" href="fix_property_atom.html"><em>fix property/atom</em></a> per-atom properties:
|
||||
i_name = custom integer vector with name
|
||||
d_name = custom integer vector with name
|
||||
i_name = custom integer vector with name
|
||||
d_name = custom integer vector with name
|
||||
</pre>
|
||||
</div>
|
||||
<div class="section" id="examples">
|
||||
|
||||
@ -138,9 +138,9 @@
|
||||
</ul>
|
||||
<pre class="literal-block">
|
||||
attributes = count, id, coord1, coord2, coord3
|
||||
count = # of atoms in chunk
|
||||
id = original chunk IDs before compression by <a class="reference internal" href="compute_chunk_atom.html"><em>compute chunk/atom</em></a>
|
||||
coord123 = coordinates for spatial bins calculated by <a class="reference internal" href="compute_chunk_atom.html"><em>compute chunk/atom</em></a>
|
||||
count = # of atoms in chunk
|
||||
id = original chunk IDs before compression by <a class="reference internal" href="compute_chunk_atom.html"><em>compute chunk/atom</em></a>
|
||||
coord123 = coordinates for spatial bins calculated by <a class="reference internal" href="compute_chunk_atom.html"><em>compute chunk/atom</em></a>
|
||||
</pre>
|
||||
</div>
|
||||
<div class="section" id="examples">
|
||||
|
||||
@ -137,25 +137,25 @@
|
||||
<li>input = one or more attributes</li>
|
||||
</ul>
|
||||
<div class="highlight-python"><div class="highlight"><pre>possible attributes = natom1 natom2 ntype1 ntype2
|
||||
patom1 patom2 ptype1 ptype2
|
||||
batom1 batom2 btype
|
||||
aatom1 aatom2 aatom3 atype
|
||||
datom1 datom2 datom3 dtype
|
||||
iatom1 iatom2 iatom3 itype
|
||||
patom1 patom2 ptype1 ptype2
|
||||
batom1 batom2 btype
|
||||
aatom1 aatom2 aatom3 atype
|
||||
datom1 datom2 datom3 dtype
|
||||
iatom1 iatom2 iatom3 itype
|
||||
</pre></div>
|
||||
</div>
|
||||
<div class="highlight-python"><div class="highlight"><pre>natom1, natom2 = IDs of 2 atoms in each pair (within neighbor cutoff)
|
||||
ntype1, ntype2 = type of 2 atoms in each pair (within neighbor cutoff)
|
||||
patom1, patom2 = IDs of 2 atoms in each pair (within force cutoff)
|
||||
ptype1, ptype2 = type of 2 atoms in each pair (within force cutoff)
|
||||
batom1, batom2 = IDs of 2 atoms in each bond
|
||||
btype = bond type of each bond
|
||||
aatom1, aatom2, aatom3 = IDs of 3 atoms in each angle
|
||||
atype = angle type of each angle
|
||||
datom1, datom2, datom3, datom4 = IDs of 4 atoms in each dihedral
|
||||
dtype = dihedral type of each dihedral
|
||||
iatom1, iatom2, iatom3, iatom4 = IDs of 4 atoms in each improper
|
||||
itype = improper type of each improper
|
||||
ntype1, ntype2 = type of 2 atoms in each pair (within neighbor cutoff)
|
||||
patom1, patom2 = IDs of 2 atoms in each pair (within force cutoff)
|
||||
ptype1, ptype2 = type of 2 atoms in each pair (within force cutoff)
|
||||
batom1, batom2 = IDs of 2 atoms in each bond
|
||||
btype = bond type of each bond
|
||||
aatom1, aatom2, aatom3 = IDs of 3 atoms in each angle
|
||||
atype = angle type of each angle
|
||||
datom1, datom2, datom3, datom4 = IDs of 4 atoms in each dihedral
|
||||
dtype = dihedral type of each dihedral
|
||||
iatom1, iatom2, iatom3, iatom4 = IDs of 4 atoms in each improper
|
||||
itype = improper type of each improper
|
||||
</pre></div>
|
||||
</div>
|
||||
</div>
|
||||
|
||||
@ -140,8 +140,8 @@
|
||||
</ul>
|
||||
<pre class="literal-block">
|
||||
<em>reduce</em> arg = none
|
||||
<em>reduce/region</em> arg = region-ID
|
||||
region-ID = ID of region to use for choosing atoms
|
||||
<em>reduce/region</em> arg = region-ID
|
||||
region-ID = ID of region to use for choosing atoms
|
||||
</pre>
|
||||
<ul class="simple">
|
||||
<li>mode = <em>sum</em> or <em>min</em> or <em>max</em> or <em>ave</em> or <em>sumsq</em> or <em>avesq</em></li>
|
||||
@ -149,11 +149,11 @@
|
||||
<li>input = x, y, z, vx, vy, vz, fx, fy, fz, c_ID, c_ID[N], f_ID, f_ID[N], v_name</li>
|
||||
</ul>
|
||||
<div class="highlight-python"><div class="highlight"><pre>x,y,z,vx,vy,vz,fx,fy,fz = atom attribute (position, velocity, force component)
|
||||
c_ID = per-atom or local vector calculated by a compute with ID
|
||||
c_ID[I] = Ith column of per-atom or local array calculated by a compute with ID
|
||||
f_ID = per-atom or local vector calculated by a fix with ID
|
||||
f_ID[I] = Ith column of per-atom or local array calculated by a fix with ID
|
||||
v_name = per-atom vector calculated by an atom-style variable with name
|
||||
c_ID = per-atom or local vector calculated by a compute with ID
|
||||
c_ID[I] = Ith column of per-atom or local array calculated by a compute with ID
|
||||
f_ID = per-atom or local vector calculated by a fix with ID
|
||||
f_ID[I] = Ith column of per-atom or local array calculated by a fix with ID
|
||||
v_name = per-atom vector calculated by an atom-style variable with name
|
||||
</pre></div>
|
||||
</div>
|
||||
<ul class="simple">
|
||||
@ -162,8 +162,8 @@
|
||||
</ul>
|
||||
<pre class="literal-block">
|
||||
<em>replace</em> args = vec1 vec2
|
||||
vec1 = reduced value from this input vector will be replaced
|
||||
vec2 = replace it with vec1[N] where N is index of max/min value from vec2
|
||||
vec1 = reduced value from this input vector will be replaced
|
||||
vec2 = replace it with vec1[N] where N is index of max/min value from vec2
|
||||
</pre>
|
||||
</div>
|
||||
<div class="section" id="examples">
|
||||
|
||||
@ -141,18 +141,18 @@
|
||||
</ul>
|
||||
<pre class="literal-block">
|
||||
<em>Kmax</em> value = Maximum distance explored from reciprocal space origin
|
||||
(inverse length units)
|
||||
<em>Zone</em> values = z1 z2 z3
|
||||
z1,z2,z3 = Zone axis of incident radiation. If z1=z2=z3=0 all
|
||||
reciprocal space will be meshed up to <em>Kmax</em>
|
||||
<em>dR_Ewald</em> value = Thickness of Ewald sphere slice intercepting
|
||||
reciprocal space (inverse length units)
|
||||
<em>c</em> values = c1 c2 c3
|
||||
c1,c2,c3 = parameters to adjust the spacing of the reciprocal
|
||||
lattice nodes in the h, k, and l directions respectively
|
||||
<em>manual</em> = flag to use manual spacing of reciprocal lattice points
|
||||
based on the values of the <em>c</em> parameters
|
||||
<em>echo</em> = flag to provide extra output for debugging purposes
|
||||
(inverse length units)
|
||||
<em>Zone</em> values = z1 z2 z3
|
||||
z1,z2,z3 = Zone axis of incident radiation. If z1=z2=z3=0 all
|
||||
reciprocal space will be meshed up to <em>Kmax</em>
|
||||
<em>dR_Ewald</em> value = Thickness of Ewald sphere slice intercepting
|
||||
reciprocal space (inverse length units)
|
||||
<em>c</em> values = c1 c2 c3
|
||||
c1,c2,c3 = parameters to adjust the spacing of the reciprocal
|
||||
lattice nodes in the h, k, and l directions respectively
|
||||
<em>manual</em> = flag to use manual spacing of reciprocal lattice points
|
||||
based on the values of the <em>c</em> parameters
|
||||
<em>echo</em> = flag to provide extra output for debugging purposes
|
||||
</pre>
|
||||
</div>
|
||||
<div class="section" id="examples">
|
||||
@ -182,7 +182,7 @@ position of each atom, fj are atomic scattering factors.</p>
|
||||
reciprocal lattice nodes. The mesh spacing is defined either (a) by
|
||||
the entire simulation domain or (b) manually using selected values as
|
||||
shown in the 2D diagram below.</p>
|
||||
<a data-lightbox="group-a30c5cda-68d7-4c1b-a26c-5316847c9103"
|
||||
<a data-lightbox="group-default"
|
||||
href="_images/saed_mesh.jpg"
|
||||
class=""
|
||||
title=""
|
||||
@ -219,7 +219,7 @@ intersecting Ewald sphere. Diffraction intensities will only be
|
||||
computed at the intersection of the reciprocal lattice mesh and a
|
||||
<em>dR_Ewald</em> thick surface of the Ewald sphere. See the example 3D
|
||||
intestiety data and the intersection of a [010] zone axis in the below image.</p>
|
||||
<a data-lightbox="group-8fb8943f-067b-4621-ab85-ef00bd1882fd"
|
||||
<a data-lightbox="group-default"
|
||||
href="_images/saed_ewald_intersect.jpg"
|
||||
class=""
|
||||
title=""
|
||||
|
||||
@ -140,9 +140,9 @@
|
||||
<li>input = c_ID, c_ID[N], f_ID, f_ID[N]</li>
|
||||
</ul>
|
||||
<div class="highlight-python"><div class="highlight"><pre>c_ID = global vector calculated by a compute with ID
|
||||
c_ID[I] = Ith column of global array calculated by a compute with ID
|
||||
f_ID = global vector calculated by a fix with ID
|
||||
f_ID[I] = Ith column of global array calculated by a fix with ID
|
||||
c_ID[I] = Ith column of global array calculated by a compute with ID
|
||||
f_ID = global vector calculated by a fix with ID
|
||||
f_ID[I] = Ith column of global array calculated by a fix with ID
|
||||
</pre></div>
|
||||
</div>
|
||||
</div>
|
||||
|
||||
@ -152,14 +152,14 @@ compute ID group-ID snav/atom ntypes rcutfac rfac0 twojmax R_1 R_2 ... w_1 w_2 .
|
||||
</ul>
|
||||
<pre class="literal-block">
|
||||
<em>diagonal</em> value = <em>0</em> or <em>1</em> or <em>2</em> or <em>3</em>
|
||||
<em>0</em> = all j1, j2, j <= twojmax, j2 <= j1
|
||||
<em>1</em> = subset satisfying j1 == j2
|
||||
<em>2</em> = subset satisfying j1 == j2 == j3
|
||||
<em>3</em> = subset satisfying j2 <= j1 <= j
|
||||
<em>rmin0</em> value = parameter in distance to angle conversion (distance units)
|
||||
<em>switchflag</em> value = <em>0</em> or <em>1</em>
|
||||
<em>0</em> = do not use switching function
|
||||
<em>1</em> = use switching function
|
||||
<em>0</em> = all j1, j2, j <= twojmax, j2 <= j1
|
||||
<em>1</em> = subset satisfying j1 == j2
|
||||
<em>2</em> = subset satisfying j1 == j2 == j3
|
||||
<em>3</em> = subset satisfying j2 <= j1 <= j
|
||||
<em>rmin0</em> value = parameter in distance to angle conversion (distance units)
|
||||
<em>switchflag</em> value = <em>0</em> or <em>1</em>
|
||||
<em>0</em> = do not use switching function
|
||||
<em>1</em> = use switching function
|
||||
</pre>
|
||||
</div>
|
||||
<div class="section" id="examples">
|
||||
|
||||
@ -139,10 +139,10 @@
|
||||
</ul>
|
||||
<pre class="literal-block">
|
||||
<em>bias</em> value = bias-ID
|
||||
bias-ID = ID of a temperature compute that removes a velocity bias
|
||||
<em>dof</em> value = <em>all</em> or <em>rotate</em>
|
||||
all = compute temperature of translational and rotational degrees of freedom
|
||||
rotate = compute temperature of just rotational degrees of freedom
|
||||
bias-ID = ID of a temperature compute that removes a velocity bias
|
||||
<em>dof</em> value = <em>all</em> or <em>rotate</em>
|
||||
all = compute temperature of translational and rotational degrees of freedom
|
||||
rotate = compute temperature of just rotational degrees of freedom
|
||||
</pre>
|
||||
</div>
|
||||
<div class="section" id="examples">
|
||||
|
||||
@ -139,8 +139,8 @@
|
||||
<li>value = <em>temp</em> or <em>kecom</em> or <em>internal</em></li>
|
||||
</ul>
|
||||
<div class="highlight-python"><div class="highlight"><pre>temp = temperature of each chunk
|
||||
kecom = kinetic energy of each chunk based on velocity of center of mass
|
||||
internal = internal kinetic energy of each chunk
|
||||
kecom = kinetic energy of each chunk based on velocity of center of mass
|
||||
internal = internal kinetic energy of each chunk
|
||||
</pre></div>
|
||||
</div>
|
||||
<ul class="simple">
|
||||
@ -149,14 +149,14 @@
|
||||
</ul>
|
||||
<pre class="literal-block">
|
||||
<em>com</em> value = <em>yes</em> or <em>no</em>
|
||||
yes = subtract center-of-mass velocity from each chunk before calculating temperature
|
||||
no = do not subtract center-of-mass velocity
|
||||
<em>bias</em> value = bias-ID
|
||||
bias-ID = ID of a temperature compute that removes a velocity bias
|
||||
<em>adof</em> value = dof_per_atom
|
||||
dof_per_atom = define this many degrees-of-freedom per atom
|
||||
<em>cdof</em> value = dof_per_chunk
|
||||
dof_per_chunk = define this many degrees-of-freedom per chunk
|
||||
yes = subtract center-of-mass velocity from each chunk before calculating temperature
|
||||
no = do not subtract center-of-mass velocity
|
||||
<em>bias</em> value = bias-ID
|
||||
bias-ID = ID of a temperature compute that removes a velocity bias
|
||||
<em>adof</em> value = dof_per_atom
|
||||
dof_per_atom = define this many degrees-of-freedom per atom
|
||||
<em>cdof</em> value = dof_per_chunk
|
||||
dof_per_chunk = define this many degrees-of-freedom per chunk
|
||||
</pre>
|
||||
</div>
|
||||
<div class="section" id="examples">
|
||||
|
||||
@ -167,7 +167,10 @@ dim = 2 or 3 = dimensionality of the simulation, N = number of atoms
|
||||
in the group, k = Boltzmann constant, and T = temperature. Note that
|
||||
the velocity of each core or shell atom used in the KE calculation is
|
||||
the velocity of the center-of-mass (COM) of the core/shell pair the
|
||||
atom is part of.</p>
|
||||
atom is part of. Note that atoms that are not core or shell particles
|
||||
are also included in the temperature calculation (if they are in the
|
||||
specified group-ID); they contribute to the total kinetic energy in
|
||||
the usual way.</p>
|
||||
<p>A kinetic energy tensor, stored as a 6-element vector, is also
|
||||
calculated by this compute for use in the computation of a pressure
|
||||
tensor. The formula for the components of the tensor is the same as
|
||||
|
||||
@ -41,7 +41,7 @@ For this compute, core and shell particles are specified by two
|
||||
respective group IDs, which can be defined using the
|
||||
"group"_group.html command. The number of atoms in the two groups
|
||||
must be the same and there should be one bond defined between a pair
|
||||
of atoms in the two groups.
|
||||
of atoms in the two groups.
|
||||
|
||||
The temperature is calculated by the formula KE = dim/2 N k T, where
|
||||
KE = total kinetic energy of the group of atoms (sum of 1/2 m v^2),
|
||||
@ -49,7 +49,10 @@ dim = 2 or 3 = dimensionality of the simulation, N = number of atoms
|
||||
in the group, k = Boltzmann constant, and T = temperature. Note that
|
||||
the velocity of each core or shell atom used in the KE calculation is
|
||||
the velocity of the center-of-mass (COM) of the core/shell pair the
|
||||
atom is part of.
|
||||
atom is part of. Note that atoms that are not core or shell particles
|
||||
are also included in the temperature calculation (if they are in the
|
||||
specified group-ID); they contribute to the total kinetic energy in
|
||||
the usual way.
|
||||
|
||||
A kinetic energy tensor, stored as a 6-element vector, is also
|
||||
calculated by this compute for use in the computation of a pressure
|
||||
|
||||
@ -139,13 +139,13 @@
|
||||
</ul>
|
||||
<pre class="literal-block">
|
||||
<em>x</em> arg = Nx
|
||||
<em>y</em> arg = Ny
|
||||
<em>z</em> arg = Nz
|
||||
<em>xy</em> args = Nx Ny
|
||||
<em>yz</em> args = Ny Nz
|
||||
<em>xz</em> args = Nx Nz
|
||||
<em>xyz</em> args = Nx Ny Nz
|
||||
Nx,Ny,Nz = number of velocity bins in x,y,z dimensions
|
||||
<em>y</em> arg = Ny
|
||||
<em>z</em> arg = Nz
|
||||
<em>xy</em> args = Nx Ny
|
||||
<em>yz</em> args = Ny Nz
|
||||
<em>xz</em> args = Nx Nz
|
||||
<em>xyz</em> args = Nx Ny Nz
|
||||
Nx,Ny,Nz = number of velocity bins in x,y,z dimensions
|
||||
</pre>
|
||||
<ul class="simple">
|
||||
<li>zero or more keyword/value pairs may be appended</li>
|
||||
|
||||
@ -139,10 +139,10 @@
|
||||
</ul>
|
||||
<pre class="literal-block">
|
||||
<em>bias</em> value = bias-ID
|
||||
bias-ID = ID of a temperature compute that removes a velocity bias
|
||||
<em>dof</em> value = <em>all</em> or <em>rotate</em>
|
||||
all = compute temperature of translational and rotational degrees of freedom
|
||||
rotate = compute temperature of just rotational degrees of freedom
|
||||
bias-ID = ID of a temperature compute that removes a velocity bias
|
||||
<em>dof</em> value = <em>all</em> or <em>rotate</em>
|
||||
all = compute temperature of translational and rotational degrees of freedom
|
||||
rotate = compute temperature of just rotational degrees of freedom
|
||||
</pre>
|
||||
</div>
|
||||
<div class="section" id="examples">
|
||||
|
||||
@ -139,17 +139,17 @@
|
||||
</ul>
|
||||
<pre class="literal-block">
|
||||
pair style args = atype v_name1 v_name2
|
||||
atype = atom type (see asterisk form below)
|
||||
v_name1 = variable with name1 that is energy scale factor and function of lambda
|
||||
v_name2 = variable with name2 that is derivative of v_name1 with respect to lambda
|
||||
<em>tail</em> args = atype v_name1 v_name2
|
||||
atype = atom type (see asterisk form below)
|
||||
v_name1 = variable with name1 that is energy tail correction scale factor and function of lambda
|
||||
v_name2 = variable with name2 that is derivative of v_name1 with respect to lambda
|
||||
<em>kspace</em> args = atype v_name1 v_name2
|
||||
atype = atom type (see asterisk form below)
|
||||
v_name1 = variable with name1 that is K-Space scale factor and function of lambda
|
||||
v_name2 = variable with name2 that is derivative of v_name1 with respect to lambda
|
||||
atype = atom type (see asterisk form below)
|
||||
v_name1 = variable with name1 that is energy scale factor and function of lambda
|
||||
v_name2 = variable with name2 that is derivative of v_name1 with respect to lambda
|
||||
<em>tail</em> args = atype v_name1 v_name2
|
||||
atype = atom type (see asterisk form below)
|
||||
v_name1 = variable with name1 that is energy tail correction scale factor and function of lambda
|
||||
v_name2 = variable with name2 that is derivative of v_name1 with respect to lambda
|
||||
<em>kspace</em> args = atype v_name1 v_name2
|
||||
atype = atom type (see asterisk form below)
|
||||
v_name1 = variable with name1 that is K-Space scale factor and function of lambda
|
||||
v_name2 = variable with name2 that is derivative of v_name1 with respect to lambda
|
||||
</pre>
|
||||
</div>
|
||||
<div class="section" id="examples">
|
||||
|
||||
@ -139,18 +139,18 @@
|
||||
</ul>
|
||||
<pre class="literal-block">
|
||||
<em>only_group</em> = no arg
|
||||
<em>occupation</em> = no arg
|
||||
<em>surface</em> arg = sgroup-ID
|
||||
sgroup-ID = compute the dividing surface between group-ID and sgroup-ID
|
||||
this keyword adds a third column to the compute output
|
||||
<em>radius</em> arg = v_r
|
||||
v_r = radius atom style variable for a poly-disperse Voronoi tessellation
|
||||
<em>edge_histo</em> arg = maxedge
|
||||
maxedge = maximum number of Voronoi cell edges to be accounted in the histogram
|
||||
<em>edge_threshold</em> arg = minlength
|
||||
minlength = minimum length for an edge to be counted
|
||||
<em>face_threshold</em> arg = minarea
|
||||
minarea = minimum area for a face to be counted
|
||||
<em>occupation</em> = no arg
|
||||
<em>surface</em> arg = sgroup-ID
|
||||
sgroup-ID = compute the dividing surface between group-ID and sgroup-ID
|
||||
this keyword adds a third column to the compute output
|
||||
<em>radius</em> arg = v_r
|
||||
v_r = radius atom style variable for a poly-disperse Voronoi tessellation
|
||||
<em>edge_histo</em> arg = maxedge
|
||||
maxedge = maximum number of Voronoi cell edges to be accounted in the histogram
|
||||
<em>edge_threshold</em> arg = minlength
|
||||
minlength = minimum length for an edge to be counted
|
||||
<em>face_threshold</em> arg = minarea
|
||||
minarea = minimum area for a face to be counted
|
||||
</pre>
|
||||
</div>
|
||||
<div class="section" id="examples">
|
||||
|
||||
@ -141,16 +141,16 @@
|
||||
</ul>
|
||||
<pre class="literal-block">
|
||||
<em>2Theta</em> values = Min2Theta Max2Theta
|
||||
Min2Theta,Max2Theta = minimum and maximum 2 theta range to explore
|
||||
(radians or degrees)
|
||||
<em>c</em> values = c1 c2 c3
|
||||
c1,c2,c3 = parameters to adjust the spacing of the reciprocal
|
||||
lattice nodes in the h, k, and l directions respectively
|
||||
<em>LP</em> value = switch to apply Lorentz-polarization factor
|
||||
0/1 = off/on
|
||||
<em>manual</em> = flag to use manual spacing of reciprocal lattice points
|
||||
based on the values of the <em>c</em> parameters
|
||||
<em>echo</em> = flag to provide extra output for debugging purposes
|
||||
Min2Theta,Max2Theta = minimum and maximum 2 theta range to explore
|
||||
(radians or degrees)
|
||||
<em>c</em> values = c1 c2 c3
|
||||
c1,c2,c3 = parameters to adjust the spacing of the reciprocal
|
||||
lattice nodes in the h, k, and l directions respectively
|
||||
<em>LP</em> value = switch to apply Lorentz-polarization factor
|
||||
0/1 = off/on
|
||||
<em>manual</em> = flag to use manual spacing of reciprocal lattice points
|
||||
based on the values of the <em>c</em> parameters
|
||||
<em>echo</em> = flag to provide extra output for debugging purposes
|
||||
</pre>
|
||||
</div>
|
||||
<div class="section" id="examples">
|
||||
@ -185,7 +185,7 @@ the optional <em>LP</em> keyword.</p>
|
||||
reciprocal lattice nodes. The mesh spacing is defined either (a)
|
||||
by the entire simulation domain or (b) manually using selected values as
|
||||
shown in the 2D diagram below.</p>
|
||||
<a data-lightbox="group-1bc3f956-68f0-4615-8f93-ca954d919b33"
|
||||
<a data-lightbox="group-default"
|
||||
href="_images/xrd_mesh.jpg"
|
||||
class=""
|
||||
title=""
|
||||
|
||||
@ -137,14 +137,14 @@
|
||||
</ul>
|
||||
<pre class="literal-block">
|
||||
<em>box</em> args = none
|
||||
<em>region</em> args = region-ID
|
||||
region-ID = particles will only be created if contained in the region
|
||||
<em>single</em> args = x y z
|
||||
x,y,z = coordinates of a single particle (distance units)
|
||||
<em>random</em> args = N seed region-ID
|
||||
N = number of particles to create
|
||||
seed = random # seed (positive integer)
|
||||
region-ID = create atoms within this region, use NULL for entire simulation box
|
||||
<em>region</em> args = region-ID
|
||||
region-ID = particles will only be created if contained in the region
|
||||
<em>single</em> args = x y z
|
||||
x,y,z = coordinates of a single particle (distance units)
|
||||
<em>random</em> args = N seed region-ID
|
||||
N = number of particles to create
|
||||
seed = random # seed (positive integer)
|
||||
region-ID = create atoms within this region, use NULL for entire simulation box
|
||||
</pre>
|
||||
<ul class="simple">
|
||||
<li>zero or more keyword/value pairs may be appended</li>
|
||||
@ -152,22 +152,22 @@
|
||||
</ul>
|
||||
<pre class="literal-block">
|
||||
<em>mol</em> value = template-ID seed
|
||||
template-ID = ID of molecule template specified in a separate <a class="reference internal" href="molecule.html"><em>molecule</em></a> command
|
||||
seed = random # seed (positive integer)
|
||||
<em>basis</em> values = M itype
|
||||
M = which basis atom
|
||||
itype = atom type (1-N) to assign to this basis atom
|
||||
<em>remap</em> value = <em>yes</em> or <em>no</em>
|
||||
<em>var</em> value = name = variable name to evaluate for test of atom creation
|
||||
<em>set</em> values = dim vname
|
||||
dim = <em>x</em> or <em>y</em> or <em>z</em>
|
||||
name = name of variable to set with x,y,z atom position
|
||||
<em>rotate</em> values = Rx Ry Rz theta
|
||||
Rx,Ry,Rz = rotation vector for single molecule
|
||||
theta = rotation angle for single molecule (degrees)
|
||||
<em>units</em> value = <em>lattice</em> or <em>box</em>
|
||||
<em>lattice</em> = the geometry is defined in lattice units
|
||||
<em>box</em> = the geometry is defined in simulation box units
|
||||
template-ID = ID of molecule template specified in a separate <a class="reference internal" href="molecule.html"><em>molecule</em></a> command
|
||||
seed = random # seed (positive integer)
|
||||
<em>basis</em> values = M itype
|
||||
M = which basis atom
|
||||
itype = atom type (1-N) to assign to this basis atom
|
||||
<em>remap</em> value = <em>yes</em> or <em>no</em>
|
||||
<em>var</em> value = name = variable name to evaluate for test of atom creation
|
||||
<em>set</em> values = dim vname
|
||||
dim = <em>x</em> or <em>y</em> or <em>z</em>
|
||||
name = name of variable to set with x,y,z atom position
|
||||
<em>rotate</em> values = Rx Ry Rz theta
|
||||
Rx,Ry,Rz = rotation vector for single molecule
|
||||
theta = rotation angle for single molecule (degrees)
|
||||
<em>units</em> value = <em>lattice</em> or <em>box</em>
|
||||
<em>lattice</em> = the geometry is defined in lattice units
|
||||
<em>box</em> = the geometry is defined in simulation box units
|
||||
</pre>
|
||||
</div>
|
||||
<div class="section" id="examples">
|
||||
@ -341,7 +341,7 @@ variable v equal "(0.2*v_y*ylat * cos(v_xx/xlat * 2.0*PI*4.0/v_x) +
|
||||
create_atoms 1 box var v set x xx set y yy
|
||||
</pre></div>
|
||||
</div>
|
||||
<a data-lightbox="group-f4ec1fe6-16bb-4d3a-9ef0-ecf4b57f6ec2"
|
||||
<a data-lightbox="group-default"
|
||||
href="_images/sinusoid.jpg"
|
||||
class=""
|
||||
title=""
|
||||
|
||||
@ -139,14 +139,14 @@
|
||||
</ul>
|
||||
<pre class="literal-block">
|
||||
<em>bond/types</em> value = # of bond types
|
||||
<em>angle/types</em> value = # of angle types
|
||||
<em>dihedral/types</em> value = # of dihedral types
|
||||
<em>improper/types</em> value = # of improper types
|
||||
<em>extra/bond/per/atom</em> value = # of bonds per atom
|
||||
<em>extra/angle/per/atom</em> value = # of angles per atom
|
||||
<em>extra/dihedral/per/atom</em> value = # of dihedrals per atom
|
||||
<em>extra/improper/per/atom</em> value = # of impropers per atom
|
||||
<em>extra/special/per/atom</em> value = # of special neighbors per atom
|
||||
<em>angle/types</em> value = # of angle types
|
||||
<em>dihedral/types</em> value = # of dihedral types
|
||||
<em>improper/types</em> value = # of improper types
|
||||
<em>extra/bond/per/atom</em> value = # of bonds per atom
|
||||
<em>extra/angle/per/atom</em> value = # of angles per atom
|
||||
<em>extra/dihedral/per/atom</em> value = # of dihedrals per atom
|
||||
<em>extra/improper/per/atom</em> value = # of impropers per atom
|
||||
<em>extra/special/per/atom</em> value = # of special neighbors per atom
|
||||
</pre>
|
||||
</div>
|
||||
<div class="section" id="examples">
|
||||
|
||||
@ -136,15 +136,15 @@
|
||||
</ul>
|
||||
<pre class="literal-block">
|
||||
<em>group</em> args = group-ID
|
||||
<em>region</em> args = region-ID
|
||||
<em>overlap</em> args = cutoff group1-ID group2-ID
|
||||
cutoff = delete one atom from pairs of atoms within the cutoff (distance units)
|
||||
group1-ID = one atom in pair must be in this group
|
||||
group2-ID = other atom in pair must be in this group
|
||||
<em>porosity</em> args = region-ID fraction seed
|
||||
region-ID = region within which to perform deletions
|
||||
fraction = delete this fraction of atoms
|
||||
seed = random number seed (positive integer)
|
||||
<em>region</em> args = region-ID
|
||||
<em>overlap</em> args = cutoff group1-ID group2-ID
|
||||
cutoff = delete one atom from pairs of atoms within the cutoff (distance units)
|
||||
group1-ID = one atom in pair must be in this group
|
||||
group2-ID = other atom in pair must be in this group
|
||||
<em>porosity</em> args = region-ID fraction seed
|
||||
region-ID = region within which to perform deletions
|
||||
fraction = delete this fraction of atoms
|
||||
seed = random number seed (positive integer)
|
||||
</pre>
|
||||
<ul class="simple">
|
||||
<li>zero or more keyword/value pairs may be appended</li>
|
||||
@ -152,8 +152,8 @@
|
||||
</ul>
|
||||
<pre class="literal-block">
|
||||
<em>compress</em> value = <em>no</em> or <em>yes</em>
|
||||
<em>bond</em> value = <em>no</em> or <em>yes</em>
|
||||
<em>mol</em> value = <em>no</em> or <em>yes</em>
|
||||
<em>bond</em> value = <em>no</em> or <em>yes</em>
|
||||
<em>mol</em> value = <em>no</em> or <em>yes</em>
|
||||
</pre>
|
||||
</div>
|
||||
<div class="section" id="examples">
|
||||
|
||||
@ -143,12 +143,12 @@
|
||||
</ul>
|
||||
<pre class="literal-block">
|
||||
<em>multi</em> arg = none
|
||||
<em>atom</em> arg = an atom type or range of types (see below)
|
||||
<em>bond</em> arg = a bond type or range of types (see below)
|
||||
<em>angle</em> arg = an angle type or range of types (see below)
|
||||
<em>dihedral</em> arg = a dihedral type or range of types (see below)
|
||||
<em>improper</em> arg = an improper type or range of types (see below)
|
||||
<em>stats</em> arg = none
|
||||
<em>atom</em> arg = an atom type or range of types (see below)
|
||||
<em>bond</em> arg = a bond type or range of types (see below)
|
||||
<em>angle</em> arg = an angle type or range of types (see below)
|
||||
<em>dihedral</em> arg = a dihedral type or range of types (see below)
|
||||
<em>improper</em> arg = an improper type or range of types (see below)
|
||||
<em>stats</em> arg = none
|
||||
</pre>
|
||||
<ul class="simple">
|
||||
<li>zero or more keywords may be appended</li>
|
||||
|
||||
@ -137,26 +137,26 @@
|
||||
</ul>
|
||||
<pre class="literal-block">
|
||||
<em>move</em> args = delx dely delz
|
||||
delx,dely,delz = distance to displace in each dimension (distance units)
|
||||
<em>ramp</em> args = ddim dlo dhi dim clo chi
|
||||
ddim = <em>x</em> or <em>y</em> or <em>z</em>
|
||||
dlo,dhi = displacement distance between dlo and dhi (distance units)
|
||||
dim = <em>x</em> or <em>y</em> or <em>z</em>
|
||||
clo,chi = lower and upper bound of domain to displace (distance units)
|
||||
<em>random</em> args = dx dy dz seed
|
||||
dx,dy,dz = random displacement magnitude in each dimension (distance units)
|
||||
seed = random # seed (positive integer)
|
||||
<em>rotate</em> args = Px Py Pz Rx Ry Rz theta
|
||||
Px,Py,Pz = origin point of axis of rotation (distance units)
|
||||
Rx,Ry,Rz = axis of rotation vector
|
||||
theta = angle of rotation (degrees)
|
||||
delx,dely,delz = distance to displace in each dimension (distance units)
|
||||
<em>ramp</em> args = ddim dlo dhi dim clo chi
|
||||
ddim = <em>x</em> or <em>y</em> or <em>z</em>
|
||||
dlo,dhi = displacement distance between dlo and dhi (distance units)
|
||||
dim = <em>x</em> or <em>y</em> or <em>z</em>
|
||||
clo,chi = lower and upper bound of domain to displace (distance units)
|
||||
<em>random</em> args = dx dy dz seed
|
||||
dx,dy,dz = random displacement magnitude in each dimension (distance units)
|
||||
seed = random # seed (positive integer)
|
||||
<em>rotate</em> args = Px Py Pz Rx Ry Rz theta
|
||||
Px,Py,Pz = origin point of axis of rotation (distance units)
|
||||
Rx,Ry,Rz = axis of rotation vector
|
||||
theta = angle of rotation (degrees)
|
||||
</pre>
|
||||
<ul class="simple">
|
||||
<li>zero or more keyword/value pairs may be appended</li>
|
||||
</ul>
|
||||
<pre class="literal-block">
|
||||
keyword = <em>units</em>
|
||||
value = <em>box</em> or <em>lattice</em>
|
||||
value = <em>box</em> or <em>lattice</em>
|
||||
</pre>
|
||||
</div>
|
||||
<div class="section" id="examples">
|
||||
|
||||
|
Before Width: | Height: | Size: 599 KiB After Width: | Height: | Size: 9.0 KiB |
|
Before Width: | Height: | Size: 610 KiB After Width: | Height: | Size: 7.0 KiB |
|
Before Width: | Height: | Size: 572 KiB After Width: | Height: | Size: 2.4 KiB |
|
Before Width: | Height: | Size: 638 KiB After Width: | Height: | Size: 26 KiB |
|
Before Width: | Height: | Size: 675 KiB After Width: | Height: | Size: 29 KiB |
|
Before Width: | Height: | Size: 604 KiB After Width: | Height: | Size: 8.9 KiB |
|
Before Width: | Height: | Size: 651 KiB After Width: | Height: | Size: 57 KiB |
|
Before Width: | Height: | Size: 594 KiB After Width: | Height: | Size: 10 KiB |
|
Before Width: | Height: | Size: 661 KiB After Width: | Height: | Size: 13 KiB |
@ -44,7 +44,7 @@ temp/rescale</A>, <A HREF = "fix_nh.html">fix npt</A>, etc.
|
||||
respective group IDs, which can be defined using the
|
||||
<A HREF = "group.html">group</A> command. The number of atoms in the two groups
|
||||
must be the same and there should be one bond defined between a pair
|
||||
of atoms in the two groups.
|
||||
of atoms in the two groups.
|
||||
</P>
|
||||
<P>The temperature is calculated by the formula KE = dim/2 N k T, where
|
||||
KE = total kinetic energy of the group of atoms (sum of 1/2 m v^2),
|
||||
@ -52,7 +52,10 @@ dim = 2 or 3 = dimensionality of the simulation, N = number of atoms
|
||||
in the group, k = Boltzmann constant, and T = temperature. Note that
|
||||
the velocity of each core or shell atom used in the KE calculation is
|
||||
the velocity of the center-of-mass (COM) of the core/shell pair the
|
||||
atom is part of.
|
||||
atom is part of. Note that atoms that are not core or shell particles
|
||||
are also included in the temperature calculation (if they are in the
|
||||
specified group-ID); they contribute to the total kinetic energy in
|
||||
the usual way.
|
||||
</P>
|
||||
<P>A kinetic energy tensor, stored as a 6-element vector, is also
|
||||
calculated by this compute for use in the computation of a pressure
|
||||
|
||||
@ -333,9 +333,9 @@ referenced by the user for each subsequent fix gcmc command.
|
||||
</P>
|
||||
<P><B>Related commands:</B>
|
||||
</P>
|
||||
<P><A HREF = "fix_atom_swap.html">fix_atom_swap</A>,
|
||||
<A HREF = "fix_nvt.html">fix_nvt</A>, <A HREF = "neighbor.html">neighbor</A>,
|
||||
<A HREF = "fix_deposit.html">fix_deposit</A>, <A HREF = "fix_evaporate.html">fix_evaporate</A>,
|
||||
<P><A HREF = "fix_atom_swap.html">fix atom/swap</A>,
|
||||
<A HREF = "fix_nvt.html">fix nvt</A>, <A HREF = "neighbor.html">neighbor</A>,
|
||||
<A HREF = "fix_deposit.html">fix deposit</A>, <A HREF = "fix_evaporate.html">fix evaporate</A>,
|
||||
<A HREF = "delete_atoms.html">delete_atoms</A>
|
||||
</P>
|
||||
<P><B>Default:</B>
|
||||
|
||||
@ -34,7 +34,7 @@ pair_coeff * * GaN_sw.polymorphic GaN
|
||||
</CENTER>
|
||||
<CENTER><IMG SRC = "Eqs/polymorphic2.jpg">
|
||||
</CENTER>
|
||||
<CENTER><IMG SRC = "Eqs/polymorphic2.jpg">
|
||||
<CENTER><IMG SRC = "Eqs/polymorphic3.jpg">
|
||||
</CENTER>
|
||||
<P>where I, J, K represent species of atoms i, j, and k, i_1, ..., i_N
|
||||
represents a list of i's neighbors, delta_ij is a Direc constant
|
||||
|
||||
@ -150,12 +150,12 @@
|
||||
</ul>
|
||||
<pre class="literal-block">
|
||||
<em>atom</em> args = none
|
||||
<em>atom/mpiio</em> args = none
|
||||
<em>cfg</em> args = same as <em>custom</em> args, see below
|
||||
<em>cfg/mpiio</em> args = same as <em>custom</em> args, see below
|
||||
<em>dcd</em> args = none
|
||||
<em>xtc</em> args = none
|
||||
<em>xyz</em> args = none
|
||||
<em>atom/mpiio</em> args = none
|
||||
<em>cfg</em> args = same as <em>custom</em> args, see below
|
||||
<em>cfg/mpiio</em> args = same as <em>custom</em> args, see below
|
||||
<em>dcd</em> args = none
|
||||
<em>xtc</em> args = none
|
||||
<em>xyz</em> args = none
|
||||
</pre>
|
||||
<pre class="literal-block">
|
||||
<em>xyz/mpiio</em> args = none
|
||||
@ -171,53 +171,53 @@
|
||||
</pre>
|
||||
<pre class="literal-block">
|
||||
<em>local</em> args = list of local attributes
|
||||
possible attributes = index, c_ID, c_ID[N], f_ID, f_ID[N]
|
||||
index = enumeration of local values
|
||||
c_ID = local vector calculated by a compute with ID
|
||||
c_ID[N] = Nth column of local array calculated by a compute with ID
|
||||
f_ID = local vector calculated by a fix with ID
|
||||
f_ID[N] = Nth column of local array calculated by a fix with ID
|
||||
possible attributes = index, c_ID, c_ID[N], f_ID, f_ID[N]
|
||||
index = enumeration of local values
|
||||
c_ID = local vector calculated by a compute with ID
|
||||
c_ID[N] = Nth column of local array calculated by a compute with ID
|
||||
f_ID = local vector calculated by a fix with ID
|
||||
f_ID[N] = Nth column of local array calculated by a fix with ID
|
||||
</pre>
|
||||
<pre class="literal-block">
|
||||
<em>custom</em> or <em>custom/mpiio</em> args = list of atom attributes
|
||||
possible attributes = id, mol, proc, procp1, type, element, mass,
|
||||
x, y, z, xs, ys, zs, xu, yu, zu,
|
||||
xsu, ysu, zsu, ix, iy, iz,
|
||||
vx, vy, vz, fx, fy, fz,
|
||||
q, mux, muy, muz, mu,
|
||||
radius, diameter, omegax, omegay, omegaz,
|
||||
angmomx, angmomy, angmomz, tqx, tqy, tqz,
|
||||
c_ID, c_ID[N], f_ID, f_ID[N], v_name
|
||||
possible attributes = id, mol, proc, procp1, type, element, mass,
|
||||
x, y, z, xs, ys, zs, xu, yu, zu,
|
||||
xsu, ysu, zsu, ix, iy, iz,
|
||||
vx, vy, vz, fx, fy, fz,
|
||||
q, mux, muy, muz, mu,
|
||||
radius, diameter, omegax, omegay, omegaz,
|
||||
angmomx, angmomy, angmomz, tqx, tqy, tqz,
|
||||
c_ID, c_ID[N], f_ID, f_ID[N], v_name
|
||||
</pre>
|
||||
<pre class="literal-block">
|
||||
id = atom ID
|
||||
mol = molecule ID
|
||||
proc = ID of processor that owns atom
|
||||
procp1 = ID+1 of processor that owns atom
|
||||
type = atom type
|
||||
element = name of atom element, as defined by <a class="reference internal" href="dump_modify.html"><em>dump_modify</em></a> command
|
||||
mass = atom mass
|
||||
x,y,z = unscaled atom coordinates
|
||||
xs,ys,zs = scaled atom coordinates
|
||||
xu,yu,zu = unwrapped atom coordinates
|
||||
xsu,ysu,zsu = scaled unwrapped atom coordinates
|
||||
ix,iy,iz = box image that the atom is in
|
||||
vx,vy,vz = atom velocities
|
||||
fx,fy,fz = forces on atoms
|
||||
q = atom charge
|
||||
mux,muy,muz = orientation of dipole moment of atom
|
||||
mu = magnitude of dipole moment of atom
|
||||
radius,diameter = radius,diameter of spherical particle
|
||||
omegax,omegay,omegaz = angular velocity of spherical particle
|
||||
angmomx,angmomy,angmomz = angular momentum of aspherical particle
|
||||
tqx,tqy,tqz = torque on finite-size particles
|
||||
c_ID = per-atom vector calculated by a compute with ID
|
||||
c_ID[N] = Nth column of per-atom array calculated by a compute with ID
|
||||
f_ID = per-atom vector calculated by a fix with ID
|
||||
f_ID[N] = Nth column of per-atom array calculated by a fix with ID
|
||||
v_name = per-atom vector calculated by an atom-style variable with name
|
||||
d_name = per-atom floating point vector with name, managed by fix property/atom
|
||||
i_name = per-atom integer vector with name, managed by fix property/atom
|
||||
mol = molecule ID
|
||||
proc = ID of processor that owns atom
|
||||
procp1 = ID+1 of processor that owns atom
|
||||
type = atom type
|
||||
element = name of atom element, as defined by <a class="reference internal" href="dump_modify.html"><em>dump_modify</em></a> command
|
||||
mass = atom mass
|
||||
x,y,z = unscaled atom coordinates
|
||||
xs,ys,zs = scaled atom coordinates
|
||||
xu,yu,zu = unwrapped atom coordinates
|
||||
xsu,ysu,zsu = scaled unwrapped atom coordinates
|
||||
ix,iy,iz = box image that the atom is in
|
||||
vx,vy,vz = atom velocities
|
||||
fx,fy,fz = forces on atoms
|
||||
q = atom charge
|
||||
mux,muy,muz = orientation of dipole moment of atom
|
||||
mu = magnitude of dipole moment of atom
|
||||
radius,diameter = radius,diameter of spherical particle
|
||||
omegax,omegay,omegaz = angular velocity of spherical particle
|
||||
angmomx,angmomy,angmomz = angular momentum of aspherical particle
|
||||
tqx,tqy,tqz = torque on finite-size particles
|
||||
c_ID = per-atom vector calculated by a compute with ID
|
||||
c_ID[N] = Nth column of per-atom array calculated by a compute with ID
|
||||
f_ID = per-atom vector calculated by a fix with ID
|
||||
f_ID[N] = Nth column of per-atom array calculated by a fix with ID
|
||||
v_name = per-atom vector calculated by an atom-style variable with name
|
||||
d_name = per-atom floating point vector with name, managed by fix property/atom
|
||||
i_name = per-atom integer vector with name, managed by fix property/atom
|
||||
</pre>
|
||||
</div>
|
||||
<div class="section" id="examples">
|
||||
|
||||
@ -147,47 +147,47 @@
|
||||
</ul>
|
||||
<pre class="literal-block">
|
||||
<em>adiam</em> value = number = numeric value for atom diameter (distance units)
|
||||
<em>atom</em> = yes/no = do or do not draw atoms
|
||||
<em>bond</em> values = color width = color and width of bonds
|
||||
color = <em>atom</em> or <em>type</em> or <em>none</em>
|
||||
width = number or <em>atom</em> or <em>type</em> or <em>none</em>
|
||||
number = numeric value for bond width (distance units)
|
||||
<em>size</em> values = width height = size of images
|
||||
width = width of image in # of pixels
|
||||
height = height of image in # of pixels
|
||||
<em>view</em> values = theta phi = view of simulation box
|
||||
theta = view angle from +z axis (degrees)
|
||||
phi = azimuthal view angle (degrees)
|
||||
theta or phi can be a variable (see below)
|
||||
<em>center</em> values = flag Cx Cy Cz = center point of image
|
||||
flag = "s" for static, "d" for dynamic
|
||||
Cx,Cy,Cz = center point of image as fraction of box dimension (0.5 = center of box)
|
||||
Cx,Cy,Cz can be variables (see below)
|
||||
<em>up</em> values = Ux Uy Uz = direction that is "up" in image
|
||||
Ux,Uy,Uz = components of up vector
|
||||
Ux,Uy,Uz can be variables (see below)
|
||||
<em>zoom</em> value = zfactor = size that simulation box appears in image
|
||||
zfactor = scale image size by factor > 1 to enlarge, factor < 1 to shrink
|
||||
zfactor can be a variable (see below)
|
||||
<em>persp</em> value = pfactor = amount of "perspective" in image
|
||||
pfactor = amount of perspective (0 = none, < 1 = some, > 1 = highly skewed)
|
||||
pfactor can be a variable (see below)
|
||||
<em>box</em> values = yes/no diam = draw outline of simulation box
|
||||
yes/no = do or do not draw simulation box lines
|
||||
diam = diameter of box lines as fraction of shortest box length
|
||||
<em>axes</em> values = yes/no length diam = draw xyz axes
|
||||
yes/no = do or do not draw xyz axes lines next to simulation box
|
||||
length = length of axes lines as fraction of respective box lengths
|
||||
diam = diameter of axes lines as fraction of shortest box length
|
||||
<em>subbox</em> values = yes/no diam = draw outline of processor sub-domains
|
||||
yes/no = do or do not draw sub-domain lines
|
||||
diam = diameter of sub-domain lines as fraction of shortest box length
|
||||
<em>shiny</em> value = sfactor = shinyness of spheres and cylinders
|
||||
sfactor = shinyness of spheres and cylinders from 0.0 to 1.0
|
||||
<em>ssao</em> value = yes/no seed dfactor = SSAO depth shading
|
||||
yes/no = turn depth shading on/off
|
||||
seed = random # seed (positive integer)
|
||||
dfactor = strength of shading from 0.0 to 1.0
|
||||
<em>atom</em> = yes/no = do or do not draw atoms
|
||||
<em>bond</em> values = color width = color and width of bonds
|
||||
color = <em>atom</em> or <em>type</em> or <em>none</em>
|
||||
width = number or <em>atom</em> or <em>type</em> or <em>none</em>
|
||||
number = numeric value for bond width (distance units)
|
||||
<em>size</em> values = width height = size of images
|
||||
width = width of image in # of pixels
|
||||
height = height of image in # of pixels
|
||||
<em>view</em> values = theta phi = view of simulation box
|
||||
theta = view angle from +z axis (degrees)
|
||||
phi = azimuthal view angle (degrees)
|
||||
theta or phi can be a variable (see below)
|
||||
<em>center</em> values = flag Cx Cy Cz = center point of image
|
||||
flag = "s" for static, "d" for dynamic
|
||||
Cx,Cy,Cz = center point of image as fraction of box dimension (0.5 = center of box)
|
||||
Cx,Cy,Cz can be variables (see below)
|
||||
<em>up</em> values = Ux Uy Uz = direction that is "up" in image
|
||||
Ux,Uy,Uz = components of up vector
|
||||
Ux,Uy,Uz can be variables (see below)
|
||||
<em>zoom</em> value = zfactor = size that simulation box appears in image
|
||||
zfactor = scale image size by factor > 1 to enlarge, factor < 1 to shrink
|
||||
zfactor can be a variable (see below)
|
||||
<em>persp</em> value = pfactor = amount of "perspective" in image
|
||||
pfactor = amount of perspective (0 = none, < 1 = some, > 1 = highly skewed)
|
||||
pfactor can be a variable (see below)
|
||||
<em>box</em> values = yes/no diam = draw outline of simulation box
|
||||
yes/no = do or do not draw simulation box lines
|
||||
diam = diameter of box lines as fraction of shortest box length
|
||||
<em>axes</em> values = yes/no length diam = draw xyz axes
|
||||
yes/no = do or do not draw xyz axes lines next to simulation box
|
||||
length = length of axes lines as fraction of respective box lengths
|
||||
diam = diameter of axes lines as fraction of shortest box length
|
||||
<em>subbox</em> values = yes/no diam = draw outline of processor sub-domains
|
||||
yes/no = do or do not draw sub-domain lines
|
||||
diam = diameter of sub-domain lines as fraction of shortest box length
|
||||
<em>shiny</em> value = sfactor = shinyness of spheres and cylinders
|
||||
sfactor = shinyness of spheres and cylinders from 0.0 to 1.0
|
||||
<em>ssao</em> value = yes/no seed dfactor = SSAO depth shading
|
||||
yes/no = turn depth shading on/off
|
||||
seed = random # seed (positive integer)
|
||||
dfactor = strength of shading from 0.0 to 1.0
|
||||
</pre>
|
||||
</div>
|
||||
<div class="section" id="examples">
|
||||
@ -220,7 +220,7 @@ from an existing dump file, and using these dump commands in the rerun
|
||||
script to generate the images/movie.</p>
|
||||
<p>Here are two sample images, rendered as 1024x1024 JPEG files. Click
|
||||
to see the full-size images:</p>
|
||||
<DIV ALIGN=center><a data-lightbox="group-24b43492-ff03-4c03-99b9-12f1dbf252d5"
|
||||
<DIV ALIGN=center><a data-lightbox="group-default"
|
||||
href="_images/dump1.jpg"
|
||||
class=""
|
||||
title=""
|
||||
@ -230,7 +230,7 @@ to see the full-size images:</p>
|
||||
width="25%"
|
||||
height="auto"
|
||||
alt=""/>
|
||||
</a><a data-lightbox="group-2540cc67-71aa-4053-b959-61451a4e4f0d"
|
||||
</a><a data-lightbox="group-default"
|
||||
href="_images/dump2.jpg"
|
||||
class=""
|
||||
title=""
|
||||
|
||||
@ -139,39 +139,39 @@
|
||||
</ul>
|
||||
<pre class="literal-block">
|
||||
<em>append</em> arg = <em>yes</em> or <em>no</em>
|
||||
<em>buffer</em> arg = <em>yes</em> or <em>no</em>
|
||||
<em>element</em> args = E1 E2 ... EN, where N = # of atom types
|
||||
E1,...,EN = element name, e.g. C or Fe or Ga
|
||||
<em>every</em> arg = N
|
||||
N = dump every this many timesteps
|
||||
N can be a variable (see below)
|
||||
<em>fileper</em> arg = Np
|
||||
Np = write one file for every this many processors
|
||||
<em>first</em> arg = <em>yes</em> or <em>no</em>
|
||||
<em>format</em> arg = C-style format string for one line of output
|
||||
<em>flush</em> arg = <em>yes</em> or <em>no</em>
|
||||
<em>image</em> arg = <em>yes</em> or <em>no</em>
|
||||
<em>label</em> arg = string
|
||||
string = character string (e.g. BONDS) to use in header of dump local file
|
||||
<em>nfile</em> arg = Nf
|
||||
Nf = write this many files, one from each of Nf processors
|
||||
<em>pad</em> arg = Nchar = # of characters to convert timestep to
|
||||
<em>precision</em> arg = power-of-10 value from 10 to 1000000
|
||||
<em>region</em> arg = region-ID or "none"
|
||||
<em>scale</em> arg = <em>yes</em> or <em>no</em>
|
||||
<em>sfactor</em> arg = coordinate scaling factor (> 0.0)
|
||||
<em>tfactor</em> arg = time scaling factor (> 0.0)
|
||||
<em>sort</em> arg = <em>off</em> or <em>id</em> or N or -N
|
||||
off = no sorting of per-atom lines within a snapshot
|
||||
id = sort per-atom lines by atom ID
|
||||
N = sort per-atom lines in ascending order by the Nth column
|
||||
-N = sort per-atom lines in descending order by the Nth column
|
||||
<em>thresh</em> args = attribute operation value
|
||||
attribute = same attributes (x,fy,etotal,sxx,etc) used by dump custom style
|
||||
operation = "<" or "<=" or ">" or ">=" or "==" or "!="
|
||||
value = numeric value to compare to
|
||||
these 3 args can be replaced by the word "none" to turn off thresholding
|
||||
<em>unwrap</em> arg = <em>yes</em> or <em>no</em>
|
||||
<em>buffer</em> arg = <em>yes</em> or <em>no</em>
|
||||
<em>element</em> args = E1 E2 ... EN, where N = # of atom types
|
||||
E1,...,EN = element name, e.g. C or Fe or Ga
|
||||
<em>every</em> arg = N
|
||||
N = dump every this many timesteps
|
||||
N can be a variable (see below)
|
||||
<em>fileper</em> arg = Np
|
||||
Np = write one file for every this many processors
|
||||
<em>first</em> arg = <em>yes</em> or <em>no</em>
|
||||
<em>format</em> arg = C-style format string for one line of output
|
||||
<em>flush</em> arg = <em>yes</em> or <em>no</em>
|
||||
<em>image</em> arg = <em>yes</em> or <em>no</em>
|
||||
<em>label</em> arg = string
|
||||
string = character string (e.g. BONDS) to use in header of dump local file
|
||||
<em>nfile</em> arg = Nf
|
||||
Nf = write this many files, one from each of Nf processors
|
||||
<em>pad</em> arg = Nchar = # of characters to convert timestep to
|
||||
<em>precision</em> arg = power-of-10 value from 10 to 1000000
|
||||
<em>region</em> arg = region-ID or "none"
|
||||
<em>scale</em> arg = <em>yes</em> or <em>no</em>
|
||||
<em>sfactor</em> arg = coordinate scaling factor (> 0.0)
|
||||
<em>tfactor</em> arg = time scaling factor (> 0.0)
|
||||
<em>sort</em> arg = <em>off</em> or <em>id</em> or N or -N
|
||||
off = no sorting of per-atom lines within a snapshot
|
||||
id = sort per-atom lines by atom ID
|
||||
N = sort per-atom lines in ascending order by the Nth column
|
||||
-N = sort per-atom lines in descending order by the Nth column
|
||||
<em>thresh</em> args = attribute operation value
|
||||
attribute = same attributes (x,fy,etotal,sxx,etc) used by dump custom style
|
||||
operation = "<" or "<=" or ">" or ">=" or "==" or "!="
|
||||
value = numeric value to compare to
|
||||
these 3 args can be replaced by the word "none" to turn off thresholding
|
||||
<em>unwrap</em> arg = <em>yes</em> or <em>no</em>
|
||||
</pre>
|
||||
<ul class="simple">
|
||||
<li>these keywords apply only to the <em>image</em> and <em>movie</em> <a class="reference internal" href="dump_image.html"><em>styles</em></a></li>
|
||||
@ -179,48 +179,48 @@
|
||||
</ul>
|
||||
<pre class="literal-block">
|
||||
<em>acolor</em> args = type color
|
||||
type = atom type or range of types (see below)
|
||||
color = name of color or color1/color2/...
|
||||
<em>adiam</em> args = type diam
|
||||
type = atom type or range of types (see below)
|
||||
diam = diameter of atoms of that type (distance units)
|
||||
<em>amap</em> args = lo hi style delta N entry1 entry2 ... entryN
|
||||
lo = number or <em>min</em> = lower bound of range of color map
|
||||
hi = number or <em>max</em> = upper bound of range of color map
|
||||
style = 2 letters = "c" or "d" or "s" plus "a" or "f"
|
||||
"c" for continuous
|
||||
"d" for discrete
|
||||
"s" for sequential
|
||||
"a" for absolute
|
||||
"f" for fractional
|
||||
delta = binsize (only used for style "s", otherwise ignored)
|
||||
binsize = range is divided into bins of this width
|
||||
N = # of subsequent entries
|
||||
entry = value color (for continuous style)
|
||||
value = number or <em>min</em> or <em>max</em> = single value within range
|
||||
color = name of color used for that value
|
||||
entry = lo hi color (for discrete style)
|
||||
lo/hi = number or <em>min</em> or <em>max</em> = lower/upper bound of subset of range
|
||||
color = name of color used for that subset of values
|
||||
entry = color (for sequential style)
|
||||
color = name of color used for a bin of values
|
||||
<em>backcolor</em> arg = color
|
||||
color = name of color for background
|
||||
<em>bcolor</em> args = type color
|
||||
type = bond type or range of types (see below)
|
||||
color = name of color or color1/color2/...
|
||||
<em>bdiam</em> args = type diam
|
||||
type = bond type or range of types (see below)
|
||||
diam = diameter of bonds of that type (distance units)
|
||||
<em>boxcolor</em> arg = color
|
||||
color = name of color for simulation box lines and processor sub-domain lines
|
||||
<em>color</em> args = name R G B
|
||||
name = name of color
|
||||
R,G,B = red/green/blue numeric values from 0.0 to 1.0
|
||||
<em>bitrate</em> arg = rate
|
||||
rate = target bitrate for movie in kbps
|
||||
<em>framerate</em> arg = fps
|
||||
fps = frames per second for movie
|
||||
type = atom type or range of types (see below)
|
||||
color = name of color or color1/color2/...
|
||||
<em>adiam</em> args = type diam
|
||||
type = atom type or range of types (see below)
|
||||
diam = diameter of atoms of that type (distance units)
|
||||
<em>amap</em> args = lo hi style delta N entry1 entry2 ... entryN
|
||||
lo = number or <em>min</em> = lower bound of range of color map
|
||||
hi = number or <em>max</em> = upper bound of range of color map
|
||||
style = 2 letters = "c" or "d" or "s" plus "a" or "f"
|
||||
"c" for continuous
|
||||
"d" for discrete
|
||||
"s" for sequential
|
||||
"a" for absolute
|
||||
"f" for fractional
|
||||
delta = binsize (only used for style "s", otherwise ignored)
|
||||
binsize = range is divided into bins of this width
|
||||
N = # of subsequent entries
|
||||
entry = value color (for continuous style)
|
||||
value = number or <em>min</em> or <em>max</em> = single value within range
|
||||
color = name of color used for that value
|
||||
entry = lo hi color (for discrete style)
|
||||
lo/hi = number or <em>min</em> or <em>max</em> = lower/upper bound of subset of range
|
||||
color = name of color used for that subset of values
|
||||
entry = color (for sequential style)
|
||||
color = name of color used for a bin of values
|
||||
<em>backcolor</em> arg = color
|
||||
color = name of color for background
|
||||
<em>bcolor</em> args = type color
|
||||
type = bond type or range of types (see below)
|
||||
color = name of color or color1/color2/...
|
||||
<em>bdiam</em> args = type diam
|
||||
type = bond type or range of types (see below)
|
||||
diam = diameter of bonds of that type (distance units)
|
||||
<em>boxcolor</em> arg = color
|
||||
color = name of color for simulation box lines and processor sub-domain lines
|
||||
<em>color</em> args = name R G B
|
||||
name = name of color
|
||||
R,G,B = red/green/blue numeric values from 0.0 to 1.0
|
||||
<em>bitrate</em> arg = rate
|
||||
rate = target bitrate for movie in kbps
|
||||
<em>framerate</em> arg = fps
|
||||
fps = frames per second for movie
|
||||
</pre>
|
||||
</div>
|
||||
<div class="section" id="examples">
|
||||
|
||||
@ -140,15 +140,15 @@
|
||||
</ul>
|
||||
<pre class="literal-block">
|
||||
<em>pair</em> args = pstyle pparam I J v_name
|
||||
pstyle = pair style name, e.g. lj/cut
|
||||
pparam = parameter to adapt over time
|
||||
I,J = type pair(s) to set parameter for
|
||||
v_name = variable with name that calculates value of pparam
|
||||
<em>kspace</em> arg = v_name
|
||||
v_name = variable with name that calculates scale factor on K-space terms
|
||||
<em>atom</em> args = aparam v_name
|
||||
aparam = parameter to adapt over time
|
||||
v_name = variable with name that calculates value of aparam
|
||||
pstyle = pair style name, e.g. lj/cut
|
||||
pparam = parameter to adapt over time
|
||||
I,J = type pair(s) to set parameter for
|
||||
v_name = variable with name that calculates value of pparam
|
||||
<em>kspace</em> arg = v_name
|
||||
v_name = variable with name that calculates scale factor on K-space terms
|
||||
<em>atom</em> args = aparam v_name
|
||||
aparam = parameter to adapt over time
|
||||
v_name = variable with name that calculates value of aparam
|
||||
</pre>
|
||||
<ul class="simple">
|
||||
<li>zero or more keyword/value pairs may be appended</li>
|
||||
@ -156,11 +156,11 @@
|
||||
</ul>
|
||||
<pre class="literal-block">
|
||||
<em>scale</em> value = <em>no</em> or <em>yes</em>
|
||||
<em>no</em> = the variable value is the new setting
|
||||
<em>yes</em> = the variable value multiplies the original setting
|
||||
<em>reset</em> value = <em>no</em> or <em>yes</em>
|
||||
<em>no</em> = values will remain altered at the end of a run
|
||||
<em>yes</em> = reset altered values to their original values at the end of a run
|
||||
<em>no</em> = the variable value is the new setting
|
||||
<em>yes</em> = the variable value multiplies the original setting
|
||||
<em>reset</em> value = <em>no</em> or <em>yes</em>
|
||||
<em>no</em> = values will remain altered at the end of a run
|
||||
<em>yes</em> = reset altered values to their original values at the end of a run
|
||||
</pre>
|
||||
</div>
|
||||
<div class="section" id="examples">
|
||||
|
||||
@ -140,16 +140,16 @@
|
||||
</ul>
|
||||
<pre class="literal-block">
|
||||
<em>pair</em> args = pstyle pparam I J v_name
|
||||
pstyle = pair style name, e.g. lj/cut
|
||||
pparam = parameter to adapt over time
|
||||
I,J = type pair(s) to set parameter for
|
||||
v_name = variable with name that calculates value of pparam
|
||||
<em>kspace</em> arg = v_name
|
||||
v_name = variable with name that calculates scale factor on K-space terms
|
||||
<em>atom</em> args = aparam v_name
|
||||
aparam = parameter to adapt over time
|
||||
I = type(s) to set parameter for
|
||||
v_name = variable with name that calculates value of aparam
|
||||
pstyle = pair style name, e.g. lj/cut
|
||||
pparam = parameter to adapt over time
|
||||
I,J = type pair(s) to set parameter for
|
||||
v_name = variable with name that calculates value of pparam
|
||||
<em>kspace</em> arg = v_name
|
||||
v_name = variable with name that calculates scale factor on K-space terms
|
||||
<em>atom</em> args = aparam v_name
|
||||
aparam = parameter to adapt over time
|
||||
I = type(s) to set parameter for
|
||||
v_name = variable with name that calculates value of aparam
|
||||
</pre>
|
||||
<ul class="simple">
|
||||
<li>zero or more keyword/value pairs may be appended</li>
|
||||
@ -157,15 +157,15 @@
|
||||
</ul>
|
||||
<pre class="literal-block">
|
||||
<em>scale</em> value = <em>no</em> or <em>yes</em>
|
||||
<em>no</em> = the variable value is the new setting
|
||||
<em>yes</em> = the variable value multiplies the original setting
|
||||
<em>reset</em> value = <em>no</em> or <em>yes</em>
|
||||
<em>no</em> = values will remain altered at the end of a run
|
||||
<em>yes</em> = reset altered values to their original values at the end
|
||||
of a run
|
||||
<em>after</em> value = <em>no</em> or <em>yes</em>
|
||||
<em>no</em> = parameters are adapted at timestep N
|
||||
<em>yes</em> = parameters are adapted one timestep after N
|
||||
<em>no</em> = the variable value is the new setting
|
||||
<em>yes</em> = the variable value multiplies the original setting
|
||||
<em>reset</em> value = <em>no</em> or <em>yes</em>
|
||||
<em>no</em> = values will remain altered at the end of a run
|
||||
<em>yes</em> = reset altered values to their original values at the end
|
||||
of a run
|
||||
<em>after</em> value = <em>no</em> or <em>yes</em>
|
||||
<em>no</em> = parameters are adapted at timestep N
|
||||
<em>yes</em> = parameters are adapted one timestep after N
|
||||
</pre>
|
||||
</div>
|
||||
<div class="section" id="examples">
|
||||
|
||||
@ -148,11 +148,11 @@
|
||||
</ul>
|
||||
<pre class="literal-block">
|
||||
<em>every</em> value = Nevery
|
||||
Nevery = add force every this many timesteps
|
||||
<em>region</em> value = region-ID
|
||||
region-ID = ID of region atoms must be in to have added force
|
||||
<em>energy</em> value = v_name
|
||||
v_name = variable with name that calculates the potential energy of each atom in the added force field
|
||||
Nevery = add force every this many timesteps
|
||||
<em>region</em> value = region-ID
|
||||
region-ID = ID of region atoms must be in to have added force
|
||||
<em>energy</em> value = v_name
|
||||
v_name = variable with name that calculates the potential energy of each atom in the added force field
|
||||
</pre>
|
||||
</div>
|
||||
<div class="section" id="examples">
|
||||
|
||||
@ -140,23 +140,23 @@
|
||||
</ul>
|
||||
<pre class="literal-block">
|
||||
<em>basis</em> values = M itype
|
||||
M = which basis atom
|
||||
itype = atom type (1-N) to assign to this basis atom
|
||||
<em>size</em> args = Lz
|
||||
Lz = z size of lattice region appended in a single event(distance units)
|
||||
<em>freq</em> args = freq
|
||||
freq = the number of timesteps between append events
|
||||
<em>temp</em> args = target damp seed extent
|
||||
target = target temperature for the region between zhi-extent and zhi (temperature units)
|
||||
damp = damping parameter (time units)
|
||||
seed = random number seed for langevin kicks
|
||||
extent = extent of thermostated region (distance units)
|
||||
<em>random</em> args = xmax ymax zmax seed
|
||||
<em>xmax</em>, <em>ymax</em>, <em>zmax</em> = maximum displacement in particular direction (distance units)
|
||||
<em>seed</em> = random number seed for random displacement
|
||||
<em>units</em> value = <em>lattice</em> or <em>box</em>
|
||||
<em>lattice</em> = the wall position is defined in lattice units
|
||||
<em>box</em> = the wall position is defined in simulation box units
|
||||
M = which basis atom
|
||||
itype = atom type (1-N) to assign to this basis atom
|
||||
<em>size</em> args = Lz
|
||||
Lz = z size of lattice region appended in a single event(distance units)
|
||||
<em>freq</em> args = freq
|
||||
freq = the number of timesteps between append events
|
||||
<em>temp</em> args = target damp seed extent
|
||||
target = target temperature for the region between zhi-extent and zhi (temperature units)
|
||||
damp = damping parameter (time units)
|
||||
seed = random number seed for langevin kicks
|
||||
extent = extent of thermostated region (distance units)
|
||||
<em>random</em> args = xmax ymax zmax seed
|
||||
<em>xmax</em>, <em>ymax</em>, <em>zmax</em> = maximum displacement in particular direction (distance units)
|
||||
<em>seed</em> = random number seed for random displacement
|
||||
<em>units</em> value = <em>lattice</em> or <em>box</em>
|
||||
<em>lattice</em> = the wall position is defined in lattice units
|
||||
<em>box</em> = the wall position is defined in simulation box units
|
||||
</pre>
|
||||
</div>
|
||||
<div class="section" id="examples">
|
||||
|
||||
@ -138,9 +138,9 @@
|
||||
</ul>
|
||||
<pre class="literal-block">
|
||||
<em>thermal</em> = thermal coupling with fields: temperature
|
||||
<em>two_temperature</em> = electron-phonon coupling with field: temperature and electron_temperature
|
||||
<em>hardy</em> = on-the-fly post-processing using kernel localization functions (see "related" section for possible fields)
|
||||
<em>field</em> = on-the-fly post-processing using mesh-based localization functions (see "related" section for possible fields)
|
||||
<em>two_temperature</em> = electron-phonon coupling with field: temperature and electron_temperature
|
||||
<em>hardy</em> = on-the-fly post-processing using kernel localization functions (see "related" section for possible fields)
|
||||
<em>field</em> = on-the-fly post-processing using mesh-based localization functions (see "related" section for possible fields)
|
||||
</pre>
|
||||
<ul class="simple">
|
||||
<li>parameter_file = name of the file with material parameters. Note: Neither hardy nor field requires a parameter file</li>
|
||||
@ -164,28 +164,28 @@ fix AtC internal atc field
|
||||
</pre></div>
|
||||
</div>
|
||||
<div class="highlight-python"><div class="highlight"><pre># initial fix to designate coupling type and group to apply it to
|
||||
# tag group physics material_file
|
||||
fix AtC internal atc thermal Ar_thermal.mat
|
||||
# tag group physics material_file
|
||||
fix AtC internal atc thermal Ar_thermal.mat
|
||||
</pre></div>
|
||||
</div>
|
||||
<div class="highlight-python"><div class="highlight"><pre># create a uniform 12 x 2 x 2 mesh that covers region contain the group
|
||||
# nx ny nz region periodicity
|
||||
fix_modify AtC mesh create 12 2 2 mdRegion f p p
|
||||
# nx ny nz region periodicity
|
||||
fix_modify AtC mesh create 12 2 2 mdRegion f p p
|
||||
</pre></div>
|
||||
</div>
|
||||
<div class="highlight-python"><div class="highlight"><pre># specify the control method for the type of coupling
|
||||
# physics control_type
|
||||
fix_modify AtC thermal control flux
|
||||
# physics control_type
|
||||
fix_modify AtC thermal control flux
|
||||
</pre></div>
|
||||
</div>
|
||||
<div class="highlight-python"><div class="highlight"><pre># specify the initial values for the empirical field "temperature"
|
||||
# field node_group value
|
||||
fix_modify AtC initial temperature all 30
|
||||
# field node_group value
|
||||
fix_modify AtC initial temperature all 30
|
||||
</pre></div>
|
||||
</div>
|
||||
<div class="highlight-python"><div class="highlight"><pre># create an output stream for nodal fields
|
||||
# filename output_frequency
|
||||
fix_modify AtC output atc_fe_output 100
|
||||
# filename output_frequency
|
||||
fix_modify AtC output atc_fe_output 100
|
||||
</pre></div>
|
||||
</div>
|
||||
<div class="highlight-python"><div class="highlight"><pre>run 1000
|
||||
@ -196,32 +196,32 @@ fix AtC internal atc field
|
||||
</pre></div>
|
||||
</div>
|
||||
<div class="highlight-python"><div class="highlight"><pre># initial fix to designate post-processing and the group to apply it to
|
||||
# no material file is allowed nor required
|
||||
fix AtC internal atc hardy
|
||||
# no material file is allowed nor required
|
||||
fix AtC internal atc hardy
|
||||
</pre></div>
|
||||
</div>
|
||||
<div class="highlight-python"><div class="highlight"><pre># for hardy fix, specific kernel function (function type and range) to # be used as a localization function
|
||||
fix AtC kernel quartic_sphere 10.0
|
||||
fix AtC kernel quartic_sphere 10.0
|
||||
</pre></div>
|
||||
</div>
|
||||
<div class="highlight-python"><div class="highlight"><pre># create a uniform 1 x 1 x 1 mesh that covers region contain the group
|
||||
# with periodicity this effectively creats a system average
|
||||
fix_modify AtC mesh create 1 1 1 box p p p
|
||||
# with periodicity this effectively creats a system average
|
||||
fix_modify AtC mesh create 1 1 1 box p p p
|
||||
</pre></div>
|
||||
</div>
|
||||
<div class="highlight-python"><div class="highlight"><pre># change from default lagrangian map to eulerian
|
||||
# refreshed every 100 steps
|
||||
fix_modify AtC atom_element_map eulerian 100
|
||||
# refreshed every 100 steps
|
||||
fix_modify AtC atom_element_map eulerian 100
|
||||
</pre></div>
|
||||
</div>
|
||||
<div class="highlight-python"><div class="highlight"><pre># start with no field defined
|
||||
# add mass density, potential energy density, stress and temperature
|
||||
fix_modify AtC fields add density energy stress temperature
|
||||
# add mass density, potential energy density, stress and temperature
|
||||
fix_modify AtC fields add density energy stress temperature
|
||||
</pre></div>
|
||||
</div>
|
||||
<div class="highlight-python"><div class="highlight"><pre># create an output stream for nodal fields
|
||||
# filename output_frequency
|
||||
fix_modify AtC output nvtFE 100 text
|
||||
# filename output_frequency
|
||||
fix_modify AtC output nvtFE 100 text
|
||||
</pre></div>
|
||||
</div>
|
||||
<div class="highlight-python"><div class="highlight"><pre>run 1000
|
||||
|
||||
@ -143,15 +143,15 @@
|
||||
</ul>
|
||||
<pre class="literal-block">
|
||||
<em>types</em> values = two or more atom types
|
||||
<em>delta_mu</em> values = number_of_types-1 relative chemical potentials (energy units)
|
||||
<em>ke</em> value = <em>no</em> or <em>yes</em>
|
||||
<em>no</em> = no conservation of kinetic energy after atom swaps
|
||||
<em>yes</em> = kinetic energy is conserved after atom swaps
|
||||
<em>semi-grand</em> value = <em>no</em> or <em>yes</em>
|
||||
<em>no</em> = particle type counts and fractions conserved
|
||||
<em>yes</em> = semi-grand canonical ensemble, particle fractions not conserved
|
||||
<em>region</em> value = region-ID
|
||||
region-ID = ID of region to use as an exchange/move volume
|
||||
<em>delta_mu</em> values = number_of_types-1 relative chemical potentials (energy units)
|
||||
<em>ke</em> value = <em>no</em> or <em>yes</em>
|
||||
<em>no</em> = no conservation of kinetic energy after atom swaps
|
||||
<em>yes</em> = kinetic energy is conserved after atom swaps
|
||||
<em>semi-grand</em> value = <em>no</em> or <em>yes</em>
|
||||
<em>no</em> = particle type counts and fractions conserved
|
||||
<em>yes</em> = semi-grand canonical ensemble, particle fractions not conserved
|
||||
<em>region</em> value = region-ID
|
||||
region-ID = ID of region to use as an exchange/move volume
|
||||
</pre>
|
||||
</div>
|
||||
<div class="section" id="examples">
|
||||
|
||||
@ -141,11 +141,11 @@ one or more input values can be listed</li>
|
||||
<li>value = x, y, z, vx, vy, vz, fx, fy, fz, c_ID, c_ID[i], f_ID, f_ID[i], v_name</li>
|
||||
</ul>
|
||||
<div class="highlight-python"><div class="highlight"><pre>x,y,z,vx,vy,vz,fx,fy,fz = atom attribute (position, velocity, force component)
|
||||
c_ID = per-atom vector calculated by a compute with ID
|
||||
c_ID[I] = Ith column of per-atom array calculated by a compute with ID
|
||||
f_ID = per-atom vector calculated by a fix with ID
|
||||
f_ID[I] = Ith column of per-atom array calculated by a fix with ID
|
||||
v_name = per-atom vector calculated by an atom-style variable with name
|
||||
c_ID = per-atom vector calculated by a compute with ID
|
||||
c_ID[I] = Ith column of per-atom array calculated by a compute with ID
|
||||
f_ID = per-atom vector calculated by a fix with ID
|
||||
f_ID[I] = Ith column of per-atom array calculated by a fix with ID
|
||||
v_name = per-atom vector calculated by an atom-style variable with name
|
||||
</pre></div>
|
||||
</div>
|
||||
</div>
|
||||
|
||||
@ -142,13 +142,13 @@
|
||||
<li>value = vx, vy, vz, fx, fy, fz, density/mass, density/number, temp, c_ID, c_ID[I], f_ID, f_ID[I], v_name</li>
|
||||
</ul>
|
||||
<div class="highlight-python"><div class="highlight"><pre>vx,vy,vz,fx,fy,fz = atom attribute (velocity, force component)
|
||||
density/number, density/mass = number or mass density
|
||||
temp = temperature
|
||||
c_ID = per-atom vector calculated by a compute with ID
|
||||
c_ID[I] = Ith column of per-atom array calculated by a compute with ID
|
||||
f_ID = per-atom vector calculated by a fix with ID
|
||||
f_ID[I] = Ith column of per-atom array calculated by a fix with ID
|
||||
v_name = per-atom vector calculated by an atom-style variable with name
|
||||
density/number, density/mass = number or mass density
|
||||
temp = temperature
|
||||
c_ID = per-atom vector calculated by a compute with ID
|
||||
c_ID[I] = Ith column of per-atom array calculated by a compute with ID
|
||||
f_ID = per-atom vector calculated by a fix with ID
|
||||
f_ID[I] = Ith column of per-atom array calculated by a fix with ID
|
||||
v_name = per-atom vector calculated by an atom-style variable with name
|
||||
</pre></div>
|
||||
</div>
|
||||
<ul class="simple">
|
||||
@ -157,28 +157,28 @@
|
||||
</ul>
|
||||
<pre class="literal-block">
|
||||
<em>norm</em> arg = <em>all</em> or <em>sample</em> or <em>none</em> = how output on <em>Nfreq</em> steps is normalized
|
||||
all = output is sum of atoms across all <em>Nrepeat</em> samples, divided by atom count
|
||||
sample = output is sum of <em>Nrepeat</em> sample averages, divided by <em>Nrepeat</em>
|
||||
none = output is sum of <em>Nrepeat</em> sums, divided by <em>Nrepeat</em>
|
||||
<em>ave</em> args = <em>one</em> or <em>running</em> or <em>window M</em>
|
||||
one = output new average value every Nfreq steps
|
||||
running = output cumulative average of all previous Nfreq steps
|
||||
window M = output average of M most recent Nfreq steps
|
||||
<em>bias</em> arg = bias-ID
|
||||
bias-ID = ID of a temperature compute that removes a velocity bias for temperature calculation
|
||||
<em>adof</em> value = dof_per_atom
|
||||
dof_per_atom = define this many degrees-of-freedom per atom for temperature calculation
|
||||
<em>cdof</em> value = dof_per_chunk
|
||||
dof_per_chunk = define this many degrees-of-freedom per chunk for temperature calculation
|
||||
<em>file</em> arg = filename
|
||||
filename = file to write results to
|
||||
<em>overwrite</em> arg = none = overwrite output file with only latest output
|
||||
<em>title1</em> arg = string
|
||||
string = text to print as 1st line of output file
|
||||
<em>title2</em> arg = string
|
||||
string = text to print as 2nd line of output file
|
||||
<em>title3</em> arg = string
|
||||
string = text to print as 3rd line of output file
|
||||
all = output is sum of atoms across all <em>Nrepeat</em> samples, divided by atom count
|
||||
sample = output is sum of <em>Nrepeat</em> sample averages, divided by <em>Nrepeat</em>
|
||||
none = output is sum of <em>Nrepeat</em> sums, divided by <em>Nrepeat</em>
|
||||
<em>ave</em> args = <em>one</em> or <em>running</em> or <em>window M</em>
|
||||
one = output new average value every Nfreq steps
|
||||
running = output cumulative average of all previous Nfreq steps
|
||||
window M = output average of M most recent Nfreq steps
|
||||
<em>bias</em> arg = bias-ID
|
||||
bias-ID = ID of a temperature compute that removes a velocity bias for temperature calculation
|
||||
<em>adof</em> value = dof_per_atom
|
||||
dof_per_atom = define this many degrees-of-freedom per atom for temperature calculation
|
||||
<em>cdof</em> value = dof_per_chunk
|
||||
dof_per_chunk = define this many degrees-of-freedom per chunk for temperature calculation
|
||||
<em>file</em> arg = filename
|
||||
filename = file to write results to
|
||||
<em>overwrite</em> arg = none = overwrite output file with only latest output
|
||||
<em>title1</em> arg = string
|
||||
string = text to print as 1st line of output file
|
||||
<em>title2</em> arg = string
|
||||
string = text to print as 2nd line of output file
|
||||
<em>title3</em> arg = string
|
||||
string = text to print as 3rd line of output file
|
||||
</pre>
|
||||
</div>
|
||||
<div class="section" id="examples">
|
||||
|
||||
@ -141,10 +141,10 @@
|
||||
<li>value = c_ID, c_ID[N], f_ID, f_ID[N], v_name</li>
|
||||
</ul>
|
||||
<div class="highlight-python"><div class="highlight"><pre>c_ID = global scalar calculated by a compute with ID
|
||||
c_ID[I] = Ith component of global vector calculated by a compute with ID
|
||||
f_ID = global scalar calculated by a fix with ID
|
||||
f_ID[I] = Ith component of global vector calculated by a fix with ID
|
||||
v_name = global value calculated by an equal-style variable with name
|
||||
c_ID[I] = Ith component of global vector calculated by a compute with ID
|
||||
f_ID = global scalar calculated by a fix with ID
|
||||
f_ID[I] = Ith component of global vector calculated by a fix with ID
|
||||
v_name = global value calculated by an equal-style variable with name
|
||||
</pre></div>
|
||||
</div>
|
||||
<ul class="simple">
|
||||
@ -153,28 +153,28 @@
|
||||
</ul>
|
||||
<pre class="literal-block">
|
||||
<em>type</em> arg = <em>auto</em> or <em>upper</em> or <em>lower</em> or <em>auto/upper</em> or <em>auto/lower</em> or <em>full</em>
|
||||
auto = correlate each value with itself
|
||||
upper = correlate each value with each succeeding value
|
||||
lower = correlate each value with each preceding value
|
||||
auto/upper = auto + upper
|
||||
auto/lower = auto + lower
|
||||
full = correlate each value with every other value, including itself = auto + upper + lower
|
||||
<em>ave</em> args = <em>one</em> or <em>running</em>
|
||||
one = zero the correlation accumulation every Nfreq steps
|
||||
running = accumulate correlations continuously
|
||||
<em>start</em> args = Nstart
|
||||
Nstart = start accumulating correlations on this timestep
|
||||
<em>prefactor</em> args = value
|
||||
value = prefactor to scale all the correlation data by
|
||||
<em>file</em> arg = filename
|
||||
filename = name of file to output correlation data to
|
||||
<em>overwrite</em> arg = none = overwrite output file with only latest output
|
||||
<em>title1</em> arg = string
|
||||
string = text to print as 1st line of output file
|
||||
<em>title2</em> arg = string
|
||||
string = text to print as 2nd line of output file
|
||||
<em>title3</em> arg = string
|
||||
string = text to print as 3rd line of output file
|
||||
auto = correlate each value with itself
|
||||
upper = correlate each value with each succeeding value
|
||||
lower = correlate each value with each preceding value
|
||||
auto/upper = auto + upper
|
||||
auto/lower = auto + lower
|
||||
full = correlate each value with every other value, including itself = auto + upper + lower
|
||||
<em>ave</em> args = <em>one</em> or <em>running</em>
|
||||
one = zero the correlation accumulation every Nfreq steps
|
||||
running = accumulate correlations continuously
|
||||
<em>start</em> args = Nstart
|
||||
Nstart = start accumulating correlations on this timestep
|
||||
<em>prefactor</em> args = value
|
||||
value = prefactor to scale all the correlation data by
|
||||
<em>file</em> arg = filename
|
||||
filename = name of file to output correlation data to
|
||||
<em>overwrite</em> arg = none = overwrite output file with only latest output
|
||||
<em>title1</em> arg = string
|
||||
string = text to print as 1st line of output file
|
||||
<em>title2</em> arg = string
|
||||
string = text to print as 2nd line of output file
|
||||
<em>title3</em> arg = string
|
||||
string = text to print as 3rd line of output file
|
||||
</pre>
|
||||
</div>
|
||||
<div class="section" id="examples">
|
||||
|
||||
@ -146,11 +146,11 @@
|
||||
<li>value = x, y, z, vx, vy, vz, fx, fy, fz, c_ID, c_ID[N], f_ID, f_ID[N], v_name</li>
|
||||
</ul>
|
||||
<div class="highlight-python"><div class="highlight"><pre>x,y,z,vx,vy,vz,fx,fy,fz = atom attribute (position, velocity, force component)
|
||||
c_ID = scalar or vector calculated by a compute with ID
|
||||
c_ID[I] = Ith component of vector or Ith column of array calculated by a compute with ID
|
||||
f_ID = scalar or vector calculated by a fix with ID
|
||||
f_ID[I] = Ith component of vector or Ith column of array calculated by a fix with ID
|
||||
v_name = value(s) calculated by an equal-style or atom-style variable with name
|
||||
c_ID = scalar or vector calculated by a compute with ID
|
||||
c_ID[I] = Ith component of vector or Ith column of array calculated by a compute with ID
|
||||
f_ID = scalar or vector calculated by a fix with ID
|
||||
f_ID[I] = Ith component of vector or Ith column of array calculated by a fix with ID
|
||||
v_name = value(s) calculated by an equal-style or atom-style variable with name
|
||||
</pre></div>
|
||||
</div>
|
||||
<ul class="simple">
|
||||
@ -159,27 +159,27 @@
|
||||
</ul>
|
||||
<pre class="literal-block">
|
||||
<em>mode</em> arg = <em>scalar</em> or <em>vector</em>
|
||||
scalar = all input values are scalars
|
||||
vector = all input values are vectors
|
||||
<em>file</em> arg = filename
|
||||
filename = name of file to output histogram(s) to
|
||||
<em>ave</em> args = <em>one</em> or <em>running</em> or <em>window</em>
|
||||
one = output a new average value every Nfreq steps
|
||||
running = output cumulative average of all previous Nfreq steps
|
||||
window M = output average of M most recent Nfreq steps
|
||||
<em>start</em> args = Nstart
|
||||
Nstart = start averaging on this timestep
|
||||
<em>beyond</em> arg = <em>ignore</em> or <em>end</em> or <em>extra</em>
|
||||
ignore = ignore values outside histogram lo/hi bounds
|
||||
end = count values outside histogram lo/hi bounds in end bins
|
||||
extra = create 2 extra bins for value outside histogram lo/hi bounds
|
||||
<em>overwrite</em> arg = none = overwrite output file with only latest output
|
||||
<em>title1</em> arg = string
|
||||
string = text to print as 1st line of output file
|
||||
<em>title2</em> arg = string
|
||||
string = text to print as 2nd line of output file
|
||||
<em>title3</em> arg = string
|
||||
string = text to print as 3rd line of output file, only for vector mode
|
||||
scalar = all input values are scalars
|
||||
vector = all input values are vectors
|
||||
<em>file</em> arg = filename
|
||||
filename = name of file to output histogram(s) to
|
||||
<em>ave</em> args = <em>one</em> or <em>running</em> or <em>window</em>
|
||||
one = output a new average value every Nfreq steps
|
||||
running = output cumulative average of all previous Nfreq steps
|
||||
window M = output average of M most recent Nfreq steps
|
||||
<em>start</em> args = Nstart
|
||||
Nstart = start averaging on this timestep
|
||||
<em>beyond</em> arg = <em>ignore</em> or <em>end</em> or <em>extra</em>
|
||||
ignore = ignore values outside histogram lo/hi bounds
|
||||
end = count values outside histogram lo/hi bounds in end bins
|
||||
extra = create 2 extra bins for value outside histogram lo/hi bounds
|
||||
<em>overwrite</em> arg = none = overwrite output file with only latest output
|
||||
<em>title1</em> arg = string
|
||||
string = text to print as 1st line of output file
|
||||
<em>title2</em> arg = string
|
||||
string = text to print as 2nd line of output file
|
||||
<em>title3</em> arg = string
|
||||
string = text to print as 3rd line of output file, only for vector mode
|
||||
</pre>
|
||||
</div>
|
||||
<div class="section" id="examples">
|
||||
|
||||
@ -141,20 +141,20 @@
|
||||
</ul>
|
||||
<pre class="literal-block">
|
||||
dim = <em>x</em> or <em>y</em> or <em>z</em>
|
||||
origin = <em>lower</em> or <em>center</em> or <em>upper</em> or coordinate value (distance units)
|
||||
delta = thickness of spatial bins in dim (distance units)
|
||||
origin = <em>lower</em> or <em>center</em> or <em>upper</em> or coordinate value (distance units)
|
||||
delta = thickness of spatial bins in dim (distance units)
|
||||
</pre>
|
||||
<ul class="simple">
|
||||
<li>one or more input values can be listed</li>
|
||||
<li>value = vx, vy, vz, fx, fy, fz, density/mass, density/number, c_ID, c_ID[I], f_ID, f_ID[I], v_name</li>
|
||||
</ul>
|
||||
<div class="highlight-python"><div class="highlight"><pre>vx,vy,vz,fx,fy,fz = atom attribute (velocity, force component)
|
||||
density/number, density/mass = number or mass density
|
||||
c_ID = per-atom vector calculated by a compute with ID
|
||||
c_ID[I] = Ith column of per-atom array calculated by a compute with ID
|
||||
f_ID = per-atom vector calculated by a fix with ID
|
||||
f_ID[I] = Ith column of per-atom array calculated by a fix with ID
|
||||
v_name = per-atom vector calculated by an atom-style variable with name
|
||||
density/number, density/mass = number or mass density
|
||||
c_ID = per-atom vector calculated by a compute with ID
|
||||
c_ID[I] = Ith column of per-atom array calculated by a compute with ID
|
||||
f_ID = per-atom vector calculated by a fix with ID
|
||||
f_ID[I] = Ith column of per-atom array calculated by a fix with ID
|
||||
v_name = per-atom vector calculated by an atom-style variable with name
|
||||
</pre></div>
|
||||
</div>
|
||||
<ul class="simple">
|
||||
@ -163,30 +163,30 @@ dim = <em>x</em> or <em>y</em> or <em>z</em>
|
||||
</ul>
|
||||
<pre class="literal-block">
|
||||
<em>region</em> arg = region-ID
|
||||
<em>bound</em> args = x/y/z lo hi
|
||||
x/y/z = <em>x</em> or <em>y</em> or <em>z</em> to bound bins in this dimension
|
||||
lo = <em>lower</em> or coordinate value (distance units)
|
||||
hi = <em>upper</em> or coordinate value (distance units)
|
||||
<em>discard</em> arg = <em>mixed</em> or <em>no</em> or <em>yes</em>
|
||||
mixed = discard atoms outside bins only if bin bounds are explicitly set
|
||||
no = always keep out-of-bounds atoms
|
||||
yes = always discard out-of-bounds atoms
|
||||
<em>norm</em> arg = <em>all</em> or <em>sample</em>
|
||||
region-ID = ID of region atoms must be in to contribute to spatial averaging
|
||||
<em>ave</em> args = <em>one</em> or <em>running</em> or <em>window M</em>
|
||||
one = output new average value every Nfreq steps
|
||||
running = output cumulative average of all previous Nfreq steps
|
||||
window M = output average of M most recent Nfreq steps
|
||||
<em>units</em> arg = <em>box</em> or <em>lattice</em> or <em>reduced</em>
|
||||
<em>file</em> arg = filename
|
||||
filename = file to write results to
|
||||
<em>overwrite</em> arg = none = overwrite output file with only latest output
|
||||
<em>title1</em> arg = string
|
||||
string = text to print as 1st line of output file
|
||||
<em>title2</em> arg = string
|
||||
string = text to print as 2nd line of output file
|
||||
<em>title3</em> arg = string
|
||||
string = text to print as 3rd line of output file
|
||||
<em>bound</em> args = x/y/z lo hi
|
||||
x/y/z = <em>x</em> or <em>y</em> or <em>z</em> to bound bins in this dimension
|
||||
lo = <em>lower</em> or coordinate value (distance units)
|
||||
hi = <em>upper</em> or coordinate value (distance units)
|
||||
<em>discard</em> arg = <em>mixed</em> or <em>no</em> or <em>yes</em>
|
||||
mixed = discard atoms outside bins only if bin bounds are explicitly set
|
||||
no = always keep out-of-bounds atoms
|
||||
yes = always discard out-of-bounds atoms
|
||||
<em>norm</em> arg = <em>all</em> or <em>sample</em>
|
||||
region-ID = ID of region atoms must be in to contribute to spatial averaging
|
||||
<em>ave</em> args = <em>one</em> or <em>running</em> or <em>window M</em>
|
||||
one = output new average value every Nfreq steps
|
||||
running = output cumulative average of all previous Nfreq steps
|
||||
window M = output average of M most recent Nfreq steps
|
||||
<em>units</em> arg = <em>box</em> or <em>lattice</em> or <em>reduced</em>
|
||||
<em>file</em> arg = filename
|
||||
filename = file to write results to
|
||||
<em>overwrite</em> arg = none = overwrite output file with only latest output
|
||||
<em>title1</em> arg = string
|
||||
string = text to print as 1st line of output file
|
||||
<em>title2</em> arg = string
|
||||
string = text to print as 2nd line of output file
|
||||
<em>title3</em> arg = string
|
||||
string = text to print as 3rd line of output file
|
||||
</pre>
|
||||
</div>
|
||||
<div class="section" id="examples">
|
||||
|
||||
@ -145,12 +145,12 @@
|
||||
<li>value = vx, vy, vz, fx, fy, fz, density/mass, density/number, c_ID, c_ID[I], f_ID, f_ID[I], v_name</li>
|
||||
</ul>
|
||||
<div class="highlight-python"><div class="highlight"><pre>vx,vy,vz,fx,fy,fz = atom attribute (velocity, force component)
|
||||
density/number, density/mass = number or mass density
|
||||
c_ID = per-atom vector calculated by a compute with ID
|
||||
c_ID[I] = Ith column of per-atom array calculated by a compute with ID
|
||||
f_ID = per-atom vector calculated by a fix with ID
|
||||
f_ID[I] = Ith column of per-atom array calculated by a fix with ID
|
||||
v_name = per-atom vector calculated by an atom-style variable with name
|
||||
density/number, density/mass = number or mass density
|
||||
c_ID = per-atom vector calculated by a compute with ID
|
||||
c_ID[I] = Ith column of per-atom array calculated by a compute with ID
|
||||
f_ID = per-atom vector calculated by a fix with ID
|
||||
f_ID[I] = Ith column of per-atom array calculated by a fix with ID
|
||||
v_name = per-atom vector calculated by an atom-style variable with name
|
||||
</pre></div>
|
||||
</div>
|
||||
<ul class="simple">
|
||||
@ -159,22 +159,22 @@
|
||||
</ul>
|
||||
<pre class="literal-block">
|
||||
<em>region</em> arg = region-ID
|
||||
region-ID = ID of region atoms must be in to contribute to spatial averaging
|
||||
<em>norm</em> arg = <em>all</em> or <em>sample</em>
|
||||
<em>units</em> arg = <em>box</em> or <em>lattice</em> or <em>reduced</em>
|
||||
<em>ave</em> args = <em>one</em> or <em>running</em> or <em>window M</em>
|
||||
one = output new average value every Nfreq steps
|
||||
running = output cumulative average of all previous Nfreq steps
|
||||
window M = output average of M most recent Nfreq steps
|
||||
<em>file</em> arg = filename
|
||||
filename = file to write results to
|
||||
<em>overwrite</em> arg = none = overwrite output file with only latest output
|
||||
<em>title1</em> arg = string
|
||||
string = text to print as 1st line of output file
|
||||
<em>title2</em> arg = string
|
||||
string = text to print as 2nd line of output file
|
||||
<em>title3</em> arg = string
|
||||
string = text to print as 3rd line of output file
|
||||
region-ID = ID of region atoms must be in to contribute to spatial averaging
|
||||
<em>norm</em> arg = <em>all</em> or <em>sample</em>
|
||||
<em>units</em> arg = <em>box</em> or <em>lattice</em> or <em>reduced</em>
|
||||
<em>ave</em> args = <em>one</em> or <em>running</em> or <em>window M</em>
|
||||
one = output new average value every Nfreq steps
|
||||
running = output cumulative average of all previous Nfreq steps
|
||||
window M = output average of M most recent Nfreq steps
|
||||
<em>file</em> arg = filename
|
||||
filename = file to write results to
|
||||
<em>overwrite</em> arg = none = overwrite output file with only latest output
|
||||
<em>title1</em> arg = string
|
||||
string = text to print as 1st line of output file
|
||||
<em>title2</em> arg = string
|
||||
string = text to print as 2nd line of output file
|
||||
<em>title3</em> arg = string
|
||||
string = text to print as 3rd line of output file
|
||||
</pre>
|
||||
</div>
|
||||
<div class="section" id="examples">
|
||||
|
||||
@ -141,10 +141,10 @@
|
||||
<li>value = c_ID, c_ID[N], f_ID, f_ID[N], v_name</li>
|
||||
</ul>
|
||||
<div class="highlight-python"><div class="highlight"><pre>c_ID = global scalar, vector, or array calculated by a compute with ID
|
||||
c_ID[I] = Ith component of global vector or Ith column of global array calculated by a compute with ID
|
||||
f_ID = global scalar, vector, or array calculated by a fix with ID
|
||||
f_ID[I] = Ith component of global vector or Ith column of global array calculated by a fix with ID
|
||||
v_name = global value calculated by an equal-style variable with name
|
||||
c_ID[I] = Ith component of global vector or Ith column of global array calculated by a compute with ID
|
||||
f_ID = global scalar, vector, or array calculated by a fix with ID
|
||||
f_ID[I] = Ith component of global vector or Ith column of global array calculated by a fix with ID
|
||||
v_name = global value calculated by an equal-style variable with name
|
||||
</pre></div>
|
||||
</div>
|
||||
<ul class="simple">
|
||||
@ -153,25 +153,25 @@
|
||||
</ul>
|
||||
<pre class="literal-block">
|
||||
<em>mode</em> arg = <em>scalar</em> or <em>vector</em>
|
||||
scalar = all input values are global scalars
|
||||
vector = all input values are global vectors or global arrays
|
||||
<em>ave</em> args = <em>one</em> or <em>running</em> or <em>window M</em>
|
||||
one = output a new average value every Nfreq steps
|
||||
running = output cummulative average of all previous Nfreq steps
|
||||
window M = output average of M most recent Nfreq steps
|
||||
<em>start</em> args = Nstart
|
||||
Nstart = start averaging on this timestep
|
||||
<em>off</em> arg = M = do not average this value
|
||||
M = value # from 1 to Nvalues
|
||||
<em>file</em> arg = filename
|
||||
filename = name of file to output time averages to
|
||||
<em>overwrite</em> arg = none = overwrite output file with only latest output
|
||||
<em>title1</em> arg = string
|
||||
string = text to print as 1st line of output file
|
||||
<em>title2</em> arg = string
|
||||
string = text to print as 2nd line of output file
|
||||
<em>title3</em> arg = string
|
||||
string = text to print as 3rd line of output file, only for vector mode
|
||||
scalar = all input values are global scalars
|
||||
vector = all input values are global vectors or global arrays
|
||||
<em>ave</em> args = <em>one</em> or <em>running</em> or <em>window M</em>
|
||||
one = output a new average value every Nfreq steps
|
||||
running = output cummulative average of all previous Nfreq steps
|
||||
window M = output average of M most recent Nfreq steps
|
||||
<em>start</em> args = Nstart
|
||||
Nstart = start averaging on this timestep
|
||||
<em>off</em> arg = M = do not average this value
|
||||
M = value # from 1 to Nvalues
|
||||
<em>file</em> arg = filename
|
||||
filename = name of file to output time averages to
|
||||
<em>overwrite</em> arg = none = overwrite output file with only latest output
|
||||
<em>title1</em> arg = string
|
||||
string = text to print as 1st line of output file
|
||||
<em>title2</em> arg = string
|
||||
string = text to print as 2nd line of output file
|
||||
<em>title3</em> arg = string
|
||||
string = text to print as 3rd line of output file, only for vector mode
|
||||
</pre>
|
||||
</div>
|
||||
<div class="section" id="examples">
|
||||
|
||||
@ -148,7 +148,7 @@
|
||||
</ul>
|
||||
<pre class="literal-block">
|
||||
<em>region</em> value = region-ID
|
||||
region-ID = ID of region atoms must be in to have added force
|
||||
region-ID = ID of region atoms must be in to have added force
|
||||
</pre>
|
||||
</div>
|
||||
<div class="section" id="examples">
|
||||
|
||||
@ -139,10 +139,10 @@
|
||||
<li>style = <em>shift</em> or <em>rcb</em></li>
|
||||
</ul>
|
||||
<div class="highlight-python"><div class="highlight"><pre>shift args = dimstr Niter stopthresh
|
||||
dimstr = sequence of letters containing "x" or "y" or "z", each not more than once
|
||||
Niter = # of times to iterate within each dimension of dimstr sequence
|
||||
stopthresh = stop balancing when this imbalance threshhold is reached
|
||||
rcb args = none
|
||||
dimstr = sequence of letters containing "x" or "y" or "z", each not more than once
|
||||
Niter = # of times to iterate within each dimension of dimstr sequence
|
||||
stopthresh = stop balancing when this imbalance threshhold is reached
|
||||
rcb args = none
|
||||
</pre></div>
|
||||
</div>
|
||||
<ul class="simple">
|
||||
@ -151,7 +151,7 @@
|
||||
</ul>
|
||||
<pre class="literal-block">
|
||||
<em>out</em> value = filename
|
||||
filename = write each processor's sub-domain to a file, at each re-balancing
|
||||
filename = write each processor's sub-domain to a file, at each re-balancing
|
||||
</pre>
|
||||
</div>
|
||||
<div class="section" id="examples">
|
||||
@ -237,7 +237,7 @@ that owns them. The leftmost diagram is the default partitioning of
|
||||
the simulation box across processors (one sub-box for each of 16
|
||||
processors); the middle diagram is after a “grid” method has been
|
||||
applied.</p>
|
||||
<a data-lightbox="group-85269663-2cc5-4073-9d2b-e7fb01db4ec3"
|
||||
<a data-lightbox="group-default"
|
||||
href="_images/balance_uniform.jpg"
|
||||
class=""
|
||||
title=""
|
||||
@ -247,7 +247,7 @@ applied.</p>
|
||||
width="25%"
|
||||
height="auto"
|
||||
alt=""/>
|
||||
</a><a data-lightbox="group-349a1bab-83d0-42b4-85c7-81050f4e3d61"
|
||||
</a><a data-lightbox="group-default"
|
||||
href="_images/balance_nonuniform.jpg"
|
||||
class=""
|
||||
title=""
|
||||
@ -257,7 +257,7 @@ applied.</p>
|
||||
width="25%"
|
||||
height="auto"
|
||||
alt=""/>
|
||||
</a><a data-lightbox="group-c6a44ae9-1930-4d73-95ea-e4e0a56b953a"
|
||||
</a><a data-lightbox="group-default"
|
||||
href="_images/balance_rcb.jpg"
|
||||
class=""
|
||||
title=""
|
||||
|
||||
@ -142,8 +142,8 @@
|
||||
</ul>
|
||||
<pre class="literal-block">
|
||||
<em>prob</em> values = fraction seed
|
||||
fraction = break a bond with this probability if otherwise eligible
|
||||
seed = random number seed (positive integer)
|
||||
fraction = break a bond with this probability if otherwise eligible
|
||||
seed = random number seed (positive integer)
|
||||
</pre>
|
||||
</div>
|
||||
<div class="section" id="examples">
|
||||
|
||||
@ -143,20 +143,20 @@
|
||||
</ul>
|
||||
<pre class="literal-block">
|
||||
<em>iparam</em> values = maxbond, newtype
|
||||
maxbond = max # of bonds of bondtype the itype atom can have
|
||||
newtype = change the itype atom to this type when maxbonds exist
|
||||
<em>jparam</em> values = maxbond, newtype
|
||||
maxbond = max # of bonds of bondtype the jtype atom can have
|
||||
newtype = change the jtype atom to this type when maxbonds exist
|
||||
<em>prob</em> values = fraction seed
|
||||
fraction = create a bond with this probability if otherwise eligible
|
||||
seed = random number seed (positive integer)
|
||||
<em>atype</em> value = angletype
|
||||
angletype = type of created angles
|
||||
<em>dtype</em> value = dihedraltype
|
||||
dihedraltype = type of created dihedrals
|
||||
<em>itype</em> value = impropertype
|
||||
impropertype = type of created impropers
|
||||
maxbond = max # of bonds of bondtype the itype atom can have
|
||||
newtype = change the itype atom to this type when maxbonds exist
|
||||
<em>jparam</em> values = maxbond, newtype
|
||||
maxbond = max # of bonds of bondtype the jtype atom can have
|
||||
newtype = change the jtype atom to this type when maxbonds exist
|
||||
<em>prob</em> values = fraction seed
|
||||
fraction = create a bond with this probability if otherwise eligible
|
||||
seed = random number seed (positive integer)
|
||||
<em>atype</em> value = angletype
|
||||
angletype = type of created angles
|
||||
<em>dtype</em> value = dihedraltype
|
||||
dihedraltype = type of created dihedrals
|
||||
<em>itype</em> value = impropertype
|
||||
impropertype = type of created impropers
|
||||
</pre>
|
||||
</div>
|
||||
<div class="section" id="examples">
|
||||
|
||||
@ -139,11 +139,11 @@
|
||||
</ul>
|
||||
<pre class="literal-block">
|
||||
<em>input</em> arg = colvars.state file name or prefix or NULL (default: NULL)
|
||||
<em>output</em> arg = output filename prefix (default: out)
|
||||
<em>seed</em> arg = seed for random number generator (default: 1966)
|
||||
<em>unwrap</em> arg = <em>yes</em> or <em>no</em>
|
||||
use unwrapped coordinates in collective variables (default: yes)
|
||||
<em>tstat</em> arg = fix id of a thermostat or NULL (default: NULL)
|
||||
<em>output</em> arg = output filename prefix (default: out)
|
||||
<em>seed</em> arg = seed for random number generator (default: 1966)
|
||||
<em>unwrap</em> arg = <em>yes</em> or <em>no</em>
|
||||
use unwrapped coordinates in collective variables (default: yes)
|
||||
<em>tstat</em> arg = fix id of a thermostat or NULL (default: NULL)
|
||||
</pre>
|
||||
</div>
|
||||
<div class="section" id="examples">
|
||||
|
||||
@ -187,14 +187,14 @@ parameter = <em>x</em> or <em>y</em> or <em>z</em> or <em>xy</em> or <em>xz</em>
|
||||
</ul>
|
||||
<pre class="literal-block">
|
||||
<em>remap</em> value = <em>x</em> or <em>v</em> or <em>none</em>
|
||||
x = remap coords of atoms in group into deforming box
|
||||
v = remap velocities of all atoms when they cross periodic boundaries
|
||||
none = no remapping of x or v
|
||||
<em>flip</em> value = <em>yes</em> or <em>no</em>
|
||||
allow or disallow box flips when it becomes highly skewed
|
||||
<em>units</em> value = <em>lattice</em> or <em>box</em>
|
||||
lattice = distances are defined in lattice units
|
||||
box = distances are defined in simulation box units
|
||||
x = remap coords of atoms in group into deforming box
|
||||
v = remap velocities of all atoms when they cross periodic boundaries
|
||||
none = no remapping of x or v
|
||||
<em>flip</em> value = <em>yes</em> or <em>no</em>
|
||||
allow or disallow box flips when it becomes highly skewed
|
||||
<em>units</em> value = <em>lattice</em> or <em>box</em>
|
||||
lattice = distances are defined in lattice units
|
||||
box = distances are defined in simulation box units
|
||||
</pre>
|
||||
</div>
|
||||
<div class="section" id="examples">
|
||||
|
||||
@ -143,40 +143,40 @@
|
||||
</ul>
|
||||
<pre class="literal-block">
|
||||
<em>region</em> value = region-ID
|
||||
region-ID = ID of region to use as insertion volume
|
||||
<em>id</em> value = <em>max</em> or <em>next</em>
|
||||
max = atom ID for new atom(s) is max ID of all current atoms plus one
|
||||
next = atom ID for new atom(s) increments by one for every deposition
|
||||
<em>global</em> values = lo hi
|
||||
lo,hi = put new atom/molecule a distance lo-hi above all other atoms (distance units)
|
||||
<em>local</em> values = lo hi delta
|
||||
lo,hi = put new atom/molecule a distance lo-hi above any nearby atom beneath it (distance units)
|
||||
delta = lateral distance within which a neighbor is considered "nearby" (distance units)
|
||||
<em>near</em> value = R
|
||||
R = only insert atom/molecule if further than R from existing particles (distance units)
|
||||
<em>attempt</em> value = Q
|
||||
Q = attempt a single insertion up to Q times
|
||||
<em>rate</em> value = V
|
||||
V = z velocity (y in 2d) at which insertion volume moves (velocity units)
|
||||
<em>vx</em> values = vxlo vxhi
|
||||
vxlo,vxhi = range of x velocities for inserted atom/molecule (velocity units)
|
||||
<em>vy</em> values = vylo vyhi
|
||||
vylo,vyhi = range of y velocities for inserted atom/molecule (velocity units)
|
||||
<em>vz</em> values = vzlo vzhi
|
||||
vzlo,vzhi = range of z velocities for inserted atom/molecule (velocity units)
|
||||
<em>target</em> values = tx ty tz
|
||||
tx,ty,tz = location of target point (distance units)
|
||||
<em>mol</em> value = template-ID
|
||||
template-ID = ID of molecule template specified in a separate <a class="reference internal" href="molecule.html"><em>molecule</em></a> command
|
||||
<em>molfrac</em> values = f1 f2 ... fN
|
||||
f1 to fN = relative probability of creating each of N molecules in template-ID
|
||||
<em>rigid</em> value = fix-ID
|
||||
fix-ID = ID of <a class="reference internal" href="fix_rigid.html"><em>fix rigid/small</em></a> command
|
||||
<em>shake</em> value = fix-ID
|
||||
fix-ID = ID of <a class="reference internal" href="fix_shake.html"><em>fix shake</em></a> command
|
||||
<em>units</em> value = <em>lattice</em> or <em>box</em>
|
||||
lattice = the geometry is defined in lattice units
|
||||
box = the geometry is defined in simulation box units
|
||||
region-ID = ID of region to use as insertion volume
|
||||
<em>id</em> value = <em>max</em> or <em>next</em>
|
||||
max = atom ID for new atom(s) is max ID of all current atoms plus one
|
||||
next = atom ID for new atom(s) increments by one for every deposition
|
||||
<em>global</em> values = lo hi
|
||||
lo,hi = put new atom/molecule a distance lo-hi above all other atoms (distance units)
|
||||
<em>local</em> values = lo hi delta
|
||||
lo,hi = put new atom/molecule a distance lo-hi above any nearby atom beneath it (distance units)
|
||||
delta = lateral distance within which a neighbor is considered "nearby" (distance units)
|
||||
<em>near</em> value = R
|
||||
R = only insert atom/molecule if further than R from existing particles (distance units)
|
||||
<em>attempt</em> value = Q
|
||||
Q = attempt a single insertion up to Q times
|
||||
<em>rate</em> value = V
|
||||
V = z velocity (y in 2d) at which insertion volume moves (velocity units)
|
||||
<em>vx</em> values = vxlo vxhi
|
||||
vxlo,vxhi = range of x velocities for inserted atom/molecule (velocity units)
|
||||
<em>vy</em> values = vylo vyhi
|
||||
vylo,vyhi = range of y velocities for inserted atom/molecule (velocity units)
|
||||
<em>vz</em> values = vzlo vzhi
|
||||
vzlo,vzhi = range of z velocities for inserted atom/molecule (velocity units)
|
||||
<em>target</em> values = tx ty tz
|
||||
tx,ty,tz = location of target point (distance units)
|
||||
<em>mol</em> value = template-ID
|
||||
template-ID = ID of molecule template specified in a separate <a class="reference internal" href="molecule.html"><em>molecule</em></a> command
|
||||
<em>molfrac</em> values = f1 f2 ... fN
|
||||
f1 to fN = relative probability of creating each of N molecules in template-ID
|
||||
<em>rigid</em> value = fix-ID
|
||||
fix-ID = ID of <a class="reference internal" href="fix_rigid.html"><em>fix rigid/small</em></a> command
|
||||
<em>shake</em> value = fix-ID
|
||||
fix-ID = ID of <a class="reference internal" href="fix_shake.html"><em>fix shake</em></a> command
|
||||
<em>units</em> value = <em>lattice</em> or <em>box</em>
|
||||
lattice = the geometry is defined in lattice units
|
||||
box = the geometry is defined in simulation box units
|
||||
</pre>
|
||||
</div>
|
||||
<div class="section" id="examples">
|
||||
|
||||
@ -143,8 +143,8 @@
|
||||
</ul>
|
||||
<pre class="literal-block">
|
||||
<em>units</em> value = <em>lattice</em> or <em>box</em>
|
||||
lattice = Xmax is defined in lattice units
|
||||
box = Xmax is defined in simulation box units
|
||||
lattice = Xmax is defined in lattice units
|
||||
box = Xmax is defined in simulation box units
|
||||
</pre>
|
||||
</div>
|
||||
<div class="section" id="examples">
|
||||
|
||||