In this work, we propose a new method to calculate molecular nonradiative electronic relaxation rates based on the numerically exact time-dependent density matrix renormalization group theory. This method could go beyond the existing frameworks under the harmonic approximation (HA) of the potential energy surface (PES) so that the anharmonic effect could be considered, which is of vital importance when the electronic energy gap is much larger than the vibrational frequency. We calculate the internal conversion (IC) rates in a two-mode model with Morse potential to investigate the validity of HA. We find that HA is unsatisfactory unless only the lowest several vibrational states of the lower electronic state are involved in the transition process when the adiabatic excitation energy is relatively low. As the excitation energy increases, HA first underestimates and then overestimates the IC rates when the excited state PES shifts toward the dissociative side of the ground state PES. On the contrary, HA slightly overestimates the IC rates when the excited state PES shifts toward the repulsive side. In both cases, a higher temperature enlarges the error of HA. As a real example to demonstrate the effectiveness and scalability of the method, we calculate the IC rates of azulene from S1 to S0 on the ab initio anharmonic PES approximated by the one-mode representation. The calculated IC rates of azulene under HA are consistent with the analytically exact results. The rates on the anharmonic PES are 30%–40% higher than the rates under HA.

Using a photonic quantum computer for boson sampling has demonstrated a tremendous advantage over classical supercomputers. It is highly desirable to develop boson sampling algorithms for realistic scientific problems. In this work, we propose a hybrid quantum-classical sampling (HQCS) algorithm to calculate the optical spectrum for complex molecules considering Duschinsky rotation effects and anharmonicity. The classical sum-over-states method for this problem has a computational complexity that exponentially increases with system size. The HQCS algorithm creates an intermediate harmonic potential energy surface (PES) to bridge the initial and final PESs. The magnitude and sign of the overlap between the initial and the intermediate state are estimated by boson sampling and classical algorithms, respectively. The overlap between the intermediate and the final state is efficiently evaluated by classical algorithms. The feasibility of HQCS is demonstrated in calculations of the emission spectrum of a Morse model as well as the pyridine molecule.

The widely known ``Energy Gap Law'' (EGL) predicts a monotonically exponential increase in the non-radiative decay rate (knr) as the energy gap narrows, which hinders the development of near-infrared (NIR) emissive molecular materials. Recently, several experiments proposed that the exciton delocalization in molecular aggregates could counteract EGL to facilitate NIR emission. In this work, the nearly exact time-dependent density matrix renormalization group (TD-DMRG) method is developed to evaluate the non-radiative decay rate for exciton-phonon coupled molecular aggregates. Systematical numerical simulations show, by increasing the excitonic coupling, knr will first decrease, then reach a minimum, and finally start to increase to follow the EGL, which is an overall result of two opposite effects of a smaller energy gap and a smaller effective electron-phonon coupling. This anomalous non-monotonic behavior is found robust in a number of models, including dimer, one-dimensional chain, and two-dimensional square lattice. The optimal excitonic coupling strength that gives the minimum knr is about half of the monomer reorganization energy and is also influenced by temperature, system size, and dimensionality.