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.
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.
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.
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.
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.
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.
A quantum dissipation theory is formulated in terms of hierarchically coupled equations of motion for an arbitrary fermionic system coupled with grand canonical fermionic bath ensembles.1 The theoretical construction starts with the second-quantization influence functional in path integral formalism, in which the fermionic creation and annihilation operators are represented by Grassmann variables. Temporal derivatives on influence functionals are then performed in a hierarchical manner, on the basis of calculus-on-path-integral algorithm. Both the multiple-frequency-dispersion and the non-Markovian reservoir parametrization schemes are considered for the desired hierarchy construction. The resulting formalism is in principle exact, applicable to interacting systems, with arbitrary time-dependent external fields. It renders an exact tool to evaluate various transient and stationary quantum transport properties of many-electron systems. At the second-tier truncation level the present theory recovers the real-time diagrammatic formalism developed by Schön and coworkers.2 For a single-particle system, the hierarchical formalism terminates at the second tier exactly, and the Landuer-Büttiker's transport current expression is readily recovered. Numerical studies will be presented to highlight the richness of transient current through both interacting and noninteracting model systems.
Recent technological advancement in scanning tunneling microscopes has enabled the measurement of spin-field and spin-spin interactions in single atomic or molecular junctions with an unprecedentedly high resolution. Theoretically, although the fermionic hierarchical equations of motion (HEOM) method has been widely applied to investigate the strongly correlated Kondo states in these junctions, the existence of low-energy spin excitations presents new challenges to numerical simulations. These include the quest for a more accurate and efficient decomposition for the non-Markovian memory of low-temperature environments and a more careful handling of errors caused by the truncation of the hierarchy. In this work, we propose several new algorithms, which significantly enhance the performance of the HEOM method, as exemplified by the calculations on systems involving various types of low-energy spin excitations. Being able to characterize both the Kondo effect and spin excitation accurately, the HEOM method offers a sophisticated and versatile theoretical tool, which is valuable for the understanding and even prediction of the fascinating quantum phenomena explored in cutting-edge experiments.