Source author record

Karsten Held

Karsten Held 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

51works
6topics
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

51 published item(s)

preprint2022arXiv

Beyond DMFT: Spin Fluctuations, Pseudogaps and Superconductivity

This Chapter of the Jülich Summer School 2022 provides for a first reading on diagrammatic extensions of dynamical-mean-field theory and their application to spin fluctuations, pseudogap physics and supercondcutivtiy. The contents is as follows: 1 Introduction 2 Dynamical vertex approximation 3 Ladder dynamical vertex approximation 4 Hubbard model, cuprates and nickelates 5 Spin fluctuations 6 Pseudogap 7 Superconductivity 8 Conclusion and outlook

preprint2022arXiv

Interplay between magnetic and superconducting fluctuations in the doped 2d Hubbard model

We study the Hubbard model on a square lattice, using the dynamical vertex approximation and the parquet approximation. These methods allow us to describe the mutual interference of spin-fluctuations in the particle-hole channel and superconducting fluctuations in the cooperon channel in an unbiased way. For small dopings we find predominant commensurable antiferromagnetic spin- and d-wave superconducting fluctuations; for larger doping incommensurate antiferromagnetic spin fluctuations are concomitant to triplet s-wave superconducting fluctuations.

preprint2022arXiv

Long-term memory magnetic correlations in the Hubbard model: A dynamical mean-field theory analysis

We investigate the onset of a not-decaying asymptotic behavior of temporal magnetic correlations in the Hubbard model in infinite dimensions. This long-term memory feature of dynamical spin correlations can be precisely quantified by computing the difference between the zero-frequency limit of the Kubo susceptibility and the corresponding static isothermal one. Here, we present a procedure for reliably evaluating this difference starting from imaginary time-axis data, and apply it to the testbed case of the Mott-Hubbard metal-insulator transition (MIT). At low temperatures, we find long-term memory effects in the entire Mott regime, abruptly ending at the first order MIT. This directly reflects the underlying local moment physics and the associated degeneracy in the many-electron spectrum. At higher temperatures, a more gradual onset of an infinitely-long time-decay of magnetic correlations occurs in the crossover regime, not too far from the Widom line emerging from the critical point of the MIT. Our work has relevant algorithmic implications for the analytical continuation of dynamical susceptibilities in strongly correlated regimes and offers a new perspective for unveiling fundamental properties of the many-particle spectrum of the problem under scrutiny.

preprint2022arXiv

Merging numerical renormalization group and intermediate representation to compactify two- and three-point correlators

The vanguard of many-body theory is nowadays dealing with the full frequency dynamics of n-point Green's functions for n higher than two. Numerically, these objects easily become a memory bottleneck, even when working with discrete imaginary-time Matsubara frequencies. Here, we use the intermediate representation (IR) to compactify the two-point Green's function and three-point Fermion-Bose vertex directly on the real frequency axis, on the basis of numerical renormalization group (NRG) data. We empirically observe an upper bound of the relative error when comparing the IR reconstructed signal with the original NRG data, and demonstrate that a IR compacification is possible.

preprint2022arXiv

Particle-hole asymmetric lifetimes promoted by spin and orbital fluctuations in SrVO$_3$ monolayers

The two-dimensional nature of engineered transition-metal ultra-thin oxide films offers a large playground of yet to be fully understood physics. Here, we study pristine SrVO$_3$ monolayers that have recently been predicted to display a variety of magnetic and orbital orders. Above all ordering temperatures, we find that the associated non-local fluctuations lead to a momentum differentiation in the self-energy, particularly in the scattering rate. In the one-band 2D Hubbard model, momentum-selectivity on the Fermi surface ("$k=k_F$") is known to lead to pseudogap physics. Here instead, in the multi-orbital case, we evidence a differentiation between momenta on the occupied ("$k<k_F$") and the unoccupied side ("$k>k_F$") of the Fermi surface. Our work, based on the dynamical vertex approximation, complements the understanding of spectral signatures of non-local fluctuations, calls to (re)examine other ultra-thin oxide films and interfaces with methods beyond dynamical mean-field theory, and may point to correlation-enhanced thermoelectric effects.

preprint2022arXiv

Phase diagram of nickelate superconductors calculated by dynamical vertex approximation

We review the electronic structure of nickelate superconductors with and without effects of electronic correlations. As a minimal model we identify the one-band Hubbard model for the Ni 3$d_{x^2-y^2}$ orbital plus a pocket around the $A$-momentum. The latter however merely acts as a decoupled electron reservoir. This reservoir makes a careful translation from {nominal} Sr-doping to the doping of the one-band Hubbard model mandatory. Our dynamical mean-field theory calculations, in part already supported by experiment, indicate that the $Γ$ pocket, Nd 4$f$ orbitals, oxygen 2$p$ and {the} other Ni 3$d$ orbitals are not relevant in the superconducting doping regime. The physics is completely different if topotactic hydrogen is present or the oxygen reduction is incomplete. Then, a two-band physics hosted by the Ni 3$d_{x^2-y^2}$ and 3$d_{3z^2-r^2}$ orbitals emerges. Based on our minimal modeling we calculated the superconducting $T_c$ vs. Sr-doping $x$ phase diagram prior to experiment using the dynamical vertex approximation. For such a notoriously difficult to determine quantity as $T_c$, the agreement with experiment is astonishingly good. The prediction that $T_c$ is enhanced with pressure or compressive strain, has been confirmed experimentally as well. This supports that the one-band Hubbard model plus an electron reservoir is the appropriate minimal model.

preprint2022arXiv

Photoexcitations in the Hubbard model -- generalized Loschmidt amplitude analysis of impact ionization in small clusters

We study photoexcitations in small Hubbard clusters of up to 12 sites, some of which show an increase of the double occupation after the electric field pulse through impact ionization. Here, the time-dependent electromagnetic field is introduced through a Peierls substitution and the time evolution is calculated by exact diagonalization with commutator-free Magnus integrators. As a tool to better analyze the out-of-equilibrium dynamics, we generalize the Loschmidt amplitude. This way, we are able to resolve which many-body energy eigenstates are responsible for impact ionization and which show pronounced changes in the double occupation and spin energy. This analysis reveals that the loss of spin energy is of little importance for impact ionization. We further demonstrate that, for one-dimensional chains, the optical conductivity has a characteristic peak structure originating solely from vertex corrections.

preprint2022arXiv

Signatures of topotactic hydrogen in nickelate superconductors

Superconductivity has entered the nickel age marked by enormous experimental and theoretical efforts. Notwithstanding, synthesizing nickelate superconductors remains extremely challenging, not least due to incomplete oxygen reduction and topotactic hydrogen. Here, we present density-functional theory calculations, identify a phonon mode as a possible indication for topotactic hydrogen and discuss the charge redistribution patterns around oxygen and hydrogen impurities.

preprint2022arXiv

Strongly correlated superconductivity with long-range spatial fluctuations

We review recent studies for superconductivity using diagrammatic extensions of dynamical mean field theory. These approaches take into account simultaneously both, the local correlation effect and spatial long-range fluctuations, which are essential to describe unconventional superconductivity in a quasi-two-dimensional plane. The results reproduce and predict the experimental phase diagrams of strongly correlated system such as cuprates and nickelates. Further studies reveal that the dynamical screening effect of the pairing interaction vertex has dramatic consequences for the transition temperature and may even support exotic mechanisms like odd-frequency pairing. We also discuss the dimensionality of layered materials and how to interpret the numerical results in two dimensions.

preprint2022arXiv

Toward Functionalized Ultrathin Oxide Films: the Impact of Surface Apical Oxygen

Thin films of transition metal oxides open up a gateway to nanoscale electronic devices beyond silicon characterized by novel electronic functionalities. While such films are commonly prepared in an oxygen atmosphere, they are typically considered to be ideally terminated with the stoichiometric composition. Using the prototypical correlated metal SrVO$_3$ as an example, it is demonstrated that this idealized description overlooks an essential ingredient: oxygen adsorbing at the surface apical sites. The oxygen adatoms, which persist even in an ultrahigh vacuum environment, are shown to severely affect the intrinsic electronic structure of a transition metal oxide film. Their presence leads to the formation of an electronically dead surface layer but also alters the band filling and the electron correlations in the thin films. These findings highlight that it is important to take into account surface apical oxygen or -- mutatis mutandis -- the specific oxygen configuration imposed by a capping layer to predict the behavior of ultrathin films of transition metal oxides near the single unit-cell limit.

preprint2022arXiv

Zoology of spin and orbital fluctuations in ultrathin oxide films

Many metallic transition-metal oxides turn insulating when grown as films that are only a few unit-cells thick. The microscopic origins of these thickness induced metal-to-insulator transitions however remain under dispute. Here, we simulate the extreme case of a monolayer of an inconspicuous correlated metal -- the strontium vanadate SrVO$_3$ -- deposited on a SrTiO$_3$ substrate. Crucially, our system can have a termination to vacuum consisting of either a SrO or a VO$_2$ top layer. While we find that both lead to Mott insulating behavior at nominal stoichiometry, the phase diagram emerging upon doping -- chemically or through an applied gate voltage -- is qualitatively different. Indeed, our many-body calculations reveal a cornucopia of nonlocal fluctuations associated with (in)commensurate antiferromagnetic, ferromagnetic, as well as stripe and checkerboard orbital ordering instabilities. Identifying that the two geometries yield crystal-field splittings of opposite signs, we elucidate the ensuing phases through the lens of the orbital degrees of freedom. Quite generally, our work highlights that interface and surface reconstruction and the deformation or severing of coordination polyhedra in ultra-thin films drive rich properties that are radically different from the material's bulk physics.

preprint2021arXiv

Correlations turn electronic structure of finite-layer nickelates upside down

Motivated by the recent discovery of superconductivity in the pentalayer nickelate Nd$_6$Ni$_5$O$_{12}$ [Nature Materials 10.1038], we calculate its electronic structure and superconducting critical temperature. We find that electronic correlations are essential for pushing Nd$_6$Ni$_5$O$_{12}$ into the superconducting doping range as they shift the electron pockets above the Fermi energy. As a consequence, Nd$_6$Ni$_5$O$_{12}$ can be described with a single $d_{x^2-y^2}$ orbital per Ni. Instead, for the bilayer nickelate Nd$_3$Ni$_2$O$_6$ we find correlations to drive the system into a three-orbital regime also involving the Ni $d_{xz,yz}$ states. We suggest, however, that single-orbital physics with optimal doping can be restored by substituting 60% of the trivalent Nd or La by tetravalent Zr.

preprint2021arXiv

Efficient Magnus-type integrators for solar energy conversion in Hubbard models

Strongly interacting electrons in solids are generically described by Hubbardtype models, and the impact of solar light can be modeled by an additional time-dependence. This yields a finite dimensional system of ordinary differential equations (ODE)s of Schrödinger type, which can be solved numerically by exponential time integrators of Magnus type. The efficiency may be enhanced by combining these with operator splittings. We will discuss several different approaches of employing exponential-based methods in conjunction with an adaptive Lanczos method for the evaluation of matrix exponentials and compare their accuracy and efficiency. For each integrator, we use defect-based local error estimators to enable adaptive time-stepping. This serves to reliably control the approximation error and reduce the computational effort

preprint2021arXiv

Hidden, one-dimensional, strongly nested, and almost half-filled Fermi surface in Ba$_2$CuO$_{3+y}$ superconductors

All previous cuprate superconductors display a set of common features: (i) vicinity to a Cu 3$d^{9}$ configuration; (ii) separated CuO$_2$ planes; (iii) superconductivity for doping $δ\sim$ 0.1$-$0.3. Recently [PNAS {\bf 24}, 12156 (2019)] challenged this picture by discovering "highly overdoped" superconducting Ba$_2$CuO$_{3+y}$. Using density-functional theory + dynamical mean-field theory, we reveal a bilayer structure of Ba$_2$CuO$_{3.2}$ of alternating quasi 2D and quasi 1D character. Correlations tune an inter-layer self-doping leading to an almost half-filled, strongly nested quasi 1D $d_{b^2-c^2}$ band, which is prone to strong antiferromagnetic fluctuations, possibly at the origin of superconductivity in Ba$_2$CuO$_{3+y}$.

preprint2021arXiv

Tiling with triangles: parquet and $GWγ$ methods unified

The parquet formalism and Hedin's $GWγ$ approach are unified into a single theory of vertex corrections, corresponding to an exact reformulation of the parquet equations in terms of boson exchange. The method has no drawbacks compared to previous parquet solvers but has the significant advantage that the vertex functions decay quickly with frequencies and with respect to distances in real space. These properties coincide with the respective separation of the length and energy scales of the two-particle correlations into long/short-ranged and high/low-energetic.

preprint2020arXiv

Electric-field-controllable high-spin SrRuO3 driven by a solid ionic junction

Controlling magnetism and spin structures in strongly correlated systems by using electric field is of fundamental importance but challenging. Here, a high-spin ruthenate phase is achieved via a solid ionic chemical junction at SrRuO3/SrTiO3 interface with distinct formation energies and diffusion barriers of oxygen vacancies, analogue to electronic band alignment in semiconductor heterojunction. Oxygen vacancies trapped within this interfacial SrRuO3 reconstruct Ru-4d electronic structure and orbital occupancy, leading to an enhanced magnetic moment. Furthermore, an interfacial magnetic phase can be switched reversibly by electric-field-rectifying oxygen migration in a solid-state ionic gating device, providing a framework for atomic design of functionalities in strongly correlated oxides using a way of solid chemistry.

preprint2020arXiv

Generic optical excitations of correlated systems: $π$-tons

The interaction of light with solids gives rise to new bosonic quasiparticles, with the exciton being---undoubtedly---the most famous of these polaritons. While excitons are the generic polaritons of semiconductors, we show that for strongly correlated systems another polariton is prevalent---originating from the dominant antiferromagnetic or charge density wave fluctuations in these systems. As these are usually associated with a wave vector ${\mathbf k}= (π,π,\ldots)$ or close to it, we propose to call the derived polaritons $π$-tons. These $π$-tons yield the leading vertex correction to the optical conductivity in all correlated models studied: the Hubbard, the extended Hubbard model, the Falicov-Kimball, and the Pariser-Parr-Pople model, both in the insulating and in the metallic phase.

preprint2020arXiv

Nickelate superconductors -- a renaissance of the one-band Hubbard model

Following the discovery of superconductivity in the cuprates and the seminal work by Anderson, the theoretical efforts to understand high-temperature superconductivity have been focusing to a large extent on a simple model: the one-band Hubbard model. However, superconducting cuprates need to be doped, and the doped holes go into the oxygen orbitals. This requires a more elaborate multi-band model such as the three-orbital Emery model. The recently discovered nickelate superconductors appear, at first glance, to be even more complicated multi-orbital systems. Here, we analyse this multi-orbital system and find that it is instead the nickelates which can be described by a one-band Hubbard model, albeit with an additional electron reservoir and only around the superconducting regime. Our calculations of the critical temperature Tc are in good agreement with experiment, and show that optimal doping is slightly below the 20% Sr-doping of Ref. 11. Even more promising than 3d nickelates are 4d palladates.

preprint2020arXiv

Parquet dual fermion approach for the Falicov-Kimball model

In the Falicov-Kimball model, a model for (annealed) disorder, we expect weak localization corrections to the optical conductivity. However, we get such weak localization effects only when employing a pp-ladder approximation in the dual fermion approach. In the full parquet approach these pp-contributions are suppressed by ph-reducible diagrams. For the optical conductivity, we find that the pht-channel yields the main contribution, even in the region where weak localization in the pp-ladder was indicated.

preprint2020arXiv

Self-consistent ladder D$Γ$A approach

We present and implement a self-consistent D$Γ$A approach for multi-orbital models and ab initio materials calculations. It is applied to the one-band Hubbard model at various interaction strengths with and without doping, to the two-band Hubbard model with two largely different bandwidths, and to SrVO$_3$. The self-energy feedback reduces critical temperatures compared to dynamical mean-field theory, even to zero temperature in two-dimensions. Compared to a one-shot, non-self-consistent calculation the non-local correlations are significantly reduced when they are strong. In case non-local correlations are weak to moderate as for SrVO$_3$, one-shot calculations are sufficient.

preprint2020arXiv

Statistical error estimates in dynamical mean-field theory and extensions thereof

We employ the jackknife algorithm to analyze the propagation of the statistical quantum Monte Carlo error through the Bethe--Salpeter equation. This allows us to estimate the error of dynamical mean-field theory calculations of the susceptibility and of dynamical vertex approximation calculations of the self-energy. We find that the different frequency components of the susceptibility are uncorrelated, whereas those of the self-energy are correlated. For improving the quality of the correlation matrix taking sufficiently many jackknife bins is key, while for reducing the standard error of the mean sufficiently many Monte Carlo measurements are necessary. We furthermore show that even in the case of the self-energy, the finite covariance does not have a sizable influence on the analytic continuation.

preprint2020arXiv

Terahertz Excitonics in Carbon Nanotubes: Exciton Autoionization and Multiplication

Excitons play major roles in optical processes in modern semiconductors, such as single-wall carbon nanotubes (SWCNTs), transition metal dichalcogenides, and 2D perovskite quantum wells. They possess extremely large binding energies (>100~meV), dominating absorption and emission spectra even at high temperatures. The large binding energies imply that they are stable, that is, hard to ionize, rendering them seemingly unsuited for optoelectronic devices that require mobile charge carriers, especially terahertz emitters and solar cells. Here, we have conducted terahertz emission and photocurrent studies on films of aligned single-chirality semiconducting SWCNTs and find that excitons autoionize, i.e., spontaneously dissociate into electrons and holes. This process naturally occurs ultrafast (<1~ps) while conserving energy and momentum. The created carriers can then be accelerated to emit a burst of terahertz radiation when a dc bias is applied, with promising efficiency in comparison to standard GaAs-based emitters. Furthermore, at high bias, the accelerated carriers acquire high enough kinetic energy to create secondary excitons through impact exciton generation, again in a fully energy and momentum conserving fashion. This exciton multiplication process leads to a nonlinear photocurrent increase as a function of bias. Our theoretical simulations based on nonequilibrium Boltzmann transport equations, taking into account all possible scattering pathways and a realistic band structure, reproduce all our experimental data semi-quantitatively. These results not only elucidate the momentum-dependent ultrafast dynamics of excitons and carriers in SWCNTs but also suggest promising routes toward terahertz excitonics despite the orders-of-magnitude mismatch between the exciton binding energies and the terahertz photon energies.

preprint2020arXiv

Topotactic hydrogen in nickelate superconductors and akin infinite-layer oxides ABO2

Superconducting nickelates appear to be difficult to synthesize. Since the chemical reduction of ABO3 (A: rare earth; B transition metal) with CaH2 may result in both, ABO2 and ABO2H, we calculate the topotactic H binding energy by density functional theory (DFT). We find intercalating H is energetically favorable for LaNiO2 but not for Sr-doped NdNiO2. This has dramatic consequences for the electronic structure as determined by DFT+dynamical mean field theory: that of 3d9 LaNiO2 is similar to (doped) cuprates, 3d8 LaNiO2H is a two-orbital Mott insulator. Topotactical H might hence explain why some nickelates are superconducting and others are not.

preprint2020arXiv

Truncated Unity Parquet Solver

We present an implementation of a truncated unity parquet solver (TUPS) which solves the parquet equations using a truncated form-factor basis for the fermionic momenta. This way fluctuations from different scattering channels are treated on an equal footing. The essentially linear scaling of computational costs in the number of untruncated bosonic momenta allows us to treat system sizes of up to 76x76 discrete lattice momenta, unprecedented by previous unbiased methods that include the frequency dependence of the vertex. With TUPS, we provide the first numerical evidence that the parquet approximation might indeed respect the Mermin-Wagner theorem and further systematically analyze the convergence with respect to the number of form factors. Using a single form factor seems to qualitatively describe the physics of the half-filled Hubbard model correctly, including the pseudogap behaviour. Quantitatively, using a single or a few form factors only is not sufficient at lower temperatures or stronger coupling.

preprint2019arXiv

The AbinitioD$Γ$A Project v1.0: Non-local correlations beyond and susceptibilities within dynamical mean-field theory

The ab initio extension of the dynamical vertex approximation (D$Γ$A) method allows for realistic materials calculations that include non-local correlations beyond $GW$ and dynamical mean-field theory. Here, we discuss the AbinitioD$Γ$A algorithm, its implementation and usage in detail, and make the program package available to the scientific community.

preprint2016arXiv

Charge self-consistency in density functional theory + dynamical mean field theory: k-space reoccupation and orbital order

We study effects of charge self-consistency within the combination of density functional theory (DFT; Wien2k) with dynamical mean field theory (DMFT; w2dynamics) in a basis of maximally localized Wannier orbitals. Using the example of two cuprates, we demonstrate that even if there is only a single Wannier orbital with fixed filling, a noteworthy charge redistribution can occur. This effect stems from a reoccupation of the Wannier orbital in k-space when going from the single, metallic DFT band to the split, insulating Hubbard bands of DMFT. We analyze another charge self-consistency effect beyond moving charge from one site to another: the correlation-enhanced orbital polarization in a freestanding layer of SrVO3.

preprint2016arXiv

Momentum structure of the self-energy and its parametrization for the two-dimensional Hubbard model

We compute the self-energy for the half-filled Hubbard model on a square lattice using lattice quantum Monte Carlo simulations and the dynamical vertex approximation. The self-energy is strongly momentum dependent, but it can be parametrized via the non-interacting energy-momentum dispersion $\varepsilon_{\mathbf{k}}$, except for pseudogap features right at the Fermi edge. That is, it can be written as $Σ(\varepsilon_{\mathbf{k}},ω)$, with two energy-like parameters ($\varepsilon$, $ω$) instead of three ($k_x$, $k_y$ and $ω$). The self-energy has two rather broad and weakly dispersing high energy features and a sharp $ω= \varepsilon_{\mathbf{k}}$ feature at high temperatures, which turns to $ω= -\varepsilon_{\mathbf{k}}$ at low temperatures. Altogether this yields a Z- and reversed-Z-like structure, respectively, for the imaginary part of $Σ(\varepsilon_{\mathbf{k}},ω)$. We attribute the change of the low energy structure to antiferromagnetic spin fluctuations.

preprint2016arXiv

Screened moments and absence of ferromagnetism in FeAl

While the stoichiometric intermetallic compound FeAl is found to be paramagnetic in experiment, standard band-theory approaches predict the material to be ferromagnetic. We show that this discrepancy can be overcome by a better treatment of electronic correlations with density functional plus dynamical mean field theory. Our results show no ferromagnetism down to 100 K and since the susceptibility is decreasing at the lowest temperatures studied we also do not expect ferromagnetism at even lower temperatures. This behavior is found to originate from temporal quantum fluctuations that screen short-lived local magnetic moments of 1.6 $μ_B$ on Fe.

preprint2016arXiv

Worm Improved Estimators in Continuous-time Quantum Monte Carlo

We derive the improved estimators for general interactions and employ these for the continuous-time quantum Monte Carlo method. Using a worm algorithm we show how measuring higher-ordered correlators leads to an improved high-frequency behavior in irreducible quantities such as the one-particle self-energy or the irreducible two-particle vertex for non-density-density interactions. A good knowledge of the asymptotics of the two-particle vertex is essential for calculating non-local electronic correlations using diagrammatic extensions to the dynamical mean field theory as well as for calculating susceptibilities. We test our algorithm against analytic results for the multi-orbital atomic-limit and the Falicov-Kimball model.

preprint2015arXiv

A route to room temperature ferromagnetic ultrathin SrRuO$_3$ films

Experimental efforts to stabilize ferromagnetism in ultrathin films of transition metal oxides have so far failed, despite expectations based on density functional theory (DFT) and DFT+U. Here, we investigate one of the most promising materials, SrRuO$_3$, and include correlation effects beyond DFT by means of dynamical mean field theory. In agreement with experiment we find an intrinsic thickness limitation for metallic ferromagnetism in SrRuO$_3$ thin films. Indeed, we demonstrate that the realization of ultrathin ferromagnetic films is out of reach of standard thin-film techniques. Proposing charge carrier doping as a new route to manipulate thin films, we predict room temperature ferromagnetism in electron-doped SrRuO$_3$ ultra thin films.

preprint2015arXiv

Continuous-time Quantum Monte Carlo using Worm Sampling

We present a worm sampling method for calculating one- and two-particle Green's functions using continuous-time quantum Monte Carlo simulations in the hybridization expansion (CT-HYB). Instead of measuring Green's functions by removing hybridization lines from partition function configurations, as in conventional CT-HYB, the worm algorithm directly samples the Green's function. We show that worm sampling is necessary to obtain general two-particle Green's functions which are not of density-density type and that it improves the sampling efficiency when approaching the atomic limit. Such two-particle Green's functions are needed to compute off-diagonal elements of susceptibilities and occur in diagrammatic extensions of the dynamical mean field theory and efficient estimators for the single-particle self-energy.

preprint2015arXiv

Efficient DMFT impurity solver using real-time dynamics with Matrix Product States

We propose to calculate spectral functions of quantum impurity models using the Time Evolving Block Decimation (TEBD) for Matrix Product States. The resolution of the spectral function is improved by a so-called linear prediction approach. We apply the method as an impurity solver within the Dynamical Mean Field Theory (DMFT) for the single- and two-band Hubbard model on the Bethe lattice. For the single-band model we observe sharp features at the inner edges of the Hubbard bands. A finite size scaling shows that they remain present in the thermodynamic limit. We analyze the real time-dependence of the double occupation after adding a single electron and observe oscillations at the same energy as the sharp feature in the Hubbard band, indicating a long-lived coherent superposition of states that correspond to the Kondo peak and the side peaks. For a two-band Hubbard model we observe an even richer structure in the Hubbard bands, which cannot be related to a multiplet structure of the impurity, in addition to sharp excitations at the band edges of a type similar to the single-band case.

preprint2015arXiv

Efficient implementation of the parquet equations -- role of the reducible vertex function and its kernel approximation

We present an efficient implementation of the parquet formalism which respects the asymptotic structure of the vertex functions at both single- and two-particle levels in momentum- and frequency-space. We identify the two-particle reducible vertex as the core function which is essential for the construction of the other vertex functions. This observation stimulates us to consider a two-level parameter-reduction for this function to simplify the solution of the parquet equations. The resulting functions, which depend on fewer arguments, are coined "kernel functions". With the use of the "kernel functions", the open boundary of various vertex functions in the Matsubara-frequency space can be faithfully satisfied. We justify our implementation by accurately reproducing the dynamical mean-field theory results from momentum-independent parquet calculations. The high-frequency asymptotics of the single-particle self-energy and the two-particle vertex are correctly reproduced, which turns out to be essential for the self-consistent determination of the parquet solutions. The current implementation is also feasible for the dynamical vertex approximation.

preprint2015arXiv

Giant switchable Rashba effect in oxide heterostructures

One of the most fundamental phenomena and a reminder of the electron's relativistic nature is the Rashba spin splitting for broken inversion symmetry. Usually this splitting is a tiny relativistic correction, hardly discernible in experiment. Interfacing a ferroelectric, BaTiO$_3$, and a heavy 5$d$ metal with a large spin-orbit coupling, Ba(Os,Ir)O$_3$, we show that giant Rashba spin splittings are indeed possible and even fully controllable by an external electric field. Based on density functional theory and a microscopic tight binding understanding, we conclude that the electric field is amplified and stored as a ferroelectric Ti-O distortion which, through the network of oxygen octahedra, also induces a large Os-O distortion. The BaTiO$_3$/BaOsO$_3$ heterostructure is hence the ideal test station for studying the fundamentals of the Rashba effect. It allows intriguing application such as the Datta-Das transistor to operate at room temperature.

preprint2014arXiv

Chebyshev expansion for Impurity Models using Matrix Product States

We improve a recently developed expansion technique for calculating real frequency spectral functions of any one-dimensional model with short-range interactions, by postprocessing computed Chebyshev moments with linear prediction. This can be achieved at virtually no cost and, in sharp contrast to existing methods based on the dampening of the moments, improves the spectral resolution rather than lowering it. We validate the method for the exactly solvable resonating level model and the single impurity Anderson model. It is capable of resolving sharp Kondo resonances, as well as peaks within the Hubbard bands when employed as an impurity solver for dynamical mean-field theory (DMFT). Our method works at zero temperature and allows for arbitrary discretization of the bath spectrum. It achieves similar precision as the dynamical density matrix renormalization group (DDMRG), at lower cost. We also propose an alternative expansion, of 1-exp(-tau H) instead of the usual H, which opens the possibility of using established methods for the time evolution of matrix product states to calculate spectral functions directly.

preprint2014arXiv

Dynamical Vertex Approximation

Dynamical vertex approximation is a Feynman diagrammatic extension of dynamical mean field theory, including non-local correlations on all time and length scales. Starting with the Dyson and the parquet equations, the lecture notes give an elementary introduction to the dynamical vertex approximation. As a benchmark, results for an exactly solvable benzene Hubbard ring are presented. Recent highlights, the calculation of the critical exponents of the Hubbard model in 3D and that long-range antiferromagnetic correlations in 2D actually shift the (paramagnetic) Mott transition to interaction U=0, are reviewed.

preprint2014arXiv

Role of impact ionization in the thermalization of photo-excited Mott insulators

We study the influence of the pulse energy and fluence on the thermalization of photo-doped Mott insulators. If the Mott gap is smaller than the width of the Hubbard bands, the kinetic energy of individual carriers can be large enough to produce doublon-hole pairs via a process analogous to impact ionization. The thermalization dynamics, which involves an adjustment of the doublon and hole densities, thus changes as a function of the energy of the photo-doped carriers and exhibits two timescales -- a fast relaxation related to impact ionization and a slower timescale associated with higher-order scattering processes. The slow dynamics depends more strongly on the gap size and the photo-doping concentration.

preprint2013arXiv

Anisotropic two-dimensional electron gas at SrTiO3(110) protected by its native overlayer

Two dimensional electron gases (2DEGs) at a oxide heterostructures are attracting considerable attention, as these might substitute conventional semiconductors for novel electronic devices [1]. Here we present a minimal set-up for such a 2DEG -the SrTiO3(110)-(4 x 1) surface, natively terminated with one monolayer of chemically-inert titania. Oxygen vacancies induced by synchrotron radiation migrate under- neath this overlayer, this leads to a confining potential and electron doping such that a 2DEG develops. Our angular resolved photoemission spectroscopy (ARPES) and theoretical results show that confinement along (110) is strikingly different from a (001) crystal orientation. In particular the quantized subbands show a surprising "semi-heavy" band, in contrast to the analogue in the bulk, and a high electronic anisotropy. This anisotropy and even the effective mass of the (110) 2DEG is tunable by doping, offering a high flexibility to engineer the properties of this system.

preprint2013arXiv

Electronics with correlated oxides: SrVO$_3$/SrTiO$_3$ as a Mott transistor

We employ density functional theory plus dynamical mean field theory and identify the physical origin of why two layers of SrVO$_3$ on a SrTiO$_3$ substrate are insulating: the thin film geometry lifts the orbital degeneracy which in turn triggers a Mott-Hubbard transition. Two layers of SrVO$_3$ are just at the verge of a Mott-Hubbard transition and hence ideally suited for technological applications of the Mott-Hubbard transition: the heterostructure is highly sensitive to strain, electric field, and temperature. A gate voltage can also turn the insulator into a metal, so that a transistor with ideal on-off (metal-insulator) switching properties is realized.

preprint2013arXiv

Importance of d-p Coulomb interaction for high T$_C$ cuprates and other oxides

Current theoretical studies of electronic correlations in transition metal oxides typically only account for the local repulsion between d-electrons even if oxygen ligand p-states are an explicit part of the effective Hamiltonian. Interatomic interactions such as Upd between d- and (ligand) p-electrons, as well as the local interaction between p-electrons, are neglected. Often, the relative d-p orbital splitting has to be adjusted "ad hoc" on the basis of the experimental evidence. By applying the merger of local density approximation and dynamical mean field theory (LDA+DMFT) to the prototypical case of the 3-band Emery dp model for the cuprates, we demonstrate that, without any "ad hoc" adjustment of the orbital splitting, the charge transfer insulating state is stabilized by the interatomic interaction Upd. Our study hence shows how to improve realistic material calculations that explicitly include the p-orbitals.

preprint2013arXiv

Oxide Heterostructures for Efficient Solar Cells

We propose an unexplored class of absorbing materials for high-efficiency solar cells: heterostructures of transition-metal oxides. In particular, LaVO_3 grown on SrTiO_3 has a direct band gap ~1.1 eV in the optimal range as well as an internal potential gradient, which can greatly help to separate the photo-generated electron-hole pairs. Furthermore, oxide heterostructures afford the flexibility to combine LaVO_3 with other materials such as LaFeO_3 in order to achieve even higher efficiencies with band-gap graded solar cells. We use density-functional theory to demonstrate these features.

preprint2013arXiv

Quantum confinement in perovskite oxide heterostructures: tight binding instead of nearly free electron picture

Most recently, orbital-selective quantum well states of $d$ electrons have been experimentally observed in SrVO$_3$ ultrathin films [K. Yoshimatsu et. al., Science 333, 319 (2011)] and SrTiO$_3$ surfaces [A. F. Santander-Syro et. al., Nature 469, 189 (2011)]. Hitherto, one tries to explain these experiments by a nearly free electron (NFE) model, an approach widely used for delocalized electrons in semiconductor heterostructures and simple metal films. We show that a tight binding (TB) model is more suitable for describing heterostructures with more localized $d$ electrons. In this paper, we construct from first principles simple TB models for perovskite oxide heterostructures and surfaces. We show that the TB model provides a simple intuitive physical picture and yields, already with only two parameters, quantitatively much more reliable results, consistent with experiment.

preprint2012arXiv

Comparing $GW$+DMFT and LDA+DMFT for the testbed material SrVO$_3$

We have implemented the $GW$+dynamical mean field theory (DMFT) approach in the Vienna ab initio simulation package. Employing the interaction values obtained from the locally unscreened random phase approximation (RPA), we compare $GW$+DMFT and LDA+DMFT against each other and against experiment for SrVO$_3$. We observed a partial compensation of stronger electronic correlations due to the reduced $GW$ bandwidth and weaker correlations due to a larger screening of the RPA interaction, so that the obtained spectra are quite similar and well agree with experiment. Noteworthily, the $GW$+DMFT better reproduces the position of the lower Hubbard side band.

preprint2012arXiv

Conserved quantities of SU(2)-invariant interactions for correlated fermions and the advantages for quantum Monte Carlo simulations

In the context of realistic calculations for strongly-correlated materials with $d$- or $f$-electrons the efficient computation of multi-orbital models is of paramount importance. Here we introduce a set of invariants for the SU(2)-symmetric Kanamori Hamiltonian which allows to massively speed up the calculation of the fermionic trace in hybridization-expansion continuous-time quantum Monte Carlo algorithms. As an application, we show that, exploiting this set of good quantum numbers, the study of the orbital-selective Mott-transition in systems with up to seven correlated orbitals becomes feasible.

preprint2012arXiv

Effective Onsite Interaction for Dynamical Mean-Field Theory

A scheme to incorporate non-local polarizations into the dynamical mean-field theory (DMFT) and a tailor-made way to determine the effective interaction for the DMFT are systematically investigated. Applying it to the two-dimensional Hubbard model, we find that non-local polarizations induce a non-trivial filling-dependent anti-screening effect for the effective interaction. The present scheme combined with density functional theory offers an ab initio way to derive effective onsite interactions for the impurity problem in DMFT. We apply it to SrVO3 and find that the anti-screening competes with the screening caused by the off-site interaction.

preprint2012arXiv

Microscopic understanding of the orbital splitting and its tuning at oxide interfaces

By means of a Wannier projection within the framework of density functional theory, we are able to identify the modified c-axis hopping and the energy mismatch between the cation bands as the main source of the $t_{2g}$ splitting around the $Γ$ point for oxide heterostructures, excluding previously proposed mechanisms such as Jahn-Teller distortions or electric field asymmetries. Interfacing LaAlO$_3$, LaVO$_3$, SrVO$_3$ and SrNbO$_3$ with SrTiO$_3$ we show how to tune this orbital splitting, designing heterostructures with more $d_{xy}$ electrons at the interface. Such an "orbital engineering" is the key for controlling the physical properties at the interface of oxide heterostructures.

preprint2012arXiv

Theory of spin-orbit coupling at LaAlO3/SrTiO3 interfaces and SrTiO3 surfaces

The theoretical understanding of the spin-orbit coupling (SOC) effects at LaAlO$_{3}$/SrTiO$_{3}$ interfaces and SrTiO$_{3}$ surfaces is still in its infancy. We perform first-principles density-functional-theory calculations and derive from these a simple tight-binding Hamiltonian, through a Wannier function projection and group theoretical analysis. We find striking differences to the standard Rashba theory for spin-orbit coupling in semiconductor heterostructures due to multi-orbital effects: by far the biggest SOC effect is at the crossing point of the $xy$ and $yz$ (or $zx$) orbitals; and around the $Γ$ point a Rashba spin splitting with a cubic dependence on the wave vector $\vec{k}$ is possible.

preprint2011arXiv

Cluster-size dependence in cellular dynamical mean-field theory

We examine the cluster-size dependence of the cellular dynamical mean-field theory (CDMFT) applied to the two-dimensional Hubbard model. Employing the continuous-time quantum Monte Carlo method as the solver for the effective cluster model, we obtain CDMFT solutions for 4-, 8-, 12-, and 16-site clusters at a low temperature. Comparing various periodization schemes, which are used to construct the infinite-lattice quantities from the cluster results, we find that the cumulant periodization yields the fastest convergence for the hole-doped Mott insulator where the most severe size dependence is expected. We also find that the convergence is much faster around (0,0) and (pi/2,pi/2) than around (pi,0) and (pi,pi). The cumulant-periodized self-energy seems to be close to its thermodynamic limit already for a 16-site cluster in the range of parameters studied. The 4-site results remarkably agree well with the 16-site results, indicating that the previous studies based on the 4-site cluster capture the essence of the physics of doped Mott insulators.

preprint2011arXiv

Effects of electronic correlations and disorder on the thermopower of NaxCoO2

For the thermoelectric properties of NaxCoO2, we analyze the effect of local Coulomb interaction and (disordered) potential differences for Co-sites with adjacent Na-ion or vacancy. The disorder potential alone increases the resistivity and reduces the thermopower, while the Coulomb interaction alone leads only to minor changes compared to the one-particle picture of the local density approximation. Only combined, these two terms give rise to a substantial increase of the thermopower: the number of (quasi-)electrons around the Fermi level is much more suppressed than that of the (quasi-)holes. Hence, there is a particle-hole imbalance acting in the same direction as a similar imbalance for the group velocities. Together, this interplay results in a large positive thermopower. Introducing a thermoelectric spectral density, we located the energies and momenta regions most relevant for the thermopower and changes thereof.

preprint2010arXiv

Dynamical vertex approximation for nanoscopic systems

With an increasing complexity of nanoscopic systems and the modeling thereof, new theoretical tools are needed for a reliable calculation of complex systems with strong electronic correlations. To this end, we propose a new approach based on the recently introduced dynamical vertex approximation. We demonstrate its reliability already on the one-particle vertex (i.e., dynamical mean field theory) level by comparison with the exact solution. Modeling a quantum point contact with 110 atoms, we show that the contact becomes insulating already before entering the tunneling regime due to a local Mott-Hubbard transition occurring on the atoms which form the point contact.

preprint2010arXiv

Wien2wannier: From linearized augmented plane waves to maximally localized Wannier functions

We present an implementaion of interface between the full-potential linearized augmented plane wave package Wien2k and the wannier90 code for the construction of maximally localized Wannier functions. The FORTRAN code and a documentation is made available and results are discussed for SrVO$_3$, Sr$_2$IrO$_4$ (including spin-orbit coupling), LaFeAsO, and FeSb$_2$.