A generalised vibronic-coupling Hamiltonian model for benzopyran (original) (raw)
Related papers
The Journal of Chemical Physics
We present a model for the lowest two potential energy surfaces (PESs) that describe the photoinduced ring-opening reaction of benzopyran taken as a model compound to study the photochromic ring-opening reaction of indolinobenzospiropyran and its evolution toward its open-chain analog. The PESs are expressed in terms of three effective rectilinear coordinates. One corresponds to the direction between the equilibrium geometry in the electronic ground state, referred to as the Franck-Condon geometry, and the minimum of conical intersection (CI), while the other two span the two-dimensional branching space at the CI. The model correctly reproduces the topography of the PESs. The ab initio calculations are performed with the extended multiconfiguration quasidegenerate perturbation theory at second order method. We demonstrate that accounting for electron dynamic correlation drastically changes the global energy landscape since some zwitterionic states become strongly stabilized. Quantum dynamics calculations using this PES model produce an absorption spectrum that matches the experimental one to a good accuracy.
A quantum dynamics study of the benzopyran ring opening guided by laser pulses
Chemical Physics, 2014
The ring-opening photoisomerization of benzopyran, which occurs via a photochemical route involving a conical intersection, has been studied with quantum dynamics calculations using the multi-configuration time-dependent Hartree method (MCTDH). We introduce a mechanistic strategy to control the conversion of benzopyran to merocyanine with laser pulses. We use a six-dimensional model developed in a previous work for the potential energy surfaces (PES) based on an extension of the vibronic-coupling Hamiltonian model (diabatization method by ansatz), which depends on the most active degrees of freedom. The main objective of these quantum dynamics simulations is to provide a set of strategies that could help experimentalists to control the photoreactivity vs. photostability ratio (selectivity). In this work we present: (i) a pump-dump technique used to control the photostability, (ii) a two-step strategy to enhance the reactivity of the system: first, a pure vibrational excitation in the electronic ground state that prepares the system and, second, an ultraviolet excitation that brings the system to the first adiabatic electronic state; (iii) finally the effect of a non-resonant pulse (Stark effect) on the dynamics.
Quantum dynamics of the photostability of pyrazine
Physical chemistry chemical physics : PCCP, 2015
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.
A generalised 17-state vibronic-coupling Hamiltonian model for ethylene
The Journal of Chemical Physics, 2012
In a previous work [B. Lasorne, M. A. Robb, H.-D. Meyer, and F. Gatti, "The electronic excited states of ethylene with large-amplitude deformations: A dynamical symmetry group investigation," Chem. Phys. 377, 30-45 (2010); 382, 132 (2011) (Erratum)], we investigated the electronic structure of ethylene (ethene, C 2 H 4 ) in terms of 17 dominant configurations selected at the multiconfiguration self-consistent field level of theory. These were shown to be sufficient to recover most of the static electron correlation among the first valence and Rydberg states at all geometries. We also devised a strategy to build a 17-quasidiabatic-state matrix representation of the electronic Hamiltonian for curvilinear coordinates using dynamical symmetry. Here, we present fitted surfaces in the form of a generalised vibronic-coupling Hamiltonian model for two nuclear coordinates, CC bond stretching and torsion. Dynamic electron correlation is included into the electronic structure to improve the energetics of the Rydberg states at the multireference configuration interaction level of theory. The chemical interpretation of the adiabatic states of interest does not change qualitatively, which validates our choice of underlying quasidiabatic states in the model. The absorption spectrum is calculated with quantum dynamics and partially assigned. This first two-dimensional model shows a surprisingly good agreement with the experimental spectrum.
Trajectory-based approaches to vibronic spectroscopy of molecular systems
2018
The goal of this thesis is to bridge the gap between the two standard theoretical approaches to vibronic spectroscopy via trajectory-based methods. As the starting point, a generalized time-correlation function is introduced and the ring-polymer molecular dynamics method is generalized to vibronic transitions, yielding an improvement over known classical approximations. Further, the vibronic spectrum is evaluated via the Matsubara dynamics. Employing an ad-hoc modification, the sign problem of the Matsubara method can be circumvented and the spectra of model systems are adequately simulated.
Molecular dynamics with restrictions derived from optical spectra
Journal of Computational Chemistry, 2009
The information about molecular structure coded in the optical spectra must often be deciphered by complicated computational procedures. A combination of spectral modeling with the molecular dynamic simulations makes the process simpler, by implicit accounting for the inhomogeneous band broadening and Boltzmann averaging of many conformations. Ideally, geometries of studied systems can be deduced by a direct confrontation of such modeling with the experiment. In this work, the comparison is enhanced by restrictions to molecular dynamics propagations based on the Raman and Raman optical activity spectra. The methodology is introduced and tested on model systems comprising idealized H2O2, H2O3 molecules, and the alanine zwitterion. An additional gradient term based on the spectral overlap smoothed by Fourier transformation is constructed and added to the molecular energy during the molecular dynamics run. For systems with one prevalent conformation the method did allow to enrich the Boltzmann ensemble by a spectroscopically favored structure. For systems with multiconformational equilibria families preferential conformations can be selected. An alternative algorithm based on the comparison of the averaged spectra with the reference enabling iterative updates of the conformer probabilities provided even more distinct distributions in shorter times. It also accounts for multiconformer equilibria and provided realistic spectra and conformer distribution for the alanine. © 2008 Wiley Periodicals, Inc. J Comput Chem, 2009
2008
Ab inito molecular dynamics, although still challenge, is becoming an available tool for the investigation of the photodynamics of aromatic heterocyclic systems. Potential energy surfaces and dynamics simulations for three particular examples and different aspects of the excited and ground state dynamics are presented and discussed. Aminopyrimidine is investigated as a model for adenine. It shows ultrafast S 1 -S 0 decay in about 400 fs. The inclusion of mass-restrictions to emulate the imidizole group increases the lifetime to about 950 fs, a value similar to the lifetime of adenine. The S 2 -S 1 deactivation, typical in the fast component of the decay of nucleobases, is investigated in pyridone. In this case, the S 2 -state lifetime is 52 fs. The hot ground-state dynamics of pyrrole starting at the puckered conical intersection is shown to produce ring-opened structures consistent with the experimental results
Review Nonadiabatic Molecular Dynamics Based on Trajectories
2013
Performing molecular dynamics in electronically excited states requires the inclusion of nonadiabatic effects to properly describe phenomena beyond the Born-Oppenheimer approximation. This article provides a survey of selected nonadiabatic methods based on quantum or classical trajectories. Among these techniques, trajectory surface hopping constitutes an interesting compromise between accuracy and efficiency for the simulation of medium-to large-scale molecular systems. This approach is, however, based on non-rigorous approximations that could compromise, in some cases, the correct description of the nonadiabatic effects under consideration and hamper a systematic improvement of the theory. With the help of an in principle exact description of nonadiabatic dynamics based on Bohmian quantum trajectories, we will investigate the origin of the main approximations in trajectory surface hopping and illustrate some of the limits of this approach by means of a few simple examples.
Nonadiabatic Excited-State Molecular Dynamics: Modeling Photophysics in Organic Conjugated Materials
Accounts of Chemical Research, 2014
CONSPECTUS: To design functional photoactive materials for a variety of technological applications, researchers need to understand their electronic properties in detail and have ways to control their photoinduced pathways. When excited by photons of light, organic conjugated materials (OCMs) show dynamics that are often characterized by large nonadiabatic (NA) couplings between multiple excited states through a breakdown of the Born−Oppenheimer (BO) approximation. Following photoexcitation, various nonradiative intraband relaxation pathways can lead to a number of complex processes. Therefore, computational simulation of nonadiabatic molecular dynamics is an indispensable tool for understanding complex photoinduced processes such as internal conversion, energy transfer, charge separation, and spatial localization of excitons. Over the years, we have developed a nonadiabatic excited-state molecular dynamics (NA-ESMD) framework that efficiently and accurately describes photoinduced phenomena in extended conjugated molecular systems. We use the fewest-switches surface hopping (FSSH) algorithm to treat quantum transitions among multiple adiabatic excited state potential energy surfaces (PESs). Extended molecular systems often contain hundreds of atoms and involve large densities of excited states that participate in the photoinduced dynamics. We can achieve an accurate description of the multiple excited states using the configuration interaction single (CIS) formalism with a semiempirical model Hamiltonian. Analytical techniques allow the trajectory to be propagated "on the fly" using the complete set of NA coupling terms and remove computational bottlenecks in the evaluation of excited-state gradients and NA couplings. Furthermore, the use of state-specific gradients for propagation of nuclei on the native excited-state PES eliminates the need for simplifications such as the classical path approximation (CPA), which only uses ground-state gradients. Thus, the NA-ESMD methodology offers a computationally tractable route for simulating hundreds of atoms on ∼10 ps time scales where multiple coupled excited states are involved. In this Account, we review recent developments in the NA-ESMD modeling of photoinduced dynamics in extended conjugated molecules involving multiple coupled electronic states. We have successfully applied the outlined NA-ESMD framework to study ultrafast conformational planarization in polyfluorenes where the rate of torsional relaxation can be controlled based on the initial excitation. With the addition of the state reassignment algorithm to identify instances of unavoided crossings between noninteracting PESs, NA-ESMD can now be used to study systems in which these so-called trivial unavoided crossings are expected to predominate. We employ this technique to analyze the energy transfer between poly(phenylene vinylene) (PPV) segments where conformational fluctuations give rise to numerous instances of unavoided crossings leading to multiple pathways and complex energy transfer dynamics that cannot be described using a simple Forster model. In addition, we have investigated the mechanism of ultrafast unidirectional energy transfer in dendrimers composed of poly(phenylene ethynylene) (PPE) chromophores and have demonstrated that differential nuclear motion favors downhill energy transfer in dendrimers. The use of native excited-state gradients allows us to observe this feature.
The Journal of Chemical Physics, 2010
We present a generalization to our previously developed quantum wavepacket ab initio molecular dynamics ͑QWAIMD͒ method by using multiple diabatic electronic reduced single particle density matrices, propagated within an extended Lagrangian paradigm. The Slater determinantal wavefunctions associated with the density matrices utilized may be orthogonal or nonorthogonal with respect to each other. This generalization directly results from an analysis of the variance in electronic structure with quantum nuclear degrees of freedom. The diabatic electronic states are treated here as classical parametric variables and propagated simultaneously along with the quantum wavepacket and classical nuclei. Each electronic density matrix is constrained to be N-representable. Consequently two sets of new methods are derived: extended Lagrangian-QWAIMD ͑xLag-QWAIMD͒ and diabatic extended Lagrangian-QWAIMD ͑DxLag-QWAIMD͒. In both cases, the instantaneous potential energy surface for the quantum nuclear degrees of freedom is constructed from the diabatic states using an on-the-fly nonorthogonal multireference formalism. By introducing generalized grid-based electronic basis functions, we eliminate the basis set dependence on the quantum nucleus. Subsequent reuse of the two-electron integrals during the on-the-fly potential energy surface computation stage yields a substantial reduction in computational costs. Specifically, both xLag-QWAIMD and DxLag-QWAIMD turn out to be about two orders of magnitude faster than our previously developed time-dependent deterministic sampling implementation of QWAIMD. Energy conservation properties, accuracy of the associated potential surfaces, and vibrational properties are analyzed for a family of hydrogen bonded systems.