Page Menu
Home
c4science
Search
Configure Global Search
Log In
Files
F81222430
dihedral_spherical.html
No One
Temporary
Actions
Download File
Edit File
Delete File
View Transforms
Subscribe
Mute Notifications
Award Token
Subscribers
None
File Metadata
Details
File Info
Storage
Attached
Created
Wed, Sep 4, 20:22
Size
10 KB
Mime Type
text/html
Expires
Fri, Sep 6, 20:22 (2 d)
Engine
blob
Format
Raw Data
Handle
20526748
Attached To
rLAMMPS lammps
dihedral_spherical.html
View Options
<!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>
dihedral_style spherical command
—
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=
"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"
>
<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
&
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
&
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>
</div>
</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"
>
<ul
class=
"wy-breadcrumbs"
>
<li><a
href=
"Manual.html"
>
Docs
</a>
»
</li>
<li>
dihedral_style spherical 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=
"dihedral-style-spherical-command"
>
<span
id=
"index-0"
></span><h1>
dihedral_style spherical command
</h1>
<div
class=
"section"
id=
"syntax"
>
<h2>
Syntax
</h2>
<div
class=
"highlight-default"
><div
class=
"highlight"
><pre><span></span><span
class=
"n"
>
dihedral_style
</span>
<span
class=
"n"
>
spherical
</span>
</pre></div>
</div>
</div>
<div
class=
"section"
id=
"examples"
>
<h2>
Examples
</h2>
<div
class=
"highlight-default"
><div
class=
"highlight"
><pre><span></span><span
class=
"n"
>
dihedral_coeff
</span>
<span
class=
"mi"
>
1
</span>
<span
class=
"mi"
>
1
</span>
<span
class=
"mf"
>
286.1
</span>
<span
class=
"mi"
>
1
</span>
<span
class=
"mi"
>
124
</span>
<span
class=
"mi"
>
1
</span>
<span
class=
"mi"
>
1
</span>
<span
class=
"mf"
>
90.0
</span>
<span
class=
"mi"
>
0
</span>
<span
class=
"mi"
>
1
</span>
<span
class=
"mf"
>
90.0
</span>
<span
class=
"mi"
>
0
</span>
<span
class=
"n"
>
dihedral_coeff
</span>
<span
class=
"mi"
>
1
</span>
<span
class=
"mi"
>
3
</span>
<span
class=
"mf"
>
286.1
</span>
<span
class=
"mi"
>
1
</span>
<span
class=
"mi"
>
114
</span>
<span
class=
"mi"
>
1
</span>
<span
class=
"mi"
>
1
</span>
<span
class=
"mi"
>
90
</span>
<span
class=
"mi"
>
0
</span>
<span
class=
"mi"
>
1
</span>
<span
class=
"mf"
>
90.0
</span>
<span
class=
"mi"
>
0
</span>
<span
class=
"o"
>
&
</span>
<span
class=
"mf"
>
17.3
</span>
<span
class=
"mi"
>
0
</span>
<span
class=
"mf"
>
0.0
</span>
<span
class=
"mi"
>
0
</span>
<span
class=
"mi"
>
1
</span>
<span
class=
"mi"
>
158
</span>
<span
class=
"mi"
>
1
</span>
<span
class=
"mi"
>
0
</span>
<span
class=
"mf"
>
0.0
</span>
<span
class=
"mi"
>
0
</span>
<span
class=
"o"
>
&
</span>
<span
class=
"mf"
>
15.1
</span>
<span
class=
"mi"
>
0
</span>
<span
class=
"mf"
>
0.0
</span>
<span
class=
"mi"
>
0
</span>
<span
class=
"mi"
>
0
</span>
<span
class=
"mf"
>
0.0
</span>
<span
class=
"mi"
>
0
</span>
<span
class=
"mi"
>
1
</span>
<span
class=
"mf"
>
167.3
</span>
<span
class=
"mi"
>
1
</span>
</pre></div>
</div>
</div>
<div
class=
"section"
id=
"description"
>
<h2>
Description
</h2>
<p>
The
<em>
spherical
</em>
dihedral style uses the potential:
</p>
<img
alt=
"_images/dihedral_spherical_angles.jpg"
class=
"align-center"
src=
"_images/dihedral_spherical_angles.jpg"
/>
<img
alt=
"_images/dihedral_spherical.jpg"
class=
"align-center"
src=
"_images/dihedral_spherical.jpg"
/>
<p>
For this dihedral style, the energy can be any function that combines the
4-body dihedral-angle (phi) and the two 3-body bond-angles (theta1, theta2).
For this reason, there is usually no need to define 3-body
“
angle
”
forces
separately for the atoms participating in these interactions.
It is probably more efficient to incorporate 3-body angle forces into
the dihedral interaction even if it requires adding additional terms to
the expansion (as was done in the second example). A careful choice of
parameters can prevent singularities that occur with traditional
force-fields whenever theta1 or theta2 approach 0 or 180 degrees.
The last example above corresponds to an interaction with a single energy
minima located at phi=114, theta1=158, theta2=167.3 degrees, and it remains
numerically stable at all angles (phi, theta1, theta2). In this example,
the coefficients 17.3, and 15.1 can be physically interpreted as the
harmonic spring constants for theta1 and theta2 around their minima.
The coefficient 286.1 is the harmonic spring constant for phi after
division by sin(158)*sin(167.3) (the minima positions for theta1 and theta2).
</p>
<p>
The following coefficients must be defined for each dihedral type via the
<a
class=
"reference internal"
href=
"dihedral_coeff.html"
><span
class=
"doc"
>
dihedral_coeff
</span></a>
command as in the example above, or in
the Dihedral Coeffs section of a data file file read by the
<a
class=
"reference internal"
href=
"read_data.html"
><span
class=
"doc"
>
read_data
</span></a>
command:
</p>
<ul
class=
"simple"
>
<li>
n (integer
>
= 1)
</li>
<li>
C1 (energy)
</li>
<li>
K1 (typically an integer)
</li>
<li>
a1 (degrees)
</li>
<li>
u1 (typically 0.0 or 1.0)
</li>
<li>
L1 (typically an integer)
</li>
<li>
b1 (degrees, typically 0.0 or 90.0)
</li>
<li>
v1 (typically 0.0 or 1.0)
</li>
<li>
M1 (typically an integer)
</li>
<li>
c1 (degrees, typically 0.0 or 90.0)
</li>
<li>
w1 (typically 0.0 or 1.0)
</li>
<li></li>
<li>
Cn (energy)
</li>
<li>
Kn (typically an integer)
</li>
<li>
an (degrees)
</li>
<li>
un (typically 0.0 or 1.0)
</li>
<li>
Ln (typically an integer)
</li>
<li>
bn (degrees, typically 0.0 or 90.0)
</li>
<li>
vn (typically 0.0 or 1.0)
</li>
<li>
Mn (typically an integer)
</li>
<li>
cn (degrees, typically 0.0 or 90.0)
</li>
<li>
wn (typically 0.0 or 1.0)
</li>
</ul>
</div>
<hr
class=
"docutils"
/>
<div
class=
"section"
id=
"restrictions"
>
<h2>
Restrictions
</h2>
<p>
This dihedral style can only be used if LAMMPS was built with the
USER_MISC 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 on packages.
</p>
</div>
<div
class=
"section"
id=
"related-commands"
>
<h2>
Related commands
</h2>
<p><a
class=
"reference internal"
href=
"dihedral_coeff.html"
><span
class=
"doc"
>
dihedral_coeff
</span></a></p>
<p><strong>
Default:
</strong>
none
</p>
</div>
</div>
</div>
</div>
<footer>
<hr/>
<div
role=
"contentinfo"
>
<p>
©
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://cdn.mathjax.org/mathjax/latest/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>
Event Timeline
Log In to Comment