doc/src/pair_mdf.txt

   1 "LAMMPS WWW Site"_lws - "LAMMPS Documentation"_ld - "LAMMPS Commands"_lc :c
   2
   3 :link(lws,http://lammps.sandia.gov)
   4 :link(ld,Manual.html)
   5 :link(lc,Section_commands.html#comm)
   6
   7 :line
   8
   9 pair_style lj/mdf command :h3
  10 pair_style buck/mdf command :h3
  11 pair_style lennard/mdf command :h3
  12
  13 [Syntax:]
  14
  15 pair_style style args :pre
  16
  17 style = {lj/mdf} or {buck/mdf} or {lennard/mdf} :ulb,l
  18 args = list of arguments for a particular style :l
  19   {lj/mdf} args = cutoff1 cutoff2
  20     cutoff1 = inner cutoff for the start of the tapering function
  21     cutoff1 = out cutoff for the end of the tapering function
  22   {buck/mdf} args = cutoff1 cutoff2
  23     cutoff1 = inner cutoff for the start of the tapering function
  24     cutoff1 = out cutoff for the end of the tapering function
  25   {lennard/mdf} args = cutoff1 cutoff2
  26     cutoff1 = inner cutoff for the start of the tapering function
  27     cutoff1 = out cutoff for the end of the tapering function :pre
  28 :ule
  29
  30 [Examples:]
  31
  32 pair_style lj/mdf 2.5 3.0
  33 pair_coeff * * 1 1
  34 pair_coeff 1 1 1 1.1 2.8 3.0 3.2 :pre
  35
  36 pair_style buck 2.5 3.0
  37 pair_coeff * * 100.0 1.5 200.0
  38 pair_coeff * * 100.0 1.5 200.0 3.0 3.5 :pre
  39
  40 pair_style lennard/mdf 2.5 3.0
  41 pair_coeff * * 1 1
  42 pair_coeff 1 1 1 1.1 2.8 3.0 3.2 :pre
  43
  44 [Description:]
  45
  46 The {lj/mdf}, {buck/mdf} and {lennard/mdf} compute the standard 12-6
  47 Lennard-Jones and Buckingham potential with the addition of a taper
  48 function that ramps the energy and force smoothly to zero between an
  49 inner and outer cutoff.
  50
  51 :c,image(Eqs/pair_mdf-1.jpg)
  52
  53 The tapering, {f(r)}, is done by using the Mei, Davenport, Fernando
  54 function "(Mei)"_#Mei.
  55
  56 :c,image(Eqs/pair_mdf-2.jpg)
  57
  58 where
  59
  60 :c,image(Eqs/pair_mdf-3.jpg)
  61
  62 Here {r_m} is the inner cutoff radius and {r_cut} is the outer cutoff
  63 radius.
  64
  65 :line
  66
  67 For the {lj/mdf} pair_style, the potential energy, {E(r)}, is the
  68 standard 12-6 Lennard-Jones written in the epsilon/sigma form:
  69
  70 :c,image(Eqs/pair_mdf-4.jpg)
  71
  72 The following coefficients must be defined for each pair of atoms
  73 types via the pair_coeff command as in the examples above, or in the
  74 data file or restart files read by the "read_data"_read_data.html or
  75 "read_restart commands"_read_restart.html, or by mixing as described
  76 below:
  77
  78 epsilon (energy units)
  79 sigma (distance units)
  80 r_m (distance units)
  81 r_{cut} (distance units) :ul
  82
  83 :line
  84
  85 For the {buck/mdf} pair_style, the potential energy, {E(r)}, is the
  86 standard Buckingham potential:
  87
  88 :c,image(Eqs/pair_mdf-5.jpg)
  89
  90 A (energy units)
  91 \rho (distance units)
  92 C (energy-distance^6 units)
  93 r_m (distance units)
  94 r_{cut}$ (distance units) :ul
  95
  96 :line
  97
  98 For the {lennard/mdf} pair_style, the potential energy, {E(r)}, is the
  99 standard 12-6 Lennard-Jones written in the $A/B$ form:
 100
 101 :c,image(Eqs/pair_mdf-6.jpg)
 102
 103 The following coefficients must be defined for each pair of atoms
 104 types via the pair_coeff command as in the examples above, or in the
 105 data file or restart files read by the read_data or read_restart
 106 commands, or by mixing as described below:
 107
 108 A (energy-distance^12 units)
 109 B (energy-distance^6 units)
 110 r_m (distance units)
 111 r_{cut} (distance units) :ul
 112
 113 :line
 114
 115 [Mixing, shift, table, tail correction, restart, rRESPA info]:
 116
 117 For atom type pairs I,J and I != J, the epsilon and sigma coefficients
 118 and cutoff distance for all of the lj/cut pair styles can be mixed.
 119 The default mix value is {geometric}.  See the "pair_modify" command
 120 for details.
 121
 122 All of the {lj/cut} pair styles support the
 123 "pair_modify"_pair_modify.html shift option for the energy of the
 124 Lennard-Jones portion of the pair interaction.
 125
 126 The {lj/cut/coul/long} and {lj/cut/tip4p/long} pair styles support the
 127 "pair_modify"_pair_modify.html table option since they can tabulate
 128 the short-range portion of the long-range Coulombic interaction.
 129
 130 All of the {lj/cut} pair styles support the
 131 "pair_modify"_pair_modify.html tail option for adding a long-range
 132 tail correction to the energy and pressure for the Lennard-Jones
 133 portion of the pair interaction.
 134
 135 All of the {lj/cut} pair styles write their information to "binary
 136 restart files"_restart.html, so pair_style and pair_coeff commands do
 137 not need to be specified in an input script that reads a restart file.
 138
 139 The {lj/cut} and {lj/cut/coul/long} pair styles support the use of the
 140 {inner}, {middle}, and {outer} keywords of the "run_style
 141 respa"_run_style.html command, meaning the pairwise forces can be
 142 partitioned by distance at different levels of the rRESPA hierarchy.
 143 The other styles only support the {pair} keyword of run_style respa.
 144 See the "run_style"_run_style.html command for details.
 145
 146 :line
 147
 148 [Restrictions:]
 149
 150 These pair styles can only be used if LAMMPS was built with the
 151 USER-MISC package.  See the "Making LAMMPS"_Section_start.html#start_3
 152 section for more info on packages.
 153
 154 [Related commands:]
 155
 156 "pair_coeff"_pair_coeff.html
 157
 158 [Default:] none
 159
 160 :line
 161
 162 :link(Mei)
 163 [(Mei)] Mei, Davenport, Fernando, Phys Rev B, 43 4653 (1991)