BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Sinisa Coh

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

ResearcherUnclaimed source record
cond-mat.mtrl-scicond-mat.mes-hallcond-mat.supr-concond-mat.str-elphysics.chem-phphysics.optics

Catalog footprint

What is connected

29works
6topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

29 published item(s)

preprint2022arXiv

Magneto-optical Kerr spectra of gold induced by spin accumulation

We report the magneto optic Kerr effect (MOKE) angle of Au magnetically excited by spin accumulation. We perform time resolved polar MOKE measurements on Au/Co heterostructures. In our experiment, the ultrafast optical excitation of the Co drives spin accumulation into an adjacent Au layer. The spin accumulation, together with spin-orbit coupling, leads to non-zero terms in the off-diagonal conductivity tensor of Au, which we measure by recording the polarization and ellipticity of light reflected from the Au surface for photon energies between 1.3 and 3.1 eV. At energies near the interband transition threshold of Au, the Kerr rotation per A/m exceeds 1 urad. Typical values for Kerr rotation per moment in transition ferromagnetic metals like Ni are < 10 nrad per A/m, while predicted values for heavy metals like Pt or W are < 13 nrad per A/m. The exceptional sensitivity of the optical properties of Au to spin magnetic moments make Au to be an exceptionally sensitive optical magnetometer, with potential applications in the development of optospintronic technologies.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Nonadiabatic Born effective charges in metals and the Drude weight

In insulators, Born effective charges describe the electrical polarization induced by the displacement of individual atomic sublattices. Such a physical property is at first sight irrelevant for metals and doped semiconductors, where the macroscopic polarization is ill-defined. Here we show that, in clean conductors, going beyond the adiabatic approximation results in nonadiabatic Born effective charges that are well defined in the low-frequency limit. In addition, we find that the sublattice sum of the nonadiabatic Born effective charges does not vanish as it does in the insulating case, but instead is proportional to the Drude weight. We demonstrate these formal results with density functional perturbation theory calculations of Al, and electron-doped SnS$_2$ and SrTiO$_3$.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Large cross-polarized Raman signal in CrI$_3$: A first-principles study

We find unusually large cross-polarized (and anti-symmetric) Raman signature of A$_{\rm g}$ phonon mode in CrI$_3$, in agreement with experiments. The signal is present only when the following three effects are considered in concert: ferromagnetism on Cr atoms, spin-orbit interaction, and resonant effects. Somewhat surprisingly, we find that the relevant spin-orbit interaction potential originates from iodine atoms, despite magnetism being mostly on chromium atoms. We analyze the Raman signature as a function of magnetic order, the direction of the magnetic moment, energy and polarization of light used for Raman scattering, as well as carrier lifetime. In addition to a strong cross-polarized Raman signal, we also find unusually strong phonon modulated magneto-optical Kerr effect (MOKE) in CrI$_3$.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Comparison of GW band structure to semi-empirical approach for an FeSe monolayer

We present the G$_0$W$_0$ band structure, core levels, and deformation potential of monolayer FeSe in the paramagnetic phase based on a starting mean field of the Kohn Sham density functional theory (DFT) with the PBE functional. We find the GW correction increases the bandwidth of the states forming the $M$ pocket near the Fermi energy, while leaving the $Γ$ pocket roughly unchanged. We then compare the G$_0$W$_0$ quasiparticle band energies with the band structure from a simple empirical +A approach, which was recently proposed to capture the renormalization of the electron-phonon interaction going beyond DFT in FeSe, when used as a starting point in density functional perturbation theory (DFPT). We show that this empirical correction succeeds in approximating the GW non-local and dynamical self energy in monolayer FeSe and reproduces the GW band structure near the Fermi surface, the core energy levels, and the deformation potential (electron-phonon coupling).

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Hydrogen plasma favored modification of anatase TiO$_2$ (001) surface with desirable water splitting performance

We show that when TiO$_2$ anatase (001) is exposed to hydrogen plasma that the pristine surface termination becomes unfavorable to another, slightly modified, surface. On this modified surface the topmost TiO$_2$ layer is intact but out of registry with the bottom layers. Nevertheless, the modified surface has significantly improved ability to split water under exposure to sunlight. We show by explicit calculation of the water splitting reaction that the energy barrier that exists on a pristine surface is not present on the modified surface. The valence band maximum of the surface is raised relative to the pristine surface, which is a favorable way of adjusting the band gap in TiO$_2$ to the solar spectrum.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Signature of multilayer graphene strain-controlled domain walls in quantum Hall effect

Domain walls, topological defects that define the frontier between regions of different stacking in multilayer graphene, have proved to host exciting physics. The ability of tuning these topological defects in-situ in an electronic transport experiment brings a wealth of possibilities in terms of fundamental understanding of domain walls as well as for electronic applications. Here, we demonstrate through a MEMS (micro-electromechanical system) actuator and magnetoresistance measurements the effect of domain walls in multilayer graphene quantum Hall effect. Reversible and controlled uniaxial strain triggers these topological defects, manifested as new quantum Hall effect plateaus as well as a discrete and reversible modulation of the current across the device. Our findings are supported by theoretical calculations and constitute the first indication of the in-situ tuning of topological defects in multilayer graphene probed through electronic transport, opening the way to the use of reversible topological defects in electronic applications.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Antiferromagnetism enables electron-phonon coupling in iron-based superconductors

It is shown that a generic form of an antiferromagnetic wave function opens strong electron-phonon coupling channels in the iron-based superconductors. In the nonmagnetic state these channels exist locally on a single iron atom, but are canceled out between the two iron atoms in the primitive unit cell. Our findings are based on symmetry and the presence of an xz/yz Fermi surface near the M point and thus should be relevant for the known iron-based pnictide or chalcogenide superconductors.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Comparative study of high-Tc superconductivity in H3S and H3P

We report on a comparative study of the electronic structure, phonon spectra, and superconducting properties for recently discovered superconducting hydrides, H3S and H3P. While the electronic structures of these two materials are similar, there are notable changes in the phonon spectra and electron-phonon coupling. The low-frequency bond-bending modes are softened in H3P and their coupling to the electrons at the Fermi surface is enhanced relative to H3S. Nevertheless, coupling to the high-frequency modes is reduced so the resulting calculated superconducting transition temperature is reduced from ~166 K in H3S to ~76 K in H3P.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Inhibiting Klein tunneling in graphene p-n junction without an external magnetic field

We study by first-principles calculations a densely packed island of organic molecules (F4TCNQ) adsorbed on graphene. We find that with electron doping the island naturally forms a p-n junction in the graphene sheet. For example, a doping level of ~ 3 * 10^{13} electrons per cm^2 results in a p-n junction with 800 meV electrostatic potential barrier. Unlike in a conventional p-n junction in graphene, in the case of the junction formed by an adsorbed organic molecular island we expect that the Klein tunneling is inhibited, even without an applied external magnetic field. Here Klein tunneling is inhibited by the ferromagnetic order that spontaneously occurs in the molecular island upon doping. We estimate that the magnetic barrier in the graphene sheet is around 10 mT.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Proposal for a bulk material based on a monolayer FeSe on SrTiO3 high-temperature superconductor

Motivated by the high superconducting transition temperature of monolayer FeSe on SrTiO3, we propose a potential three-dimensional high-temperature superconductor superlattice FeSe-SrTiO3 and study its structural stability and electronic structure using density functional theory. We find that the binding energy between the FeSe and SrTiO3 layers is about 0.7 eV per (Fe2Se2) unit and that it saturates already within a single TiO2 atomic layer of SrTiO3. In addition we analyzed the dynamical stability of the superlattice and compared it to the case of bulk SrTiO3.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Why do Collapsed Carbon Nanotubes Twist?

We study the collapsing and subsequent spontaneous twisting of a carbon nanotube by in-situ transmission electron microscopy. A custom-sized nanotube is first created in the microscope by selectively extracting shells from a parent multi-wall tube. The few-wall, large-diameter daughter nanotube is driven to collapse via mechanical stimulation, after which the ribbon-like collapsed tube spontaneously twists along its long axis. In-situ diffraction experiments fully characterize the uncollapsed and collapsed tubes. From the experimental observations and associated theoretical analysis, the origin of the twisting is determined to be compressive strain due to charge imbalance.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Automated construction of maximally localized Wannier functions: Optimized projection functions method

Maximally localized Wannier functions are widely used in electronic structure theory for analyses of bonding, electric polarization, orbital magnetization, and for interpolation. The state of the art method for their construction is based on the method of Marzari and Vanderbilt. One of the practical difficulties of this method is guessing functions (initial projections) that approximate the final Wannier functions. Here we present an approach based on optimized projection functions that can construct maximally localized Wannier functions without a guess. We describe and demonstrate this approach on several realistic examples.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Calculation of the specific heat of optimally K-doped BaFe$_2$As$_2$

The calculated specific heat of optimally K-doped BaFe$_2$As$_2$ in density functional theory is about five times smaller than that found in the experiment. We report that by adjusting the potential on the iron atom to be slightly more repulsive for electrons improves the calculated heat capacity as well as the structural, magnetic, and electronic properties of Ba$_{0.6}$K$_{0.4}$Fe$_2$As$_2$. Applying the same correction to the antiferromagnetic state, we find that the electron-phonon coupling is strongly enhanced.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Intrinsic surface dipole in topological insulators

We calculate the local density of states of two prototypical topological insulators (Bi$_2$Se$_3$ and Bi$_2$Te$_2$Se) as a function of distance from the surface within density functional theory. We find that, in the absence of disorder or doping, there is a 2 nm thick surface dipole the origin of which is the occupation of the topological surface states above the Dirac point. As a consequence, the bottom of the conduction band is bent upward by about 75 meV near the surface, and there is a hump-like feature associated with the top of the valence band. We expect that band bending will occur in all pristine topological insulators as long as the Fermi level does not cross the Dirac point. Our results show that topological insulators are intrinsic Schottky barrier solar cells.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Large electron-phonon interactions from FeSe phonons in a monolayer

We show that electron-phonon coupling can induce strong electron pairing in an FeSe monolayer on a SrTiO3 substrate (experimental indications for superconducting Tc are between 65 and 109K). The role of the SrTiO3 substrate in increasing the coupling is two-fold. First, the interaction of the FeSe and TiO2 terminated face of SrTiO3 prevents the FeSe monolayer from undergoing a shear-type (orthorhombic, nematic) structural phase transition. Second, the substrate allows an anti-ferromagnetic ground state of FeSe which opens electron-phonon coupling channels within the monolayer that are prevented by symmetry in the non-magnetic phase. The spectral function for the electron-phonon coupling (alpha2F) in our calculations agrees well with inelastic tunneling data.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Atomically perfect torn graphene edges and their reversible reconstruction

The atomic structure of graphene edges is critical in determining the electrical, magnetic, and chemical properties of truncated graphene structures, notably nanoribbons. Unfortunately, graphene edges are typically far from ideal and suffer from atomic-scale defects, structural distortion, and unintended chemical functionalization, leading to unpredictable properties. Here we report that graphene edges fabricated by electron-beam-initiated mechanical rupture or tearing in high vacuum are clean and largely atomically perfect, oriented in either the armchair or zigzag direction. Via aberration-corrected transmission electron microscopy, we demonstrate reversible and extended pentagon-heptagon (5-7) reconstruction at zigzag edges, and explore experimentally and theoretically the dynamics of the transitions between configuration states. Good theoretical-experimental agreement is found for the flipping rates between 5-7 and 6-6 zigzag edge states. Our study demonstrates that simple ripping is remarkably effective in producing atomically clean, ideal terminations thus providing a valuable tool for realizing atomically-tailored graphene and facilitating meaningful experimental study.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Canonical magnetic insulators with isotropic magnetoelectric coupling

We have performed an exhaustive representation-theory-based search for the simplest structures allowing isotropic magnetoelectric coupling. We find 30 such structures, all sharing a common pattern of atomic displacements in the direction of atomic magnetic moments. We focus on one of these 30 canonical structures and find that it is generically realized in a class of fractionally substituted pyrochlore compounds with an all-in-all-out magnetic order. Furthermore, we find that these substituted pyrochlore compounds have a substantial Chern-Simons orbital magnetoelectric component (θ=0.1--0.2). While this component is also formally present in strong Z_2 topological insulators (θ=π), its effects are observable there only if time-reversal symmetry is broken at the surface.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Corner states of topological fullerenes

The unusual electronic properties of the quantum spin Hall or Chern insulator become manifest in the form of robust edge states when samples with boundaries are studied. In this work, we ask if and how the topologically non-trivial electronic structure of these two-dimensional systems can be passed on to their zero-dimensional relatives, namely fullerenes or other closed-cage molecules. To address this question, we study Haldane's honeycomb lattice model on polyhedral nano-surfaces. We find that for sufficiently large surfaces characteristic corner states appear for parameters for which the planar model displays a quantized Hall effect. In the electronic structure, these corner states show up as in-gap modes which are well separated from the quasi-continuum of states. We discuss the role of finite size effects and how the coupling between the corner states lifts the degeneracy in a characteristic way determined by the combined Berry phases which leads to an effective magnetic monopole of charge 2 at the center of the nano-surface. Experimental implications for fullerenes in the large spin-orbit regime are also pointed out.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Imprint of transition metal d-orbitals on graphene Dirac cone

We investigate the influence of SiO2, Au, Ag, Cu, and Pt substrates on the Raman spectrum of graphene. Experiments reveal particularly strong modifications to the intensity, position, width, and shape of the Raman signal of graphene on platinum, compared to that of suspended graphene. The modifications strongly depend on the relative orientation of the graphene and platinum lattices. These observations are theoretically investigated and shown to originate from hybridization of electronic states in graphene and d-orbitals in platinum. It is expected that, quite generally, hybridization between graphene and any material with d-orbitals near the Fermi level will result in an imprint on the graphene Dirac cone which depends sensitively on the relative orientation of the respective lattices.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Intermolecular interactions and substrate effects for an adamantane monolayer on the Au(111) surface

We study theoretically and experimentally the infrared (IR) spectrum of an adamantane monolayer on a Au(111) surface. Using a new STM-based IR spectroscopy technique (IRSTM) we are able to measure both the nanoscale structure of an adamantane monolayer on Au(111) as well as its infrared spectrum, while DFT-based ab initio calculations allow us to interpret the microscopic vibrational dynamics revealed by our measurements. We find that the IR spectrum of an adamantane monolayer on Au(111) is substantially modified with respect to the gas-phase IR spectrum. The first modification is caused by the adamantane--adamantane interaction due to monolayer packing and it reduces the IR intensity of the 2912 cm$^{-1}$ peak (gas phase) by a factor of 3.5. The second modification originates from the adamantane--gold interaction and it increases the IR intensity of the 2938 cm$^{-1}$ peak (gas phase) by a factor of 2.6, and reduces its frequency by 276 cm$^{-1}$. We expect that the techniques described here can be used for an independent estimate of substrate effects and intermolecular interactions in other diamondoid molecules, and for other metallic substrates.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Surface Atom Motion to Move Iron Nanocrystals through Constrictions in Carbon Nanotubes under the Action of an Electric Current

Under the application of electrical currents, metal nanocrystals inside carbon nanotubes can be bodily transported. We examine experimentally and theoretically how an iron nanocrystal can pass through a constriction in the carbon nanotube with a smaller cross-sectional area than the nanocrystal itself. Remarkably, through in situ transmission electron imaging and diffraction, we find that, while passing through a constriction, the nanocrystal remains largely solid and crystalline and the carbon nanotube is unaffected. We account for this behavior by a pattern of iron atom motion and rearrangement on the surface of the nanocrystal. The nanocrystal motion can be described with a model whose parameters are nearly independent of the nanocrystal length, area, temperature, and electromigration force magnitude. We predict that metal nanocrystals can move through complex geometries and constrictions, with implications for both nanomechanics and tunable synthesis of metal nanoparticles.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Theoretical study of solid iron nanocrystal movement inside a carbon nanotube

We use a first-principles based kinetic Monte Carlo simulation to study the movement of a solid iron nanocrystal inside a carbon nanotube driven by the electrical current. The origin of the iron nanocrystal movement is the electromigration force. Even though the iron nanocrystal appears to be moving as a whole, we find that the core atoms of the nanocrystal is completely stationary, and only the surface atoms are moving. Movement in the contact region with the carbon nanotube is driven by electromigration forces, and the movement on the remaining surfaces is driven by diffusion. Results of our calculations also provide a simple model which can predict the center of mass speed of the iron nanocrystal over a wide range of parameters. We find both qualitative and quantitative agreement of the iron nanocrystal center of mass speed with experimental data.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Theory of the Raman spectrum of rotated double-layer graphene

We study theoretically the Raman spectrum of the rotated double-layer graphene, consisting of two graphene layers rotated with respect to each other by an arbitrary angle θ. We find a relatively simple dependence of the Raman G peak intensity on the angle θ. On the other hand, the Raman 2D peak position, intensity, and width show a much more complicated dependence on the angle θ. We account for all of these effects, including dependence on the incoming photon energy, in good agreement with the experimental data. We find that it is sufficient to include the interaction between the graphene layers on the electronic degrees of freedom (resulting in the occurrence of van Hove singularities in the density of states). We assume that the phonon degrees of freedom are unaffected by the interaction between the layers. Furthermore, we decompose the Raman 2D peak into two components having very different linewidths; these widths are almost independent of the angle θ. The change in the intensity and the peak position of one of these two components gives insight into the dependence of the overall Raman 2D peak features as a function of the angle θ. Additionally, we show regions in the phonon spectrum giving rise to the Raman 2D peak signal. This work provides an insight into the interplay between the mechanical degree of freedom (angle θ) and the electronic degrees of freedom (singularities in the density of states) in rotated double-layer graphene. Additionally, this work provides a way to establish experimentally the value of the rotation angle θusing Raman spectroscopy measurement. This procedure becomes even more robust if one repeats the Raman spectroscopy measurement with a different incoming photon energy.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Full magnetoelectric response of Cr2O3 from first principles

The linear magnetoelectric response of Cr2O3 at zero temperature is calculated from first principles by tracking the change in magnetization under a macroscopic electric field. Both the spin and the orbital contributions to the induced magnetization are computed, and in each case the response is decomposed into lattice and electronic parts. We find that the transverse response is dominated by the spin-lattice and spin-electronic contributions, whose calculated values are consistent with static and optical magnetoelectric measurements. In the case of the longitudinal response, orbital contributions dominate over spin contributions, but the net calculated longitudinal response remains much smaller than the experimentally measured one at low temperatures. We also discuss the absolute sign of the magnetoelectric coupling in the two time-reversed magnetic domains of Cr2O3.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Raman spectroscopy study of rotated double-layer graphene: Misorientation-angle dependence of electronic structure

We present a systematic Raman study of unconventionally-stacked double-layer graphene, and find that the spectrum strongly depends on the relative rotation angle between layers. Rotation-dependent trends in the position, width and intensity of graphene 2D and G peaks are experimentally established and accounted for theoretically. Our theoretical analysis reveals that changes in electronic band structure due to the interlayer interaction, such as rotational-angle dependent Van Hove singularities, are responsible for the observed spectra features. Our combined experimental and theoretical study provides a deeper understanding of the electronic band structure of rotated double-layer graphene, and leads to a practical way to identify and analyze rotation angles of misoriented double-layer graphene.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Large isosymmetric reorientation of oxygen octahedra rotation axes in epitaxially strained perovskites

Using first-principles density functional theory calculations, we discover an anomalously large bi-axial strain-induced octahedral rotation axis reorientation in orthorhombic perovskites with tendency towards rhombohedral symmetry. The transition between crystallographically equivalent (isosymmetric) structures with different octahedral rotation magnitudes originates from strong strain--octahedral rotation coupling available to perovskites and the energetic hierarchy among competing octahedral tilt patterns. By elucidating these criteria, we suggest many functional perovskites would exhibit the transition in thin film form, thus offering a new landscape in which to tailor highly anisotropic electronic responses.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Chern-Simons orbital magnetoelectric coupling in generic insulators

We present a Wannier-based method to calculate the Chern-Simons orbital magnetoelectric coupling in the framework of first-principles density-functional theory. In view of recent developments in connection with strong Z2 topological insulators, we anticipate that the Chern-Simons contribution to the magnetoelectric coupling could, in special cases, be as large or larger than the total magnetoelectric coupling in known magnetoelectrics like Cr2O3. The results of our calculations for the ordinary magnetoelectrics Cr2O3, BiFeO3 and GdAlO3 confirm that the Chern-Simons contribution is quite small in these cases. On the other hand, we show that if the spatial inversion and time-reversal symmetries of the Z2 topological insulator Bi2Se3 are broken by hand, large induced changes appear in the Chern-Simons magnetoelectric coupling.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Si-compatible candidates for high-K dielectrics with the Pbnm perovskite structure

We analyze both experimentally (where possible) and theoretically from first-principles the dielectric tensor components and crystal structure of five classes of Pbnm perovskites. All of these materials are believed to be stable on silicon and are therefore promising candidates for high-K dielectrics. We also analyze the structure of these materials with various simple models, decompose the lattice contribution to the dielectric tensor into force constant matrix eigenmode contributions, explore a peculiar correlation between structural and dielectric anisotropies in these compounds and give phonon frequencies and infrared activities of those modes that are infrared-active. We find that CaZrO_3, SrZrO_3, LaHoO_3, and LaYO_3 are among the most promising candidates for high-K dielectrics among the compounds we considered.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Theory of orbital magnetoelectric response

We extend the recently-developed theory of bulk orbital magnetization to finite electric fields, and use it to calculate the orbital magnetoelectric response of periodic insulators. Working in the independent-particle framework, we find that the finite-field orbital magnetization can be written as a sum of three gauge-invariant contributions, one of which has no counterpart at zero field. The extra contribution is collinear with and explicitly dependent on the electric field. The expression for the orbital magnetization is suitable for first-principles implementations, allowing to calculate the magnetoelectric response coefficients by numerical differentiation. Alternatively, perturbation-theory techniques may be used, and for that purpose we derive an expression directly for the linear magnetoelectric tensor by taking the first field-derivative analytically. Two types of terms are obtained. One, the `Chern-Simons' term, depends only on the unperturbed occupied orbitals and is purely isotropic. The other, `Kubo' term, involves the first-order change in the orbitals and gives isotropic as well as anisotropic contributions to the response. In ordinary magnetoelectric insulators both terms are generally present, while in strong Z2 topological insulators only the former type is allowed, and is quantized. In order to validate the theory we have calculated under periodic boundary conditions the linear magnetoelectric coupling for a 3-D tight-binding model of an ordinary magnetoelectric insulator, using both the finite-field and perturbation-theory expressions. The results are in excellent agreement with calculations on bounded samples.

0 citations0 reviewsNo code linkedOpen access