The simulation of quantum transport in a realistic, many-particle system is a nontrivial problem with no quantitatively satisfactory solution. While real-time propagation has the potential to overcome the shortcomings of conventional transport methods, this approach is prone to finite size effects that are associated with modeling an open system on a closed spatial domain. Using a master equation framework, we exploit an equivalence between the superoperators coupling an open system to external particle reservoirs and non-Hermitian terms defined at the periphery of a quantum device. By taking the mean-field limit, the equation of motion for the single-particle reduced density matrix becomes equivalent to real-time time-dependent density functional theory in the presence of imaginary source and sink potentials. This method may be used to converge nonequilibrium steady states for a many-body quantum system using a previously reported constraint algorithm.
The simulation of quantum transport in a realistic, many-particle system is a nontrivial problem with no quantitatively satisfactory solution. While real-time propagation has the potential to overcome the shortcomings of conventional transport methods, this approach is prone to finite size effects that are associated with modeling an open system on a closed spatial domain. Using a master equation framework, we exploit an equivalence between the superoperators coupling an open system to external particle reservoirs and non-Hermitian terms defined at the periphery of a quantum device. By taking the mean-field limit, the equation of motion for the single-particle reduced density matrix becomes equivalent to real-time time-dependent density functional theory in the presence of imaginary source and sink potentials. This method may be used to converge nonequilibrium steady states for a many-body quantum system using a previously reported constraint algorithm.
The accurate simulation of real--time quantum transport is notoriously difficult, requiring a consistent scheme to treat incoming and outgoing fluxes at the boundary of an open system. We demonstrate a method to converge non--equilibrium steady states using non--Hermitian source and sink potentials alongside the application of dynamic density constraints during wavefunction propagation. This scheme adds negligible cost to existing computational methods while exhibiting exceptional stability and numerical accuracy.
Basing on our hierarchical equations of motion for time-dependent quantum transport [X. Zheng, G. H. Chen, Y. Mo, S. K. Koo, H. Tian, C. Y. Yam, and Y. J. Yan, J. Chem. Phys. 133, 114101 (2010)], we develop an efficient and accurate numerical algorithm to solve the Liouville-von-Neumann equation. We solve the real-time evolution of the reduced single-electron density matrix at the tight-binding level. Calculations are carried out to simulate the transient current through a linear chain of atoms, with each represented by a single orbital. The self-energy matrix is expanded in terms of multiple Lorentzian functions, and the Fermi distribution function is evaluated via the Padè spectrum decomposition. This Lorentzian-Padè decomposition scheme is employed to simulate the transient current. With sufficient Lorentzian functions used to fit the self-energy matrices, we show that the lead spectral function and the dynamics response can be treated accurately. Compared to the conventional master equation approaches, our method is much more efficient as the computational time scales cubically with the system size and linearly with the simulation time. As a result, the simulations of the transient currents through systems containing up to one hundred of atoms have been carried out. As density functional theory is also an effective one-particle theory, the Lorentzian-Padè decomposition scheme developed here can be generalized for first-principles simulation of realistic systems.
Time-dependent [current]-density functional theory for open quantum systems (OQS) has emerged as a formalism that can incorporate dissipative effects in the dynamics of many-body quantum systems. Here, we review and clarify some formal aspects of these theories that have been recently questioned in the literature. In particular, we provide theoretical support for the following conclusions: (1) contrary to what we and others had stated before, within the master equation framework, there is in fact a one-to-one mapping between vector potentials and current densities for fixed initial state, particle-particle interaction, and memory kernel; (2) regardless of the first conclusion, all of our recently suggested Kohn-Sham (KS) schemes to reproduce the current and particle densities of the original OQS, and in particular, the use of a KS closed driven system, remains formally valid; (3) the Lindblad master equation maintains the positivity of the density matrix regardless of the time-dependence of the Hamiltonian or the dissipation operators; (4) within the stochastic Schrödinger equation picture, a one-to-one mapping from stochastic vector potential to stochastic current density for individual trajectories has not been proven so far, except in the case where the vector potential is the same for every member of the ensemble, in which case, it reduces to the Lindblad master equation picture; (5) master equations may violate certain desired properties of the density matrix, such as positivity, but they remain as one of the most useful constructs to study OQS when the environment is not easily incorporated explicitly in the calculation. The conclusions support our previous work as formally rigorous, offer new insights into it, and provide a common ground to discuss related theories.
We propose a simple, yet feasible, model for quantum transport of fermionic carriers across tight-binding chain connecting two reservoirs maintained at arbitrary temperatures and chemical potentials. The model allows for elementary derivation of the master equation for the reduced single particle density matrix (SPDM) in a closed form in both Markov and Born approximations. In the Markov approximation the master equation is solved analytically, whereas in the Born approximation the problem is reduced to an algebraic equation for the SPDM in the Redfield form. The non-Markovian equation is shown to lead to resonant transport similar to Landauer's conductance. It is shown that in the deep non-Markovian regime the transport current can be matched with that obtained by the non-equilibrium Green's function method.
We derive a non-Markovian master equation for the evolution of a class of open quantum systems consisting of quadratic fermionic systems coupled to wideband reservoirs. This is done by providing an explicit correspondence between master equations and nonequilibrium Green's function approaches. Our findings permit us to address non-Markovian regimes characterized by negative decoherence rates and to characterize the dynamics with respect to a recently proposed measurement of ``non-Markovianity.'' We study the real-time dynamics and the steady-state solution of two illustrative models: a tight-binding electronic model and XY spin chains. The rich set of nonequilibrium steady-state phases encountered extends previous results to the non-Markovian regime.
Open quantum systems have drawn attention over decades due to its applicability in the foundation of theoretical physics, e.g. statistical mechanics, quantum mechanics and condensed matter physics. The dynamics of open quantum systems has been described as separate entities from their surrounding environment that consists of a very large number of modes, somehow coupled to the mode of the system. Even though the exact solution of the dynamical behavior of the system is impossible to calculate, we obtain a tentative solution using the crucial Markov approximation. The input-output formalism of the quantum Langevin equation (QLE) has been considered as a useful tool which provides a semi-classical description of the dynamics of the system, whereas the master equation provides a complete picture of the dynamics of the system expressed in terms of density matrix. While studying the dynamics of nonlinear system/environment coupling using QLE, we see, for a small value of external field, that the steady state system field does not change much from the steady state field obtained in the absence of nonlinear dissipation. However, in a case where a stronger external field is applied, we see that the deviation becomes substantial from the solution of linear system. We also see that the nonlinear coupling introduces significant difference in the cavity fluctuation spectrum. The description, therefore, provides a potential explanation of parametric effects in terms of nonlinear dissipation phenomena associated with the nonlinear coupling. Even though the theories developed in the context of open quantum systems have proven to be powerful tools, they do not provide a satisfactory platform to be implemented on non-linear Hamiltonians. We often approximate it by linearizing over nonlinear steady state field amplitude, and therefore, the interesting effects are often overlooked. The limitation of the analytics provokes us to simulate open quantum dynamics numerically. The numerical method consists of transformation of the environmental degrees of freedom to a one-dimensional many-body chain, and the computational technique includes numerical diagonalization and renormalization process. The time-adaptive density matrix renormalisation group (t-DMRG) is known as one of the most powerful techniques for the simulation of strongly-correlated many-body quantum systems. In this thesis, along with the theoretical modeling, we implement DMRG numerical scheme for the simulation of canonical S/B model by mapping it to one-dimensional harmonic chain with nearest neighbor interactions, and use the method to investigate the dynamics of the free dissipative system. The thermalization of open quantum systems is also studied by generating minimally entangled typical thermal states (METTS) through imaginary time evolution, and real-time evolving an empty system in the presence of the thermal bath. Further, we simulate coherently driven free dissipative Kerr nonlinear system numerically using Euler's method by solving Heisenberg equation of motion and t-DMRG algorithm, and demonstrate how the numerical results are analogous to classical bistability. By comparing with analytics, we see that the DMRG numerics is analogous to the quantum-mechanical exact solution obtained by mapping the equation of motion of the density matrix of the system to a Fokker-Plank equation. The comparison between two different numerical techniques shows that the semi-classical Euler's method determines the dynamics of the system field of one among two coherent branches, whereas DMRG numerics gives the superposition of both of them. Hence, DMRG-determined time dynamics undergoes generating non-classical states. Our approach of dealing with nonlinearity represents an important contribution in the developments of technique to study the dynamical and steady-state behavior of open quantum systems, which is a fundamental aspect of quantum physics.