Hybridization-expansion Continuous-time Quantum Research Articles

Improved Hubbard-I approximation impurity solver for quantum impurity models

We develop a fast impurity solver which is based on the combination of Hubbard-I approximation and hybridization expansion continuous-time quantum Monte Carlo algorithm. This solver inherits the advantages of both algorithms. In order to demonstrate the power and usefulness of this solver, we use it to study the magnetic phase transitions of single-band and two-band Hubbard models in the framework of single-site dynamical mean-field theory. The calculated results agree well with those obtained by hybridization expansion quantum impurity solver. It is suggested that this solver is very suitable to solve the multi-orbital quantum impurity models efficiently and can be used to study more realistic systems with magnetic long-range order in the future.

Formula omitted]QIST v0.7: An open source continuous-time quantum Monte Carlo impurity solver toolkit

formula omitted]QIST v0.7: An open source continuous-time quantum Monte Carlo impurity solver toolkit

Detecting phase transitions and crossovers in Hubbard models using the fidelity susceptibility

A generalized version of the fidelity susceptibility of single-band and multi-orbital Hubbard models is systematically studied using single-site dynamical mean-field theory in combination with a hybridization expansion continuous-time quantum Monte Carlo impurity solver. We find that the fidelity susceptibility is extremely sensitive to changes in the state of the system. It can be used as a numerically inexpensive tool to detect and characterize a broad range of phase transitions and crossovers in Hubbard models, including (orbital-selective) Mott metal-insulator transitions, high-spin to low-spin transitions, Fermi-liquid to non-Fermi-liquid crossovers, and spin-freezing crossovers.

Kernel polynomial representation for imaginary-time Green’s functions in continuous-time quantum Monte Carlo impurity solver**Project supported by the National Natural Science Foundation of China (Grant No. 11504340).

Inspired by the recently proposed Legendre orthogonal polynomial representation for imaginary-time Green’s functions G(τ), we develop an alternate and superior representation for G(τ) and implement it in the hybridization expansion continuous-time quantum Monte Carlo impurity solver. This representation is based on the kernel polynomial method, which introduces some integral kernel functions to filter the numerical fluctuations caused by the explicit truncations of polynomial expansion series and can improve the computational precision significantly. As an illustration of the new representation, we re-examine the imaginary-time Green’s functions of the single-band Hubbard model in the framework of dynamical mean-field theory. The calculated results suggest that with carefully chosen integral kernel functions, whether the system is metallic or insulating, the Gibbs oscillations found in the previous Legendre orthogonal polynomial representation have been vastly suppressed and remarkable corrections to the measured Green’s functions have been obtained.

First-principles investigation of cubicBaRuO3: A Hund's metal

A first-principles investigation of cubic-BaRuO$_3$, by combining density functional theory with dynamical mean-field theory and a hybridization expansion continuous time quantum Monte-Carlo solver, has been carried out. Non-magnetic calculations with appropriately chosen on-site Coulomb repulsion, $U$ and Hund's exchange, $J$, for single-particle dynamics and static susceptibility show that cubic-BaRuO$_3$ is in a spin-frozen state at temperatures above the ferromagnetic transition point. A strong red shift with increasing $J$ of the peak in the real frequency dynamical susceptibility indicates a dramatic suppression of the Fermi liquid coherence scale as compared to the bare parameters in cubic-BaRuO$_3$. The self-energy also shows clear deviation from Fermi liquid behaviour that manifests in the single-particle spectrum. Such a clean separation of energy scales in this system provides scope for an incoherent spin-frozen (SF) phase, that extends over a wide temperature range, to manifest in non-Fermi liquid behaviour and to be the precursor for the magnetically ordered ground state.

Formula omitted]: 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 iQIST). 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 iQIST are also explained. Finally, several simple tutorials are presented in order to demonstrate the usage and power of iQIST. Program summaryProgram title:iQISTCatalogue identifier: AEWQ_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEWQ_v1_0.htmlProgram obtainable from: CPC Program Library, Queen’s University, Belfast, N. IrelandLicensing provisions: GNU General Public License, version 3No. of lines in distributed program, including test data, etc.: 226270No. of bytes in distributed program, including test data, etc.: 5263144Distribution format: tar.gzProgramming language: Fortran 2008 and Python.Computer: Desktop PC, laptop, high performance computing cluster.Operating system: Unix, Linux, Mac OS X, Windows.Has the code been vectorized or parallelized?: Yes, it is parallelized by MPI and OpenMPRAM: Depends on the complexity of the problemClassification: 7.3.External routines: BLAS, LAPACK, Latex is required to build the user manual.Nature of problem:Quantum impurity models were originally proposed to describe magnetic impurities in metallic hosts. In these models, the Coulomb interaction acts between electrons occupying the orbitals of the impurity atom. Electrons can hop between the impurity and the host, and in an action formulation, this hopping is described by a time-dependent hybridization function. Nowadays quantum impurity models have a broad range of applications, from the description of heavy fermion systems, and Kondo insulators, to quantum dots in nano-science. They also play an important role as auxiliary problems in dynamical mean-field theory and its diagrammatic extensions [1–3], where an interacting lattice model is mapped onto a quantum impurity model in a self-consistent manner. Thus, the accurate and efficient solution of quantum impurity models becomes an essential task.Solution method:The quantum impurity model can be solved by the numerically exact continuous-time quantum Monte Carlo method, which is the most efficient and powerful impurity solver for finite temperature simulations. In the iQIST software package, we implemented the hybridization expansion version of continuous-time quantum Monte Carlo algorithm. Both the segment representation and general matrix formalism are supported. The key idea of this algorithm is to expand the partition function diagrammatically in powers of the impurity-bath hybridization, and to stochastically sample these diagrams to all relevant orders using the Metropolis Monte Carlo algorithm. For a detailed review of the continuous-time quantum Monte Carlo algorithms, please refer to [4].Running time:Depends on the complexity of the problem. The sample run supplied in the distribution takes about 1.5 minutes.

Odd-parity triplet superconducting phase in multiorbital materials with a strong spin-orbit coupling: application to doped Sr₂IrO₄.

We explore possible superconducting states in t(2g) multiorbital correlated electron systems with strong spin-orbit coupling (SOC). In order to study such systems in a controlled manner, we employ large-scale dynamical mean-field theory (DMFT) simulations with the hybridization expansion continuous-time quantum Monte Carlo (CTQMC) impurity solver. To determine the pairing symmetry, we go beyond the local DMFT formalism using parquet equations to introduce the momentum dependence in the two-particle vertex and correlation functions. In the strong SOC limit, a singlet, d-wave pairing state in the electron-doped side of the phase diagram is observed at weak Hund's coupling, which is triggered by antiferromagnetic fluctuations. When the Hund's coupling is comparable to SOC, a twofold degenerate, triplet p-wave pairing state with relatively high transition temperature emerges in the hole-doped side of the phase diagram, which is associated with enhanced charge fluctuations. Experimental implications to doped Sr2IrO4 are discussed.

Improved Green’s function measurement for hybridization expansion quantum Monte Carlo

We present an algorithm for measurement of the Green’s function in the hybridization expansion continuous-time quantum Monte-Carlo method based on continuous estimators. Compared to the standard method, the present algorithm has a similar or better accuracy with an improvement notable especially at low perturbation orders and high frequencies. The resulting statistical noise is weakly correlated so high-accuracy data for the numerical analytical continuation can be produced.

Metal-insulator transition in three-band Hubbard model with strong spin-orbit interaction

Recent investigations suggest that both spin-orbit coupling and electron correlation play very crucial roles in the 5d transition metal oxides. By using the generalized Gutzwiller variational method and dynamical mean-field theory with the hybridization expansion continuous time quantum Monte Carlo as impurity solver, the three-band Hubbard model with full Hund’s rule coupling and spin-orbit interaction terms, which contains the essential physics of partially filled t 2g sub-shell of 5d materials, is studied systematically. The calculated phase diagram of this model exhibits three distinct phase regions, which are metal, band insulator and Mott insulator. We find that the spin-orbit coupling term tends to greatly enhance the tendency towards the Mott insulator phase. Furthermore, the influence of the electron-electron interaction on the effective strength of spin-orbit coupling in the metallic phase is studied in detail. We conclude that the electron correlation effect on the effective spin-orbit coupling is far beyond the mean-field treatment even in the intermediate coupling region.

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.

Magnetic properties of the two-band Hubbard model with different bandwidths

Abstract We investigate magnetic properties of the two-band Hubbard model with different bandwidths. By means of dynamical mean-field theory combined with the hybridization expansion continuous time quantum Monte Carlo method, we discuss how the difference in the bandwidths affects the magnetism at low temperatures. It is found that the nature of the magnetic phase transition is not affected by the difference in the bandwidths at half filling and thus a single phase transition occurs at the same critical temperature. By contrast, the introduction of holes induces an orbital-selective insulator–metal transition. The characteristic properties of the intermediate phase with one orbital localized and the other itinerant are also discussed.