jianmu
/
lammps
forked from lijiext/lammps
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

git-svn-id: svn://svn.icms.temple.edu/lammps-ro/trunk@13733 f3b2605a-c512-4ea7-a41b-209d697bcdaa

This commit is contained in:
sjplimp 2015-07-27 18:05:49 +00:00
parent 059f719d62
commit f524105378
2 changed files with 259 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

133
doc/pair_thole.html Normal file
View File

@ -0,0 +1,133 @@
<HTML>
<script type="text/javascript"
src="https://cdn.mathjax.org/mathjax/latest/MathJax.js?config=TeX-AMS-MML_HTMLorMML">
</script>
<script type="text/x-mathjax-config">
MathJax.Hub.Config({ TeX: { equationNumbers: {autoNumber: "AMS"} } });
</script>
<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>
</CENTER>
<HR>
<H3>pair_style thole command
</H3>
<P><B>Syntax:</B>
</P>
<PRE>pair_style thole damp cutoff
</PRE>
<UL><LI>thole = style name
<LI>damp = global damping parameter
<LI>cutoff = global cutoff
</UL>
<P><B>Examples:</B>
</P>
<PRE>pair_style hybrid/overlay ... thole 2.6 12.0
pair_coeff 1 1 thole 1.0
pair_coeff 1 2 thole 1.0 2.6 10.0
pair_coeff * 2 thole 1.0 2.6
</PRE>
<P><B>Description:</B>
</P>
<P>The <I>thole</I> pair style is meant to be used with force fields that
include explicit polarization through Drude dipoles. This link
describes how to use the <A HREF = "tutorial_drude.html">thermalized Drude oscillator
model</A> in LAMMPS and polarizable models in LAMMPS
are discussed in <A HREF = "Section_howto.html#howto_25">this Section</A>.
</P>
<P>The <I>thole</I> pair style should be used as a sub-style within in the
<A HREF = "pair_hybrid.html">pair_hybrid/overlay</A> command, in conjunction with a
main pair style including Coulomb interactions, i.e. any pair style
containing <I>coul/cut</I> or <I>coul/long</I> in its style name.
</P>
<P>The <I>thole</I> pair style computes the Coulomb interaction damped at
short distances by a function
</P>
<P>\begin{equation} T_{ij}(r_{ij}) = 1 - \left( 1 +
\frac{s_{ij} r_{ij} }{2} \right)
\exp \left( - s_{ij} r_{ij} \right) \end{equation}
</P>
<P>This function results from an adaptation to point charges
<A HREF = "#Noskov">(Noskov)</A> of the dipole screening scheme originally proposed
by <A HREF = "#Thole">Thole</A>. The scaling coefficient \(s_{ij} \) is determined
by the polarizability of the atoms, \( \alpha_i \), and by a Thole
damping parameter \( a \). This Thole damping parameter usually takes
a value of 2.6, but in certain force fields the value can depend upon
the atom types. The mixing rule for Thole damping parameters is the
arithmetic average, and for polarizabilities the geometric average
between the atom-specific values.
</P>
<P>\begin{equation} s_{ij} = \frac{ a_{ij} }{
(\alpha_{ij})^{1/3} } = \frac{ (a_i + a_j)/2 }{
[(\alpha_i\alpha_j)^{1/2}]^{1/3} } \end{equation}
</P>
<P>The damping function is only applied to the interactions between the
point charges representing the induced dipoles on polarizable sites,
that is, charges on Drude particles, \( q_{D,i} \), and opposite
charges, \( -q_{D,i} \), located on the respective core particles
(to which each Drude particle is bonded). Therefore, Thole screening
is not applied to the full charge of the core particle \( q_i \), but
only to the \( -q_{D,i} \) part of it.
</P>
<P>The interactions between core charges are subject to the weighting
factors set by the <A HREF = "special_bonds.html">special_bonds</A> command. The
interactions between Drude particles and core charges or
non-polarizable atoms are also subject to these weighting factors. The
Drude particles inherit the 1-2, 1-3 and 1-4 neighbor relations from
their respective cores.
</P>
<P>The following coefficients must be defined for each pair of atoms
types via the <A HREF = "pair_coeff.html">pair_coeff</A> command as in the example
above.
</P>
<UL><LI>alpha (distance units^3)
<LI>damp
<LI>cutoff (distance units)
</UL>
<P>The last two coefficients are optional. If not specified the global
Thole damping parameter or global cutoff specified in the pair_style
command are used. In order to specify a cutoff (third argument) a damp
parameter (second argument) must also be specified.
</P>
<P><B>Mixing</B>:
</P>
<P>The <I>thole</I> pair style does not support mixing. Thus, coefficients
for all I,J pairs must be specified explicitly.
</P>
<P><B>Restrictions:</B>
</P>
<P>These pair styles are part of the USER-DRUDE package. They are only
enabled if LAMMPS was built with that package. See the <A HREF = "Section_start.html#start_3">Making
LAMMPS</A> section for more info.
</P>
<P>This pair_style should currently not be used with the <A HREF = "dihedral_charmm.html">charmm dihedral
style</A> if the latter non-zero 1-4 weighting
factors. This is because the <I>thole</I> pair style does not know which
pairs are 1-4 partners of which dihedrals.
</P>
<P><B>Related commands:</B>
</P>
<P><A HREF = "fix_drude.html">fix drude</A>, <A HREF = "fix_langevin_drude.html">fix
langevin/drude</A>, <A HREF = "fix_drude_transform.html">fix
drude/transform</A>, <A HREF = "compute_temp_drude.html">compute
temp/drude</A>
</P>
<P><B>Default:</B> none
</P>
<HR>
<A NAME = "Noskov"></A>
<P><B>(Noskov)</B> Noskov, Lamoureux and Roux, J Phys Chem B, 109, 6705 (2005).
</P>
<A NAME = "Thole"></A>
<P><B>(Thole)</B> Chem Phys, 59, 341 (1981).
</P>
</HTML>

126
doc/pair_thole.txt Normal file
View File

@ -0,0 +1,126 @@
<script type="text/javascript"
src="https://cdn.mathjax.org/mathjax/latest/MathJax.js?config=TeX-AMS-MML_HTMLorMML">
</script>
<script type="text/x-mathjax-config">
MathJax.Hub.Config({ TeX: { equationNumbers: {autoNumber: "AMS"} } });
</script>
"LAMMPS WWW Site"_lws - "LAMMPS Documentation"_ld - "LAMMPS Commands"_lc :c
:link(lws,http://lammps.sandia.gov)
:link(ld,Manual.html)
:link(lc,Section_commands.html#comm)
:line
pair_style thole command :h3
[Syntax:]
pair_style thole damp cutoff :pre
thole = style name
damp = global damping parameter
cutoff = global cutoff :ul
[Examples:]
pair_style hybrid/overlay ... thole 2.6 12.0
pair_coeff 1 1 thole 1.0
pair_coeff 1 2 thole 1.0 2.6 10.0
pair_coeff * 2 thole 1.0 2.6 :pre
[Description:]
The {thole} pair style is meant to be used with force fields that
include explicit polarization through Drude dipoles. This link
describes how to use the "thermalized Drude oscillator
model"_tutorial_drude.html in LAMMPS and polarizable models in LAMMPS
are discussed in "this Section"_Section_howto.html#howto_25.
The {thole} pair style should be used as a sub-style within in the
"pair_hybrid/overlay"_pair_hybrid.html command, in conjunction with a
main pair style including Coulomb interactions, i.e. any pair style
containing {coul/cut} or {coul/long} in its style name.
The {thole} pair style computes the Coulomb interaction damped at
short distances by a function
\begin\{equation\} T_\{ij\}(r_\{ij\}) = 1 - \left( 1 +
\frac\{s_\{ij\} r_\{ij\} \}\{2\} \right)
\exp \left( - s_\{ij\} r_\{ij\} \right) \end\{equation\}
This function results from an adaptation to point charges
"(Noskov)"_#Noskov of the dipole screening scheme originally proposed
by "Thole"_#Thole. The scaling coefficient \(s_\{ij\} \) is determined
by the polarizability of the atoms, \( \alpha_i \), and by a Thole
damping parameter \( a \). This Thole damping parameter usually takes
a value of 2.6, but in certain force fields the value can depend upon
the atom types. The mixing rule for Thole damping parameters is the
arithmetic average, and for polarizabilities the geometric average
between the atom-specific values.
\begin\{equation\} s_\{ij\} = \frac\{ a_\{ij\} \}\{
(\alpha_\{ij\})^\{1/3\} \} = \frac\{ (a_i + a_j)/2 \}\{
\[(\alpha_i\alpha_j)^\{1/2\}\]^\{1/3\} \} \end\{equation\}
The damping function is only applied to the interactions between the
point charges representing the induced dipoles on polarizable sites,
that is, charges on Drude particles, \( q_\{D,i\} \), and opposite
charges, \( -q_\{D,i\} \), located on the respective core particles
(to which each Drude particle is bonded). Therefore, Thole screening
is not applied to the full charge of the core particle \( q_i \), but
only to the \( -q_\{D,i\} \) part of it.
The interactions between core charges are subject to the weighting
factors set by the "special_bonds"_special_bonds.html command. The
interactions between Drude particles and core charges or
non-polarizable atoms are also subject to these weighting factors. The
Drude particles inherit the 1-2, 1-3 and 1-4 neighbor relations from
their respective cores.
The following coefficients must be defined for each pair of atoms
types via the "pair_coeff"_pair_coeff.html command as in the example
above.
alpha (distance units^3)
damp
cutoff (distance units) :ul
The last two coefficients are optional. If not specified the global
Thole damping parameter or global cutoff specified in the pair_style
command are used. In order to specify a cutoff (third argument) a damp
parameter (second argument) must also be specified.
[Mixing]:
The {thole} pair style does not support mixing. Thus, coefficients
for all I,J pairs must be specified explicitly.
[Restrictions:]
These pair styles are part of the USER-DRUDE package. They are only
enabled if LAMMPS was built with that package. See the "Making
LAMMPS"_Section_start.html#start_3 section for more info.
This pair_style should currently not be used with the "charmm dihedral
style"_dihedral_charmm.html if the latter non-zero 1-4 weighting
factors. This is because the {thole} pair style does not know which
pairs are 1-4 partners of which dihedrals.
[Related commands:]
"fix drude"_fix_drude.html, "fix
langevin/drude"_fix_langevin_drude.html, "fix
drude/transform"_fix_drude_transform.html, "compute
temp/drude"_compute_temp_drude.html
[Default:] none
:line
:link(Noskov)
[(Noskov)] Noskov, Lamoureux and Roux, J Phys Chem B, 109, 6705 (2005).
:link(Thole)
[(Thole)] Chem Phys, 59, 341 (1981).