2
<CENTER><A HREF = "http://lammps.sandia.gov">LAMMPS WWW Site</A> - <A HREF = "Manual.html">LAMMPS Documentation</A> - <A HREF = "Section_commands.html#comm">LAMMPS Commands</A>
12
<H3>compute basal/atom command
16
<PRE>compute ID group-ID basal/atom
18
<UL><LI>ID, group-ID are documented in <A HREF = "compute.html">compute</A> command
19
<LI>basal/atom = style name of this compute command
23
<PRE>compute 1 all basal/atom
25
<P><B>Description:</B>
27
<P>Defines a computation that calculates the hexagonal close-packed "c"
28
lattice vector for each atom in the group. It does this by
29
calculating the normal unit vector to the basal plane for each atom.
30
The results enable efficient identification and characterization of
31
twins and grains in hexagonal close-packed structures.
33
<P>The output of the compute is thus the 3 components of a unit vector
34
associdate with each atom. The components are set to 0.0 for
35
atoms not in the group.
37
<P>Details of the calculation are given in <A HREF = "#Barrett">(Barrett)</A>.
39
<P>The neighbor list needed to compute this quantity is constructed each
40
time the calculation is performed (i.e. each time a snapshot of atoms
41
is dumped). Thus it can be inefficient to compute/dump this quantity
42
too frequently or to have multiple compute/dump commands, each of
43
which computes this quantity.
45
<P>An example input script that uses this compute is provided
46
in examples/USER/misc/basal.
48
<P><B>Output info:</B>
50
<P>This compute calculates a per-atom array with 3 columns, which can be
51
accessed by indices 1-3 by any command that uses per-atom values from
52
a compute as input. See <A HREF = "Section_howto.html#howto_15">Section_howto
53
15</A> for an overview of LAMMPS output
56
<P>The per-atom vector values are unitless since the 3 columns represent
57
components of a unit vector.
59
<P><B>Restrictions:</B>
61
<P>This compute is part of the USER-MISC package. It is only enabled if
62
LAMMPS was built with that package. See the <A HREF = "Section_start.html#start_3">Making
63
LAMMPS</A> section for more info.
65
<P>The output of this compute will be meaningless unless the atoms are on
66
(or near) hcp lattice sites, since the calculation assumes a
67
well-defined basal plane.
69
<P><B>Related commands:</B>
71
<P><A HREF = "compute_centro_atom.html">compute centro/atom</A>, <A HREF = "compute_ackland_atom.html">compute
74
<P><B>Default:</B> none
78
<A NAME = "Barrett"></A>
80
<P><B>(Barrett)</B> Barrett, Tschopp, El Kadiri, Scripta Mat. 66, p.666 (2012).