# How To calculate formation energy of Ni doped ZnO system in Quantum ESPRESSO?

I have successfully completed the SCF of doped system Ni doped ZnO using GGA+U. I have checked the SCF.out file and found the energy value : E= -6112.0389 Ry.

Hwo to calculate the formation energy of doped system using this $$formula^1$$ ?

$$E_f = E(doped_{ZnO}) -E_{ZnO}+2(\mu_{Zn}+\mu_{Ni})$$

Where I can find the chemical value of Zn and Ni ($$\mu_{Zn},\mu_{Ni}$$) in my outputfil.

Also I would like to know what is the minimum formation energy for a crystal to be stable?

Based on reference paper given above, the reference system of interest here is $$\ce{Zn_{22}M_{2}O_{24}}$$. As two $$\ce{Zn}$$ atoms were substituted by transition metal ($$M=\ce{V}, \ce{Cr}, \ce{Mn}, \ce{Fe}, \ce{Co}, \ce{Ni}).$$ Formation energy for doping is given as $$E_f = E(doped_{\ce{ZnO}}) -E_{\ce{ZnO}}+2(\mu_{\ce{Zn}}+\mu_{\ce{Ni}})$$ where $$E(doped_{\ce{ZnO}})$$ is total energy of $$\ce{Zn_{22}M_{2}O_{24}}$$ and $$E_{\ce{ZnO}}$$ is energy of $$\ce{Zn_{24}O_{24}}$$. The chemical potentials of the $$\ce{V}$$, $$\ce{Cr}$$, $$\ce{Mn}$$, or $$\ce{Fe}$$ are energies of single atom in the body centered cubic structure, where as in case of $$\ce{Co}$$ and $$\ce{Ni}$$, hexagonal crystal structure and face centered cubic structure were considered respectively. For $$\ce{Zn}$$ they have considered considered two state, one $$\ce{Zn}$$ as its elementary substance and second as $$\mu_{\ce{Zn}}=E_{\ce{ZnO}}-\mu_{\ce{O}}$$ where $$\mu_{\ce{O}}$$ is energy of oxygen as elementary substance. For stability of any structure you can do mechanical(elastic constants) and dynamic stability(phonon) calculations. Based on formation energy one can say that if it is negative, such doping are more favorable.
• +1 But please take a look at my edit where I fixed the molecular formula for $\ce{Zn22M2O24}$ and try yourself to fix the rest of the atoms/chemical-formulas! May 12 at 17:00