A hierarchical equations of motion based numerical approach is developed for accurate and efficient evaluation of dynamical observables of strongly correlated quantum impurity systems. This approach is capable of describing quantitatively Kondo resonance and Fermi-liquid characteristics, achieving the accuracy of the latest high-level numerical renormalization group approach, as demonstrated on single-impurity Anderson model systems. Its application to a two-impurity Anderson model results in differential conductance versus external bias, which correctly reproduces the continuous transition from Kondo states of individual impurity to singlet spin states formed between two impurities. The outstanding performance on characterizing both equilibrium and nonequilibrium properties of quantum impurity systems makes the hierarchical equations of motion approach potentially useful for addressing strongly correlated lattice systems in the framework of dynamical mean-field theory.
A hierarchical equations of motion based numerical approach is developed for accurate and efficient evaluation of dynamical observables of strongly correlated quantum impurity systems. This approach is capable of describing quantitatively Kondo resonance and Fermi-liquid characteristics, achieving the accuracy of the latest high-level numerical renormalization group approach, as demonstrated on single-impurity Anderson model systems. Its application to a two-impurity Anderson model results in differential conductance versus external bias, which correctly reproduces the continuous transition from Kondo states of individual impurity to singlet spin states formed between two impurities. The outstanding performance on characterizing both equilibrium and nonequilibrium properties of quantum impurity systems makes the hierarchical equations of motion approach potentially useful for addressing strongly correlated lattice systems in the framework of dynamical mean-field theory.
A nonperturbative quantum impurity solver is proposed based on a formally exact hierarchical equations of motion (HEOM) formalism for open quantum systems. It leads to quantitatively accurate evaluation of physical properties of strongly correlated electronic systems, in the framework of dynamical mean-field theory (DMFT). The HEOM method is also numerically convenient to achieve the same level of accuracy as that using the state-of-the-art numerical renormalization group impurity solver at finite temperatures. The practicality of the HEOM+DMFT method is demonstrated by its applications to the Hubbard models with Bethe and hypercubic lattice structures. We investigate the metal-insulator transition phenomena, and address the effects of temperature on the properties of strongly correlated lattice systems.
We propose a hierarchical dynamics approach for evaluation of nonequilibrium dynamic response properties of quantum impurity systems. It is based on a hierarchical equations of motion formalism, in conjunction with a linear response theory established in the hierarchical Liouville space. This provides an accurate and universal tool for characterization of a variety of response and correlation functions of local impurities, as well as transport related response properties. The practicality of our proposed approach is demonstrated via the evaluation of various dynamical properties of a single-impurity Anderson model. These include the impurity spectral density, local charge fluctuation, local magnetic susceptibility, and current-voltage admittance, in both equilibrium and nonequilibrium situations. The numerical results are considered to be quantitatively accurate, as long as they converge quantitatively with respect to the truncation of the hierarchy.
In recent years, nano-sized systems involving local charge or spin states have received wide interests because of their potential application in emerging fields such as quantum information and quantum computation. For instance, organometallic molecular complexes may serve as building blocks of quantum storage devices, because the spin-unpaired d or f electrons at transition metal centers may be employed to construct spin qubits. Moreover, if the system contains strong electron-electron interaction, the involving local quantum states are subject to prominent electron correlation effects (such as the Kondo effect). Theoretically, spatially confined nano-systems are often described by quantum impurity models. Thus the accurate prediction of the intrinsic properties of quantum impurity systems and the deep understanding on the response and evolution of local quantum states under external fields or in dissipative environment are fundamentally important for the design and fabrication of quantum devices. The accurate characterization of quantum coherence, correlation, and entanglement in quantum impurity systems remains a great challenge. Enormous efforts have been made to achieve this goal. A variety of theoretical methods have been developed, including the numerical renormalization group method, the quantum Monte Carlo method, and many others. However, all the existing methods are subject to certain limitations regarding accuracy or efficiency. Therefore, we choose to view this problem from a new perspective—the perspective of open quantum systems. Over the past decade, we have developed a formally exact quantum dissipation theory, the hierarchical equations of motion (HEOM) theory, for fermionic open systems. The HEOM theory captures the combined effects of system-environment dissipation, many-body interaction, and non-Markovian memory in a nonperturbative manner. It is capable of addressing static and dynamic responses of system observables in both equilibrium and nonequilibrium situations. We have implemented the HEOM method in our self-designed computer program HEOM-QUICK. We have also devised a series of advanced algorithms which substantially enhance the numerical efficiency of HEOM. The HEOM-QUICK program thus provides an accurate, efficient, and versatile theoretical tool for the investigation of strongly correlated quantum impurity systems. We have applied the HEOM method to study a variety of problems associated with quantum impurity systems. These include the intrinsic properties of quantum impurity models and lattice models, the transient current response of quantum dots to ac voltages, the thermopower and local heating effect in nonequilibrium quantum dots, and the tuning of local spin states in adsorbed molecular magnets. In particular, the HEOM method has been combined with density functional theory (DFT) method, and thus allows for first-principles-based simulation on the precise tuning of local spin states in adsorbed magnetic molecules. Numerical simulations have discovered or reproduced many important quantum phenomena that are essential for relevant experiments, including the Kondo effect, Mott metal-insulator transition, quantum memristive effect, etc. This paper gives a brief overview of our development of the HEOM method for fermionic open systems, as well as its applications to various strongly correlated quantum impurity systems.
Investigations of strongly correlated quantum impurity systems (QIS), which exhibit diversified novel and intriguing quantum phenomena, have become a highly concerning subject in recent years. The hierarchical equations of motion (HEOM) method is one of the most popular numerical methods to characterize QIS linearly coupled to the environment. This review provides a comprehensive account of a formally rigorous and numerical convergent HEOM method, including a modeling description of the QIS and an overview of the fermionic HEOM formalism. Moreover, a variety of spectrum decomposition schemes and hierarchal terminators have been proposed and developed, which significantly improve the accuracy and efficiency of the HEOM method, especially in cryogenic temperature regimes. The practicality and usefulness of the HEOM method to tackle strongly correlated issues are exemplified by numerical simulations for the characterization of nonequilibrium quantum transport and strongly correlated Kondo states as well as the investigation of nonequilibrium quantum thermodynamics.
We present a numerical method for the study of correlated quantum impurity problems out of equilibrium, which is particularly suited to address steady-state properties within dynamical mean field theory. The approach, recently introduced by Arrigoni et al. [Phys. Rev. Lett. 110, 086403 (2013)], is based upon a mapping of the original impurity problem onto an auxiliary open quantum system, consisting of the interacting impurity coupled to bath sites as well as to a Markovian environment. The dynamics of the auxiliary system is governed by a Lindblad master equation whose parameters are used to optimize the mapping. The accuracy of the results can be readily estimated and systematically improved by increasing the number of auxiliary bath sites, or by introducing a linear correction. Here, we focus on a detailed discussion of the proposed approach including technical remarks. To solve for the Green's functions of the auxiliary impurity problem, a non-Hermitian Lanczos diagonalization is applied. As a benchmark, results for the steady-state current-voltage characteristics of the single-impurity Anderson model are presented. Furthermore, the bias dependence of the single-particle spectral function and the splitting of the Kondo resonance are discussed. In its present form, the method is fast, efficient, and features a controlled accuracy.
Understanding the time evolution of strongly correlated open quantum systems (OQSs) in response to perturbations (quenches) is of fundamental importance to the precise control of quantum devices. It is, however, rather challenging in multi-impurity quantum systems because such evolution often involves multiple intricate dynamical processes. In this work, we apply the numerically exact hierarchical equations of motion approach to explore the influence of two different types of perturbations, i.e., sudden swapping of the energy levels of impurity systems and activating the inter-impurity spin-exchange interaction, on the dissipation dynamics of the Kondo-correlated two-impurity Anderson model over a wide range of energetic parameters. By evaluating the time-dependent impurity spectral function and other system properties, we analyze the time evolution of the Kondo state in detail and conclude a phenomenologically scaling relation for Kondo dynamics driven by these perturbations. The evolutionary scaling relationship is not only related to the Kondo characteristic energy TK but also significantly affected by the simultaneous non-Kondo dynamic characteristic energy. We expect these results will inspire subsequent theoretical studies on the dynamics of strongly correlated OQSs.