Photodissociation of methyl iodide embedded in a host-guest complex: A full dimensional (189D) quantum dynamics study of CH3I@resorc[4]arene

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.

Quantum rates in dissipative systems with spatially varying friction

We investigate whether making the friction spatially dependent on the reaction coordinate introduces quantum effects into the thermal reaction rates for dissipative reactions. Quantum rates are calculated using the numerically exact multi-configuration time-dependent Hartree method, as well as the approximate ring-polymer molecular dynamics (RPMD), ring-polymer instanton methods, and classical molecular dynamics. By conducting simulations across a wide range of temperatures and friction strengths, we can identify the various regimes that govern the reactive dynamics. At high temperatures, in addition to the spatial-diffusion and energy-diffusion regimes predicted by Kramer's rate theory, a (coherent) tunneling-dominated regime is identified at low friction. At low temperatures, incoherent tunneling dominates most of Kramer's curve, except at very low friction, when coherent tunneling becomes dominant. Unlike in classical mechanics, the bath's influence changes the equilibrium time-independent properties of the system, leading to a complex interplay between spatially dependent friction and nuclear quantum effects even at high temperatures. More specifically, a realistic friction profile can lead to an increase (or decrease) of the quantum (classical) rates with friction within the spatial-diffusion regime, showing that classical and quantum rates display qualitatively different behaviors. Except at very low frictions, we find that RPMD captures most of the quantum effects in the thermal reaction rates.

On the multiphoton ionisation photoelectron spectra of phenol

The phenol molecule is a prototype for non-adiabatic dynamics and the excited-state photochemistry of biomolecules. In this article, we report a joint theoretical and experimental investigation on the resonance enhanced multiphoton ionisation photoelectron (REMPI) spectra of the two lowest ionisation bands of phenol. The focus is on the theoretical interpretation of the measured spectra using quantum dynamics simulations. These were performed by numerically solving the time-dependent Schrödinger equation using the multi-layer variant of the multiconfiguration time-dependent Hartree algorithm together with a vibronic coupling Hamiltonian model. The ionising laser pulse is modelled explicitly within the ionisation continuum model to simulate experimental femtosecond 1+1 REMPI photoelectron spectra. These measured spectra are sensitive to very short lived electronically excited states, providing a rigorous benchmark for our theoretical methods. The match between experiment and theory allows for an interpretation of the features of the spectra at different wavelengths and shows that there are features due to both 'direct' and 'indirect' ionisation, resulting from non-resonant and resonant excitation by the pump pulse.

Optimal Mode Combination in the Multiconfiguration Time-Dependent Hartree Method through Multivariate Statistics: Factor Analysis and Hierarchical Clustering

The multiconfiguration time-dependent Hartree (MCTDH) method and its multilayer extension (ML-MCTDH) are powerful algorithms for the efficient computation of nuclear quantum dynamics in high-dimensional systems. By providing time-dependent variational orbitals and an optimal choice of layered effective degrees of freedom, one is able to reduce the computational cost to an amenable number of configurations. However, choices related to selecting properly the mode grouping and tensor tree are strongly system dependent and, thus far, subjectively based on intuition and/or experience. Therefore, herein we detail a new protocol based on multivariate statistics─more specifically, factor analysis and hierarchical clustering─for a reliable and convenient guiding in the optimal design of such complex "system-of-systems" tensor-network decompositions. The advantages of employing the new algorithm and its applicability are tested on water and two floppy protonated water clusters with large amplitude motions.

Development of a fully coupled diabatic spin-orbit model for the photodissociation of phenyl iodide

The theoretical treatment of the quantum dynamics of the phenyl iodide photodissociation requires an accurate analytical potential energy surface (PES) model. This model must also account for spin-orbit (SO) coupling. This study is the first step to construct accurate SO coupled PESs, namely, for the C-I dissociation coordinate. The model is based on the Effective Relativistic Coupling by Asymptotic Representation (ERCAR) method developed over the past ten years. The SO-free Hamiltonian is represented in an asymptotic diabatic basis and then combined with an atomic effective relativistic coupling operator determined analytically. In contrast to the previously studied cases (HI, CH₃I), the diabatic basis states are due to excitations in the phenyl fragment rather than the iodine atom. An accurate analytical model of the ab initio reference data is determined in two steps. The first step is a simple reference model describing the data qualitatively. This reference model is corrected through a trained artificial neural-network to achieve high accuracy. The SO-free and the fine structure states resulting from this ERCAR model are discussed extensively in the context of the photodissociation.

Comparison of the multi-layer multi-configuration time-dependent Hartree (ML-MCTDH) method and the density matrix renormalization group (DMRG) for ground state properties of linear rotor chains

We demonstrate the applicability of the Multi-Layer Multi-Configuration Time-Dependent Hartree (ML-MCTDH) method to the problem of computing ground states of one-dimensional chains of linear rotors with dipolar interactions. Specifically, we successfully obtain energies, entanglement entropies, and orientational correlations that are in agreement with the Density Matrix Renormalization Group (DMRG), which has been previously used for this system. We find that the entropies calculated by ML-MCTDH for larger system sizes contain nonmonotonicity, as expected in the vicinity of a second-order quantum phase transition between ordered and disordered rotor states. We observe that this effect remains when all couplings besides nearest-neighbor are omitted from the Hamiltonian, which suggests that it is not sensitive to the rate of decay of the interactions. In contrast to DMRG, which is tailored to the one-dimensional case, ML-MCTDH (as implemented in the Heidelberg MCTDH package) requires more computational time and memory, although the requirements are still within reach of commodity hardware. The numerical convergence and computational demand of two practical implementations of ML-MCTDH and DMRG are presented in detail for various combinations of system parameters.

Formation of excited triplet states in naphthalene diimide and perylene diimide derivatives: A detailed theoretical study

Mechanistic details of the excited triplet state formation upon photoexcitation to the low-lying singlet manifold in naphthalene diimide and perylene diimide derivatives are explored theoretically. Static and dynamic aspects of two singlets (S₁ and S₂) and six triplets (T₁-T₆) of these molecules are investigated. Suitable vibronic Hamiltonians are constructed to investigate the internal conversion dynamics in both the singlet and triplet manifolds. Computed singlet-triplet energetics, spin-orbit coupling matrix elements, and intersystem crossing rates strongly suggest an efficient intersystem crossing process involving higher triplet states (T₆, T₅, and T₄). Separate full dimensional quantum wavepacket simulations of singlet and triplet manifolds in the approximate linear vibronic model by assuming initial Franck-Condon conditions are carried out to unravel the internal conversion decay dynamics in the respective manifolds. The obtained diabatic electronic populations and nuclear densities are analyzed to illustrate the triplet generation pathways involving higher triplet states in these molecules.

Excited-state intramolecular proton transfer driven by conical intersection in hydroxychromones

Nonradiative decay pathways associated with vibronically coupled S₁ (ππ*)-S₂ (nπ*) potential energy surfaces of 3- and 5-hydroxychromones are investigated by employing the linear vibronic coupling approach. The presence of a conical intersection close to the Franck-Condon point is identified based on the critical examination of computed energetics and structural parameters of stationary points. We show that very minimal displacements of relevant atoms of intramolecular proton transfer geometry are adequate to drive the molecule toward the conical intersection nuclear configuration. The evolving wavepacket on S₁ (ππ*) bifurcates at the conical intersection: a part of the wavepacket moves to S₂ (nπ*) within a few femtoseconds while the other decays to S₁ minimum. Our findings indicate the possibility of forming the proton transfer tautomer product via S₂ (nπ*), competing with the traditional pathway occurring on S₁ (ππ*).

Probing photodissociation dynamics using ring polymer molecular dynamics

The performance of the ring polymer molecular dynamics (RPMD) approach to simulate typical photodissociation processes is assessed. The correct description of photodissociation requires the calculation of correlation functions or expectation values associated with non-equilibrium initial conditions, which was shown to be possible with RPMD very recently [J. Chem. Phys. 145, 204118 (2016)]. This approach is combined with treatment of the nonadiabatic dynamics employing the ring polymer surface hopping approach (RPSH), which is based on Tully's fewest switches surface hopping (FSSH) approach. The performance is tested using one-dimensional photodissociation models. It is found that RPSH with non-equilibrium initial conditions can well reproduce the time-dependent dissociation probability, and adiabatic and diabatic populations for cases where the crossing point is below and above the Franck-Condon point, respectively, while standard FSSH fails to reproduce the exact quantum dynamics in the latter case. Thus, it is shown that RPSH is an efficient and accurate alternative to standard FSSH, which is one of the most widely employed approaches to study photochemistry.

Fast and slow excited-state intramolecular proton transfer in 3-hydroxychromone: a two-state story?

The photodynamics of 3-hydroxychromone in its first-excited singlet electronic state (bright state of ππ* character) is investigated with special emphasis given to two types of reaction pathways: the excited-state intramolecular-proton-transfer coordinate and the hydrogen-torsion coordinate linking the excited cis and trans isomers. A newly-found conical intersection with the second-excited singlet electronic state (dark state of nπ* character) is suspected to be, to some extent, the reason for the slower rate constant. This hypothesis based on quantum-chemistry calculations is supported by quantum-dynamics simulations in full dimensionality. They show significant transfer of electronic population and provide consistently a vibronic interpretation for the forbidden band in the UV absorption spectrum.

Wavepacket dynamics and the multi-configurational time-dependent Hartree approach

Multi-configurational time-dependent Hartree (MCTDH) based approaches are efficient, accurate, and versatile methods for high-dimensional quantum dynamics simulations. Applications range from detailed investigations of polyatomic reaction processes in the gas phase to high-dimensional simulations studying the dynamics of condensed phase systems described by typical solid state physics model Hamiltonians. The present article presents an overview of the different areas of application and provides a comprehensive review of the underlying theory. The concepts and guiding ideas underlying the MCTDH approach and its multi-mode and multi-layer extensions are discussed in detail. The general structure of the equations of motion is highlighted. The representation of the Hamiltonian and the correlated discrete variable representation (CDVR), which provides an efficient multi-dimensional quadrature in MCTDH calculations, are discussed. Methods which facilitate the calculation of eigenstates, the evaluation of correlation functions, and the efficient representation of thermal ensembles in MCTDH calculations are described. Different schemes for the treatment of indistinguishable particles in MCTDH calculations and recent developments towards a unified multi-layer MCTDH theory for systems including bosons and fermions are discussed.

Multilayer Multiconfiguration Time-Dependent Hartree Theory

Multilayer multiconfiguration time-dependent Hartree (ML-MCTDH) theory is a rigorous and powerful method to simulate quantum dynamics in complex many-body systems. This approach extends the original MCTDH theory of Meyer, Manthe, and Cederbaum to include dynamically contracted layers in a recursive way, within which the equations of motion are determined from the Dirac-Frenkel variational principle. This paper presents the general derivation of the theory and analyzes the important features that make the ML-MCTDH method numerically efficient. Furthermore, we discuss the generalization of the theory to treat many-body identical particles (fermions or bosons) as well as calculating energy eigenstates via the improved relaxation method.

A multilayer MCTDH study on the full dimensional vibronic dynamics of naphthalene and anthracene cations

Employing the multilayer multiconfiguration time-dependent Hartree (ML-MCTDH) method in conjunction with the multistate multimode vibronic coupling Hamiltonian model, we perform a full dimensional quantum dynamical study on the naphthalene (48D) and anthracene (66D) radical cations in their six lowest-lying doublet electronic states. For easily comparing results of full and reduced dimensionalities, MCTDH simulations based on larger sizes of primitive basis functions and single-particle functions than the previous ones [S. Ghanta, V. S. Reddy, and S. Mahapatra, Phys. Chem. Chem. Phys. 13, 14531 (2011)], are also performed. Extensive ML-MCTDH test calculations are performed to find appropriate ML separations of the wave functions (so-called ML-trees), and the convergence of the dynamical calculations are carefully checked. The ML-MCTDH method was developed for efficiently simulating quantum dynamics of large systems, and in fact the full dimensional ML-MCTDH calculations save a considerable amount of CPU-time in comparison with corresponding reduced dimensional MCTDH simulations. On basis of the present full and reduced dimensional simulations, the photoelectron (PE) spectra of these two cations are simulated and compared with corresponding experimental spectra. The agreement between theoretical and experimental PE spectra is good. Both full and reduced dimensional simulations give all main bands in the PE spectra. The vibronic energy-level positions from both ML-MCTDH and MCTDH calculations agree with corresponding experimental results. These quantum dynamical studies also complement the observations on diffuse interstellar bands with the wavelength of ~7088, ~6707, ~6490, ~6120, and ~5959 Å measured by astronomers as well as laboratory experimentalists.