Source author record

Nicola A. Spaldin

Nicola A. Spaldin appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

53works
8topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

53 published item(s)

preprint2023arXiv

Chiral phonons probed by X rays

The concept of chirality is of great relevance in nature, from chiral molecules such as sugar to parity transformations in particle physics. In condensed matter physics, recent studies have demonstrated chiral fermions and their relevance in emergent phenomena closely related to topology. The experimental verification of chiral phonons (bosons) remains challenging, however, despite their expected strong impact on fundamental physical properties. Here we show experimental proof of chiral phonons using resonant inelastic X-ray scattering with circularly polarized X rays. Using the prototypical chiral material, quartz, we demonstrate that circularly polarized X rays, which are intrinsically chiral, couple to chiral phonons at specific positions in reciprocal space, allowing us to determine the chiral dispersion of the lattice modes. Our experimental proof of chiral phonons demonstrates a new degree of freedom in condensed matter that is both of fundamental importance and opens the door to exploration of novel emergent phenomena based on chiral bosons

preprint2022arXiv

Accelerated search for new ferroelectric materials

We report the development of a combined machine-learning and high-throughput density functional theory (DFT) framework to accelerate the search for new ferroelectric materials. The framework can predict potential ferroelectric compounds using only elemental composition as input. A series of machine-learning algorithms initially predict the possible stable and insulating stoichiometries with a polar crystal structure, necessary for ferroelectricity, within a given chemical composition space. A classification model then predicts the point groups of these stoichiometries. A subsequent series of high-throughput DFT calculations finds the lowest energy crystal structure within the point group. As a final step, non-polar parent structures are identified using group theory considerations, and the values of the spontaneous polarization are calculated using DFT. By predicting the crystal structures and the polarization values, this method provides a powerful tool to explore new ferroelectric materials beyond those in existing databases.

preprint2022arXiv

Charge multipoles correlations and order in Cs$_2$TaCl$_6$

We examine the role of charge, structural, and spin degrees of freedom in the previously poorly understood phase transition in the 5$d^1$ transition-metal double perovskite Cs$_2$TaCl$_6$, using a combination of computational and experimental techniques. Our heat capacity measurements of single-crystalline Cs$_2$TaCl$_6$, reveal a clear anomaly at the transition temperature, $T_\mathrm{Q}$, which was not previously observed in polycrystalline samples. Density functional calculations indicate the emergence of local charge quadrupoles in the cubic phase, mediated by the paramagnetic spins or local structural distortions which then develop into long-range ordered charge quadrupoles in the tetragonal phase. Our resonant elastic x-ray scattering on Cs$_2$TaCl$_6$, single crystals lend support to our calculations. Our work provides new insight into the phase transition in Cs$_2$TaCl$_6$, at $T_\mathrm{Q}$, and demonstrates the utility of this combination of techniques in understanding the complex physics of hidden orders in paramagnetic spin-orbit-entangled compounds.

preprint2022arXiv

Crystal responses to general dark matter-electron interactions

We develop a formalism to describe the scattering of dark matter (DM) particles by electrons bound in crystals for a general form of the underlying DM-electron interaction. Such a description is relevant for direct-detection experiments of DM particles lighter than a nucleon, which might be observed in operating DM experiments via electron excitations in semiconductor crystal detectors. Our formalism is based on an effective theory approach to general non-relativistic DM-electron interactions, including the anapole, and magnetic and electric dipole couplings, combined with crystal response functions defined in terms of electron wave function overlap integrals. Our main finding is that, for the usual simplification of the velocity integral, the rate of DM-induced electronic transitions in a semiconductor material depends on at most five independent crystal response functions, four of which were not known previously. We identify these crystal responses, and evaluate them using density functional theory for crystalline silicon and germanium, which are used in operating DM direct detection experiments. Our calculations allow us to set 90% confidence level limits on the strength of DM-electron interactions from data reported by the SENSEI and EDELWEISS experiments. The novel crystal response functions discovered in this work encode properties of crystalline solids that do not interact with conventional experimental probes, suggesting the use of the DM wind as a probe to reveal new kinds of hidden order in materials.

preprint2022arXiv

Influence of the triangular Mn-O breathing mode on the magnetic ordering in multiferroic hexagonal manganites

We use a combination of symmetry analysis, phenomenological modelling and first-principles density functional theory to explore the interplay between the magnetic ground state and the detailed atomic structure in the hexagonal rare-earth manganites. We find that the magnetic ordering is sensitive to a breathing mode distortion of the Mn and O ions in the $ab$ plane, which is described by the K\textsubscript{1} mode of the high-symmetry structure. Our density functional calculations of the magnetic interactions indicate that this mode particularly affects the single-ion anisotropy and the inter-planar symmetric exchanges. By extracting the parameters of a magnetic model Hamiltonian from our first-principles results, we develop a phase diagram to describe the magnetic structure as a function of the anisotropy and exchange interactions. This in turn allows us to explain the dependence of the magnetic ground state on the identity of the rare-earth ion and on the K\textsubscript{1} mode.

preprint2022arXiv

Leggett Modes Accompanying Crystallographic Phase Transitions

Higgs and Goldstone modes, well known in high energy physics, have been realized in a number of condensed matter physics contexts, including superconductivity and magnetism. The Goldstone-Higgs concept is also applicable to and gives rise to new insights on structural phase transitions. Here, we show that the Leggett mode, a collective mode observed in multi-band superconductors, also has an analog in crystallographic phase transitions. Such structural Leggett modes can occur in the phase channel as in the original work of Leggett, \href{https://doi.org/10.1143/PTP.36.901}{Prog.\ Theor.\ Phys.\ \textbf{36}, 901 (1966)}. That is, they are antiphase Goldstone modes (anti-phasons). In addition, a new collective mode can also occur in the amplitude channel, an out-of phase (antiphase) Higgs mode, that should be observable in multi-band superconductors as well. We illustrate the existence and properties of these structural Leggett modes using the example of the pyrochlore relaxor ferroelectric, Cd$_2$Nb$_2$O$_7$.

preprint2022arXiv

Local structure and its implications for the relaxor ferroelectric Cd$_2$Nb$_2$O$_7$

The relaxor ferroelectric transition in Cd$_2$Nb$_2$O$_7$ is thought to be described by the unusual condensation of two $Γ$-centered phonon modes, $Γ_4^-$ and $Γ_5^-$. However, their respective roles have proven to be ambiguous, with disagreement between $\textit{ab initio}$ studies, which favor $Γ_4^-$ as the primary mode, and global crystal refinements, which point to $Γ_5^-$ instead. Here, we resolve this issue by demonstrating from x-ray pair distribution function measurements that locally, $Γ_4^-$ dominates, but globally, $Γ_5^-$ dominates. This behavior is consistent with the near degeneracy of the energy surfaces associated with these two distortion modes found in our own $\textit{ab initio}$ simulations. Our first-principles calculations also show that these energy surfaces are almost isotropic, providing an explanation for the numerous structural transitions found in Cd$_2$Nb$_2$O$_7$, as well as its relaxor behavior. Our results point to several candidate descriptions of the local structure, some of which demonstrate two-in/two-out behavior for Nb displacements within a given Nb tetrahedron. Although this suggests the possibility of a charge analog of spin ice in Cd$_2$Nb$_2$O$_7$, our results are more consistent with a Heisenberg-like description for dipolar fluctuations rather than an Ising one. We hope this encourages future experimental investigations of the Nb and Cd dipolar fluctuations, along with their associated mode dynamics.

preprint2022arXiv

Magnetoelectric classification of skyrmions

We develop a general theory to classify magnetic skyrmions and related spin textures in terms of their magnetoelectric multipoles. Since magnetic skyrmions are now established in insulating materials, where the magnetoelectric multipoles govern the linear magnetoelectric response, our classification provides a recipe for manipulating the magnetic properties of skyrmions using applied electric fields. We apply our formalism to skyrmions and anti-skyrmions of different helicities, as well as to magnetic bimerons, which are topologically, but not geometrically, equivalent to skyrmions. We show that the non-zero components of the magnetoelectric multipole and magnetoelectric response tensors are uniquely determined by the topology, helicity and geometry of the spin texture. Therefore, we propose straightforward linear magnetoelectric response measurements as an alternative to Lorentz microscopy for characterizing insulating skyrmionic textures.

preprint2022arXiv

Magnetoelectric effect in hydrogen harvesting: magnetic field as a trigger of catalytic reactions

Magnetic fields have been regarded as an additional stimulus for electro- and photocatalytic reactions, but not as a direct trigger for catalytic processes. Multiferroic/magnetoelectric materials, whose electrical polarization and surface charges can be magnetically altered, are especially suitable for triggering and control of catalytic reactions solely with magnetic fields. Here, we demonstrate that magnetic fields can be employed as an independent input energy source for hydrogen harvesting by means of the magnetoelectric effect. Composite multiferroic CoFe2O4-BiFeO3 core-shell nanoparticles act as catalysts for the hydrogen evolution reaction (HER) that is triggered when an alternating magnetic field is applied to an aqueous dispersion of the magnetoelectric nanocatalysts. Based on density functional calculations, we propose that the hydrogen evolution is driven by changes in the ferroelectric polarization direction of BiFeO3 caused by the magnetoelectric coupling. We believe our findings will open new avenues towards magnetically induced renewable energy harvesting.

preprint2022arXiv

Nonvolatile Electric-Field Control of Inversion Symmetry

In condensed-matter systems, competition between ground states at phase boundaries can lead to significant changes in material properties under external stimuli, particularly when these ground states have different crystal symmetries. A key scientific and technological challenge is to stabilize and control coexistence of symmetry-distinct phases with external stimuli. Using BiFeO3 (BFO) layers confined between layers of the dielectric TbScO3 as a model system, we stabilize the mixed-phase coexistence of centrosymmetric and non-centrosymmetric BFO phases with antipolar, insulating and polar, semiconducting behavior, respectively at room temperature. Application of in-plane electric (polar) fields can both remove and introduce centrosymmetry from the system resulting in reversible, nonvolatile interconversion between the two phases. This interconversion between the centrosymmetric insulating and non-centrosymmetric semiconducting phases coincides with simultaneous changes in the non-linear optical response of over three orders of magnitude, a change in resistivity of over five orders of magnitude, and a change in the polar order. Our work establishes a materials platform allowing for novel cross-functional devices which take advantage of changes in optical, electrical, and ferroic responses.

preprint2021arXiv

Ferroelectric, quantum paraelectric or paraelectric? Calculating the evolution from BaTiO$_3$ to SrTiO$_3$ to KTaO$_3$ using a single-particle quantum-mechanical description of the ions

We present an inexpensive first-principles approach for describing quantum paraelectricity that combines density functional theory (DFT) treatment of the electronic subsystem with quantum mechanical treatment of the ions through solution of the single-particle Schrödinger equation with the DFT-calculated potential. Using BaTiO$_3$, SrTiO$_3$ and KTaO$_3$ as model systems, we show that the approach can straightforwardly distinguish between ferroelectric, paraelectric and quantum paraelectric materials, based on simple quantities extracted from standard density functional and density functional perturbation theories. We calculate the influence of isotope substitution and strain on the quantum paraelectric behavior, and find that, while complete replacement of oxygen-16 by oxygen-18 has a surprisingly small effect, experimentally accessible strains can induce large changes. Finally, we collect the various choices for the phonon mass that have been introduced in the literature. We identify those that are most physically meaningful by comparing with our results that avoid such a choice through the use of mass-weighted coordinates.

preprint2020arXiv

Atomic responses to general dark matter-electron interactions

In the leading paradigm of modern cosmology, about 80% of our Universe's matter content is in the form of hypothetical, as yet undetected particles. These do not emit or absorb radiation at any observable wavelengths, and therefore constitute the so-called Dark Matter (DM) component of the Universe. Detecting the particles forming the Milky Way DM component is one of the main challenges for astroparticle physics and basic science in general. One promising way to achieve this goal is to search for rare DM-electron interactions in low-background deep underground detectors. Key to the interpretation of this search is the response of detectors' materials to elementary DM-electron interactions defined in terms of electron wave functions' overlap integrals. In this work, we compute the response of atomic argon and xenon targets used in operating DM search experiments to general, so far unexplored DM-electron interactions. We find that the rate at which atoms can be ionized via DM-electron scattering can in general be expressed in terms of four independent atomic responses, three of which we identify here for the first time. We find our new atomic responses to be numerically important in a variety of cases, which we identify and investigate thoroughly using effective theory methods. We then use our atomic responses to set 90% confidence level (C.L.) exclusion limits on the strength of a wide range of DM-electron interactions from the null result of DM search experiments using argon and xenon targets.

preprint2020arXiv

Comparison of coherent phonon generation by electronic and ionic Raman scattering in LaAlO$_3$

In ionic Raman scattering, infrared-active phonons mediate a scattering process that results in the creation or destruction of a Raman-active phonon. This mechanism relies on nonlinear interactions between phonons and has in recent years been associated with a variety of emergent lattice-driven phenomena in complex transition-metal oxides, but the underlying mechanism is often obscured by the presence of multiple coupled order parameters in play. Here, we use time-resolved spectroscopy to compare coherent phonons generated by ionic Raman scattering with those created by more conventional electronic Raman scattering on the nonmagnetic and non-strongly-correlated wide band-gap insulator LaAlO$_3$. We find that the oscillatory amplitude of the low-frequency Raman-active $E_g$ mode exhibits a sharp peak when we tune our pump frequency into resonance with the high-frequency infrared-active $E_u$ mode, consistent with first-principles calculations. Our results suggest that ionic Raman scattering can strongly dominate electronic Raman scattering in wide band-gap insulating materials. We also see evidence of competing scattering channels at fluences above 28~mJ/cm$^2$ that alter the measured amplitude of the coherent phonon response.

preprint2020arXiv

On the happiness of ferroelectric surfaces and its role in water dissociation: the example of bismuth ferrite

We investigate, using density functional theory, how the interaction between the ferroelectric polarization and the chemical structure of the (001) surfaces of bismuth ferrite influences the surface properties and reactivity of this material. A precise understanding ofthe surface behavior of ferroelectrics is necessary for their use in surface science applications such as catalysis as well as for their incorporation in microelectronic devices. Using the (001) surface of bismuth ferrite as a model system we show that the most energetically favoured surface geometries are combinations of surface termination and polarization direction that lead to uncharged, stable surfaces. On the unfavorable charged surfaces, we explore the compensation mechanisms of surface charges provided by the introduction of point defects and adsorbates, such as water. Finally, we propose that the special surface properties of bismuth ferrite (001) could be used to produce an effective water splitting cycle through cyclic polarization switching.

preprint2019arXiv

Origin and evolution of ferroelectricity in the layered rare-earth-titanate, $R_2$Ti$_2$O$_{7}$, Lichtenberg phases

We report a systematic first-principles study based on density functional theory (DFT) of the structural and ferroelectric properties of the $R_2$Ti$_2$O$_{7}$ perovskite-related oxides with $R=$ La, Ce, Pr, and Nd. We show that, in all cases, the ferroelectric distortion to the ground-state polar $P2_1$ structure from its parent centrosymmetric $P2_1/m$ phase is driven by a single polar soft mode consisting of rotations and tilts of the TiO$_6$ octahedra combined with displacements of the $R$ ions. A secondary centrosymmetric distortion, which is stable in the parent structure, contributes substantially to the energy lowering of the ground state. We evaluate the trends in structure and polarization across the series as a function of $R$ and reconcile discrepancies in reported values of polarization in the literature. Our results confirm that the family of $R_2$Ti$_2$O$_{7}$ materials belong to the class of proper geometric ferroelectrics.

preprint2018arXiv

Magnetoelectric multipoles in metals

We demonstrate computationally the existence of magnetoelectric multipoles, arising from the second order term in the multipole expansion of a magnetization density in a magnetic field, in noncentrosymmetric magnetic metals. While magnetoelectric multipoles have long been discussed in the context of the magnetoelectric effect in noncentrosymmetric magnetic insulators, they have not previosuly been identified in metallic systems, in which the mobile carriers screen any electrical polarization. Using first-principles density functional calculations we explore three specific systems: First, a conventional centrosymmetric magnetic metal, Fe, in which we break inversion symmetry by introducing a surface, which both generates magnetoelectric monopoles and allows a perpendicular magnetoelectric response. Next, the hypothetical cation-ordered perovskite, SrCaRu$_2$O$_6$, in which we study the interplay between the magnitude of the polar symmetry breaking and that of the magnetic dipoles and multipoles, finding that both scale proportionally to the structural distortion. Finally, we identify a hidden antiferromultipolar order in the noncentrosymmetric, antiferromagnetic metal Ca$_3$Ru$_2$O$_7$, and show that, while its competing magnetic phases have similar magnetic dipolar structures, their magnetoelectric multipolar structures are distinctly different, reflecting the strong differences in transport properties.

preprint2017arXiv

On the relationship between topological and geometric defects

The study of topology in solids is undergoing a renaissance following renewed interest in the properties of ferroic domain walls as well as recent discoveries regarding topological insulators and skyrmionic lattices. Each of these systems possess a property that is `protected' in a symmetry sense, and is defined rigorously using a branch of mathematics known as topology. In this article we review the formal definition of topological defects as they are classified in terms of homotopy theory, and discuss the precise symmetry-breaking conditions that lead to their formation. We distinguish topological defects from geometric defects, which arise from the details of the stacking or structure of the material but are not protected by symmetry. We provide simple material examples of both topological and geometric defect types, and discuss the implications of the classification on the resulting material properties.

preprint2016arXiv

Coupling and competition between ferroelectricity, magnetism, strain and oxygen vacancies in AMnO3 perovskites

We use first-principles calculations based on density functional theory to investigate the interplay between oxygen vacancies, A-site cation size / tolerance factor, epitaxial strain, ferroelectricity and magnetism in the perovskite manganite series, AMnO3 (A=Ca2+, Sr2+, Ba2+). We find that, as expected, increasing the volume through either chemical pressure or tensile strain generally lowers the formation energy of neutral oxygen vacancies consistent with their established tendency to expand the lattice. Increased volume also favors polar distortions, both because competing rotations of the oxygen octahedra are suppressed and because coulomb repulsion associated with cation off-centering is reduced. Interestingly, the presence of ferroelectric polarization favors ferromagnetic over antiferromagnetic ordering due to suppressed antiferromagnetic superexchange as the polar distortion bends the Mn-O-Mn bond angles away from the optimal 180deg. Intriguingly, we find that polar distortions compete with the formation of oxygen vacancies, which have a higher formation energy in the polar phases; conversely the presence of oxygen vacancies suppresses the onset of polarization. In contrast, oxygen vacancy formation energies are lower for ferromagnetic than antiferromagnetic orderings of the same structure type. Our findings suggest a rich and complex phase diagram, in which defect chemistry, polarization, structure and magnetism can be modified using chemical potential, strain, and electric or magnetic fields.

preprint2016arXiv

Effects of intense optical phonon pumping on the structure and electronic properties of yttrium barium copper oxide

We investigate the structural modulations induced by optical excitation of a polar phonon mode in YBa$_2$Cu$_3$O$_7$, using first-principles calculations based on density functional theory. We focus on the intense-excitation regime in which we expect that fourth-order phonon-phonon coupling terms dominate, and model the structural modulations induced by pulses of such intensity. Our calculations of the phonon-phonon anharmonicities confirm that the cubic coupling between modes, shown in earlier work to cause a quasi-static change in the apical O - Cu distance and a buckling of the CuO$_2$ planes, is the leading contribution at moderate pump strengths. At higher pump strengths ($\sim$10 MV/cm) the previously neglected quartic couplings become relevant and produce an additional shearing of the CuO$_2$ planes. Finally, we analyze the changes in the electronic and magnetic properties associated with the induced structural changes.

preprint2016arXiv

First-principles calculation of the bulk magnetoelectric monopolization: Berry phase and Wannier function approaches

We present a formalism to calculate the macroscopic magnetoelectric monopolization from first principles within the density functional theory framework. An expression for the monopolization in the case of insulating collinear magnetism is derived first in terms of spin-polarized Wannier functions then recast as a Berry phase. We propose an extension to the general, non-collinear case which we implement computationally in the Wannier function form and use to calculate the magnetoelectric monopolizations of LiMnPO$_4$ and Cr$_2$O$_3$. We find that, while the former is well approximated by a summation over the formal local spin moments, the latter shows significant deviations from this approximation. We suggest that equating the Berry phase value with a sum over local moments provides an unambiguous route to defining the size of the local magnetic moment in magnetoelectric antiferromagnets containing only one type of magnetic ion.

preprint2016arXiv

Interplay between strain, defect charge state and functionality in complex oxides

We use first-principles calculations to investigate the interplay between strain and the charge state of point defect impurities in complex oxides. Our work is motivated by recent interest in using defects as active elements to provide novel functionality in coherent epitaxial films. Using oxygen vacancies as model point defects, and CaMnO$_3$ and MnO as model materials, we calculate the changes in internal strain caused by changing the charge state of the vacancies, and conversely the effect of strain on charge-state stability. Our results show that the charge state is a degree of freedom that can be used to control the interaction of defects with strain and hence the concentration and location of defects in epitaxial films. We propose the use of field-effect gating to reversibly change the charge state of defects and hence the internal strain and corresponding strain-induced functionalities.

preprint2016arXiv

Strain-induced structural instability in FeRh

We perform density functional calculations to investigate the structure of the inter-metallic alloy FeRh under epitaxial strain. Bulk FeRh exhibits a metamagnetic transition from a low-temperature antiferromagnetic (AFM) phase to a ferromagnetic (FM) phase at 350K, and its strain dependence is of interest for tuning the transition temperature to the room-temperature operating conditions of typical memory devices. We find an unusually strong dependence of the structural energetics on the choice of exchange-correlation functional, with the usual local density approximation (LDA) yielding the wrong ground-state structure, and generalized gradient (GGA) extensions being in better agreement with the bulk experimental structure. Using the GGA we show the existence of a metastable face-centered-cubic (fcc)-like AFM structure that is reached from the ground state body-centered-cubic (bcc) AFM structure by following the epitaxial Bain path. We predict that this metastable fcc-like structure has a significantly higher conductivity than the bcc AFM phase. We show that the behavior is well described using non-linear elasticity theory, which captures the softening and eventual sign change of the orthorhombic shear modulus under compressive strain, consistent with this structural instability. Finally, we predict the existence of an additional unit-cell-doubling lattice instability, which should be observable at low temperature.

preprint2015arXiv

Ab initio study of the ferroelectric strain dependence and 180$^\circ$ domain walls in the barium metal fluorides BaMgF$_4$ and BaZnF$_4$

We investigate the strain dependence of the ferroelectric polarization and the structure of the ferroelectric domain walls in the layered perovskite-related barium fluorides, BaMF$_4$ (M=Mg, Zn). The unusual "geometric ferroelectricity" in these materials is driven by the softening of a single polar phonon mode consisting of rotations of the MF$_6$ octahedra accompanied by polar displacements of the Ba cations, and in contrast to conventional ferroelectrics involves minimal electronic rehybridization. We therefore anticipate a different strain dependence of the polarization, and alternative domain wall structures compared with those found in conventional ferroelectric materials. Using first-principles calculations based on density functional theory (DFT) within the general gradient approximation (GGA), we calculate the variation of the crystal structure and the ferroelectric polarization under both compressive and tensile strain. We perform structural relaxations of neutral domain walls between oppositely oriented directions of the ferroelectric polarization and calculate their corresponding energies to determine which are most likely to form. We compare our results to literature values for conventional perovskite oxides to provide a source of comparison for understanding the ferroelectric properties of alternative non-oxide materials such as the barium fluorides.

preprint2015arXiv

Effect of epitaxial strain on cation and anion vacancy formation in MnO

Biaxial strain in coherent epitaxial thin films can have a pronounced effect on the point defect profile in the film material. Detailed fundamental knowledge of the interaction of strain with point defects is crucial in understanding the stoichiometry and resulting properties of strained thin films. Here we investigate the effect of biaxial strain on the formation energy of cation and anion vacancies using MnO as a model system. Our density functional theory calculations show that, as expected from local volume arguments, compressive strain favours the formation of cation vacancies. Interestingly, we find that small compressive and tensile strains lead to ordering of the resulting holes along the in-plane and normal direction respectively, which should manifest in different anisotropic properties in the two strain states.

preprint2015arXiv

Quantum critical origin of the superconducting dome in SrTiO$_3$

We investigate the origin of superconductivity in doped SrTiO$_3$ (STO) using a combination of density functional and strong coupling theories within the framework of quantum criticality. Our density functional calculations of the ferroelectric soft mode frequency as a function of doping reveal a crossover from quantum paraelectric to ferroelectric behavior at a doping level coincident with the experimentally observed top of the superconducting dome. Based on this finding, we explore a model in which the superconductivity in STO is enabled by its proximity to the ferroelectric quantum critical point and the soft mode fluctuations provide the pairing interaction on introduction of carriers. Within our model, the low doping limit of the superconducting dome is explained by the emergence of the Fermi surface, and the high doping limit by departure from the quantum critical regime. We predict that the highest critical temperature will increase and shift to lower carrier doping with increasing $^{18}$O isotope substitution, a scenario that is experimentally verifiable.

preprint2015arXiv

Strain-Induced Magnetic Anisotropy in Epitaxial Thin Films of the Spinel CoCr$_2$O$_4$

We show that the magnetic anisotropy in spinel-structure CoCr$_2$O$_4$ thin films exhibits a strain dependence in which compressive strain induces an out-of-plane magnetic easy axis and tensile strain an in-plane easy axis, exactly opposite to the behavior reported for the related compound CoFe$_2$O$_4$. We use density functional theory calculations within the LSDA+U approximation to reproduce and explain the observed behavior. Using second-order perturbation theory, we analyse the anisotropy tensor of the Co$^{2+}$ ions in both octahedral and tetrahedral coordination, allowing us to extend our results to spinels with general arrangements of Co$^{2+}$ ions.

preprint2014arXiv

Biquadratic and ring exchange interactions in orthorhombic perovskite manganites

We use ab initio electronic structure calculations within the generalized gradient approximation (GGA+U) to density functional theory (DFT) to determine the microscopic exchange interactions in the series of orthorhombic rare-earth manganites, o-$R$MnO$_3$. Our motivation is to construct a model Hamiltonian (excluding effects due to spin-orbit coupling), which can provide an accurate description of the magnetism in these materials. First, we consider TbMnO$_3$, which exhibits a spiral magnetic order at low temperatures. We map the exchange couplings in this compound onto a Heisenberg Hamiltonian and observe a clear deviation from the Heisenberg-like behavior. We consider first the coupling between magnetic and orbital degrees of freedom as a potential source of non-Heisenberg behavior in TbMnO$_3$, but conclude that it does not explain the observed deviation. We find that higher order magnetic interactions (biquadratic and four-spin ring couplings) should be taken into account for a proper treatment of the magnetism in TbMnO$_3$ as well as in the other representatives of the o-$R$MnO$_3$ series with small radii of the $R$ cation.

preprint2014arXiv

The magnetic field generated by a charge in a uniaxial magnetoelectric material

We revisit the description of the magnetic field around antiferromagnetic magnetoelectrics in the context of recent developments regarding magnetoelectric monopoles. Using Maxwell's equations, we calculate the magnetic and electric fields associated with a free charge in a bulk uniaxial magnetoelectric, as well as in a finite sphere of magnetoelectric material. We show that a charge in the prototypical magnetoelectric Cr$_2$O$_3$, which is uniaxial with a diagonal magnetoelectric response, induces an internal magnetic field with both monopolar and quadrupolar components, but that only the quadrupolar contribution extends beyond the sample surface. We discuss the behavior of the external quadrupolar field and compare its magnitude to those of magnetic fields from other sources.

preprint2013arXiv

Competition and cooperation between antiferrodistortive and ferroelectric instabilities in SrTiO$_3$

We use density functional theory to explore the interplay between octahedral rotations and ferroelectricity in the model compound SrTiO$_3$. We find that over the experimentally relevant range octahedral rotations suppress ferroelectricity as is generally assumed in the literature. Somewhat surprisingly we observe that at larger angles the previously weakened ferroelectric instability strengthens significantly. By analysing geometry changes, energetics, force constants and charges, we explain the mechanisms behind this transition from competition to cooperation with increasing octahedral rotation angle.

preprint2013arXiv

Strain controlled oxygen vacancy formation and ordering in CaMnO$_3$

We use first-principles calculations to investigate the stability of bi-axially strained \textit{Pnma} perovskite CaMnO$_3$ towards the formation of oxygen vacancies. Our motivation is provided by promising indications that novel material properties can be engineered by application of strain through coherent heteroepitaxy in thin films. While it is usually assumed that such epitaxial strain is accommodated primarily by changes in intrinsic lattice constants, point defect formation is also a likely strain relaxation mechanism. This is particularly true at the large strain magnitudes ($>$4%) which first-principles calculations often suggest are required to induce new functionalities. We find a strong dependence of oxygen vacancy defect formation energy on strain, with tensile strain lowering the formation energy consistent with the increasing molar volume with increasing oxygen deficiency. In addition, we find that strain differentiates the formation energy for different lattice sites, suggesting its use as a route to engineering vacancy ordering in epitaxial thin films.

preprint2013arXiv

Strong coupling of Jahn-Teller distortion to oxygen-octahedron rotation and functional properties in epitaxially-strained orthorhombic LaMnO$_3$

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.

preprint2013arXiv

The macroscopic monopolization in diagonal magnetoelectrics

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

preprint2012arXiv

First-principles prediction of oxygen octahedral rotations in perovskite-structure EuTiO3

We present a systematic first-principles study of the structural and vibrational properties of perovskite-structure EuTiO3. Our calculated phonon spectrum of the high-symmetry cubic structural prototype shows strong M- and R-point instabilities, indicating a tendency to symmetry-lowering structural deformations composed of rotations and tilts of the oxygen octahedra. Subsequent explicit study of 14 different octahedral tilt-patterns showed that the I4/mcm, Imma, and R\bar{3}c structures, all with antiferrodistortive rotations of the octahedra, have significantly lower total energy than the prototype Pm\bar{3}m structure. We discuss the dynamical stability of these structures, and the influence of the antiferrodistortive structural distortions on the vibrational, optical, and magnetic properties of EuTiO3, in the context of recent unexplained experimental observations.

preprint2012arXiv

Interplay between strain and oxygen vacancies in lanthanum aluminate

We evaluate the interplay between epitaxial strain and oxygen vacancy formation in the perovskite-structure oxide lanthanum aluminate, LaAlO$_3$. Using density functional theory within the GGA$+U$ approximation we calculate the dependence of the oxygen vacancy formation energy on the biaxial strain conditions. We find that the change in formation energy with strain is negligible over the range of strain values usually accessible through coherent epitaxial growth. Our findings suggest that experimental reports of different oxygen vacancy concentrations in strained films result from extrinsic factors, or from other impurities such as defect complexes.

preprint2012arXiv

Linear Magnetoelectric Effect by Orbital Magnetism

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 ($μ_{(L)}\approx 0.3 μ_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}$.

preprint2012arXiv

Observation of persistent centrosymmetricity in the hexagonal manganite family

The controversy regarding the ferroelectric behavior of hexagonal InMnO$_3$ is resolved by using a combination of x-ray diffraction (XRD), piezoresponse force microscopy (PFM), second harmonic generation (SHG), and density functional theory (DFT). While XRD data show a symmetry-lowering unit-cell tripling, which is also found in the multiferroic hexagonal manganites of $P6_3cm$ symmetry, PFM and SHG do not detect ferroelectricity at ambient or low temperature, in striking contrast to the behavior in the multiferroic counterparts. We propose instead a centrosymmetric $P\bar{3}c$ phase as the ground state structure. Our DFT calculations reveal that the relative energy of the ferroelectric and nonferroelectric structures is determined by a competition between electrostatics and oxygen-$R$-site covalency, with an absence of covalency favoring the ferroelectric phase.

preprint2012arXiv

Structural Domain Walls in Polar Hexagonal Manganites

We present a microscopic first-principles study of the neutral structural domain walls (DWs) in the multiferroic hexagonal manganites, which have been shown to exhibit cross-couplings between ferroelectricity and struc- tural antiphase. We find that, in contradiction with previously proposed models, the interlocked antiphase- ferroelectric domain walls have approximately zero width, and their energy is lower than that of antiphase-only or ferroelectric-only domain walls. Furthermore, we show that the ferroelectric-only and antiphase-only DWs are superpositions of interlocked antiphase ferroelectric (AP+FE) DWs and intermediate domains inevitably exist through the DWs. Our results shed light on the question of why only AP+FE DWs are observed and how the topological defects emerge in polar hexagonal manganites.

preprint2011arXiv

Ab initio investigation of FeAs/GaAs heterostructures for potential spintronic and superconducting applications

Ultra-thin FeAs is of interest both as the active component in the newly identified pnictide superconductors, and in spintronic applications at the interface between ferromagnetic Fe and semiconducting GaAs. Here we use first-principles density functional theory to investigate the properties of FeAs/GaAs heterostructures. We find that the Fermi surface is modified from that characteristic of the pnictide superconductors by interactions between the FeAs layer and the As atoms in the GaAs layers. Regardless of the number of FeAs layers, the Fe to As ratio, or the strain state, the lowest energy magnetic ordering is always antiferromagnetic, suggesting that such heterostructures are not promising spintronic systems, and offering an explanation for the failure of spin injection across Fe/GaAs interfaces.

preprint2011arXiv

Anisotropic conductance at improper ferroelectric domain walls

Transition metal oxides hold great potential for the development of new device paradigms because of the field-tunable functionalities driven by their strong electronic correlations, combined with their earth abundance and environmental friendliness. Recently, the interfaces between transition-metal oxides have revealed striking phenomena such as insulator-metal transitions, magnetism, magnetoresistance, and superconductivity. Such oxide interfaces are usually produced by sophisticated layer-by-layer growth techniques, which can yield high quality, epitaxial interfaces with almost monolayer control of atomic positions. The resulting interfaces, however, are fixed in space by the arrangement of the atoms. Here we demonstrate a route to overcoming this geometric limitation. We show that the electrical conductance at the interfacial ferroelectric domain walls in hexagonal ErMnO3 is a continuous function of the domain wall orientation, with a range of an order of magnitude. We explain the observed behaviour using first-principles density functional and phenomenological theories, and relate it to the unexpected stability of head-to-head and tail-to-tail domain walls in ErMnO3 and related hexagonal manganites. Since the domain wall orientation in ferroelectrics is tunable using modest external electric fields, our finding opens a degree of freedom that is not accessible to spatially fixed interfaces.

preprint2011arXiv

Band alignment at metal/ferroelectric interfaces: insights and artifacts from first principles

Based on recent advances in first-principles theory, we develop a general model of the band offset at metal/ferroelectric interfaces. We show that, depending on the polarization of the film, a pathological regime might occur where the metallic carriers populate the energy bands of the insulator, making it metallic. As the most common approximations of density functional theory are affected by a systematic underestimation of the fundamental band gap of insulators, this scenario is likely to be an artifact of the simulation. We provide a number of rigorous criteria, together with extensive practical examples, to systematically identify this problematic situation in the calculated electronic and structural properties of ferroelectric systems. We discuss our findings in the context of earlier literature studies, where the issues described in this work have often been overlooked. We also discuss formal analogies to the physics of polarity compensation at LaAlO3/SrTiO3 interfaces, and suggest promising avenues for future research.

preprint2011arXiv

High-temperature multiferroicity and strong magnetocrystalline anisotropy in 3d-5d double perovskites

Using density functional calculations we explore the properties of as-yet-unsynthesized $3d - 5d$ ordered double perovskites ($A_2BB'$O$_6$) with highly polarizable Bi$^{3+}$ ions on the $A$ site. We find that the Bi$_2$NiReO$_6$ and Bi$_2$MnReO$_6$ compounds are insulating and exhibit a robust net magnetization that persists above room temperature. When the in-plane lattice vectors of the pseudocubic unit cell are constrained to be orthogonal (for example, by coherent heteroepitaxy), the ground states are ferroelectric with large polarization and a very large uniaxial magnetocrystalline anisotropy with easy axis along the ferroelectric polarization direction. Our results suggest a route to multiferroism and electrically controlled magnetization orientation at room temperature.

preprint2010arXiv

Chemical control of polar behavior in bicomponent short-period superlattices

Using first-principles density functional calculations, we study the interplay of ferroelectricity and polar discontinuities in a range of 1-1 oxide superlattices, built out of ferroelectric and paraelectric components. Studies have been carried out for a varied choice of chemical composition of the components. We find that, when polar interfaces are present, the polar discontinuities induce off- centric movements in the ferroelectric layers, even though the ferroelectric is only one unit cell thick. The distortions yield non-switchable polarizations, with magnitudes comparable to those of the corresponding bulk ferroelectrics. In contrast, in superlattices with no polar discontinuity at the interfaces, the off-centric movements in the ferroelectric layer are usually suppressed. The details of the behavior and functional properties are, however, found to be sensitive to epitaxial strain, rotational instabilities and second-order Jahn-Teller activity, and are therefore strongly in uenced by the chemical composition of the paraelectric layer.

preprint2010arXiv

Electric and magnetic polarizabilities of hexagonal Ln2CuTiO6 (Ln=Y, Dy, Ho, Er and Yb)

We investigated the rare-earth transition metal oxide series, Ln2CuTiO6 (Ln=Y, Dy, Ho, Er and Yb), crystallizing in the hexagonal structure with non-centrosymmetric P63cm space group for possible occurrences of multiferroic properties. Our results show that while these compounds, except Ln=Y, exhibit a low temperature antiferromagnetic transition due to the ordering of the rare-earth moments, the expected ferroelectric transition is frustrated by the large size difference between Cu and Ti at the B-site. Interestingly, this leads these compounds to attain a rare and unique combination of desirable paraelectric properties with high dielectric constants, low losses and weak temperature and frequency dependencies. First-principles calculations establish these exceptional properties result from a combination of two effects. A significant difference in the MO5 polyhedral sizes for M = Cu and M = Ti suppress the expected co-operative tilt pattern of these polyhedra, required for the ferroelectric transition, leading to relatively large values of the dielectric constant for every compound investigated in this series. Additionally, it is shown that the majority contribution to the dielectric constant arises from intermediate-frequency polar vibrational modes, making it relatively stable against any temperature variation. Changes in the temperature stability of the dielectric constant amongst different members of this series are shown to arise from changes in relative contributions from soft polar modes.

preprint2010arXiv

Structural phases of strained LaAlO3 driven by octahedral tilt instabilities

We investigate the effect of epitaxial strain on [001]-oriented LaAlO3 using first-principles density functional calculations. We find a series of structural phase transitions between states characterized by distinct patterns of tilting of the AlO6 octahedra. By tuning the biaxial strain from compressive to tensile, we induce an evolution in the crystal structure in which the tilt axis changes from out-of-plane to in-plane, corresponding to space groups I4/mcm and Imma. We also study the effect of uniaxial relaxation of the usual biaxial epitaxial constraint and explore this as a mechanism for selectively stabilizing different patterns of octahedral tilts.

preprint2009arXiv

Electron-lattice instabilities suppress cuprate-like electronic structures in SrFeO$_3$/SrTiO$_3$ superlattices

Using {\it ab initio} density functional theory we explore the behavior of thin layers of metallic $d^4$ SrFeO$_3$ confined between the $d^0$ dielectric SrTiO$_3$ in a superlattice geometry. We find the presence of insulating SrTiO$_3$ spacer layers strongly affects the electronic properties of SrFeO$_3$: For single SrFeO$_3$ layers constrained to their bulk cubic structure, the Fermi surface is two-dimensional, nested and resembles the hole-doped superconducting cuprates. A Jahn-Teller instability couples to an octahedral tilt mode, however, to remove this degeneracy resulting in insulating superlattices.

preprint2009arXiv

First-principles modeling of ferroelectric capacitors via constrained-D calculations

First-principles modeling of ferroelectric capacitors presents several technical challenges, due to the coexistence of metallic electrodes, long-range electrostatic forces and short-range interface chemistry. Here we show how these aspects can be efficiently and accurately rationalized by using a finite-field density-functional theory formalism in which the fundamental electrical variable is the displacement field D. By performing calculations on model Pt/BaTiO3/Pt and Au/BaZrO3/Au capacitors we demonstrate how the interface-specific and bulk-specific properties can be identified and rigorously separated. Then, we show how the electrical properties of capacitors of arbitrary thickness and geometry (symmetric or asymmetric) can be readily reconstructed by using such information. Finally, we show how useful observables such as polarization and dielectric, piezoelectric and electrostrictive coefficients are easily evaluated as a byproduct of the above procedure. We apply this methodology to elucidate the relationship between chemical bonding, Schottky barriers and ferroelectric polarization at simple-metal/oxide interfaces. We find that BO2-electrode interfaces behave analogously to a layer of linear dielectric put in series with a bulk-like perovskite film, while a significant non-linear effect occurs at AO-electrode interfaces.

preprint2008arXiv

Electric displacement as the fundamental variable in electronic-structure calculations

Finite-field calculations in periodic insulators are technically and conceptually challenging, due to fundamental problems in defining polarization in extended solids. While significant progress has been made recently with the establishment of techniques to fix the electric field E or the macroscopic polarization P in first-principles calculations, both methods lack the ease of use and conceptual clarity of standard zero-field calculations. Here we develop a new formalism in which the electric displacement D, rather than E or P, is the fundamental electrical variable. Fixing D has the intuitive interpretation of imposing open-circuit electrical boundary conditions, which is particularly useful in studying ferroelectric systems. Furthermore, the analogy to open-circuit capacitors suggests an appealing reformulation in terms of free charges and potentials, which dramatically simplifies the treatment of stresses and strains. Using PbTiO3 as an example, we show that our technique allows full control over the electrical variables within the density functional formalism.

preprint2008arXiv

Enhancement of ferroelectricity at metal/oxide interfaces

By performing first-principles calculations on four capacitor structures based on BaTiO3 and PbTiO3, we determine the intrinsic interfacial effects that are responsible for the destabilization of the polar state in thin-film ferroelectric devices. We show that, contrary to the established interpretation, the interfacial capacitance depends crucially on the local chemical environment of the interface through the force constants of the metal-oxide bonds, and is not necessarily related to the bulk screening properties of the electrode material. In particular, in the case of interfaces between AO-terminated perovskites and simple metals, we demonstrate a novel mechanism of "interfacial ferroelectricity" that produces an overall enhancement of the ferroelectric instability of the film, rather than a suppression as is usually assumed. The resulting "negative dead layer" suggests a route to thin-film ferroelectric devices that are free of deleterious size effects.

preprint2007arXiv

Self-interaction correction with Wannier functions

We describe the behavior of the Perdew-Zunger self-interaction-corrected local density approximation (SIC-LDA) functional when implemented in a plane-wave pseudopotential formalism with Wannier functions. Prototypical semiconductors and wide-bandgap oxides show a large overcorrection of the LDA bandgap. Application to transition-metal oxides and elements with d electrons is hindered by a serious breaking of the spherical symmetry, which appears even in a closed shell free atom. Our results indicate that, when all spherical approximations are lifted, the general applicability of orbital-dependent potentials is very limited and should be reconsidered in favor of rotationally invariant functionals.

preprint2007arXiv

Tri-layer superlattices: A route to magnetoelectric multiferroics?

We explore computationally the formation of tri-layer superlattices as an alternative approach for combining ferroelectricity with magnetism to form magnetoelectric multiferroics. We find that the contribution to the superlattice polarization from tri-layering is small compared to typical polarizations in conventionalferroelectrics, and the switchable ferroelectric component is negligible. In contrast, we show that epitaxial strain and ``negative pressure'' can yield large, switchable polarizations that are compatible with the coexistence of magnetism, even in materials with no active ferroelectric ions.

preprint2005arXiv

Accurate polarization within a unified Wannier function formalism

We present an alternative formalism for calculating the maximally localized Wannier functions in crystalline solids, obtaining an expression which is extremely simple and general. In particular, our scheme is exactly invariant under Brillouin zone folding, and therefore it extends trivially to the Gamma-point case. We study the convergence properties of the Wannier functions, their quadratic spread and centers as obtained by our simplified technique. We show how this convergence can be drastically improved by a simple and inexpensive ``refinement'' step, which allows for very efficient and accurate calculations of the polarization in zero external field.

preprint2005arXiv

Recent progress in first-principles studies of magnetoelectric multiferroics

Materials that combine magnetic and ferroelectric properties have generated increasing interest over the last few years, due to both their diverse properties and their potential utility in new types of magnetoelectric device applications. In this review we discuss recent progress in the study of such magnetoelectric multiferroics which has been achieved using computational first-principles methods based on density functional theory. In particular, we show how first-principles methods have been successfully used to explain various properties of multiferroic materials and to predict novel effects and new systems that exhibit multiferroic properties.

preprint2003arXiv

A search for ferromagnetism in transition-metal-doped piezoelectric ZnO

We present the results of a computational study of ZnO in the presence of Co and Mn substitutional impurities. The goal of our work is to identify potential ferromagnetic ground states within the (Zn,Co)O or (Zn,Mn)O material systems that are also good candidates for piezoelectricity. We find that, in contrast to previous results, robust ferromagnetism is not obtained by substitution of Co or Mn on the Zn site, unless additional carriers (holes) are also incorporated. We propose a practical scheme for achieving such $p$-type doping in ZnO.