Quantum transport in molecular junctions has attracted great attention. The charge motion in a molecular junction can cause geometric deformation, leading to strong electron phonon coupling, which was often overlooked. We have formulated a nearly exact method to assess the time-dependent current and occupation number in the molecular junction modeled by the electron-phonon coupled bridge state using the time-dependent density matrix renormalization group (TD-DMRG) method. The oscillation period and amplitude of the current are found to be dependent on the electron phonon coupling strength and energy level alignment with the electrodes. In an attempt to better understand these phenomena, we have devised a new approximation that explains the bistability phenomenon and the behavior of steady currents in the strong electron-phonon coupling regime. Comparisons have been made with the multilayer-multiconfiguration time-dependent Hartree (ML-MCTDH) method and the analytical result in the purely electronic limit. Furthermore, we explore the entropy of different orderings, extending to the electron phonon model problems. Regarding finite temperature, the thermal Bogoliubov transformation of both fermions and bosons is used and compared with imaginary time evolution results.
Quantum transport in molecular junctions has attracted great attention. The charge motion in a molecular junction can cause geometric deformation, leading to strong electron phonon coupling, which was often overlooked. We have formulated a nearly exact method to assess the time-dependent current and occupation number in the molecular junction modeled by the electron-phonon coupled bridge state using the time-dependent density matrix renormalization group (TD-DMRG) method. The oscillation period and amplitude of the current are found to be dependent on the electron phonon coupling strength and energy level alignment with the electrodes. In an attempt to better understand these phenomena, we have devised a new approximation that explains the bistability phenomenon and the behavior of steady currents in the strong electron-phonon coupling regime. Comparisons have been made with the multilayer-multiconfiguration time-dependent Hartree (ML-MCTDH) method and the analytical result in the purely electronic limit. Furthermore, we explore the entropy of different orderings, extending to the electron phonon model problems. Regarding finite temperature, the thermal Bogoliubov transformation of both fermions and bosons is used and compared with imaginary time evolution results.
Recent advances in experimental techniques have enabled precise characterization of fundamental nonadiabatic and even electronic dynamics in molecules, on ultrafast time scales reaching femtosecond to attosecond resolution. However, accurate theoretical simulation of these ultrafast chemical dynamics processes in large systems remains challenging, largely due to the overwhelming number of degrees of freedom (DoFs) and the pronounced many-body correlations. In recent years, by leveraging efficient decomposition schemes for high-dimensional wave function and operator tensors, the time-dependent density matrix renormalization group (TD-DMRG) has emerged as a powerful and accurate quantum dynamics method for simulating nonadiabatic and electron dynamics in large chemical systems, in conjunction with realistic electron/exciton-vibration/phonon models or ab initio quantum chemistry many-electron Hamiltonians. This review outlines the fundamentals of TD-DMRG for chemical dynamics, covering matrix product state/operator (MPS/MPO) frameworks and algorithms from ground-state calculations to time evolution. We discuss thermal/environmental effects and compare TD-DMRG with other tensor network methods such as multiconfiguration time-dependent Hartree (MCTDH) and multilayer MCTDH (ML-MCTDH). Demonstrated applications include simulations of pyrazine absorption, singlet fission in rubrene crystal, and charge migration in chloroacetylene cation. These show TD-DMRG’s capability for modeling complex ultrafast processes from femtoseconds to attoseconds with controlled accuracy.
The problem of quantum dynamics in open systems has gained attention in recent decades and not the least due to the advances made in quantum transport in molecular systems. The main motivation behind quantum transport and molecular electronics is the futuristic goal to be able at some point to replace, or to complement, the silicon-based technology and to make the electronic devices faster. On a fundamental level, one has to deal with time-dependent processes where electron-electron or electron- phonon interactions are of great importance, and they can cause profound quantitative and qualitative changes on the physical and dynamical properties of electronic systems compared to the non-interacting case. Most of the studies of quantum transport have been focused on the steady-state description while neglecting the short-time dynamics. However, the dynamical effects are of great importance since fast- switching processes play a pivotal role in the operation of future devices. We studied the problem of time-dependent electron transport through the Anderson impurity model by using many-body perturbation theory (MBPT) together with Keldysh Green’s functions as well as with time-dependent density functional theory (TDDFT). These methods were compared with numerically exact time-dependent density-matrix renormalization group (tDMRG) method. We found that the many-body perturbation theory results beyond Hartree-Fock approximation were in close agreement with tDMRG results. In addition we studied the possibility of multistablity in the density and current of an interacting nanoscale junction as well as how to reversibly switch between the multiple solutions in time domain. An accurate theoretical treatment of electron correlation even in as simple model as an interacting electron gas at metallic densities still continues to be a challenge; especially description of features in the photoemission spectra due to electron correlations provides a theoretical challenge. The many-body perturbation theory yields a systematic way to study electron-electron (electron-phonon) correlations in various systems. One of the widely used approximations in MBPT is the GW approximation in which the bare interaction line is replaced with screened interaction line in the first order exchange diagram. The GW approximation gives good estimates for the band gap values close to experimental ones but especially the self-consistent GW approximation has a number of deficiencies like washing out of plasmon features and overestimation of bandwidths compared to experiment. One way to improve GW calculations is to include vertex corrections. Unfortunately, the straightforward inclusion of vertex corrections yields negative spectra in some frequency regions. We developed a diagrammatic approach to construct self-energy approximations with positive spectral properties. Our approach consists of expressing a self-energy of response diagram as a product of half-diagrams after which a minimal set of additional diagrams is identified to construct a perfect square. We applied this method to study vertex corrections in a homogeneous electron gas. In addition we analyzed the diagrammatic content of photocurrent with density functional theory. The expression for the photocurrent was obtained as an integral over the Kohn-Sham spectral function renormalized by effective potentials that depend on the exchange correlation kernel of current density functional theory. The expression for the photocurrent gives us the angular dependence of the photocurrent but it does not provide a direct access to the kinetic energy distribution of the photoelectrons.
Organic/polymeric materials are of emerging importance for thermoelectric conversion. The soft nature of these materials implies strong electron-phonon coupling, often leading to carrier localization. This poses great challenges for the conventional Boltzmann transport description based on relaxation time approximation and band structure calculations. In this work, combining Kubo formula with finite-temperature time-dependent density matrix renormalization group (FT-TD-DMRG) in the grand canonical ensemble, we developed a nearly exact algorithm to calculate the thermoelectric power factor PF=α^2 σ, where α is the Seebeck coefficient and σ is electrical conductivity, and apply the algorithm to organic materials modeled by Holstein Hamiltonian with electron-phonon coupling. Our algorithm can provide a unified description covering the weak coupling limit described by bandlike Boltzmann transport to the strong coupling hopping limit.
In this work we develop a new numerical method, the adaptive time-dependent density-matrix renormalization-group (adaptive t-DMRG), which turns out to be very well suited to investigate time-dependent phenomena in one-dimensional strongly correlated systems. The method is based on the original density-matrix renormalization group method by White and the time evolving block-decimation procedure by Vidal. We further show the applicability of the new method to different physical systems: bosonic, fermionic, and spin chains. We first discuss the consequences of the presence of an external trapping potential on the static properties of trapped cold bosonic atoms in an optical lattice using the Bose-Hubbard model. It is shown that a properly rescaled one-particle density matrix characterizes superfluid versus insulating states just as in the homogeneous system. We note that the superfluid to Mott-insulating transition is seen most directly in the half width of the interference peak. Then we turn to the evolution of density perturbations in these systems. In particular, we discuss the dependence of the velocity and the decay of the amplitude on density, interaction strength, and the extent and height of the perturbation in a numerically exact way, covering a wide range of interaction and perturbation strengths. By comparing our results for the sound velocity to theoretical predictions, we determine the limits of a Gross-Pitaevskii or Bogoliubov type description and the regime where repulsive one-dimensional Bose gases exhibit fermionic behaviour. As a second application of the adaptive t-DMRG, we investigate the phenomenon of spin-charge separation in the Hubbard model and propose an experimental setup for its observation in cold Fermi gases. We show the robustness of this separation beyond the low-energy regime by studying the time evolution of density wave packets of finite strength and at length scales down to a few lattice spacings. A striking signature of spin-charge separation is found in 1D cold Fermi gases in a harmonic trap using the different propagation properties of the liquid and Mott-insulating phases. The third application relates to the area of spintronics. A simplified model for the magnetization transport between two coupled reservoirs with opposite spin polarization is studied. We investigate whether a simple long time limit for the spin transport exists and if so, how it depends on the properties of the system such as integrability and criticality. Time-scales accessible to us are of the order of 100 units of time measured in 1/J while maintaining insignificant error in the observables.
Electronic and/or vibronic coherence has been found by recent ultrafast spectroscopy experiments in many chemical, biological, and material systems. This indicates that there are strong and complicated interactions between electronic states and vibration modes in realistic chemical systems. Therefore, simulations of quantum dynamics with a large number of electronic and vibrational degrees of freedom are highly desirable. Due to the efficient compression and localized representation of quantum states in the matrix-product state (MPS) formulation, time-evolution methods based on the MPS framework, which we summarily refer to as tDMRG (time-dependent density-matrix renormalization group) methods, are considered to be promising candidates to study the quantum dynamics of realistic chemical systems. In this work, we benchmark the performances of four different tDMRG methods, including global Taylor, global Krylov, and local one-site and two-site time-dependent variational principles (1TDVP and 2TDVP), with a comparison to multiconfiguration time-dependent Hartree and experimental results. Two typical chemical systems of internal conversion and singlet fission are investigated: one containing strong and high-order local and nonlocal electron-vibration couplings and the other exhibiting a continuous phonon bath. The comparison shows that the tDMRG methods (particularly, the 2TDVP method) can describe the full quantum dynamics in large chemical systems accurately and efficiently. Several key parameters in the tDMRG calculation including the truncation error threshold, time interval, and ordering of local sites were also investigated to strike the balance between efficiency and accuracy of results.
A procedure based on the recently developed ``adaptive'' time-dependent density-matrix-renormalization-group (DMRG) technique is presented to calculate the zero temperature conductance of nanostructures, such as quantum dots (QDs) or molecular conductors, when represented by a small number of active levels. The leads are modeled using noninteracting tight-binding Hamiltonians. The ground state at time zero is calculated at zero bias. Then, a small bias is applied between the two leads, the wave function is DMRG evolved in time, and currents are measured as a function of time. Typically, the current is expected to present periodicities over long times, involving intermediate well-defined plateaus that resemble steady states. The conductance can be obtained from those steady-state-like currents. To test this approach, several cases of interacting and noninteracting systems have been studied. Our results show excellent agreement with exact results in the noninteracting case. More importantly, the technique also reproduces quantitatively well-established results for the conductance and local density of states in both the cases of one and two coupled interacting QDs. The technique also works at finite bias voltages, and it can be extended to include interactions in the leads.