Source author record

Philipp Werner

Philipp Werner 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

116works
11topics
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

116 published item(s)

preprint2023arXiv

Non-thermal superconductivity in photo-doped multi-orbital Hubbard systems

Superconductivity in laser-excited correlated electron systems has attracted considerable interest due to reports of light-induced superconducting-like states. Here we explore the possibility of non-thermal superconducting order in strongly interacting multi-orbital Hubbard systems, using non-equilibrium dynamical mean field theory. We find that a staggered $η$-type superconducting phase can be realized on a bipartite lattice in the high photo-doping regime, if the effective temperature of the photo-carriers is sufficiently low. The $η$ superconducting state is stabilized by Hund coupling - a positive Hund coupling favors orbital-singlet spin-triplet $η$ pairing, whereas a negative Hund coupling stabilizes spin-singlet orbital-triplet $η$ pairing.

preprint2023arXiv

Quaternary borocarbides: a testbed for DFT for superconductors

Using ab-initio density functional theory for superconductors (SCDFT), we systematically study the quaternary borocarbides $RM_2$B$_2$C. Treating the retarded (frequency-dependent) interaction $W(ω)$ within the random-phase approximation (RPA), we find good agreement with experiments for the calculated superconducting critical temperature $T_c$ in the nonmagnetic Ni- and Pd-based compounds. We argue that the problem of accurately placing the $f$-bands within DFT, and possibly the lack of an explicit magnetic pair-breaking mechanism, explain the difficulties of SCDFT in reproducing $T_c$ in members of the magnetic $R$Ni$_2$B$_2$C series ($R$ rare-earth with partially filled $4f$). While the calculated $T_c$ is overestimated, SCDFT qualitatively captures the experimentally observed trend along the rare-earth series, which indicates that the electron-phonon couplings and dynamically screened interactions have a significant effect on $T_c$.

preprint2022arXiv

$GW$+EDMFT investigation of Pr$_{1-x}$Sr$_x$NiO$_2$ under pressure

Motivated by the recent experimental observation of a large pressure effect on $T_c$ in Pr$_{1-x}$Sr$_x$NiO$_2$, we study the electronic properties of this compound as a function of pressure for $x=0$ and $0.2$ doping using self-consistent $GW$+EDMFT. Our numerical results demonstrate a non-trivial interplay between chemical doping and physical pressure, and small but systematic changes in the orbital occupations, local level energies, and interaction parameters with increasing pressure. The proper treatment of correlation effects, beyond density function theory, is shown to play an important role in revealing these trends. While the pressure dependent changes in the electronic structure of the undoped compound suggest a more single-band-like behavior in the high-pressure regime, a qualitatively different behavior is found in the doped system. We also point out that the fluctuations in the orbital occupations and spin states are not consistent with a single-band picture, and that at least a two-band model is necessary to reproduce the full result. This multi-orbital nature manifests itself most clearly in the doped compound.

preprint2022arXiv

Collective theory for an interacting solid in a single-mode cavity

We investigate the control of interacting matter through strong coupling to a single electromagnetic mode, such as the photon mode in a Fabry-Perot or split-ring cavity. For this purpose, we analyze the exact effective theory for the collective light-matter hybrid modes of a generic system of $N$ transition dipoles within an interacting solid. The approach allows to predict properties of the coupled light-matter system from the nonlinear response functions of the uncoupled matter ``outside the cavity''. The limit of large $N$ corresponds to a conventional macroscopic description based on the polarizability of matter. In this limit, the cavity does not affect the static ferroelectric response. Corrections, which are needed to understand finite size systems and to obtain the nonlinear light-matter response, can be obtained from the non-linear susceptibilities of the matter outside the cavity. The theory is benchmarked for the Dicke model, and for a quantum Ising model which serves as a minimal mean-field model for a quantum paraelectric material like SrTiO3.

preprint2022arXiv

Doping asymmetry and layer-selective metal-insulator transition in trilayer K$_{3+x}$C$_{60}$

Thin films provide a versatile platform to tune electron correlations and explore new physics in strongly correlated materials. Epitaxially grown thin films of the alkali-doped fulleride K$_{3+x}$C$_{60}$, for example, exhibit %various intriguing phenomena, including Mott transitions and superconductivity, depending on dimensionality and doping. Surprisingly, in the trilayer case, a strong electron-hole doping asymmetry has been observed in the superconducting phase, which is absent in the three-dimensional bulk limit. Using density-functional theory plus dynamical mean-field theory, we show that this doping asymmetry results from a substantial charge reshuffling from the top layer to the middle layer. While the nominal filling per fullerene is close to $n=3$, the top layer rapidly switches to an $n=2$ insulating state upon hole doping, which \textcolor{black}{implies a} doping asymmetry of the superconducting gap. The interlayer charge transfer and layer-selective metal-insulator transition result from the interplay between crystal field splittings, strong Coulomb interactions, and an effectively negative Hund coupling. This peculiar charge reshuffling is absent in the monolayer system, which is an $n=3$ Mott insulator, as expected from the nominal filling.

preprint2022arXiv

Eliashberg theory of the Jahn-Teller-Hubbard model

We study a multiorbital Hubbard model coupled to local Jahn-Teller phonons to investigate the superconducting state realized in fullerides. A weak-coupling approach is employed in combination with a local self-energy approximation. In addition to the normal and anomalous self-energies of the electrons, we consider the phonon self-energy, which allows a self-consistent treatment of the energetics. The frequency dependence of the self-energies and their characteristic coefficients, such as renormalization factors and dampings, are investigated in detail using numerical calculations. It is clarified that the anisotropic phonons play an important role in the stabilization of the superconducting state. By comparing the full results to those without phonon self-energies, we show that the superconductivity is stabilized by the softening of the phonon frequency. The effects of electronic fluctuations are also considered, which leads to the coupling to orbitons, an analog of plasmons in the electron gas. This additional contribution further stabilizes the superconducting state.

preprint2022arXiv

Entropy-cooled nonequilibrium states of the Hubbard model

We show that the recently proposed cooling-by-doping mechanism allows to efficiently prepare interesting nonequilibrium states of the Hubbard model. Using nonequilibrium dynamical mean field theory and a particle-hole symmetric setup with dipolar excitations to full and empty bands we produce cold photo-doped Mott insulating states with a sharp Drude peak in the optical conductivity, a superconducting state in the repulsive Hubbard model with an inverted population, and $η$-paired states in systems with a large density of doublons and holons. The reshuffling of entropy into full and empty bands not only provides an efficient cooling mechanism, it also allows to overcome thermalization bottlenecks and slow dynamics that have been observed in systems cooled by the coupling to boson baths.

preprint2022arXiv

Mott versus hybridization gap in the low-temperature phase of $1T$-TaS$_2$

We compute the correlated electronic structure of stacked $1T$-TaS$_2$ bilayers using the $GW$ + EDMFT method. Depending on the surface termination, the semi-infinite uncorrelated system is either band-insulating or exhibits a metallic surface state. For realistic values of the onsite and intersite interactions, a Mott gap opens in the surface state, but this gap is smaller than the gap originating from the bilayer structure. Our results are consistent with recent scanning tunneling spectroscopy measurements for different terminating layers, and with our own photoemission measurements, which indicate the coexistence of spatial regions with different gaps in the electronic spectrum.

preprint2022arXiv

Pseudoparticle vertex solver for quantum impurity models

We present a quantum impurity solver based on a pseudo-particle framework, which combines diagrammatic resummations for a three-point vertex with diagrammatic Monte Carlo sampling of a four-point vertex. This recently proposed approach [A. J. Kim et al., arXiv:2112.15549] is generalized here to fermionic impurity problems and we discuss the technical details of the implementation, including the time-stepping approach, the Monte Carlo updates, and the routines for checking the two-particle irreducibility of the four-point vertex. We also explain how the vertex information can be efficiently stored using a Dubiner basis representation. The convergence properties of the algorithm are demonstrated with applications to exactly solvable impurity models and dynamical mean field theory simulations of the single-orbital Hubbard model. It is furthermore shown that the algorithm can handle a two-orbital problem with off-diagonal hybridizations, which would cause a severe sign problem in standard hybridization-expansion Monte Carlo simulations. Since the vertex-based algorithm successfully handles sign oscillating integrals in equilibrium and samples only connected diagrams, it may be a promising approach for real-time simulations.

preprint2022arXiv

Superconductivity in black phosphorus and the role of dynamical screening

Simple cubic phosphorus exhibits superconductivity with a maximum $T_c$ of up to 12 K under pressure. The pressure dependence of $T_c$ cannot be consistently explained with a simple electron-phonon mechanism, which has stimulated investigations into the role of electronic correlations and plasmonic contributions. Here, we solve the gap equation of density functional theory for superconductors using different electron-electron and electron-phonon contributions to the kernel. We find that the phonon contribution alone yields an overestimation of $T_c$, while the addition of the static electronic contribution results in an underestimation. Taking into account the full frequency dependence of the screened interaction, the one-shot $GW$ approximation predicts $T_c$ values in good agreement with the experiments in the pressure range appropriate for the cubic phase. We also explore the use of quasi-particle bands in the calculation of the electronic and phononic kernels, and show that this modification significantly improves $T_c$ in the high-pressure region.

preprint2021arXiv

Exploring nonequilibrium phases of photo-doped Mott insulators with Generalized Gibbs ensembles

Photo-excited strongly correlated systems can exhibit intriguing non-thermal phases, but the theoretical investigation of them poses significant challenges. In this work, we introduce a generalized Gibbs ensemble type description for long-lived photo-doped states in Mott insulators. This framework enables systematic studies of photo-induced phases based on equilibrium methods, as demonstrated here for the one-dimensional extended Hubbard model. We determine the nonequilibrium phase diagram, which features $η$-pairing and charge density wave phases in a wide doping range, and reveal physical properties of these phases. We show that the peculiar kinematics of photo-doped carriers, and the interaction between them, play an essential role in the formation of the non-thermal phases, and we clarify the differences between photo-doped Mott insulators, chemically-doped Mott insulators and photo-doped semiconductors. Our results demonstrate a new path for the systematic exploration of nonequilibrium strongly correlated systems and show that photo-doped Mott insulators host different phases than conventional semiconductors.

preprint2021arXiv

Memory truncated Kadanoff-Baym equations

The Keldysh formalism for nonequilibrium Green's functions is a powerful theoretical framework for the description of the electronic structure, spectroscopy, and dynamics of strongly correlated systems. However, the underlying Kadanoff-Baym equations (KBE) for the two-time Keldysh Green's functions involve a memory kernel which results in a high computational cost for long simulation times $t_\text{max}$, with a cubic scaling of the computation time with $t_\text{max}$. Truncation of the memory kernel can reduce the computational cost to linear scaling with $t_\text{max}$, but the required memory times will depend on the model and the diagrammatic approximation to the self-energy. We explain how a truncation of the memory kernel can be incorporated into the time-propagation algorithm to solve the KBE, and investigate the systematic truncation of the memory kernel for the Hubbard model in different parameter regimes, and for different diagrammatic approximations. The truncation is easier to control within dynamical mean-field solutions, where it is applied to a momentum-independent self-energy. Here, simulation times up to two orders of magnitude longer are accessible both in the weak and strong coupling regime, allowing for a study of long-time phenomena such as the crossover between pre-thermalization and thermalization dynamics.

preprint2021arXiv

Vertex-based Diagrammatic Treatment of Light-Matter-Coupled Systems

We propose a diagrammatic Monte Carlo approach for general spin-boson models, which can be regarded as a generalization of the strong-coupling expansion for fermionic impurity models. The algorithm is based on a self-consistently computed three-point vertex and a stochastically sampled four-point vertex, and achieves convergence to the numerically exact result in a wide parameter regime. The performance of the algorithm is demonstrated with applications to a spin-boson model representing an emitter in a waveguide. As a function of the coupling strength, the spin exhibits a delocalization-localization crossover at low temperatures, signaling a qualitative change in the real-time relaxation. In certain parameter regimes, the response functions of the emitter coupled to the electromagnetic continuum can be described by an effective Rabi model with appropriately defined parameters. We also discuss the spatial distribution of the photon density around the emitter.

preprint2020arXiv

Alleviating the Sign Problem in Quantum Monte Carlo Simulations of Spin-Orbit-Coupled Multi-Orbital Hubbard Models

We present a strategy to alleviate the sign problem in continuous-time quantum Monte Carlo (CTQMC) simulations of the dynamical-mean-field-theory (DMFT) equations for the spin-orbit-coupled multiorbital Hubbard model. We first identify the combinations of rotationally invariant Hund coupling terms present in the relativistic basis which lead to a severe sign problem. Exploiting the fact that the average sign in CTQMC depends on the choice of single-particle basis, we propose a bonding-antibonding basis $V_{j3/2\mathrm{BA}}$ which shows an improved average sign compared to the widely used relativistic basis for most parameter sets investigated. We then generalize this procedure by introducing a stochastic optimization algorithm that exploits the space of single-particle bases and show that $V_{j3/2\mathrm{BA}}$ is very close to optimal within the parameter space investigated. Our findings enable more efficient DMFT simulations of materials with strong spin-orbit coupling.

preprint2020arXiv

Collective Modes in Excitonic Insulators: Effects of Electron-Phonon Coupling and Signatures in Optical Response

We consider a two-band spinless model describing an excitonic insulator (EI) on the two-dimensional square lattice with anisotropic hopping parameters and electron-phonon (el-ph) coupling, inspired by the EI candidate Ta$_2$NiSe$_5$. We systematically study the nature of the collective excitations in the ordered phase which originates from the interband Coulomb interaction and the el-ph coupling. When the ordered phase is stabilized only by the Coulomb interaction (pure EI phase), its collective response exhibits a massless phase mode in addition to the amplitude mode. We show that in the BEC regime, the signal of the amplitude mode becomes less prominent and that the anisotropy in the phase mode velocities is relaxed compared to the model bandstructure. Through coupling to the lattice, the phase mode acquires a mass and the signal of the amplitude mode becomes less prominent. Importantly, character of the softening mode at the boundary between the normal semiconductor phase and the ordered phase depends on the parameter condition. In particular, we point out that even for el-ph coupling smaller than the Coulomb interaction the mode that softens to zero at the boundary can have a phonon character. We also discuss how the collective modes can be observed in the optical conductivity. Furthermore, we study the effects of nonlocal interactions on the collective modes and show the possibility of realizing a coexistence of an in-gap mode and an above-gap mode split off from the single amplitude mode in the system with the local interaction only.

preprint2020arXiv

Effects of frustration on the nonequilibrium dynamics of photo-excited lattice systems

We theoretically investigate the effects of the lattice geometry on the nonequilibrium dynamics of photo-excited carriers in a half-filled two-dimensional Hubbard model. Using a nonequilibrium generalization of the dynamical cluster approximation, we compare the relaxation dynamics in lattices which interpolate between the triangular lattice and square lattice configuration and thus reveal the role of the geometric frustration in these strongly correlated nonequilibrium systems. In particular, we show that the cooling effect resulting from the disordering of the spin background is less effective in the triangular case because of the frustration. This manifests itself in a longer relaxation time of the photo-doped population, as measured by the time-resolved photo-emission signal, and a higher effective temperature of the photo-doped carriers in the non-thermal steady state after the intra-Hubbard-band thermalization.

preprint2020arXiv

Entropy and specific heat of the infinite-dimensional three-orbital Hubbard model

The Hund coupling in multiorbital Hubbard systems induces spin freezing and associated Hund metal behavior. Using dynamical mean field theory, we explore the effect of local moment formation, spin and charge excitations on the entropy and specific heat of the three-orbital model. In particular, we demonstrate a substantial enhancement of the entropy in the spin-frozen metal phase at low temperatures, and peaks in the specific heat associated with the activation of spin and charge fluctuations at high temperature. We also clarify how these temperature scales depend on the interaction parameters and filling.

preprint2020arXiv

High-harmonic generation in one-dimensional Mott insulator

We study high-harmonic generation (HHG) in the one-dimensional Hubbard model in order to understand its relation to elementary excitations as well as the similarities and differences to semiconductors. The simulations are based on the infinite time-evolving block decimation (iTEBD) method and exact diagonalization. We clarify that the HHG originates from the doublon-holon recombination, and the scaling of the cutoff frequency is consistent with a linear dependence on the external field. We demonstrate that the subcycle features of the HHG can be reasonably described by a phenomenological three step model for a doublon-holon pair. We argue that the HHG in the one-dimensional Mott insulator is closely related to the dispersion of the doublon-holon pair with respect to its relative momentum, which is not necessarily captured by the single-particle spectrum due to the many-body nature of the elementary excitations. For the comparison to semiconductors, we introduce effective models obtained from the Schrieffer-Wolff transformation, i.e. a strong-coupling expansion, which allows us to disentangle the different processes involved in the Hubbard model: intraband dynamics of doublons and holons, interband dipole excitations, and spin exchanges. These demonstrate the formal similarity of the Mott system to the semiconductor models in the dipole gauge, and reveal that the spin dynamics, which does not directly affect the charge dynamics, can reduce the HHG intensity. We also show that the long-range component of the intraband dipole moment has a substantial effect on the HHG intensity, while the correlated hopping terms for the doublons and holons essentially determine the shape of the HHG spectrum. A new numerical method to evaluate single-particle spectra within the iTEBD method is also introduced.

preprint2020arXiv

High-harmonic generation in spin-orbit coupled systems

We study high-harmonic generation in two-dimensional electron systems with Rashba and Dresselhaus spin-orbit coupling and derive harmonic generation selection rules with the help of group theory. Based on the bandstructures of these minimal models and explicit simulations we reveal how the spin-orbit parameters control the cutoff energy in the high-harmonic spectrum. We also show that the magnetic field and polarization dependence of this spectrum provides information on the magnitude of the Rashba and Dresselhaus spin-orbit coupling parameters. The shape of the Fermi surface can be deduced at least qualitatively and if only one type of spin-orbit coupling is present, the coupling strength can be determined.

preprint2020arXiv

Nonthermal excitonic condensation near a spin-state transition

We consider a two-orbital Hubbard model with Hund coupling and crystal-field splitting and show that in the vicinity of the high-spin/low-spin transition, crystal-field quenches can induce an excitonic condensation at initial temperatures above the highest ordering temperature in equilibrium. This condensation is the effect of an increase in the spin entropy and an associated cooling of the effective electronic temperature. We identify a dynamical phase transition and show that such quenches can result in long-lived nonthermal excitonic condensates which have no analogue in the equilibrium phase diagram. The results are interpreted by means of an effective pseudo-spin model.

preprint2020arXiv

Novel $J_{\rm{eff}}$=3/2 Metallic Phase and Unconventional Superconductivity in GaTa$_4$Se$_8$

By means of density functional theory plus dynamical mean-field theory (DFT+DMFT) calculations and resonant inelastic x-ray scattering (RIXS) experiments, we investigate the high-pressure phases of the spin-orbit-coupled $J_{\rm{eff}}=3/2$ insulator GaTa$_4$Se$_8$. Its metallic phase, derived from the Mott state by applying pressure, is found to carry $J_{\rm{eff}}=3/2$ moments. The characteristic excitation peak in the RIXS spectrum maintains its destructive quantum interference of $J_{\rm{eff}}$ at the Ta $L_2$-edge up to 10.4 GPa. Our exact diagonalization based DFT+DMFT calculations including spin-orbit coupling also reveal that the $J_{\rm{eff}}=3/2$ character can be clearly identified under high pressure. These results establish the intriguing nature of the correlated metallic magnetic phase, which represents the first confirmed example of $J_{\rm{eff}}$=3/2 moments residing in a metal. They also indicate that the pressure-induced superconductivity is likely unconventional and influenced by these $J_{\rm{eff}}=3/2$ moments. Based on a self-energy analysis, we furthermore propose the possibility of doping-induced superconductivity related to a spin-freezing crossover.

preprint2020arXiv

Signatures of bosonic excitations in high-harmonic spectra of Mott insulators

The high harmonic spectrum of the Mott insulating Hubbard model has recently been shown to exhibit plateau structures with cutoff energies determined by $n$th nearest neighbor doublon-holon recombination processes. The spectrum thus allows to extract the on-site repulsion $U$. Here, we consider generalizations of the single-band Hubbard model and discuss the signatures of bosonic excitations in high harmonic spectra. Specifically, we study an electron-plasmon model which captures the essential aspects of the dynamically screened Coulomb interaction in solids and a multi-orbital Hubbard model with Hund coupling which allows to analyze the effect of local spin excitations. For the electron-plasmon model, we show that the high harmonic spectrum can reveal information about the screened and bare onsite interaction, the boson frequency, as well as the relation between boson coupling strength and boson frequency. In the multi-orbital case, string states formed by local spin excitations result in an increase of the radiation intensity and cutoff energy associated with higher order recombination processes.

preprint2020arXiv

Superconducting optical response of photodoped Mott insulators

Ultrafast laser pulses can redistribute charges in Mott insulators on extremely short time scales, leading to the fast generation of photocarriers. It has recently been demonstrated that these photocarriers can form a novel $η$--paired condensate at low temperatures, featuring a staggered superconducting pairing field. In this conference paper, we discuss the origin of the $η$--paired hidden phase and its optical response which may be detected in a pump-probe experiment. The hidden phase may be relevant for possible light-induced superconductivity in Mott insulators.

preprint2020arXiv

Ultrafast nonequilibrium evolution of excitonic modes in semiconductors

We study the time evolution of excitonic states after photo-excitation in the one-dimensional spin-less extended Falicov-Kimball model. Several numerical methods are employed and benchmarked against each other: time-dependent mean-field simulations, the second-Born approximation (2BA) within the Kadanoff-Baym formalism, the generalized Kadanoff-Baym Ansatz (GKBA) implemented with the 2BA and the infinite time-evolving block decimation (iTEBD) method. It is found that the GKBA gives the best agreement with iTEBD and captures the relevant physics. We find that excitations to the particle-hole continuum and resonant excitations of the equilibrium exciton result in a qualitatively different dynamics. In the former case, the exciton binding energy remains positive and the frequency of the corresponding coherent oscillations is smaller than the band gap. On the other hand, resonant excitations trigger a collective mode whose frequency is larger than the band gap. We discuss the origin of these different behaviors by evaluating the nonequilibrium susceptibility using the nonthermal distribution and a random phase approximation. The peculiar mode with frequency larger than the band gap is associated with a partial population inversion with a sharp energy cutoff. We also discuss the effects of the cooling by a phonon bath. We demonstrate the real-time development of coherence in the polarization, which indicates excitonic condensation out of equilibrium.

preprint2019arXiv

Ferromagnetic spin correlations in the two-dimensional Hubbard model

We analyze the dynamical nearest-neighbor and next-nearest-neighbor spin correlations in the 4-site and 8-site dynamical cluster approximation to the two-dimensional Hubbard model. Focusing on the robustness of these correlations at long imaginary times, we reveal enhanced ferromagnetic correlations on the lattice diagonal, consistent with the emergence of composite spin-1 moments at a temperature scale that essentially coincides with the pseudo-gap temperature $T^*$. We discuss these results in the context of the spin-freezing theory of unconventional superconductivity.

preprint2019arXiv

Floquet DMFT Analysis of High Harmonic Generation in Mott Insulators

We study the high harmonic generation (HHG) in Mott insulators using Floquet dynamical mean-field theory (DMFT). We show that the main origin of the HHG in Mott insulators is the doublon-holon recombination, and that the character of the HHG spectrum differs depending on the field strength. In the weaker-field regime, the HHG spectrum shows a single plateau as in the HHG from gases, and its cut-off energy $ε_{\rm cut}$ scales linearly with the field strength $E_0$ as $ε_{\rm cut}=Δ_{\rm gap} + αE_0$, where $Δ_{\rm gap}$ is the Mott gap. On the other hand, in the stronger-field regime, multiple plateaus emerge and the $m$-th cut-off scales as $ε_{\rm cut,m}=U + m E_0$. We show that this difference originates from the different dynamics of the doublons and holons in the weak- and strong-field regimes. We also comment on the similarities and differences between HHG from Mott insulators and from semiconductors. This proceedings paper complements our recent work, Phys. Rev. Lett. 121, 057405 (2018), with additional results and analyses.

preprint2019arXiv

Local Berry curvature signatures in dichroic angle-resolved photoelectron spectroscopy

Topologically nontrivial two-dimensional materials hold great promise for next-generation optoelectronic applications. However, measuring the Hall or spin-Hall response is often a challenge and practically limited to the ground state. An experimental technique for tracing the topological character in a differential fashion would provide useful insights. In this work, we show that circular dichroism angle-resolved photoelectron spectroscopy (ARPES) provides a powerful tool which can resolve the topological and quantum-geometrical character in momentum space. In particular, we investigate how to map out the signatures of the local Berry curvature by exploiting its intimate connection to the orbital angular momentum. A spin-resolved detection of the photoelectrons allows to extend the approach to spin-Chern insulators. Our predictions are corroborated by state-of-the art \emph{ab initio} simulations employing time-dependent density functional theory, complemented with model calculations. The present proposal can be extended to address topological properties in materials out of equilibrium in a time-resolved fashion.

preprint2019arXiv

NESSi: The Non-Equilibrium Systems Simulation package

The nonequilibrium dynamics of correlated many-particle systems is of interest in connection with pump-probe experiments on molecular systems and solids, as well as theoretical investigations of transport properties and relaxation processes. Nonequilibrium Green's functions are a powerful tool to study interaction effects in quantum many-particle systems out of equilibrium, and to extract physically relevant information for the interpretation of experiments. We present the open-source software package NESSi (The Non-Equilibrium Systems Simulation package) which allows to perform many-body dynamics simulations based on Green's functions on the L-shaped Kadanoff-Baym contour. NESSi contains the library libcntr which implements tools for basic operations on these nonequilibrium Green's functions, for constructing Feynman diagrams, and for the solution of integral and integro-differential equations involving contour Green's functions. The library employs a discretization of the Kadanoff-Baym contour into time $N$ points and a high-order implementation of integration routines. The total integrated error scales up to $\mathcal{O}(N^{-7})$, which is important since the numerical effort increases at least cubically with the simulation time. A distributed-memory parallelization over reciprocal space allows large-scale simulations of lattice systems. We provide a collection of example programs ranging from dynamics in simple two-level systems to problems relevant in contemporary condensed matter physics, including Hubbard clusters and Hubbard or Holstein lattice models. The libcntr library is the basis of a follow-up software package for nonequilibrium dynamical mean-field theory calculations based on strong-coupling perturbative impurity solvers.

preprint2019arXiv

Nickelate superconductors $--$ Multiorbital nature and spin freezing

Superconductivity with a remarkably high $T_c$ has recently been found in Sr-doped NdNiO$_2$ thin films. While this system bears strong similarities to the cuprates, some differences, such as a weaker antiferromagnetic exchange coupling and possible high-spin moments on the doped Ni sites have been pointed out. Here, we investigate the effect of Hund coupling and crystal field splitting in a simple model system and argue that a multiorbital description of nickelate superconductors is warranted, especially in the strongly hole-doped regime. We then look at this system from the viewpoint of the spin-freezing theory of unconventional superconductivity, which provides a unified understanding of unconventional superconductivity in a broad range of compounds. Sr$_{0.2}$Nd$_{0.8}$NiO$_2$ falls into a parameter regime influenced by two spin-freezing crossovers, one related to the emergent multi-orbital nature in the strongly doped regime and the other related to the single-band character and square lattice geometry in the weakly doped regime.

preprint2019arXiv

Photo-enhanced excitonic correlations in a Mott insulator with nonlocal interactions

We investigate the effect of nonlocal interactions on the photo-doped Mott insulating state of the two-dimensional Hubbard model using a nonequilibrium generalization of the dynamical cluster approximation. In particular, we compare the situation where the excitonic states are lying within the continuum of doublon-holon excitations to a set-up where the excitons appear within the Mott gap. In the first case, the creation of nearest-neighbor doublon-holon pairs by excitations across the Mott gap results in enhanced excitonic correlations, but these excitons quickly decay into uncorrelated doublons and holons. In the second case, photo-excitation results in long-lived excitonic states. While in a low-temperature equilibrium state, excitonic features are usually not evident in single-particle observables such as the photoemission spectrum, we show that the photo-excited nonequilibrium system can exhibit in-gap states associated with the excitons. The comparison with exact-diagonalization results for small clusters allows us to identify the signatures of the excitons in the photo-emission spectrum.

preprint2019arXiv

Revealing Hund's multiplets in Mott insulators under strong electric fields

We investigate the strong-field dynamics of a paramagnetic two-band Mott insulator using real-time dynamical mean-field theory. A dielectric breakdown occurs due to many-body Landau-Zener tunnelling, with a threshold field determined by the gap. For a large range of fields, however, we predict that the tunnelling currents are small enough to allow the observation of field-induced localization of electrons, which becomes most strikingly evident in atomic-like local spin multiplets determined by the Hund's coupling $J$. This field-induced localization might provide a way of measuring the value of $J$ in correlated materials. It should be observable in transition metal oxides using time-resolved photo-emission spectroscopy or optical measurements in the presence of strong THz field transients.

preprint2019arXiv

Screening from $e_g$ states and antiferromagnetic correlations in $d^{(1,2,3)}$ perovskites: A $GW$+EDMFT investigation

We perform a systematic {\it ab initio} study of the electronic structure of Sr(V,Mo,Mn)O$_3$ perovskites, using the parameter-free $GW$+EDMFT method. This approach self-consistently calculates effective interaction parameters, taking into account screening effects due to nonlocal charge fluctuations. Comparing the results of a 3-band ($t_{2g}$) description to those of a 5-band ($t_{2g}$+$e_g$) model, it is shown that the $e_g$ states have little effect on the low-energy properties and the plasmonic features for the first two compounds but play a more active role in SrMnO$_3$. In the case of SrMnO$_3$ paramagnetic $GW$+EDMFT yields a metallic low-temperature solution on the verge of a Mott transition, while antiferromagnetic $GW$+EDMFT produces an insulating solution with the correct gap size. We discuss the possible implications of this result for the nature of the insulating state above the Néel temperature, and the reliability of the $GW$+EDMFT scheme.

preprint2019arXiv

Ultrafast coupled charge and spin dynamics in strongly correlated NiO

Charge excitations across an electronic band gap play an important role in opto-electronics and light harvesting. In contrast to conventional semiconductors, studies of above-band-gap photoexcitations in strongly correlated materials are still in their infancy. Here we reveal the ultrafast dynamics controlled by Hund's physics in strongly correlated photo-excited NiO. By combining time-resolved two-photon photoemission experiments with state-of-the-art numerical calculations, an ultrafast ($\lesssim$ 10\,fs) relaxation due to Hund excitations and related photo-induced in-gap states are identified. Remarkably, the weight of these in-gap states displays long-lived coherent THz oscillations up to 2\,ps at low temperature. The frequency of these oscillations corresponds to the strength of the antiferromagnetic superexchange interaction in NiO and their lifetime vanishes as the Néel temperature is approached. Numerical simulations of a two-band $t$-$J$ model reveal that the THz oscillations originate from the interplay between local many-body excitations and long-range antiferromagnetic order.

preprint2016arXiv

A spin-freezing perspective on cuprates

The high-temperature superconducting state in cuprates appears if charge carriers are doped into a Mott insulating parent compound. An unresolved puzzle is the unconventional nature of the normal state above the superconducting dome, and its connection to the superconducting instability. At weak hole-doping, a "pseudo-gap" metal state with signatures of time-reversal symmetry breaking is observed, which near optimal doping changes into a "strange metal" with non-Fermi liquid properties. Qualitatively similar phase diagrams are found in multi-orbital systems, such as pnictides, where the unconventional metal states arise from a Hund coupling induced spin-freezing. Here, we show that the relevant model for cuprates, the single-orbital Hubbard model on the square lattice, can be mapped onto an effective multi-orbital problem with strong ferromagnetic Hund coupling. The spin-freezing physics of this multi-orbital system explains the phenomenology of cuprates, including the pseudo-gap, the strange metal, and the d-wave superconducting instability. Our analysis suggests that spin-freezing is the universal mechanism which controls the properties of unconventional superconductors.

preprint2016arXiv

Bosonic self-energy functional theory

We derive the self-energy functional theory for bosonic lattice systems with broken $U(1)$ symmetry by parametrizing the bosonic Baym-Kadanoff effective action in terms of one- and two-point self-energies. The formalism goes beyond other approximate methods such as the pseudoparticle variational cluster approximation, the cluster composite boson mapping, and the Bogoliubov+U theory. It simplifies to bosonic dynamical-mean field theory when constraining to local fields, whereas when neglecting kinetic contributions of non-condensed bosons it reduces to the static mean-field approximation. To benchmark the theory we study the Bose-Hubbard model on the two- and three-dimensional cubic lattice, comparing with exact results from path integral quantum Monte Carlo. We also study the frustrated square lattice with next-nearest neighbor hopping, which is beyond the reach of Monte Carlo simulations. A reference system comprising a single bosonic state, corresponding to three variational parameters, is sufficient to quantitatively describe phase-boundaries, and thermodynamical observables, while qualitatively capturing the spectral functions, as well as the enhancement of kinetic fluctuations in the frustrated case. On the basis of these findings we propose self-energy functional theory as the omnibus framework for treating bosonic lattice models, in particular, in cases where path integral quantum Monte Carlo methods suffer from severe sign problems (e.g. in the presence of non-trivial gauge fields or frustration). Self-energy functional theory enables the construction of diagrammatically sound approximations that are quantitatively precise and controlled in the number of optimization parameters, but nevertheless remain computable by modest means.

preprint2016arXiv

Dynamical screening in correlated electron systems - from lattice models to realistic materials

Recent progress in treating the dynamically screened nature of the Coulomb interaction in strongly correlated lattice models and materials is reviewed with a focus on computational schemes based on the dynamical mean field approximation. We discuss approximate and exact methods for the solution of impurity models with retarded interactions, and explain how these models appear as auxiliary problems in various extensions of the dynamical mean field formalism. The current state of the field is illustrated with results from recent applications of these schemes to U-V Hubbard models and correlated materials.

preprint2016arXiv

Electronic orders in multi-orbital Hubbard models with lifted orbital degeneracy

We study the symmetry-broken phases in two- and three-orbital Hubbard models with lifted orbital degeneracy using dynamical mean field theory. On the technical level, we explain how symmetry relations can be exploited to measure the four-point correlation functions needed for the calculation of the lattice susceptibilities. In the half-filled two-orbital model with crystal field splitting, we find an instability of the metallic phase to spin-orbital order with neither spin nor orbital moment. This ordered phase is shown to be related to the recently discovered fluctuating-moment induced spin-triplet superconducting state in the orbitally degenerate model with shifted chemical potential. In the three-orbital case, we consider the effect of a crystal field splitting on the spin-triplet superconducting state in the model with positive Hund coupling, and the spin-singlet superconducting state in the case of negative Hund coupling. It is demonstrated that for certain crystal field splittings the higher energy orbitals instead of the lower ones are relevant for superconductivity, and that $T_c$ can be enhanced by the crystal field effect. We comment on the implications of our results for the superconductivity in strontium ruthenates, and for the recently reported light-enhanced superconducting state in alkali-doped fullerides.

preprint2016arXiv

Long-range orders, spin- and orbital-freezing in the two-band Hubbard model

We solve the orbitally degenerate two-band Hubbard model within dynamical mean field theory and map out the instabilities to various symmetry-broken phases based on an analysis of the corresponding lattice susceptibilities. Phase diagrams as a function of the Hund coupling parameter J are obtained both for the model with rotationally invariant interaction and for the model with Ising-type anisotropy. For negative J, an intra-orbital spin-singlet superconducting phase appears at low temperatures, while the normal state properties are characterized by an orbital freezing phenomenon. This is the negative-J analogue of the recently discovered fluctuating-moment induced s-wave spin-triplet superconductivity in the spin-freezing regime of multi-orbital models with J>0.

preprint2016arXiv

Nonequilibrium Dynamical Cluster Approximation study of the Falicov-Kimball Model

We use a nonequilibrium implementation of the dynamical cluster approximation (DCA) to study the effect of short-range correlations on the dynamics of the two-dimensional Falicov-Kimball model after an interaction quench. As in the case of single-site dynamical mean field theory, thermalization is absent in DCA simulations, and for quenches across the metal-insulator boundary, nearest-neighbor charge correlations in the nonthermal steady state are found to be larger than in the thermal state with identical energy. We investigate to what extent it is possible to define an effective temperature of the trapped state after a quench. Based on the ratio between the lesser and retarded Green's function we conclude that a roughly thermal distribution is reached within the energy intervals corresponding to the momentum-patch dependent subbands of the spectral function. The effectively different chemical potentials of these distributions however lead to a very hot, or even negative, effective temperature in the energy intervals between these subbands.

preprint2016arXiv

Photo-induced gap closure in an excitonic insulator

We study the dynamical phase transition out of an excitonic insulator phase after photo-excitation using a time-dependent extension of the selfconsistent GW method. We connect the evolution of the photoemission spectra to the dynamics of the excitonic order parameter and identify two dynamical phase transition points marked by a slowdown in the relaxation: one critical point is connected with the trapping in a nonthermal state with reduced exciton density and the second corresponds to the thermal phase transition. The transfer of kinetic energy from the photoexcited carriers to the exciton condensate is shown to be the main mechanism for the gap melting. We analyze the low energy dynamics of screening, which strongly depends on the presence of the excitonic gap, and argue that it is difficult to interpret the static component of the screened interaction as the effective interaction of some low energy model. Instead we propose a phenomenological measure for the effective interaction which indicates that screening has minor effects on the low energy dynamics.

preprint2016arXiv

When strong correlations become weak: Consistent merging of $GW$ and DMFT

The cubic perovskite SrVO$_3$ is generally considered to be a prototype strongly correlated metal with a characteristic three-peak structure of the $d$-electron spectral function, featuring a renormalized quasiparticle band in between pronounced Hubbard sidebands. Here we show that this interpretation, which has been supported by numerous "ab-initio" simulations, has to be reconsidered. Using a fully self-consistent $GW$+extended dynamical mean-field theory calculation we find that the screening from nonlocal Coulomb interactions substantially reduces the effective local Coulomb repulsion, and at the same time leads to strong plasmonic effects. The resulting effective local interactions are too weak to produce pronounced Hubbard bands in the local spectral function, while prominent plasmon satellites appear at energies which agree with those of the experimentally observed sidebands. Our results demonstrate the important role of nonlocal interactions and dynamical screening in determining the effective interaction strength of correlated compounds.

preprint2015arXiv

$i$QIST: An open source continuous-time quantum Monte Carlo impurity solver toolkit

Quantum impurity solvers have a broad range of applications in theoretical studies of strongly correlated electron systems. Especially, they play a key role in dynamical mean-field theory calculations of correlated lattice models and realistic materials. Therefore, the development and implementation of efficient quantum impurity solvers is an important task. In this paper, we present an open source interacting quantum impurity solver toolkit (dubbed $i$QIST). This package contains several highly optimized quantum impurity solvers which are based on the hybridization expansion continuous-time quantum Monte Carlo algorithm, as well as some essential pre- and post-processing tools. We first introduce the basic principle of continuous-time quantum Monte Carlo algorithm and then discuss the implementation details and optimization strategies. The software framework, major features, and installation procedure for $i$QIST are also explained. Finally, several simple tutorials are presented in order to demonstrate the usage and power of $i$QIST.

preprint2015arXiv

Accuracy of downfolding based on the constrained random phase approximation

We study the reliability of the constrained random phase approximation (cRPA) method for the calculation of low-energy effective Hamiltonians by considering multi-orbital lattice models with one strongly correlated "target" band and two weakly correlated "screening" bands. The full multi-orbital system and the effective model are solved within dynamical mean field theory (DMFT) in a consistent way. By comparing the quasi-particle weights for the correlated bands, we examine how accurately the effective model describes the low-energy properties of the multi-band system. We show that the violation of the Pauli principle in the cRPA method leads to overscreening effects when the inter-orbital interaction is small. This problem can be overcome by using a variant of the cRPA method which restores the Pauli principle.

preprint2015arXiv

Beyond the Hubbard bands in strongly correlated lattice bosons

We investigate features in the single-particle spectral function beyond the Hubbard bands in the strongly correlated normal phase of the Bose-Hubbard model. There are two distinct classes of additional peaks generated by the bosonic statistics. The first type is thermally activated Hubbard "sidebands", with the same physical origin as the zero-temperature Hubbard bands, but generated by excitations from thermally activated local occupation number states. The second class are two-particle fluctuation resonances driven by the lattice dynamics. In the unity filling Mott insulator, this takes the form of a localized triplon combined with a dispersing holon. Both types of resonances also manifest themselves in the structure factor and the interaction modulation spectra obtained from nonequilibrium bosonic dynamical mean-field theory calculations. Our findings explain experimental lattice modulation and Bragg spectroscopy results, and they predict a strong temperature dependence of the first sideband, thereby opening the door to precise thermometry of strongly correlated lattice bosons.

preprint2015arXiv

Double-expansion impurity solver for multiorbital models with dynamically screened U and J

We present a continuous-time Monte Carlo impurity solver for multiorbital impurity models which combines a strong-coupling hybridization expansion and a weak-coupling expansion in the Hund's coupling parameter J. This double-expansion approach allows to treat the dominant density-density interactions U within the efficient segment representation. We test the approach for a two-orbital model with static interactions, and then explain how the double-expansion allows to simulate models with frequency dependent U(ω) and J(ω). The method is used to investigate spin state transitions in a toy model for fullerides, with repulsive bare J but attractive screened J.

preprint2015arXiv

Dynamics of screening in photo-doped Mott insulators

We use a nonequilibrium implementation of extended dynamical mean field theory to study the effect of dynamical screening in photo-excited Mott insulators. The insertion of doublons and holes adds low-energy screening modes and leads to a reduction of the Mott gap. The coupling to low-energy bosonic modes further- more opens new relaxation channels and significantly speeds up the thermalization process. We also consider the effect of the energy distribution of the doped carriers on the screening.

preprint2015arXiv

Effective doublon and hole temperatures in the photo-doped dynamic Hubbard model

Hirsch's dynamic Hubbard model describes the effect of orbital expansion with occupancy by coupling the doublon operator to an auxiliary boson. We use the nonequilibrium dynamical mean field method to study the properties of doublon and hole carriers in this model in the strongly correlated regime. In particular, we discuss how photodoping leads to doublon and hole populations with different effective temperatures, and we analyze the relaxation behavior as a function of the boson coupling and boson energy. In the polaronic regime, the nontrivial energy exchange between doublons, holes and bosons can result in a negative temperature distribution for the holes.

preprint2015arXiv

First order dynamical phase transitions

Recently, dynamical phase transitions have been identified based on the non-analytic behavior of the Loschmidt echo in the thermodynamic limit [Heyl et al., Phys.~Rev.~Lett.~{\bf 110}, 135704 (2013)]. By introducing conditional probability amplitudes, we show how dynamical phase transitions can be further classified, both mathematically, and potentially in experiment. This leads to the definition of first-order dynamical phase transitions. Furthermore, we develop a generalized Keldysh formalism which allows to use nonequilibrium dynamical mean-field theory to study the Loschmidt echo and dynamical phase transitions in high-dimensional, non-integrable models. We find dynamical phase transitions of first order in the Falicov-Kimball model and in the Hubbard model.

preprint2015arXiv

Negative sign problem in continuous-time quantum Monte Carlo: optimal choice of single-particle basis for impurity problems

The negative sign problem in quantum Monte Carlo (QMC) simulations of cluster impurity problems is the major bottleneck in cluster dynamical mean field calculations. In this paper we systematically investigate the dependence of the sign problem on the single-particle basis. We explore both the hybridization-expansion and the interaction-expansion variants of continuous-time QMC for three-site and four-site impurity models with baths that are diagonal in the orbital degrees of freedom. We find that the sign problem in these models can be substantially reduced by using a non-trivial single-particle basis. Such bases can be generated by diagonalizing a subset of the intracluster hoppings.

preprint2015arXiv

Non-thermal melting of Neel order in the Hubbard model

We study the unitary time evolution of antiferromagnetic order in the Hubbard model after a quench starting from the perfect Néel state. In this setup, which is well suited for experiments with cold atoms, one can distinguish fundamentally different pathways for melting of long-range order at weak and strong interaction. In the Mott insulating regime, melting of long-range order occurs due to the ultra-fast transfer of energy from charge excitations to the spin background, while local magnetic moments and their exchange coupling persist during the process. The latter can be demonstrated by a local spin-precession experiment. At weak interaction, local moments decay along with the long-range order. The dynamics is governed by residual quasiparticles, which are reflected in oscillations of the off-diagonal components of the momentum distribution. Such oscillations provide an alternative route to study the prethermalization phenomenon and its influence on the dynamics away from the integrable (noninteracting) limit. The Hubbard model is solved within nonequilibrium dynamical mean-field theory, using the density matrix-renormalization group as an impurity solver.

preprint2015arXiv

Nonequilibrium dynamical cluster theory

We study the effect of spatially nonlocal correlations on the nonequilibrium dynamics of interacting fermions by constructing the nonequilibrium dynamical cluster theory, a cluster generalization of the nonequilibrium dynamical mean-field theory (DMFT). The formalism is applied to interaction quenches in the Hubbard model in one and two dimensions, and the results are compared with data from single-site DMFT, the time-dependent density matrix renormalization group, and lattice perturbation theory. Both in one and two dimensions the double occupancy quickly thermalizes, while the momentum distribution relaxes only on much longer time scales. For the two-dimensional square lattice we find a strongly momentum-dependent evolution of the momentum distribution around the Fermi energy, with a much faster relaxation near the momenta $(0,π)$ and $(π,0)$ than near $(π/2,π/2)$. This result is interpreted as reflecting the momentum-anisotropic quasiparticle lifetime of the marginal Fermi liquid. The method is further applied to the two-dimensional Hubbard model driven by a dc electric field, where the damping of the Bloch oscillation of the current is found to be less effective than predicted by DMFT and lattice perturbation theory.

preprint2015arXiv

Nonequilibrium dynamical mean-field theory for bosonic lattice models

We develop the nonequilibrium extension of bosonic dynamical mean field theory (BDMFT) and a Nambu real-time strong-coupling perturbative impurity solver. In contrast to Gutzwiller mean-field theory and strong coupling perturbative approaches, nonequilibrium BDMFT captures not only dynamical transitions, but also damping and thermalization effects at finite temperature. We apply the formalism to quenches in the Bose-Hubbard model, starting both from the normal and Bose-condensed phases. Depending on the parameter regime, one observes qualitatively different dynamical properties, such as rapid thermalization, trapping in metastable superfluid or normal states, as well as long-lived or strongly damped amplitude oscillations. We summarize our results in non-equilibrium "phase diagrams" which map out the different dynamical regimes.

preprint2015arXiv

On the dangers of partial diagrammatic summations: Benchmarks for the two-dimensional Hubbard model in the weak-coupling regime

We study the two-dimensional Hubbard model in the weak-coupling regime and compare the self-energy obtained from various approximate diagrammatic schemes to the result of diagrammatic Monte Carlo simulations, which sum up all weak-coupling diagrams up to a given order. While dynamical mean-field theory provides a good approximation for the local part of the self-energy, including its frequency dependence, the partial summation of bubble and/or ladder diagrams typically yields worse results than second order perturbation theory. Even widely used self-consistent schemes such as GW or the fluctuation-exchange approximation (FLEX) are found to be unreliable. Combining the dynamical mean-field self-energy with the nonlocal component of GW in GW+DMFT yields improved results for the local self-energy and nonlocal self-energies of the correct order of magnitude, but here, too, a more reliable scheme is obtained by restricting the nonlocal contribution to the second order diagram. FLEX+DMFT is found to give accurate results in the low-density regime, but even worse results than FLEX near half-filling.

preprint2015arXiv

Phase diagram of pyrochlore iridates: all-in-all-out magnetic ordering and non-Fermi liquid properties

We study the prototype 5d pyrochlore iridate Y2Ir2O7 from first principles using the local density approximation and dynamical mean-field theory (LDA+DMFT). We map out the phase diagram in the space of temperature, onsite Coulomb repulsion U, and filling. Consistent with experiments, we find that an all-in-all-out ordered insulating phase is stable for realistic values of U. The trigonal crystal field enhances the magnetic moments of the jeff=1/2 doublet, and strong inter-band correlations are induced by the Coulomb interaction, which indicates that a three-band description is important. We demonstrate a substantial band narrowing in the paramagnetic metallic phase and non-Fermi liquid behavior in the electron/hole doped system originating from long-lived quasi-spin moments induced by nearly flat bands.

preprint2015arXiv

Superconductivity from emerging magnetic moments

Multi-orbital Hubbard models are shown to exhibit a spatially isotropic spin-triplet superconducting phase, where equal-spin electrons in different local orbitals are paired. This superconducting state is stabilized in the spin-freezing crossover regime, where local moments emerge in the metal phase, and the pairing is substantially assisted by spin anisotropy. The phase diagram features a superconducting dome below a non-Fermi liquid metallic region and next to a magnetically ordered phase. We suggest that this type of fluctuating-moment induced superconductivity, which is not originating from fluctuations near a quantum critical point, may be realized in spin-triplet superconductors such as strontium ruthenates and uranium compounds.

preprint2015arXiv

Superconductivity in the two-band Hubbard model

We study the two-band Hubbard model in infinite dimensions by solving the dynamical mean-field equations with a strong coupling continuous-time quantum Monte Carlo method and show that an $s$-wave superconducting state can be stabilized in the repulsively interacting case. We discuss how this superconducting state competes with the metallic and paired Mott states. The effects of the Hund coupling and crystalline electric field are also addressed.

preprint2014arXiv

Downfolding electron-phonon Hamiltonians from ab-initio calculations: application to K$_3$Picene

We propose an electron-phonon parameterization which reliably reproduces the geometry and harmonic frequencies of a real system. With respect to standard electron-phonon models, it adds a "double-counting" correction, which takes into account the lattice deformation as the system is dressed by low-energy electron-phonon processes. We show the importance of this correction by studying potassium-doped picene (K$_3$Picene), recently claimed to be a superconductor with a $T_c$ of up to 18 K. The Hamiltonian parameters are derived from ab-initio density functional theory, and the lattice model is solved by dynamical mean-field theory. Our calculations include the effects of electron-electron interactions and local electron-phonon couplings. Even with the inclusion of a strongly coupled molecular phonon, the Hubbard repulsion prevails and the system is an insulator with a small Mott gap of $\approx$ 0.2 eV.

preprint2014arXiv

Dynamical screening in La2CuO4

We show that the dynamical screening of the Coulomb interaction among Cu-d electrons in high-Tc cuprates is very strong and that a proper treatment of this effect is essential for a consistent description of the electronic structure. In particular, we find that ab-initio calculations for undoped La2CuO4 yield an insulator only if the frequency dependence of the Coulomb interaction is taken into account. We also identify a collective excitation in the screened interaction at 9 eV which is rather localized on the copper site, and which is responsible for a satellite structure at energy -13 eV, located below the p bands.

preprint2014arXiv

Electronic excitation spectra of the five-orbital Anderson impurity model: From the atomic limit to itinerant atomic magnetism

We study the competition of Coulomb interaction and hybridization effects in the five-orbital Anderson impurity model by means of continuous time quantum Monte Carlo, exact diagonalization, and Hartree Fock calculations. The dependence of the electronic excitation spectra and thermodynamic ground state properties on hybridization strength and the form of the Coulomb interaction is systematically investigated for impurity occupation number $N \approx 6$. With increasing hybridization strength, a Kondo resonance emerges, broadens and merges with some of the upper and lower Hubbard peaks. Concomitantly, there is an increase of charge fluctuations at the impurity site. In contrast to the single orbital model, some atomic multiplet peaks and exchange split satellites persist despite strong charge fluctuations. We find that Hund's coupling leads to a state that may be characterized as an itinerant single atom magnet. As the filling is increased, the magnetic moment decreases, but the spin freezing phenomenon persists up to $N\approx 8$. When the hybridization is weak, the positions of atomic ionization peaks are rather sensitive to shifts of the impurity on-site energies. This allows to distinguish atomic ionization peaks from quasiparticle peaks or satellites in the electronic excitation spectra. On the methodological side we show that a comparison between the spectra obtained from Monte Carlo and exact diagonalization calculations is possible if the charge fluctuations are properly adjusted.

preprint2014arXiv

Extended dynamical mean-field study of the Hubbard model with long range interactions

Using extended dynamical mean-field theory and its combination with the $GW$ approximation, we compute the phase diagrams and local spectral functions of the single-band extended Hubbard model on the square and simple cubic lattices, considering long range interactions up to the third nearest neighbors. The longer range interactions shift the boundaries between the metallic, charge-ordered insulating and Mott insulating phases, and lead to characteristic changes in the screening modes and local spectral functions. Momentum-dependent self-energy contributions enhance the correlation effects and thus compete with the additional screening effect from longer range Coulomb interactions. Our results suggest that the influence of longer range intersite interactions is significant, and that these effects deserve attention in realistic studies of correlated materials.

preprint2014arXiv

Extended Mean Field study of complex $ϕ^4$-theory at finite density and temperature

We review the Extended Mean Field Theory (EMFT) approximation and apply it to complex, scalar $ϕ^4$-theory on the lattice. We study the critical properties of the Bose condensation driven by a nonzero chemical potential $μ$ at both zero and nonzero temperature and determine the $(T,μ)$ phase diagram. The results are in very good agreement with recent Monte Carlo data for all parameter values considered. EMFT can be formulated directly in the thermodynamic limit which allows us to study lattice spacings for which Monte Carlo studies are not feasible with present techniques. We find that the EMFT approximation accurately reproduces many known phenomena of the exact solution, like the "Silver Blaze" behavior at zero temperature and dimensional reduction at finite temperature.

preprint2014arXiv

Field-induced polaron formation in the Holstein-Hubbard model

We study the effect of strong DC and pulsed electric fields on a Mott insulating system with coupling to optical phonons. A DC field of the order of the gap induces a metallic state characterized by polaronic features in the gap region and a partially inverted population. In this quasi-steady state, the field-induced doublon-hole production is balanced by a phonon-enhanced doublon-hole recombination. The photo-excitation of carriers by a pulsed field leads to similar modifications of the electronic structure in the gap region, and an associated reduction of the doublon life-time. We demonstrate that the field-induced localization of electrons effectively enhances the phonon coupling, an effect which should be measureable with time-resolved photoemission spectroscopy.

preprint2014arXiv

Hybridization expansion Monte Carlo simulation of multi-orbital quantum impurity problems: matrix product formalism and improved Monte Carlo sampling

We explore two complementary modifications of the hybridization-expansion continuous-time Monte Carlo method, aiming at large multi-orbital quantum impurity problems. One idea is to compute the imaginary-time propagation using a matrix product states representation. We show that bond dimensions considerably smaller than the dimension of the Hilbert space are sufficient to obtain accurate results, and that this approach scales polynomially, rather than exponentially with the number of orbitals. Based on scaling analyses, we conclude that a matrix product state implementation will outperform the exact-diagonalization based method for quantum impurity problems with more than 12 orbitals. The second idea is an improved Monte Carlo sampling scheme which is applicable to all variants of the hybridization expansion method. We show that this so-called sliding window sampling scheme speeds up the simulation by at least an order of magnitude for a broad range of model parameters, with the largest improvements at low temperature.

preprint2014arXiv

Interaction quench in the Holstein model: Thermalization crossover from electron- to phonon-dominated relaxation

We study the relaxation of the Holstein model after a sudden switch-on of the interaction by means of the nonequilibrium dynamical mean field theory, with the self-consistent Migdal approximation as an impurity solver. We show that there exists a qualitative change in the thermalization dynamics as the interaction is varied in the weak-coupling regime. On the weaker interaction side of this crossover, the phonon oscillations are damped more rapidly than the electron thermalization timescale, as determined from the relaxation of the electron momentum distribution function. On the stronger interaction side, the relaxation of the electrons becomes faster than the phonon damping. In this regime, despite long-lived phonon oscillations, a thermalized momentum distribution is realized temporarily. The origin of the thermalization crossover found here is traced back to different behaviors of the electron and phonon self-energies as a function of the electron-phonon coupling. In addition, the importance of the phonon dynamics is demonstrated by comparing the self-consistent Migdal results with those obtained with a simpler Hatree-Fock impurity solver that neglects the phonon self-energy. The latter scheme does not properly describe the evolution and thermalization of isolated electron-phonon systems.

preprint2014arXiv

Nonequilibrium dynamical mean-field theory and its applications

The study of nonequilibrium phenomena in correlated lattice systems has developed into an active and exciting branch of condensed matter physics. This research field provides rich new insights that could not be obtained from the study of equilibrium situations, and the theoretical understanding of the physics often requires the development of new concepts and methods. On the experimental side, ultra-fast pump-probe spectroscopies enable studies of excitation and relaxation phenomena in correlated electron systems, while ultra-cold atoms in optical lattices provide a new way to control and measure the time-evolution of interacting lattice systems with a vastly different characteristic timescale compared to electron systems. A theoretical description of these phenomena is challenging because, firstly, we have to compute the quantum-mechanical time-evolution of many-body systems out of equilibrium, and secondly, deal with strong-correlation effects which can be of nonperturbative nature. In this review, we discuss the nonequilibrium extension of the dynamical mean field theory (DMFT), which treats quantum fluctuations in the time domain and works directly in the thermodynamic limit. The method reduces the complexity of the calculation via a mapping to a self-consistent impurity problem. Particular emphasis is placed on a detailed derivation of the formalism, and on a discussion of numerical techniques, which enable solutions of the effective nonequilibrium DMFT impurity problem. We summarize the insights gained into the properties of the infinite-dimensional Hubbard model under strong non-equilibrium conditions. These examples illustrate the current ability of the theoretical framework to reproduce and understand fundamental nonequilibrium phenomena, such as the dielectric breakdown of Mott insulators, photo-doping, and collapse-and-revival oscillations in quenched systems.

preprint2014arXiv

Nonequilibrium dynamical mean-field theory based on weak-coupling perturbation expansions: Application to dynamical symmetry breaking in the Hubbard model

We discuss the general formalism and validity of weak-coupling perturbation theory as an impurity solver for nonequilibrium dynamical mean-field theory. The method is implemented and tested in the Hubbard model, using expansions up to fourth order for the paramagnetic phase at half filling and third order for the antiferromagnetic and paramagnetic phase away from half filling. We explore various types of weak-coupling expansions and examine the accuracy and applicability of the methods for equilibrium and nonequilibrium problems. We find that in most cases an expansion of local self-energy diagrams including all the tadpole diagrams with respect to the Weiss Green's function (bare-diagram expansion) gives more accurate results than other schemes such as self-consistent perturbation theory using the fully interacting Green's function (bold-diagram expansion). In the paramagnetic phase at half filling, the fourth-order bare expansion improves the result of the second-order expansion in the weak-coupling regime, while both expansions suddenly fail at some intermediate interaction strength. The higher-order bare perturbation is especially advantageous in the antiferromagnetic phase near half filling. We use the third-order bare perturbation expansion within the nonequilibrium dynamical mean-field theory to study dynamical symmetry breaking from the paramagnetic to the antiferromagnetic phase induced by an interaction ramp in the Hubbard model. The results show that the order parameter, after an initial exponential growth, exhibits an amplitude oscillation around a nonthermal value followed by a slow drift toward the thermal value. The transient dynamics seems to be governed by a nonthermal critical point, associated with a nonthermal universality class, which is distinct from the conventional Ginzburg-Landau theory.

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.

preprint2014arXiv

Supersolid Phase Accompanied by a Quantum Critical Point in the Intermediate Coupling Regime of the Holstein Model

We reveal that electron-phonon systems described by the Holstein model on a bipartite lattice exhibit, away from half filling, a supersolid (SS) phase characterized by coexisting charge order (CO) and superconductivity (SC), and an accompanying quantum critical point (QCP). The SS, demonstrated by the dynamical mean-field theory with a quantum Monte Carlo impurity solver, emerges in the intermediate-coupling regime, where the peak of the Tc dome is located and the metal-insulator crossover occurs. On the other hand, in the weak- and strong-coupling regimes the CO-SC boundary is of first order with no intervening SS phases. The QCP is associated with the continuous transition from SS to SC and characterized by a reentrant behavior of the SS around it. We further show that the SS-SC transition is hallmarked by diverging charge fluctuations and a kink (peak) in the superfluid density.

preprint2014arXiv

Ultra-fast photo-carrier relaxation in Mott insulators with short-range spin correlations

We compute the time-resolved photoemission spectrum after photo-doping in a two-dimensional Mott-Hubbard insulator. We find that the relaxation rate of high-energy photo-doped electrons in the paramagnetic phase scales with the strength of the nearest-neighbor spin correlations, which implies a pronounced increase of the relaxation times with temperature and excitation density. Finite doping, in contrast, opens additional scattering channels and leads to a faster relaxation. To obtain our results we have implemented a nonequilibrium version of the dynamical cluster approximation (DCA), which, in contrast to single-site dynamical mean-field theory, captures the effect of short-range correlations.

preprint2014arXiv

Ultrafast separation of photo-doped carriers in Mott antiferromagnets

We use inhomogeneous nonequilibrium dynamical mean-field theory to investigate the spreading of photo-excited carriers in Mott insulating heterostructures with strong internal fields. Antiferromagnetic correlations are found to affect the carrier dynamics in a crucial manner: An antiferromagnetic spin background can absorb energy from photo-excited carriers on an ultrafast timescale, thus enabling fast transport between different layers and the separation of electron and hole-like carriers, whereas in the paramagnetic state, carriers become localized in strong fields. This interplay between charge and spin degrees of freedom can be exploited to control the functionality of devices based on Mott insulating heterostructures with polar layers, e.g., for photovoltaic applications.

preprint2013arXiv

Dynamical Mean Field Approximation Applied to Quantum Field Theory

We apply the Dynamical Mean Field (DMFT) approximation to the real, scalar phi^4 quantum field theory. By comparing to lattice Monte Carlo calculations, perturbation theory and standard mean field theory, we test the quality of the approximation in two, three, four and five dimensions. The quantities considered in these tests are the critical coupling for the transition to the ordered phase and the associated critical exponents nu and beta. We also map out the phase diagram in four dimensions. In two and three dimensions, DMFT incorrectly predicts a first order phase transition for all bare quartic couplings, which is problematic, because the second order nature of the phase transition of lattice phi^4-theory is crucial for taking the continuum limit. Nevertheless, by extrapolating the behaviour away from the phase transition, one can obtain critical couplings and critical exponents. They differ from those of mean field theory and are much closer to the correct values. In four dimensions the transition is second order for small quartic couplings and turns weakly first order as the coupling increases beyond a tricritical value. In dimensions five and higher, DMFT gives qualitatively correct results, predicts reasonable values for the critical exponents and considerably more accurate critical couplings than standard mean field theory. The approximation works best for small values of the quartic coupling. We investigate the change from first to second order transition in the local limit of DMFT which is computationally much less intensive. We also discuss technical issues related to the convergence of the non-linear self-consistency equation solver and the solution of the effective single-site model using Fourier-space Monte Carlo updates in the presence of a phi^4-interaction.

preprint2013arXiv

Dynamical mean-field analysis of ordered phases in the half-filled Holstein-Hubbard model

We investigate ordered phases (s-wave superconducting (SC), antiferromagnetic (AF) and charge ordered (CO) phases) of the half-filled Holstein-Hubbard model using the dynamical mean-field theory in combination with a continuous time quantum Monte Carlo impurity solver. We particularly focus on the parameter regime where the electron-electron Coulomb interaction and the phonon mediated retarded attractive interaction are both comparable to the bare electron bandwidth to understand the interplay of these interactions and the retardation effects due to finite phonon frequency ω_0. Results for the transition temperatures for SC and AF as a function of the effective interaction U_eff(=U-λwith U the Coulomb repulsion and λthe phonon-mediated attraction) show deviations from the the corresponding results in the anti-adiabatic limit. These are analyzed by comparing with an effective static model in the polaron picture. We also determine the phase diagram around U_eff=0 at low, but nonzero temperatures to show that paramagnetic phases appear between AF and CO when U and λare small or large, while AF and CO directly compete with each other at moderate U and λ.

preprint2013arXiv

Efficient implementation of the continuous-time hybridization expansion quantum impurity solver

Strongly correlated quantum impurity problems appear in a wide variety of contexts ranging from nanoscience and surface physics to material science and the theory of strongly correlated lattice models, where they appear as auxiliary systems within dynamical mean-field theory. Accurate and unbiased solutions must usually be obtained numerically, and continuous-time quantum Monte Carlo algorithms, a family of algorithms based on the stochastic sampling of partition function expansions, perform well for such systems. With the present paper we provide an efficient and generic implementation of the hybridization expansion quantum impurity solver, based on the segment representation. We provide a complete implementation featuring most of the recently developed extensions and optimizations. Our implementation allows one to treat retarded interactions and provides generalized measurement routines based on improved estimators for the self-energy and for vertex functions. The solver is embedded in the ALPS-DMFT application package.

preprint2013arXiv

Gauge theory of topological phases of matter

We study the response of quantum many-body systems to coupling some of their degrees of freedom to external gauge fields. This serves to understand the current Green functions and transport properties of interacting many-body systems. Our analysis leads to a "gauge theory of states of matter" complementary to the well known Landau theory of order parameters. We illustrate the power of our approach by deriving and interpreting the gauge-invariant effective actions of (topological) superconductors, 2D electron gases exhibiting the quantized Hall- and spin-Hall effect, 3D topological insulators, as well as axion electrodynamics. We also use the theory to elucidate the structure of surface modes in these systems.

preprint2013arXiv

Long-range Coulomb interactions in surface systems: a first principles description within self-consistently combined GW and dynamical mean field theory

Systems of adatoms on semiconductor surfaces display competing ground states and exotic spectral properties typical of two-dimensional correlated electron materials which are dominated by a complex interplay of spin and charge degrees of freedom. We report a fully ab initio derivation of low energy Hamiltonians for the adatom systems Si(111):X, with X=Sn, Si, C, Pb, that we solve within self-consistent combined GW and dynamical mean field theory ("GW+DMFT"). Calculated photoemission spectra are in agreement with available experimental data. We rationalize experimentally observed tendencies from Mott physics towards charge-ordering along the series as resulting from substantial long-range interactions.

preprint2013arXiv

Nonequilibrium dynamical mean field simulation of inhomogeneous systems

We extend the nonequilibrium dynamical mean field (DMFT) formalism to inhomogeneous systems by adapting the "real-space" DMFT method to Keldysh Green's functions. Solving the coupled impurity problems using strong-coupling perturbation theory, we apply the formalism to homogeneous and inhomogeneous layered systems with strong local interactions and up to 39 layers. We study the diffusion of doublons and holes created by photo-excitation in a Mott insulating system, the time-dependent build-up of the polarization and the current induced by a linear voltage bias across a multi-layer structure, and the photo-induced current in a Mott insulator under bias.

preprint2013arXiv

Nonthermal antiferromagnetic order and nonequilibrium criticality in the Hubbard model

We study dynamical phase transitions from antiferromagnetic to paramagnetic states driven by an interaction quench in the fermionic Hubbard model using the nonequilibrium dynamical mean-field theory. We identify two dynamical transition points where the relaxation behavior qualitatively changes: one corresponds to the thermal phase transition at which the order parameter decays critically slowly in a power law $\propto t^{-1/2}$, and the other is connected to the existence of nonthermal antiferromagnetic order in systems with effective temperature above the thermal critical temperature. The frequency of the amplitude mode extrapolates to zero as one approaches the nonthermal (quasi)critical point, and thermalization is significantly delayed by the trapping in the nonthermal state. A slow relaxation of the nonthermal order is followed by a faster thermalization process.

preprint2013arXiv

Ordered phases in the Holstein-Hubbard model: Interplay of strong Coulomb interaction and electron-phonon coupling

We study the Holstein-Hubbard model at half filling to explore ordered phases including su- perconductivity (SC), antiferromagnetism (AF), and charge order (CO) in situations where the electron-electron and electron-phonon interactions are strong (~ electronic bandwidth). The model is solved in the dynamical mean-field approximation using a continuous-time quantum Monte Carlo impurity solver. We determine the superconducting transition temperature (Tc) and the SC order parameter and show that the phonon-induced retardation or the strong Coulomb interaction leads to a significant reduction and shift of the Tc dome if one interprets the system as having an effective static interaction Ueff given by the Hubbard U reduced by the phonon-mediated attraction in the static limit. This behavior is analyzed by comparison to an effective static model in the polaron representation with a reduced bandwidth. We also determine the finite-temperature phase diagram including AF and CO. In the moderate-coupling regime, there is a hysteretic region of AF and CO around Ueff = 0, while the two phases are separated by a paramagnetic metal in the weak-coupling regime and a paramagnetic insulator in the strong coupling regime.

preprint2013arXiv

Phonon-enhanced relaxation and excitation in the Holstein-Hubbard model

We study quenches of the interaction and electron-phonon coupling parameter in the Hubbard-Holstein model, using nonequilibrium dynamical mean field theory. The calculations are based on a generalized Lang-Firsov scheme for time-dependent interactions or externally driven phonons, and an approximate strong-coupling impurity solver. The interaction quench calculations reveal the phonon-assisted decay of excess doublons, while the quenches of the electron-phonon coupling lead to persistent oscillations of the phonons and to a phonon-enhanced doublon production.

preprint2012arXiv

Dielectric breakdown of Mott insulators -- doublon production and doublon heating

Using dynamical mean-field theory and the non-crossing approximation as impurity solver, we study the response of a Mott insulator to strong dc electric fields. The breakdown of the Mott insulating state is triggered by field-induced creation of doublon-hole pairs. In a previous investigation, Ref. [1], it was found that the system approaches a long-lived quasi-steady state in which the current is time-independent although the number of carriers constantly increases. Here we investigate and clarify the nature of this state, which exists only because thermalization is slow in the Hubbard model at strong coupling. The current is time-independent because doublons and holes have an infinite temperature distribution. Evidence for this fact is obtained from spectral functions and by comparing the electric current with the field-induced doublon-hole creation rate. Implications to real experiments, in systems with energy dissipation, are discussed.

preprint2012arXiv

Extracting spectral properties from Keldysh Green functions

We investigate the possibility to assist the numerically ill-posed calculation of spectral properties of interacting quantum systems in thermal equilibrium by extending the imaginary-time simulation to a finite Schwinger-Keldysh contour. The effect of this extension is tested within the standard Maximum Entropy approach to analytic continuation. We find that the inclusion of real-time data improves the resolution of structures at high energy, while the imaginary-time data are needed to correctly reproduce low-frequency features such as quasi-particle peaks. As a nonequilibrium application, we consider the calculation of time-dependent spectral functions from retarded Green function data on a finite time-interval, and compare the Maximum Entropy approach to direct Fourier transformation and a method based on Padé approximants.

preprint2012arXiv

From the Cooper problem to canted supersolids in Bose-Fermi mixtures

We calculate the phase diagram of the Bose-Fermi Hubbard model on the 3d cubic lattice at fermionic half filling and bosonic unit filling by means of single-site dynamical mean-field theory. For fast bosons, this is equivalent to the Cooper problem in which the bosons can induce s-wave pairing between the fermions. We also find miscible superfluid and canted supersolid phases depending on the interspecies coupling strength. In contrast, slow bosons favor fermionic charge density wave structures for attractive fermionic interactions. These competing instabilities lead to a rich phase diagram within reach of cold gas experiments.

preprint2012arXiv

Improved Estimators for the Self-Energy and Vertex Function in Hybridization Expansion Continuous-Time Quantum Monte Carlo Simulations

We propose efficient measurement procedures for the self-energy and vertex function of the Anderson impurity model within the hybridization expansion continuous-time quantum Monte Carlo algorithm. The method is based on the measurement of higher-order correlation functions related to the quantities being sought through the equation of motion, a technique previously introduced in the NRG context. For the case of interactions of density-density type, the additional correlators can be obtained at essentially no additional computational cost. In combination with a recently introduced method for filtering the Monte Carlo noise using a representation in terms of orthogonal polynomials, we obtain data with unprecedented accuracy. This leads to an enhanced stability in analytical continuations of the self-energy or in two-particle based theories such as the dual fermion approach. As an illustration of the method we reexamine the previously reported spin-freezing and high-spin to low-spin transitions in a two-orbital model with density-density interactions. In both cases, the vertex function undergoes significant changes, which suggests significant corrections to the dynamical mean-field solutions in dual fermion calculations.

preprint2012arXiv

Nonthermal symmetry broken states in the strongly interacting Hubbard model

We study the time evolution of the antiferromagnetic order parameter after interaction quenches in the Hubbard model. Using the nonequilibrium dynamical mean field formalism, we show that the system, after a quench from intermediate to strong interaction, is trapped in a nonthermal state which is reminiscent of a photo-doped state and protected by the slow decay of doublons. If the effective doping of this state is low enough, it exhibits robust antiferromagnetic order, even if the system is highly excited and the thermal state thus expected to be paramagnetic. We comment on the implication of our findings for the stability of nonthermal superconducting states.

preprint2012arXiv

Photo-induced states in a Mott insulator

We investigate the properties of the metallic state obtained by photo-doping carriers into a Mott insulator. In a strongly interacting system, these carriers have a long life-time, so that they can dissipate their kinetic energy to a phonon bath. In the relaxed state, the scattering rate saturates at a non-zero temperature-independent value, and the momentum-resolved spectral function features broad bands which differ from the well-defined quasi-particle bands of a chemically doped system. Our results indicate that a photo-doped Mott insulator behaves as a bad metal, in which strong scattering between doublons and holes inhibits Fermi-liquid behavior down to low temperature.

preprint2012arXiv

Relaxation dynamics of the Kondo lattice model

We study the relaxation properties of the Kondo lattice model using the nonequilibrium dynamical mean field formalism in combination with the non-crossing approximation. The system is driven out of equilibrium either by a magnetic field pulse which perturbs the local singlets, or by a sudden quench of the Kondo coupling. For relaxation processes close to thermal equilibrium (after a weak perturbation), the relaxation time increases substantially as one crosses from the local moment regime into the heavy Fermi liquid. A strong perturbation, which injects a large amount of energy, can rapidly transform the heavy Fermi liquid into a local moment state. Upon cooling, the heavy Fermi liquid reappears in a two-stage relaxation, where the first step opens the Kondo gap and the second step corresponds to a slow approach of the equilibrium state via a nonthermal pathway.

preprint2012arXiv

Repulsion-to-attraction transition in correlated electron systems triggered by a monocycle pulse

We study the time evolution of the Hubbard model driven by a half-cycle or monocycle pulsed electric field F(t) using the nonequilibrium dynamical mean-field theory. We find that for properly chosen pulse shapes the electron-electron interaction can be effectively and permanently switched from repulsive to attractive if there is no energy dissipation. The physics behind the interaction conversion is a nonadiabatic shift $δ$ of the population in momentum space. When $δ\simπ$, the shifted population relaxes to a negative-temperature state, which leads to the interaction switching. Due to electron correlation effects $δ$ deviates from the dynamical phase $ϕ=\int dt F(t)$, which enables the seemingly counterintuitive repulsion-to-attraction transition by a monocycle pulse with $ϕ=0$.

preprint2012arXiv

Rubidium superoxide: a p-electron Mott insulator

Rubidium superoxide, RbO_2, is a rare example of a solid with partially-filled electronic p states, which allows to study the interplay of spin and orbital order and other effects of strong electronic correlations in a material that is quite different from the conventional d or f electron systems. Here we show, using a combination of density functional theory (DFT) and dynamical mean-field theory, that at room temperature RbO_2 is indeed a paramagnetic Mott insulator. We construct the metal-insulator phase diagram as a function of temperature and Hubbard interaction parameters U and J. Due to the strong particle-hole asymmetry of the RbO_2 band-structure, we find strong differences compared to a simple semi-elliptical density of states, which is often used to study the multiband Hubbard model. In agreement with our previous DFT study, we also find indications for complex spin and orbital order at low temperatures.

preprint2012arXiv

Screening and Non-local Correlations in the Extended Hubbard Model from Self-Consistent Combined GW and Dynamical Mean Field Theory

We describe a recent implementation of the combined GW and dynamical mean field (DMFT) method "GW+DMFT" for the two-dimensional Hubbard model with on-site and nearest-neighbor repulsion. We clarify the relation of the GW+DMFT scheme to alternative approaches in the literature, and discuss the corresponding approximations to the free energy functional of the model. Furthermore, we describe a numerically exact technique for the self-consistent solution of the GW+DMFT equations, namely the hybridization expansion continuous-time algorithm for solving the dynamical impurity model that arises within the GW+DMFT scheme. We compute the low-temperature phase diagram of the extended Hubbard model, addressing the metal-insulator transition at small intersite interactions and the transition to a charge-ordered state for stronger intersite repulsions. Within the GW+DMFT framework, as in extended DMFT, the intersite repulsion translates into a frequency-dependence of the local effective interaction. We analyze this dependence by extracting a charateristic plasma frequency and show how it affects the local spectral function.

preprint2012arXiv

Spectral Properties of Correlated Materials: Local Vertex and Non-Local Two-Particle Correlations from Combined GW and Dynamical Mean Field Theory

We present a fully self-consistent combined GW and dynamical mean field (GW+DMFT) study of the spectral properties of the extended two-dimensional Hubbard model. The inclusion of the local dynamical vertex stemming from the DMFT self-energy and polarization is shown to cure the problems of self-consistent GW in the description of spectral properties. We calculate the momentum-resolved spectral functions, the two-particle polarization and electron loss spectra, and show that the inclusion of GW in extended DMFT leads to a narrowing of the quasi-particle width and more pronounced Hubbard bands in the metallic regime as one approaches the charge-ordering transition. Finally, the momentum-dependence introduced by GW into the extended DMFT description of collective modes is found to affect their shape, giving rise to dispersive plasmon-like long-wavelength and stripe modes.

preprint2011arXiv

Charge and spin criticality for the continuous Mott transition in a two-dimensional organic conductor

We study the continuous bandwidth-controlled Mott transition in the two-dimensional single-band Hubbard model with a focus on the critical scaling behavior of charge and spin degrees of freedom. Using plaquette cluster dynamical mean-field theory, we find charge and spin criticality consistent with experimental results for organic conductors. In particular, the charge degree of freedom measured via the local density of states at the Fermi level shows a smoother transition than expected for the Ising universality class and in single-site dynamical mean-field theory, revealing the importance of short-ranged nonlocal correlations in two spatial dimensions. The spin criticality measured via the local spin susceptibility agrees quantitatively with nuclear magnetic resonance measurements of the spin-lattice relaxation rate.

preprint2011arXiv

Continuous-time Monte Carlo methods for quantum impurity models

Quantum impurity models describe an atom or molecule embedded in a host material with which it can exchange electrons. They are basic to nanoscience as representations of quantum dots and molecular conductors and play an increasingly important role in the theory of "correlated electron" materials as auxiliary problems whose solution gives the "dynamical mean field" approximation to the self energy and local correlation functions. These applications require a method of solution which provides access to both high and low energy scales and is effective for wide classes of physically realistic models. The continuous-time quantum Monte Carlo algorithms reviewed in this article meet this challenge. We present derivations and descriptions of the algorithms in enough detail to allow other workers to write their own implementations, discuss the strengths and weaknesses of the methods, summarize the problems to which the new methods have been successfully applied and outline prospects for future applications.

preprint2011arXiv

Damping of Bloch oscillations in the Hubbard model

Using nonequilibrium dynamical mean-field theory, we study the isolated Hubbard model in a static electric field in the limit of weak interactions. Linear response behavior is established at long times, but only if the interaction exceeds a critical value, below which the system exhibits an AC-type response with Bloch oscillations. The transition from AC to DC response is defined in terms of the universal long-time behavior of the system, which does not depend on the initial condition.

preprint2011arXiv

Dynamical band flipping in fermionic lattice systems: An ac-field-driven change of the interaction from repulsive to attractive

We show theoretically that the sudden application of an appropriate ac field to correlated lattice fermions flips the band structure and effectively switches the interaction from repulsive to attractive. The nonadiabatically driven system is characterized by a negative temperature with a population inversion. We numerically demonstrate the converted interaction in an ac-driven Hubbard model with the nonequilibrium dynamical mean-field theory solved by the continuous-time quantum Monte Carlo method. Based on this, we propose an efficient ramp-up protocol for ac fields that can suppress heating, which leads to an effectively attractive Hubbard model with a temperature below the superconducting transition temperature of the equilibrium system.

preprint2011arXiv

Dynamical mean-field theory for bosons

We discuss the recently developed bosonic dynamical mean-field (B-DMFT) framework, which maps a bosonic lattice model onto the selfconsistent solution of a bosonic impurity model with coupling to a reservoir of normal and condensed bosons. The effective impurity action is derived in several ways: (i) as an approximation to the kinetic energy functional of the lattice problem, (ii) using a cavity approach, and (iii) by using an effective medium approach based on adding a one-loop correction to the selfconsistently defined condensate. To solve the impurity problem, we use a continuous-time Monte Carlo algorithm based on a sampling of a perturbation expansion in the hybridization functions and the condensate wave function. As applications of the formalism we present finite temperature B-DMFT phase diagrams for the bosonic Hubbard model on a 3d cubic and 2d square lattice, the condensate order parameter as a function of chemical potential, critical exponents for the condensate, the approach to the weakly interacting Bose gas regime for weak repulsions, and the kinetic energy as a function of temperature.

preprint2011arXiv

Low temperature properties of the infinite-dimensional attractive Hubbard model

We investigate the attractive Hubbard model in infinite spatial dimensions by combining dynamical mean-field theory with a strong-coupling continuous-time quantum Monte Carlo method. By calculating the superfluid order parameter and the density of states, we discuss the stability of the superfluid state. In the intermediate coupling region above the critical temperature, the density of states exhibits a heavy fermion behavior with a quasi-particle peak in the dense system, while a dip structure appears in the dilute system. The formation of the superfluid gap is also addressed.

preprint2011arXiv

Multi-orbital Kondo physics of Co in Cu hosts

We investigate the electronic structure of cobalt atoms on a copper surface and in a copper host by combining density functional calculations with a numerically exact continuous-time quantum Monte Carlo treatment of the five-orbital impurity problem. In both cases we find low energy resonances in the density of states of all five Co $d$-orbitals. The corresponding self-energies indicate the formation of a Fermi liquid state at low temperatures. Our calculations yield the characteristic energy scale -- the Kondo temperature -- for both systems in good agreement with experiments. We quantify the charge fluctuations in both geometries and suggest that Co in Cu must be described by an Anderson impurity model rather than by a model assuming frozen impurity valency at low energies. We show that fluctuations of the orbital degrees of freedom are crucial for explaining the Kondo temperatures obtained in our calculations and measured in experiments.

preprint2011arXiv

Optimal ramp shapes for the fermionic Hubbard model in infinite dimensions

We use non-equilibrium dynamical mean field theory and a real-time diagrammatic impurity solver to study the heating associated with time-dependent changes of the interaction in a fermionic Hubbard model. Optimal ramp shapes U(t) which minimize the excitation energy are determined for an infinitesimal change. For ramp times of a few inverse hoppings, these optimal U(t) are strongly oscillating with a frequency determined by the bandwidth. We show that the scaled versions of the optimized ramps yield substantially lower temperatures than linear ramps even far outside the perturbative regime.

preprint2011arXiv

Satellites and large doping- and temperature-dependence of electronic properties in hole-doped BaFe2As2

Over the last years, superconductivity has been discovered in several families of iron-based compounds. Despite intense research, even basic electronic properties of these materials, such as Fermi surfaces, effective electron masses, or orbital characters are still subject to debate. Here, we address an issue that has not been considered before, namely the consequences of dynamical screening of the Coulomb interactions among Fe-d electrons. We demonstrate its importance not only for correlation satellites seen in photoemission spectroscopy, but also for the low-energy electronic structure. From our analysis of the normal phase of BaFe2As2 emerges the picture of a strongly correlated compound with strongly doping- and temperature-dependent properties. In the hole overdoped regime, an incoherent metal is found, while Fermi-liquid behavior is recovered in the undoped compound. At optimal doping, the self-energy exhibits an unusual square-root energy dependence which leads to strong band renormalizations near the Fermi level.

preprint2011arXiv

Superconductivity and Pairing Fluctuations in the Half-Filled Two-Dimensional Hubbard Model

The two-dimensional Hubbard model exhibits superconductivity with d-wave symmetry even at half-filling in the presence of next-nearest neighbor hopping. Using plaquette cluster dynamical mean-field theory with a continuous-time quantum Monte Carlo impurity solver, we reveal the non-Fermi liquid character of the metallic phase in proximity to the superconducting state. Specifically, the low-frequency scattering rate for momenta near (π, 0) varies non-monotonously at low temperatures, and the dc conductivity is T-linear at elevated temperatures with an upturn upon cooling. Evidence is provided that pairing fluctuations dominate the normal-conducting state even considerably above the superconducting transition temperature.

preprint2011arXiv

Thermalization of a pump-excited Mott insulator

We use nonequilibrium dynamical mean-field theory in combination with a recently implemented strong-coupling impurity solver to investigate the relaxation of a Mott insulator after a laser excitation with frequency comparable to the Hubbard gap. The time evolution of the double occupancy exhibits a crossover from a strongly damped transient at short times towards an exponential thermalization at long times. In the limit of strong interactions, the thermalization time is consistent with the exponentially small decay rate for artificially created doublons, which was measured in ultracold atomic gases. When the interaction is comparable to the bandwidth, on the other hand, the double occupancy thermalizes within a few times the inverse bandwidth along a rapid thermalization path in which the exponential tail is absent. Similar behavior can be observed in time-resolved photoemission spectroscopy. Our results show that a simple quasi-equilibrium description of the electronic state breaks down for pump-excited Mott insulators characterized by strong interactions.

preprint2010arXiv

Continuous-Time Quantum Monte Carlo and Maximum Entropy Approach to an Imaginary-Time Formulation of Strongly Correlated Steady-State Transport

Recently Han and Heary proposed an approach to steady-state quantum transport through mesoscopic structures, which maps the non-equilibrium problem onto a family of auxiliary quantum impurity systems subject to imaginary voltages. We employ continuous-time quantum Monte-Carlo solvers to calculate accurate imaginary time data for the auxiliary models. The spectral function is obtained from a maximum entropy analytical continuation in both Matsubara frequency and complexified voltage. To enable the analytical continuation we construct a kernel which is compatible with the analytical structure of the theory. While it remains a formidable task to extract reliable spectral functions from this unbiased procedure, particularly for large voltages, our results indicate that the method in principle yields results in agreement with those obtained by other methods.

preprint2010arXiv

Dielectric breakdown of Mott insulators in dynamical mean-field theory

Using nonequilibrium dynamical mean-field theory, we compute the time evolution of the current in a Mott insulator after a strong electric field is turned on. We observe the formation of a quasistationary state in which the current is almost time-independent although the system is constantly excited. At moderately strong fields this state is stable for quite long times. The stationary current exhibits a threshold behavior as a function of the field, in which the threshold increases with the Coulomb interaction and vanishes as the metal-insulator transition is approached.

preprint2010arXiv

Dynamical mean field solution of the Bose-Hubbard model

We present the effective action and self-consistency equations for the bosonic dynamical mean field (B-DMFT) approximation to the bosonic Hubbard model and show that it provides remarkably accurate phase diagrams and correlation functions. To solve the bosonic dynamical mean field equations we use a continuous-time Monte Carlo method for bosonic impurity models based on a diagrammatic expansion in the hybridization and condensate coupling. This method is readily generalized to bosonic mixtures, spinful bosons, and Bose-Fermi mixtures.

preprint2010arXiv

Dynamical Screening in Correlated Electron Materials

We present an efficient method for incorporating the dynamical effects of the screening of the Hubbard U by electronic degrees of freedom in the solid into the single site dynamical mean field approximation. The formalism is illustrated by model system calculations which capture the essential features of the frequency dependent interactions proposed for Gd, Ni, SrVO_3 and other compounds. Screening leads to shifts in the metal-insulator phase boundary, changes in the spectral function near the Mott-Hubbard gap edge and to a renormalization of the quasiparticle weight. Hubbard bands are generically neither separated by the screened nor the unscreened interaction energy, implying that the common practice of extracting the Hubbard U from the energies of features in photoemission and inverse photoemission spectra requires reexamination.

preprint2010arXiv

Nonequilibrium dynamical mean-field calculations based on the non-crossing approximation and its generalizations

We solve the impurity problem which arises within nonequilibrium dynamical mean-field theory for the Hubbard model by means of a self-consistent perturbation expansion around the atomic limit. While the lowest order, known as the non-crossing approximation (NCA), is reliable only when the interaction U is much larger than the bandwidth, low-order corrections to the NCA turn out to be sufficient to reproduce numerically exact Monte Carlo results in a wide parameter range that covers the insulating phase and the metal-insulator crossover regime at not too low temperatures. As an application of the perturbative strong-coupling impurity solver we investigate the response of the double occupancy in the Mott insulating phase of the Hubbard model to a dynamical change of the interaction or the hopping, a technique which has been used as a probe of the Mott insulating state in ultracold fermionic gases.

preprint2010arXiv

Polarized Superfluidity in the imbalanced attractive Hubbard model

We investigate the attractive Hubbard model in infinite spatial dimensions by means of dynamical mean-field theory. Using a continuous-time Monte Carlo algorithm in the Nambu formalism as an impurity solver, we directly deal with the superfluid phase in the population imbalanced system. By calculating the superfluid order parameter, the magnetization, and the density of states, we discuss how the polarized superfluid state is realized in the attractive Hubbard model at quarter filling. We find that a drastic change in the density of states is induced by spin imbalanced populations in the superfluid state.

preprint2010arXiv

Superfluid state in the periodic Anderson model with attractive interactions

We investigate the periodic Anderson model with attractive interactions by means of dynamical mean-field theory (DMFT). Using a continuous-time quantum Monte Carlo impurity solver, we study the competition between the superfluid state and the paramagnetic Kondo insulating state, and determine the phase diagram. At the chemical potential-induced phase transition from the Kondo insulating state to the superfluid state, a low-energy peak characteristic of the superfluid state appears inside the hybridization gap. We also address the effect of the confining potential in optical lattice systems by means of real-space DMFT calculations.

preprint2009arXiv

Interaction quench in the Hubbard model: Relaxation of the spectral function and the optical conductivity

We use non-equilibrium dynamical mean-field theory in combination with a recently developed Quantum Monte Carlo impurity solver to study the real-time dynamics of a Hubbard model which is driven out of equilibrium by a sudden increase in the on-site repulsion U. We discuss the implementation of the self-consistency procedure and some important technical improvements of the QMC method. The exact numerical solution is compared to iterated perturbation theory, which is found to produce accurate results only for weak interaction or short times. Furthermore we calculate the spectral functions and the optical conductivity from a Fourier transform on the finite Keldysh contour, for which the numerically accessible timescales allow to resolve the formation of Hubbard bands and a gap in the strongly interacting regime. The spectral function, and all one-particle quantities that can be calculated from it, thermalize rapidly at the transition between qualitatively different weak- and strong-coupling relaxation regimes.

preprint2009arXiv

Magnetism and orbital-ordering in an interacting three band model: a dynamical mean field study

Single-site dynamical mean field theory is used to determine the magnetic and orbital-ordering phase diagram for a model of electrons moving on a lattice with three orbital states per site and with the fully rotationally invariant Slater-Kanamori on-site interactions. The model captures important aspects of the physics of transition metal oxides with partially filled $t_{2g}$ shells, and of electron-doped C$_{60}$. We introduce an unbiased, computationally simple and inexpensive method for estimating the presence of two sublattice order, determine the regimes in which spatially uniform and two-sublattice spin and orbital orderings are present and give physical arguments for the origins of the different phases. Guidelines are determined for optimizing the presence of ferromagnetism, which may be desirable in applications

preprint2009arXiv

Weak-coupling quantum Monte Carlo calculations on the Keldysh contour: theory and application to the current-voltage characteristics of the Anderson model

We present optimized implementations of the weak-coupling continuous-time Monte Carlo method defined for nonequilibrium problems on the Keldysh contour. We describe and compare two methods of preparing the system before beginning the real-time calculation: the "interaction quench" and the "voltage quench", which are found to be suitable for large and small voltage biasses, respectively. We also discuss technical optimizations which increase the efficiency of the real-time measurements. The methods allow the accurate simulation of transport through quantum dots over wider interaction ranges and longer times than have heretofore been possible. The current-voltage characteristics of the particle-hole symmetric Anderson impurity model is presented for interactions U up to 10 times the intrinsic level width Gamma. We compare the Monte Carlo results to fourth order perturbation theory, finding that perturbation theory begins to fail at U/Gamma>4. Within the parameter range studied we find no evidence for a splitting of the Kondo resonance due to the applied voltage. The interplay of voltage and temperature and the Coulomb blockade conductance regime are studied.

preprint2007arXiv

Doping-driven Mott transition in the one-band Hubbard model

A powerful new impurity solver is shown to permit a systematic study of the doping driven Mott transition in a one-band Hubbard model within the framework of single-site dynamical mean field theory. At small dopings and large interaction strengths we are able access low enough temperatures that a reliable extrapolation to temperature T=0 may be performed, and ground state energies of insulating and metallic states may be compared. We find that the T=0 doping-driven transition is of second order and is characterized by an interaction-strength dependent electronic compressibility, which vanishes at the critical interaction strength of the half filled model. Over wide parameter ranges the compressibility is substantially reduced relative to the non-interacting system. The metal insulator transition is characterized by the appearance of in-gap states, but these are relevant only for very low dopings of less than 3%.

preprint2007arXiv

Performance analysis of continuous-time solvers for quantum impurity models

Impurity solvers play an essential role in the numerical investigation of strongly correlated electrons systems within the "dynamical mean field" approximation. Recently, a new class of continuous-time solvers has been developed, based on a diagrammatic expansion of the partition function in either the interactions or the impurity-bath hybridization. We investigate the performance of these two complementary approaches and compare them to the well-established Hirsch-Fye method. The results show that the continuous-time methods, and in particular the version which expands in the hybridization, provide substantial gains in computational efficiency.

preprint2007arXiv

Resistively shunted Josephson junctions: QFT predictions versus MC results

During the last fourteen years several exact results were obtained for the so-called boundary sine-Gordon model. In the case of a conformal bulk this 2D boundary quantum field theory describes the universal scaling behavior of the Caldeira-Leggett model of resistively shunted Josephson junctions. In this work, we use a recently developed Monte Carlo technique to test some of the analytical predictions.

preprint2004arXiv

Quantum spin chains with site dissipation

We use Monte Carlo simulations to study chains of Ising- and XY-spins with dissipation coupling to the site variables. The phase diagram and critical exponents of the dissipative Ising chain in a transverse magnetic field have been computed previously, and here we consider a universal ratio of susceptibilities. We furthermore present the phase diagram and exponents of the dissipative XY-chain, which exhibits a second order phase transition. All our results compare well with the predictions from a dissipative $ϕ^4$ field theory.