Skip to content
Snippets Groups Projects
Unverified Commit efbee4a2 authored by Michael DeLyser's avatar Michael DeLyser Committed by GitHub
Browse files

Delete man_fix_bocs.html

parent 2e2f1977
No related branches found
No related tags found
No related merge requests found
<!DOCTYPE html>
<!--[if IE 8]><html class="no-js lt-ie9" lang="en" > <![endif]-->
<!--[if gt IE 8]><!--> <html class="no-js" lang="en" > <!--<![endif]-->
<head>
<meta charset="utf-8">
<meta name="viewport" content="width=device-width, initial-scale=1.0">
<title>fix bocs command &mdash; LAMMPS documentation</title>
<link rel="stylesheet" href="_static/css/theme.css" type="text/css" />
<link rel="stylesheet" href="_static/sphinxcontrib-images/LightBox2/lightbox2/css/lightbox.css" type="text/css" />
<link rel="index" title="Index"
href="genindex.html"/>
<link rel="search" title="Search" href="search.html"/>
<link rel="top" title="LAMMPS documentation" href="index.html"/>
<script src="_static/js/modernizr.min.js"></script>
</head>
<body class="wy-body-for-nav" role="document">
<div class="wy-grid-for-nav">
<nav data-toggle="wy-nav-shift" class="wy-nav-side">
<div class="wy-side-nav-search">
<a href="Manual.html" class="icon icon-home"> LAMMPS
</a>
<div role="search">
<form id="rtd-search-form" class="wy-form" action="search.html" method="get">
<input type="text" name="q" placeholder="Search docs" />
<input type="hidden" name="check_keywords" value="yes" />
<input type="hidden" name="area" value="default" />
</form>
</div>
</div>
<div class="wy-menu wy-menu-vertical" data-spy="affix" role="navigation" aria-label="main navigation">
<p class="caption"><span class="caption-text">User Documentation</span></p>
<ul>
<li class="toctree-l1"><a class="reference internal" href="Section_intro.html">1. Introduction</a></li>
<li class="toctree-l1"><a class="reference internal" href="Section_start.html">2. Getting Started</a></li>
<li class="toctree-l1"><a class="reference internal" href="Section_commands.html">3. Commands</a></li>
<li class="toctree-l1"><a class="reference internal" href="Section_packages.html">4. Packages</a></li>
<li class="toctree-l1"><a class="reference internal" href="Section_accelerate.html">5. Accelerating LAMMPS performance</a></li>
<li class="toctree-l1"><a class="reference internal" href="Section_howto.html">6. How-to discussions</a></li>
<li class="toctree-l1"><a class="reference internal" href="Section_example.html">7. Example problems</a></li>
<li class="toctree-l1"><a class="reference internal" href="Section_perf.html">8. Performance &amp; scalability</a></li>
<li class="toctree-l1"><a class="reference internal" href="Section_tools.html">9. Additional tools</a></li>
<li class="toctree-l1"><a class="reference internal" href="Section_modify.html">10. Modifying &amp; extending LAMMPS</a></li>
<li class="toctree-l1"><a class="reference internal" href="Section_python.html">11. Python interface to LAMMPS</a></li>
<li class="toctree-l1"><a class="reference internal" href="Section_errors.html">12. Errors</a></li>
<li class="toctree-l1"><a class="reference internal" href="Section_history.html">13. Future and history</a></li>
</ul>
<p class="caption"><span class="caption-text">Index</span></p>
<ul>
<li class="toctree-l1"><a class="reference internal" href="tutorials.html">Tutorials</a></li>
<li class="toctree-l1"><a class="reference internal" href="commands.html">Commands</a></li>
<li class="toctree-l1"><a class="reference internal" href="fixes.html">Fixes</a></li>
<li class="toctree-l1"><a class="reference internal" href="computes.html">Computes</a></li>
<li class="toctree-l1"><a class="reference internal" href="pairs.html">Pair Styles</a></li>
<li class="toctree-l1"><a class="reference internal" href="bonds.html">Bond Styles</a></li>
<li class="toctree-l1"><a class="reference internal" href="angles.html">Angle Styles</a></li>
<li class="toctree-l1"><a class="reference internal" href="dihedrals.html">Dihedral Styles</a></li>
<li class="toctree-l1"><a class="reference internal" href="impropers.html">Improper Styles</a></li>
</ul>
</div>
&nbsp;
</nav>
<section data-toggle="wy-nav-shift" class="wy-nav-content-wrap">
<nav class="wy-nav-top" role="navigation" aria-label="top navigation">
<i data-toggle="wy-nav-top" class="fa fa-bars"></i>
<a href="Manual.html">LAMMPS</a>
</nav>
<div class="wy-nav-content">
<div class="rst-content">
<div role="navigation" aria-label="breadcrumbs navigation">
<div style="text-align: center; margin-bottom: -1.5em; display: block"><b>LAMMPS</b> 11 Aug 2017</div>
<ul class="wy-breadcrumbs">
<li><a href="Manual.html">Docs</a> &raquo;</li>
<li>fix bocs command</li>
<li class="wy-breadcrumbs-aside">
<a href="http://lammps.sandia.gov">Website</a>
<a href="Section_commands.html#comm">Commands</a>
</li>
</ul>
<hr/>
</div>
<div role="main" class="document" itemscope="itemscope" itemtype="http://schema.org/Article">
<div itemprop="articleBody">
<div class="section" id="fix-bocs-command">
<span id="index-0"></span><h1>fix bocs command</h1>
<div class="section" id="syntax">
<h2>Syntax</h2>
<pre class="literal-block">
fix ID group-ID bocs keyword values ...
keyword = <em>temp</em> or <em>cgiso</em> or <em>analytic</em> or <em>linear_spline</em> or <em>cubic_spline</em>
<em>temp</em> values = Tstart Tstop Tdamp
<em>cgiso</em> values = Pstart Pstop Pdamp
<em>basis set</em>
<em>analytic</em> values = V_avg N_particles N_coeff Coeff_1 Coeff_2 ... Coeff_N
<em>linear_spline</em> values = input_filename
<em>cubic_spline</em> values = input_filename
</pre>
</div>
<div class="section" id="examples">
<h2>Examples</h2>
<div class="highlight-default"><div class="highlight"><pre><span></span><span class="n">fix</span> <span class="mi">1</span> <span class="nb">all</span> <span class="n">bocs</span> <span class="n">temp</span> <span class="mf">300.0</span> <span class="mf">300.0</span> <span class="mf">100.0</span> <span class="n">cgiso</span> <span class="mf">0.986</span> <span class="mf">0.986</span> <span class="mf">1000.0</span> <span class="n">analytic</span> <span class="mf">66476.015</span> <span class="mi">968</span> <span class="mi">2</span> <span class="mf">245030.10</span> <span class="mf">8962.20</span>
<span class="n">fix</span> <span class="mi">1</span> <span class="nb">all</span> <span class="n">bocs</span> <span class="n">temp</span> <span class="mf">300.0</span> <span class="mf">300.0</span> <span class="mf">100.0</span> <span class="n">cgiso</span> <span class="mf">0.986</span> <span class="mf">0.986</span> <span class="mf">1000.0</span> <span class="n">cubic_spline</span> <span class="n">input_Fv</span><span class="o">.</span><span class="n">dat</span>
<span class="n">compute</span> <span class="n">bocsPress</span> <span class="nb">all</span> <span class="n">pressureBocs</span> <span class="n">thermo_temp</span>
<span class="n">thermo_modify</span> <span class="n">press</span> <span class="n">bocsPress</span>
<span class="n">fix_modify</span> <span class="mi">1</span> <span class="n">press</span> <span class="n">bocsPress</span>
</pre></div>
</div>
</div>
<div class="section" id="description">
<h2>Description</h2>
<p>These commands incorporate a pressure correction as described by
Dunn and Noid in <a class="reference internal" href="#bocs-dunn1"><span class="std std-ref">(Dunn1)</span></a> to the standard MTTK
barostat by Martyna et. al. in <a class="reference internal" href="#bocs-martyna"><span class="std std-ref">(Martyna)</span></a> .
The first half of the command mimics a standard fix npt command:</p>
<div class="highlight-default"><div class="highlight"><pre><span></span><span class="n">fix</span> <span class="mi">1</span> <span class="nb">all</span> <span class="n">bocs</span> <span class="n">temp</span> <span class="n">Tstart</span> <span class="n">Tstop</span> <span class="n">Tcoupl</span> <span class="n">cgiso</span> <span class="n">Pstart</span> <span class="n">Pstop</span> <span class="n">Pdamp</span>
</pre></div>
</div>
<p>The two differences are replacing <em>npt</em> with <em>bocs</em>, and replacing
<em>iso</em>/<em>aniso</em>/<em>etc</em> with <em>cgiso</em>.
The rest of the command details what form you would like to use for
the pressure correction equation. The choices are: <em>analytic</em>, <em>linear_spline</em>,
or <em>cubic_spline</em>.</p>
<p>With either spline method, the only argument that needs to follow it
is the name of a file that contains the desired pressure correction
as a function of volume. The file should be formatted so each line has:</p>
<div class="highlight-default"><div class="highlight"><pre><span></span><span class="n">Volume_i</span><span class="p">,</span> <span class="n">PressureCorrection_i</span>
</pre></div>
</div>
<p>Note both the COMMA and the SPACE separating the volume&#8217;s
value and its corresponding pressure correction. The volumes in the file
should be uniformly spaced. Both the volumes and the pressure corrections
should be provided in the proper units, e.g. if you are using <em>units real</em>,
the volumes should all be in cubic angstroms, and the pressure corrections
should all be in atomspheres. Furthermore, the table should start/end at a
volume considerably smaller/larger than you expect your system to sample
during the simulation. If the system ever reaches a volume outside of the
range provided, the simulation will stop.</p>
<p>With the <em>analytic</em> option, the arguments are as follows:</p>
<div class="highlight-default"><div class="highlight"><pre><span></span><span class="o">...</span> <span class="n">analytic</span> <span class="n">V_avg</span> <span class="n">N_particles</span> <span class="n">N_coeff</span> <span class="n">Coeff_1</span> <span class="n">Coeff_2</span> <span class="o">...</span> <span class="n">Coeff_N</span>
</pre></div>
</div>
<p>Note that <em>V_avg</em> and <em>Coeff_i</em> should all be in the proper units, e.g. if you
are using <em>units real</em>, <em>V_avg</em> should be in cubic angstroms, and the
coefficients should all be in atmospheres * cubic angstroms.</p>
</div>
<div class="section" id="restrictions">
<h2>Restrictions</h2>
<p>As this is computing a (modified) pressure, group-ID should be <em>all</em>.</p>
<p>The pressure correction has only been tested for use with an isotropic
pressure coupling in 3 dimensions.</p>
<p>There are three additional commands that must be supplied along with fix
bocs. They are given in the examples section. The name <em>bocsPress</em> can be
changed for any name of your choosing, provided it remains consistent across
all three commands. Additionally, the ID <em>1</em> in the <em>fix_modify</em> command must
match up with the <em>ID</em> specified in the original <em>fix bocs</em> command. Everything else
in those three commands must be unchanged. The first two commands can be
specified before the <em>fix bocs</em> command, but the <em>fix_modify</em> command must be
given after the original <em>fix bocs</em> command.</p>
<p>The <em>compute</em> command tells LAMMPS to compute the pressure using the modified
barostat.</p>
<p>The <em>thermo_modify</em> command tells LAMMPS to report the pressure from the modified
barostat instead of the default pressure, i.e. thermo_press.</p>
<p>The <em>fix_modify</em> command tells LAMMPS to use the pressure from the modified
barostat for the fix, instead of using the default presure, i.e. thermo_press.</p>
<p><strong>Related:</strong></p>
<p>For more details about the pressure correction and the entire BOCS software
package, visit the <a class="reference external" href="https://github.com/noid-group/BOCS">BOCS package on github</a> and read the release
paper by Dunn et. al. <a class="reference internal" href="#bocs-dunn2"><span class="std std-ref">(Dunn2)</span></a> .</p>
<hr class="docutils" />
<p id="bocs-dunn1"><strong>(Dunn1)</strong> Dunn and Noid, J Chem Phys, 143, 243148 (2015).</p>
<p id="bocs-martyna"><strong>(Martyna)</strong> Martyna, Tobias, and Klein, J Chem Phys, 101, 4177 (1994).</p>
<p id="bocs-dunn2"><strong>(Dunn2)</strong> Dunn, Lebold, DeLyser, Rudzinski, and Noid, J. Phys. Chem. B, 122, 3363 (2018).</p>
</div>
</div>
</div>
</div>
<footer>
<hr/>
<div role="contentinfo">
<p>
&copy; Copyright 2013 Sandia Corporation.
</p>
</div>
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>.
</footer>
</div>
</div>
</section>
</div>
<script type="text/javascript">
var DOCUMENTATION_OPTIONS = {
URL_ROOT:'./',
VERSION:'',
COLLAPSE_INDEX:false,
FILE_SUFFIX:'.html',
HAS_SOURCE: true
};
</script>
<script type="text/javascript" src="_static/jquery.js"></script>
<script type="text/javascript" src="_static/underscore.js"></script>
<script type="text/javascript" src="_static/doctools.js"></script>
<script type="text/javascript" src="https://cdnjs.cloudflare.com/ajax/libs/mathjax/2.7.0/MathJax.js?config=TeX-AMS-MML_HTMLorMML"></script>
<script type="text/javascript" src="_static/sphinxcontrib-images/LightBox2/lightbox2/js/jquery-1.11.0.min.js"></script>
<script type="text/javascript" src="_static/sphinxcontrib-images/LightBox2/lightbox2/js/lightbox.min.js"></script>
<script type="text/javascript" src="_static/sphinxcontrib-images/LightBox2/lightbox2-customize/jquery-noconflict.js"></script>
<script type="text/javascript" src="_static/js/theme.js"></script>
<script type="text/javascript">
jQuery(function () {
SphinxRtdTheme.StickyNav.enable();
});
</script>
</body>
</html>
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment