1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
|
"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 zbl command :h3
pair_style zbl/omp command :h3
[Syntax:]
pair_style zbl inner outer :pre
inner = distance where switching function begins
outer = global cutoff for ZBL interaction :ul
[Examples:]
pair_style zbl 3.0 4.0
pair_coeff * * 73.0
pair_coeff 1 1 14.0 :pre
[Description:]
Style {zbl} computes the Ziegler-Biersack-Littmark (ZBL) screened nuclear
repulsion for describing high-energy collisions between atoms.
"(Ziegler)"_#Ziegler. It includes an additional switching function
that ramps the energy, force, and curvature smoothly to zero
between an inner and outer cutoff. The potential
energy due to a pair of atoms at a distance r_ij is given by:
:c,image(Eqs/pair_zbl.jpg)
where e is the electron charge, epsilon_0 is the electrical
permittivity of vacuum, and Z_i and Z_j are the nuclear charges of the
two atoms. The switching function S(r) is identical to that used by
"pair_style lj/gromacs"_pair_gromacs.html. Here, the inner and outer
cutoff are the same for all pairs of atom types.
The following coefficient must be defined for each pair of atom types
via the "pair_coeff"_pair_coeff.html command as in the examples above,
or in the LAMMPS data file. Z can not be specified for two different
atoms types. Therefore the lists of atom types I and atom types J
must match.
Z (multiples of proton charge, e.g. 13.0 for aluminum) :ul
Although Z must be defined for all atom type pairs I,J, it is only
stored for individual atom types, i.e. when I = J. Z is normally equal
to the atomic number of the atom type.
IMPORTANT NOTE: The numerical values of the exponential decay
constants in the screening function depend on the unit of distance. In
the above equation they are given for units of angstroms. LAMMPS will
automatically convert these values to the distance unit of the
specified LAMMPS "units"_units.html setting. The values of Z should
always be given as multiples of a proton's charge, e.g. 29.0 for
copper.
:line
Styles with a {cuda}, {gpu}, {omp}, or {opt} suffix are functionally
the same as the corresponding style without the suffix. They have
been optimized to run faster, depending on your available hardware, as
discussed in "Section_accelerate"_Section_accelerate.html of the
manual. The accelerated styles take the same arguments and should
produce the same results, except for round-off and precision issues.
These accelerated styles are part of the USER-CUDA, GPU, USER-OMP and OPT
packages, respectively. They are only enabled if LAMMPS was built with
those packages. See the "Making LAMMPS"_Section_start.html#start_3
section for more info.
You can specify the accelerated styles explicitly in your input script
by including their suffix, or you can use the "-suffix command-line
switch"_Section_start.html#start_7 when you invoke LAMMPS, or you can
use the "suffix"_suffix.html command in your input script.
See "Section_accelerate"_Section_accelerate.html of the manual for
more instructions on how to use the accelerated styles effectively.
:line
[Mixing, shift, table, tail correction, restart, rRESPA info]:
Mixing is not relevant for this pair style, since as explained above,
Z values are stored on a per-type basis, and both Zi and Zj are used
explicitly in the ZBL formula.
The ZBL pair style does not support the "pair_modify"_pair_modify.html
shift option, since the ZBL interaction is already smoothed to 0.0 at
the cutoff.
The "pair_modify"_pair_modify.html table option is not relevant for
this pair style.
This pair style does not support the "pair_modify"_pair_modify.html
tail option for adding long-range tail corrections to energy and
pressure, since there are no corrections for a potential that goes to
0.0 at the cutoff.
This pair style does not write information to "binary restart
files"_restart.html, so pair_style and pair_coeff commands must be
specified in an input script that reads a restart file.
This pair style can only be used via the {pair} keyword of the
"run_style respa"_run_style.html command. It does not support the
{inner}, {middle}, {outer} keywords.
:line
[Restrictions:] none
[Related commands:]
"pair_coeff"_pair_coeff.html
[Default:] none
:line
:link(Ziegler)
[(Ziegler)] J.F. Ziegler, J. P. Biersack and U. Littmark, "The
Stopping and Range of Ions in Matter," Volume 1, Pergamon, 1985.
|
