Researcher profile

Mark van Schilfgaarde

Mark van Schilfgaarde contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
14works
0followers
8topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

14 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.