First-principles study of phase stability of Gd-doped EuO and EuS

J. M. An, S. V. Barabash, V. Ozolins, M. van Schilfgaarde, and K. D. Belashchenko

Phys. Rev. B 83, 064105 – Published 17 February 2011

Abstract

Phase diagrams of isoelectronic Eu $_{1 - x}$ Gd $_{x}$ O and Eu $_{1 - x}$ Gd $_{x}$ S quasibinary alloy systems are constructed using first-principles calculations combined with a standard cluster-expansion approach and Monte Carlo simulations. The oxide system has a wide miscibility gap on the Gd-rich side but forms ordered compounds on the Eu-rich side, exhibiting a deep asymmetric convex hull in the formation enthalpy diagram. The sulfide system has no stable compounds. The large difference in the formation enthalpies of the oxide and sulfide compounds is attributed to the contribution of local lattice relaxation, which is sensitive to the anion size. The solubility of Gd in both EuO and EuS is in the range of 10%–20% at room temperature and quickly increases at higher temperatures, indicating that highly doped disordered solid solutions can be produced without the precipitation of secondary phases. We also predict that rocksalt GdO can be stabilized under appropriate experimental conditions.

Received 22 May 2010

DOI:https://doi.org/10.1103/PhysRevB.83.064105

Authors & Affiliations

J. M. An^1,*, S. V. Barabash², V. Ozolins², M. van Schilfgaarde³, and K. D. Belashchenko¹

¹Department of Physics and Astronomy, Nebraska Center for Materials and Nanoscience, University of Nebraska, Lincoln, Nebraska 68588, USA
²Department of Materials Science and Engineering, University of California, Los Angeles, California 90095-1595, USA
³Arizona State University, Tempe, Arizona 85284, USA

^*anj_777@msn.com

Article Text (Subscription Required)

Click to Expand

References (Subscription Required)

Click to Expand

Issue

Vol. 83, Iss. 6 — 1 February 2011

Reuse & Permissions

Access Options

Article Available via CHORUS

Download Accepted Manuscript

Author publication services for translation and copyediting assistance advertisement

Images

Figure 1
Electronic band structure of EuO calculated by (a) $GGA + U$ and (b) QS $GW$ , and of ferromagnetic phase GdO, calculated by (c) $GGA + U$ and (d) QS $GW$ . Black (red) curves correspond to the majority (minority) spin bands. Energy is referenced either from the valence band maximum (VBM) or from the Fermi level ( $E$ $_{F}$ ). In GdO, states near $E_{F}$ at $Γ$ are of Gd $d$ character; the band whose value is $-$ 1.5 eV at $Γ$ is of $s$ character. If correlations were strong, the $d$ band would become narrow and possibly shift relative to the $s$ band. As can be seen, the GGA and QS $GW$ results are very similar for these bands: the GGA and QS $GW$ results differ mainly in the positions of the O $2 p$ bands, at around $-$ 6 eV. That O $2 p$ states shift downward relative to results in the GGA (or LDA) seems to be a universal property of oxide insulators.Reuse & Permissions
Figure 2
Density of states (DOS) of FM EuO and EuS, AFM type II GdO, and AFM type II GdS calculated in the $GGA + U$ . For the AFM II phases, the solid line shows the total DOS including both the cation and the anion contributions, whereas the dashed lines show the partial majority- and minority-spin DOSs from the Gd cations. Energy is referenced from the Fermi level $E_{F}$ .Reuse & Permissions
Figure 3
Three different types of ECI parameters as functions of effective radius in units of the nearest-cation-neighbor distance $r_{NN}$ in the cluster expansions for Eu $_{x}$ Gd $_{1 - x}$ O, denoted by red crosses and Eu $_{x}$ Gd $_{1 - x}$ S, denoted by blue circles.Reuse & Permissions
Figure 4
Formation enthalpy per cation vs composition $x$ for all distinct cation orderings within rocksalt structure for (a) Eu $_{x}$ Gd $_{1 - x}$ O with up to 30 atoms per unit cell and (b) Eu $_{x}$ Gd $_{1 - x}$ S with up to 20 atoms per unit cell. The black open circles are the first-principles inputs, the red crosses are the fitted CE values for the input structures, and the gray dots are the predicted $Δ H_{CE}$ for all other structures.Reuse & Permissions
Figure 5
Calculated phase diagrams obtained using Monte Carlo simulations from the cluster expansions for (a) Eu $_{x}$ Gd $_{1 - x}$ O and(b) Eu $_{x}$ Gd $_{1 - x}$ S alloys. The phase labels indicate the ordering within the cation sublattice. In (a), the phases in the region bounded by the dot-dashed line have not been identified. In (b), the solid lines correspond to the paramagnetic high-temperature phases, while the dotted lines show the prediction of the $T = 0$ ferromagnetic cluster expansion.Reuse & Permissions
Figure 6
Monte Carlo heating simulation of the $L$ 1 $_{1}$ -to-disordered phase transition for Eu $_{0.5}$ Gd $_{0.5}$ O. Formation enthalpy and heat capacity are shown by black and red lines, respectively. The inset shows the temperature dependence of the composition $x$ .Reuse & Permissions

Physical Review B

covering condensed matter and materials physics