Dynamical Mean-field Theory Simulations Research Articles

Kondo frustration via charge fluctuations: a route to Mott localisation

We propose a minimal effective impurity model that captures the phenomenology of the Mott–Hubbard metal–insulator transition of the half-filled Hubbard model on the Bethe lattice in infinite dimensions as observed by dynamical mean field theory (DMFT). This involves extending the standard Anderson impurity model Hamiltonian to include an explicit Kondo coupling J, as well as a local on-site correlation Ub on the conduction bath site connected directly to the impurity. For the case of attractive local bath correlations (), the extended Anderson impurity model (e-SIAM) sheds new light on several aspects of the DMFT phase diagram. For example, the T = 0 metal-to-insulator quantum phase transition (QPT) is preceded by an excited state QPT (ESQPT) where the local moment eigenstates are emergent in the low-lying spectrum. Long-ranged fluctuations are observed near both the QPT and ESQPT, suggesting that they are the origin of the quantum critical scaling observed recently at high temperatures in DMFT simulations. The T = 0 gapless excitations at the quantum critical point display particle-hole interconversion processes, and exhibit power-law behaviour in self-energies and two-particle correlations. These are signatures of non-Fermi liquid behaviour that emerge from the partial breakdown of the Kondo screening.

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.

Simulation of time-dependent resonant inelastic x-ray scattering using nonequilibrium dynamical mean-field theory

We develop a framework to evaluate the time-dependent resonant inelastic X-ray scattering (RIXS) signal with the use of non-equilibrium dynamical mean field theory simulations. The approach is based on the solution of a time-dependent impurity model which explicitly incorporates the probe pulse. It avoids the need to compute four-point correlation functions, and can in principle be combined with different impurity solvers. This opens a path to study time-resolved RIXS processes in multi-orbital systems. The approach is exemplified with a study of the RIXS signal of a melting Mott antiferromagnet.

Alleviating the sign problem in quantum Monte Carlo simulations of spin-orbit-coupled multiorbital 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.

Dynamical mean field theory simulations with the adaptive sampling configuration interaction method

In the pursuit of accurate descriptions of strongly correlated quantum many-body systems, dynamical mean-field theory (DMFT) has been an invaluable tool for elucidating the spectral properties and quantum phases of both phenomenological models and ab initio descriptions of real materials. Key to the DMFT process is the self-consistent map of the original system into an Anderson impurity model, the ground state of which is computed using an impurity solver. The power of the method is thus limited by the complexity of the impurity model the solver can handle. Simulating realistic systems generally requires many correlated sites. By adapting the recently proposed adaptive sampling configuration interaction (ASCI) method as an impurity solver, we enable much more efficient zero temperature DMFT simulations. The key feature of the ASCI method is that it selects only the most relevant Hilbert space degrees of freedom to describe the ground state. This reduces the numerical complexity of the calculation, which will allow us to pursue future DMFT simulations with more correlated impurity sites than in previous works. Here we present the ASCI-DMFT method and example calculations on the one-dimensional and two-dimensional Hubbard models that exemplify its efficient convergence and timing properties. We show that the ASCI approach is several orders of magnitude faster than the current best published ground state DMFT simulations, which allows us to study the bath discretization error in simulations with small clusters, as well as to address cluster sizes beyond the current state of the art. Our approach can also be adapted for other embedding methods such as density matrix embedding theory and self-energy embedding theory.

Quantum Boltzmann equation for strongly correlated systems: Comparison to dynamical mean field theory

We investigate the potential of a quantum Boltzmann equation without momentum conservation for description of strongly correlated electron systems out of equilibrium. In a spirit similar to dynamical mean field theory (DMFT), the momentum conservation of the electron-electron scattering is neglected, which yields a time-dependent occupation function for the equilibrium spectral function, even in cases where well-defined quasiparticles do not exist. The main assumption of this method is that the spectral function remains sufficiently rigid under the non-equilibrium evolution. We compare the result of the quantum Boltzmann equation to non-equilibrium DMFT simulations for the case of photo-carrier relaxation in Mott insulators, where processes on very different timescales emerge, i.e., impact ionization, intra-Hubbard-band thermalization, and full thermalization. Since quantum Boltzmann simulations without momentum conservation are computationally cheaper than non-equilibrium DMFT, this method allows the simulation of more complicated systems or devices, and to access much longer times.

Truncating the memory time in nonequilibrium dynamical mean field theory calculations

The nonequilibrium Green's functions (NEGF) approach is a versatile theoretical tool, which allows to describe the electronic structure, spectroscopy and dynamics of strongly correlated systems. The applicability of this method is, however, limited by its considerable computational cost. Due to the treatment of the full two-time dependence of the NEGF the underlying equations of motion involve a long-lasting non-Markovian memory kernel that results in at least a $N^3_t$ scaling in the number of time points $N_t$. The system's memory time is, however, reduced in the presence of a thermalizing bath. In particular, dynamical mean-field theory (DMFT) -- one of the most successful approaches to strongly correlated lattice systems - maps extended systems to an effective impurity coupled to a bath. In this work, we systematically investigate how the memory time can be truncated in nonequilibrium DMFT simulations of the Hubbard and Hubbard-Holstein model. We show that suitable truncation schemes, which substantially reduce the computational cost, result in excellent approximations to the full time evolution. This approach enables the propagation to longer times, making fundamental processes like prethermalization and the final stages of thermalization accessible to nonequilibrium DMFT.

Evidence for spin-triplet odd-parity superconductivity close to type-II van Hove singularities

Searching for unconventional Cooper pairing states has been at the heart of superconductivity research since the discovery of BCS superconductors. In particular, spin-triplet odd-parity pairing states were recently revisited due to the possibility of tuning towards topological superconductors. In this context, it is interesting to note a recent proposal that such a spin-triplet pairing instability occurs when the band filling is near van Hove singularities (vHS) associated with momenta away from time-reversal invariant momenta named type-II vHS. However, this result was obtained within a weak coupling renormalization group with Fermi surface patch approximation. To explore superconducting instabilities beyond this weak coupling Fermi surface patch approximation, we perform systematic study on Hubbard model in a two-dimensional square lattice using three different methods: random phase approximation, large-scale dynamical mean field theory simulations with continuous time quantum Monte Carlo (CTQMC) impurity solver, and large-scale dynamical cluster simulations with CTQMC cluster solver. We find, in a wide doping range centered around the type-II van Hove filling, a two-fold degenerate, spin-triplet, odd-parity $p$-wave pairing state emerges due to repulsive interaction, when the Fermi surface is not sufficiently nested. Possible relevance of our findings to the recently discovered superconductors LaO$_{1-x}$F$_{x}$BiS$_{2}$, Ir$_{1-x}$Pt$_{x}$Te$_{2}$ and proposed doped BC$_{3}$ are also discussed.