Theory of vibronic coupling in linear molecules

Resolving Photoinduced Femtosecond Three-Dimensional Solute-Solvent Dynamics through Surface Hopping Simulations

Photoinduced dynamics in solution is governed by mutual solute-solvent interactions, which give rise to phenomena like solvatochromism, the Stokes shift, dual fluorescence, or charge transfer. Understanding these phenomena requires simulating the solute's photoinduced dynamics and simultaneously resolving the three-dimensional solvent distribution dynamics. If using trajectory surface hopping (TSH) to this aim, thousands of trajectories are required to adequately sample the time-dependent three-dimensional solvent distribution functions, and thus resolve the solvent dynamics with sub-Ångstrom and femtosecond accuracy and sufficiently low noise levels. Unfortunately, simulating thousands of trajectories with TSH in the framework of hybrid quantum mechanical/molecular mechanical (QM/MM) can be prohibitively expensive when employing ab initio electronic structure methods. To tackle this challenge, we recently introduced a computationally efficient approach that combines efficient linear vibronic coupling models with molecular mechanics (LVC/MM) via electrostatic embedding [Polonius et al., JCTC 2023, 19, 7171-7186]. This method provides solvent-embedded, nonadiabatically coupled potential energy surfaces while scaling similarly to MM force fields. Here, we employ TSH with LVC/MM to unravel the photoinduced dynamics of two small thiocarbonyl compounds solvated in water. We describe how to estimate the number of trajectories required to produce nearly noise-free three-dimensional solvent distribution functions and present an analysis based on approximately 10,000 trajectories propagated for 3 ps. In the electronic ground state, both molecules exhibit in-plane hydrogen bonds to the sulfur atom. Shortly after excitation, these bonds are broken and reform perpendicular to the molecular plane on timescales that differ by an order of magnitude due to steric effects. We also show that the solvent relaxation dynamics is coupled to the electronic dynamics, including intersystem crossing. These findings are relevant to advance the understanding of the coupled solute-solvent dynamics of solvated photoexcited molecules, e.g., biologically relevant thio-nucleobases.

Time dependent vibrational electronic coupled cluster (VECC) theory for non-adiabatic nuclear dynamics

A time-dependent vibrational electronic coupled-cluster (VECC) approach is proposed to simulate photo-electron/UV-VIS absorption spectra as well as time-dependent properties for non-adiabatic vibronic models, going beyond the Born-Oppenheimer approximation. A detailed derivation of the equations of motion and a motivation for the ansatz are presented. The VECC method employs second-quantized bosonic construction operators and a mixed linear and exponential ansatz to form a compact representation of the time-dependent wave-function. Importantly, the method does not require a basis set, has only a few user-defined inputs, and has a classical (polynomial) scaling with respect to the number of degrees of freedom (of the vibronic model), resulting in a favorable computational cost. In benchmark applications to small models and molecules, the VECC method provides accurate results compared to multi-configurational time-dependent Hartree calculations when predicting short-time dynamical properties (i.e., photo-electron/UV-VIS absorption spectra) for non-adiabatic vibronic models. To illustrate the capabilities, the VECC method is also successfully applied to a large vibronic model for hexahelicene with 14 electronic states and 63 normal modes, developed in the group by Aranda and Santoro [J. Chem. Theory Comput. 17, 1691, (2021)].

LVC/MM: A Hybrid Linear Vibronic Coupling/Molecular Mechanics Model with Distributed Multipole-Based Electrostatic Embedding for Highly Efficient Surface Hopping Dynamics in Solution

We present a theoretical framework for a hybrid linear vibronic coupling model electrostatically embedded into a molecular mechanics environment, termed the linear vibronic coupling/molecular mechanics (LVC/MM) method, for the surface hopping including arbitrary coupling (SHARC) molecular dynamics package. Electrostatic embedding is realized through the computation of interactions between environment point charges and distributed multipole expansions (DMEs, up to quadrupoles) that represent each electronic state and transition densities in the diabatic basis. The DME parameters are obtained through a restrained electrostatic potential (RESP) fit, which we extended to yield higher-order multipoles. We also implemented in SHARC a scheme for achieving roto-translational invariance of LVC models as well as a general quantum mechanics/molecular mechanics (QM/MM) interface, an OpenMM interface, and restraining potentials for simulating liquid droplets. Using thioformaldehyde in water as a test case, we demonstrate that LVC/MM can accurately reproduce the solvation structure and energetics of rigid solutes, with errors on the order of 1-2 kcal/mol compared to a BP86/MM reference. The implementation in SHARC is shown to be very efficient, enabling the simulation of trajectories on the nanosecond time scale in a matter of days.

A vibrational action spectroscopic study of the Renner-Teller- and spin-orbit-affected cyanoacetylene radical cation HC3N. J Chem Phys 2023;158:084305. [PMID: 36859081 DOI: 10.1063/5.0135000]

The linear radical cation of cyanoacetylene, HC₃N⁺ (²Π), is not only of astrophysical interest, being the, so far undetected, cationic counterpart of the abundant cyanoaceteylene, but also of fundamental spectroscopic interest due to its strong spin-orbit and Renner-Teller interactions. Here, we present the first broadband vibrational action spectroscopic investigation of this ion through the infrared pre-dissociation (IRPD) method using a Ne tag. Experiments have been performed using the FELion cryogenic ion-trap instrument in combination with the FELIX free-electron lasers and a Laservision optical parametric oscillator/optical parametric amplifier system. The vibronic splitting patterns of the three interacting bending modes (ν₅, ν₆, ν₇), ranging from 180 to 1600 cm^-1, could be fully resolved revealing several bands that were previously unobserved. The associated Renner-Teller and intermode coupling constants have been determined by fitting an effective Hamiltonian to the experimental data, and the obtained spectroscopic constants are in reasonable agreement with previous photoelectron spectroscopy (PES) studies and ab initio calculations on the HC₃N⁺ ion. The influence of the attached Ne atom on the infrared spectrum has been investigated by ab initio calculations at the RCCSD(T)-F12a level of theory, which strongly indicates that the discrepancies between the IRPD and PES data are a result of the effects of the Ne attachment.

Vibronic Coupling through the Continuum in the e+CO_{2} System

We report two-dimensional electron energy-loss spectra of CO_{2}. The high-resolution experiment reveals a counterintuitive fine structure at energy losses where CO_{2} states form a vibrational pseudocontinuum. Guided by the symmetry of the system, we constructed a four-dimensional nonlocal model for the vibronic dynamics involving two shape resonances (forming a Renner-Teller Π_{u} doublet at the equilibrium geometry) coupled to a virtual Σ_{g}^{+} state. The model elucidates the extremely non-Born-Oppenheimer dynamics of the coupled nuclear motion and explains the origin of the observed structures. It is a prototype of the vibronic coupling of metastable states in continuum.

Vibronic Coupling in Spherically Encapsulated, Diatomic Molecules: Prediction of a Renner-Teller-like Effect for Endofullerenes

In the year 1933, Herzberg and Teller realized that the potential energy surface of a triatomic, linear molecule splits into two as soon as the molecule is bent. The phenomenon, later dubbed the Renner–Teller effect due to the detailed follow-up work of Renner on the subject, describes the coupling of a symmetry-reducing molecular vibration with degenerate electronic states. In this article, we show that a very similar type of nonadiabatic coupling can occur for certain translational degrees of freedom of diatomic, electronically degenerate molecules when trapped in a nearly spherical or cylindrical quantum confinement, e.g., realized through electromagnetic fields or molecular encapsulation. We illustrate this on the example of fullerene-encapsulated nitric oxide, and provide a prediction of its interesting, perturbed vibronic spectrum.

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.

How important are the residual nonadiabatic couplings for an accurate simulation of nonadiabatic quantum dynamics in a quasidiabatic representation?

Diabatization of the molecular Hamiltonian is a standard approach to remove the singularities of nonadiabatic couplings at conical intersections of adiabatic potential energy surfaces. In general, it is impossible to eliminate the nonadiabatic couplings entirely-the resulting "quasidiabatic" states are still coupled by smaller but nonvanishing residual nonadiabatic couplings, which are typically neglected. Here, we propose a general method for assessing the validity of this potentially drastic approximation by comparing quantum dynamics simulated either with or without the residual couplings. To make the numerical errors negligible to the errors due to neglecting the residual couplings, we use the highly accurate and general eighth-order composition of the implicit midpoint method. The usefulness of the proposed method is demonstrated on nonadiabatic simulations in the cubic Jahn-Teller model of nitrogen trioxide and in the induced Renner-Teller model of hydrogen cyanide. We find that, depending on the system, initial state, and employed quasidiabatization scheme, neglecting the residual couplings can result in wrong dynamics. In contrast, simulations with the exact quasidiabatic Hamiltonian, which contains the residual couplings, always yield accurate results.

On the generation of electronic ring currents under vibronic coupling effects

We study the generation of electronic ring currents in the presence of nonadiabatic coupling using circularly polarized light. For this, we introduce a solvable model consisting of an electron and a nucleus rotating around a common center and subject to their mutual Coulomb interaction. The simplicity of the model brings to the forefront the non-trivial properties of electronic ring currents in the presence of coupling to the nuclear coordinates and enables the characterization of various limiting situations transparently. Employing this model, we show that vibronic coupling effects play a crucial role even when a single E degenerate eigenstate of the system supports the current. The maximum current of a degenerate eigenstate depends on the strength of the nonadiabatic interactions. In the limit of large nuclear to electronic masses, in which the Born-Oppenheimer approximation becomes exact, constant ring currents and time-averaged oscillatory currents necessarily vanish.

Stability and Dissociation of Ethylenedione (OCCO)

This work examines the electronic structure and apparent instability of ethylenedione (OCCO), including an analysis of the singlet and triplet potential energy surfaces along the bending vibrations. While the singlet state is inherently unstable due to the Renner-Teller effect, theory predicts the triplet state to have a stable minimum on the potential energy surface. The stability of the triplet state is examined in detail, taking into account spin-orbit interactions. Using multireference quantum chemical methods, the lifetime of the triplet state is estimated to be in the picosecond range, significantly lower than previously computed. A quasi-atomic molecular orbital (QUAO) analysis is also used to elucidate the nature of bonding along the potential energy surface in both the singlet and triplet states. These calculations confirm the transient nature of the OCCO molecule, although they do not fully explain the lack of experimental detection via spectroscopy, which is known have the capability to probe even shorter lifetimes.

Signatures of vibronic coupling in two-dimensional electronic-vibrational and vibrational-electronic spectroscopies

Two-Dimensional Electronic-Vibrational (2D EV) spectroscopy and Two-Dimensional Vibrational-Electronic (2D VE) spectroscopy are new coherent four-wave mixing spectroscopies that utilize both electronically resonant and vibrationally resonant field-matter interactions to elucidate couplings between electronic and vibrational degrees of freedom. A system Hamiltonian is developed here to lay a foundation for interpreting the 2D EV and 2D VE signals that arise from a vibronically coupled molecular system in the condensed phase. A molecular system consisting of one anharmonic vibration and two electronic states is modeled. Equilibrium displacement of the vibrational coordinate and vibrational frequency shifts upon excitation to the first electronic excited state are included in our Hamiltonian through linear and quadratic vibronic coupling terms. We explicitly consider the nuclear dependence of the electronic transition dipole moment and demonstrate that these spectroscopies are sensitive to non-Condon effects. A series of simulations of 2D EV and 2D VE spectra obtained by varying parameters of the system, system-bath, and interaction Hamiltonians demonstrate that one of the following conditions must be met to observe signals: (1) non-zero linear and/or quadratic vibronic coupling in the electronic excited state, (2) vibrational-coordinate dependence of the electronic transition dipole moment, or (3) electronic-state-dependent vibrational dephasing dynamics. We explore how these vibronic interactions are manifested in the positions, amplitudes, and line shapes of the peaks in 2D EV and 2D VE spectroscopies.

Practical approximation of the non-adiabatic coupling terms for same-symmetry interstate crossings by using adiabatic potential energies only

A very simple equation, F_ij^App=[(∂²(V_i^a-V_j^a)/∂Q²)/(V_i^a-V_j^a)]^1/2/2, giving a reliable magnitude of non-adiabatic coupling terms (NACTs, F_ij's) based on adiabatic potential energies only (V_i^a and V_j^a) was discovered, and its reliability was tested for several prototypes of same-symmetry interstate crossings in LiF, C₂, NH₃Cl, and C₆H₅SH molecules. Our theoretical derivation starts from the analysis of the relationship between the Lorentzian dependence of NACTs along a diabatization coordinate and the well-established linear vibronic coupling scheme. This analysis results in a very simple equation, α=2κ/Δ_c, enabling the evaluation of the Lorentz function α parameter in terms of the coupling constant κ and the energy gap Δ_c (Δ_c=|V_i^a-V_j^a|_Qc ) between adiabatic states at the crossing point Q_C. Subsequently, it was shown that Q_C corresponds to the point where F_ij^App exhibit maximum values if we set the coupling parameter as κ=[(V_i^a-V_j^a)⋅(∂²(V_i^a-V_j^a)/∂Q²)]_Qc^1/2/2. Finally, we conjectured that this relation could give reasonable values of NACTs not only at the crossing point but also at other geometries near Q_C. In this final approximation, the pre-defined crossing point Q_C is not required. The results of our test demonstrate that the approximation works much better than initially expected. The present new method does not depend on the selection of an ab initio method for adiabatic electronic states but is currently limited to local non-adiabatic regions where only two electronic states are dominantly involved within a nuclear degree of freedom.

A benchmark study of electronic excitation energies, transition moments, and excited-state energy gradients on the nicotine molecule

In this work, we report a comparative study of computed excitation energies, oscillator strengths, and excited-state energy gradients of (S)-nicotine, chosen as a test case, using multireference methods, coupled cluster singles and doubles, and methods based on time-dependent density functional theory. This system was chosen because its apparent simplicity hides a complex electronic structure, as several different types of valence excitations are possible, including n-π(*), π-π(*), and charge-transfer states, and in order to simulate its spectrum it is necessary to describe all of them consistently well by the chosen method.

Probing the electronic structure and Au-C chemical bonding in AuC2(-) and AuC2 using high-resolution photoelectron spectroscopy

We report photoelectron spectroscopy (PES) and high-resolution PE imaging of AuC2(-) at a wide range of photon energies. The ground state of AuC2(-) is found to be linear (C∞v, (1)Σ(+)) with a …8π(4)4δ(4)17σ(2)9π(4)18σ(2) valence configuration. Detachments from all the five valence orbitals of the ground state of AuC2(-) are observed at 193 nm. High-resolution PE images are obtained in the energy range from 830 to 330 nm, revealing complicated vibronic structures from electron detachment of the 18σ, 9π, and 17σ orbitals. Detachment from the 18σ orbital results in the (2)Σ(+) ground state of neutral AuC2, which, however, is bent due to strong vibronic coupling with the nearby (2)Π state from detachment of a 9π electron. The (2)Σ(+)-(2)Π vibronic and spin-orbit coupling results in complicated vibronic structures for the (2)Σ(+) and (2)Π3/2 states with extensive bending excitations. The electron affinity of AuC2 is measured accurately to be 3.2192(7) eV with a ground state bending frequency of 195(6) cm(-1). The first excited state ((2)A') of AuC2, corresponding to the (2)Π3/2 state at the linear geometry, is only 0.0021 eV above the ground state ((2)A') and has a bending frequency of 207(6) cm(-1). The (2)Π1/2 state, 0.2291 eV above the ground state, is linear with little geometry change relative to the anion ground state. The detachment of the 17σ orbital also results in complicated vibronic structures, suggesting again a bent state due to possible vibronic coupling with the lower (2)Π state. The spectrum at 193 nm shows the presence of a minor species with less than 2% intensity relative to the ground state of AuC2(-). High-resolution data of the minor species reveal several vibrational progressions in the Au-C stretching mode, which are assigned to be from the metastable (3)Π2,1,0 spin-orbit excited states of AuC2(-) to the (2)Π3/2,1/2 spin-orbit states of neutral AuC2. The spin-orbit splittings of the (3)Π and (2)Π states are accurately measured at the linear geometry. The current study provides a wealth of electronic structure information about AuC2(-) and AuC2, which are ideal systems to investigate the strong Σ-Π and spin-orbit vibronic couplings.

Dissociative recombination of highly symmetric polyatomic ions

A general first-principles theory of dissociative recombination is developed for highly symmetric molecular ions and applied to H(3)O(+) and CH(3)(+), which play an important role in astrophysical, combustion, and laboratory plasma environments. The theoretical cross sections obtained for the dissociative recombination of the two ions are in good agreement with existing experimental data from storage ring experiments.

Spin-orbit vibronic coupling in Π3 states of linear triatomic molecules

The Renner-Teller vibronic-coupling problem of a 3Pi electronic state of a linear molecule is analyzed with the inclusion of the spin-orbit coupling of the 3Pi electronic state, employing the microscopic (Breit-Pauli) spin-orbit coupling operator for the two unpaired electrons. The 6x6 Hamiltonian matrix in a diabatic spin-electronic basis is obtained by an expansion of the molecular Hamiltonian in powers of the bending amplitude. The symmetry properties of the Hamiltonian with respect to the time-reversal operator and the relativistic vibronic angular momentum operator are analyzed. It is shown that there exists a linear vibronic-coupling term of spin-orbit origin, which has not been considered so far in the Renner-Teller theory of 3Pi electronic states. While two of the six adiabatic electronic wave functions do not exhibit a geometric phase, the other four carry nontrivial topological phases which depend on the radius of the integration contour. The spectroscopic effects of the linear spin-orbit vibronic-coupling mechanism have been analyzed by numerical calculations of the vibronic spectrum for selected model examples.

Calculation of the vibronic structure of the photodetachment spectra of CCCl- and CCBr-

The vibronic structure of the closely spaced and strongly coupled X 2Sigma+ and A 2Pi states in the photodetachment spectra of CCCl- and CCBr- has been calculated by considering Sigma-Pi vibronic coupling together with spin-orbit coupling. The stretching modes are treated within the so-called linear-vibronic-coupling model. The vibronic and spin-orbit parameters have been determined by accurate ab initio electronic-structure calculations. While the nonrelativistic vibronic-coupling parameters are of approximately equal strength in CCCl and CCBr, the vibronic-coupling parameters of spin-orbit origin are found to be larger in the latter. The calculated photodetachment spectra of both systems are shown to exhibit a complicated vibronic structure due to strong Sigma-Pi vibronic coupling. The spectral envelopes of the calculated photodetachment spectra exhibit a double-hump reminiscent of strongly coupled Exe Jahn-Teller systems.

Spectroscopic effects of first-order relativistic vibronic coupling in linear triatomic molecules

It has recently been shown that there exists, in addition to the well-known nonrelativistic Renner-Teller coupling, a linear (that is, of the first order in the bending distortion) vibronic-coupling mechanism of relativistic (that is, spin-orbit) origin in 2II electronic states of linear molecules [L. V. Poluyanov and W. Domcke, Chem. Phys. 301, 111 (2004)]. The generic aspects of the relativistic linear vibronic-coupling mechanism have been analyzed in the present work by numerical calculations of the vibronic spectrum for appropriate models. The vibronic and spin-orbit parameters have been determined by accurate ab initio electronic-structure calculations for the X 2II states of a series of triatomic radicals and radical cations. It is shown for the example of GeCH that the relativistic linear vibronic-coupling mechanism provides a quantitative explanation of the pronounced perturbations in the vibronic spectrum of the X 2II state of GeCH, which previously have been termed "Sears resonances" [S.-G. He, H. Li, T. C. Smith, D. J. Clouthier, and A. J. Merer, J. Chem. Phys. 119, 10115 (2003)]. The X 2II vibronic spectra of the series BS2, CS2+, OCS+, and OBS illustrate the interplay of nonrelativistic and relativistic vibronic-coupling mechanisms in Renner-Teller systems.