Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userKrylov implementation of the hybridization expansion impurity solver and application to 5-orbital models
242
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
We present an implementation of the hybridization expansion impurity solver which employs sparse matrix exact-diagonalization techniques to compute the time evolution of the local Hamiltonian. This method avoids computationally expensive matrix-matrix multiplications and becomes advantageous over the conventional implementation for models with 5 or more orbitals. In particular, this method will allow the systematic investigation of 7-orbital systems (lanthanide and actinide compounds) within single-site dynamical mean field theory. We illustrate the power and usefulness of our approach with dynamical mean field results for a 5-orbital model which captures some aspects of the physics of the iron based superconductors.
rate research
Read More
Orbital and spin ordering phenomena in strongly correlated systems are commonly studied using the local-density approximation + dynamical mean-field theory approach. Typically, however, such simulations are restricted to simplified models (density-density Coulomb interactions, high symmetry couplings and few-band models). In this work we implement an efficient general hybridization-expansion continuous-time quantum Monte Carlo impurity solver (Krylov approach) which allows us to investigate orbital and spin ordering in a more realistic setting, including interactions that are often neglected (e.g., spin-flip and pair-hopping terms), enlarged basis sets (full d versus eg), low-symmetry distortions, and reaching the very low-temperature (experimental) regime. We use this solver to study ordering phenomena in a selection of exemplary low-symmetry transition-metal oxides: LaMnO3 and rare-earth manganites as well as the perovskites CaVO3 and YTiO3. We show that spin-flip and pair hopping terms do not affect the Kugel-Khomskii orbital-order melting transition in rare-earth manganites, or the suppression of orbital fluctuations driven by crystal field and Coulomb repulsion. For the Mott insulator YTiO3 we find a ferromagnetic transition temperature 50 K, in remarkably good agreement with experiments. For LaMnO3 we show that the classical t2g-spin approximation, commonly adopted for studying manganites, yields indeed an occupied eg orbital in very good agreement with that obtained for the full d 5-orbital Hubbard model, while the spin-spin e_g-t_{2g} correlation function calculated from the full d model is 0.74, very close to the value expected for aligned eg and t2g spins; the eg spectral function matrix is also well reproduced. Finally, we show that the t2g screening reduces the eg-eg Coulomb repulsion by about 10%
We present a very efficient solver for the general Anderson impurity problem. It is based on the perturbation around a solution obtained from exact diagonalization using a small number of bath sites. We formulate a perturbation theory which is valid for both weak and strong coupling and interpolates between these limits. Good agreement with numerically exact quantum Monte-Carlo results is found for a single bath site over a wide range of parameters. In particular, the Kondo resonance in the intermediate coupling regime is well reproduced for a single bath site and the lowest order correction. The method is particularly suited for low temperatures and alleviates analytical continuation of imaginary time data due to the absence of statistical noise compared to quantum Monte-Carlo impurity solvers.
We extend a previously proposed rotation and truncation scheme to optimize quantum Anderson impurity calculations with exact diagonalization [PRB 90, 085102 (2014)] to density-matrix renormalization group (DMRG) calculations. The method reduces the solution of a full impurity problem with virtually unlimited bath sites to that of a small subsystem based on a natural impurity orbital basis set. The later is solved by DMRG in combination with a restricted-active-space truncation scheme. The method allows one to compute Greens functions directly on the real frequency or time axis. We critically test the convergence of the truncation scheme using a one-band Hubbard model solved in the dynamical mean-field theory. The projection is exact in the limit of both infinitely large and small Coulomb interactions. For all parameter ranges the accuracy of the projected solution converges exponentially to the exact solution with increasing subsystem size.
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.
We have developed a new efficient and accurate impurity solver for the single impurity Anderson model (SIAM), which is based on a non-perturbative recursion technique in a space of operators and involves expanding the self-energy as a continued fraction. The method has no special occupation number or temperature restrictions; the only approximation is the number of levels of the continued fraction retained in the expansion. We also show how this approach can be used as a new approach to Dynamical Mean Field Theory (DMTF) and illustrate this with the Hubbard model. The three lowest orders of recursion give the Hartree-Fock, Hubbard I, and Hubbard III approximations. A higher level of recursion is able to reproduce the expected 3-peak structure in the spectral function and Fermi liquid behavior.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
