MEAM Potential for Au-Si

MEAM Potential for Au-Si

Adriano Satana and Wei Cai

Created Aug, 2015, Last modified Sep, 2015

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

---

Potential for Pure Elements

MEAM Potential for Au

We use the 'Au' potential whose parameters are originally given in M. I. Baskes, Phys. Rev. B 46, 2727 (1992), and later modified by S. Ryu, et al. Model. Simul. Mater. Sci. Eng. 17, 075008 (2009).

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 'Au' 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   
'Au' 'fcc' 12.     79     196.967 6.34090112  5.449   2.20    6     2.20  

$(R_{i}^{0})$ $E_{i}^{0}$   $A_{i}$   $t_{i}^{(0)}$   $t_{i}^{(1)}$        $t_{i}^{(2)}$         $t_{i}^{(3)}$     $\rho _{0}^{\rm {Au}}$
 alat  esub  asub t0     t1          t2           t3     rozero  ibar
 4.07  3.93  1.04 1.0  1.58956328 1.50776392  2.60609758    1.      3

Note that the nearest neighbor distance $R_{i}^{0}$ = alat / ${\sqrt {2}}$.

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

ibar is a setting used in the equation of state (EOS), and will be explained later.

MEAM Potential for Si

We use the 'Si4' potential whose parameters are originally given in M. I. Baskes, Phys. Rev. B 46, 2727 (1992), and later modified by S. Ryu, et al. Model. Simul. Mater. Sci. Eng. 17, 075008 (2009).

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 'Siz' 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   
'Si4' 'dia' 4.     14     28.086     4.87     4.4     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)}$     Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \rho_0^{\rm Si}}

 alat  esub  asub t0     t1          t2           t3     rozero  ibar
 5.431 4.63  1.  1.0    3.13        4.47          -1.8      1.48    0

Note that the nearest neighbor distance Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle R_i^0 } = alat Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \times \sqrt{3}/4} for the diamond cubic structure.

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \rho_0^{\rm Si}} = rozero will be important only for cross-potential.

ibar is a setting used in the equation of state (EOS), and will be explained later.

Cross Potential between Au and Si

ABC