Linear and nonlinear optical response of polyenes: A density matrix renormalization group study

Optimal tree tensor network operators for tensor network simulations: Applications to open quantum systems

Tree tensor network states (TTNS) decompose the system wavefunction to the product of low-rank tensors based on the tree topology, serving as the foundation of the multi-layer multi-configuration time-dependent Hartree method. In this work, we present an algorithm that automatically constructs the optimal and exact tree tensor network operators (TTNO) for any sum-of-product symbolic quantum operator. The construction is based on the minimum vertex cover of a bipartite graph. With the optimal TTNO, we simulate open quantum systems, such as spin relaxation dynamics in the spin-boson model and charge transport in molecular junctions. In these simulations, the environment is treated as discrete modes and its wavefunction is evolved on equal footing with the system. We employ the Cole-Davidson spectral density to model the glassy phonon environment and incorporate temperature effects via thermo-field dynamics. Our results show that the computational cost scales linearly with the number of discretized modes, demonstrating the efficiency of our approach.

On the fly swapping algorithm for ordering of degrees of freedom in density matrix renormalization group

Density matrix renormalization group (DMRG) and its time-dependent variants have found widespread applications in quantum chemistry, includingab initioelectronic structure of complex bio-molecules, spectroscopy for molecular aggregates, and charge transport in bulk organic semiconductors. The underlying wavefunction ansatz for DMRG, matrix product state (MPS), requires mapping degrees of freedom (DOF) into a one-dimensional topology. DOF ordering becomes a crucial factor for DMRG accuracy. In this work, we propose swapping neighboring DOFs during the DMRG sweeps for DOF ordering, which we term 'on the fly swapping' (OFS) algorithm. We show that OFS is universal for both static and time-dependent DMRG with minimum computational overhead. Examples are given for one dimensional antiferromagnetic Heisenberg model,ab initioelectronic structure of N₂molecule, and the S₁/S₂internal conversion dynamics of pyrazine molecule. It is found that OFS can indeed improve accuracy by finding better DOF ordering in all cases.

Finite Temperature Dynamical Density Matrix Renormalization Group for Spectroscopy in Frequency Domain

We present a novel methodology through casting the dynamical density matrix renormalization group (DDMRG) into the matrix product state (MPS) formulation to calculate the spectroscopy at finite temperature for molecular aggregates. The frequency domain algorithm can avoid the time evolution accumulation of error and is naturally suitable for parallelization, in addition to facile graphic processing unit (GPU) acceleration. The high accuracy is demonstrated by simulating the optical spectra of vibronic model systems ranging from an exactly solvable dimer model to a more complex real-world perylene bisimide (PBI) J-aggregate. The relationship between the 0-0 emission strength and the exciton thermal coherent length is discussed for linearly stacked aggregates. The computing performance largely boosted by GPU demonstrates that DDMRG emerges as a promising method to study dynamical properties for complex systems.

Numerical assessment for accuracy and GPU acceleration of TD-DMRG time evolution schemes

The time dependent density matrix renormalization group (TD-DMRG) has become one of the cutting edge methods of quantum dynamics for complex systems. In this paper, we comparatively study the accuracy of three time evolution schemes in the TD-DMRG, the global propagation and compression method with the Runge-Kutta algorithm (P&C-RK), the time dependent variational principle based methods with the matrix unfolding algorithm (TDVP-MU), and with the projector-splitting algorithm (TDVP-PS), by performing benchmarks on the exciton dynamics of the Fenna-Matthews-Olson complex. We show that TDVP-MU and TDVP-PS yield the same result when the time step size is converged and they are more accurate than P&C-RK4, while TDVP-PS tolerates a larger time step size than TDVP-MU. We further adopt the graphical processing units to accelerate the heavy tensor contractions in the TD-DMRG, and it is able to speed up the TDVP-MU and TDVP-PS schemes by up to 73 times.

Organic light-emitting diodes: theoretical understanding of highly efficient materials and development of computational methodology

Theoretical understanding of organic light-emitting diodes started from the quest to the nature of the primary excitation in organic molecular and polymeric materials. We found the electron correlation strength, bond-length alternation as well as the conjugation extent have strong influences on the orderings of the lowest lying excited states through the first application of density matrix renormalization group theory to quantum chemistry. The electro-injected free carriers (with spin 1/2) can form both singlet and triplet bound states. We found that the singlet exciton formation ratio can exceed the conventional 25% spin statistics limit. We proposed a vibration correlation function formalism to evaluate the excited-state decay rates, which is shown to not only give reasonable estimations for the quantum efficiency but also a quantitative account for the aggregation-induced emission (AIE). It is suggested to unravel the AIE mechanism through resonance Raman spectroscopy.

HEISENBERG MODEL AND ITS APPLICATIONS TO π-CONJUGATED SYSTEMS

This paper briefly reviewed the Heisenberg model and its improvement by including higher order corrections, and their applications to bond lengths, stability, and reactivity of non-benzenoids and the low-lying excitation spectra of conjugated systems. Two efficient computational methods, the Lanczos method and the density-matrix renormalization group (DMRG), for exactly and approximately solving various Heisenberg models, respectively, were briefly introduced.

COUPLED-CLUSTER EQUATION OF MOTION STUDY FOR THE ELECTRONIC AND OPTICAL PROPERTIES OF CONJUGATED SYSTEMS

We have implemented the Coupled-Cluster Equation of Motion (CC-EOM) with single and double excitations coupled with semiempirical parameterizations, Pariser–Parr–Pople model and INDO Hamiltonians, to investigate the optical and nonlinear optical properties, electronic structures and the excited states properties for the conjugated polymers. The semiempirical parameters allow us to study the conjugated systems with extensive sizes. Firstly, by comparing with the quasi-exact Density Matrix Renormalization Group theory for the 1D conjugated chain, we find that the CC-EOM approach can give satisfactory results for both the ground state and the excited states energies. We demonstrate that our approach can be adopted to evaluate linear and nonlinear response. We find that both the real and imaginary parts of the third-order polarizability can be solved either for static and dynamic responses. Then we apply the CC-EOM approach to study the optical signatures for the polarons in conjugated polymers. We have established a solid relationship between the rigidity of a polymer and its optical signature of the polarons.

Excitation energy calculation of conjugated hydrocarbons: a new Pariser-Parr-Pople model parameterization approaching CASPT2 accuracy

A new parameterization for the Pariser-Parr-Pople (PPP) model for conjugated hydrocarbons is proposed in this work. The distance-dependence of PPP parameters are obtained from CASPT2 ground state and low-lying excited state energies of ethylene and its cation at various C-C single bond lengths and are fitted to a set of carefully chosen mathematical functions. Our new PPP model is applied to the calculation of vertical singlet-triplet energy gaps and the excitation energies for low-lying π→π(*) valence excitations in various π-conjugated molecules. Results with the new PPP model are consistently better than the standard PPP model in use. It often surpasses density functional theory and single-reference excited state methods such as configuration interaction singles or time-dependent density functional theory in terms of its accuracy and agrees reasonably well with high-level theories or experiments.

Dynamical simulations of charged soliton transport in conjugated polymers with the inclusion of electron-electron interactions

We present numerical studies of the transport dynamics of a charged soliton in conjugated polymers under the influence of an external time-dependent electric field. All relevant electron-phonon and electron-electron interactions are nearly fully taken into account by simulating the monomer displacements with classical molecular dynamics and evolving the wave function for the pi electrons by virtue of the adaptive time-dependent density matrix renormalization group simultaneously and nonadiabatically. It is found that after a smooth turn on of the external electric field the charged soliton is accelerated at first up to a stationary constant velocity as one entity consisting of both the charge and the lattice deformation. An Ohmic region (6 mV/A</=E(0)</=12 mV/A) where the stationary velocity increases linearly with the electric field strength is observed. The relationship between electron-electron interactions and charged soliton transport is also investigated in detail. We find that the dependence of the stationary velocity of a charged soliton on the on-site Coulomb interactions U and the nearest-neighbor interactions V is due to the extent of delocalization of the charged soliton defect.