Abstract
Multiferroics in which giant ferroelectric polarization and magnetism coexist are of tremendous potential for engineering disruptive applications in information storage and energy conversion. Yet the functional properties of multiferroics are thought to be affected detrimentally by the presence of point defects, which may be abundant due to the volatile nature of some constituent atoms and the high temperatures involved in the synthesis of materials. Here, we demonstrate with theoretical methods that oxygen vacancies may enhance the functionality of multiferroics by radically changing their magnetic interactions in thin films. Specifically, oxygen vacancies may restore missing magnetic super-exchange interactions in large axial ratio phases, leading to full antiferromagnetic spin ordering, and induce the stabilization of ferrimagnetic states with considerable net magnetizations. Our theoretical study should help to clarify the origins of long-standing controversies in bismuth ferrite and improve the design of technological applications based on multiferroics.
Similar content being viewed by others
Introduction
Finding multiferroics in which ferroelectricity and magnetism coexist and influence each other is of great fundamental and applied interests1,2. Salient technological features of multiferroics include the control of the magnetization with electric fields to design efficient logic and memory devices3,4, and the realization of large piezomagnetic coefficients to facilitate the miniaturization of antennas and sensors5,6. Furthermore, competition between phases displaying distinct electric polarization and magnetic ordering offers also encouraging prospects for energy conversion applications like photovoltaics and solid-state cooling7,8,9,10.
But, unfortunately, multiferroics are rare in nature, typically present weak magnetoelectric coupling (BiFeO3)11, and require extreme thermodynamic conditions for their synthesis (PbVO3 and BiCoO3)12,13,14. In addition, magnetoelectric multiferroics mostly are antiferromagnetic hence potential applications based on external magnetic bias generally are frustrated due to their little influence on antiparallel magnetic spins15. Common strategies employed to synthesize bettered multiferroic materials include doping16,17, solid solutions18,19, and strain engineering1,20.
Through epitaxial strain, it is actually possible to create new multiferroic materials in the laboratory that exhibit giant electric polarization and unexpected magnetic spin ordering21,22,23. An illustrative example is given by BiFeO3 (BFO) thin films, in which large spontaneous polarization and ferromagnetism (FM) have been observed under moderate compressive biaxial strains at room temperature24,25. The origins of the net magnetization reported for BFO thin films, however, are not clear yet and from a technological point of view it is crucial to understand them correctly.
Most magnetic ferroelectrics with chemical formula AMO3 and perovskite-like structure present antiferromagnetic spin ordering along the three pseudo-Cartesian directions (AFM-G), due to the dominant role of oxygen-mediated super-exchange interactions between neighbouring transition metal atoms M26,27,28,29. In large axial ratio structures (c/a > 1), also known as “super-tetragonal”, the covalency of M–O bonds oriented along the electric polarization direction is significantly reduced and consequently the magnetic exchange interactions, which favour parallel magnetic spins, dominate in that direction. The coexistence of “in-plane” antiferromagnetism and “out-of-plane” ferromagnetism in super-tetragonal phases (AFM-C, Fig. 1a, b) leads anyway to null crystal magnetization when small spin canting effects are neglected. Therefore, intrinsic and robust FM in principle is not expected to occur in BFO or any other similar multiferroic30,31,32,33.
a Representation of the super-tetragonal phase (\({\mathcal{T}}\), space group P4mm) characteristic of bulk BCO and other multiferroics exhibiting giant electric polarization. b Magnetic structure in the \({\mathcal{T}}\) phase rendering C-type antiferromagnetism (AFM-C); green arrows represent atomic magnetic moments and their orientation. c Examples of substrates in which to grow BCO thin films displaying the effects predicted in this study (bottom); other well-known materials exhibiting super-tetragonal phases are shown for comparison (top)21,24. “SE” stands for magnetic super-exchange interactions. Bi, Co, and O atoms are represented with magenta, blue, and red spheres, respectively.
A plausible explanation for the appearance of FM in AMO3 perovskite oxide thin films is based on extrinsic causes like point defects25,34. The volatile nature of bismuth and the high temperatures involved in the preparation of samples, for instance, make the presence of oxygen vacancies (VO) almost inevitable in Bi-based multiferroics17,19. In fact, oxygen defects may modify significantly the structural and functional properties of perovskite thin films via changes in the M oxidation states and their coupling with the lattice strain35,36,37,38,39,40. However, a number of theoretical works based on first-principles methods have agreed in that the combined action of VO and lattice strain may not affect considerably the magnetic properties of Bi-based multiferroics41,42,43.
Here, we present new theoretical evidence showing that the presence of VO may in fact change radically the magnetic properties of multiferroic thin films via previously overlooked electro-structural mechanisms. We select BiCoO3 (BCO) as the model multiferroic in which to perform first-principles calculations based on density functional theory (DFT) because (i) this material already exhibits a large axial ratio in the absence of any strain, (ii) the magnetic effects that we predict could be realized on substrates that are commonly employed for growth of epitaxial perovskite thin films (Fig. 1c), and (iii) BCO naturally exhibits a high concentration of oxygen vacancies (i.e., it decomposes easily at temperatures around 700 K and needs to be synthesised at high pressures13,14). In particular, it is found that oxygen vacancies occupying specific lattice positions can induce the stabilization of full antiferromagnetic (AFM-G) super-tetragonal and ferrimagnetic (FiM) monoclinic polar phases, depending on the lattice strain. As a consequence, phase competition is enriched and magnetic functionalities further enhanced in comparison to perfectly stoichiometric thin films. We show that most of the results obtained in BCO thin films can be generalized to BFO and other Bi-based multiferroics, hence our conclusions are of broad applicability and significance to the field of multifunctional materials.
Results
We have performed first-principles simulations of stoichiometric and non-stoichiometric BCO thin films based on density functional theory (DFT) techniques (see “Methods” section). The VO-induced effects explained in this section have been observed for compositions BiCoO2.75 and BiCoO2.875. Detailed discussions on the influence of the selected DFT exchange-correlation energy functional and other technical parameters on our general conclusions can be found in the “Methods” and Supplementary Discussion sections. Overall, the results presented next are robust and qualitatively accurate.
V O-induced effects on phase competition and functionality
It is known both from experiments and first-principles calculations that bulk BCO presents a polar tetragonal \({\mathcal{T}}\) phase (space group P4mm) with a large axial ratio of c/a ≈ 1.3 and relatively small lattice paramater a = 3.76 Å9,14,31,44 (Fig. 1a). The competing structures are a non-polar orthorhombic \({\mathcal{O}}\) phase (space group Pnma) and a polar monoclinic \({\mathcal{M}}\) phase (space group Pc) (Supplementary Fig. 1, Supplementary Tables 1–3)9,14,31,44; both competing phases have cells that are slightly distorted versions of the ideal cubic perovskite structure with c/a ≈ 1.
The polar phases in BCO present spontaneous polarizations along quite different crystallographic directions, namely, pseudocubic [001]pc in \({\mathcal{T}}\) and ~ [111]pc in \({\mathcal{M}}\). As regards magnetism, recent first- principles calculations have shown that phases \({\mathcal{O}}\) and \({\mathcal{M}}\) exhibit G-type antiferromagnetism (AFM-G) with a quite high Néel temperature, TN ≈ 500 K, whereas the \({\mathcal{T}}\) phase presents C-type antiferromagnetism (AFM-C) with a relatively low TN of ≈ 310 K31. In stoichiometric BCO thin films, and by completely neglecting temperature effects, a multiferroic \({\mathcal{T}}\to {\mathcal{M}}\) phase transition involving large structural, polar, and magnetic changes has been predicted to occur at in-plane parameter ain = 3.91 Å9 (Fig. 2a–c).
Results were obtained for composition BiCoO2.75. a Zero-temperature energy of competing phases expressed as a function of in-plane lattice parameter. Metastable phases are indicated by grey curves and strain-induced phase transitions by vertical dashed lines; the oxygen vacancy positions leading to the lowest energies, either apical “Ap” or equatorial “Eq”, are indicated within parentheses. “G” stands for G-type antiferromagnetism, “C” for C-type antiferromagnetism, and “FiM” for ferrimagnetism. b Electric polarization of stoichiometric and non-stoichiometric ground-state phases. c Magnetic transition temperature of stoichiometric and non-stoichiometric ground-state phases. d Phase transition sequence occurring in non-stoichiometric BCO thin films under increasing ain; the green arrows represent the magnetic spin ordering in each phase. e Change in volume, ΔΩ, change in electric polarization orientation, Δα, and change in magnetic moment per formula unit, Mi, associated with the multiferroic \({\mathcal{T}}\)-G → \({\mathcal{M}}\)-FiM phase transition.
Figure 2a–c show the influence of neutral oxygen vacancies (see “Methods” section), VO, on the structural, ferroelectric, and magnetic properties of BCO thin films (the accompanying changes in atomic lattice positions and energy band gap are reported in Supplementary Tables 4–7 and Supplementary Fig. 2). For the smallest in-plane lattice parameters, a \({\mathcal{T}}\)-C phase (magnetic spin ordering is indicated along with the structure symmetry) containing oxygen vacancies in equatorial (Eq) positions (Fig. 1a) renders the lowest energy. The electric polarization and Néel temperature calculated for \({\mathcal{T}}\)-C(Eq) are significantly lower than for the analogous stoichiometric phase, in particular, we estimate differences of ΔP ≈ −75 μC cm−2 and ΔTN ≈ −50 K when considering identical in-plane parameters.
At ain = 3.77 Å, an unusual magnetic phase transition from AFM-C to AFM-G spin ordering occurs that is accompanied by the appearance of a small in-plane electric polarization (Pxy ~ 10 μC cm−2) and change in VO symmetry. The Néel temperature of the \({\mathcal{T}}\)-G(Ap) phase is lower than that of \({\mathcal{T}}\)-C(Eq) by approximately 50 K. Furthermore, at ain ≥ 3.91 Å the system adopts a monoclinic ferrimagnetic (FiM) phase with oxygen vacancies preferably in equatorial positions, \({\mathcal{M}}\)-FiM(Eq), and a considerable net magnetization (i.e., 0.5 μB per formula unit for BiCoO2.75—Fig. 2d—and 0.25 μB for BiCoO2.875). The Néel temperature of the \({\mathcal{M}}\)-FiM(Eq) phase is larger than that of \({\mathcal{T}}\)-G(Ap) and remains close to room temperature almost independently of ain.
Interestingly, the formation energy of neutral oxygen vacancies in BCO thin films (see “Methods” section) is found to decrease steadily for ain ≥ 3.77 Å as compared to the bulk equilibrium system (Supplementary Fig. 3). For instance, at ain = 3.95 Å the formation energy of the most stable type of VO is found to be 0.25 eV smaller than obtained at 3.76 Å (per formula unit). The difference in formation energy between Ap and Eq oxygen vacancies can be deduced directly from the energy curves shown in Fig. 2a (see “Methods” section). For instance, in the \({\mathcal{M}}\) phase the creation of an equatorial VO is ~0.1 eV more favourable than an apical vacancy, almost independently of the in-plane lattice parameter.
The physical mechanisms responsible for the two multiferroic phase transitions represented in Fig. 2d will be explained in detail in the next subsections. Let us now comment briefly on the functionality enhancement deriving from the \({\mathcal{T}}\)-G(Ap) → \({\mathcal{M}}\)-FiM(Eq) transformation. First, a large change in the electric polarization orientation involving a rotation of ≈ 60∘ is observed during the transition (Fig. 2e); as a consequence, and in analogy to what has been observed in Pb(Zr1−xTix)O3 alloys45,46 and Bi(Fe1−xCox)O3 thin films47, it should be possible to realize large piezoelectric responses under small electric bias for ain ≈ 3.91 Å. Second, the sizeable changes in electric polarization and total magnetization in principle should allow for control of the polarization with magnetic fields and vice versa, which hints at the likely existence of large magnetoelectric couplings1,2. And third, the thin film dimension along the direction perpendicular to the epitaxial plane shrinks by an impressive ≈11% (Fig. 2e) hence there is the possibility of realizing giant piezomagnetic responses5,6 and multicaloric effects8,9,10 through the application of external bias. As it will be shown later, similar magnetic phenomena than found for BCO are likely to occur also in other Bi-based multiferroic thin films, including BiFeO3.
V O-induced magnetic super-exchange interactions in the \({\mathcal{T}}\) phase
Figure 3 summarizes the electronic, structural, and magnetic properties of stoichiometric and non-stoichiometric \({\mathcal{T}}\) BCO thin films. In the stoichiometric \({\mathcal{T}}\)-C phase (Fig. 3a), the square-pyramidal O5 crystal field splits the electronic Co d levels into nondegenerate b2g (dxy), doubly degenerate eg (dxz, dyz), and nondegenerate a1g (\({d}_{{z}^{2}}\)) and b1g (\({d}_{{x}^{2}-{y}^{2}}\)). Our first-principles calculations render a high-spin Co state characterised by the electronic occupation \({b}_{2g}^{2}{e}_{g}^{2}{a}_{1g}^{1}{b}_{1g}^{1}\) and magnetic moment of 3.1 μB (in absolute value), in good agreement with the available experimental data14.
Results were obtained for composition BiCoO2.75. a Stoichiometric \({\mathcal{T}}\) phase. The first moment of the square of the d densities of states (DOS) are indicated with thick vertical lines; from those moments the represented electronic orbital occupation sketch follows (idem for the following subfigures). Spin-up densities are represented in the top DOS panel and spin-down in the bottom DOS panel (idem for the following subfigures). b Non-stoichiometric \({\mathcal{T}}\) thin films with neutral oxygen vacancies (VO) in Eq positions and AFM-C spin ordering. The orange plane indicates the two Co ions that are reduced as a consequence of creating such VO. The electronic orbital occupation sketch corresponds to the Co ion highlighted with blue colour in the transparency. The purple arrows in the d-orbitals occupation sketch indicate electronic delocalization. c Non-stoichiometric \({\mathcal{T}}\) thin films with neutral VO's in Ap positions and AFM-G spin ordering. The orange plane indicates the two Co ions that are reduced as a consequence of creating such VO. The electronic orbital occupations correspond to the Co ions highlighted with blue colour in the transparencies. The red arrows in the d-orbitals occupation sketch indicate differences with respect to the stoichiometric case. d Charge density surface plot for the stoichiometric \({\mathcal{T}}\)-C phase; the surface over which the charge density is calculated is indicated by a grey plane in the structural ball-stick representation. Isovalue paths are represented with black and coloured lines. e Charge density surface plot for non-stoichiometric \({\mathcal{T}}\)-G(Ap). Regions of interest describing Co-O bonds are indicated with yellow arrows.
In the non-stoichiometric \({\mathcal{T}}\)-C(Eq) phase (Fig. 3b), the relative splitting of electronic d levels remains invariant with respect to the stoichiometric case and the electronic occupation of the two cobalt ions nearest to the neutral VO, which are reduced and are electronically equivalent (indicated with an orange plane in Fig. 3b), changes slightly to \({b}_{2g}^{2}{e}_{g}^{3}{a}_{1g}^{1}{b}_{1g}^{1}\). Interestingly, when VO is created in an apical position the magnetic spin ordering of the \({\mathcal{T}}\) phase changes to AFM-G at specific ain’s (Fig. 2a). Figure 3c shows the electronic orbital occupation deduced from our DFT calculations for the two reduced cobalt ions in the \({\mathcal{T}}\)-G(Ap) phase (indicated with an orange plane in Fig. 3c), which in this case turn out to be electronically inequivalent. In particular, the doubly degenerate eg (dxz, dyz) orbitals in the cobalt ion closest to the apical VO (Co3, as labelled in Fig. 3d, e) undergo a significant energy reduction and become fully populated rendering the occupation state \({e}_{g}^{4}{b}_{2g}^{1}{a}_{1g}^{1}{b}_{1g}^{1}\), while the other reduced metal ion (Co1, as labelled in Fig. 3d, e) exhibits the more usual distribution \({b}_{2g}^{2}{a}_{1g}^{2}{e}_{g}^{2}{b}_{1g}^{1}\).
These drastic electronic rearrangements are correlated with the appearance of a strong structural distortion in the system that pushes Co3 towards the oxygen atom underneath of it (see inversion of the corresponding O5 square-pyramid in Fig. 3c, e), which tends to restore (partially) the missing magnetic super-exchange interactions along the out-of-plane direction. This super-exchange restoration mechanism leading to the stabilization of AFM-G spin ordering, which is accompanied by a covalency increase in the out-of-plane Co1–O–Co3 bonds, is clearly observed in the electronic density plots shown in Fig. 3d–e (see yellow arrows therein). We have explicitly checked that the stabilization of the \({\mathcal{T}}\)-G(Ap) phase is robust with respect to the concentration of oxygen vacancies and presence of finite-size bias in our simulations (see “Methods” and Supplementary Discussion).
The VO-induced AFM-C → AFM-G phase transition disclosed in large axial ratio BCO thin films may shed some light on uncomprehended experimental observations of antiferromagnetic spin ordering in other super-tetragonal multiferroic systems. For instance, first-principles works agree in predicting AFM-C spin ordering for \({\mathcal{T}}\) BFO thin films48,49,50 whereas in most experimental studies AFM-G is reported for the same system51,52. As it will be shown later, by explicitly considering the presence of oxygen vacancies in BFO thin films such disagreements between theoretical and experimental works can be reconciled.
V O-induced stabilization of a ferrimagnetic \({\mathcal{M}}\) phase
Figure 4 summarizes the electronic, structural, and magnetic properties of stoichiometric and non-stoichiometric \({\mathcal{M}}\) BCO thin films. In the stoichiometric \({\mathcal{M}}\)-G phase (Fig. 4a), the octahedral O6 crystal field splits the electronic Co d levels into doubly degenerate eg (\({d}_{{z}^{2}}\), \({d}_{{x}^{2}-{y}^{2}}\)) and triply degenerate t2g (dxy, dyz, dxz). A strong Jahn-Teller distortion rendering a large Q2 value of 0.37 Å53 lifts further the degeneracy in the eg and t2g manifolds54, promoting the electronic occupation state \({b}_{2g}^{2}{e}_{g}^{2}{a}_{1g}^{1}{b}_{1g}^{1}\) (Fig. 4a).
Results were obtained for composition BiCoO2.75. a Stoichiometric \({\mathcal{M}}\)-G phase. The first moment of the square of the d densities of states (DOS) are indicated with thick vertical lines; from those moments the represented electronic orbital occupation sketch follows (idem for the following subfigures). Spin-up densities are represented in the top DOS panel and spin-down in the bottom DOS panel (idem for the following subfigures). b Non-stoichiometric \({\mathcal{M}}\) thin films with neutral oxygen vacancies (VO) in Eq positions and FiM spin ordering. The orange plane indicates the two equivalent Co ions that are reduced as a consequence of creating such VO. The electronic orbital occupation sketch corresponds to the Co ions highlighted with blue color in the transparency. The red arrows in the d-orbitals occupation sketch indicate the difference with respect to the stoichiometric case. c Spin-up (dark green) and spin-down (light green) electronic charge densities calculated for stoichiometric \({\mathcal{M}}\) BCO thin films considering two perpendicular planes; highly magnetized oxygen atoms are indicated with red squares. Stabilization of FiM spin ordering in the \({\mathcal{M}}\)-FiM(Eq) phase is achieved when a highly magnetized equatorial oxygen atom is removed from the stoichiometric system. d Charge density surface plot for the stoichiometric \({\mathcal{M}}\)-G phase; the surface over which the charge density is calculated is indicated by a grey plane in the structural ball-stick representation. Isovalue paths are represented with black and coloured lines. e Charge density surface plot for non-stoichiometric \({\mathcal{M}}\)-FiM(Eq). Regions of interest describing structural distortions and Co-O bonds are indicated with yellow arrows.
Remarkably, when specific VO’s are created in equatorial positions the lowest-energy magnetic spin ordering becomes FiM. In particular, the net magnetization per formula unit associated with the \({\mathcal{M}}\)-FiM(Eq) phase is considerable: 0.5 and 0.25 μB for the off-stoichiometries BiCoO2.75 and BiCoO2.875, respectively. The two Co ions that are reduced by a neutral vacancy present the same magnetic moment orientation, equivalent electronic orbital occupancy, and sit within the [111]pc plane (indicated with an orange plane in Fig. 4b). Interestingly, the \({\mathcal{M}}\)-FiM(Eq) phase appears only when highly magnetized oxygen atoms (0.1–0.2 μB) occupying equatorial positions in the stoichiometric \({\mathcal{M}}\)-G crystal are removed (Fig. 4c).
For the pair of distant ions Co1 and Co4 to get reduced (Co3+ → Co2+, Fig. 4c), the crystal needs to undergo sizable structural distortions involving the Bi and O atoms that are close to the equatorial VO (Supplementary Tables 6–7). For example, the non-magnetic oxygen in the apical position above Co2 is displaced significantly out of Co2–Co4 vertical (Fig. 4c). Such a concerted electro-structural process has a clear imprint on the electronic density calculated for the plane oriented perpendicular to the substrate that contains Co1, Co2, Co3, and Co4 (Fig. 4d, e, see yellow arrows therein). As can be appreciated in Fig. 4d, e, the covalency of the Co4–O bond in the \({\mathcal{M}}\)-FiM(Eq) phase is reduced significantly as compared to that in the stoichiometric system. The vertical tilt of the Co2–O bond, which brings the apical oxygen closer to the equatorial VO, is also evidenced in Fig. 4e. We have carried out several numerical tests to assess the influence of the dimensions of the simulation cell and concentration of oxygen vacancies on the findings just described. Our results indicate that the revealed VO-induced stabilization of the \({\mathcal{M}}\)-FiM(Eq) phase is robust and general (see “Methods” and Supplementary Discussion).
The discovery of FiM spin ordering in the non-stoichiometric \({\mathcal{M}}\) phase motivated us to search for similar magnetic states, even if metastable, in the \({\mathcal{T}}\) and \({\mathcal{O}}\) BCO thin film geometries. The presence of highly magnetized O atoms was acknowledged in both stoichiometric phases (Supplementary Fig. 7), however, upon removal of those oxygens no net magnetization emerged (neglecting small spin canting effects). These results suggest that the structural distortions just described are key to the stabilization of FiM spin ordering. For instance, in the \({\mathcal{T}}\) phase the non-magnetic oxygens appear in equatorial positions that are fairly clamped to the substrate (Supplementary Fig. 7), hence they hardly can move off-symmetry. Meanwhile, in the \({\mathcal{O}}\) phase the lack of polar order and high dielectric permittivity appears to impose high energy penalties to any symmetry-lowering distortion (Supplementary Table 7), thus frustrating FiM spin ordering.
Based on the previous outcomes, we hypothesize that the following three conditions may be necessary for the stabilization of VO-induced FiM spin ordering in magnetic AMO3 oxide perovskites: (1) lack of inversion symmetry leading to polar order and structural deformation ease, (2) moderate axial ratio structures allowing for the reduction of second-nearest transition metal ions, and (3) the existence of highly magnetized oxygen ions.
Discussion
The results obtained in BCO thin films raise the natural question: are there any other multiferroic materials in which similar VO-induced magnetic phenomena may occur? To answer this question we investigated the special case of BiFeO3 (BFO) and other Bi-based multiferroic (BiMnO3 and BiCrO3) thin films. Figure 5 encloses the energy, magnetic, and vibrational properties of non-stoichiometric BFO thin films. In order to be consistent with the notation employed heretofore, we label the usual rhombohedral-like monoclinic phase of BFO as \({\mathcal{M}}\) and the large axial ratio tetragonal-like phase as \({\mathcal{T}}\) (in spite of the fact that the space groups corresponding to those structures are different in BFO and BCO)50,55. Our zero-temperature calculations (Fig. 5a) predict a ground-state \({\mathcal{M}}\)-G(Ap) phase at ain ≥ 4.01 Å, followed by \({\mathcal{T}}\)-G(Eq) at 3.91 ≤ ain ≤ 4.01 Å, and \({\mathcal{T}}\)-C(Eq) at ain ≤ 3.91 Å. Indeed, when the likely existence of VO is explicitly considered in the simulations a broad ain region appears in which AFM-G spin ordering is stable for the \({\mathcal{T}}\) phase. (We recall that AFM-G spin ordering is metastable in the corresponding stoichiometric system48,50). The causes of the stabilization of the \({\mathcal{T}}\)-G(Eq) phase in BFO are very similar to those explained previously for \({\mathcal{T}}\)-G(Ap) in BCO thin films (Fig. 3). As mentioned earlier, these results may shed new light on the origins of some unresolved discrepancies between theory and experiments regarding the identification of antiferromagnetic ordering in \({\mathcal{T}}\) BFO thin films48,49,50,51,52.
a Zero-temperature energy of competing phases expressed as a function of in-plane lattice parameter. b Structural representation of non-stoichiometric \({\mathcal{M}}\) phases considering different magnetic orderings. c Γ-point phonon spectrum calculated for non-stoichiometric monoclinic BFO thin films. For equivalent eigenmodes, the phonon frequencies of the \({\mathcal{M}}\)-G(Ap) phase in general are higher than of the \({\mathcal{M}}\)-FiM(Ap) phase hence the vibrational free energy of the latter phase is lower. d Spin-up (dark green) and spin-down (light green) electronic charge densities calculated for stoichiometric BCO and BFO thin films considering two perpendicular planes; highly magnetized oxygen atoms are indicated with red squares.
By creating an apical VO in the stoichiometric \({\mathcal{M}}\)-C phase, we found a thus far neglected ferrimagnetic phase in monoclinic BFO thin films, \({\mathcal{M}}\)-FiM(Ap). In this case the net magnetization per formula unit amounts also to 0.5 μB for BiFeO2.75. Nevertheless, the \({\mathcal{M}}\)-FiM(Ap) phase is metastable at zero temperature due to a small energy difference of 25–30 meV per formula unit with respect to the ground-state phase \({\mathcal{M}}\)-G(Ap) (Fig. 5a). The atomic structure of the metastable \({\mathcal{M}}\)-FiM(Ap) and ground-state \({\mathcal{M}}\)-G(Ap) phases are highly distorted and surprisingly very similar (Fig. 5b and Supplementary Tables 8–9). Analysis of the Γ-point phonon modes (Fig. 5c), however, indicates that the \({\mathcal{M}}\)-FiM(Ap) phase is vibrationally softer than \({\mathcal{M}}\)-G(Ap). Consequently, the \({\mathcal{M}}\)-FiM(Ap) phase could be entropically stabilized over the \({\mathcal{M}}\)-G(Ap) under increasing temperature due to the fact that the zero-temperature energies of the two states are similar and the vibrational free energy of the former phase appears to be more favourable9,31. Our results, therefore, can be interpreted as theoretical evidence showing that the observation of “ferromagnetic” behaviour in BFO thin films24 may be caused by the presence of oxygen vacancies (as it was suggested previously by other authors25).
Figure 5d shows spin-up spin-down surface charge densities calculated for non-stoichiometric \({\mathcal{M}}\) and \({\mathcal{T}}\) BFO thin films considering in-plane and out-of-plane orientations. The reason why we could only find one FiM solution when generating neutral VO’s in stoichiometric \({\mathcal{M}}\) phases is now clear: all oxygen atoms in the \({\mathcal{M}}\)-G phase are non-magnetic whereas all apical O in the \({\mathcal{M}}\)-C phase are highly magnetized (Fig. 5d). The stoichiometric \({\mathcal{T}}\)-C phase also displays highly magnetized oxygen atoms in apical positions however, as we have explained before, it seems difficult to observe FiM spin ordering in c/a >1 geometries. Furthermore, we repeated the same type of calculations and analysis for monoclinic-like BiMnO3 and BiCrO3 thin films, which fulfill conditions (1) and (2) above. We found that the appearance of FiM spin ordering is also correlated with the presence of highly magnetized oxygen atoms (Supplementary Fig. 8). In particular, a \({\mathcal{M}}\)-FiM(Ap) phase displaying a net magnetization of ≈0.2 μB is found in BiCrO2.75. These last results confirm that the simple rules provided in this work for understanding the appearance of FiM phases in non-stoichiometric multiferroic thin films may be general and robust.
In summary, by using first-principles calculations we have disclosed a number of previously overlooked electro-structural mechanisms induced by the presence of oxygen vacancies that facilitate the stabilization of unexpected magnetic states in multiferroic thin films. In particular, AFM-G and FiM spin orderings may naturally appear in large axial ratio and monoclinic phases under specific lattice strain conditions. Our theoretical results may clarify the origins of some long-standing controversies in BFO, the paradigm of single-phase multiferroics and one of the most intensively studied multifunctional materials. We provide general and simple rules to fundamentally understand and systematically predict FiM phases in non-stoichiometric multiferroics, thus offering new approaches for the rational engineering of bettered multifunctional materials. The present work shows that oxygen vacancies should be considered as a design opportunity to create new funcionalities in multiferroic thin films, especially to what concerns magnetism.
Methods
Density functional theory calculations
First-principles spin-polarized calculations based on density functional theory (DFT) are performed with the generalized gradient approximation proposed by Perdew, Burke and Ernzerhof (GGA-PBE) as implemented in the VASP package56,57. We employ the “Hubbard-U” scheme derived by Dudarev et al. to deal with the 3d electrons in Co (Fe) atoms and, as done in previous works, a U value of 6 (4) eV is adopted9,30,31. We use the “projected augmented wave” method58 considering the following electronic states as valence: Co’s 4s13d8, Fe’s 3p64s13d7, Mn’s 4s13d6, Cr’s 4s13d5, Bi’s 6s25d106p3, and O’s 2s22p4. The energy cut-off is truncated at 650 eV and we employ a Γ-centered k-point grid of 4 × 6 × 6 for a \(2\times \sqrt{2}\times \sqrt{2}\) supercell containing 20 atoms (that is, four formula units) 55. The reported net magnetizations correspond to sums over the local magnetic moments, which are accurate to within 0.01 μB per formula unit.
Thin film geometry relaxations are carried out by using a conjugated gradient algorithm that allows to change the simulation-cell volume and atomic positions while constraining the length and orientation of the two in-plane lattice vectors (that is, ∣a∣ = ∣b∣ and γ = 90°). Periodic boundary conditions are applied along the three lattice-vector directions, thus the influence of surface and interface effects are systematically disregarded in our simulations. This type of calculations are known as “strained-bulk” geometry relaxations and typically are considered to be a good approximation for thin films with thicknesses of at least few nanometers48,55 (i.e., for which surface and interface effects can be safely neglected). In our simulations, we assume systems that are elastically coupled to a substrate hence we disregard the existence of possible stress relaxation mechanisms in thin films. The geometry relaxations are stopped when the forces on the ions are smaller than 0.01 eV/Å. We have checked the vibrational stability of every phase by estimating the lattice phonons at the Γ-point with the small-displacement method59 and considering central differences for the calculation of atomic forces derivatives. Full phonon spectra calculations have not been pursued in this work due to the huge computational expense associated with large and low-symmetry off-stoichiometric unit cells.
In this study, most non-stoichiometric configurations are generated by removing one oxygen atom from an apical or equatorial position in a 20-atoms simulation cell, thus rendering the chemical composition BiCoO2.75. Apical and equatorial VO positions are investigated systematically by generating all inequivalent configurations in all competing phases and considering all possible magnetic spin orderings (FM, AFM-G, AFM-G, and AFM-A31). Additional simulations in a 40-atoms simulation cell have been performed to assess the influence of size effects and VO concentration in our main conclusions (Supplementary Fig. 10 and Supplementary Discussion). The results of our tests indicate that our findings obtained for BiCoO2.75 are general and qualitatively correct also for BiCoO2.875 thin films. The influence of possible kinetic effects on the distribution of oxygen vacancies has been totally neglected in this study.
The results presented in the main text apply to neutral oxygen vacancies since we have found that the creation of neutral VO’s in energetically more favourable than that of charged oxygen vacancies (Supplementary Fig. 9, Supplementary Table 10, and Supplementary Discussion). In particular, to estimate the ranking of VO formation energies as a function of charge, q, we have employed the following well-established formula60:
where \(E[{V}_{O}^{q}]\) is the energy of the non-stoichiometric system containing the oxygen vacancy, Estoi the energy of the corresponding stoichiometric system, \({E}_{corr}^{q}\) a finite-size supercell correction, nO the number of created VO’s (typically equal to 1 in our calculations), μO the chemical potential of oxygen atoms, ϵF the Fermi energy in the non-stoichiometric system, ϵv the top energy in the valence band of the non-stoichiometric system, and ΔV a term used for aligning the electrostatic potentials of the stoichiometric and defective supercells. In order to calculate \({E}_{corr}^{q}\) and ΔV, we have followed the methods explained in work60. According to our estimations, neutral VO’s (q = 0) are energetically more favourable than charged vacancies (q = +2 e) by about ~1 eV per formula unit (Supplementary Fig. 9, Supplementary Table 10, and Supplementary Discussion).
We have performed several tests to assess the influence of the adopted DFT exchange-correlation functional and U value on our theoretical predictions (Supplementary Figs. 4–6, Supplementary Tables 11–13, and Supplementary Discussion). Specifically, we repeated most calculations by considering the PBEsol functional61 and 2 ≤ U ≤ 6 eV values. It is found that the main conclusions presented in the main text are not affected qualitatively by the choice of the U parameter or exchange-correlation functional. At the quantitative level, the ain parameters at which the phase transitions occur and the energy ranking of the competing phases change in some cases (Supplementary Figs. 4–6, Supplementary Tables 11–13, and Supplementary Discussion). For a detailed discussion on these technical aspects, see Supplementary Discussion. It is also worth noting that lattice thermal excitations have been neglected systematically throughout this study35. Nevertheless, we have performed a test to assess the influence of this approximation in our main findings (Supplementary Fig. 11 and Supplementary Discussion) and found that the conclusions reported in this study appear to be robust and physically meaningful also for T ≠ 0 conditions.
Regarding the estimation of the electric polarization, P, we started by employing the Berry phase formalism62. However, the presence of oxygen vacancies induces a notable reduction in the energy band gap of the system (Supplementary Fig. 2) that in some cases frustrates the determination of the corresponding Berry phase (due to the appearance of intermediate metallic phases). In order to overcome such a limitation, we opted for calculating the electric polarization perturbatively. Specifically, we estimate P with the formula55:
where Ω is the volume of the cell, κ runs over all the atoms, α, β = x, y, z represent Cartesian directions, uκ is the displacement vector of the κ-th atom as referred to a non-polar reference phase, and \({{\boldsymbol{Z}}}_{\kappa }^{* }\) the Born effective charge tensor calculated in a non-polar reference state (\({{\boldsymbol{Z}}}_{\kappa }^{* }\) has been calculated for each defective system). In the stoichiometric phases we do not find the technical limitations just explained for non-stoichiometric systems, hence in that case we have been able to compare the P values obtained with the Berry phase approach (exact) and Eq. (2) (approximate). According to our estimations, the electric polarizations calculated perturbatively are accurate to within ~10% of, and systematically larger than, the P values calculated with the Berry phase method. It is reasonable to assume then a similar level of accuracy for the P values estimated for non-stoichiometric thin films that are reported in Fig. 2.
Heisenberg model Monte Carlo simulations
To simulate the effects of thermal excitations on the magnetic ordering of \({\mathcal{T}}\), \({\mathcal{O}}\), and \({\mathcal{M}}\) BCO thin films, both stoichiometric and non-stoichiometric, we construct several spin Heisenberg models of the form \(\hat{H}=\frac{1}{2}{\sum }_{ij}{J}_{ij}{S}_{i}{S}_{j}\), in which the value of the involved exchange constants are obtained from zero-temperature DFT calculations (following analogous methods to those explained in detail in previous works of ours9,30,31).
For the \({\mathcal{T}}\) phase (c > a), we consider spin couplings between first nearest in-plane and out-of-plane neighbours (Ja and Jc, respectively) and between second nearest neighbours (Jac). The value of such magnetic exchange constants are obtained with the formulas9,30,31:
where E represents energy and superscripts “FM”, “G”, “A”, and “C” ferromagnetic, antiferromagnetic G-type, antiferromagnetic A-type, and antiferromagnetic C-type spin arrangements, respectively. For the \({\mathcal{O}}\) and \({\mathcal{M}}\) phases (c ≈ a), we consider only spin couplings between first nearest neighbours and obtain their value with the expression9,30,31:
The value of some characteristic magnetic exchange constants estimated in the manner just described are reported in Supplementary Table 14.
We use the Heisenberg models described above to perform Monte Carlo (MC) simulations in a periodically-repeated simulation box of 20 × 20 × 20 spins; thermal averages are computed from runs of 50,000 MC sweeps after equilibration. These simulations allow us to monitor the T-dependence of the magnetic ordering through the computation of the AFM-C (in the \({\mathcal{T}}\) phase) and AFM-G (in the \({\mathcal{O}}\) and \({\mathcal{M}}\) phases) order parameters, namely, \({S}^{{\rm{C}}}\equiv \frac{1}{N}{\sum }_{i}{(-1)}^{{n}_{ix}+{n}_{iy}}{S}_{iz}\) and \({S}^{{\rm{G}}}\equiv \frac{1}{N}{\sum }_{i}{(-1)}^{{n}_{ix}+{n}_{iy}+{n}_{iz}}{S}_{iz}\). Here, nix, niy, and niz are the three integers locating the i-th lattice cell, and N is the total number of spins in the simulation box. For the calculation of SC and SG, we considered only the z component of the spins because a small symmetry-breaking magnetic anisotropy was introduced in the Hamiltonian in order to facilitate the numerical analysis9,30,31.
Data availability
The data that support the findings of this study are available from the corresponding author (C.C.) upon reasonable request.
References
Ramesh, R. & Spaldin, N. A. Multiferroics: progress and prospects in thin films. Nat. Mater. 6, 21 (2007).
Spaldin, N. A. & Ramesh, R. Advances in magnetoelectric multiferroics. Nat. Mater. 18, 203 (2019).
Heron, J. T. et al. Deterministic switching of ferromagnetism at room temperature using an electric field. Nature 516, 370 (2014).
Allibe, J. et al. Room temperature electrical manipulation of giant magnetoresistance in spin valves exchange-biased with BiFeO3. Nano Lett. 12, 1141 (2012).
Domann, J. P. & Carman, G. P. Strain powered antennas. J. Appl. Phys. 121, 044905 (2017).
Nan, T. et al. Acoustically actuated ultra-compact NEMS magnetoelectric antennas. Nat. Commun. 8, 296 (2017).
Huang, W. et al. Multiferroic Bi2 FeCrO6 based p-i-n heterojunction photovoltaic devices. J. Mater. Chem. A 5, 10355 (2017).
Stern-Taulats, E. et al. Multicaloric materials and effects. MRS Bull. 43, 295 (2018).
Cazorla, C. & Íñiguez, J. Giant direct and inverse electrocaloric effects in multiferroic thin films. Phys. Rev. B 98, 174105 (2018).
Cazorla, C. Novel mechanocaloric materials for solid-state cooling applications. Appl. Phys. Rev. 6, 041316 (2019).
Bertinshaw, J. et al. Direct evidence for the spin cycloid in strained nanoscale bismuth ferrite thin films. Nat. Commun. 7, 12664 (2016).
Belik, A. A., Azuma, M., Saito, T., Shimakawa, Y. & Takano, M. Crystallographic features and tetragonal phase stability of PbVO3, a new member of PbTiO3 family. Chem. Mater. 17, 269 (2005).
Belik, A. A. et al. Neutron powder diffraction study on the crystal and magnetic structures of BiCoO3. Chem. Mater. 18, 798 (2006).
Oka, K. et al. Pressure-induced spin-state transition in BiCoO3. J. Am. Chem. Soc. 132, 9438 (2010).
Wang, X. et al. E-field control of the RKKY interaction in FeCoB/Ru/FeCoB/PMN-PT (011) multiferroic heterostructures. Adv. Mater. 30, 1803612 (2018).
Das, R. & Mandal, K. Magnetic, ferroelectric, and magnetoelectric properties of Ba-doped BiFeO3. J. Magn. Magn. Mater. 324, 1913 (2012).
Das, R., Sharma, S. & Mandal, K. Aliovalent Ba2+ doping: A way to reduce oxygen vancancy in multiferroic BiFeO3. J. Magn. Magn. Mater. 401, 129 (2016).
Sakai, H. et al. Displacement-type ferroelectricity with off-center magnetic ions in perovskite Sr1−x Bax MnO3. Phys. Rev. Lett. 107, 137601 (2011).
Hojo, H. et al. Development of bismuth ferrite as a piezoelectric and multiferroic material by cobalt substitution. Adv. Mater. 30, 1705665 (2018).
Martin, L. W. et al. Multiferroics and magnetoelectrics: thin films and nanostructures. J. Phys.: Condens. Matter 20, 434220 (2008).
Zhang, L. et al. Giant polarization in super-tetragonal thin films through interphase strain. Science 361, 494 (2018).
Varga, T. et al. Coupled lattice polarization and ferromagnetism in multiferroic NiTiO3 thin films. ACS Appl. Mater. Interfaces 9, 21879 (2017).
Töpfer, J. & Goodenough, J. G. LaMnO3+δ revisited. J. Solid State Chem. 130, 117 (1997).
Wang, J. et al. Epitaxial BiFeO3 multiferroic thin film heterostructures. Science 299, 1719 (2003).
Eerenstein, W. et al. Comment on “Epitaxial BiFeO3 multiferroic thin film heterostructures”. Science 307, 1203 (2005).
Goodenough, J. B. Theory of the role of covalence in the perovskite-type manganites (La,M (II))MnO3. Phys. Rev. 100, 564 (1955).
Hill, N. A. & Filipetti, A. Why are there any magnetic ferroelectrics? J. Magn. Magn. Mater. 242, 976 (2002).
Spaldin, N.A. Magnetic Materials: Fundamentals and Applications. 2nd edn (Cambridge University Press), https://doi.org/10.1017/CBO9780511781599.021 (2012).
Bhattacharya, A. & May, S. Magnetic oxide heterostructures. Annu. Rev. Mater. Res. 44, 65 (2014).
Cazorla, C. & Íñiguez, J. Insights into the phase diagram of bismuth ferrite from quasiharmonic free-energy calculations. Phys. Rev. B 88, 214430 (2013).
Cazorla, C., Diéguez, O. & Íñiguez, J. Multiple structural transitions driven by spin-phonon couplings in a perovskite oxide. Sci. Adv. 3, e1700288 (2017).
Singh, D. J. Electronic structure and bond competition in the polar magnet PbVO3. Phys. Rev. B 73, 094102 (2006).
Solovyev, I. V. Magnetic structure of the noncentrosymmetric perovskites PbVO3 and BiCoO3: theoretical analysis. Phys. Rev. B 85, 054420 (2012).
Niu, W. et al. Direct demonstration of the emergent magnetism resulting from the multivalence Mn in a LaMnO3 epitaxial thin film system. Adv. Electron. Mater. 4, 1800055 (2018).
Cazorla, C. Lattice effects on the formation of oxygen vacancies in perovskite thin films. Phys. Rev. Appl. 7, 044025 (2017).
Herklotz, A. et al. Strain coupling of oxygen non-stoichiometry in perovskite thin films. J. Phys.: Condens. Matter 29, 493001 (2017).
Mangalam, R. V. K., Ray, N., Waghmare, U. V., Sundaresana, A. & Rao, C. N. R. Multiferroic properties of nanocrystalline BaTiO3. Solid State Commun. 149, 1 (2009).
Shimada, T. et al. Multiferroic vacancies at ferroelectric PbTiO3 surfaces. Phys. Rev. Lett. 115, 107202 (2015).
Shimada, T. et al. Multiferroic nature of intrinsic point defects in BiFeO3 : a hybrid Hartree-Fock density functional study. Phys. Rev. B 93, 174107 (2016).
Li, W. et al. Oxygen-vacancy-induced antiferromagnetism to ferromagnetism transformation in Eu0.5 Ba0.5 TiO3−δ multiferroic thin films. Sci. Rep. 3, 2618 (2013).
Ederer, C. & Spaldin, N. A. Influence of strain and oxygen vacancies on the magnetoelectric properties of multiferroic bismuth ferrite. Phys. Rev. B 71, 224103 (2005).
Paudel, T. R., Jaswal, S. S. & Tsymbal, E. Y. Intrinsic defects in multiferroic BiFeO3 and their effect on magnetism. Phys. Rev. B 85, 104409 (2012).
Chen, X.-Y. et al. Tuning the polarization and magnetism in BiCoO3 by strain an oxygen vacancy effect: a first-principle study. J. Appl. Phys. 111, 013901 (2012).
Diéguez, O. & Íñiguez, J. First-principles investigation of morphotropic transitions and phase-change functional responses in BiFeO3 -BiCoO3 multiferroic solid solutions. Phys. Rev. Lett. 107, 057601 (2011).
Jaffe, B., Roth, R. S. & Marzullo, S. Piezoelectric properties of lead zirconate-lead titanate solid-solution ceramics. J. Appl. Phys. 25, 809 (1954).
Bellaiche, L., García, A. & Vanderbilt, D. Finite-temperature properties of Pb(Zr1−x Tix)O3 alloys from first principles. Phys. Rev. Lett. 84, 5427 (2000).
Shimizu, K., Hojo, H., Ikuhara, Y. & Azuma, M. Enhanced piezoelectric response due to polarization rotation in cobalt-substituted BiFeO3 epitaxial thin films. Adv. Mater. 28, 8639 (2016).
Hatt, A. J., Spaldin, N. A. & Ederer, C. Strain-induced isosymmetric phase transition in BiFeO3. Phys. Rev. B 81, 054109 (2010).
Diéguez, O., González-Vázquez, O. E., Wojdel, J. C. & Íñiguez, J. First-principles predictions of low-energy phases of multiferroic BiFeO3. Phys. Rev. B 83, 094105 (2011).
Heo, Y. et al. Impact of isovalent and aliovalent doping on mechanical properties of mixed phase BiFeO3. ACS Nano 11, 2805 (2017).
Béa, H. et al. Evidence for room-temperature multiferroicity in a compound with a giant axial ratio. Phys. Rev. Lett. 102, 217603 (2009).
MacDougall, G. J. et al. Antiferromagnetic transitions in tetragonal-like BiFeO3. Phys. Rev. B 85, 100406 (2012).
Rivero, P. & Cazorla, C. Revisiting the zero-temperature phase diagram of stoichiometric SrCoO3 with first-principles methods. Phys. Chem. Chem. Phys. 18, 30686 (2016).
Halcrow, M. A. Jahn-Teller distortions in transition metal compounds, and their importance in functional molecular and inorganic materials. Chem. Soc. Rev. 42, 1784 (2013).
Cazorla, C. & Stengel, M. Electrostatic engineering of strained ferroelectric perovskites from first principles. Phys. Rev. B 92, 214108 (2015).
Kresse, G. & Furthmüller, J. Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set. Phys. Rev. B 54, 11169 (1996).
Perdew, J. P., Burke, K. & Ernzerhof, M. Generalized gradient approximation made simple. Phys. Rev. Lett. 77, 3865 (1996).
Blöchl, P. E. Projector augmented-wave method. Phys. Rev. B 50, 17953 (1994).
Cazorla, C. & Boronat, J. Simulation and understanding of atomic and molecular quantum crystals. Rev. Mod. Phys. 89, 035003 (2017).
Komsa, H.-P., Rantala, P. P. & Pasquarello, A. Finite-size supercell correction schemes for charged defect calculations. Phys. Rev. B 86, 045112 (2012).
Perdew, J. P. et al. Restoring the density-gradient expansion for exchange in solids and surfaces. Phys. Rev. Lett. 100, 136406 (2008).
King-Smith, R. D. & Vanderbilt, D. Theory of polarization of crystalline solids. Phys. Rev. B 47, 1651(R) (1993).
Acknowledgements
Computational resources and technical assistance were provided by the Australian Government and the Government of Western Australia through the National Computational Infrastructure (NCI) and Magnus under the National Computational Merit Allocation Scheme and The Pawsey Supercomputing Centre.
Author information
Authors and Affiliations
Contributions
C.C. conceived the study and planned the research. C.M. and C.C. performed the theoretical calculations. Results were discussed by C.M., D.C., and C.C. The paper was written by C.M., D.C., and C.C.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary information
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Menéndez, C., Chu, D. & Cazorla, C. Oxygen-vacancy induced magnetic phase transitions in multiferroic thin films. npj Comput Mater 6, 76 (2020). https://doi.org/10.1038/s41524-020-0344-3
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41524-020-0344-3
This article is cited by
-
Giant multiphononic effects in a perovskite oxide
npj Computational Materials (2023)
-
Defect graph neural networks for materials discovery in high-temperature clean-energy applications
Nature Computational Science (2023)
-
Mechanisms for point defect-induced functionality in complex perovskite oxides
Applied Physics A (2022)
-
Ab initio description of oxygen vacancies in epitaxially strained \(\hbox {SrTiO}_{{3}}\) at finite temperatures
Scientific Reports (2021)