Search

Polarons can naturally form in materials from the interaction of extra charge carriers with the atomic lattice. Ubiquitous, they are central to various topics and phenomena such as high-T$_c$ superconductivity, electrochromism, photovoltaics, photocatalysis or ion batteries. However, polaron formation remains poorly understood and mostly relies on few historical models such as Landau-Pekar, Fr\"olich, Holstein or Jahn-Teller polarons. Here, from advanced first-principles calculations, we show that the formation of intriguing medium-size polarons in WO$_3$ does not fit with traditional models but instead arises from the undoing of distortive atomic motions inherent to the pristine phase, which lowers the bandgap through dynamical covalency effects. We so introduce the innovative concept of {\it anti-distortive} polaron and rationalize it from a quantum-dot model. We demonstrate that anti-distortive polarons are generic to different families of compounds and clarify how this new concept opens concrete perspectives for a better control of the polaronic state and related properties., Comment: 9 pages, 2 figures

We compute from density functional theory calculations the dynamical magnetic effective charges of the Jahn-Teller ground state of KCoF$_3$ and we unexpectedly obtain the largest up-to-date values reported in real materials (close to 200$\times 10^{-2}\mu_{\text{B}}/\text{\r{A}}$). We determine that those large values are due to both spin and orbital contributions, as expected from the presence of unquenched Co-d orbital magnetic moment, even though we show that the Co dynamical magnetic effective charges values are surprisingly zero by symmetry. Hence, we conclude that even if KCoF$_3$ is not multiferroic, it is antimagnetoelectric with a compensating but large individual magnetic sublattice magnetoelectric response of 210 ps/m per spin-channel. Our study shows that the dynamical magnetic effective charges might be an important property of magnetic materials beyond the magnetoelectric ones.

Chirality refers to the asymmetry of objects that cannot be superimposed on their mirror image. It is a concept that exists in various scientific fields and has profound consequences. Although these are perhaps most widely recognized within biology, chemistry, and pharmacology, recent advances in chiral phonons, topological systems, crystal enantiomorphic materials, and magneto-chiral materials have brought this topic to the forefront of condensed matter physics research. Our review discusses the symmetry requirements and the features associated with structural chirality in inorganic materials. This allows us to explore the nature of phase transitions in these systems, the coupling between order parameters, and their impact on the material's physical properties. We highlight essential contributions to the field, particularly recent progress in the study of chiral phonons, altermagnetism, magnetochirality between others. Despite the rarity of naturally occurring inorganic chiral crystals, this review also highlights a significant knowledge gap, presenting challenges and opportunities for structural chirality mostly at the fundamental level, e.g., chiral displacive phase transitions and ferrochirality, possibilities of tuning and switching structural chirality by external means (electric, magnetic, or strain fields), whether chirality could be an independent order parameter, and whether structural chirality could be quantified, etc. Beyond simply summarising this field of research, this review aims to inspire further research in materials science by addressing future challenges, encouraging the exploration of chirality beyond traditional boundaries, and seeking the development of innovative materials with superior or new properties.

Using first-principles calculations, we explore the magnetoelectric properties of the room-temperature multiferroic crystal BiCoO$_3$. We use both applied magnetic field and finite-difference techniques to show that BiCoO$_3$ is anti-magnetoelectric at the linear level. The calculation of the dynamical effective charges reveals that the total magnetoelectric response is zero due to the compensating non-zero magnetoelectric response of each magnetic sublattice. This calculation also highlights that the the orbital contribution to the response is remarkably larger than the spin one and that each sublattice has a rather large total magnetoelectric response of 85 ps/m. Furthermore, we provide an intuitive recipe to visualize the dynamical magnetic effective charge, allowing to examine its multipolar nature which we confirm by means of ab initio calculations. Given the large value of the local response, we investigate the ferromagnetic phase as well, which gives a giant magnetoelectric response of about 1000 ps/m and coming mainly from the spin contribution this time. Finally, we discuss the possible reasons for such a large magnetoelectric response in BiCoO3 and propose possible strategies to unveil this potentially large response.

We compare the various chirality measures most widely used in the literature to quantify chiral symmetry in extended solids, i.e., the continuous chirality measure, the Hausdorff distance, and the angular momentum. By studying these functions in an algebraically tractable case, we can evaluate their strengths and weaknesses when applied to more complex crystals. Going beyond those classical calculations, we propose a new method to quantify the handedness of a crystal based on a pseudoscalar function, i.e., the helicity. This quantity, borrowed from hydrodynamics, can be computed from the eigenvector carrying the system from the high-symmetry non-chiral phase to the low-symmetry chiral phase. Different model systems like K$_3$NiO$_2$, CsCuCl$_3$ and MgTi$_2$O$_4$ are used as test cases where we show the superior interest of using helicity to quantify chirality together with the handedness distinction.

We demonstrate through first-principles calculations applied to the K$_{3}$NiO$_{2}$ crystal that a structural phase transition from an achiral to a chiral phase can be mediated by a degenerate soft phonon mode and controlled by pressure and epitaxial strain. Breaking such degeneracy with an electric field generates a competition between enantiomorphic, polar, and orthorhombic displacements. Originated by the interaction between spontaneous chiral and induced polar and axial modes, an optimal parameter window for converting the handedness of the system into its opposite kind is observed.

We study from first-principles calculations the ferroelectric structural phase transition of Pb$_5$Ge$_3$O$_{11}$ crystal. The calculations of phonons and Born effective charges of the paraelectric phase allow us to identify a polar instability that is unstable in both transverse-optic and longitudinal-optic versions, giving rise to an entire branch of instability along a propagation vector parallel to the mode polarization (the hexagonal axe). This is the hint of hyperferroelectricity and the stable head-to-head and tail-to-tail domain, as recently reported from both experiments and theory. Then, our analysis of the ferroelectric phase shows that the polarization of Pb$_5$Ge$_3$O$_{11}$ is uniaxial along the hexagonal axes and with small in-plane components due to a piezoelectric effect. The symmetry-adapted mode analysis shows that the total ferroelectric ground state distortion comes mainly from polar distortions of the unstable polar phonon mode but also from an invariant, cooperative mode that amplifies the polar deformation. We also build a phenomenological model that highlights how the coupling between these modes is at play and helps us understand how to reproduce the second-order phase transition. At last, we also quantify the structural chirality through the continuous symmetry measure method and trace its origin to the polar unstable mode itself. By extending our approach to the phonon states we further show that the chirality is poorly affected by the relaxation but could also be enhanced by activating high frequency modes with polar symmetry. Finally we study the phonon angular momentum (AM) distribution in both phases and identify trends in the AM behaviour across the Brillouin zone.

We establish the sign of the linear magnetoelectric (ME) coefficient, $\alpha$, in chromia, Cr$_2$O$_3$. Cr$_2$O$_3$ is the prototypical linear ME material, in which an electric (magnetic) field induces a linearly proportional magnetization (polarization), and a single magnetic domain can be selected by annealing in combined magnetic (H) and electric (E) fields. Opposite antiferromagnetic domains have opposite ME responses, and which antiferromagnetic domain corresponds to which sign of response has previously been unclear. We use density functional theory (DFT) to calculate the magnetic response of a single antiferromagnetic domain of Cr$_2$O$_3$ to an applied in-plane electric field at 0 K. We find that the domain with nearest neighbor magnetic moments oriented away from (towards) each other has a negative (positive) in-plane ME coefficient, $\alpha_{\perp}$, at 0 K. We show that this sign is consistent with all other DFT calculations in the literature that specified the domain orientation, independent of the choice of DFT code or functional, the method used to apply the field, and whether the direct (magnetic field) or inverse (electric field) ME response was calculated. Next, we reanalyze our previously published spherical neutron polarimetry data to determine the antiferromagnetic domain produced by annealing in combined E and H fields oriented along the crystallographic symmetry axis at room temperature. We find that the antiferromagnetic domain with nearest-neighbor magnetic moments oriented away from (towards) each other is produced by annealing in (anti-)parallel E and H fields, corresponding to a positive (negative) axial ME coefficient, $\alpha_{\parallel}$, at room temperature. Since $\alpha_{\perp}$ at 0 K and $\alpha_{\parallel}$ at room temperature are known to be of opposite sign, our computational and experimental results are consistent., Comment: 11 pages, 5 figures

We explore from first-principles calculations the ferroelectric material Pb5Ge3O11 as a model for controlling the spin-orbit interaction (SOC) in crystalline solids. The SOC has a surprisingly strong effect on the structural energy landscape by deepening the ferroelectric double well. We observe that this effect comes from a specific Pb Wyckoff site that lies on the verge of a natural cavity channel of the crystal. We also find that a unique cavity state is formed by the empty 6p states of another Pb site at the edge of the cavity channel. This cavity state exhibits a sizeable spin splitting with a mixed Rashba-Weyl character and a topologically protected crossing of the related bands. We also show that the ferroelectric properties and the significant SOC effects are exceptionally robust against n-doping up to several electrons per unit cell. We trace the provenance of these original effects to the unique combination of the structural cavity channel and the chemistry of the Pb atoms with 6p orbitals localizing inside the channel.

Out-of-plane polar domain structures have recently been discovered in strained and twisted bilayers of inversion symmetry broken systems such as hexagonal boron nitride. Here we show that this symmetry breaking also gives rise to an in-plane component of polarization, and the form of the total polarization is determined purely from symmetry considerations. The in-plane component of the polarization makes the polar domains in strained and twisted bilayers topologically non-trivial, forming a network of merons and antimerons (half-skyrmions and half-antiskyrmions). For twisted systems, the merons are of Bloch type whereas for strained systems they are of N\'eel type. We propose that the polar domains in strained or twisted bilayers may serve as a platform for exploring topological physics in layered materials and discuss how control over topological phases and phase transitions may be achieved in such systems., Comment: Accepted version

We present a general relation between the electromechanical couplings of infinitesimal strain and electric field to arbitrary order, measured at fixed voltage and at fixed electric field. We show that the improper response at fixed field can be written as the strain derivative of the $n^{\text{th}}$ order susceptibility tensor, and the proper response at fixed voltage drop can be written as the response at fixed field plus corrections for dilations and 90$^{\circ}$ rotations induced by strain. Our theory correctly reproduces the proper piezoelectric response and we go beyond with the electrostrictive response. We present first-principles calculations of the improper electrostrictive response at fixed field, and illustrate how the correction is used to obtain the proper response at fixed voltage. This distinction is of high importance given the recent interest in giant electrostrictors exhibiting electromechanical responses as large as the piezoelectric ones., Comment: Accepted version

Lead vanadate PbVO3 is a polar crystal with a P4mm space group at ambient conditions. It is isostructural with the model soft-mode driven ferroelectric PbTiO3, but differs from it by the so-called 'supertetragonal' elongation of its unit cell. In this paper, we report a combined study of the lattice dynamics of PbVO3 by Raman spectroscopy at room temperature and first-principle calculations. All zone-center transverse optical (TO) phonon modes are identified by polarized, angle-dependent Raman spectroscopy and assigned as follows: E modes at 136, 269, 374 and 508 cm-1, A1 modes at 188, 429 and 874 cm-1 and B1 mode at 319 cm-1. The calculations confirm the experimental symmetry assignment and allow to obtain the longitudinal (LO) phonons wavenumbers. Besides, we analyze the mode eigenvectors in detail, in order to identify the atomic displacements associated with each mode and compare them with PbTiO3. In spite of their differences in chemistry and strain, the phonon eigenvectors are found to be remarkably comparable in both compounds. We discuss the position of the ferroelectric soft mode in PbVO3 as compared to PbTiO3. A sizeable splitting of the B1+E modes appears as a characteristic feature of supertetragonal phases. The peculiarity of the vanadyl V-O bond frequency in PbVO3 is also addressed., Comment: 11 pages, 8 figures, 6 Tables, Supplemental Information 3 figures and 6 Tables, under review

Epitaxial SrTiO$_3$ (STO) on Si is nowadays the benchmark initial platform for the further addition of functional oxides on Si. Starting the growth of STO on a Sr-passivated Si substrate with 1/2 monolayer (ML) Sr coverage and a (1 $\times$ 2) reconstructed Si surface with rows of Si dimers, the final STO/Sr/Si stack exhibits 1 ML Sr coverage and a (1 $\times$ 1) Si surface without dimer. Using first-principles density functional theory calculations, we investigate how the interface evolves from 1/2 ML to 1 ML Sr coverage, concluding that the latter is indeed most stable and that the reconstruction of the interface takes place during the early stage of the layer-by-layer deposition. Going further, we determine the band alignment of the final stable interface and assess its potential interest as photocathode for water reduction., Comment: 9 pages, 6 figures

We investigate the electrostrictive response across a ferroelectric phase transition from first-principles calculations and refute the prevailing view of constant electrostriction across the ferroelectric phase boundary. We take as a case study the epitaxial strain-induced transition from para- to feroelectricity of \ce{KTaO3}. We show that the magnitude of the electrostriction diverges with the permitivity at the transition, hence exhibiting giant responses through a calculation of both the M and Q electrostrictive tensors. We explain the origin of this giant electrostrictive response in \ce{KTaO3} using a microscopic decomposition of the electrostriction coefficients, and use this understanding to propose design rules for the development of future giant electrostrictors for electromechanical applications. Finally, we introduce a further means to calculate electrostriction, specific to ferroelectrics, and not yet utilised in the literature.

Multiferroics are a unique class of materials where magnetic and ferroelectric orders coexist. The research on multiferroics contributes significantly to the fundamental understanding of the strong correlations between different material degrees of freedom and provides an energy-efficient route toward the electrical control of magnetism. While multiple ABO3 oxide perovskites have been identified as being multiferroic, their magnetoelectric coupling strength is often weak, necessitating the material search in different compounds. Here, we report the observation of room-temperature multiferroic orders in multi-anion SrNbO3-xNx thin films. In these samples, the multi-anion state enables the room-temperature ferromagnetic ordering of the Nb d-electrons. Simultaneously, we find ferroelectric responses that originate from the structural symmetry breaking associated with both the off-center displacements of Nb and the geometric displacements of Sr, depending on the relative O-N arrangements within the Nb-centered octahedra. Our findings not only diversify the available multiferroic material pool but also demonstrate a new multiferroism design strategy via multi-anion engineering.

We report a theoretical study of the non-linear magnetoelectric response of GdFeO$_3$ through an analytical approach combined with a Heisenberg model which is fitted against first-principles calculations. Our theory reproduces the non-linear change of polarization under applied magnetic field reported experimentally such that it allows to analyze the origin of the large responses in the different directions. We show that the non-linear character of the response in these materials originates from the fact that the antiferromagnetic order of Gd atoms changes non-linearly with respect to the applied magnetic field. Our model can be generalized to other materials in which the antiferromagnetic ordering breaks inversion symmetry.

Raman spectroscopy is an advantageous method for studying the local structure of materials, but the interpretation of measured spectra is complicated by the presence of oblique phonons in polycrystals of polar materials. Whilst group theory considerations and standard ab initio calculations are helpful, they are often valid only for single crystals. In this paper, we introduce a method for computing Raman spectra of polycrystalline materials from first principles. We start from the standard approach based on the (Placzek) rotation invariants of the Raman tensors and extend it to include the effect of the coupling between the lattice vibrations and the induced electric field, and the electro-optic contribution, relevant for polar materials like ferroelectrics. As exemplified by applying the method to rhombohedral BaTiO3, AlN, and LiNbO3, such an extension brings the simulated Raman spectrum to a much better correspondence with the experimental one. Additional advantages of the method are that it is general, permits automation, and thus can be used in high-throughput fashion., Comment: 23 pages, 4 figures

Spin reorientation and magnetisation reversal are two important features of the rare-earth orthorhombic provskites ($RM$O$_{3}$'s) that have attracted a lot of attention, though their exact microscopic origin has eluded researchers. Here, using density functional theory and classical atomistic spin dynamics we build a general Heisenberg magnetic model that allows to explore the whole phase diagram of the chromite and ferrite compounds and to scrutinize the microscopic mechanism responsible for spin reorientations and magnetisation reversals. We show that the occurrence of a magnetization reversal transition depends on the relative strength and sign of two interactions between rare-earth and transition-metal atoms: superexchange and Dzyaloshinsky-Moriya. We also conclude that the presence of a smooth spin reorientation transition between the so-called $\Gamma_4$ and the $\Gamma_2$ phases through a coexisting region, and the temperature range in which it occurs, depends on subtle balance of metal--metal (superexchange and Dzyaloshinsky-Moriya) and metal--rare-earth (Dzyaloshinsky-Moriya) couplings. In particular, we show that the intermediate coexistence region occurs because the spin sublattices rotate at different rates.

In this work we present a new method for the calculation of the electrostrictive properties of materials using density functional theory. The method relies on the thermodynamical equivalence, in a dielectric, of the quadratic mechanical responses (stress or strain) to applied electric stimulus (electric or polarisation fields) to the strain or stress dependence of its dielectric susceptibility or stiffness tensors. Comparing with current finite-field methodologies for the calculation of electrostriction, we demonstrate that our presented methodology offers significant advantages of efficiency, robustness, and ease of use. These advantages render tractable the highthroughput theoretical investigation into the largely unknown electrostrictive properties of materials.

Third-generation photovoltaic (PV) materials combine many advantageous properties, including a high optical absorption together with a large charge carrier mobility, facilitated by small effective masses. Halide perovskites (ABX3, where X is I, Br or Cl) appear to be the most promising third-generation PV materials at present. Their opto-electronic properties are governed by the B-X bond. A quantum-chemical bond analysis reveals that this bond differs significantly from ionic, metallic or covalent bonds. Instead, it is better regarded as metavalent, since it shares approximately one p-electron between adjacent atoms. The resulting sigma-bond is half-filled, which causes pronounced optical absorption. Electron transfer and lattice distortions open a moderate band gap, resulting in charge carriers with small effective masses. Hence metavalent bonding explains the favorable PV properties of halide perovskites. This is summarized in a map for different bond types, which provides a blueprint to design third-generation PV materials.

Combining ferroelectricity with other properties such as visible light absorption or long-range magnetic order requires the discovery of new families of ferroelectric materials. Here, through the analysis of a high-throughput database of phonon band structures, we identify a new structural family of anti-Ruddlesden-Popper phases A$_4$X$_2$O (A=Ca, Sr, Ba, Eu, X=Sb, P, As, Bi) showing ferroelectric and anti-ferroelectric behaviors. The discovered ferroelectrics belong to the new class of hyperferroelectrics which polarize even under open-circuit boundary conditions. The polar distortion involves the movement of O anions against apical A cations and is driven by geometric effects resulting from internal chemical strains. Within this new structural family, we show that Eu$_4$Sb$_2$O combines coupled ferromagnetic and ferroelectric order at the same atomic site, a very rare occurrence in materials physics., Comment: 21 pages, 4 figures with supplemental information attached

We present TB2J, a Python package for the automatic computation of magnetic interactions, including exchange and Dzyaloshinskii-Moriya interactions, between atoms of magnetic crystals from the results of density functional calculations. The program is based on the Green's function method with the local rigid spin rotation treated as a perturbation. As input,the package uses the output of either Wannier90, which is interfaced with many density functional theory packages,or of codes based on localized orbitals. A minimal user input is needed, which allows for easy integration into high-throughput workflows. The package is open source under BSD 2-Clause license, available at https://github.com/mailhexu/TB2J.

Here, we present the theoretical analysis of the structural and electronic degrees of freedom of two different oxide/fluoride perovskite superlattices, KTaO$_3$/KZnF$_3$ and KTaO$_3$/KNiF$_3$. Using first-principles calculations, we found the appearance of a two-dimensional electron, 2DEG, and hole, 2DHG, gases as a function of the number of layers of the different pristine materials. We demonstrate that the phonon-dynamics at the KTO/K$M$F superlattices play a crucial role in the appearance of these effects. Additionally, our results reveal a rather sizeable Rashba-type spin-splitting at these interfaces in comparison to another oxide/oxide counterparts.

The PyProcar Python package plots the band structure and the Fermi surface as a function of site and/or s,p,d,f - projected wavefunctions obtained for each $k$-point in the Brillouin zone and band in an electronic structure calculation. This can be performed on top of any electronic structure code, as long as the band and projection information is written in the PROCAR format, as done by the VASP and ABINIT codes. PyProcar can be easily modified to read other formats as well. This package is particularly suitable for understanding atomic effects into the band structure, Fermi surface, spin texture, etc. PyProcar can be conveniently used in a command line mode, where each one of the parameters define a plot property. In the case of Fermi-surfaces, the package is able to plot the surface with colors depending on other properties such as the electron velocity or spin projection. The mesh used to calculate the property does not need to be the same as the one used to obtain the Fermi surface. A file with a specific property evaluated for each $k$-point in a $k-$mesh and for each band can be used to project other properties such as electron-phonon mean path, Fermi velocity, electron effective mass, etc. Another existing feature refers to the band unfolding of supercell calculations into predefined unit cells., Comment: 32 pages, 14 figures, Submitted to "Computer Physics Communications", project url: "https://github.com/romerogroup/pyprocar"

Ferroelectric Rashba semiconductors (FERSC), in which Rashba spin-splitting can be controlled and reversed by an electric field, have recently emerged as a new class of functional materials useful for spintronic applications. The development of concrete devices based on such materials is, however, still hampered by the lack of robust FERSC compounds. Here, we show that the coexistence of large spontaneous polarisation and sizeable spin-orbit coupling is not sufficient to have strong Rashba effects and clarify why simple ferroelectric oxide perovskites with transition metal at the B-site are typically not suitable FERSC candidates. By rationalizing how this limitation can be by-passed through band engineering of the electronic structure in layered perovskites, we identify the Bi$_2$WO$_6$ Aurivillius crystal as the first robust ferroelectric with large and reversible Rashba spin-splitting, that can even be substantially doped without losing its ferroelectric properties. Importantly, we highlight that a unidirectional spin-orbit field arises in layered Bi$_2$WO$_6$, resulting in a protection against spin-decoherence.We highlight moreover that a unidirectional spin-orbit field arises in Bi$_2$WO$_6$, in which the spin-texture is so protected against spin-decoherence.

We extend the density functional perturbation theory formalism to the case of non-collinear magnetism. The main problem comes with the exchange-correlation (XC) potential derivatives, which are the only ones that are affected by the non-collinearity of the system. Most of the present XC functionals are constructed at the collinear level, such that the off-diagonal (containing magnetization densities along $x$ and $y$ directions) derivatives cannot be calculated simply in the non-collinear framework. To solve this problem, we consider here possibilities to transform the non-collinear XC derivatives to a local collinear basis, where the $z$ axis is aligned with the local magnetization at each point. The two methods we explore are i) expanding the spin rotation matrix as a Taylor series, ii) evaluating explicitly the XC for the local density approximation through an analytical expression of the expansion terms. We compare the two methods and describe their practical implementation. We show their application for atomic displacement and electric field perturbations at the second order, within the norm-conserving pseudopotential methods.

This version withdrawn by arXiv administrators because the author did not have the right to agree to our license at the time of submission, Comment: This version withdrawn by arXiv administrators because the author did not have the right to agree to our license at the time of submission

Entropy is a fundamental thermodynamic quantity that is a measure of the accessible microstates available to a system, with the stability of a system determined by the magnitude of the total entropy of the system. This is valid across truly mind boggling length scales - from nanoparticles to galaxies. However, quantitative measurements of entropy change using calorimetry are predominantly macroscopic, with direct atomic scale measurements being exceedingly rare. Here for the first time, we experimentally quantify the polar configurational entropy (in meV/K) using sub-\r{a}ngstr\"{o}m resolution aberration corrected scanning transmission electron microscopy. This is performed in a single crystal of the prototypical ferroelectric $\mathsf{LiNbO_3}$ through the quantification of the niobium and oxygen atom column deviations from their paraelectric positions. Significant excursions of the niobium - oxygen polar displacement away from its symmetry constrained direction is seen in single domain regions which increases in the proximity of domain walls. Combined with first principles theory plus mean field effective Hamiltonian methods, we demonstrate the variability in the polar order parameter, which is stabilized by an increase in the magnitude of the configurational entropy. This study presents a powerful tool to quantify entropy from atomic displacements and demonstrates its dominant role in local symmetry breaking at finite temperatures in classic, nominally Ising ferroelectrics., Comment: 23 pages, 21 figures (5 main, 16 supplemental)

Despite oxides and some fluorides perovskites have emerged as prototypes of multiferroic and magnetoelectric materials, they have not impacted real devices. Unfortunately, their working temperatures are very low and the magnetoelectric coupling has been reported to be rather small. Herewith, we report from first-principles calculations an ideal magnetization reversal through polarization switching in BaCuF$_4$ which, according to our results, could be achieved close to room temperature. We also show that this ideal coupling is driven by a soft mode that combines both, polarization and octahedral rotation. The later being directly coupled to the weak ferromagnetism of BaCuF$_4$. This, added to its strong Jahn-Teller distortion and its orbital ordering, makes this material as a very appealing prototype for crystals in the $ABX_4$ family in multifunctional applications.

Despite of the tremendous efforts spent on the oxygen vacancy migration in determining the property optimization of oxygen-vacancy enrichment transition metal oxides, few has focused on their dynamic behaviors non-equilibrium states. In this work, we performed multi-timescale ultrafast X-ray diffraction measurements by using picosecond synchrotron X-ray pulses and femtosecond table-top X-ray pulses to monitor the structural dynamics in the oxygen-vacancy ordered SrCoO2.5 thin films. A giant photoinduced strain ({\Delta}c/c > 1%) was observed, whose distinct correlation with the pump photon energy indicates a non-thermal origin of the photoinduced strain. The sub-picosecond resolution X-ray diffraction reveals the formation and propagation of the coherent acoustic phonons inside the film. We also simulate the effect of photoexcited electron-hole pairs and the resulting lattice changes using the Density Function Theory method to obtain further insight on the microscopic mechanism of the measured photostriction effect. Comparable photostrictive responses and the strong dependence on excitation wavelength are predicted, revealing a bonding to anti-bonding charge transfer or high spin to intermediate spin crossover induced lattice expansion in the oxygen-vacancy films., Comment: 12 pages, 4 figures, support material

The effect of epitaxial strain on structural, magnetic and electronic properties of BaFeO 3 per- ovskite oxide are investigated from first principles calculations, using the Density Functional The- ory (DFT) plus the Hubbard approach (DFT+U) within the Generalized Gradient Approximation (GGA). Hybrid functional calculations, based on mixed exact Hartree-Fock (HF) and DFT exchange energy functionals, are also performed. For the ground state calculations,the DFT+U is found more suitable to describe the half metallic and ferromagnetic state of cubic BaFeO 3 . The possible oc- curence of oxygen octahedra rotations, Jahn-Teller distortions and charge orderings through biaxial strain are explored. The obtained results reveal that the Jahn-teller distortion is induced under tensile biaxial strain while the oxygen octahedra rotations and breathing are unusualy not observed. Then, the strained BaFeO 3 is considered as a particular Jahn-Teller distorted perovskite with excep- tional properties when compared to CaFeO 3 and SrFeO 3 . These findings lead to a strain engineering of the JT distortions in BaFeO 3 , and thus for high fundamental and technological interests.

Using first-principles calculations, we analyze the structural properties of tungsten trioxide WO3 . Our calculations rely on density functional theory and the use of the B1-WC hybrid functional, which provides very good agreement with experimental data. The hypothetical high-symmetry cubic reference structure combines several ferroelectric and antiferrodistortive (antipolar cation motions, rotations and tilts of oxygen octahedra) structural instabilities. The instability related to antipolar W motions combines with those associated to oxygen rotations and tilts to produce the biggest energy reduction, yielding a P21/c ground state. This non-polar P21/c phase is only different from the experimentally reported Pc ground state by the absence of a very tiny additional ferroelectric distortion. The calculations performed on a stoichiometric compound so suggest that the low temperature phase of WO3 is not intrinsically ferroelectric and that the reported ferroelectric character might arise from extrinsic defects such as oxygen vacancies. Independently, we also identify never observed R3m and R3c ferroelectric phases with large polarizations and low energies close to the P 21/c ground state, which makes WO3 a potential antiferroelectric material. The relative stability of various phases is discussed in terms of the couplings between different structural distortions, highlighting a very complex interplay involving improper-like couplings up to fourth order in the energy expansion in the cubic phase., Comment: 12 pages, 7 figures

Using first-principles calculations, we estimated the impact of large applied electric E fields on the structural, dielectric, and ferroelectric properties of typical ferroelectrics. At large fields, the structural parameters change significantly, decreasing the strain between the different structural phases. This effect favours a polarization rotation model for ferroelectric switching in which the electronic polarization rotates between the directions of tetragonal, rhombohedral and orthorhombic phases. We estimate coercive fields E_c ~31 MV/m and ~52 MV/m at zero temperature for bulk ferroelectric monodomains of BaTiO3 and PbTiO3, respectively. The dielectric permittivity and tunability of BaTiO3 are the least affected at large fields, making this material attractive for applications in electronics and energy storage., Comment: 6 pages, 6 figures, and 1 table

The amazing possibility of magnetoelectric crystals to cross couple electric and magnetic properties without the need of time-dependent Maxwell's equations has attracted a lot of interest in material science. This enthusiasm has re-emerged during the last decade where magnetoelectric and multiferroic crystals have captivated a tremendous number of studies, mostly driven by the quest of low-power-consumption spintronic devices. While several new candidates have been discovered, the desirable magnetoelectric coupling at room temperature is still sparse and calls for new promising candidates. Here, we show from first-principles studies that the troilite phase of the iron sulfide based compounds, one of the most common mineral of Earth, Moon, Mars or meteors, is magnetoelectric up to temperatures as high as 415 K.

Interest in first-principles calculations within the multiferroic community has been rapidly on the rise over the last decade. Initially considered as a powerful support to explain experimentally observed behaviours, the trend has evolved and, nowadays, density functional theory calculations has become also an essential predicting tool for identifying original rules to achieve multiferroism and design new magneto-electric compounds. This chapter aims to highlight the key advances in the field of multiferroics for which first-principles methods have contributed significantly. The essential theoretical developments that made this search possible are also briefly presented.

Perovskite oxides are already widely used in industry and have huge potential for novel device applications thanks to the rich physical behaviour displayed in these materials. The key to the functional electronic properties exhibited by perovskites is often the so-called Jahn-Teller distortion. For applications, an electrical control of the Jahn-Teller distortions, which is so far out of reach, would therefore be highly desirable. Based on universal symmetry arguments, we determine new lattice mode couplings that can provide exactly this paradigm, and exemplify the effect from first-principles calculations. The proposed mechanism is completely general, however for illustrative purposes, we demonstrate the concept on vanadium based perovskites where we reveal an unprecedented orbital ordering and Jahn-Teller induced ferroelectricity. Thanks to the intimate coupling between Jahn-Teller distortions and electronic degrees of freedom, the electric field control of Jahn-Teller distortions is of general relevance and may find broad interest in various functional devices., Comment: The first two authors contributed equally

First-principles calculations reveal a large cooperative coupling of Jahn-Teller (JT) distortion to oxygen-octahedron rotations in perovskite LaMnO$_3$. The combination of the two distortions is responsible for stabilizing the strongly orthorhombic $A$-AFM insulating ($I$) $Pbnm$ ground state relative to a metallic ferromagnetic (FM-$M$) phase. However, epitaxial strain due to coherent matching to a crystalline substrate can change the relative stability of the two states. In particular, coherent matching to a square-lattice substrate favors the less orthorhombic FM-$M$ phase, with the $A$-AFM phase stabilized at higher values of tensile epitaxial strain due to its larger volume per formula unit, resulting in a coupled magnetic and metal-insulator transition at a critical strain close to 1%. At the phase boundary, colossal magneto-resistance is expected. Tensile epitaxial strain enhances the JT distortion and opens the band gap in the $A$-AFM-$I$ $c$-$Pbnm$ phase, offering the opportunity for band-gap engineering. Compressive epitaxial strain induces an orientational transition within the FM-$M$ phase from $c$-$Pbnm$ to $ab$-$Pbnm$ with a change in the direction of the magnetic easy axis relative to the substrate, yielding strain-controlled magnetization at the phase boundary. The strong couplings between the JT distortion, the oxygen-octahedron rotations and the magnetic and electronic properties, and associated functional behavior, motivate interest in other orthorhombic $Pbnm$ perovskites with large JT distortions, which should also exhibit a rich variety of coupled magnetic, structural and electronic phase transitions driven by epitaxial strain., Comment: 5 figures

We develop the formalism of the macroscopic monopolization -- that is the monopole moment per unit volume -- in periodic solids, and discuss its relationship to the diagonal magnetoelectric effect. For the series of lithium transition metal phosphate compounds we use first-principles density functional theory to calculate the contributions to the macroscopic monopolization from the global distribution of magnetic moments within the unit cell, as well as from the distribution of magnetization around the atomic sites. We find one example within the series (LiMnPO$_4$) that shows a macroscopic monopolization corresponding to a ferromonopolar ordering consistent with its diagonal magnetoelectric response. The other members of the series (LiMPO$_4$, with M = Co, Fe and Ni) have zero net monopolization but have antiferromonopolar orderings that should lead to $q$-dependent diagonal magnetoelectric effects

We use symmetry analysis and first principles calculations to show that the linear magnetoelectric effect can originate from the response of orbital magnetic moments to the polar distortions induced by an applied electric field. Using LiFePO4 as a model compound we show that spin-orbit coupling partially lifts the quenching of the 3d orbitals and causes small orbital magnetic moments ($\mu_{(L)}\approx 0.3 \mu_B$) parallel to the spins of the Fe$^{2+}$ ions. An applied electric field $\mathbf{E}$ modifies the size of these orbital magnetic moments inducing a net magnetization linear in $\mathbf{E}$., Comment: 5 pages and 2 figures, plus 3 pages and 1 figure of supplementary material; 1 section added in supplementary material, 1 additional author, 1 reference updated

The link between the crystal distortions of the perovskite structure and the magnetic exchange interaction, the single-ion anisotropy (SIA) and the Dzyaloshinsky-Moriya (DM) interaction are investigated by means of density-functional calculations. Using BiFeO$_3$ and LaFeO$_3$ as model systems, we quantify the relationship between the oxygen octahedra rotations, the ferroelectricity and the weak ferromagnetism (wFM). We recover the fact that the wFM is due to the DM interaction induced by the oxygen octahedra rotations. We find a simple relationship between the wFM, the oxygen rotation amplitude and the ratio between the DM vector and the exchange parameter such as the wFM increases with the oxygen octahedra rotation when the SIA does not compete with the DM forces induced on the spins. Unexpectedly, we also find that, in spite of the $d^5$ electronic configuration of Fe$^{3+}$, the SIA is very large in some structures and is surprisingly strongly sensitive to the chemistry of the $A$-site cation of the $A$BO$_3$ perovskite. In the ground $R3c$ state phase we show that the SIA shape induced by the ferroelectricity and the oxygen octahedra rotations are in competition such as it is possible to tune the wFM "on" and "off" through the relative size of the two types of distortion.

In order to better understand the reconstructive ferroelectric-paraelectric transition of Bi2WO6, which is unusual within the Aurivillius family of compounds, we performed first principles calculations of the dielectric and dynamical properties on two possible high-temperature paraelectic structures: the monoclinic phase of A2/m symmetry observed experimentally and the tetragonal phase of I4/mmm symmetry, common to most Aurivillius phase components. Both paraelectric structures exhibits various unstable modes, which after their condensation bring the system toward more stable structures of lower symmetry. The calculations confirms that, starting from the paraelectric A2/m phase at high temperature, the system must undergo a reconstructive transition to reach the P2_1ab ferroelectric ground state., Comment: added Appendix and two tables

The antiferrodistortive (AFD) structural transitions of calcium titanate (CaTiO3) at ambient pressure have been extensively studied during the last few years. It is found none of the AFD polymorphs is polar or ferroelectric. However, it was recently shown theoretically and later experimentally confirmed that a ferroelectric transition in CaTiO3 can be induced by tensile strains. The ferroelectric instability is believed to be strongly coupled to the AFD soft modes. In this article, we present a complete thermodynamic potential for describing the coupling between the AFD and ferroelectric phase transitions. We analyzed the dependence of transition temperatures on stress and strain condition. Based on this potential, a (001) CaTiO3 thin film diagram was constructed. The results show good agreement with available experimental observations. The strong suppression of ferroelectric transition by the AFD transition is discussed.

We use symmetry analysis to show that the G, C and A-type antiferromagnetic Pnma perovskites can exhibit magnetoelectric (ME) responses when a ferroelectric instability is induced with epitaxial strain. Using first-principles calculations we compute the values of the allowed ME response in strained CaMnO3 as a model system. Our results show that large linear and non-linear ME responses are present and can diverge when close to the ferroelectric phase transition. By decomposing the electronic and ionic contributions, we explore the detailed mechanism of the ME response.

We re-examine the commonly used density functional theory plus Hubbard \textit{U} (DFT$+U$) method for the case of non-collinear magnets. While many studies neglect to explicitly include the exchange correction parameter \textit{J}, or consider its exact value to be unimportant, here we show that in the case of non-collinear magnetism calculations the \textit{J} parameter can strongly affect the magnetic ground state. We illustrate the strong \textit{J}-dependence of magnetic canting and magnetocrystalline anisotropy by calculating trends in the magnetic lithium orthophosphate family LiMPO$_4$ (M = Fe and Ni) and difluorite family MF$_2$ (M = Mn, Fe, Co and Ni). Our results can be readily understood by expanding the usual DFT$+U$ equations within the spinor scheme, in which the \textit{J} parameter acts directly on the off-diagonal components which determine the spin canting.

We report a first-principles study of (BaTiO${_3}$)$_m$/(BaO)$_n$ superlattices for a wide range of periodicities $m/n$. We show that such a system develops a polar zone-center instability for sufficiently large \textit{m}/\textit{n} ratio, which can be understood, at least qualitatively, from a simple electrostatic model and should lead to a ferroelectric ground-state. However, the analysis of the phonon dispersion curves also points out the appearance of stronger antiferroelectric instabilities at the zone boundaries around $m=4$, before the critical ratio for ferroelectricity is reached and which still dominate beyond it. The dominant character of the anti-ferroelectric instability is explained from the depolarizing field which hardens the ferroelectric mode. This analysis allows us to predict that, (BaTiO${_3}$)$_m$/(BaO)$_n$ superlattices should present an antiferroelectric ground state for $m$ larger than 4, which should smoothly evolve to a multidomain structure for increasing $m$ values and only become ferroelectric for large $m$.

We demonstrate that the electronic contribution to the linear magnetoelectric response, usually omitted in first-principles studies, can be comparable in magnitude to that mediated by lattice distortions, even for materials in which responses are strong. Using a self-consistent Zeeman response to an applied magnetic field for noncollinear electron spins, we show how electric polarization emerges in linear magnetoelectrics through both electronic- and lattice-mediated components -- in analogy with the high- and low-frequency dielectric response to an electric field. The approach we use is conceptually and computationally simple, and can be applied to study both linear and non-linear responses to magnetic fields., Comment: 5 pages, 3 figures

The alkaline earth binary oxides adopt a simple rocksalt structure and form an important family of compounds because of their large presence in the earth's mantle and their potential use in microelectronic devices. In comparison to the class of multifunctional ferroelectric perovskite oxides, however, their practical applications remain limited and the emergence of ferroelectricity and related functional properties in simple binary oxides seems so unlikely that it was never previously considered. Here, we show using first-principles density functional calculations that ferroelectricity can be easily induced in simple alkaline earth binary oxides such as barium oxide (BaO) using appropriate epitaxial strains. Going beyond the fundamental discovery, we show that the functional properties (polarization, dielectric constant and piezoelectric response) of such strained binary oxides are comparable in magnitude to those of typical ferroelectric perovskite oxides, so making them of direct interest for applications. Finally, we show that magnetic binary oxides such as EuO, with the same rocksalt structure, behave similarly to the alkaline earth oxides, suggesting a route to new multiferroics combining ferroelectric and magnetic properties.

From first-principles calculations, we investigate the structural instabilities of CaMnO$_3$. We point out that, on top of a strong antiferrodistortive instability responsible for its orthorhombic ground-state, the cubic perovskite structure of CaMnO$_3$ also exhibit a weak ferroelectric instability. Although ferroelectricity is suppressed by antiferrodistortive oxygen motions, we show that it can be favored using strain or chemical engineering in order to make CaMnO$_3$ multiferroic. We finally highlight that the FE instability of CaMnO$_3$ is Mn-dominated. This illustrates that, contrary to the common believe, ferroelectricity and magnetism are not necessarily exclusive but can be driven by the same cation.

The structural, dielectric and dynamical properties of the low temperature antiferromagnetic orthorhombic phase of CaMnO$_3$ have been computed from first principles using a density functional theory approach within the local spin density approximation. The theoretical structural parameters are in good agreement with experiment. The full set of infrared and Raman zone-center phonons is reported and compared to experimental data. It is shown that coherently with the anomalous Born effective charges and the presence of low frequency polar modes, the static dielectric constant is very large and highly anisotropic., Comment: Sumbitted to Phys. Rev. B

The size dependence of the ferroelectric properties of BaTiO$_3$ nanowires is studied from first-principles. We show that the ferroelectric distortion along the wire axis disappears below a critical diameter of about 1.2 nm. This disappearance is related to a global contraction of the unit cell resulting from low atomic coordinations at the wire surface. It is shown that a ferroelectric distortion can be recovered under appropriate tensile strain conditions., Comment: 3 pages, 3 figures. Accepted in Applied Physics Letters