No Arabic abstract
We present an efficient ab initio dynamical mean-field theory (DMFT) implementation for quantitative simulations in solids. Our DMFT scheme employs ab initio Hamiltonians defined for impurities comprising the full unit cell or a supercell of atoms and for realistic quantum chemical basis sets. We avoid double counting errors by using Hartree-Fock as the low-level theory. Intrinsic and projected atomic orbitals (IAO+PAO) are chosen as the local embedding basis, facilitating numerical bath truncation. Using an efficient integral transformation and coupled-cluster Greens function (CCGF) impurity solvers, we are able to handle embedded impurity problems with several hundred orbitals. We apply our ab initio DMFT approach to study a hexagonal boron nitride monolayer, crystalline silicon, and nickel oxide in the antiferromagnetic phase, with up to 104 and 78 impurity orbitals in spin-restricted and unrestricted cluster DMFT calculations and over 100 bath orbitals. We show that our scheme produces accurate spectral functions compared to both benchmark periodic coupled-cluster computations and experimental spectra.
We describe the use of coupled-cluster theory as an impurity solver in dynamical mean-field theory (DMFT) and its cluster extensions. We present numerical results at the level of coupled-cluster theory with single and double excitations (CCSD) for the density of states and self-energies of cluster impurity problems in the one- and two-dimensional Hubbard models. Comparison to exact diagonalization shows that CCSD produces accurate density of states and self-energies at a variety of values of $U/t$ and filling fractions. However, the low cost allows for the use of many bath sites, which we define by a discretization of the hybridization directly on the real frequency axis. We observe convergence of dynamical quantities using approximately 30 bath sites per impurity site, with our largest 4-site cluster DMFT calculation using 120 bath sites. We suggest coupled cluster impurity solvers will be attractive in ab initio formulations of dynamical mean-field theory.
We present here two alternative schemes designed to correct the high-frequency truncation errors in the numerical treatment of the Bethe-Salpeter equations. The schemes are applicable to all Bethe-Salpeter calculations with a local two-particle irreducible local, which is relevant, e.g., for the dynamical mean-field theory (DMFT) and its diagrammatic extensions. In particular, within a purely diagrammatic framework, we could extend existing algorithms for treating the static case in the particle-hole sector to more general procedures applicable to all bosonic frequencies and all channels. After illustrating the derivation and the theoretical interrelation of the two proposed schemes, these have been applied to the Bethe-Salpeter equations for the auxiliary Anderson impurity models of selected DMFT calculations, where results can be compared against a numerically exact solution. The successful performance of the proposed schemes suggests that their implementation can significantly improve the accuracy of dynamical mean-field theory (DMFT) calculations at the two-particle level, in particular for more realistic multi-orbital calculations where the large number of degrees of freedom substantially restricts the actual frequency range for numerical calculations, as well as -on a broader perspective- of the diagrammatic extensions of DMFT.
We present an efficient and numerically stable algorithm for calculation of two-particle response functions within the dynamical mean-field theory. The technique is based on inferring the high frequency asymptotic behavior of the irreducible vertex function from the local dynamical susceptibility. The algorithm is tested on several examples. In all cases rapid convergence of the vertex function towards its asymptotic form is observed.
We study layered systems and heterostructures of s-wave superconductors by means of a suitable generalization of Dynamical Mean-Field Theory. In order to reduce the computational effort, we consider an embedding scheme in which a relatively small number of active layers is embedded in an effective potential accounting for the effect of the rest of the system. We introduce a feedback of the active layers on the embedding potential that improves on previous approaches and essentially eliminates the effects of the finiteness of the active slab allowing for cheap computation of very large systems. We extend the method to the superconducting state, and we benchmark the approach by means of simple paradigmatic examples showing some examples on how an interface affects the superconducting properties. As examples, we show that superconductivity can penetrate from an intermediate coupling superconductor into a weaker coupling one for around ten layers, and that the first two layers of a system with repulsive interaction can turn superconducting by proximity effects even when charge redistribution is inhibited.
The dynamical mean-field theory (DMFT) is a widely applicable approximation scheme for the investigation of correlated quantum many-particle systems on a lattice, e.g., electrons in solids and cold atoms in optical lattices. In particular, the combination of the DMFT with conventional methods for the calculation of electronic band structures has led to a powerful numerical approach which allows one to explore the properties of correlated materials. In this introductory article we discuss the foundations of the DMFT, derive the underlying self-consistency equations, and present several applications which have provided important insights into the properties of correlated matter.