git-svn-id: svn://svn.icms.temple.edu/lammps-ro/trunk@16053 f3b2605a-c512-4ea7-a41b...
[lammps.git] / doc / src / pair_mie.txt
blobfec00ed44c32bf856d8769d4eba75a22598920f3
1 "LAMMPS WWW Site"_lws - "LAMMPS Documentation"_ld - "LAMMPS Commands"_lc :c
3 :link(lws,http://lammps.sandia.gov)
4 :link(ld,Manual.html)
5 :link(lc,Section_commands.html#comm)
7 :line
9 pair_style mie/cut command :h3
10 pair_style mie/cut/gpu command :h3
12 [Syntax:]
14 pair_style mie/cut cutoff :pre
16 cutoff = global cutoff for mie/cut interactions (distance units) :ul
18 [Examples:]
20 pair_style mie/cut 10.0
21 pair_coeff 1 1 0.72 3.40 23.00 6.66
22 pair_coeff 2 2 0.30 3.55 12.65 6.00
23 pair_coeff 1 2 0.46 3.32 16.90 6.31 :pre
25 [Description:]
27 The {mie/cut} style computes the Mie potential, given by
29 :c,image(Eqs/pair_mie.jpg)
31 Rc is the cutoff and C is a function that depends on the repulsive and
32 attractive exponents, given by:
34 :c,image(Eqs/pair_mie2.jpg)
36 Note that for 12/6 exponents, C is equal to 4 and the formula is the
37 same as the standard Lennard-Jones potential.
39 The following coefficients must be defined for each pair of atoms
40 types via the "pair_coeff"_pair_coeff.html command as in the examples
41 above, or in the data file or restart files read by the
42 "read_data"_read_data.html or "read_restart"_read_restart.html
43 commands, or by mixing as described below:
45 epsilon (energy units)
46 sigma (distance units)
47 gammaR
48 gammaA
49 cutoff (distance units) :ul
51 The last coefficient is optional.  If not specified, the global
52 cutoff specified in the pair_style command is used.
54 :line
56 [Mixing, shift, table, tail correction, restart, rRESPA info]:
58 For atom type pairs I,J and I != J, the epsilon and sigma coefficients
59 and cutoff distance for all of the mie/cut pair styles can be mixed.
60 If not explicity defined, both the repulsive and attractive gamma
61 exponents for different atoms will be calculated following the same
62 mixing rule defined for distances.  The default mix value is
63 {geometric}. See the "pair_modify" command for details.
65 This pair style supports the "pair_modify"_pair_modify.html shift
66 option for the energy of the pair interaction.
68 This pair style supports the "pair_modify"_pair_modify.html tail
69 option for adding a long-range tail correction to the energy and
70 pressure of the pair interaction.
72 This pair style writes its information to "binary restart
73 files"_restart.html, so pair_style and pair_coeff commands do not need
74 to be specified in an input script that reads a restart file.
76 This pair style supports the use of the {inner}, {middle}, and {outer}
77 keywords of the "run_style respa"_run_style.html command, meaning the
78 pairwise forces can be partitioned by distance at different levels of
79 the rRESPA hierarchy.  See the "run_style"_run_style.html command for
80 details.
82 :line
84 [Restrictions:] none
86 [Related commands:]
88 "pair_coeff"_pair_coeff.html
90 [Default:] none
92 :line
94 :link(Mie)
95 [(Mie)] G. Mie, Ann Phys, 316, 657 (1903).
97 :link(Avendano)
98 [(Avendano)] C. Avendano, T. Lafitte, A. Galindo, C. S. Adjiman,
99 G. Jackson, E. Muller, J Phys Chem B, 115, 11154 (2011).