Vibronic Coupling Model Research Articles

Jahn-Teller and pseudo-Jahn-Teller effects on the vibronic structure of the photoionized spectrum of cyanopropyne.

Nonadiabatic quantum dynamics are carried out to illustrate the photoionized spectrum of the cyanopropyne (CH3-C≡C-C≡N) as reported in recent experimental measurements [Lamarre et al., J. Mol. Spectrosc. 315, 206 (2015)]. A detailed electronic structure calculation is performed to analyze the topographical details of the first five ionized states, of which three are degenerate states (X̃2E, B̃2E, and C̃2E) and two are non-degenerate states (Ã2A1 and D̃2A1). The degenerate E states of the C3V symmetry molecule are prone to Jahn-Teller (JT) instability, and in addition, symmetry allowed A1 - E vibronic coupling, i.e., pseudo-Jahn-Teller (PJT), effects are expected to have a significant impact in the detailed vibronic structure of these electronic states. The JT splittings of X̃2E and B̃2E degenerate states are small, whereas it is quite large at three high frequencies in the C̃2E electronic states. The large energy separation of X̃2E from the other states and the non-zero PJT coupling of the B̃2E state with the close-lying Ã2A1 state indicate the uncoupled nature of the X̃, Ã, and B̃ vibronic bands of C4H3N. The intersection minima of B̃ and C̃ states with the D̃ state nearly coincide with the energetic minimum of D̃ state. Therefore, the PJT couplings among these states will lead to a strong vibronic interaction to shape the respective band structure. To completely understand the JT and PJT interactions in the photoionized spectrum of C4H3N, the vibronic coupling model Hamiltonian was constructed to perform nuclear dynamics studies for these electronic states. The vibrational progressions in each vibronic band are identified and compared with the available experimental data in the literature. The impacts of JT and PJT effects in the first five ionized states of cyanopropyne are investigated and discussed in detail.

Photoionization of aziridine: Nonadiabatic dynamics of the first six low-lying electronic states of the aziridine radical cation.

In this article, the theoretical photoionization spectroscopy of the aziridine (C2H5N) molecule is investigated. To start with, we have optimized the geometry of this molecule at the neutral electronic ground state at the density functional theory/augmented correlation-consistent polarized valence triple zeta level of theory using the G09 program. The electronic structure calculations were restricted to the first six low-lying electronic states in order to account for the experimental photoelectron spectrum of the C2H5N molecule. The first six low-lying electronic states (X̃2A', Ã2A', B̃2A″, C̃2A″, D̃2A', and Ẽ2A') of the potential energy surfaces (PESs) are calculated by both equation of motion-ionization potential-coupled cluster singles and doubles and multi-configuration quasi-degenerate perturbation theory abinitio quantum chemistry methods along the dimensionless normal displacement coordinates in which multiple conical intersections were established among the considered electronic states. A (6 × 6) model vibronic Hamiltonian is constructed on a diabatic electronic basis, using the symmetry selection rules and Taylor series expansion. The Cs symmetry point group of the aziridine molecule leads to electronic states symmetry of either A' or A″, and these states are close in energy, due to which the same symmetry electronic states avoid each other. To get a smooth diabatic PES, a fourfold diabatization scheme is used, which is implemented in the General Atomic and Molecular Electronic Structure Systems suite of programs. All the parameters used in the diabatic vibronic coupling model Hamiltonian are calculated in terms of the normal modes of vibrational coordinates. Finally, the vibronic model Hamiltonian constructed for the coupled six electronic states is used to solve both time-independent and time-dependent Schrödinger equations using the multi-configuration time-dependent Hartree program module to obtain the dynamical observables. The theoretical vibronic band structure is found to be in good accord with the available experimental results.

Understanding molecular geometric phase effects with exact effective force: case study of a model system

In this work, molecular geometric phase effects are studied using the idea of exact factorization (EF) (Abediet al2010Phys. Rev. Lett.105123002) and exact effective force (Liet al2022Phys. Rev. Lett.128113001). In particular, we performed dynamics simulations for a two-state vibronic coupling model, and interpreted the results in three different perspectives: the Born-Huang expansion, the exact time-dependent potential energy surface (TDPES) and the exact effective force. We find that (i) at particular moment, while the vanishing nuclear density that occurs periodically in space is conventionally attributed to destructive interference of the nuclear wave packet owing to the geometric phase, such phenomenon can be equally well interpreted through the energy perspective, as manifested in the exact TDPES in the EF scheme; (ii) when combined with trajectory-based classical dynamics, the exact effective force obtained through EF qualitatively reproduces the correct nuclear density, while the adiabatic force gives the wrong density, particularly in the interference region. Our results suggest that the exact effective force is a potential starting point for making approximations and improving trajectory-based computational methods towards an accurate description of geometric phase effects.

Electronic dynamics through conical intersections via non-Markovian stochastic Schrödinger equation with complex modes.

Conical intersections (CIs) play a crucial role in photochemical reactions, offering an efficient channel for ultrafast non-adiabatic relaxation of excited states. This significantly influences the reaction pathways and the resulting products. In this work, we utilize the non-Markovian stochastic Schrödinger equation with complex modes method to explore the dynamics of electronic transitions through conical intersections (CIs) in pyrazine. The linear vibronic coupling model serves as the foundational framework, incorporating both intra-state and inter-state electron-vibrational interactions. The dynamics of the excited electronic transitions are analyzed across varying strengths of system-bath coupling and different bath relaxation times. The accuracy of this method is demonstrated by comparing its predictions with those from the hierarchical equations of motion method.

Ultrafast Excited-State Energy Transfer in Phenylene Ethynylene Dendrimer: Quantum Dynamics with the Tensor Network Method.

Photoinduced excited-state energy transfer (EET) processes play an important role in solar energy conversions. Owing to their excellent photoharvesting and exciton-transport properties, phenylene ethynylene (PE) dendrimers display great potential for improving the efficiency of solar cells. In this work, we investigated the intramolecular EET dynamics in a dendrimer composed of two linear PE units (2-ring and 3-ring) using a fully quantum description based on the tensor network method. We first constructed a diabatic model Hamiltonian based on the electronic structure calculations. Using this diabatic vibronic coupling model, we tried to obtain the main features of the EET dynamics in terms of the several diabatic models with different numbers of vibrational modes (from 4 modes to 129 modes) and to explore the corresponding vibronic coupling interactions. The results show that the EET in this PE dendrimer is ultrafast. Four modes of A' symmetry play dominant roles in the dynamics; the remaining 86 modes of A' symmetry can dampen the electronic coherence; and the modes of A″ symmetry do not exhibit significant influence on the EET process. Overall, the first-order intrastate vibronic coupling terms show the dominant role in the EET dynamics, while the second-order intrastate vibronic coupling terms cause damping of the electronic coherence and slow down the overall EET process. This work provides a microscopic understanding of the EET dynamics in PE dendrimers.

Benchmarking various nonadiabatic semiclassical mapping dynamics methods with tensor-train thermo-field dynamics.

Accurate quantum dynamics simulations of nonadiabatic processes are important for studies of electron transfer, energy transfer, and photochemical reactions in complex systems. In this comparative study, we benchmark various approximate nonadiabatic dynamics methods with mapping variables against numerically exact calculations based on the tensor-train (TT) representation of high-dimensional arrays, including TT-KSL for zero-temperature dynamics and TT-thermofield dynamics for finite-temperature dynamics. The approximate nonadiabatic dynamics methods investigated include mixed quantum-classical Ehrenfest mean-field and fewest-switches surface hopping, linearized semiclassical mapping dynamics, symmetrized quasiclassical dynamics, the spin-mapping method, and extended classical mapping models. Different model systems were evaluated, including the spin-boson model for nonadiabatic dynamics in the condensed phase, the linear vibronic coupling model for electronic transition through conical intersections, the photoisomerization model of retinal, and Tully's one-dimensional scattering models. Our calculations show that the optimal choice of approximate dynamical method is system-specific, and the accuracy is sensitively dependent on the zero-point-energy parameter and the initial sampling strategy for the mapping variables.

Mean-Field Ring Polymer Rates Using a Population Dividing Surface.

Mean-field ring polymer molecular dynamics offers a computationally efficient method for the simulation of reaction rates in multilevel systems. Previous work has established that, to model a nonadiabatic state-to-state reaction accurately, it is necessary to ensure reactive trajectories form kinked ring polymer configurations at the dividing surface. Building on this idea, we introduce a population difference coordinate and a reactive flux expression modified to only include contributions from kinked configurations. We test the accuracy of the resulting mean-field rate theory on a series of linear vibronic coupling model systems. We demonstrate that this new kMF-RP rate approach is efficient to implement and quantitatively accurate for models over a wide range of driving forces, coupling strengths, and temperatures.

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.

Prediction through quantum dynamics simulations: Photo-excited cyclobutanone.

Quantum dynamics simulations are becoming a standard tool for simulating photo-excited molecular systems involving a manifold of coupled states, known as non-adiabatic dynamics. While these simulations have had many successes in explaining experiments and giving details of non-adiabatic transitions, the question remains as to their predictive power. In this work, we present a set of quantum dynamics simulations on cyclobutanone using both grid-based multi-configuration time-dependent Hartree and direct dynamics variational multi-configuration Gaussian methods. The former used a parameterized vibronic coupling model Hamiltonian, and the latter generated the potential energy surfaces on the fly. The results give a picture of the non-adiabatic behavior of this molecule and were used to calculate the signal from a gas-phase ultrafast electron diffraction (GUED) experiment. Corresponding experimental results will be obtained and presented at a later stage for comparison to test the predictive power of the methods. The results show that over the first 500fs after photo-excitation to the S2 state, cyclobutanone relaxes quickly to the S1 state, but only a small population relaxes further to the S0 state. No significant transfer of population to the triplet manifold is found. It is predicted that the GUED experiments over this time scale will see signals related mostly to the C-O stretch motion and elongation of the molecular ring along the C-C-O axis.

Linear Vibronic Coupling Approach for Surface-Enhanced Raman Scattering: Quantifying the Charge-Transfer Enhancement Mechanism.

The outstanding amplification observed in surface-enhanced Raman scattering (SERS) is due to several enhancement mechanisms, and standing out among them are the plasmonic (PL) and charge-transfer (CT) mechanisms. The theoretical estimation of the enhancement factors of the CT mechanism is challenging because the excited-state coupling between bright plasmons and dark CT states must be properly introduced into the model to obtain reliable intensities. In this work, we aim at simulating electrochemical SERS spectra, considering models of pyridine on silver clusters subjected to an external electric field E⃗ that represents the effect of an electrode potential Vel. The method adopts quantum dynamical propagations of nuclear wavepackets on the coupled PL and CT states described with linear vibronic coupling models parametrized for each E⃗ through a fragment-based maximum-overlap diabatization. By presenting results at different values of E⃗, we show that indeed there is a relation between the population transferred to the CT states and the total scattered intensity. The tuning and detuning processes of the CT states with the bright PLs as a function of the electric field are in good agreement with those observed in experiments. Finally, our estimations for the CT enhancement factors predict values in the order of 105 to 106, meaning that when the CT and PL states are both in resonance with the excitation wavelength, the CT and PL enhancements are comparable, and vibrational bands whose intensity is amplified by different mechanisms can be observed together, in agreement with what was measured by typical experiments on silver electrodes.

Exact-Factorization-Based Surface Hopping without Velocity Adjustment.

While surface hopping has emerged as a powerful method for simulating non-adiabatic dynamics in large molecules, the ad hoc nature of the necessary velocity adjustments and decoherence corrections in the algorithm somewhat reduces its reliability. Here we propose a new scheme that eliminates these aspects by combining the nuclear equation from the quantum-trajectory surface-hopping approach with the electronic equation derived from the exact-factorization approach. The resulting method, denoted QTSH-XF, yields a surface-hopping method on firmer ground than previous and is shown to successfully capture dynamics in Tully models and in a linear vibronic coupling model of the photoexcited uracil cation.

Studies of nonadiabatic dynamics in the singlet fission processes of pentacene dimer via tensor network method

Singlet fission (SF) is a very significant photophysical phenomenon and possesses potential applications. In this work, we try to give a rather detailed theoretical investigation of the SF process in the stacked polyacene dimer by combining the high-level quantum chemistry calculations and the quantum dynamics simulations based on the tensor network method. Starting with the construction of the linear vibronic coupling model, we explore the pure electronic dynamics and the vibronic dynamics in the SF processes. The role of vibrational modes in nonadiabatic dynamics is addressed. The results show that the super-exchange mechanism mediated by the charge-transfer state is found in both pure electronic dynamics and the nonadiabatic dynamics. Particularly the vibrational modes with the frequencies resonance with the adiabatic energy gap play very import roles in the SF dynamics. This work not only provides a deep and detailed understanding of the SF process but also verifies the efficiency of the tensor network method with the train structure that can serve as the reference dynamics method to explore the dynamics behaviors of complex systems.

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.

Vibrational excitation in the e + CO2 system: Analysis of the two-dimensional energy-loss spectrum

We present a detailed analysis of the two-dimensional electron energy-loss spectrum of ${\mathrm{CO}}_{2}$, which extends our recent Letter [Phys. Rev. Lett. 129, 013401 (2022)]. We show that our vibronic coupling model [Phys. Rev. A 105, 062821 (2022)] captures primary features of the multidimensional dynamics of the temporary molecular anion, and the calculations qualitatively reproduce the spectrum. The shape of the spectrum is given by two overlapping contributions that originate in excitation of vibrational states within ${\mathrm{\ensuremath{\Sigma}}}_{g}^{+}$ and ${\mathrm{\ensuremath{\Pi}}}_{u}$ Fermi polyads. Propensity rules in terms of scattered and vibrational wave functions are also discussed to clarify the selectivity of states from the vibrational pseudocontinuum that is responsible for the observed fine structure.

Construction of Beyond Born-Oppenheimer Based Diabatic Surfaces and Generation of Photoabsorption Spectra: The Touchstone Pyrazine (C4 N2 H4 ).

We construct theoretically "exact" and numerically "accurate" Beyond Born-Oppenheimer (BBO) based diabatic potential energy surfaces (PESs) of pyrazine (C4 N2 H4 ) molecule involving lowest four excited adiabatic PESs (S1 to S4 ) and nonadiabatic coupling terms (NACTs) among those surfaces as functions of nonadiabatically active normal modes (Q1 , Q6a , Q9a and Q10a ) to compute its photoabsorption (PA) spectra. Those adiabatic PESs are calculated using CASSCF as well as MRCI based methodologies, where NACTs are obtained from CP-MCSCF approach. Employing ab initio quantities (adiabatic PESs and NACTs), it is possible to depict the conical intersections (CIs) and develop matrices of diabatic PESs over six normal mode planes. Once single-valued, smooth, symmetric and continuous 2×2 and 4×4 diabatic surface matrices are in hand for the first time, such matrices are used to perform multi-state multi-mode nuclear dynamics with the aid of Time-Dependent Discrete Variable Representation (TDDVR) methodology initializing the product type wavefunction on (S1 ) and (S2 ) states to obtain the corresponding PA spectra. TDDVR calculated spectra for those states (S1 and S2 ) obtained from BBO based 2×2 and 4×4 diabatic surface matrices show good and better agreement with the experimental results, respectively. Both of these calculated results depict better peak progression over the existing profiles of Multi-Configuration Time-Dependent Hartree (MCTDH) dynamics over 2×2 Vibronic Coupling Model (VCM) Hamiltonian.

Protocols for the in-silico screening of the perceived color of industrial dyes: Anthraquinones and indigos as study cases

Being able to accurately predict the perceived color of industrially relevant dyes is still a challenge for current theoretical approaches, especially when aiming at keeping a low computational cost and user input. In this work focusing on time-dependent density functional theory as the underlying quantum method of choice, we have analyzed the performance of different integrated approaches to reproduce the perceived color of selected anthraquinones and indigo derivatives. To reconstruct absorption spectra band shape both fully empirical approaches and different vibronic coupling approaches have been tested. The results obtained show how empirical broadening methods, and particularly the pseudo-Voigt model, provide good agreement with experiment in perceived colors but need parametrization and validation for each family of molecules. On the other hand, vibronic coupling models represent a possible powerful alternative, especially in the case of simplified models, though intrinsically limited to rigid dyes.

Multi-state and Multi-mode Vibronic Coupling Effects in the Photoionization Spectroscopy of Acetaldehyde.

Multi-state and multi-mode vibronic dynamics in the seven energetically low-lying (, , , , , , and ) electronic states of the acetaldehyde radical cation is theoretically studied in this article. Adiabatic energies of these electronic states are calculated by ab initio quantum chemistry methods. A vibronic coupling model of seven electronic states is constructed in a diabatic electronic basis to carry out the first-principles nuclear dynamics study. The vibronic spectrum is calculated and compared with the experimental findings reported in the literature. The progressions of vibrational modes found in the spectrum are assigned. The findings reveal that the and electronic states are energetically well-separated from the other electronic states and the remaining states ( to ) are energetically very close or even quasi-degenerate at the equilibrium geometry of the reference electronic ground state of acetaldehyde. The energetic proximity of to electronic states results in multiple multi-state conical intersections. The impact of electronic nonadiabatic interactions due to conical intersections on the vibronic structure of the photoionization band and nonradiative internal conversion dynamics is discussed.

Nonadiabatic Vibrational Resonance Raman Spectra from Quantum Dynamics Propagations with LVC Models. Application to Thymine.

We present a viable protocol to compute vibrational resonance Raman (vRR) spectra for systems with several close-lying and potentially coupled electronic states. It is based on the parametrization of linear vibronic coupling (LVC) models from time-dependent density functional theory (TD-DFT) calculations and quantum dynamics propagations of vibronic wavepackets with the multilayer version of the multiconfiguration time-dependent Hartree (ML-MCTDH) method. Our approach is applied to thymine considering seven coupled electronic states, comprising the three lowest bright states, and all vibrational coordinates. Computed vRR at different excitation wavelengths are in good agreement with the available experimental data. Up to 250 nm the signal is dominated by the lowest HOMO → LUMO transition, whereas at 233 nm, in the valley between the two lowest energy absorption bands, the contributions of all the three bright states, and their interferences and couplings, are important. Inclusion of solvent (water) effects improves the agreement with experiment, reproducing the coalescence of vibrational bands due to CC and C=O stretchings. With our approach we disentangle and assess the effect of interferences between the contribution of different quasi-resonant states to the transition polarizability and the effect of interstate couplings. Our findings strongly suggest that in cases of close-lying and potentially coupled states a simple inclusion of interference effects is not sufficient, and a fully nonadiabatic computation should instead be performed. We also document that for systems with strong couplings and quasi-degenerate states, the use of HT perturbative approach, not designed for these cases, may lead to large artifacts.

Velocity map imaging spectroscopy of C2H- and C2D-: A benchmark study of vibronic coupling interactions.

High-resolution velocity-map imaged photoelectron spectra of the ethynyl anions C2H- and C2D- are measured at photon wavelengths between 355 and 266 nm to investigate the complex interactions between the closely lying X̃2Σ+ and Ã2Π electronic states. An indicative kinetic energy resolution of 0.4%, together with the full angular dependence of the fast electrons, provides a detailed description of the vibronically coupled structure. It is demonstrated that a modest quadratic vibronic coupling model, parameterized by the quasidiabatic ansatz, is sufficient to accurately recreate all the observed vibronic interactions. Simulated spectra are shown to be in excellent agreement with the experimental data, verifying the proposed model and providing a framework that may be used to accurately simulate spectra of larger C2nH monohydride carbon chains. New spectral assignments are supported by experimental electron anisotropy measurements and Dyson orbital calculations.

Controlling energy conservation in quantum dynamics with independently moving basis functions: Application to multi-configuration Ehrenfest.

Application of the time-dependent variational principle to a linear combination of frozen-width Gaussians describing the nuclear wavefunction provides a formalism where the total energy is conserved. The computational downside of this formalism is that trajectories of individual Gaussians are solutions of a coupled system of differential equations, limiting implementation to serial propagation algorithms. To allow for parallelization and acceleration of the computation, independent trajectories based on simplified equations of motion were suggested. Unfortunately, within practical realizations involving finite Gaussian bases, this simplification leads to breaking the energy conservation. We offer a solution for this problem by using Lagrange multipliers to ensure the energy and normconservation regardless of basis function trajectories or basis completeness. We illustrate our approach within the multi-configurational Ehrenfest method considering a linear vibronic coupling model.