A generalized quantum master equation theory that governs the exact, nonperturbative quantum dissipation and quantum transport is formulated in terms of hierarchically coupled equations of motion for an arbitrary electronic system in contact with electrodes under either a stationary or a nonstationary electrochemical potential bias. The theoretical construction starts with the influence functional in path integral, in which the electron creation and annihilation operators are Grassmann variables. Time derivatives on the influence functionals are then performed in a hierarchical manner. Both the multiple-frequency dispersion and the non-Markovian reservoir parametrization schemes are considered for the desired hierarchy construction. The resulting hierarchical equations of motion formalism is in principle exact and applicable to arbitrary electronic systems, including Coulomb interactions, under the influence of arbitrary time-dependent applied bias voltage and external fields. Both the conventional quantum master equation and the real-time diagrammatic formalism of Schon and co-workers can be readily obtained at well defined limits of the present theory. We also show that for a noninteracting electron system, the present hierarchical equations of motion formalism terminates at the second tier exactly, and the Landuer-Buttiker transport current expression is recovered. The present theory renders an exact and numerically tractable tool to evaluate various transient and stationary quantum transport properties of many-electron systems, together with the involving nonperturbative dissipative dynamics.

An exact quantum master equation formalism is constructed for the efficient evaluation of quantum non-Markovian dissipation beyond the weak system-bath interaction regime in the presence of time-dependent external field. A novel truncation scheme is further proposed and compared with other approaches to close the resulting hierarchically coupled equations of motion. The interplay between system-bath interaction strength, non-Markovian property, and required level of hierarchy is also demonstrated with the aid of simple spin-boson systems.

We propose an efficient method to propagate the hierarchical quantum master equations based on a reformulation of the original formalism and the incorporation of a filtering algorithm that automatically truncates the hierarchy with a preselected tolerance. The new method is applied to calculate electron transfer dynamics in a spin-boson model and the absorption spectra of an excitonic dimmer. The proposed method significantly reduces the number of auxiliary density operators used in the hierarchical equation approach and thus provides an efficient way capable of studying real time dynamics of non-Markovian quantum dissipative systems in strong system-bath coupling and low temperature regimes.

Quantum dissipation involves both energy relaxation and decoherence, leading toward quantum thermal equilibrium. There are several theoretical prescriptions of quantum dissipation but none of them is simple enough to be treated exactly in real applications. As a result, formulations in different prescriptions are practically used with different approximation schemes. This review examines both theoretical and application aspects on various perturbative formulations, especially those that are exact up to second-order but nonequivalent in high-order system-bath coupling contributions. Discrimination is made in favor of an unconventional formulation that in a sense combines the merits of both the conventional time-local and memory-kernel prescriptions, where the latter is least favorite in terms of the applicability range of parameters for system-bath coupling, non-Markovian, and temperature. Also highlighted is the importance of correlated driving and dissipation effects, not only on the dynamics under strong external field driving, but also in the calculation of field-free correlation and response functions.

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.

Padé approximant is exploited for an efficient sum-over-poles decomposition of Fermi and Bose functions. The resulting poles are all pure imaginary and can therefore be used to define Padé frequencies, in analogy with the celebrated Matsubara frequencies. The proposed Padé spectrum decomposition is shown to be equivalent to a truncated continued fraction. It converges significantly faster than other schemes such as the Matsubara expansion at all temperatures. By introducing the characteristic validity length as the measure of approximant, we analyze the convergence properties of different schemes thoroughly. Our results qualify the present scheme the best among all sum-over-poles approaches. Thus, it is of great value in efficient numerical evaluations of integrals involving Fermi/Bose function in various condensed-phase matter problems.

Padé spectrum decomposition is an optimal sum-over-poles expansion scheme of Fermi function and Bose function [J. Hu, R. X. Xu, and Y. J. Yan, J. Chem. Phys. 133, 101106 (2010)]. In this work, we report two additional members to this family, from which the best among all sum-over-poles methods could be chosen for different cases of application. Methods are developed for determining these three Padé spectrum decomposition expansions at machine precision via simple algorithms. We exemplify the applications of present development with optimal construction of hierarchical equations-of-motion formulations for nonperturbative quantum dissipation and quantum transport dynamics. Numerical demonstrations are given for two systems. One is the transient transport current to an interacting quantum-dots system, together with the involved high-order co-tunneling dynamics. Another is the non-Markovian dynamics of a spin-boson system.

Using theory to guide the choice of pulse shape, we have synthesized frequency-chirped laser pulses and used them to control the evolution of vibrational wave packets on the 8 excited state of iodine.A negatively chirped pulse produces a wave packet at the The spread of the target in position and momentum is given by the formulas Ax = Qh/2m' and Ap = /mesh/2, where m is the reduced mass of the iodine molecule and the parameter co (units of frequency) was chosen to correspond to 250 cm '. The target is centered about a chosen position (0.372 nm) and a chosen momentum (corresponding to a kinetic energy of 50 meV).In this case, the momentum is negative, meaning that the iodine atoms are moving toward each other. The target is achieved by using a tailored, ultrashort light field to create an initial wave packet that evolves under the influ- state at a later time. The wave packet evolution is monitored by exciting the wave packet to the higher-lying e state with a second ultrashort pulse which opens up a "window" in internuclear distance (shown schematically by the striped box).The experimental signal is the total laser-induced fluorescence (LIF) from the e state as a function of delay between the pump and probe pulses, allowing detection of the wave packet density in a region of internuclear distance selected by the probe wavelength.ence of the molecular Hamiltonian to yield a wave packet that has maximum overlap with the target distribution at a later time. To reach the target in a reasonable amount of time, the electric field is allowed to be nonzero only during a temporal interval of our choice (550 fs, in the present case). Thus, the task of theory is to find the electric field F(t) that best steers the system to our chosen target at the chosen time. We note that to achieve this target, the optimal field must overcome the wave packet spreading which ordinarily occurs in anharmonic potentials. 33600031-9007/95/74(17)/3360(4)$06. 00