Orbital magnetization in periodic insulators

Orbital Hall Effect Accompanying Quantum Hall Effect: Landau Levels Cause Orbital Polarized Edge Currents

The quantum Hall effect emerges when two-dimensional samples are subjected to strong magnetic fields at low temperatures: Topologically protected edge states cause a quantized Hall conductivity in multiples of e^{2}/h. Here we show that the quantum Hall effect is accompanied by an orbital Hall effect. Our quantum mechanical calculations fit well the semiclassical interpretation in terms of "skipping orbits." The chiral edge states of a quantum Hall system are orbital polarized akin to a hypothetical orbital version of the quantum anomalous Hall effect in magnetic systems. The orbital Hall resistivity scales quadratically with the magnetic field, making it the dominant effect at high fields.

Theory of spin and orbital Edelstein effects

In systems with broken spatial inversion symmetry, such as surfaces, interfaces, or bulk systems lacking an inversion center, the application of a charge current can generate finite spin and orbital densities associated with a nonequilibrium magnetization, which is known as spin and orbital Edelstein effect (SEE and OEE), respectively. Early reports on this current-induced magnetization focus on two-dimensional Rashba systems, in which an in-plane nonequilibrium spin density is generated perpendicular to the applied charge current. However, until today, a large variety of materials have been theoretically predicted and experimentally demonstrated to exhibit a sizeable Edelstein effect, which comprises contributions from the spin as well as the orbital degrees of freedom, and whose associated magnetization may be out of plane, nonorthogonal, and even parallel to the applied charge current, depending on the system's particular symmetries. In this review, we give an overview on the most commonly used theoretical approaches for the discussion and prediction of the SEE and OEE. Further, we introduce a selection of the most intensely discussed materials exhibiting a finite Edelstein effect, and give a brief summary of common experimental techniques.

Magnon Orbital Nernst Effect in Honeycomb Antiferromagnets without Spin-Orbit Coupling

Recently, topological responses of magnons have emerged as a central theme in magnetism and spintronics. However, resulting Hall responses are typically weak and infrequent, since, according to present understanding, they arise from effective spin-orbit couplings, which are weaker compared to the exchange energy. Here, by investigating transport properties of magnon orbital moments, we predict that the magnon orbital Nernst effect is an intrinsic characteristic of the honeycomb antiferromagnet and therefore, it manifests even in the absence of spin-orbit coupling. For the electric detection, we propose an experimental scheme based on the magnetoelectric effect. Our results break the conventional wisdom that the Hall transport of magnons requires spin-orbit coupling by predicting the magnon orbital Nernst effect in a system without it, which leads us to envision that our work initiates the intensive search for various magnon Hall effects in generic magnetic systems with no reliance on spin-orbit coupling.

High spin axion insulator

Axion insulators possess a quantized axion field θ = π protected by combined lattice and time-reversal symmetry, holding great potential for device applications in layertronics and quantum computing. Here, we propose a high-spin axion insulator (HSAI) defined in large spin-s representation, which maintains the same inherent symmetry but possesses a notable axion field θ = (s + 1/2)2π. Such distinct axion field is confirmed independently by the direct calculation of the axion term using hybrid Wannier functions, layer-resolved Chern numbers, as well as the topological magneto-electric effect. We show that the guaranteed gapless quasi-particle excitation is absent at the boundary of the HSAI despite its integer surface Chern number, hinting an unusual quantum anomaly violating the conventional bulk-boundary correspondence. Furthermore, we ascertain that the axion field θ can be precisely tuned through an external magnetic field, enabling the manipulation of bonded transport properties. The HSAI proposed here can be experimentally verified in ultra-cold atoms by the quantized non-reciprocal conductance or topological magnetoelectric response. Our work enriches the understanding of axion insulators in condensed matter physics, paving the way for future device applications.

Chiral Spectroscopy of Bulk Systems with Propagated Localized Orbitals

We present approaches for the simulation of electronic circular dichroism, Raman, and Raman optical activity (ROA) spectra for isolated and periodic systems as well as subsystem analysis thereof. The method is based on the use of time-dependent maximally localized Wannier functions in the CP2K package and accounts for origin dependencies inherent to the Gaussian and plane wave with pseudopotentials approach as well as the origin dependence of the magnetic dipole and electric quadrupole operators. Tests on the H-bonded enantiomers of alanine by harmonic normal-mode analysis and on an aqueous solution of l-alanine by ab initio molecular dynamics obeying periodic boundary conditions (PBCs) are presented as total and subsystem-resolved spectra. To our knowledge, this is the first instance of an ROA spectrum derived from real-time propagation obeying PBCs and the first ROA simulation considering off-, pre-, and on-resonance effects within PBCs.

Topological magnetoelectric response in ferromagnetic axion insulators

The topological magnetoelectric effect (TME) is a hallmark response of the topological field theory, which provides a paradigm shift in the study of emergent topological phenomena. However, its direct observation is yet to be realized due to the demanding magnetic configuration required to gap all surface states. Here, we theoretically propose that axion insulators with a simple ferromagnetic configuration, such as the MnBi2Te4/(Bi2Te3)n family, provide an ideal playground to realize the TME. In the designed triangular prism geometry, all the surface states are magnetically gapped. Under a vertical electric field, the surface Hall currents give rise to a nearly half-quantized orbital moment, accompanied by a gapless chiral hinge mode circulating in parallel. Thus, the orbital magnetization from the two topological origins can be easily distinguished by reversing the electric field. Our work paves the way for direct observation of the TME in realistic axion-insulator materials.

Monopole-like orbital-momentum locking and the induced orbital transport in topological chiral semimetals

The interplay between chirality and topology nurtures many exotic electronic properties. For instance, topological chiral semimetals display multifold chiral fermions that manifest nontrivial topological charge and spin texture. They are an ideal playground for exploring chirality-driven exotic physical phenomena. In this work, we reveal a monopole-like orbital-momentum locking texture on the three-dimensional Fermi surfaces of topological chiral semimetals with B20 structures (e.g., RhSi and PdGa). This orbital texture enables a large orbital Hall effect (OHE) and a giant orbital magnetoelectric (OME) effect in the presence of current flow. Different enantiomers exhibit the same OHE which can be converted to the spin Hall effect by spin-orbit coupling in materials. In contrast, the OME effect is chirality-dependent and much larger than its spin counterpart. Our work reveals the crucial role of orbital texture for understanding OHE and OME effects in topological chiral semimetals and paves the path for applications in orbitronics, spintronics, and enantiomer recognition.

A quantum ruler for orbital magnetism in moiré quantum matter

For almost a century, magnetic oscillations have been a powerful "quantum ruler" for measuring Fermi surface topology. In this study, we used Landau-level spectroscopy to unravel the energy-resolved valley-contrasting orbital magnetism and large orbital magnetic susceptibility that contribute to the energies of Landau levels of twisted double-bilayer graphene. These orbital magnetism effects led to substantial deviations from the standard Onsager relation, which manifested as a breakdown in scaling of Landau-level orbits. These substantial magnetic responses emerged from the nontrivial quantum geometry of the electronic structure and the large length scale of the moiré lattice potential. Going beyond traditional measurements, Landau-level spectroscopy performed with a scanning tunneling microscope offers a complete quantum ruler that resolves the full energy dependence of orbital magnetic properties in moiré quantum matter.

Transport of Spin Magnetic Multipole Moments Carried by Bloch Quasiparticles

Magnetic ordering beyond the standard dipolar order has attracted significant attention in recent years, but it remains an open question how to effectively manipulate such nontrivial order parameters using external perturbations such as electric currents or fields. In particular, it is desirable to have a conceptual tool similar to nonequilibrium spin currents in spintronics to describe the creation and transport of multipole moments. In this context, we present a theory for Cartesian spin magnetic multipole moments of Bloch quasiparticles and their transport based on a general gauge-invariant formula obtained using the wave packet approach. As a concrete example, we point out that the low-energy Hamiltonian of phosphorene subject to a perpendicular electric field has a valley structure that hosts magnetic octupole moments. The magnetic octupole moments can be exhibited by an in-plane electric current and lead to accumulation of staggered spin densities at the corners of a rectangular sample. Our Letter paves the way for systematically seeking and utilizing quasiparticles with higher-order magnetic multipole moments in crystal materials towards the emergence of multipoletronics.

Time-Reversal-Even Nonlinear Current Induced Spin Polarization

We propose a time-reversal-even spin generation in second order of electric fields, which dominates the current induced spin polarization in a wide class of centrosymmetric nonmagnetic materials, and leads to a novel nonlinear spin-orbit torque in magnets. We reveal a quantum origin of this effect from the momentum space dipole of the anomalous spin polarizability. First-principles calculations predict sizable spin generations in several nonmagnetic hcp metals, in monolayer TiTe_{2}, and in ferromagnetic monolayer MnSe_{2}, which can be detected in experiment. Our work opens up the broad vista of nonlinear spintronics in both nonmagnetic and magnetic systems.

Theoretical methods for structural phase transitions in elemental solids at extreme conditions: statics and dynamics

In recent years, theoretical studies have moved from a traditionally supporting role to a more proactive role in the research of phase transitions at high pressures. In many cases, theoretical prediction leads the experimental exploration. This is largely owing to the rapid progress of computer power and theoretical methods, particularly the structure prediction methods tailored for high-pressure applications. This review introduces commonly used structure searching techniques based on static and dynamic approaches, their applicability in studying phase transitions at high pressure, and new developments made toward predicting complex crystalline phases. Successful landmark studies for each method are discussed, with an emphasis on elemental solids and their behaviors under high pressure. The review concludes with a perspective on outstanding challenges and opportunities in the field.

Colossal Orbital Edelstein Effect in Noncentrosymmetric Superconductors

In superconductors that lack inversion symmetry, the flow of supercurrent can induce a nonvanishing magnetization, a phenomenon which is at the heart of nondissipative magnetoelectric effects, also known as Edelstein effects. For electrons carrying spin and orbital moments, a question of fundamental relevance deals with the orbital nature of magnetoelectric effects in conventional spin-singlet superconductors with Rashba coupling. Remarkably, we find that the supercurrent-induced orbital magnetization is more than 1 order of magnitude greater than that due to the spin, giving rise to a colossal magnetoelectric effect. The induced orbital magnetization is shown to be sign tunable, with the sign change occurring for the Fermi level lying in proximity of avoiding crossing points in the Brillouin zone. In the presence of superconducting phase inhomogeneities, a modulation of the Edelstein signal on the scale of the superconducting coherence length appears, leading to domains with opposite orbital moment orientations. These hallmarks are robust to real-space self-consistent treatment of the superconducting order parameter. The orbital-dominated magnetoelectric phenomena, hence, have clear-cut marks for detection both in the bulk and at the edge of the system and are expected to be a general feature of multiorbital superconductors with inversion symmetry breaking.

Intelligent infrared sensing enabled by tunable moiré quantum geometry

Quantum geometric properties of Bloch wave functions in solids, that is, Berry curvature and the quantum metric, are known to significantly influence the ground- and excited-state behaviour of electrons^1-5. The bulk photovoltaic effect (BPVE), a nonlinear phenomenon depending on the polarization of excitation light, is largely governed by the quantum geometric properties in optical transitions^6-10. Infrared BPVE has yet to be observed in graphene or moiré systems, although exciting strongly correlated phenomena related to quantum geometry have been reported in this emergent platform^11-14. Here we report the observation of tunable mid-infrared BPVE at 5 µm and 7.7 µm in twisted double bilayer graphene (TDBG), arising from the moiré-induced strong symmetry breaking and quantum geometric contribution. The photoresponse depends substantially on the polarization state of the excitation light and is highly tunable by external electric fields. This wide tunability in quantum geometric properties enables us to use a convolutional neural network^15,16 to achieve full-Stokes polarimetry together with wavelength detection simultaneously, using only one single TDBG device with a subwavelength footprint of merely 3 × 3 µm². Our work not only reveals the unique role of moiré engineered quantum geometry in tunable nonlinear light-matter interactions but also identifies a pathway for future intelligent sensing technologies in an extremely compact, on-chip manner.

Phonon Magnetic Moment from Electronic Topological Magnetization

The traditional theory of magnetic moments for chiral phonons is based on the picture of the circular motion of the Born effective charge, typically yielding a small fractional value of the nuclear magneton. Here we investigate the adiabatic evolution of electronic states induced by the lattice vibration of a chiral phonon and obtain an electronic orbital magnetization in the form of a topological second Chern form. We find that the traditional theory needs to be refined by introducing a k resolved Born effective charge, and identify another contribution from the phonon-modified electronic energy together with the momentum-space Berry curvature. The second Chern form can diverge when there is a Yang's monopole near the parameter space of interest as illustrated by considering a phonon at the Brillouin zone corner in a gapped graphene model. We also find large magnetic moments for the optical phonon in bulk topological materials where nontopological contribution is also important. Our results agree with recent observations in experiments.

Size and shape effects in the orbital magnetization of TMDs monolayers

The intrinsic orbital magnetization of a TMD monolayer is usually calculated for periodic crystals without mentioning the geometrical shape of samples and boundary conditions (BCs) for the electron wave functions. Such approaches, based on Bloch's theorem, involves a contribution from the Berry curvature, also in the case when the system is described by the two-band minimal model (Xiaoet al2012Phys. Rev. Lett.108196802). In the present paper, we show that the geometrical and topological properties of the specimen, as well as the BCs, play an important role in the problem of valley orbital magnetization even for a macroscopic sample.

Absence of Hall effect due to Berry curvature in phase space

Transverse current due to Berry curvature in phase space is formulated based on the Boltzmann equations with the semiclassical equations of motion for an electron wave packet. It is shown that the Hall effect due to the phase space Berry curvature is absent because the contributions from “anomalous velocity” and “effective Lorentz force” are completely cancelled out.

Electronic-structure methods for materials design

The accuracy and efficiency of electronic-structure methods to understand, predict and design the properties of materials has driven a new paradigm in research. Simulations can greatly accelerate the identification, characterization and optimization of materials, with this acceleration driven by continuous progress in theory, algorithms and hardware, and by adaptation of concepts and tools from computer science. Nevertheless, the capability to identify and characterize materials relies on the predictive accuracy of the underlying physical descriptions, and on the ability to capture the complexity of realistic systems. We provide here an overview of electronic-structure methods, of their application to the prediction of materials properties, and of the different strategies employed towards the broader goals of materials design and discovery.

Electron g-factor in nanostructures: continuum media and atomistic approach

We report studies of [Formula: see text]-dependent Landé g-factor, performed by both continuous media approximation [Formula: see text] method, and atomistic tight-binding [Formula: see text] approach. We propose an effective, mesoscopic model for InAs that we are able to successfully compare with atomistic calculations, for both very small and very large nanostructures, with a number of atoms reaching over 60 million. Finally, for nanostructure dimensions corresponding to near-zero g-factor we report electron spin states anti-crossing as a function of system size, despite no shape-anisotropy nor strain effects included, and merely due to breaking of atomistic symmetry of cation/anion planes constituting the system.

Voltage-Controlled Magnetic Reversal in Orbital Chern Insulators

Chern insulator ferromagnets are characterized by a quantized anomalous Hall effect and have so far been identified experimentally in magnetically doped topological insulator thin films and in bilayer graphene moiré superlattices. We classify Chern insulator ferromagnets as either spin or orbital, depending on whether the orbital magnetization results from spontaneous spin polarization combined with spin-orbit interactions, as in the magnetically doped topological insulator case, or directly from spontaneous orbital currents, as in the moiré superlattice case. We argue that, in a given magnetic state, characterized, for example, by the sign of the anomalous Hall effect, the magnetization of an orbital Chern insulator will often have opposite signs for weak n and weak p electrostatic or chemical doping. This property enables pure electrical switching of a magnetic state in the presence of a fixed magnetic field.

Ab Initio Studies of Exciton g Factors: Monolayer Transition Metal Dichalcogenides in Magnetic Fields

The effect of a magnetic field on the optical absorption in semiconductors has been measured experimentally and modeled theoretically for various systems in previous decades. We present a new first-principles approach to systematically determine the response of excitons to magnetic fields, i.e., exciton g factors. By utilizing the GW-Bethe-Salpeter equation methodology we show that g factors extracted from the Zeeman shift of electronic bands are strongly renormalized by many-body effects which we trace back to the extent of the excitons in reciprocal space. We apply our approach to monolayers of transition metal dichalcogenides (MoS_{2}, MoSe_{2}, MoTe_{2}, WS_{2}, and WSe_{2}) with strongly bound excitons for which g factors are weakened by about 30%.

Numerical study of disorder on the orbital magnetization in two dimensions

The modern theory of orbital magnetization (OM) was developed by using Wannier function method, which has a formalism similar with the Berry phase. In this manuscript, we perform a numerical study on the fate of the OM under disorder, by using this method on the Haldane model in two dimensions, which can be tuned between a normal insulator or a Chern insulator at half filling. The effects of increasing disorder on OM for both cases are simulated. Energy renormalization shifts are observed in the weak disorder regime and topologically trivial case, which was predicted by a self-consistentT-matrix approximation. Besides this, two other phenomena can be seen. One is the localization trend of the band orbital magnetization. The other is the remarkable contribution from topological chiral states arising from nonzero Chern number or large value of integrated Berry curvature. If the fermi energy is fixed at the gap center of the clean system, there is an enhancement of |M| at the intermediate disorder, for both cases of normal and Chern insulators, which can be attributed to the disorder induced topological metal state before localization.

Universal Behavior of the Thermal Hall Conductivity

We report theoretical and experimental analyses of the thermal Hall conductivity in correlated systems. For both fermionic and bosonic excitations with nontrivial topology, we show that at "intermediate" temperatures, the thermal Hall conductivity exhibits an unexpected universal scaling with a simple exponential form. At low temperatures, it behaves differently and reflects the spectral properties of underlying excitations. Our predictions are examined as examples in two prototype compounds, the quantum paraelectric SrTiO_{3} and the spin-liquid compound RuCl_{3}. The experimental data can be largely covered by our proposed minimal phenomenological model independent of microscopic details, revealing dominant bosonic contributions in SrTiO_{3} and gapped fermionic excitations in RuCl_{3}. Our work establishes a phenomenological link between microscopic models and experimental data and provides a unified basis for analyzing the thermal Hall conductivity in correlated systems over a wide temperature region.

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

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

Berry Phase Effects in Dipole Density and the Mott Relation

We provide a unified semiclassical theory for thermoelectric responses of any observable represented by an operator θ[over ^] that is well defined in periodic crystals. The Einstein and Mott relations are established generally in the presence of Berry phase effects for various physical realizations of θ[over ^] in electronic systems, including the familiar case of the electric current as well as the currently controversial cases of the spin polarization and spin current. The magnetization current, which has been proven indispensable in the thermoelectric response of electric current, is generalized to the cases of various θ[over ^]. In our theory the dipole density of a physical quantity emerges and plays a vital role, which contains not only the statistical sum of the dipole moment of θ[over ^] but also a Berry phase correction.

Theory of Current-Induced Angular Momentum Transfer Dynamics in Spin-Orbit Coupled Systems

Motivated by the importance of understanding various competing mechanisms to the current-induced spin-orbit torque on magnetization in complex magnets, we develop a theory of current-induced spin-orbital coupled dynamics in magnetic heterostructures. The theory describes angular momentum transfer between different degrees of freedom in solids, e.g., the electron orbital and spin, the crystal lattice, and the magnetic order parameter. Based on the continuity equations for the spin and orbital angular momenta, we derive equations of motion that relate spin and orbital current fluxes and torques describing the transfer of angular momentum between different degrees of freedom, achieved in a steady state under an applied external electric field. We then propose a classification scheme for the mechanisms of the current-induced torque in magnetic bilayers. We evaluate the sources of torque using density functional theory, effectively capturing the impact of the electronic structure on these quantities. We apply our formalism to two different magnetic bilayers, Fe/W(110) and Ni/W(110), which are chosen such that the orbital and spin Hall effects in W have opposite sign and the resulting spin- and orbital-mediated torques can compete with each other. We find that while the spin torque arising from the spin Hall effect of W is the dominant mechanism of the current-induced torque in Fe/W(110), the dominant mechanism in Ni/W(110) is the orbital torque originating in the orbital Hall effect of the non-magnetic substrate. Thus the effective spin Hall angles for the total torque are negative and positive in the two systems. Our prediction can be experimentally identified in moderately clean samples, where intrinsic contributions dominate. This clearly demonstrates that our formalism is ideal for studying the angular momentum transfer dynamics in spin-orbit coupled systems as it goes beyond the "spin current picture" by naturally incorporating the spin and orbital degrees of freedom on an equal footing. Our calculations reveal that, in addition to the spin and orbital torque, other contributions such as the interfacial torque and self-induced anomalous torque within the ferromagnet are not negligible in both material systems.

Mixed topological semimetals driven by orbital complexity in two-dimensional ferromagnets

The concepts of Weyl fermions and topological semimetals emerging in three-dimensional momentum space are extensively explored owing to the vast variety of exotic properties that they give rise to. On the other hand, very little is known about semimetallic states emerging in two-dimensional magnetic materials, which present the foundation for both present and future information technology. Here, we demonstrate that including the magnetization direction into the topological analysis allows for a natural classification of topological semimetallic states that manifest in two-dimensional ferromagnets as a result of the interplay between spin-orbit and exchange interactions. We explore the emergence and stability of such mixed topological semimetals in realistic materials, and point out the perspectives of mixed topological states for current-induced orbital magnetism and current-induced domain wall motion. Our findings pave the way to understanding, engineering and utilizing topological semimetallic states in two-dimensional spin-orbit ferromagnets.

Quantized angular momentum in topological optical systems

The Chern index characterizes the topological phases of nonreciprocal photonic systems. Unlike in electronics, the photonic Chern number has no clear physical meaning, except that it determines the number of unidirectional edge states supported by an interface with a trivial mirror. Here, we fill in this gap by demonstrating that the photonic Chern number can be understood as the quantum of the light-angular momentum in a photonic insulator cavity. It is proven that for a large cavity, the thermal fluctuation-induced angular momentum is precisely quantized in the band-gaps of the bulk states. The nontrivial expectation of the light angular momentum is due to a circulation of thermal energy in closed orbits. Remarkably, this result can be extended to systems without a topological classification, and in such a case the “quantum” of the angular momentum density is determined by the net number of unidirectional edge states supported by the cavity walls.

While the Chern number is used to characterize topological states in both electronic and photonic systems, its direct physical meaning in photonics has not yet been established. Here, Silveirinha shows that the photonic Chern number can be understood as a quantum of the optical angular momentum.

Ultrathin Films of Superconducting Metals as a Platform for Topological Superconductivity

The ingredients normally required to achieve topological superconductivity (TSC) are Cooper pairing, broken inversion symmetry, and broken time-reversal symmetry. We present a theoretical exploration of the possibility of using ultrathin films of superconducting metals as a platform for TSC. Because they necessarily break inversion symmetry when prepared on a substrate and have intrinsic Cooper pairing, they can be TSCs when time-reversal symmetry is broken by an external magnetic field. Using microscopic density functional theory calculations we show that, for ultrathin Pb and β-Sn superconductors, the position of the Fermi level can be tuned to quasi-2D band extrema energies using strain, and that the g factors of states at time-reversal invariant momenta can be extremely large, enhancing the influence of external magnetic fields.

Drude weight and superconducting weight

The Drude weight D (also called charge stiffness) measures the effective electron density contributing to dc conductivity; it vanishes in insulators. It is a general concept, which applies to any metal, including cases with disorder and electron-electron interaction. We provide a thorough analysis of D and of its meaning, both in the general case and in the special case of band metals, where we also show that D has a close relationship to orbital magnetization. The superconducting weight D _s measures instead the superfluid density accounting for the Meissner effect. The two quantities D and D _s are the main criteria to discriminate between insulators, metals, and superconductors.

Electronic orbital response of regular extended and infinite periodic systems to magnetic fields. I. Theoretical foundations for static case

A theoretical treatment for the orbital response of an infinite, periodic system to a static, homogeneous, magnetic field is presented. It is assumed that the system of interest has an energy gap separating occupied and unoccupied orbitals and a zero Chern number. In contrast to earlier studies, we do not utilize a perturbation expansion, although we do assume the field is sufficiently weak that the occurrence of Landau levels can be ignored. The theory is developed by analyzing results for large, finite systems and also by comparing with the analogous treatment of an electrostatic field. The resulting many-electron Hamilton operator is forced to be hermitian, but hermiticity is not preserved, in general, for the subsequently derived single-particle operators that determine the electronic orbitals. However, we demonstrate that when focusing on the canonical solutions to the single-particle equations, hermiticity is preserved. The issue of gauge-origin dependence of approximate solutions is addressed. Our approach is compared with several previously proposed treatments, whereby limitations in some of the latter are identified.

Orbital Edelstein Effect as a Condensed-Matter Analog of Solenoids

We theoretically study current-induced orbital magnetization in a chiral crystal. This phenomenon is an orbital version of the Edelstein effect. We propose an analogy between the current-induced orbital magnetization and an Ampère field in a solenoid in classical electrodynamics. To quantify this effect, we define a dimensionless parameter from the response coefficients relating a current density with an orbital magnetization. This dimensionless parameter can be regarded as a number of turns within a unit cell when the crystal is regarded as a solenoid, and it represents how "chiral" the crystal is. By focusing on the dimensionless parameter, one can design a band structure that realizes the induction of large orbital magnetization. In particular, a Weyl semimetal with all of the Weyl nodes close to the Fermi energy can have a large value for this dimensionless parameter, which can exceed that of a classical solenoid.

Advanced capabilities for materials modelling with Quantum ESPRESSO

Quantum EXPRESSO is an integrated suite of open-source computer codes for quantum simulations of materials using state-of-the-art electronic-structure techniques, based on density-functional theory, density-functional perturbation theory, and many-body perturbation theory, within the plane-wave pseudopotential and projector-augmented-wave approaches. Quantum EXPRESSO owes its popularity to the wide variety of properties and processes it allows to simulate, to its performance on an increasingly broad array of hardware architectures, and to a community of researchers that rely on its capabilities as a core open-source development platform to implement their ideas. In this paper we describe recent extensions and improvements, covering new methodologies and property calculators, improved parallelization, code modularization, and extended interoperability both within the distribution and with external software.

Toward surface orbitronics: giant orbital magnetism from the orbital Rashba effect at the surface of sp-metals

As the inversion symmetry is broken at a surface, spin-orbit interaction gives rise to spin-dependent energy shifts - a phenomenon which is known as the spin Rashba effect. Recently, it has been recognized that an orbital counterpart of the spin Rashba effect - the orbital Rashba effect - can be realized at surfaces even without spin-orbit coupling. Here, we propose a mechanism for the orbital Rashba effect based on sp orbital hybridization, which ultimately leads to the electric polarization of surface states. For the experimentally well-studied system of a BiAg₂ monolayer, as a proof of principle, we show from first principles that this effect leads to chiral orbital textures in k-space. In predicting the magnitude of the orbital moment arising from the orbital Rashba effect, we demonstrate the crucial role played by the Berry phase theory for the magnitude and variation of the orbital textures. As a result, we predict a pronounced manifestation of various orbital effects at surfaces, and proclaim the orbital Rashba effect to be a key platform for surface orbitronics.

Prototypical topological orbital ferromagnet γ-FeMn

We predict from first principles an entirely topological orbital magnetization in the noncoplanar bulk antiferromagnet γ-FeMn originating in the nontrivial topology of the underlying spin structure, without any reference to spin-orbit interaction. Studying the influence of strain, composition ratio, and spin texture on the topological orbital magnetization and the accompanying topological Hall effect, we promote the scalar spin chirality as key mechanism lifting the orbital degeneracy. The system is thus a prototypical topological orbital ferromagnet, the macroscopic orbital magnetization of which is prominent even without spin-orbit coupling. One of the remarkable features of γ-FeMn is the possibility for pronounced orbital magnetostriction mediated by the complex spin topology in real space.

The inverse thermal spin-orbit torque and the relation of the Dzyaloshinskii-Moriya interaction to ground-state energy currents

Using the Kubo linear-response formalism we derive expressions to calculate the electronic contribution to the heat current generated by magnetization dynamics in ferromagnetic metals with broken inversion symmetry and spin-orbit interaction (SOI). The effect of producing heat currents by magnetization dynamics constitutes the Onsager reciprocal of the thermal spin-orbit torque (TSOT), i.e. the generation of torques on the magnetization due to temperature gradients. We find that the energy current driven by magnetization dynamics contains a contribution from the Dzyaloshinskii-Moriya interaction (DMI), which needs to be subtracted from the Kubo linear response of the energy current in order to extract the heat current. We show that the expressions of the DMI coefficient can be derived elegantly from the DMI energy current. Guided by formal analogies between the Berry phase theory of DMI on the one hand and the modern theory of orbital magnetization on the other hand we are led to an interpretation of the latter in terms of energy currents as well. Based on ab initio calculations we investigate the electronic contribution to the heat current driven by magnetization dynamics in Mn/W(0 0 1) magnetic bilayers. We predict that fast domain walls drive strong heat currents.

Irrelevance of the Boundary on the Magnetization of Metals

The macroscopic current density responsible for the mean magnetization M of a uniformly magnetized bounded sample is localized near its surface. In order to evaluate M one needs the current distribution in the whole sample: bulk and boundary. In recent years it has been shown that the boundary has no effect on M in insulators: therein, M admits an alternative expression not based on currents. M can be expressed in terms of the bulk electron distribution only, which is "nearsighted" (exponentially localized); this virtue is not shared by metals, having a qualitatively different electron distribution. We show, by means of simulations on paradigmatic model systems, that even in metals the M value can be retrieved in terms of the bulk electron distribution only.

Current response of a topological insulator to a static Zeeman field

We study the magnetoelectric coupling at the surface of a topological insulator. We are particularly interested in the surface current induced by a static Zeeman/exchange field. This surface current can be related to the orbital magnetization of the system. For an insulator with zero Chern number, the orbital magnetization is independent of the details at the boundary. With the appearance of surface states in the topological insulator, it is not immediately obvious if the response is affected by the conditions at the surface. We investigate this question using exact diagonalization to a lattice model. By applying a time-reversal symmetry-breaking term near the boundary, even if the surface states are gapped out, we still find no change in the surface current. This arises from cancelations between Pauli and Van Vleck contributions between surface and bulk scattering states. We also show that the surface current response is independent of the chemical potential when it is within the bulk gap. Our results are consistent with the claim that orbital magnetization is a bulk property.

Quantum Quenches in Chern Insulators

We explore the nonequilibrium response of Chern insulators. Focusing on the Haldane model, we study the dynamics induced by quantum quenches between topological and nontopological phases. A notable feature is that the Chern number, calculated for an infinite system, is unchanged under the dynamics following such a quench. However, in finite geometries, the initial and final Hamiltonians are distinguished by the presence or absence of edge modes. We study the edge excitations and describe their impact on the experimentally observable edge currents and magnetization. We show that, following a quantum quench, the edge currents relax towards new equilibrium values, and that there is light-cone spreading of the currents into the interior of the sample.

The Calculation of NMR Chemical Shifts in Periodic Systems Based on Gauge Including Atomic Orbitals and Density Functional Theory

We present here a method that can calculate NMR shielding tensors from first principles for systems with translational invariance. Our approach is based on Kohn-Sham density functional theory and gauge-including atomic orbitals. Our scheme determines the shielding tensor as the second derivative of the total electronic energy with respect to an external magnetic field and a nuclear magnetic moment. The induced current density due to a periodic perturbation from nuclear magnetic moments is obtained through numerical differentiation, whereas the influence of the responding perturbation in terms of the external magnetic field is evaluated analytically. The method is implemented into the periodic program BAND. It employs a Bloch basis set made up of Slater-type or numeric atomic orbitals and represents the Kohn-Sham potential fully without the use of effective core potentials. Results from calculations of NMR shielding constants based on the present approach are presented for isolated molecules as well as systems with one-, two- and three-dimensional periodicity. The reported values are compared to experiment and results from calculations on cluster models.

Current-induced Orbital and Spin Magnetizations in Crystals with Helical Structure

We theoretically show that in a crystal with a helical lattice structure, orbital and spin magnetizations along a helical axis are induced by an electric current along the helical axis. We propose a simple tight-binding model for calculations, and the results can be generalized to any helical crystals. The induced magnetizations are opposite for right-handed and left-handed helices. The current-induced spin magnetization along the helical axis comes from a radial spin texture on the Fermi surface. This is in sharp contrast to Rashba systems where the induced spin magnetization is perpendicular to the applied current.

Spin-orbit-induced circulating currents in a semiconductor nanostructure

Circulating orbital currents produced by the spin-orbit interaction for a single electron spin in a quantum dot are explicitly evaluated at zero magnetic field, along with their effect on the total magnetic moment (spin and orbital) of the electron spin. The currents are dominated by coherent superpositions of the conduction and valence envelope functions of the electronic state, are smoothly varying within the quantum dot, and are peaked roughly halfway between the dot center and edge. Thus the spatial structure of the spin contribution to the magnetic moment (which is peaked at the dot center) differs greatly from the spatial structure of the orbital contribution. Even when the spin and orbital magnetic moments cancel (for g=0) the spin can interact strongly with local magnetic fields, e.g., from other spins, which has implications for spin lifetimes and spin manipulation.

Berry phase theory of Dzyaloshinskii-Moriya interaction and spin-orbit torques

Recent experiments on current-induced domain-wall motion in chiral domain walls reveal important contributions both from spin-orbit torques (SOTs) and from the Dzyaloshinskii-Moriya interaction (DMI). We derive a Berry phase expression for the DMI and show that within this Berry phase theory DMI and SOTs are intimately related, in a way formally analogous to the relation between orbital magnetization (OM) and anomalous Hall effect (AHE). We introduce the concept of the twist torque moment, which probes the internal twist of wavepackets in chiral magnets in a similar way as the orbital moment probes the wavepacket's internal self-rotation. We propose to interpret the Berry phase theory of DMI as a theory of spiralization in analogy to the modern theory of OM. We show that the twist torque moment and the spiralization together give rise to a Berry phase governing the response of the SOT to thermal gradients, in analogy to the intrinsic anomalous Nernst effect. The Berry phase theory of DMI is computationally very efficient because it only needs the electronic structure of the collinear magnetic system as input. As an application of the formalism we compute the DMI in Co/Pt(111), O/Co/Pt(111) and Al/Co/Pt(111) magnetic bi- and trilayers and show that the DMI is highly anisotropic in these systems.

Topological magnetoelectric effect decay

We address the influence of realistic disorder and finite doping on the effective magnetic monopole induced near the surface of an ideal topological insulator (TI) by currents that flow in response to a suddenly introduced external electric charge. We show that when the longitudinal conductivity σ(xx)=g(e(2)/h)≠0, the apparent position of a magnetic monopole initially retreats from the TI surface at speed v(M)=αcg, where α is the fine structure constant and c is the speed of light. For the particular case of TI surface states described by a massive Dirac model, we further find that the temperature T=0 Hall currents vanish when the external potential is screened.

Orbital magnetization as a local property

The modern expressions for polarization P and orbital magnetization M are k-space integrals. But a genuine bulk property should also be expressible in r space, as an unambiguous function of the ground-state density matrix, "nearsighted" in insulators, independently of the boundary conditions--either periodic or open. While P--owing to its "quantum" indeterminacy--is not a bulk property in this sense, M is. We provide its r-space expression for any insulator, even with a nonzero Chern invariant. Simulations on a model Hamiltonian validate our theory.

Controllable chirality-induced geometrical Hall effect in a frustrated highly correlated metal

A current of electrons traversing a landscape of localized spins possessing non-coplanar magnetic order gains a geometrical (Berry) phase, which can lead to a Hall voltage independent of the spin-orbit coupling within the material-a geometrical Hall effect. Here we show that the highly correlated metal UCu(5) possesses an unusually large controllable geometrical Hall effect at T<1.2 K due to its frustration-induced magnetic order. The magnitude of the Hall response exceeds 20% of the ν=1 quantum Hall effect per atomic layer, which translates into an effective magnetic field of several hundred Tesla acting on the electrons. The existence of such a large geometric Hall response in UCu(5) opens a new field of enquiry into the importance of the role of frustration in highly correlated electron materials.

Theoretical prediction of a rotating magnon wave packet in ferromagnets

We theoretically show that the magnon wave packet has a rotational motion in two ways: a self-rotation and a motion along the boundary of the sample (edge current). They are similar to the cyclotron motion of electrons, but unlike electrons the magnons have no charge and the rotation is not due to the Lorentz force. These rotational motions are caused by the Berry phase in momentum space from the magnon band structure. Furthermore, the rotational motion of the magnon gives an additional correction term to the magnon Hall effect. We also discuss the Berry curvature effect in the classical limit of long-wavelength magnetostatic spin waves having macroscopic coherence length.