Assessing the performance of trajectory surface hopping methods: Ultrafast internal conversion in pyrazine.

Trajectory surface hopping (TSH) methods have been widely used to study photoinduced nonadiabatic processes. In the present study, nonadiabatic dynamics simulations with the widely used Tully's fewest switches surface hopping (FSSH) algorithm and a Landau-Zener-type TSH (LZSH) algorithm have been performed for the internal conversion dynamics of pyrazine. The accuracy of the two TSH algorithms has been critically evaluated by a direct comparison with exact quantum dynamics calculations for a model of pyrazine. The model comprises the three lowest excited electronic states (B3u(nπ*), A1u(nπ*), and B2u(ππ*)) and the nine most relevant vibrational degrees of freedom. Considering photoexcitation to the diabatic B2u(ππ*) state, we examined the time-dependent diabatic and adiabatic electronic population dynamics. It is found that the diabatic populations obtained with both TSH methods are in good agreement with the exact quantum results. Fast population oscillations between the B3u(nπ*) and A1u(nπ*) states, which reflect nonadiabatic electronic transitions driven by coherent dynamics in the normal mode Q8a, are qualitatively reproduced by both TSH methods. In addition to the model study, the TSH methods have been interfaced with the second-order algebraic diagrammatic construction ab initio electronic-structure method to perform full-dimensional on-the-fly nonadiabatic dynamics simulations for pyrazine. It is found that the electronic population dynamics obtained with the LZSH method is in excellent agreement with that obtained by the FSSH method using a local diabatization algorithm. Moreover, the electronic populations of the full-dimensional on-the-fly calculations are in excellent agreement with the populations of the three-state nine-mode model, which confirms that the internal conversion dynamics of pyrazine is accurately represented by this reduced-dimensional model on the time scale under consideration (200 fs). The original FSSH method, in which the electronic wave function is propagated in the adiabatic representation, yields less accurate results. The oscillations in the populations of the diabatic B3u(nπ*) and A1u(nπ*) states driven by the mode Q8a are also observed in the full-dimensional dynamics simulations.

Imaging large amplitude out-of-plane motion in photoexcited pentafluorobenzene using time-resolved photoelectron spectroscopy: a computational study.

Excited-state dynamics of pentafluorobenzene is studied in detail for a quartic vibronic coupling model including the six b1 vibrational modes of the molecule and the two lowest excited electronic states. The study analyzes the influence of the large-amplitude out-of-plane vibrational motion on the electronic dynamics and extends to the simulation of the emerging time-resolved photoelectron spectra. The mapping of coherent non-separable electron-nuclear motion into oscillatory photoelectron signals is discussed.

Quantum dynamics of the photostability of pyrazine.

We investigate the radiationless decay of photoexcited pyrazine to its ground electronic state using multireference electronic structure and quantum dynamics calculations. We construct a quadratic vibronic coupling Hamiltonian, including the four lowest electronic states and ten vibrational modes, by fitting to more than 5000 ab initio points. We then use this model to simulate the non-adiabatic excited state dynamics of the molecule using the multi-configuration time-dependent Hartree method. On the basis of these calculations, we propose a new mechanism for this decay process involving a conical intersection between the Au(nπ*) state and the ground state. After excitation to the B2u(ππ*) state, the molecule decays to both the B3u(nπ*) and Au(nπ*) states on an ultrashort timescale of approximately 20 fs. The radiationless decay to the ground state then occurs from the Au(nπ*) state on a much longer timescale.

Comparing the electronic relaxation dynamics of aniline and d(7)-aniline following excitation at 272-238 nm.

Femtosecond time-resolved photoelectron spectroscopy experiments have been used to compare the electronic relaxation dynamics of aniline and d7-aniline following photoexcitation in the range 272-238 nm. Together with the results of recent theoretical investigations of the potential energy landscape [M. Sala, O. M. Kirkby, S. Guérin and H. H. Fielding, Phys. Chem. Chem. Phys., 2014, 16, 3122], these experiments allow us to resolve a number of unanswered questions surrounding the nonradiative relaxation mechanism. We find that tunnelling does not play a role in the electronic relaxation dynamics, which is surprising given that tunnelling plays an important role in the electronic relaxation of isoelectronic phenol and in pyrrole. We confirm the existence of two time constants associated with dynamics on the 1(1)πσ* surface that we attribute to relaxation through a conical intersection between the 1(1)πσ* and 1(1)ππ* states and motion on the 1(1)πσ* surface. We also present what we believe is the first report of an experimental signature of a 3-state conical intersection involving the 2(1)ππ*, 1(1)πσ* and 1(1)ππ* states.

Coherent destruction of tunneling in a six-dimensional model of NHD2: a computational study using the multi-configuration time-dependent Hartree method.

We investigate the phenomenon of coherent destruction of tunneling in a six-dimensional model of the NHD2 molecule. Two regimes are considered for the frequency of the laser field. A non-resonant regime where the frequency of the laser field is high with respect to the ground vibrational state tunneling splitting but smaller than the transition frequencies between the ground and excited vibrational states; and a quasi-resonant regime where the frequency of the laser field is close to the transition frequency between the ground and first excited vibrational states. In each case, we study the laser driven dynamics in the framework of the Floquet formalism and derive simple analytical formulas that explain the shape of the quasienergy curves associated with the two tunneling components of the ground vibrational state. This analysis allows us to obtain the parameters (frequency and amplitude) of the laser field that lead to the coherent destruction of tunneling. The multi-configuration time-dependent Hartree method is then used to solve the time-dependent Schrödinger equation for a six-dimensional model of the molecule in interaction with an adiabatically turned on monochromatic laser field, in order to confirm the results obtained from this analysis.

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

We present a full quantum-mechanical study of the laser control of the radiationless decay between the B3u(nπ(*)) and B2u(ππ(*)) states of pyrazine using the dynamic Stark effect. In contrast to our previous study [Sala et al., J. Chem. Phys. 140, 194309 (2014)], where a four-dimensional model was used, all the 24 degrees of freedom are now included in order to test the robustness of the strategy of control. Using a vibronic coupling Hamiltonian model in a diabatic representation, the multi-layer version of the multi-configuration time-dependent Hartree method is exploited to propagate the corresponding wave packets. We still observe a trapping of the wavepacket on the B2u(ππ(*)) potential energy surface due to the Stark effect for a longer time than the "non-resonant field-free" B2u(ππ(*)) lifetime.

The role of the low-lying dark nπ* states in the photophysics of pyrazine: a quantum dynamics study.

The excited state dynamics of pyrazine has attracted considerable attention in the last three decades. It has long been recognized that after UV excitation, the dynamics of the molecule is impacted by strong non-adiabatic effects due to the existence of a conical intersection between the B2u(ππ*) and B3u(nπ*) electronic states. However, a recent study based on trajectory surface hopping dynamics simulations suggested the participation of the Au(nπ*) and B2g(nπ*) low-lying dark electronic states in the ultrafast radiationless decay of the molecule after excitation to the B2u(ππ*) state. The purpose of this work was to pursue the investigation of the role of the Au(nπ*) and B2g(nπ*) states in the photophysics of pyrazine. A linear vibronic coupling model hamiltonian including the four lowest excited electronic states and the sixteen most relevant vibrational degrees of freedom was constructed using high level XMCQDPT2 electronic structure calculations. Wavepacket propagations using the MCTDH method were then performed and used to simulate the absorption spectrum and the electronic state population dynamics of the system. Our results show that the Au(nπ*) state plays an important role in the photophysics of pyrazine.

Laser control of the radiationless decay in pyrazine using the dynamic Stark effect.

The laser control of the radiationless decay between the B(3u)(nπ*) and B(2u)(ππ*) states of pyrazine using the dynamic Stark effect has been investigated. A vibronic coupling model Hamiltonian in diabatic representation, including potential energy, transition dipole, and static polarizability surfaces as a function of the four most important vibrational modes of the molecule has been parametrized using multi-reference electronic structure calculations. The interaction of the molecule with a strong non-resonant laser pulse has been analyzed in terms of dressed potential energy surfaces. Because of the large polarizability difference between the vibronically coupled B(3u)(nπ*) and B(2u)(ππ*) states, the Stark effect induced by the non-resonant laser pulse shifts the conical intersection away from the Franck-Condon region. We have shown, by solving the time-dependent Schrödinger equation for the molecule interacting with a relatively weak pump pulse driving the electronic excitation from the ground state to the B(2u)(ππ*) state, and a strong non-resonant control pulse, that this control mechanism can be used to trap the wavepacket on the B(2u)(ππ*) potential energy surface for a much longer time than the natural B(2u)(ππ*) lifetime.

New insight into the potential energy landscape and relaxation pathways of photoexcited aniline from CASSCF and XMCQDPT2 electronic structure calculations.

There have been a number of recent experimental investigations of the nonadiabatic relaxation dynamics of aniline following excitation to the first three singlet excited states, 1(1)ππ*, 1(1)π3s/πσ* and 2(1)ππ*. Motivated by differences between the interpretations of experimental observations, we have employed CASSCF and XMCQDPT2 calculations to explore the potential energy landscape and relaxation pathways of photoexcited aniline. We find a new prefulvene-like MECI connecting the 1(1)ππ* state with the GS in which the carbon-atom carrying the amino group is distorted out-of-plane. This suggests that excitation above the 1(1)π3s/πσ* vertical excitation energy could be followed by electronic relaxation from the 1(1)ππ* state to the ground-electronic state through this MECI. We find a MECI connecting the 1(1)π3s/πσ* and 1(1)ππ* states close to the local minimum on 1(1)π3s/πσ* which suggests that photoexcitation to the 1(1)π3s/πσ* state could be followed by relaxation to the 1(1)ππ* state and to the dissociative component of the 1(1)π3s/πσ* state. We also find evidence for a new pathway from the 2(1)ππ* state to the ground electronic state that is likely to pass through a three-state conical intersection involving the 2(1)ππ*, 1(1)π3s/πσ* and 1(1)ππ* states.

Laser-induced enhancement of tunneling in NHD2.

We apply and explore techniques aiming at enhancing the tunneling by laser fields, originally developed for a one-dimensional model, to a complete six-dimensional vibrational model of the inversion motion in NHD(2). The computational study is performed with the multi-configuration time-dependent Hartree method. Assuming an ideal three-dimensional alignment we obtain a driven tunneling time twenty times smaller than the natural one, in rather good agreement with an oversimplified three-state model. In the case of one-dimensional alignment, a linearly polarized field leads to a poor enhancement of the tunneling probability, after averaging over the rotation about the alignment axis, whereas a circularly polarized field improves the rotationally averaged tunneling probability at the end of the pulse.

Effect of the overall rotation on the cis-trans isomerization of HONO induced by an external field.

Rovibrational eigenenergies of HONO are computed and compared to experimental energies available in the literature. For their computation, we use a previously developed potential energy surface (PES) and a newly derived exact kinetic energy operator (KEO) including the overall rotation for a tetra-atomic molecule in non-orthogonal coordinates. In addition, we use the Heidelberg Multi-Configuration Time-Dependent Hartree (MCTDH) package. We compare the experimental rovibrational eigenvalues of HONO available in the literature with those obtained with MCTDH and a previously developed potential energy surface (PES) [F. Richter et al., J. Chem. Phys., 2004, 120, 1306.] for the cis geometry. The effect of the overall rotation on the process studied in our previous work on HONO [F. Richter et al., J. Chem. Phys., 2007, 127, 164315.] leading to the cis→trans isomerization of HONO is investigated. This effect on this process is found to be weak.