Revisiting properties of ferroelectric and multiferroic thin films under tensile strain from first principles

Lone-pair activated ferroelectricity and stable charged domain wall in Bi monolayer

Ferroelectricity has been predicted in two-dimensional Group-Va elemental materials and confirmed in high-quality Bi monolayers by a recent experiment. The origin of such elemental ferroelectricity is related to the spontaneous lattice distortion with atomic layer buckling. A surprising observation in experiment is the abundance of charged 180° head-to-head/tail-to-tail domain walls, distinct from conventional ferroelectrics, where the naturally occurring ferroelectric domain walls are mostly charge neutral. Here, we clarify the origin of this phenomenon. We find that distinct from conventional ferroelectrics, in such single-element ferroelectric monolayers, it is the strain energy rather than the electrostatic energy that dominates the energetics. This leads to intrinsically stable 180° charged domain walls. The orbital interaction and the lone-pair activation mechanism play a key role in this picture. We further predict and confirm experimentally that the most stable domain wall type changes from charged to neutral ones under small applied strain. Our work reveals a mechanism to generate polarization and stabilize intrinsic charged domain walls, which will shed light on potential applications of ferroelectronics based on charged domain walls.

Ultrafast Optically Induced Perturbation of Oxygen Octahedral Rotations in Multiferroic BiFeO3 Thin Films

The functional properties of complex oxides, including magnetism and ferroelectricity, are closely linked to subtle structural distortions. Ultrafast optical excitations provide the means to manipulate structural features and ultimately to affect the functional properties of complex oxides with picosecond-scale precision. We report that the lattice expansion of multiferroic BiFeO3 following above-bandgap optical excitation leads to distortion of the oxygen octahedral rotation (OOR) pattern. The continuous coupling between OOR and strain was probed using time-resolved X-ray free-electron laser diffraction with femtosecond time resolution. Density functional theory calculations predict a relationship between the OOR and the elastic strain consistent with the experiment, demonstrating a route to employing this approach in a wider range of systems. Ultrafast control of the functional properties of BiFeO3 thin films is enabled by this approach because the OOR phenomena are related to ferroelectricity, and via the Fe-O-Fe bond angles, the superexchange interaction between Fe atoms.

Unexpected Competition between Ferroelectricity and Rashba Effects in Epitaxially Strained SrTiO_{3}

The Rashba parameter α_{R} is usually assumed to scale linearly with the amplitude of polar displacements by construction of the spin-orbit interaction. On the basis of first-principles simulations, ferroelectric phases of SrTiO_{3} reached under epitaxial compressive strain are characterized by (i) large Rashba effects at the bottom of the conduction band near the paraelectric-ferroelectric boundary and (ii) an unexpected suppression of the phenomena when the amplitude of polar displacements increases. This peculiar behavior is ascribed to the inverse dependence of the Rashba parameter with the crystal field Δ_{CF} induced by the polar displacements that alleviates the degeneracy of Ti t_{2g} states and annihilates the Rashba effects. Although α_{R} has intrinsically a linear dependance on polar displacements, the latter becomes antagonist to Rashba phenomena at large polar mode amplitude. Thus, the Rashba coefficient may be bound to an upper value.

Tunable and parabolic piezoelectricity in hafnia under epitaxial strain

Piezoelectrics are a class of functional materials that have been extensively used for application in modern electro-mechanical and mechatronics technologies. The sign of longitudinal piezoelectric coefficients is typically positive but recently a few ferroelectrics, such as ferroelectric polymer poly(vinylidene fluoride) and van der Waals ferroelectric CuInP2S6, were experimentally found to have negative piezoelectricity. Here, using first-principles calculation and measurements, we show that the sign of the longitudinal linear piezoelectric coefficient of HfO2 can be tuned from positive to negative via epitaxial strain. Nonlinear and even parabolic piezoelectric behaviors are further found at tensile epitaxial strain. This parabolic piezoelectric behavior implies that the polarization decreases when increasing the magnitude of either compressive or tensile longitudinal strain, or, equivalently, that the strain increases when increasing the magnitude of electric field being either parallel or antiparallel to the direction of polarization. The unusual piezoelectric effects are from the chemical coordination of the active oxygen atoms. These striking piezoelectric features of positive and negative sign, as well as linear and parabolical behaviors, expand the current knowledge in piezoelectricity and broaden the potential of piezoelectric applications towards electro-mechanical and communications technology.

Strain-Induced Structural Phase Transitions in Epitaxial (001) BiCoO3 Films: A First-Principles Study

We have simulated BiCoO3 films epitaxially grown along (001) direction with density functional theory computations. Leading candidates for the lowest-energy phases have been identified. The tensile strains induce magnetic phase transition in the ground state (P4mm symmetry) from a C-type antiferromagnetic order to a G-type order for the in-plane lattice parameter above 3.922 Å. The G-type antiferromagnetic order will be maintained with larger tensile strains; however, a continuous structural phase transition will occur, combining the ferroelectric and antiferrodistortive modes. In particular, the larger tensile strain allows an isostructural transition, the so-called Cowley's ''Type Zero'' phase transitions, from Cc-(I) to Cc-(II), with a slight volume collapse. The orientation of ferroelectric polarization changes from the out-of-plane direction in the P4mm to the in-plane direction in the Pmc21 state under epitaxial tensile strain; meanwhile, the magnetic ordering temperature TN can be strikingly affected by the variation of misfit strain.

Polar Magnetism Above 600 K with High Adaptability in Perovskite Oxides

High magnetic order temperature, sustainable polar insulating state, and tolerance to device integrations are substantial advantages for applications in next-generation spintronics. However, engineering such functionality in a single-phase system remains a challenge owing to the contradicted chemical and electronic requirements for polar nature and magnetism, especially with an ordering state highly above room temperature. Perovskite-related oxides with unique flexibility allow electron-unpaired subsystems to merge into the polar lattice to induce magnetic interactions, combined with their inherent asymmetry, thereby promising polar magnet design. Herein, by atomic-level composition assembly, a family of Ti/Fe co-occupied perovskite oxide films Pb(Ti_1-x,Fe_x)O₃ (PFT(x)) with a Ruddlesden-Popper superstructure are successfully synthesized on several different substrates, demonstrating exceptional adaptability to different integration conditions. Furthermore, second-harmonic generation measurements convince the symmetry-breaking polar character. Notably, a ferromagnetic ground state up to 600 K and a steady insulating state far beyond room temperature were achieved simultaneously in these films. This strategy of constructing layered modular superlattices in perovskite oxides could be extended to other strongly correlated systems for triggering nontrivial quantum physical phenomena.

Physics-Guided Descriptors for Prediction of Structural Polymorphs

We develop a method combining machine learning (ML) and density functional theory (DFT) to predict low-energy polymorphs by introducing physics-guided descriptors based on structural distortion modes. We systematically generate crystal structures utilizing the distortion modes and compute their energies with single-point DFT calculations. We then train a ML model to identify low-energy configurations on the material's high-dimensional potential energy surface. Here, we use BiFeO₃ as a case study and explore its phase space by tuning the amplitudes of linear combinations of a finite set of distinct distortion modes. Our procedure is validated by rediscovering several known metastable phases of BiFeO₃ with complex crystal structures, and its efficiency is proved by identifying 21 new low-energy polymorphs. This approach proposes a new avenue toward accelerating the prediction of low-energy polymorphs in solid-state materials.

Comparative density functional studies of pristine and doped bismuth ferrite polymorphs by GGA+U and meta-GGA SCAN+U

We analyzed the effect of nonmagnetic dopants Al, Ga, Sc, and In at the Fe-site on the phase stability, structural, and electronic properties of different bismuth ferrite (BiFeO₃) polymorphs in the framework of density functional theory with the Hubbard U correction (DFT+U). We started our consideration from the determination of the magnitude of the U parameter. First, we calculated the structural, electronic, and magnetic properties of the rhombohedral R3c-G phase of BiFeO₃ within the generalized gradient approximation (GGA) and strongly constrained and appropriately normed (SCAN) meta-GGA for different values of the U. Next, we compared these results with those obtained within the parameter-free hybrid functional. After determining the optimal values of the Hubbard U parameter we analyzed the total energies between the selected BiFeO₃ polymorphs without and with dopants within both GGA+U and SCAN+U. For all dopants the concentration was 12.5% which was close to their solubility limit in BiFeO₃ under ambient conditions. We found that none of these dopants led to the structural phase transition. However, DFT+U calculations revealed that the doping of BiFeO₃ with Al and Ga reduced the energy barrier between R3c-G and Cm-C phases whereas for Sc and In the energy difference between both phases increased. For the orthorhombic phases the considered dopants do not affect the energy barrier between them and the rhombohedral phase. In addition, the ferroelectric polarization does not change after replacing the Fe atom by the dopant for the all considered BiFeO₃ polymorphs.

Dzyaloshinskii-Moriya-like interaction in ferroelectrics and antiferroelectrics

The Dzyaloshinskii-Moriya interaction (DMI) between two magnetic moments m_i and m_j is of the form [Formula: see text]. It originates from spin-orbit coupling, and is at the heart of fascinating phenomena involving non-collinear magnetism, such as magnetic topological defects (for example, skyrmions) as well as spin-orbit torques and magnetically driven ferroelectricity, that are of significant fundamental and technological interest. In sharp contrast, its electric counterpart, which is an electric DMI characterized by its [Formula: see text] strength and describing an interaction between two polar displacements u_i and u_j, has rarely been considered, despite the striking possibility that it could also generate new features associated with non-collinear patterns of electric dipoles. Here we report first-principles simulations combined with group theoretical symmetry analysis which not only demonstrate that electric DMI does exist and has a one-to-one correspondence with its magnetic analogue, but also reveals a physical source for it. These findings can be used to explain and/or design phenomena of possible technological importance in ferroelectrics and multiferroics.

Strain-induced structural phase transition, electric polarization and unusual electric properties in photovoltaic materials CsMI3 (M = Pb, Sn)

The structural phase transition, ferroelectric polarization, and electric properties have been investigated for photovoltaic films CsMI3 (M = Pb, Sn) epitaxially grown along (001) direction based on the density functional theory. The calculated results indicate that the phase diagrams of two epitaxial CsPbI3 and CsSnI3 films are almost identical, except critical transition strains varying slightly. The epitaxial tensile strains induce two ferroelectric phases Pmc21, and Pmn21, while the compressive strains drive two paraelectric phases P212121, P21212. The larger compressive strain enhances the ferroelectric instability in these two films, eventually rendering them another ferroelectric state Pc. Whether CsPbI3 or CsSnI3, the total polarization of Pmn21 phase comes from the main contribution of B-position cations (Pb or Sn), whereas, for Pmc21 phase, the main contributor is the I ion. Moreover, the epitaxial strain effects on antiferrodistortive vector, polarization and band gap of CsMI3 (M = Pb, Sn) are further discussed. Unusual electronic properties under epitaxial strains are also revealed and interpreted.

Magneto-electric multiferroics: designing new materials from first-principles calculations

In parallel with the revival of interest for magneto-electric multiferroic materials in the beginning of the century, first-principles simulations have grown incredibly in efficiency during the last two decades. Density functional theory calculations, in particular, have so become a must-have tool for physicists and chemists in the multiferroic community. While these calculations were originally used to support and explain experimental behaviour, their interest has progressively moved to the design of novel magneto-electric multiferroic materials. In this article, we mainly focus on oxide perovskites, an important class of multifunctional material, and review some significant advances to which contributed first-principles calculations. We also briefly introduce the various theoretical developments that were at the core of all these advances.

Frustrated Dipole Order Induces Noncollinear Proper Ferrielectricity in Two Dimensions

Within Landau theory, magnetism and polarity are homotopic, displaying a one-to-one correspondence between most physical characteristics. However, despite widely reported noncollinear magnetism, spontaneous noncollinear electric dipole order as a ground state is rare. Here, a dioxydihalides family is predicted to display noncollinear ferrielectricity, induced by competing ferroelectric and antiferroelectric soft modes. This intrinsic of dipoles generates unique physical properties, such as Z_{2}×Z_{2} topological domains, atomic-scale dipole vortices, and negative piezoelectricity.

A group-theoretical approach to enumerating magnetoelectric and multiferroic couplings in perovskites

A group-theoretical approach is used to enumerate the possible couplings between magnetism and ferroelectric polarization in the parent Pm{\overline 3}m perovskite structure. It is shown that third-order magnetoelectric coupling terms must always involve magnetic ordering at the A and B sites which either transforms both as R-point or both as X-point time-odd irreducible representations (irreps). For fourth-order couplings it is demonstrated that this criterion may be relaxed allowing couplings involving irreps at X-, M- and R-points which collectively conserve crystal momentum, producing a magnetoelectric effect arising from only B-site magnetic order. In this case, exactly two of the three irreps entering the order parameter must be time-odd irreps and either one or all must be odd with respect to inversion symmetry. It is possible to show that the time-even irreps in this triad must transform as one of: X1+, M3,5- or R5+, corresponding to A-site cation order, A-site antipolar displacements or anion rocksalt ordering, respectively. This greatly reduces the search space for type-II multiferroic perovskites. Similar arguments are used to demonstrate how weak ferromagnetism may be engineered and a variety of schemes are proposed for coupling this to ferroelectric polarization. The approach is illustrated with density functional theory calculations on magnetoelectric couplings and, by considering the literature, suggestions are given of which avenues of research are likely to be most promising in the design of novel magnetoelectric materials.

Growth and lattice dynamics of ultrathin BaO films on Pt(0 0 1)

We report on the surface phonons of long-range ordered BaO thin films grown on Pt(0 0 1). In the thickness range between 4 and 28 ML, we find unstrained BaO(0 0 1)-([Formula: see text]) bulk-like terminated films which coincide with a Pt(0 0 1)-c([Formula: see text]) lattice periodicity. The dipole-active lattice vibrations were determined using high-resolution electron energy loss spectroscopy. For all BaO film thicknesses, a single Fuchs-Kliewer phonon-polariton is observed. Its intensity increases and its frequency softens with increasing film thickness. These thickness-dependent properties and the spectral shape are quantitatively discussed on the basis of dielectric theory, where a proper modeling requires three components: the dielectric response of the BaO film itself, the plasmonic response of the metallic substrate, and a weak damping due to a defect-induced doping within the oxide film. For a full description, also the quantization of the phonon wavevector due to the confinement within the film of finite thickness has to be taken into account.

Strain-induced ferroelectricity and lattice coupling in BaSnO3 and SrSnO3

Perovskite stannates such as BaSnO₃ and SrSnO₃ exhibit promising photovoltaic properties, and hold promise for application in solar cell devices. However, the lack of ferroelectricity and the wide band gap in these materials limit their potential for photovoltaic applications. Here, by first-principles calculations, we demonstrate the realization of a primary ferroelectric polarization in non-ferroelectric BaSnO₃ and SrSnO₃ through strain engineering. In addition to the appearance of polarization, the band gaps of the materials are greatly narrowed when the paraelectric to ferroelectric phase transition takes place under compressive strain. Furthermore, an intriguing Q₂ mode triggered by lattice coupling with the polar mode is found in the stannates subjected to a sufficient tensile strain and this mode has a significant effect on the band gap, which suggests another pathway to narrow the band gap through the electric field control of the Q₂ mode. The fruitful electronic, structural, and energetic properties are discussed in detail to achieve a fundamental understanding of the strain-induced ferroelectricity, tunable band gap, and lattice couplings between the Q₂ mode and different polar/rotational distortions in the perovskite stannates.

First-principles study on the electric structure and ferroelectricity in epitaxial CsSnI3 films

Epitaxial strain can induce ferroelectricity and enhance the polarization in the inorganic hybrid perovskite material CsSnI₃ for energy applications.

Electric Field Control of Jahn-Teller Distortions in Bulk Perovskites

The Jahn-Teller distortion, by its very nature, is often at the heart of the various electronic properties displayed by perovskites and related materials. Despite the Jahn-Teller mode being nonpolar, we devise and demonstrate, in the present Letter, an electric field control of Jahn-Teller distortions in bulk perovskites. The electric field control is enabled through an anharmonic lattice mode coupling between the Jahn-Teller distortion and a polar mode. We confirm this coupling and quantify it through first-principles calculations. The coupling will always exist within the Pb2_{1}m space group, which is found to be the favored ground state for various perovskites under sufficient tensile epitaxial strain. Intriguingly, the calculations reveal that this mechanism is not only restricted to Jahn-Teller active systems, promising a general route to tune or induce novel electronic functionality in perovskites as a whole.

Structural phase transition and spin reorientation of LaFeO3 films under epitaxial strain

Structural phase transition and spin reorientation of orthoferrites LaFeO₃ epitaxially grown along the pseudocubic (001) direction are investigated based on first-principles calculation.

Large polarization and dielectric response in epitaxial SrZrO3 films

First-principles calculations are performed to investigate the ferroelectric and dielectric properties of (001) epitaxial SrZrO₃ thin films under misfit strain.

Large Elasto-Optic Effect in Epitaxial PbTiO(3) Films

First-principles calculations are performed to investigate the elasto-optic properties of four different structural phases in (001) epitaxial PbTiO(3) films under tensile strain: a tetragonal (T) phase and an orthorhombic (O) phase, which are the ground states for small and large strain, respectively, and two low-symmetry, monoclinic phases of Cm and Pm symmetries that have low total energy in the intermediate strain range. It is found that the refractive indices of the T and O phases respond differently to epitaxial strain, evidenced by a change of sign of their effective elasto-optic coefficients, and as a result of presently discovered correlations between refractive index, axial ratio, and magnitude of the ferroelectric polarization. The difference in refractive indices between T and O and the existence of such correlations naturally lead to large elasto-optic coefficients in the Cm and Pm states in the intermediate strain range, because Cm structurally bridges the T and O phases (via polarization rotation and a rapid change of its axial ratio) and Pm adopts a similar axial ratio and polarization magnitude to Cm. The present results therefore broaden the palette of functionalities of ferroelectric materials, and suggest new routes to generate systems with unprecedentedly large elasto-optic conversion.

Predicted energetics and properties of rare-earth ferrites films grown on cubic (111)- and hexagonal (0001)-oriented substrates

First-principles calculations are performed to compare the energetics of several phases, including hexagonal polar P6(3)cm and perovskite non-polar Pbnm-like states, of epitaxial RFeO3 films (with R  =  Ce, Pr, Nd, Sm, Gd, Tb, Dy, Ho, Er and Lu) grown on different cubic (1 1 1)- and hexagonal (0 0 0 1)-oriented substrates. The P63cm phase is found to be the ground state for large enough in-plane lattice parameters in all investigated RFeO3 films, and its polarization is tunable by the amount of epitaxial strain. Series of available substrates allowing the growth of hexagonal polar RFeO3 films, as well as other phenomena of fundamental and technological importance (e.g. different ground states and coexistence between several phases) are also predicted.

BiFeO3 epitaxial thin films and devices: past, present and future

The celebrated renaissance of the multiferroics family over the past ten years has also been that of its most paradigmatic member, bismuth ferrite (BiFeO3). Known since the 1960s to be a high temperature antiferromagnet and since the 1970s to be ferroelectric, BiFeO3 only had its bulk ferroic properties clarified in the mid-2000s. It is however the fabrication of BiFeO3 thin films and their integration into epitaxial oxide heterostructures that have fully revealed its extraordinarily broad palette of functionalities. Here we review the first decade of research on BiFeO3 films, restricting ourselves to epitaxial structures. We discuss how thickness and epitaxial strain influence not only the unit cell parameters, but also the crystal structure, illustrated for instance by the discovery of the so-called T-like phase of BiFeO3. We then present its ferroelectric and piezoelectric properties and their evolution near morphotropic phase boundaries. Magnetic properties and their modification by thickness and strain effects, as well as optical parameters, are covered. Finally, we highlight various types of devices based on BiFeO3 in electronics, spintronics, and optics, and provide perspectives for the development of further multifunctional devices for information technology and energy harvesting.

Creating multiferroics with large tunable electrical polarization from paraelectric rare-earth orthoferrites

The quest for materials possessing both a magnetic ordering temperature above room temperature and a large electrical polarization is an important research direction in order to design novel spintronic and memory devices. Up to now, BiFeO3 and related systems are the only known compounds simultaneously possessing such characteristics. Here, first-principles calculations predict that another family of materials, namely epitaxial films made of rare-earth orthoferrites (RFeO3), can also exhibit such desired features. As a matter of fact, applying a large enough strain to these compounds, which are nominally paraelectric and have a high magnetic transition temperature, is predicted to render them ferroelectric, and thus multiferroic. At high compressive strain, the resulting ferroelectric phase of RFeO3 systems having large rare-earth ions is even a tetragonal state characterized by a giant polarization and axial ratio. For large tensile strain, two striking inhomogenous ferroelectric phases--including one never observed before in any perovskite--are further predicted as having significant polarization. A multiphase boundary also occurs, which may lead to optimization of properties or unusual features. Finally, many quantities, including electrical polarization and magnetic ordering temperature, are tunable by varying the epitaxial strain and/or chemical pressure.

Structural instability of epitaxial (001) BiFeO₃ thin films under tensile strain

We explore BiFeO3 under tensile strain using first-principles calculations. We find that the actual structures are more complex than what had been previously thought, and that there is a strong shear deformation type structural instability which modifies the properties. Specifically, we find that normal tensile strain leads to structural instabilities with a large induced shear deformation in (001) BiFeO3 thin films. These induced shear deformations in (001) BiFeO3 thin films under tension stabilize the (001) BiFeO3 thin films and lead to Cc and Ima2 phases that are more stable than the Pmc21 phase at high tensile strain. The induced shear deformation shifts the Cc to Ima2 phase transition towards lower tensile strain region (~1% less), prevents monoclinic tilt and oxygen octahedral tilts, and increases the ferroelectric polarization. The induced shear deformation also strongly affects the electronic structure. The results are discussed in relation to growth of BiFeO3 thin films on cubic and tetragonal substrates involving high levels of tensile strain.

First-principles study on the phase diagram and multiferroic properties of (SrCoO3)1/(SrTiO3)1 superlattices

To design a multiferroic material at atomic scale, strong spin-lattice and charge-lattice couplings play crucial roles. Our first-principles calculation on (SrCoO₃)₁/(SrTiO₃)₁ superlattices, with above coupling properties, yields a rich physical phase diagram as a function of epitaxial strain. In particular, a robust ferroelectric ferromagnetic insulator of Pc symmetry is stabilized at tensile strain Δa/a₀ = 0.86%–5.53%. The polarization can be as large as 36 μC/cm² and magnetic moment can reach 6μ_B per unit cell. The magnetocrystalline anisotropy energy (0.16 meV/Co in (001) plane, 0.6 meV/Co in (100) plane) is comparable with that of TbMnO₃ compound and the magnetoelectric constant α (1.44 × 10⁻³ Gaussian unit) is comparable with that of Co₃B₇O₁₃Br compound. Our study suggests that epitaxially strained (SrCoO₃)₁/(SrTiO₃)₁ superlattices not only offer an excellent candidate for multiferroic materials, but also demonstrate the half-metal and ferromagnetic insulator properties with potential application in spintronic devices.

Prediction of a novel magnetoelectric switching mechanism in multiferroics

We report a first-principles study of the recently predicted Pmc21 phase of the multiferroic BiFeO3 material, revealing a novel magnetoelectric effect that makes it possible to control magnetism with an electric field. The effect can be viewed as a two-step process: Switching the polarization first results in the change of the sense of the rotation of the oxygen octahedra, which in turn induces the switching of the secondary magnetic order parameter. The first step is governed by an original trilinear-coupling energy between polarization, octahedral tilting, and an antiferroelectric distortion. The second step is controlled by another trilinear coupling, this one involving the predominant and secondary magnetic orders as well as the oxygen octahedral tilting. In contrast with other trilinear-coupling effects in the literature, the present ones occur in a simple ABO3 perovskite and involve a large polarization.

BiFeO3 domain wall energies and structures: a combined experimental and density functional theory+U study

We determined the atomic structures and energies of 109°, 180°, and 71° domain walls in BiFeO3, combining density functional theory+U calculations and aberration-corrected transmission electron microscopy images. We find a substantial Bi sublattice shift and a rather uniform Fe sublattice across the walls. The calculated wall energies (γ) follow the sequence γ109<γ180<γ71 for the 109°, 180°, and 71° walls. We attribute the high 71° wall energy to an opposite tilting rotation of the oxygen octahedra and the low 109° wall energy to the opposite twisting rotation of the oxygen octahedra across the domain walls.

Ferroelectric domains in multiferroic BiFeO3 films under epitaxial strains

First-principles calculations are performed to investigate energetic and atomistic characteristics of ferroelectric domains walls (DWs) of BiFeO(3) (BFO) films subject to compressive strain. Significantly lower DW energies than those previously reported, and a different energetic hierarchy between the various DW types, are found for small strains. In all investigated cases (corresponding to ideal angles of 71°, 109°, and 180° formed by the domain polarizations), the DW energy reaches its maximum value for misfit strains that are around the critical strain at which the transition between the R-like and T-like phases occurs in single-domain BFO. Near these strains, several quantities depend strongly on the type of domain wall; such distinct behavior is associated with an elastic difference and a large out-of-plane polarization at the DW in the 180° case. A further increase of the magnitude of the strain leads to (i) a change of hierarchy of the DW energies, (ii) large out-of-plane polarizations inside each up and down domain, and (iii) novel atomic arrangements at the domain walls. Our study can thus initiate a new research direction, namely strain engineering of domain-wall functionalities.

Understanding and revisiting properties of EuTiO3 bulk material and films from first principles

Ab initio computations are performed to investigate properties of bulk material and epitaxial films made of EuTiO3 (ETO). A whole family of nanoscale twinned phases, that present complex oxygen octahedra tilting (OOT) and unusual antiferroelectricity, is found to be degenerate in energy with simpler phases (all possessing typical OOT) in bulk ETO. Such degeneracy provides a successful explanation of recently observed anomalous phenomena. The calculations also lead to revisiting the (rich) phase diagram of ETO films.