diff --git a/doc/Eqs/pair_tersoff_zbl.jpg b/doc/Eqs/pair_tersoff_zbl.jpg index ac80bba44f..20d60d2256 100644 Binary files a/doc/Eqs/pair_tersoff_zbl.jpg and b/doc/Eqs/pair_tersoff_zbl.jpg differ diff --git a/doc/Eqs/pair_tersoff_zbl.tex b/doc/Eqs/pair_tersoff_zbl.tex index 4accac098e..902819aa1b 100644 --- a/doc/Eqs/pair_tersoff_zbl.tex +++ b/doc/Eqs/pair_tersoff_zbl.tex @@ -24,7 +24,7 @@ f_F(r_{ij}) & = & \frac{1}{1 + e^{-A_F(r_{ij} - r_C)}}\\ f_A(r) & = & -B \exp (-\lambda_2 r) \\ b_{ij} & = & \left( 1 + \beta^n {\zeta_{ij}}^n \right)^{-\frac{1}{2n}} \\ \zeta_{ij} & = & \sum_{k \neq i,j} f_C(r_{ik}) g(\theta_{ijk}) - \exp \left[ {\lambda_3}^3 (r_{ij} - r_{ik})^m \right] \\ + \exp \left[ {\lambda_3}^m (r_{ij} - r_{ik})^m \right] \\ g(\theta) & = & \gamma_{ijk} \left( 1 + \frac{c^2}{d^2} - \frac{c^2}{\left[ d^2 + (\cos \theta - \cos \theta_0)^2\right]} \right) diff --git a/doc/pair_tersoff_zbl.html b/doc/pair_tersoff_zbl.html index f0a6d381dc..5b59312f67 100644 --- a/doc/pair_tersoff_zbl.html +++ b/doc/pair_tersoff_zbl.html @@ -136,14 +136,20 @@ three-body parameters for SiCSi and SiSiC entries will not, in general, be the same. The parameters used for the two-body interaction come from the entry where the 2nd element is repeated. Thus the two-body parameters for Si interacting with C, comes from the -SiCC entry. By symmetry, the twobody parameters in the SiCC and CSiSi -entries should thus be the same. The parameters used for a particular +SiCC entry. +

+

The parameters used for a particular three-body interaction come from the entry with the corresponding three elements. The parameters used only for two-body interactions (n, beta, lambda2, B, lambda1, and A) in entries whose 2nd and 3rd element are different (e.g. SiCSi) are not used for anything and can be set to 0.0 if desired.

+

Note that the twobody parameters in entries such as SiCC and CSiSi +are often the same, due to the common use of symmetric mixing rules, +but this is not always the case. For example, the beta and n parameters in +Tersoff_2 (Tersoff_2) are not symmetric. +

We chose the above form so as to enable users to define all commonly used variants of the Tersoff portion of the potential. In particular, our form reduces to the original Tersoff form when m = 3 and gamma = @@ -154,14 +160,13 @@ different but equivalent form for alloys, which we will refer to as Tersoff_2 potential (Tersoff_2).

LAMMPS parameter values for Tersoff_2 can be obtained as follows: -gamma = 1, just as for Tersoff_1, but now lambda3 = 0 and the value of +gamma = omega_ijk, lambda3 = 0 and the value of m has no effect. The parameters for species i and j can be calculated using the Tersoff_2 mixing rules:

-

Values not shown are determined by the first atom type. Finally, the -Tersoff_2 parameters R and S must be converted to the LAMMPS +

Tersoff_2 parameters R and S must be converted to the LAMMPS parameters R and D (R is different in both forms), using the following relations: R=(R'+S')/2 and D=(S'-R')/2, where the primes indicate the Tersoff_2 parameters. @@ -260,11 +265,11 @@ of Ions in Matter' Vol 1, 1985, Pergamon Press.

-

(Albe) J. Nord, K. Albe, P. Erhartand K. Nordlund, J. Phys.: +

(Albe) J. Nord, K. Albe, P. Erhart and K. Nordlund, J. Phys.: Condens. Matter, 15, 5649(2003).

-

(Tersoff_2) J. Tersoff, Phys Rev B, 39, 5566 (1989) +

(Tersoff_2) J. Tersoff, Phys Rev B, 39, 5566 (1989); errata (PRB 41, 3248)

diff --git a/doc/pair_tersoff_zbl.txt b/doc/pair_tersoff_zbl.txt index 01e13ffeff..0ec5b69a24 100644 --- a/doc/pair_tersoff_zbl.txt +++ b/doc/pair_tersoff_zbl.txt @@ -132,14 +132,20 @@ three-body parameters for SiCSi and SiSiC entries will not, in general, be the same. The parameters used for the two-body interaction come from the entry where the 2nd element is repeated. Thus the two-body parameters for Si interacting with C, comes from the -SiCC entry. By symmetry, the twobody parameters in the SiCC and CSiSi -entries should thus be the same. The parameters used for a particular +SiCC entry. + +The parameters used for a particular three-body interaction come from the entry with the corresponding three elements. The parameters used only for two-body interactions (n, beta, lambda2, B, lambda1, and A) in entries whose 2nd and 3rd element are different (e.g. SiCSi) are not used for anything and can be set to 0.0 if desired. +Note that the twobody parameters in entries such as SiCC and CSiSi +are often the same, due to the common use of symmetric mixing rules, +but this is not always the case. For example, the beta and n parameters in +Tersoff_2 "(Tersoff_2)"_#Tersoff_2 are not symmetric. + We chose the above form so as to enable users to define all commonly used variants of the Tersoff portion of the potential. In particular, our form reduces to the original Tersoff form when m = 3 and gamma = @@ -150,13 +156,12 @@ different but equivalent form for alloys, which we will refer to as Tersoff_2 potential "(Tersoff_2)"_#Tersoff_2. LAMMPS parameter values for Tersoff_2 can be obtained as follows: -gamma = 1, just as for Tersoff_1, but now lambda3 = 0 and the value of +gamma = omega_ijk, lambda3 = 0 and the value of m has no effect. The parameters for species i and j can be calculated using the Tersoff_2 mixing rules: :c,image(Eqs/pair_tersoff_2.jpg) -Values not shown are determined by the first atom type. Finally, the Tersoff_2 parameters R and S must be converted to the LAMMPS parameters R and D (R is different in both forms), using the following relations: R=(R'+S')/2 and D=(S'-R')/2, where the primes indicate the @@ -253,8 +258,8 @@ units. of Ions in Matter' Vol 1, 1985, Pergamon Press. :link(Albe) -[(Albe)] J. Nord, K. Albe, P. Erhartand K. Nordlund, J. Phys.: +[(Albe)] J. Nord, K. Albe, P. Erhart and K. Nordlund, J. Phys.: Condens. Matter, 15, 5649(2003). :link(Tersoff_2) -[(Tersoff_2)] J. Tersoff, Phys Rev B, 39, 5566 (1989) +[(Tersoff_2)] J. Tersoff, Phys Rev B, 39, 5566 (1989); errata (PRB 41, 3248)