6533b834fe1ef96bd129d680

RESEARCH PRODUCT

Assessment of the Potential Energy Hypersurfaces in Thymine within Multiconfigurational Theory: CASSCF vs. CASPT2

Manuela MerchánAntonio Francés-monerrisJavier Segarra‐martíDaniel Roca-sanjuán

subject

Chemistry Multidisciplinary2-DIMENSIONAL ELECTRONIC SPECTROSCOPYPharmaceutical Sciencephotostability0305 Organic Chemistry01 natural sciencesLOWEST TRIPLET-STATEAnalytical ChemistryInterpretation (model theory)Molecular dynamicschemistry.chemical_compoundComputational chemistryDrug DiscoveryComplete active spacePerturbation theoryRETINAL CHROMOPHORE MODELComputingMilieux_MISCELLANEOUSAB-INITIOphotochemistry010304 chemical physicsBasis (linear algebra)ChemistryCOUPLED-CLUSTER METHODSPhotochemical ProcessesPotential energy[CHIM.THEO]Chemical Sciences/Theoretical and/or physical chemistryChemistryChemistry (miscellaneous)2ND-ORDER PERTURBATION-THEORYPhysical SciencesANO BASIS-SETSCASSCF/CASPT2Molecular MedicineThermodynamicsLife Sciences & BiomedicineBiochemistry & Molecular BiologyField (physics)INITIO MOLECULAR-DYNAMICSMolecular Dynamics Simulation010402 general chemistryMolecular physicsArticlelcsh:QD241-441lcsh:Organic chemistryCASSCF/CASPT2; photochemistry; DNA; thymine; photostability0103 physical sciencesthyminePhysical and Theoretical ChemistryULTRAFAST INTERNAL-CONVERSIONScience & TechnologyOrganic ChemistryDNAEXCITED-STATE DYNAMICS0104 chemical sciencesThymineModels Chemical

description

The present study provides new insights into the topography of the potential energy hypersurfaces (PEHs) of the thymine nucleobase in order to rationalize its main ultrafast photochemical decay paths by employing two methodologies based on the complete active space self-consistent field (CASSCF) and the complete active space second-order perturbation theory (CASPT2) methods: (i) CASSCF optimized structures and energies corrected with the CASPT2 method at the CASSCF geometries and (ii) CASPT2 optimized geometries and energies. A direct comparison between these strategies is drawn, yielding qualitatively similar results within a static framework. A number of analyses are performed to assess the accuracy of these different computational strategies under study based on a variety of numerical thresholds and optimization methods. Several basis sets and active spaces have also been calibrated to understand to what extent they can influence the resulting geometries and subsequent interpretation of the photochemical decay channels. The study shows small discrepancies between CASSCF and CASPT2 PEHs, displaying a shallow planar or twisted 1(ππ*) minimum, respectively, and thus featuring a qualitatively similar scenario for supporting the ultrafast bi-exponential deactivation registered in thymine upon UV-light exposure. A deeper knowledge of the PEHs at different levels of theory provides useful insight into its correct characterization and subsequent interpretation of the experimental observations. The discrepancies displayed by the different methods studied here are then discussed and framed within their potential consequences in on-the-fly non-adiabatic molecular dynamics simulations, where qualitatively diverse outcomes are expected.

10.3390/molecules21121666https://dx.doi.org/10.3390/molecules21121666