Full-dimensional control of the radiationless decay in pyrazine using the dynamic Stark effect

A multi-layer multi-configurational time-dependent Hartree approach to lattice models beyond one dimension

The multi-layer multi-configurational time-dependent Hartree (MCTDH) approach is an efficient method to study quantum dynamics in real and imaginary time. The present work explores its potential to describe quantum fluids. The multi-layer MCTDH approach in second quantization representation is used to study lattice models beyond one dimension at finite temperatures. A scheme to map the lattice sites onto the MCTDH tree representation for multi-dimensional lattice models is proposed. A statistical sampling scheme previously used in MCTDH calculations is adapted to facilitate an efficient description of the thermal ensemble. As example, a two-dimensional hard-core Bose-Hubbard model is studied considering up to 64 × 64 lattice sites. The single particle function basis set size required to obtain converged results is found to not increase with the lattice size. The numerical results properly simulate the finite temperature Berezinskii-Kosterlitz-Thouless phase transition.

Light-induced photodissociation in the lowest three electronic states of the NaH molecule

It has been known that electronic conical intersections in a molecular system can also be created by laser light even in diatomics. The direct consequence of these light-induced degeneracies is the appearance of a strong mixing between the electronic and vibrational motions, which has a strong fingerprint on the ultrafast nuclear dynamics. In the present work, pump and probe numerical simulations are performed with the NaH molecule involving the first three singlet electronic states (X1Σ+(X), A1Σ+(A) and B1Π(B)) and several light-induced degeneracies in the numerical description. To demonstrate the impact of the multiple light-induced non-adiabatic effects together with the molecular rotation on the dynamical properties of the molecule, the dissociation probabilities, kinetic energy release spectra (KER) and the angular distributions of the photofragments were calculated by discussing the role of the permanent dipole moment as well.

Eigenstate calculation in the state-averaged (multi-layer) multi-configurational time-dependent Hartree approach

A new approach for the calculation of eigenstates with the state-averaged (multi-layer) multi-configurational time-dependent Hartree (MCTDH) approach is presented. The approach is inspired by the recent work of Larsson [J. Chem. Phys. 151, 204102 (2019)]. It employs local optimization of the basis sets at each node of the multi-layer MCTDH tree and successive downward and upward sweeps to obtain a globally converged result. At the top node, the Hamiltonian represented in the basis of the single-particle functions (SPFs) of the first layer is diagonalized. Here p wavefunctions corresponding to the p lowest eigenvalues are computed by a block Lanczos approach. At all other nodes, a non-linear operator consisting of the respective mean-field Hamiltonian matrix and a projector onto the space spanned by the respective SPFs is considered. Here, the eigenstate corresponding to the lowest eigenvalue is computed using a short iterative Lanczos scheme. Two different examples are studied to illustrate the new approach: the calculation of the vibrational states of methyl and acetonitrile. The calculations for methyl employ the single-layer MCTDH approach, a general potential energy surface, and the correlation discrete variable representation. A five-layer MCTDH representation and a sum of product-type Hamiltonian are used in the acetonitrile calculations. Very fast convergence and order of magnitude reductions in the numerical effort compared to the previously used block relaxation scheme are found. Furthermore, a detailed comparison with the results of Avila and Carrington [J. Chem. Phys. 134, 054126 (2011)] for acetonitrile highlights the potential problems of convergence tests for high-dimensional systems.

A non-hierarchical correlation discrete variable representation

The correlation discrete variable representation (CDVR) facilitates (multi-layer) multi-configurational time-dependent Hartree (MCTDH) calculations with general potentials. It employs a layered grid representation to efficiently evaluate all potential matrix elements appearing in the MCTDH equations of motion. The original CDVR approach and its multi-layer extension show a hierarchical structure: the size of the grids employed at the different layers increases when moving from an upper layer to a lower one. In this work, a non-hierarchical CDVR approach, which uses identically structured quadratures at all layers of the MCTDH wavefunction representation, is introduced. The non-hierarchical CDVR approach crucially reduces the number of grid points required, compared to the hierarchical CDVR, shows superior scaling properties, and yields identical results for all three representations showing the same topology. Numerical tests studying the photodissociation of NOCl and the vibrational states of CH₃ demonstrate the accuracy of the non-hierarchical CDVR approach.

Identification of an ultrafast internal conversion pathway of pyrazine by time-resolved vacuum ultraviolet photoelectron spectrum simulations

The internal conversion from the optically bright S₂ (¹B_2u, ππ*) state to the dark S₁ (¹B_3u, nπ*) state in pyrazine is a standard benchmark for experimental and theoretical studies on ultrafast radiationless decay. Since 2008, a few theoretical groups have suggested significant contributions of other dark states S₃ (¹A_u, nπ*) and S₄ (¹B_2g, nπ*) to the decay of S₂. We have previously reported the results of nuclear wave packet simulations [Kanno et al., Phys. Chem. Chem. Phys. 17, 2012 (2015)] and photoelectron spectrum calculations [Mignolet et al., Chem. Phys. 515, 704 (2018)] that support the conventional two-state picture. In this article, the two different approaches, i.e., wave packet simulation and photoelectron spectrum calculation, are combined: We computed the time-resolved vacuum ultraviolet photoelectron spectrum and photoelectron angular distribution for the ionization of the wave packet transferred from S₂ to S₁. The present results reproduce almost all the characteristic features of the corresponding experimental time-resolved spectrum [Horio et al., J. Chem. Phys. 145, 044306 (2016)], such as a rapid change from a three-band to two-band structure. This further supports the existence and character of the widely accepted pathway (S₂ → S₁) of ultrafast internal conversion in pyrazine.

Symmetries in the multi-configurational time-dependent Hartree wavefunction representation and propagation

In multi-configurational time-dependent Hartree (MCTDH) approaches, different multi-layered wavefunction representations can be used to represent the same physical wavefunction. Transformations between different equivalent representations of a physical wavefunction that alter the tree structure used in the multi-layer MCTDH wavefunction representation interchange the role of single-particle functions (SPFs) and single-hole functions (SHFs) in the MCTDH formalism. While the physical wavefunction is invariant under these transformations, this invariance does not hold for the standard multi-layer MCTDH equations of motion. Introducing transformed SPFs, which obey normalization conditions typically associated with SHFs, revised equations of motion are derived. These equations do not show the singularities resulting from the inverse single-particle density matrix and are invariant under tree transformations. Based on the revised equations of motion, a new integration scheme is introduced. The scheme combines the advantages of the constant mean-field approach of Beck and Meyer [Z. Phys. D 42, 113 (1997)] and the singularity-free integrator suggested by Lubich [Appl. Math. Res. Express 2015, 311]. Numerical calculations studying the spin boson model in high dimensionality confirm the favorable properties of the new integration scheme.

Application of optimal control simulation to selective photodissociation of IBr by non-resonant dynamic Stark effects

We apply nonlinear optimal control simulation to design a non-resonant control pulse that maximizes the probability of specified photodissociation of IBr by utilizing the non-resonant dynamic Stark effect in the presence of a predetermined pump pulse. The optimal pulses are always composed of several subpulses that increase the target probability considerably depending on the wavelength of the pump pulse. Focusing on the cases of high target probabilities, we systematically examine how the subpulses cooperate with each other on the basis of pulse-partitioning analyses. We show that the subpulses largely cooperate with the pump pulse, which can explain their irradiation timings. On the other hand, the cooperation between the subpulses is mainly expressed as the sum of the contribution from each subpulse.

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.

Quantum dynamics of the photostability of pyrazine

We propose a new mechanism for the radiationless decay of photoexcited pyrazine to its ground electronic state involving a conical intersection between the dark A_u(nπ) state and the ground state.