Fully quantal treatment of nonadiabatic molecular photodynamics: General considerations and application to the benzene cation

An Exhaustive Quantum-Classical Study of C6F6+ Using the Newly Formulated Parallel TDDVR Method

We recently implemented our parallelized quantum-classical dynamical approach, known as the Time-Dependent Discrete Variable Representation (TDDVR) method, which is applied to the spectroscopically important hexafluorobenzene (HFBz) radical cation, where several conical intersections exist in their seven lowest excited electronic states (S11B2u, S21E1g, S31B1u, S41E1u, and S51A2u) considering degeneracy among potential energy surfaces (PESs), to demonstrate their various dynamical aspects. This new parallel version shows almost linear scalability with increasing number of computing processors. To get photoelectron (PE) spectra, Mass-Analyzed Threshold Ionization (MATI) spectra, population dynamics, and many other dynamical observables, the first-principles dynamics is applied at the state-of-the-art level to the corresponding Hamiltonian, where the Jahn-Teller (JT) and pseudo-Jahn-Teller (PJT) type interactions are involved in those coupled seven electronic states. The quantum-classical method is used to thoroughly analyze the effects of these couplings on the nuclear dynamics of the involved electronic states, and the findings are compared with those observables obtained from experiments. Intrinsic dynamical properties are explained using the reduced densities of the wave packet (WP) in a coupled electronic manifold. The PE and MATI spectra of HFBz computed using TDDVR are found to be in good agreement with earlier experimental data and other theoretically simulated spectra.

Surface Hopping Dynamics on Vibronic Coupling Models

The simulation of photoinduced non-adiabatic dynamics is of great relevance in many scientific disciplines, ranging from physics and materials science to chemistry and biology. Upon light irradiation, different relaxation processes take place in which electronic and nuclear motion are intimately coupled. These are best described by the time-dependent molecular Schrödinger equation, but its solution poses fundamental practical challenges to contemporary theoretical chemistry. Two widely used and complementary approaches to this problem are multiconfigurational time-dependent Hartree (MCTDH) and trajectory surface hopping (SH). MCTDH is an accurate fully quantum-mechanical technique but often is feasible only in reduced dimensionality, in combination with approximate vibronic coupling (VC) Hamiltonians, or both (i.e., reduced-dimensional VC potentials). In contrast, SH is a quantum–classical technique that neglects most nuclear quantum effects but allows nuclear dynamics in full dimensionality by calculating potential energy surfaces on the fly. If nuclear quantum effects do not play a central role and a linear VC (LVC) Hamiltonian is appropriate—e.g., for stiff molecules that generally keep their conformation in the excited state—then it seems advantageous to combine the efficient LVC and SH techniques. In this Account, we describe how surface hopping based on an LVC Hamiltonian (SH/LVC)—as recently implemented in the SHARC surface hopping package—can provide an economical and automated approach to simulate non-adiabatic dynamics. First, we illustrate the potential of SH/LVC in a number of showcases, including intersystem crossing in SO₂, intra-Rydberg dynamics in acetone, and several photophysical studies on large transition-metal complexes, which would be much more demanding or impossible to perform with other methods. While all of the applications provide very useful insights into light-induced phenomena, they also hint at difficulties faced by the SH/LVC methodology that need to be addressed in the future. Second, we contend that the SH/LVC approach can be useful to benchmark SH itself. By the use of the same (LVC) potentials as MCTDH calculations have employed for decades and by relying on the efficiency of SH/LVC, it is possible to directly compare multiple SH test calculations with a MCTDH reference and ponder the accuracy of various correction algorithms behind the SH methodology, such as decoherence corrections or momentum rescaling schemes. Third, we demonstrate how the efficiency of SH/LVC can also be exploited to identify essential nuclear and electronic degrees of freedom to be employed in more accurate MCTDH calculations. Lastly, we show that SH/LVC is able to advance the development of SH protocols that can describe nuclear dynamics including explicit laser fields—a very challenging endeavor for trajectory-based schemes. To end, this Account compiles the typical costs of contemporary SH simulations, evidencing the great advantages of using parametrized potentials. The LVC model is a sleeping beauty that, kissed by SH, is fueling the field of excited-state molecular dynamics. We hope that this Account will stimulate future research in this direction, leveraging the advantages of the SH/VC schemes to larger extents and extending their applicability to uncharted territories.

Exploring vibronic coupling in the benzene radical cation and anion with different levels of the GW approximation

The linear vibronic coupling constants of the benzene radical cation and anion have been obtained with different levels of the GW approximation, including G₀W₀, eigenvalue self-consistent GW, and quasiparticle self-consistent GW, as well as DFT with the following exchange-correlation functionals: BLYP, B3LYP, CAM-B3LYP, tuned CAM-B3LYP, and an IP-tuned CAM-B3LYP functional. The vibronic coupling constants were calculated numerically using the gradients of the eigenvalues of the degenerate HOMOs and LUMOs of the neutral benzene molecule for DFT, while the numerical gradients of the quasiparticle energies were used in the case of GW. The results were evaluated against those of high level wave function methods in the literature, and the approximate self-consistent GW methods and G₀W₀ with long-range corrected functionals were found to yield the best results on the whole.

Accurate prediction and measurement of vibronic branching ratios for laser cooling linear polyatomic molecules

We report a generally applicable computational and experimental approach to determine vibronic branching ratios in linear polyatomic molecules to the 10^-5 level, including for nominally symmetry-forbidden transitions. These methods are demonstrated in CaOH and YbOH, showing approximately two orders of magnitude improved sensitivity compared with the previous state of the art. Knowledge of branching ratios at this level is needed for the successful deep laser cooling of a broad range of molecular species.

Table-Top X-ray Spectroscopy of Benzene Radical Cation

Ultrafast table-top X-ray spectroscopy at the carbon K-edge is used to measure the X-ray spectral features of benzene radical cations (Bz⁺). The ground state of the cation is prepared selectively by two-photon ionization of neutral benzene, and the X-ray spectra are probed at early times after the ionization by transient absorption using X-rays produced by high harmonic generation (HHG). Bz⁺ is well-known to undergo Jahn-Teller distortion, leading to a lower symmetry and splitting of the π orbitals. Comparison of the X-ray absorption spectra of the neutral and the cation reveals a splitting of the two degenerate π* orbitals as well as an appearance of a new peak due to excitation to the partially occupied π-subshell. The π* orbital splitting of the cation, elucidated on the basis of high-level calculations in a companion theoretical paper [Vidal et al. J. Phys. Chem. A. http://dx.doi.org/10.1021/acs.jpca.0c08732], is discovered to be due to both the symmetry distortion and even more dominant spin coupling of the unpaired electron in the partially vacant π orbital (from ionization) with the unpaired electrons resulting from the transition from the 1s_C core orbital to the fully vacant π* orbitals.

Interplay of Open-Shell Spin-Coupling and Jahn-Teller Distortion in Benzene Radical Cation Probed by X-ray Spectroscopy

We report a theoretical investigation and elucidation of the X-ray absorption spectra of neutral benzene and of the benzene cation. The generation of the cation by multiphoton ultraviolet (UV) ionization and the measurement of the carbon K-edge spectra of both species using a table-top high-harmonic generation source are described in the companion experimental paper [Epshtein, M.; et al. J. Phys. Chem. A http://dx.doi.org/10.1021/acs.jpca.0c08736]. We show that the 1s_C → π transition serves as a sensitive signature of the transient cation formation, as it occurs outside of the spectral window of the parent neutral species. Moreover, the presence of the unpaired (spectator) electron in the π-subshell of the cation and the high symmetry of the system result in significant differences relative to neutral benzene in the spectral features associated with the 1s_C → π* transitions. High-level calculations using equation-of-motion coupled-cluster theory provide the interpretation of the experimental spectra and insight into the electronic structure of benzene and its cation. The prominent split structure of the 1s_C → π* band of the cation is attributed to the interplay between the coupling of the core → π* excitation with the unpaired electron in the π-subshell and the Jahn-Teller distortion. The calculations attribute most of the splitting (∼1-1.2 eV) to the spin coupling, which is visible already at the Franck-Condon structure, and we estimate the additional splitting due to structural relaxation to be around ∼0.1-0.2 eV. These results suggest that X-ray absorption with increased resolution might be able to disentangle electronic and structural aspects of the Jahn-Teller effect in the benzene cation.