MEAM Potential for Au-Ge

MEAM Potential for Au-Ge

Adriano Santana and Wei Cai

Created Aug, 2015, Last modified Sep, 2015

This tutorial explains how to specify the parameters for the Au-Ge MEAM potential in MD++. It starts with the parameters in pure Au and pure Ge potentials, then talks about the Au-Ge cross potential.

---

Potential for Pure Elements

Original MEAM Potential for Au

The details for the original 'Au' potential can be found here JUST REFER TO OTHER PAGE....!

UNDER CONSTRUCTION!!!!!

MEAM Potential for Ge

START HERE...!!!!

We use the 'Ge' potential whose parameters are originally given in M. I. Baskes, The main parameters in the MEAM potential is specified in the meamf file. (In MD++, this file is in the potentials/MEAMDATA folder.) The lines correspond to 'Ge' is given below. Most of these parameters correspond to Table III of Baskes PRB (1992), as shown below.

                                      $\alpha _{i}$       $\beta _{i}^{(0)}$     $\beta _{i}^{(1)}$   $\beta _{i}^{(2)}$  $\beta _{i}^{(3)}$
elt  lat   z    ielement   atwt      alpha    b0       b1     b2    b3      
'Ge' 'dia'  4.     32       72.64     4.98   4.55      5.5    5.5  5.5 

$(R_{i}^{0})$    $E_{i}^{0}$   $A_{i}$   $t_{i}^{(0)}$   $t_{i}^{(1)}$         $t_{i}^{(2)}$          $t_{i}^{(3)}$       $\rho _{0}^{\rm {Ge}}$
 alat  esub  asub t0     t1          t2           t3     rozero  ibar
 5.6575 3.85  1.0 1.0   4.02       5.23          -1.6     1.35     0

Note that the nearest neighbor distance $R_{i}^{0}$ = alat $\times {\sqrt {3}}/4$ for the diamond cubic structure.

$\rho _{0}^{\rm {Ge}}$ = rozero will be important only for cross-potential.

ibar is a setting used in the equation of state (EOS). ibar selects the G(gamma) function in Eq (4) and (5) of the paper by BJ Lee, PRB 68, 144112 (2003).

While the functional form is quite different, the modulus is almost not affected by the choice of ibar.

The parameters for the cross potential are specified in the AuGe2nn.meam file. The variables in Eq.(A.1) of Ryu and Cai JPCM (2010) are given in the parenthesis.

erose_form = 3
rc = 4.5
attrac(2,2) = -0.36 ($\gamma$)
repuls(2,2) = 16.0  ($\lambda$)
Cmin(2,2,2) = 1.85  ($C_{\rm {min}}$)

Note that we label the atomic species of Si as 2.

Cross Potential between Au and Ge

The parameters for the cross potential are specified in AuGe2nn.meam file. The lines relevant for the cross potential (i.e. between species 1 and 2) are shown below. The values are to be published soon. THey are calculated from VASP LDA/PAW.

The values correspond to Table 3 of Ryu and Cai, J. Phys. Condens. Matter, 22, 055401 (2010).

RECALCULATE THIS VALUES FOR MY POTENTIAL AU-GE!!!

re(1,2) = 2.6495 new!    ($r_{e}$)
delta(1,2) = 0.048 new!  (related to $E_{c}$, see below)
lattce(1,2) = b1
alpha(1,2) = 4.927 new!     ($\alpha$)
attrac(1,2) = 0.0      
repuls(1,2) = 0.26     ($\gamma$)
Cmin(1,1,2) = 1.9     ($C_{\min }(1,1,2)$)
Cmin(1,2,1) = 0.95     ($C_{\min }(1,2,1)$)
Cmin(1,2,2) = 1.85     ($C_{\min }(1,2,2)$)
Cmin(2,2,1) = 1.0     ($C_{\min }(2,2,1)$)

The values for $E_{c}({\rm {AuGe}})=4.734$. This value is related to delta(1,2) through

$E_{c}({\rm {AuGe}})=0.5*[E_{c}({\rm {Au}})+E_{c}({\rm {Ge}})]-{\rm {delta}}(1,2)=0.5*(4.394+5.17)-0.048=4.734$.

$\rho _{0}^{\rm {Ge}}/\rho _{0}^{\rm {Au}}$ = 1.35 because of the $\rho _{0}^{\rm {Ge}}$ and $\rho _{0}^{\rm {Au}}$ values specified above.

Cmax = 2.8 is the default value.