Source author record

Mark van Schilfgaarde

Mark van Schilfgaarde 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

27works
9topics
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

27 published item(s)

preprint2025arXiv

Interplay of Quasi-Quantum Hall Effect and Coulomb Disorder in Semimetals

Low carrier densities in topological semimetals (TSMs) enable the exploration of novel magnetotransport in the quantum limit (QL). Recent findings consistent with 3D quasi-quantum Hall effect (QQHE) have positioned TSMs as promising platforms for exploring 3D quantum Hall transport, but the lack of tunability in the Fermi level has thus far limited the ability to observe a QQHE signal. Here, we tune the defect concentrations in the Dirac semimetal Cd${}_3$As${}_2$ to achieve ultra-low carrier concentrations at 2 K around $2.9\times10^{16}$cm${}^{-3}$, giving way to QQHE signal at modest fields near 10 T. At low carrier densities, where QQHE is most accessible, we find that clear QQHE is obscured by a carrier density dependent background originating from Coulomb disorder from charged point defects and Landau level broadening. Our results highlight the interplay between QQHE and Coulomb disorder, demonstrating that clear observation of QQHE in TSMs intricately depends on Fermi level and disorder magnitudes. We find that Coulomb disorder, as theoretically predicted, is an essential ingredient for understanding the magnetoresistivity for a spectrum of Fermi levels in Cd${}_3$As${}_2$, anchoring the role of defects and charged disorder in TSM applications. We discuss future constraints and opportunities in exploring 3D QQHE and quantum Hall effects in TSMs.

preprint2022arXiv

Electronic and optical properties of crystalline nitrogen versus black phosphorus: A comparative first-principles study

Crystalline black nitrogen (BN) is an allotrope of nitrogen with the black phosphorus (BP) structure recently synthesized at high pressure by two independent research groups [Ji et al., Sci. Adv. 6, eaba9206 (2020); Laniel et al., Phys. Rev. Lett. 124, 216001 (2020)]. Here, we present a systematic study of the electronic and optical properties of BN focusing on its comparison with BP. To this end, we use the state-of-the-art quasiparticle self-consistent $GW$ approach with vertex corrections in both the electronic and optical channels. Despite many similarities, the properties of BN are found to be considerably different. Unlike BP, BN exhibits a larger optical gap (2.5 vs 0.26 eV), making BN transparent in the visible spectral region with a highly anisotropic optical response. This difference can be primarily attributed to a considerably reduced dielectric screening in BN, leading to enhancement of the effective Coulomb interaction. Despite relatively strong Coulomb interaction, exciton formation is largely suppressed in both materials. Our analysis of the elastic properties shows exceptionally high stiffness of BN, comparable to that of diamond.

preprint2022arXiv

Real space representation of the quasiparticle self-consistent $GW$ self-energy and its application to defect calculations

The quasiparticle self-consistent QS$GW$ approach incorporates the corrections of the quasiparticle energies from their Kohn-Sham density functional theory (DFT) eigenvalues by means of an energy independent and Hermitian self-energy matrix usually given in the basis set of the DFT eigenstates. By expanding these into an atom-centered basis set (specifically here the linearized muffin-tin orbitals) a real space representation of the self-energy corrections becomes possible. We show that this representation is relatively short-ranged. This offers new opportunities to construct the self-energy of a complex system from parts of the system by a cut-and-paste method. Specifically for a point defect, represented in a large supercell, the self-eneregy can be constructed from those of the host and a smaller defect containing cell. The self-energy of the periodic host can be constructed simply from a $GW$ calculation for the primitive cell. We show for the case of the As$_\mathrm{Ga}$ in GaAs that the defect part can already be well represented by a minimal 8 atom cell and allows us to construct the self-energy for a 64 cell in good agreement with direct QS$GW$ calculations for the large cell. Using this approach to an even larger 216 atom cell shows the defect band approaches an isolated defect level. The calculations also allow to identify a second defect band which appears as a resonance near the conduction band minimum. The results on the extracted defect levels agree well with Green's function calculations for an isolated defect and with experimental data.

preprint2022arXiv

Role of nematicity in controlling spin fluctuations and superconducting Tc in bulk FeSe

FeSe undergoes a transition from a tetragonal to a slightly orthorhombic phase at 90\,K, and becomes a superconductor below 8\,K. The orthorhombic phase is sometimes called a nematic phase because quantum oscillation, neutron, and other measurements detect a significant asymmetry in $x$ and $y$. How nematicity affects superconductivity has recently become a matter of intense speculation. Here we employ an advanced \emph{ab-initio} Green's function description of superconductivity and show that bulk tetragonal FeSe would, in principle, superconduct with almost the same T$_{c}$ as the nematic phase. The mechanism driving the observed nematicity is not yet understood. Since the present theory underestimates nematicity, we simulate the full nematic asymmetry by artificially enhancing the orthorhombic distortion. For benchmarking, we compare theoretical spin susceptibilities against experimentally observed data over all energies and relevant momenta. When the orthorhombic distortion is adjusted to correlate with observed nematicity in spin susceptibility, the enhanced nematicity causes spectral weight redistribution in the Fe-3d$_{xz}$ and d$_{yz}$ orbitals, but it leads to at most 10-15$\%$ increment in T$_{c}$. This is because the d$_{xy}$ orbital always remains the most strongly correlated and provides most of the source of the superconducting glue. Nematicity suppresses the density of states at Fermi level; nevertheless T$_{c}$ increases, in contradiction to both BCS and BEC theories. We show how the increase is connected to the structure of the particle-particle vertex. Our results suggest while nematicity may be intrinsic property of bulk FeSe, is not the primary force driving the superconducting pairing.

preprint2021arXiv

Excitons in Bulk and Layered Chromium Tri-Halides: From Frenkel to the Wannier-Mott Limit

Excitons with large binding energies $\sim$2-3 eV in CrX$_{3}$ are historically characterized as being localized (Frenkel) excitons that emerge from the atomic $d{-}d$ transitions between the Cr-3$d$-$t_{2g}$ and $e_{g}$ orbitals. The argument has gathered strength in recent years as the excitons in recently made monolayers are found at almost the same energies as the bulk. The Laporte rule, which restricts such parity forbidden atomic transitions, can relax if, at least, one element is present: spin-orbit coupling, odd-parity phonons or Jahn-Teller distortion. While what can be classified as a purely Frenkel exciton is a matter of definition, we show using an advanced first principles parameter-free approach that these excitons in CrX$_{3}$, in both its bulk and monolayer variants, have band-origin and do not require the relaxation of Laporte rule as a fundamental principle. We show that, the character of these excitons is mostly determined by the Cr-$d$ orbital manifold, nevertheless, they appear only as a consequence of X-p states hybridizing with the Cr-$d$. The hybridization enhances as the halogen atom becomes heavier, bringing the X-$p$ states closer to the Cr-$d$ states in the sequence Cl{\textrightarrow}Br{\textrightarrow}I, with an attendant increase in exciton intensity and decrease in binding energy. By applying a range of different kinds of perturbations, we show that, moderate changes to the two-particle Hamiltonian that essentially modifies the Cr-$d$-X-$p$ hybridization, can alter both the intensities and positions of the exciton peaks. A detailed analysis of several deep lying excitons, with and without strain, reveals that the exciton is most Frenkel like in CrCl$_{3}$ and acquires mixed Frenkel-Wannier character in CrI$_{3}$.

preprint2021arXiv

First principles vs second principles: Role of charge self-consistency in strongly correlated systems

First principles approaches have been successful in solving many-body Hamiltonians for real materials to an extent when correlations are weak or moderate. As the electronic correlations become stronger often embedding methods based on first principles approaches are used to better treat the correlations by solving a suitably chosen many-body Hamiltonian with a higher level theory. Such combined methods are often referred to as second principles approaches. At such level of the theory the self energy, i.e. the functional that embodies the stronger electronic correlations, is either a function of energy or momentum or both. The success of such theories is commonly measured by the quality of the self energy functional. However, self-consistency in the self-energy should, in principle, also change the real space charge distribution in a correlated material and be able to modify the electronic eigenfunctions, which is often undermined in second principles approaches. Here we study the impact of charge self-consistency within two example cases: TiSe$_{2}$, a three-dimensional charge-density-wave candidate material, and CrBr$_{3}$, a two-dimensional ferromagnet, and show how real space charge re-distribution due to correlation effects taken into account within a first principles Green's function based many-body perturbative approach is key in driving qualitative changes to the final electronic structure of these materials.

preprint2021arXiv

Supercurrent decay in ballistic magnetic Josephson junctions

We investigate transport properties of ballistic magnetic Josephson junctions and establish that suppression of supercurrent is an intrinsic property of the junctions, even in absence of disorder. By studying the role of ferromagnet thickness, magnetization, and crystal orientation we show how the supercurrent decays exponentially with thickness and identify two mechanisms responsible for the effect: (i) large exchange splitting may gap out minority or majority carriers leading to the suppression of Andreev reflection in the junction, (ii) loss of synchronization between different modes due to the significant dispersion of the quasiparticle velocity with the transverse momentum. Our results for Nb/Ni/Nb junctions are in good agreement with recent experimental studies. Our approach combines density functional theory and Bogoliubov-de Gennes model and opens a path for material composition optimization in magnetic Josephson junctions and superconducting magnetic spin valves.

preprint2020arXiv

Controlling T$_c$ through band structure and correlation engineering in collapsed and uncollapsed phases of iron arsenides

Recent observations of selective emergence (suppression) of superconductivity in the uncollapsed (collapsed) tetragonal phase of LaFe$_2$As$_2$ has rekindled interest in understanding what features of the band structure control the superconducting T$_c$. We show that the proximity of the narrow Fe-d$_{xy}$ state to the Fermi energy emerges as the primary factor. In the uncollapsed phase this state is at the Fermi energy, and is most strongly correlated and source of enhanced scattering in both single and two particle channels. The resulting intense and broad low energy spin fluctuations suppress magnetic ordering and simultaneously provide glue for Cooper pair formation. In the collapsed tetragonal phase, the d$_{xy}$ state is driven far below the Fermi energy, which suppresses the low-energy scattering and blocks superconductivity. A similar source of broad spin excitation appears in uncollapsed and collapsed phases of CaFe$_{2}$As$_{2}$. This suggests controlling coherence provides a way to engineer T$_c$ in unconventional superconductors primarily mediated through spin fluctuations.

preprint2020arXiv

Electron-Phonon-Driven Three-Dimensional Metallicity in an Insulating Cuprate

The role of the crystal lattice for the electronic properties of cuprates and other high-temperature superconductors remains controversial despite decades of theoretical and experimental efforts. While the paradigm of strong electronic correlations suggests a purely electronic mechanism behind the insulator-to-metal transition, recently the mutual enhancement of the electron-electron and the electron-phonon interaction and its relevance to the formation of the ordered phases have also been emphasized. Here, we combine polarization-resolved ultrafast optical spectroscopy and state-of-the-art dynamical mean-field theory to show the importance of the crystal lattice in the breakdown of the correlated insulating state in an archetypal undoped cuprate. We identify signatures of electron-phonon coupling to specific fully-symmetric optical modes during the build-up of a three-dimensional metallic state that follows charge photodoping. Calculations for coherently displaced crystal structures along the relevant phonon coordinates indicate that the insulating state is remarkably unstable toward metallization despite the seemingly large charge-transfer energy scale. This hitherto unobserved insulator-to-metal transition mediated by fully-symmetric lattice modes can find extensive application in a plethora of correlated solids.

preprint2020arXiv

Electronic structure correspondence of singlet-triplet scale separation in strained Sr2RuO4

At a temperature of roughly 1\,K, \ce{Sr2RuO4} undergoes a transition from a normal Fermi liquid to a superconducting phase. Even while the former is relatively simple and well understood, the superconducting state is not even after 25 years of study. More recently it has been found that critical temperatures can be enhanced by application of uniaxial strain, up to a critical strain, after which it falls off. In this work, we take an `instability' approach and seek for divergences in susceptibilities. This provides an unbiased way to distinguish tendencies to competing ground states. We show that in the unstrained compound the singlet and triplet instabilities of the normal Fermi liquid phase are closely spaced. Under uniaxial strain electrons residing on all orbitals contributing to the Fermiology become more coherent while the electrons of Ru-$d_{xy}$ character become heavier and electrons of Ru-$d_{xz,yz}$ characters become lighter. In the process, Im\,$χ(\mathbf{q},ω)$ increases rapidly around the incommensurate vector $\mathbf{q}{=}(0.3,0.3,0)2π/a$ while it gets suppressed at all other commensurate vectors, in particular at $q{=}0$, which is essential for spin-triplet superconductivity. Thus the triplet superconducting instability remains the lagging instability of the system and the singlet instability enhances under strain, leading to a large energy-scale separation between these competing instabilities. At large strain an instability to a spin density wave overtakes the superconducting one. The analysis relies on a high-fidelity, \emph{ab initio} description of the one-particle properties and two-particle susceptibilities, based on the Quasiparticle Self-Consistent \emph{GW} approximation augmented by Dynamical Mean Field theory. This approach is described and its high fidelity confirmed by comparing to observed one- and two-particle properties.

preprint2020arXiv

Interplay between band structure and Hund's correlation to increase T$_{c}$ in FeSe

FeSe is classed as a Hund's metal, with a multiplicity of $d$ bands near the Fermi level. Correlations in Hund's metals mostly originate from the exchange parameter \emph{J}, which can drive a strong orbital selectivity in the correlations. The Fe-chalcogens are the most strongly correlated of the Fe-based superconductors, with $d_{xy}$ the most correlated orbital. Yet little is understood whether and how such correlations directly affect the superconducting instability in Hund's systems. By applying a recently developed high-fidelity \emph{ab initio} theory, we show explicitly the connections between correlations in $d_{xy}$ and the superconducting critical temperature $T_{c}$. Starting from the \emph{ab initio} results as a reference, we consider various kinds of excursions in parameter space around the reference to determine what controls $T_{c}$. We show small excursions in $J$ can cause colossal changes in $T_{c}$. Additionally we consider changes in hopping by varying the Fe-Se bond length in bulk, in the free standing monolayer M-FeSe, and M-FeSe on a SrTiO$_{3}$ substrate (M-FeSe/STO). The twin conditions of proximity of the $d_{xy}$ state to the Fermi energy, and the strength of $J$ emerge as the primary criteria for incoherent spectral response and enhanced single- and two-particle scattering that in turn controls $T_{c}$. Using constrained RPA, we show further that FeSe in monolayer form (M-FeSe) provides a natural mechanism to enhance $J$. We explain why M-FeSe/STO has a high $T_{c}$, whereas M-FeSe in isolation should not. Our study opens a paradigm for a unified understanding what controls $T_{c}$ in bulk, layers, and interfaces of Hund's metals by hole pocket and electron screening cloud engineering.

preprint2019arXiv

Direct comparison of many-body methods for realistic electronic Hamiltonians

A large collaboration carefully benchmarks 20 first principles many-body electronic structure methods on a test set of 7 transition metal atoms, and their ions and monoxides. Good agreement is attained between the 3 systematically converged methods, resulting in experiment-free reference values. These reference values are used to assess the accuracy of modern emerging and scalable approaches to the many-electron problem. The most accurate methods obtain energies indistinguishable from experimental results, with the agreement mainly limited by the experimental uncertainties. Comparison between methods enables a unique perspective on calculations of many-body systems of electrons.

preprint2019arXiv

Possible phonon-induced electronic bi-stability in VO$_2$ for ultrafast memory at room temperature

VO$_{2}$ is a model material system which exhibits a metal to insulator transition at 67$^\circ$C. This holds potential for future ultrafast switching in memory devices, but typically requires a purely electronic process to avoid the slow lattice response. The role of lattice vibrations is thus important, but it is not well understood and it has been a long-standing source of controversy. We use a combination of ultrafast spectroscopy and ab initio quantum calculations to unveil the mechanism responsible for the transition. We identify an atypical Peierls vibrational mode which acts as a trigger for the transition. This rules out the long standing paradigm of a purely electronic Mott transition in VO$_{2}$; however, we found a new electron-phonon pathway for a purely reversible electronic transition in a true bi-stable fashion under specific conditions. This transition is very atypical, as it involves purely charge-like excitations and requires only small nuclear displacement. Our findings will prompt the design of future ultrafast electro-resistive non-volatile memory devices.

preprint2019arXiv

Questaal: a package of electronic structure methods based on the linear muffin-tin orbital technique

This paper summarises the theory and functionality behind Questaal, an open-source suite of codes for calculating the electronic structure and related properties of materials from first principles. The formalism of the linearised muffin-tin orbital (LMTO) method is revisited in detail and developed further by the introduction of short-ranged tight-binding basis functions for full-potential calculations. The LMTO method is presented in both Green's function and wave function formulations for bulk and layered systems. The suite's full-potential LMTO code uses a sophisticated basis and augmentation method that allows an efficient and precise solution to the band problem at different levels of theory, most importantly density functional theory, LDA+U, quasi-particle self-consistent GW and combinations of these with dynamical mean field theory. This paper details the technical and theoretical bases of these methods, their implementation in Questaal, and provides an overview of the code's design and capabilities.

preprint2016arXiv

A photon ratchet route to high-efficiency hybrid halide perovskite intermediate band solar cells

The spin-split indirect bandgap in hybrid-halide perovskites provides a momentum-space realisation of a photon-ratchet intermediate band. Excited electrons thermalise to recombination-protected Rashba pockets offset in momentum space, building up the charge density to have sufficient flux to the higher lying conduction band. This effect could be used to form an intrinsic intermediate band solar cell with efficiencies beyond the Shockley-Queisser limit if a selective low-electron affinity contact can be made to the higher conduction state. This concept is supported by analysis of the many-body electronic structure. Production of above-bandgap voltages under illumination would affirm the physical mechanism proposed here.

preprint2016arXiv

Band-filling effect on magnetic anisotropy using a Green's function method

We use an analytical model to describe the magnetocrystalline anisotropy energy (MAE) in solids as a function of band filling. The MAE is evaluated in second-order perturbation theory, which makes it possible to decompose the MAE into a sum of transitions between occupied and unoccupied pairs. The model enables us to characterize the MAE as a sum of contributions from different, often competing terms. The nitridometalates Li$_{2}$[(Li$_{1-x}$T$_{x}$)N], with $T$=Mn, Fe, Co, Ni, provide a system where the model is very effective because atomic like orbital characters are preserved and the decomposition is fairly clean. Model results are also compared against MAE evaluated directly from first-principles calculations for this system. Good qualitative agreement is found.

preprint2016arXiv

First-principles treatment of Mott insulators: linearized QSGW+DMFT approach

The theoretical understanding of emergent phenomena in quantum materials is one of the greatest challenges in condensed matter physics. In contrast to simple materials such as noble metals and semiconductors, macroscopic properties of quantum materials cannot be predicted by the properties of individual electrons. One of the examples of scientific importance is strongly correlated electron system. Neither localized nor itinerant behaviors of electrons in partially filled 3d, 4f, and 5f orbitals give rise to rich physics such as Mott insulators, high-temperature superconductors, and superior thermoelectricity, but hinder quantitative understanding of low-lying excitation spectrum. Here we present a new first-principles approach to strongly correlated solids. It is Q4 based on a combination of the quasiparticle self-consistent GW approximation and the dynamical mean-field theory. The sole input in this method is the projector to the set of correlated orbitals for which all local Feynman graphs are being evaluated. For that purpose, we choose very localized quasiatomic orbitals spanning large energy window, which contains most strongly hybridized bands, as well as upper and lower Hubbard bands. The self-consistency is carried out on the Matsubara axis. This method enables the first-principles study of Mott insulators in both their paramagnetic and antiferromagnetic phases. We illustrate the method on the archetypical charge transfer correlated insulators La2CuO4 and NiO, and obtain spectral properties and magnetic moments in good agreement with experiments.

preprint2016arXiv

Relativistic origin of slow electron-hole recombination in hybrid halide perovskite solar cells

The hybrid perovskite CH3NH3PbI3 (MAPI) exhibits long minority-carrier lifetimes and diffusion lengths. We show that slow recombination originates from a spin-split indirect-gap. Large internal electric fields act on spin-orbit-coupled band extrema, shifting band-edges to inequivalent wavevectors, making the fundamental gap indirect. From a description of photoluminescence within the quasiparticle self-consistent GW (QSGW) approximation for MAPI, CdTe and GaAs, we predict carrier lifetime as a function of light intensity and temperature. At operating conditions we find radiative recombination in MAPI is reduced by a factor of more than 350 compared to direct gap behavior. The indirect gap is retained with dynamic disorder.

preprint2015arXiv

Direct Method for Calculating Temperature-Dependent Transport Properties

We show how temperature-induced disorder can be combined in a direct way with first-principles scattering theory to study diffusive transport in real materials. Excellent (good) agreement with experiment is found for the resistivity of Cu, Pd, Pt (and Fe) when lattice (and spin) disorder are calculated from first principles. For Fe, the agreement with experiment is limited by how well the magnetization (of itinerant ferromagnets) can be calculated as a function of temperature. By introducing a simple Debye-like model of spin disorder parameterized to reproduce the experimental magnetization, the temperature dependence of the average resistivity, the anisotropic magnetoresistance and the spin polarization of a Ni$_{80}$Fe$_{20}$ alloy are calculated and found to be in good agreement with existing data. Extension of the method to complex, inhomogeneous materials as well as to the calculation of other finite-temperature physical properties within the adiabatic approximation is straightforward.

preprint2014arXiv

Atomistic origins of high-performance in hybrid halide perovskite solar cells

The performance of organometallic perovskite solar cells has rapidly surpassed that of both conventional dye-sensitised and organic photovoltaics. High power conversion efficiency can be realised in both mesoporous and thin-film device architectures. We address the origin of this success in the context of the materials chemistry and physics of the bulk perovskite as described by electronic structure calculations. In addition to the basic optoelectronic properties essential for an efficient photovoltaic device (spectrally suitable band gap, high optical absorption, low carrier effective masses), the materials are structurally and compositionally flexible. As we show, hybrid perovskites exhibit spontaneous electric polarisation; we also suggest ways in which this can be tuned through judicious choice of the organic cation. The presence of ferroelectric domains will result in internal junctions that may aid separation of photoexcited electron and hole pairs, and reduction of recombination through segregation of charge carriers. The combination of high dielectric constant and low effective mass promotes both Wannier-Mott exciton separation and effective ionisation of donor and acceptor defects. The photoferroic effect could be exploited in nanostructured films to generate a higher open circuit voltage and may contribute to the current-voltage hysteresis observed in perovskite solar cells.

preprint2014arXiv

Relativistic quasiparticle self-consistent electronic structure of hybrid halide perovskite photovoltaic absorbers

Solar cells based on a light absorbing layer of the organometal halide perovskite CH$_3$NH$_3$PbI$_3$ have recently reached 15% conversion efficiency, though how these materials work remains largely unknown. We analyse the electronic structure and optical properties within the quasiparticle self-consistent GW approximation. While this compound bears some similarity to conventional sp semiconductors, it also displays unique features. Quasiparticle self-consistency is essential for an accurate description of the band structure: bandgaps are much larger than what is predicted by the local density approximation (LDA) or GW based on the LDA. Valence band dispersions are modified in a very unusual manner. In addition, spin orbit coupling strongly modifies the band structure and gives rise to unconventional dispersion relations and a Dresselhaus splitting at the band edges. The average hole mass is small, which accounts for the long diffusion lengths recently observed. The surface ionisation potential (workfunction) is calculated to be 5.7 eV with respect to the vacuum level, explaining efficient carrier transfer to TiO$_2$ and Au electrical contacts.

preprint2013arXiv

Effects of alloying and strain on the magnetic properties of Fe$_{16}$N$_2$

The electronic structure and magnetic properties of pure and doped {Fe$_{16}$N$_2$} systems have been studied in the local-density (LDA) and quasiparticle self-consistent {\emph{GW}} approximations. The {\emph{GW}} magnetic moment of pure {Fe$_{16}$N$_2$} is somewhat larger compared to LDA but not anomalously large. The effects of doping on magnetic moment and exchange coupling were analyzed using the coherent potential approximation. The theoretical Curie temperature in pure {Fe$_{16}$N$_2$} is significantly higher than the measured value, which is attributed to the quality of available samples and the interpretation of experimental results. We found that different Fe sites contribute very differently to the magnetocrystalline anisotropy energy (MAE), which offers a way to increase MAE by small additions of Co or Ti. MAE also increases under tetragonal strain.

preprint2012arXiv

Doping-induced metallicity and coexistence of magnetic subsystems in {K$_{2}$Fe$_{4+x}$Se$_{5}$}

Electronic structure calculations are used to analyze the electronic and magnetic properties in {K$_{2}$Fe$_{4+x}$Se$_{5}$}. Fe atoms can be divided into two distinct groups. The $x{=}0$ (parent) compound forms an insulating, collinear, local moment phase with high N{é}el temperature. We show that large biquadratic exchange coupling and exchange-elastic interactions stabilize the magnetic order. For $x{>}0$ the additional Fe atoms fill vacancy sites. They form impurity bands for small $x$, which broaden as $x$ increases. They determine the states at the Fermi level and may be viewed as a magnetic subsystem separate from the host. Spin fluctuations are prevalent because magnetic interactions between the `defect' and the `parent' atoms are relatively weak, while chemical fluctuations are prevalent for low $x$. Fluctuations of either type leads to the formation of a weakly metallic state. The unusual coexistence of the two magnetic subsystems offers a new perspective as to how superconductivity and strong antiferromagnetism can coexist. We argue that spin fluctuations of the impurity subsystem share common features with the Fe-pnictide superconductors.

preprint2012arXiv

Effect of electron correlations on (001) Fe/MgO interfaces

We developed a parametrization of transmission probability that reliably captures essential elements of the tunneling process in magnetic tunnel junctions. The electronic structure of Fe/MgO system is calculated within the quasiparticle self-consistent GW approximation and used to evaluate transmission probability across (001) Fe/MgO interface. The transmission has a peak at +0.12 V, in excellent agreement with recent differential conductance measurements for electrodes with antiparallel spin. These findings confirm that the observed current-voltage characteristics are intrinsic to well defined (001) Fe/MgO interfaces, in contrast to previous predictions based on the local spin-density approximation, and also that many-body effects are important to realistically describe electron transport across well defined metal-insulator interfaces.

preprint2012arXiv

Spin excitations in K$_{2}$Fe$_{4+x}$Se$_{5}$: linear response approach

Using \emph{ab initio} linear response techniques we calculate spin wave spectra in K$_{2}$Fe$_{4+x}$Se$_{5}$, and find it to be in excellent agreement with a recent experiment. The spectrum can be alternatively described rather well by localized spin Hamiltonian restricted to first and second nearest neighbor couplings. We confirm that exchange coupling between nearest neighbor Fe magnetic moments is strongly anisotropic, and show directly that in the ideal system this anisotropy has itinerant nature which can be imitated by introducing higher order terms in effective localized spin Hamiltonian (biquadratic coupling). In the real system, structural relaxation provides an additional source of the exchange anisotropy of approximately the same magnitude. The dependence of spin wave spectra on filling of Fe vacancy sites is also discussed.

preprint2009arXiv

Prediction of large linear-in-k spin splitting for holes in the 2D GaAs/AlAs system

The spin-orbit interaction generally leads to spin splitting (SS) of electron and hole energy states in solids, a splitting that is characterized by a scaling with the wavevector $\bf k$. Whereas for {\it 3D bulk zincblende} solids the electron (heavy hole) SS exhibits a cubic (linear) scaling with $k$, in {\it 2D quantum-wells} the electron (heavy hole) SS is currently believed to have a mostly linear (cubic) scaling. Such expectations are based on using a small 3D envelope function basis set to describe 2D physics. By treating instead the 2D system explicitly in a multi-band many-body approach we discover a large linear scaling of hole states in 2D. This scaling emerges from hole bands coupling that would be unsuspected by the standard model that judges coupling by energy proximity. This discovery of a linear Dresselhaus k-scaling for holes in 2D implies a different understanding of hole-physics in low-dimensions.

preprint2008arXiv

Impact ionization rates for Si, GaAs, InAs, ZnS, and GaN in the $GW$ approximation

We present first-principles calculations of the impact ionization rate (IIR) in the $GW$ approximation ($GW$A) for semiconductors. The IIR is calculated from the quasiparticle (QP) width in the $GW$A, since it can be identified as the decay rate of a QP into lower energy QP plus an independent electron-hole pair. The quasiparticle self-consistent $GW$ method was used to generate the noninteracting hamiltonian the $GW$A requires as input. Small empirical corrections were added so as to reproduce experimental band gaps. Our results are in reasonable agreement with previous work, though we observe some discrepancy. In particular we find high IIR at low energy in the narrow gap semiconductor InAs.