HEOM-QUICK: a program for accurate, efficient, and universal characterization of strongly correlated quantum impurity systems

Hierarchical Equations of Motion for Quantum Chemical Dynamics: Recent Methodology Developments and Applications

ConspectusQuantum effects are critical to understanding many chemical dynamical processes in condensed phases, where interactions between molecules and their environment are usually strong and non-Markovian. In this Account, we review recent progress from our group in development and application of the hierarchical equations of motion (HEOM) method, highlighting its ability to address some challenging problems in quantum chemical dynamics.In the HEOM method, the bath degrees of freedom are represented using effective modes from exponential decomposition of the bath correlation function. Complex spectral densities and low temperature simulations often require a larger number of modes, making the simulations very expensive. Recent advances, such as the barycentric spectral decomposition (BSD) technique, can significantly reduce the number of effective modes, allowing to handle complex spectral densities and enabling simulations at very low temperatures, including near-zero temperature dynamics.Another key improvement in the computational efficiency is the use of tensor network methods like matrix product states and hierarchical tensor networks. These techniques allow for efficient HEOM propagation with thousands of effective modes, crucial for simulating large molecular systems interacting with multiple baths. This combination enables simulations of excitation energy transfer (EET) in systems like the Fenna-Matthews-Olson (FMO) complex and even larger systems with experimentally determined spectral densities.The versatility of the HEOM method is demonstrated through applications to a wide range of chemical dynamics problems. Simulations of EET and related ultrafast spectroscopy are first briefly covered. Applications of the HEOM to quantum tunneling effects in proton transfer reactions are then presented. Early works have studied the non-Kramers dependence of the rate constant as a function of bath friction due to deep tunneling and revealed vibrationally nonadiabatic dynamics within the so-called nontraditional view of proton transfer reactions. A recent work on the large kinetic isotope effects in soybean lipoxygenase also indicated that many quantum correction approximations to classical transition-state theory may fall short in describing deep tunneling effects.Charge transport and separation dynamics in organic semiconductors are another area where the HEOM method has been instrumental. We first demonstrate that the HEOM provides a unified description of both band-like and thermally assisted charge carrier transport in organic materials. The effect of non-nearest neighbor transitions is then investigated by combining generalized master equations with exact memory kernels. The HEOM method also enables simulation of charge separation in organic photovoltaics (OPVs) and reveals how factors such as external electric fields, entropy, and charge delocalization influence the charge separation barrier and dynamics.Moreover, HEOM has been applied to investigate hydrogen atom scattering on the Au(111) surface and vibrational energy relaxation at molecule-metal interfaces. These studies provide deeper insights into how electron-hole pair excitations and temporary charge transfer states influence the nuclear motion, offering a new framework for simulating nonadiabatic dynamics on metal surfaces.In summary, the HEOM method has developed into a robust tool for simulating quantum effects in condensed phases. Future developments in algorithm efficiency and computational power will likely expand its applicability to even more complex systems.

Nonadiabatic dynamics of molecules interacting with metal surfaces: Extending the hierarchical equations of motion and Langevin dynamics approach to position-dependent metal-molecule couplings

Electronic friction and Langevin dynamics is a popular mixed quantum-classical method for simulating the nonadiabatic dynamics of molecules interacting with metal surfaces, as it can be computationally more efficient than fully quantum approaches. In this work, we extend the theory of electronic friction within the hierarchical equations of motion formalism to models with a position-dependent metal-molecule coupling. We show that the addition of a position-dependent metal-molecule coupling adds new contributions to the electronic friction and other forces, which are highly relevant for many physical processes. Our expressions for the electronic forces within the Langevin equation are valid both in and out of equilibrium and for molecular models containing strong interactions. We demonstrate the approach by applying it to different models of interest.

Differences in Kondo Splitting of Surface Quantum Systems Induced by Two Distinct Magnetic Tips: A Joint Method of DFT and HEOM

The interactions between a magnetic tip and local spin impurities initiate unconventional Kondo phenomena, such as asymmetric suppression or even splitting of the Kondo peak. However, a lack of realistic theoretical models and comprehensive explanations for this phenomenon persists due to the complexity of the interactions. This research employs a joint method of density functional theory (DFT) and hierarchical equation of motion (HEOM) to simulate and contrast the modulation of the spin state and Kondo behavior in the Fe/Cu(100) system with two distinct magnetic tips. A cobalt tip, possessing a larger magnetic moment, incites greater atomic displacement of the iron atom, more notable alterations in electronic structure, and enhanced charge transfer with the environment compared with the control process utilizing a nickel tip. Furthermore, the Kondo resonance undergoes asymmetric splitting as a result of the ferromagnetic correlation between the iron atom and the magnetic tip. The Co tip's higher spin polarization results in a wider spacing between the splitting peaks. This investigation underscores the precision of the DFT + HEOM approach in predicting complex quantum phenomena and explaining the underlying physical principles. This provides valuable theoretical support for developing more sophisticated quantum regulation experiments.

Unraveling the Nature of Spin Coupling in a Metal-Free Diradical: Theoretical Distinction of Ferromagnetic and Antiferromagnetic Interactions

Metal-free diradicals based on polycyclic aromatic hydrocarbons are promising candidates for organic spintronics due to their stable magnetism and tunable spin coupling. However, distinguishing and elucidating the origins of ferromagnetic and antiferromagnetic interactions in these systems remain challenging. Here, we investigate the 2-OS diradical molecule sandwiched between gold electrodes using a combined density functional theory and hierarchical equations of motion approach. We find that the dihedral angle between the radical moieties controls the nature and strength of the intramolecular spin coupling, transitioning smoothly from antiferromagnetic to ferromagnetic as the angle increases. Distinct features in the inelastic electron tunneling spectra are identified that can discern the two coupling regimes, including spin excitation steps whose energies directly reveal the exchange coupling constant. Mechanical stretching of the junction is predicted to modulate the spectral line shapes by adjusting the hybridization of the molecular radicals with the electrodes. Our work elucidates the electronic origin of tunable spin interactions in 2-OS and provides spectroscopic fingerprints for characterizing magnetism in metal-free diradicals.

Seeking a quantum advantage with trapped-ion quantum simulations of condensed-phase chemical dynamics

Simulating the quantum dynamics of molecules in the condensed phase represents a longstanding challenge in chemistry. Trapped-ion quantum systems may serve as a platform for the analog-quantum simulation of chemical dynamics that is beyond the reach of current classical-digital simulation. To identify a 'quantum advantage' for these simulations, performance analysis of both analog-quantum simulation on noisy hardware and classical-digital algorithms is needed. In this Review, we make a comparison between a noisy analog trapped-ion simulator and a few choice classical-digital methods on simulating the dynamics of a model molecular Hamiltonian with linear vibronic coupling. We describe several simple Hamiltonians that are commonly used to model molecular systems, which can be simulated with existing or emerging trapped-ion hardware. These Hamiltonians may serve as stepping stones towards the use of trapped-ion simulators for systems beyond the reach of classical-digital methods. Finally, we identify dynamical regimes in which classical-digital simulations seem to have the weakest performance with respect to analog-quantum simulations. These regimes may provide the lowest hanging fruit to make the most of potential quantum advantages.

Open quantum systems with nonlinear environmental backactions: Extended dissipaton theory vs core-system hierarchy construction

In this paper, we present a comprehensive account of quantum dissipation theories with the quadratic environment couplings. The theoretical development includes the Brownian solvation mode embedded hierarchical quantum master equations, a core-system hierarchy construction that verifies the extended dissipaton equation of motion (DEOM) formalism [R. X. Xu et al., J. Chem. Phys. 148, 114103 (2018)]. Developed are also the quadratic imaginary-time DEOM for equilibrium and the λ(t)-DEOM for nonequilibrium thermodynamics problems. Both the celebrated Jarzynski equality and Crooks relation are accurately reproduced, which, in turn, confirms the rigorousness of the extended DEOM theories. While the extended DEOM is more numerically efficient, the core-system hierarchy quantum master equation is favorable for "visualizing" the correlated solvation dynamics.

Unveiling the Decisive Factor for the Sharp Transition in the Scanning Tunneling Spectroscopy of a Single Nickelocene Molecule

Scanning tunneling microscopy (STM) has been utilized to realize the precise measurement and control of local spin states. Experiments have demonstrated that when a nickelocene (Nc) molecule is attached to the apex of an STM tip, the dI/dV spectra exhibit a sharp or a smooth transition when the tip is displaced toward the substrate. However, what leads to the two distinct types of transitions remains unclear, and more intriguingly, the physical origin of the abrupt change in the line shape of dI/dV spectra remains unclear. To clarify these intriguing issues, we perform first-principles-based simulations on the STM tip control process for the Cu tip/Nc/Cu(100) junction. In particular, we find that the suddenly enhanced hybridization between the d orbitals on the Ni ion and the metallic bands in the substrate leads to Kondo correlation overwhelming spin excitation, which is the main cause of the sharp transition in the dI/dV spectra observed experimentally.

Nonequilibrium reaction rate theory: Formulation and implementation within the hierarchical equations of motion approach

The study of chemical reactions in environments under nonequilibrium conditions has been of interest recently in a variety of contexts, including current-induced reactions in molecular junctions and scanning tunneling microscopy experiments. In this work, we outline a fully quantum mechanical, numerically exact approach to describe chemical reaction rates in such nonequilibrium situations. The approach is based on an extension of the flux correlation function formalism to nonequilibrium conditions and uses a mixed real and imaginary time hierarchical equations of motion approach for the calculation of rate constants. As a specific example, we investigate current-induced intramolecular proton transfer reactions in a molecular junction for different applied bias voltages and molecule-lead coupling strengths.

Universal time-domain Prony fitting decomposition for optimized hierarchical quantum master equations

In this work, we propose the time-domain Prony fitting decomposition (t-PFD) as an accurate and effcient exponential series method, applicable to arbitrary bath correlation functions. The resulting numerical effciency of hierarchical equations of motion (HEOM) formalism is greatly optimized, especially in low temperature regimes that would be inaccessible with other methods. For demonstration, we calibrate the present t-PFD against the celebrated Padé spectrum decomposition (PSD) method, followed by converged HEOM evaluations on the single-impurity Anderson model system.

Hierarchical equations of motion approach to hybrid fermionic and bosonic environments: Matrix product state formulation in twin space

We extend the twin-space formulation of the hierarchical equations of motion approach in combination with the matrix product state representation (introduced in J. Chem. Phys. 150, 234102, [2019]) to nonequilibrium scenarios where the open quantum system is coupled to a hybrid fermionic and bosonic environment. The key ideas used in the extension are a reformulation of the hierarchical equations of motion for the auxiliary density matrices into a time-dependent Schrödinger-like equation for an augmented multi-dimensional wave function as well as a tensor decomposition into a product of low-rank matrices. The new approach facilitates accurate simulations of non-equilibrium quantum dynamics in larger and more complex open quantum systems. The performance of the method is demonstrated for a model of a molecular junction exhibiting current-induced mode-selective vibrational excitation.

Study the mixed valence problem in asymmetric Anderson model: Fano-Kondo resonance around Fermi level

We numerically calculate the local density of states (LDOS) in asymmetric Anderson model in mixed valence regime using hierarchical equations of motion approach. Based on the idea that the asymmetric line shape of LDOS around Fermi level stems from the interference between the single particle resonance and the Kondo resonance, we perform a fitting. From the fitting results, we obtain the Kondo temperatures and the Fano factors with changing the single particle energy. The tendency of Kondo temperature agrees with the previous analytic expressions and the Fano factors are in an expected variation of Fano resonance. Our study shows that the Fano-Kondo resonance can reasonably explain the asymmetric line shape of the LDOS around the Fermi level.

Origin of Asymmetric Splitting of Kondo Peak in Spin-Polarized Scanning Tunneling Spectroscopy: Insights from First-Principles-Based Simulations

The spin-polarized scanning tunneling microscope (SP-STM) has served as a versatile tool for probing and manipulating the spintronic properties of atomic and molecular devices with high precision. The interplay between the local spin state and its surrounding magnetic environment significantly affects the transport behavior of the device. Particularly, in the contact regime, the strong hybridization between the SP-STM tip and the magnetic atom or molecule could give rise to unconventional Kondo resonance signatures in the differential conductance (dI/dV) spectra. This poses challenges for the simulation of a realistic tip control process. By combining the density functional theory and the hierarchical equations of motion methods, we achieve first-principles-based simulation of the control of a Ni-tip/Co/Cu(100) junction in both the tunneling and contact regimes. The calculated dI/dV spectra reproduce faithfully the experimental data. A cotunneling mechanism is proposed to elucidate the physical origin of the observed unconventional Kondo signatures.

Nonequilibrium Kondo regime current noise spectrum of quantum dot systems with the single impurity Anderson model

We investigate the nonequilibrium current noise spectrum of single impurity Anderson model quantum dot systems on the basis of the accurate dissipation equation of motion evaluations. By comparing between the equilibrium and nonequilibrium cases and between the non-Kondo and Kondo regimes, we identify the current noise spectrum of the nonequilibrium Kondo features that actually appear in the entire region of ω ∈ [-eV, eV]. It is well known that the primary Kondo characteristics at ω = ±eV = ±(μ_L - μ_R) display asymmetrical upturns and remarkable peaks in S(ω) and dS(ω)/dω, respectively. These features are originated from the Rabi interference of the transport current dynamics, with the Kondo oscillation frequency of |eV|. Moreover, we also identify the minor but very distinguishable inflections, crossing over from ω = -eV to ω = +eV. This uncovered feature would be related to the interference between two Kondo resonance channels.

Correlated vibration-solvent effects on the non-Condon exciton spectroscopy

Excitation energy transfer is crucially involved in a variety of systems. During the process, the non-Condon vibronic coupling and the surrounding solvent interaction may synergetically play important roles. In this work, we study the correlated vibration-solvent influences on the non-Condon exciton spectroscopy. Statistical analysis is elaborated for the overall vibration-plus-solvent environmental effects. Analytic solutions are derived for the linear absorption of monomer systems. General simulations are accurately carried out via the dissipaton-equation-of-motion approach. The resulted spectra in either the linear absorption or strong field regime clearly demonstrate the coherence enhancement due to the synergetic vibration-solvent correlation.

Unraveling current-induced dissociation mechanisms in single-molecule junctions

Understanding current-induced bond rupture in single-molecule junctions is both of fundamental interest and a prerequisite for the design of molecular junctions, which are stable at higher-bias voltages. In this work, we use a fully quantum mechanical method based on the hierarchical quantum master equation approach to analyze the dissociation mechanisms in molecular junctions. Considering a wide range of transport regimes, from off-resonant to resonant, non-adiabatic to adiabatic transport, and weak to strong vibronic coupling, our systematic study identifies three dissociation mechanisms. In the weak and intermediate vibronic coupling regime, the dominant dissociation mechanism is stepwise vibrational ladder climbing. For strong vibronic coupling, dissociation is induced via multi-quantum vibrational excitations triggered either by a single electronic transition at high bias voltages or by multiple electronic transitions at low biases. Furthermore, the influence of vibrational relaxation on the dissociation dynamics is analyzed and strategies for improving the stability of molecular junctions are discussed.

Efficient Steady-State Solver for the Hierarchical Equations of Motion Approach: Formulation and Application to Charge Transport through Nanosystems

An iterative approach is introduced, which allows the efficient solution of the hierarchical equations of motion (HEOM) for the steady-state of open quantum systems. The approach combines the method of matrix equations with an efficient preconditioning technique to reduce the numerical effort of solving the HEOM. Illustrative applications to simulate nonequilibrium charge transport in single-molecule junctions demonstrate the performance of the method.

Numerically "exact" approach to open quantum dynamics: The hierarchical equations of motion (HEOM)

An open quantum system refers to a system that is further coupled to a bath system consisting of surrounding radiation fields, atoms, molecules, or proteins. The bath system is typically modeled by an infinite number of harmonic oscillators. This system-bath model can describe the time-irreversible dynamics through which the system evolves toward a thermal equilibrium state at finite temperature. In nuclear magnetic resonance and atomic spectroscopy, dynamics can be studied easily by using simple quantum master equations under the assumption that the system-bath interaction is weak (perturbative approximation) and the bath fluctuations are very fast (Markovian approximation). However, such approximations cannot be applied in chemical physics and biochemical physics problems, where environmental materials are complex and strongly coupled with environments. The hierarchical equations of motion (HEOM) can describe the numerically "exact" dynamics of a reduced system under nonperturbative and non-Markovian system-bath interactions, which has been verified on the basis of exact analytical solutions (non-Markovian tests) with any desired numerical accuracy. The HEOM theory has been used to treat systems of practical interest, in particular, to account for various linear and nonlinear spectra in molecular and solid state materials, to evaluate charge and exciton transfer rates in biological systems, to simulate resonant tunneling and quantum ratchet processes in nanodevices, and to explore quantum entanglement states in quantum information theories. This article presents an overview of the HEOM theory, focusing on its theoretical background and applications, to help further the development of the study of open quantum dynamics.

Molecular molds for regularizing Kondo states at atom/metal interfaces

Adsorption of magnetic transition metal atoms on a metal surface leads to the formation of Kondo states at the atom/metal interfaces. However, the significant influence of surrounding environment presents challenges for potential applications. In this work, we realize a novel strategy to regularize the Kondo states by moving a CoPc molecular mold on an Au(111) surface to capture the dispersed Co adatoms. The symmetric and ordered structures of the atom-mold complexes, as well as the strong dπ-π bonding between the Co adatoms and conjugated isoindole units, result in highly robust and uniform Kondo states at the Co/Au(111) interfaces. Even more remarkably, the CoPc further enables a fine tuning of Kondo states through the molecular-mold-mediated superexchange interactions between Co adatoms separated by more than 12 Å. Being highly precise, efficient and reproducible, the proposed molecular mold strategy may open a new horizon for the construction and control of nano-sized quantum devices.

Current-induced effective Dzyaloshinskii-Moriya interaction and its Kondo enhancement in double quantum dot

We studied the nonequilibrium transport of serially coupled double quantum dots connected to ferromagnetic electrodes. We demonstrated that the nonadiabatic part of the spin gauge field resulted in a current-induced Dzyaloshinskii-Moriya (DM) interaction effect in a double quantum dot and numerically confirmed this observation through the hierarchical equations of motion approach. We report that the spin current and the effective DM interaction are enhanced in the Kondo regime. We demonstrate that this enhancement occurs because the Kondo resonance, which is supposed to be suppressed by the local ferromagnetic exchange, is enhanced by the inter-dot coupling. This additional Kondo resonance channel increases the spin current. In addition, the impact of the spin-spin interaction and the Kondo effect on tunnel magnetoresistance is discussed. Our results offer a new approach for controlling the non-collinear spin interaction in double quantum dot devices.

Hierarchical equations of motion method based on Fano spectrum decomposition for low temperature environments

The hierarchical equations of motion (HEOM) method has become one of the most popular methods for the studies of the open quantum system. However, its applicability to systems at ultra-low temperatures is largely restrained by the enormous computational cost, which is caused by the numerous exponential functions required to accurately characterize the non-Markovian memory of the reservoir environment. To overcome this problem, a Fano spectrum decomposition (FSD) scheme has been proposed recently [Cui et al., J. Chem. Phys. 151, 024110 (2019)], which expands the reservoir correlation functions using polynomial-exponential functions and hence greatly reduces the size of the memory basis set. In this work, we explicitly establish the FSD-based HEOM formalisms for both bosonic and fermionic environments. The accuracy and efficiency of the FSD-based HEOM are exemplified by the calculated low-temperature dissipative dynamics of a spin-boson model and the dynamic and static properties of a single-orbital Anderson impurity model in the Kondo regime. The encouraging numerical results highlight the practicality and usefulness of the FSD-based HEOM method for general open systems at ultra-low temperatures.

Extending the hierarchical quantum master equation approach to low temperatures and realistic band structures

The hierarchical quantum master equation (HQME) approach is an accurate method to describe quantum transport in interacting nanosystems. It generalizes perturbative master equation approaches by including higher-order contributions as well as non-Markovian memory and allows for the systematic convergence to the numerically exact result. As the HQME method relies on a decomposition of the bath correlation function in terms of exponentials, however, its application to systems at low temperatures coupled to baths with complexer band structures has been a challenge. In this publication, we outline an extension of the HQME approach, which uses re-summation over poles and can be applied to calculate transient currents at a numerical cost that is independent of temperature and band structure of the baths. We demonstrate the performance of the extended HQME approach for noninteracting tight-binding model systems of increasing complexity as well as for the spinless Anderson-Holstein model.

On the exact truncation tier of fermionic hierarchical equations of motion

The hierarchical equations of motion (HEOM) theory is in principle exact for describing the dissipative dynamics of quantum systems linearly coupled to Gaussian environments. In practice, the hierarchy needs to be truncated at a finite tier. We demonstrate that, for general systems described by the fermionic HEOM, the (n+L̃)th-tier truncation with L̃=2N_σN_ν yields the exact density operators up to the nth tier. Here, N_σ = 2 for fermionic systems and N_ν is the system degrees of freedom. For noninteracting systems, L̃ is further reduced by half. Such an exact termination pattern originates from the Pauli exclusion principle for fermions, and it holds true regardless of the system-environment coupling strength, the number of coupling reservoirs, or the specific scheme employed to unravel the environment memory contents. The relatively small L̃ emphasizes the nonperturbative nature of the HEOM theory. We also propose a simplified HEOM approach to further reduce the memory cost for practical calculations.

A unified stochastic formulation of dissipative quantum dynamics. I. Generalized hierarchical equations

We extend a standard stochastic theory to study open quantum systems coupled to a generic quantum environment. We exemplify the general framework by studying a two-level quantum system coupled bilinearly to the three fundamental classes of non-interacting particles: bosons, fermions, and spins. In this unified stochastic approach, the generalized stochastic Liouville equation (SLE) formally captures the exact quantum dissipations when noise variables with appropriate statistics for different bath models are applied. Anharmonic effects of a non-Gaussian bath are precisely encoded in the bath multi-time correlation functions that noise variables have to satisfy. Starting from the SLE, we devise a family of generalized hierarchical equations by averaging out the noise variables and expand bath multi-time correlation functions in a complete basis of orthonormal functions. The general hierarchical equations constitute systems of linear equations that provide numerically exact simulations of quantum dynamics. For bosonic bath models, our general hierarchical equation of motion reduces exactly to an extended version of hierarchical equation of motion which allows efficient simulation for arbitrary spectral densities and temperature regimes. Similar efficiency and flexibility can be achieved for the fermionic bath models within our formalism. The spin bath models can be simulated with two complementary approaches in the present formalism. (I) They can be viewed as an example of non-Gaussian bath models and be directly handled with the general hierarchical equation approach given their multi-time correlation functions. (II) Alternatively, each bath spin can be first mapped onto a pair of fermions and be treated as fermionic environments within the present formalism.

Precise Control of Local Spin States in an Adsorbed Magnetic Molecule with an STM Tip: Theoretical Insights from First-Principles-Based Simulation

The precise tuning of local spin states in adsorbed organometallic molecules by a mechanically controlled scanning tunneling microscope (STM) tip has become a focus of recent experiments. However, the underlying mechanisms remain somewhat unclear. We investigate theoretically the STM tip control of local spin states in a single iron(II) porphyrin molecule adsorbed on the Pb(111) substrate. A combined density functional theory and hierarchical equations of motion approach is employed to simulate the tip tuning process in conjunction with the complete active space self-consistent field method for accurate computation of magnetic anisotropy. Our first-principles-based simulation accurately reproduces the tuning of magnetic anisotropy realized in experiment. Moreover, we elucidate the evolution of geometric and electronic structures of the composite junction and disclose the delicate competition between the Kondo resonance and local spin excitation. The understanding and insight provided by the first-principles-based simulation may help to realize more fascinating quantum state manipulations.

Manipulation of spin and magnetic anisotropy in bilayer magnetic molecular junctions

The Kondo effect and magnetic anisotropy in bilayer TMPc/TMPc/Pb(111) junctions can be actively tuned by changing the intermediate decoupling layer.

Ferromagnetic Phase in Nonequilibrium Quantum Dots

By nonperturbatively solving the nonequilibrium Anderson two-impurity model with the hierarchical equations of motion approach, we report a robust ferromagnetic (FM) phase in series-coupled double quantum dots, which can suppress the antiferromagnetic (AFM) phase and dominate the phase diagram at finite bias and detuning energy in the strongly correlated limit. The FM exchange interaction origins from the passive parallel spin arrangement caused by the Pauli exclusion principle during the electrons transport. At very low temperature, the Kondo screening of the magnetic moment in the FM phase induces some nonequilibrium Kondo effects in magnetic susceptibility, spectral functions and current. In the weakly correlated limit, the AFM phase is found still stable, therefore, a magnetic-field-free internal control of spin states can be expected through the continuous FM–AFM phase transition.

Many-body Tunneling and Nonequilibrium Dynamics of Doublons in Strongly Correlated Quantum Dots

Quantum tunneling dominates coherent transport at low temperatures in many systems of great interest. In this work we report a many–body tunneling (MBT), by nonperturbatively solving the Anderson multi-impurity model, and identify it a fundamental tunneling process on top of the well–acknowledged sequential tunneling and cotunneling. We show that the MBT involves the dynamics of doublons in strongly correlated systems. Proportional to the numbers of dynamical doublons, the MBT can dominate the off–resonant transport in the strongly correlated regime. A T^3/2–dependence of the MBT current on temperature is uncovered and can be identified as a fingerprint of the MBT in experiments. We also prove that the MBT can support the coherent long–range tunneling of doublons, which is well consistent with recent experiments on ultracold atoms. As a fundamental physical process, the MBT is expected to play important roles in general quantum systems.

Corrected Kondo temperature beyond the conventional Kondo scaling limit

In the Kondo systems such as the magnetic impurity screened by the conduction electrons in a metal host, as well as the quantum dots connected by the leads, the low energy behaviors have universal dependence on the [Formula: see text] or [Formula: see text], where [Formula: see text] is the conventional Kondo temperature. However, it was shown that this scaling behavior is only valid at low-energy; this is called the Kondo scaling limit. Here we explore the extention of the scaling parameter range by introducing the corrected Kondo temperature T _K, which may depend on the temperature and bias, as well as the other external parameters. We define the corrected Kondo temperature by scaling the local density of states near the Fermi level, obtained by accurate hierarchy of equations of motion approach at finite temperature and finite bias, and thus obtain a phenomenological expression of the corrected Kondo temperature. By using the corrected Kondo temperature as a characteristic energy scale, the conductance of the quantum dot can be well scaled in a wide parameter range, even two orders beyond the conventional scaling parameter range. Our work indicates that the Kondo scaling, although dominated by the conventional Kondo temperature in the low-energy of the Kondo system, could be extended to a higher energy regime, which is useful for analyzing the physics of the Kondo transport in non-equilibrium or high temperature cases.

Fokker-Planck quantum master equation for mixed quantum-semiclassical dynamics

We revisit Caldeira-Leggett's quantum master equation representing mixed quantum-classical theory, but with limited applications. Proposed is a Fokker-Planck quantum master equation theory, with a generic bi-exponential correlation function description on semiclassical Brownian oscillators' environments. The new theory has caustic terms that bridge between the quantum description on primary systems and the semiclassical or quasi-classical description on environments. Various parametrization schemes, both analytical and numerical, for the generic bi-exponential environment bath correlation functions are proposed and scrutinized. The Fokker-Planck quantum master equation theory is of the same numerical cost as the original Caldeira-Leggett's approach but acquires a significantly broadened validity and accuracy range, as illustrated against the exact dynamics on model systems in quantum Brownian oscillators' environments, at moderately low temperatures.