<p>which are fully explained in <a class="reference internal" href="#janssens"><span class="std std-ref">(Janssens)</span></a>. For fcc crystals
this order parameter Xi for atom I in equation (1) is a sum over the
12 nearest neighbors of atom I. For bcc crystals it is the
corresponding sum of the 8 nearest neighbors. Rj is the vector from
atom I to its neighbor J, and RIj is a vector in the reference
(perfect) crystal. That is, if dir = 0/1, then RIj is a vector to an
atom coord from file 0/1. Equation (2) gives the expected value of
the order parameter XiIJ in the other grain. Hi and lo cutoffs are
defined in equations (3) and (4), using the input parameters <em>cutlo</em>
and <em>cuthi</em> as thresholds to avoid adding grain boundary energy when
the deviation in the order parameter from 0 or 1 is small (e.g. due to
thermal fluctuations in a perfect crystal). The added potential
energy Ui for atom I is given in equation (6) where it is interpolated
between 0 and dE using the two threshold Xi values and the Wi value of
equation (5).</p>
<p>The derivative of this energy expression gives the force on each atom
which thus depends on the orientation of its neighbors relative to the
2 grain orientations. Only atoms near the grain boundary feel a net
force which tends to drive them to one of the two grain orientations.</p>
<p>In equation (1), the reference vector used for each neighbor is the
reference vector closest to the actual neighbor position. This means
it is possible two different neighbors will use the same reference
vector. In such cases, the atom in question is far from a perfect
orientation and will likely receive the full dE addition, so the
effect of duplicate reference vector usage is small.</p>
<p>The <em>dir</em> parameter determines which grain wants to grow at the
expense of the other. A value of 0 means the first grain will shrink;
a value of 1 means it will grow. This assumes that <em>dE</em> is positive.
The reverse will be true if <em>dE</em> is negative.</p>
<p>The <em>alat</em> parameter is the cubic lattice constant for the fcc or bcc
material and is only used to compute a cutoff distance of 1.57 * alat
/ sqrt(2) for finding the 12 or 8 nearest neighbors of each atom
(which should be valid for an fcc or bcc crystal). A longer/shorter
cutoff can be imposed by adjusting <em>alat</em>. If a particular atom has
less than 12 or 8 neighbors within the cutoff, the order parameter of
equation (1) is effectively multiplied by 12 or 8 divided by the
actual number of neighbors within the cutoff.</p>
<p>The <em>dE</em> parameter is the maximum amount of additional energy added to
each atom in the grain which wants to shrink.</p>
<p>The <em>cutlo</em> and <em>cuthi</em> parameters are used to reduce the force added
to bulk atoms in each grain far away from the boundary. An atom in
the bulk surrounded by neighbors at the ideal grain orientation would
compute an order parameter of 0 or 1 and have no force added.
However, thermal vibrations in the solid will cause the order
parameters to be greater than 0 or less than 1. The cutoff parameters
mask this effect, allowing forces to only be added to atoms with
order-parameters between the cutoff values.</p>
<p><em>File0</em> and <em>file1</em> are filenames for the two grains which each
contain 6 vectors (6 lines with 3 values per line) which specify the
grain orientations. Each vector is a displacement from a central atom
(0,0,0) to a nearest neighbor atom in an fcc lattice at the proper
orientation. The vector lengths should all be identical since an fcc
lattice has a coordination number of 12. Only 6 are listed due to
symmetry, so the list must include one from each pair of
equal-and-opposite neighbors. A pair of orientation files for a
Sigma=5 tilt boundary are shown below. A tutorial that can help for
writing the orientation files is given in <a class="reference internal" href="#wicaksono2"><span class="std std-ref">(Wicaksono2)</span></a></p>
<p><strong>Restart, fix_modify, output, run start/stop, minimize info:</strong></p>
<p>No information about this fix is written to <a class="reference internal" href="restart.html"><span class="doc">binary restart files</span></a>.</p>
<p>The <a class="reference internal" href="fix_modify.html"><span class="doc">fix_modify</span></a> <em>energy</em> option is supported by this
fix to add the potential energy of atom interactions with the grain
boundary driving force to the system’s potential energy as part of
<p>The <a class="reference internal" href="fix_modify.html"><span class="doc">fix_modify</span></a> <em>respa</em> option is supported by these
fixes. This allows to set at which level of the <a class="reference internal" href="run_style.html"><span class="doc">r-RESPA</span></a>
integrator a fix is adding its forces. Default is the outermost level.</p>
<p>This fix calculates a global scalar which can be accessed by various
<a class="reference internal" href="Section_howto.html#howto-15"><span class="std std-ref">output commands</span></a>. The scalar is the
potential energy change due to this fix. The scalar value calculated
by this fix is “extensive”.</p>
<p>This fix also calculates a per-atom array which can be accessed by
various <a class="reference internal" href="Section_howto.html#howto-15"><span class="std std-ref">output commands</span></a>. The array
stores the order parameter Xi and normalized order parameter (0 to 1)
for each atom. The per-atom values can be accessed on any timestep.</p>
<p>No parameter of this fix can be used with the <em>start/stop</em> keywords of
the <a class="reference internal" href="run.html"><span class="doc">run</span></a> command. This fix is not invoked during <a class="reference internal" href="minimize.html"><span class="doc">energy minimization</span></a>.</p>
</div>
<div class="section" id="restrictions">
<h2>Restrictions</h2>
<p>This fix is part of the MISC package. It is only enabled if LAMMPS
was built with that package. See the <a class="reference internal" href="Section_start.html#start-3"><span class="std std-ref">Making LAMMPS</span></a> section for more info.</p>
<p>This fix should only be used with fcc or bcc lattices.</p>
Built with <a href="http://sphinx-doc.org/">Sphinx</a> using a <a href="https://github.com/snide/sphinx_rtd_theme">theme</a> provided by <a href="https://readthedocs.org">Read the Docs</a>.