<spanid="index-0"></span><h1>compute omega/chunk command<aclass="headerlink"href="#compute-omega-chunk-command"title="Permalink to this headline">¶</a></h1>
<divclass="section"id="syntax">
<h2>Syntax<aclass="headerlink"href="#syntax"title="Permalink to this headline">¶</a></h2>
<divclass="highlight-python"><divclass="highlight"><pre>compute ID group-ID omega/chunk chunkID
</pre></div>
</div>
<ulclass="simple">
<li>ID, group-ID are documented in <aclass="reference internal"href="compute.html"><em>compute</em></a> command</li>
<li>omega/chunk = style name of this compute command</li>
<li>chunkID = ID of <aclass="reference internal"href="compute_chunk_atom.html"><em>compute chunk/atom</em></a> command</li>
</ul>
</div>
<divclass="section"id="examples">
<h2>Examples<aclass="headerlink"href="#examples"title="Permalink to this headline">¶</a></h2>
<h2>Description<aclass="headerlink"href="#description"title="Permalink to this headline">¶</a></h2>
<p>Define a computation that calculates the angular velocity (omega) of
multiple chunks of atoms.</p>
<p>In LAMMPS, chunks are collections of atoms defined by a <aclass="reference internal"href="compute_chunk_atom.html"><em>compute chunk/atom</em></a> command, which assigns each atom
to a single chunk (or no chunk). The ID for this command is specified
as chunkID. For example, a single chunk could be the atoms in a
molecule or atoms in a spatial bin. See the <aclass="reference internal"href="compute_chunk_atom.html"><em>compute chunk/atom</em></a> doc page and “<aclass="reference internal"href="Section_howto.html#howto-23"><span>Section_howto 23</span></a> for details of how chunks can be
defined and examples of how they can be used to measure properties of
a system.</p>
<p>This compute calculates the 3 components of the angular velocity
vector for each chunk, via the formula L = Iw where L is the angular
momentum vector of the chunk, I is its moment of inertia tensor, and w
is omega = angular velocity of the chunk. The calculation includes
all effects due to atoms passing thru periodic boundaries.</p>
<p>Note that only atoms in the specified group contribute to the
calculation. The <aclass="reference internal"href="compute_chunk_atom.html"><em>compute chunk/atom</em></a> command
defines its own group; atoms will have a chunk ID = 0 if they are not
in that group, signifying they are not assigned to a chunk, and will
thus also not contribute to this calculation. You can specify the
“all” group for this command if you simply want to include atoms with
non-zero chunk IDs.</p>
<divclass="admonition warning">
<pclass="first admonition-title">Warning</p>
<pclass="last">The coordinates of an atom contribute to the chunk’s
angular velocity in “unwrapped” form, by using the image flags
associated with each atom. See the <aclass="reference internal"href="dump.html"><em>dump custom</em></a> command
for a discussion of “unwrapped” coordinates. See the Atoms section of
the <aclass="reference internal"href="read_data.html"><em>read_data</em></a> command for a discussion of image flags
and how they are set for each atom. You can reset the image flags
(e.g. to 0) before invoking this compute by using the <aclass="reference internal"href="set.html"><em>set image</em></a> command.</p>
</div>
<p>The simplest way to output the results of the compute omega/chunk
calculation to a file is to use the <aclass="reference internal"href="fix_ave_time.html"><em>fix ave/time</em></a>
command, for example:</p>
<divclass="highlight-python"><divclass="highlight"><pre>compute cc1 all chunk/atom molecule
<p>This compute calculates a global array where the number of rows = the
number of chunks <em>Nchunk</em> as calculated by the specified <aclass="reference internal"href="compute_chunk_atom.html"><em>compute chunk/atom</em></a> command. The number of columns =
3 for the 3 xyz components of the angular velocity for each chunk.
These values can be accessed by any command that uses global array
values from a compute as input. See <aclass="reference internal"href="Section_howto.html#howto-15"><span>Section_howto 15</span></a> for an overview of LAMMPS output
options.</p>
<p>The array values are “intensive”. The array values will be in
Built with <ahref="http://sphinx-doc.org/">Sphinx</a> using a <ahref="https://github.com/snide/sphinx_rtd_theme">theme</a> provided by <ahref="https://readthedocs.org">Read the Docs</a>.