Comparative Study of Uracil Excited-State Photophysics in Water and Acetonitrile via RMS-CASPT2-Driven Quantum-Classical Trajectories

We present a nonadiabatic molecular dynamics study of the ultrafast processes occurring in uracil upon UV light absorption, leading to electronic excitation and subsequent nonradiative decay. Previous studies have indicated that the mechanistic details of this process are drastically different depending on whether the process takes place in the gas phase, acetonitrile, or water. However, such results have been produced using quantum chemical methods that did not incorporate both static and dynamic electron correlation. In order to assess the previously proposed mechanisms, we simulate the photodynamics of uracil in the three environments mentioned above using quantum-classical trajectories and, for solvated uracil, hybrid quantum mechanics/molecular mechanics (QM/MM) models driven by the rotated multistate complete active space second-order perturbation (RMS-CASPT2) method. To do so, we exploit the gradient recently made available in OpenMolcas and compare the results to those obtained using the complete active space self-consistent field (CASSCF) method only accounting for static electron correlation. We show that RMS-CASPT2 produces, in general, a mechanistic picture different from the one obtained at the CASSCF level but confirms the hypothesis advanced on the basis of previous ROKS and TDDFT studies thus highlighting the importance of incorporating dynamic electron correlation in the investigation of ultrafast electronic deactivation processes.

Robust ΔSCF calculations with direct energy functional minimization methods and STEP for molecules and materials

The direct energy functional minimization method using the orbital transformation (OT) scheme in the program package CP2K has been employed for Δ self-consistent field (ΔSCF) calculations. The OT method for non-uniform molecular orbitals occupations allows us to apply the ΔSCF method for various kinds of molecules and periodic systems. Vertical excitation energies of heteroaromatic molecules and condensed phase systems, such as solvated ethylene and solvated uracil obeying periodic boundary conditions, are reported using the ΔSCF method. In addition, a Re-phosphate molecule attached to the surface of anatase (TiO₂) has been investigated. Additionally, we have implemented a recently proposed state-targeted energy projection ΔSCF algorithm [K. Carter-Fenk and J. M. Herbert, J. Chem. Theory Comput. 16(8), 5067-5082 (2020)] for diagonalization based SCF in CP2K. It is found that the OT scheme provides a smooth and robust SCF convergence for all investigated excitation energies and (non-)periodic systems.

Coupling 6-chloro-3-methyluracil with copper: structural features, theoretical analysis, and biofunctional properties

As nucleobases in RNA and DNA, uracil and 5-methyluracil represent a recognized class of bioactive molecules and versatile ligands for coordination compounds with various biofunctional properties. In this study, 6-chloro-3-methyluracil (Hcmu) was used as an unexplored building block for the self-assembly generation of a new bioactive copper(II) complex, [Cu(cmu)₂(H₂O)₂]·4H₂O (1). This compound was isolated as a stable crystalline solid and fully characterized in solution and solid state by a variety of spectroscopic methods (UV-vis, EPR, fluorescence spectroscopy), cyclic voltammetry, X-ray diffraction, and DFT calculations. The structural, topological, H-bonding, and Hirshfeld surface features of 1 were also analyzed in detail. The compound 1 shows a distorted octahedral {CuN₂O₄} coordination environment with two trans cmu^- ligands adopting a bidentate N,O-coordination mode. The monocopper(II) molecular units participate in strong H-bonding interactions with water molecules of crystallization, leading to structural 0D → 3D extension into a 3D H-bonded network with a tfz-d topology. Molecular docking and ADME analysis as well as antibacterial and antioxidant activity studies were performed to assess the bioactivity of 1. In particular, this compound exhibits a prominent antibacterial effect against Gram negative (E. coli, P. aeruginosa) and positive (S. aureus, B. cereus) bacteria. The obtained copper(II) complex also represents the first structurally characterized coordination compound derived from 6-chloro-3-methyluracil, thus introducing this bioactive building block into a family of uracil metal complexes with notable biofunctional properties.

IR-VUV spectroscopy of pyridine dimers, trimers and pyridine-ammonia complexes in a supersonic jet

The infrared spectra of the C-H stretching vibrations of (pyridine)m, m = 1-3, and the N-H stretching vibrations of (pyridine)m-(NH3)n, m = 1, 2; n = 1-4, complexes were investigated by infrared (IR)-vacuum ultraviolet (VUV) spectroscopy under jet-cooled conditions. The ionization potential (IP0) of the pyridine monomer was determined to be 74 546 cm-1 (9.242 eV), while its complexes showed only smooth curves of the ionization thresholds at ∼9 eV, indicating large structural changes in the ionic form. The pyridine monomer exhibits five main features with several satellite bands in the C-H stretching region at 3000-3200 cm-1. Anharmonic calculations including Fermi-resonance were carried out to analyze the candidates of the overtone and combination bands which can couple to the C-H stretching fundamentals. For (pyridine)2 and (pyridine)3, most C-H bands are blue-shifted by 3-5 cm-1 from those of the monomer. The structures revealed by random searching algorithms with density functional methods indicate that the π-stacked structure is most stable for (pyridine)2, while (pyridine)3 prefers the structures stabilized by dipole-dipole and C-Hπ interactions. For the (pyridine)m-(NH3)n complexes, the mass spectrum exhibited a wide range distribution of the complexes. The observed IR spectra in the N-H stretching vibrations of the complexes showed four main bands in the 3200-3450 cm-1 region. These features are very similar to those of (NH3)n complexes, and the bands are assigned to the anti-symmetric N-H stretching band (ν3), the symmetric N-H stretching (ν1) band, and the first overtone bands of the N-H bending vibrations (2ν4). The anharmonic calculations including the Fermi-resonance between ν1 and 2ν4 well reproduced the observed spectra.

How intermolecular interactions influence electronic absorption spectra: insights from the molecular packing of uracil in condensed phases

Intermolecular interactions in terms of molecular packing are crucial for the investigation of the absorption spectra of uracil in different environments.

Solvent effects on excitation energies obtained using the state-specific TD-DFT method with a polarizable continuum model based on constrained equilibrium thermodynamics

Nonequilibrium solvation effects need to be treated properly in the study of electronic absorption processes of solutes since solvent polarization is not in equilibrium with the excited-state charge density of the solute. In this work, we developed a state specific (SS) method based on the novel nonequilibrium solvation model with constrained equilibrium manipulation to account for solvation effects in electronic absorption processes. Time-dependent density functional theory (TD-DFT) is adopted to calculate electronic excitation energies and a polarizable continuum model is employed in the treatment of bulk solvent effects on both the ground and excited electronic states. The equations based on this novel nonequilibrium solvation model in the framework of TDDFT to calculate vertical excitation energy are presented and implemented in the Q-Chem package. The implementation is validated by comparing reorganization energies for charge transfer excitations between two atoms obtained from Q-Chem and those obtained using a two-sphere model. Solvent effects on electronic transitions of coumarin 153 (C153), acetone, pyridine, (2E)-3-(3,4-dimethoxyphenyl)-1-(2-hydroxyphenyl)prop-2-en-1-one (DMHP), and uracil in different solvents are investigated using the newly developed code. Our results show that the obtained vertical excitation energies as well as spectral shifts generally agree better with the available experimental values than those obtained using the traditional nonequlibrium solvation model. This new model is thus appropriate to study nonequilibrium excitation processes in solution.

Electrostatic Interaction-Induced Room-Temperature Phosphorescence in Pure Organic Molecules from QM/MM Calculations

Room temperature phosphorescence (RTP) from pure organic material is rare due to the low phosphorescence quantum efficiency. That is why the recent discovery of crystallization induced RTP for several organic molecules aroused strong interests. Through a combined quantum and molecular mechanics CASPT2/AMBER scheme taking terephthalic acid (TPA) as example, we found that electrostatic interaction not only can induce an enhanced radiative decay T1 → S0 through the dipole-allowed S1 intermediate state, but also can hinder the nonradiative decay process upon crystallization. From gas phase to crystal, the nature of S1 state is converted to (1)(π,π*) from (1)(n,π*) character, enhancing transition dipole moment and serving as an efficient intermediate radiative pathway for T1 → S0 transition, and eventually leading to a boosted RTP. The intermolecular packing also blocks the nonradiative decay channel of the high-frequency C═O stretching vibration with large vibronic coupling, rather than the conventional low-frequency aromatic rotation in crystal. This mechanism also holds for other organic compounds that contain both ketones and aromatic rings.

Quantum Mechanical Studies on the Photophysics and the Photochemistry of Nucleic Acids and Nucleobases

The photophysics and photochemistry of DNA is of great importance due to the potential damage of the genetic code by UV light. Quantum mechanical studies have played a key role in interpretating the results of modern time-resolved pump-probe spectroscopy, and in elucidating the main photoactivated reactive paths. This review provides a concise, complete picture of the computational studies carried out, approximately, in the past decade. We start with an overview of the photophysics of the nucleobases in the gas phase and in solution. We discuss the proposed mechanisms for ultrafast decay to the ground state, that involve conical intersections, consider the role of triplet states, and analyze how the solvent modulates the photophysics. Then we move to larger systems, from dinucleotides to single- and double-stranded oligonucleotides. We focus on the possible role of charge transfer and delocalized or excitonic states in the photophysics of these systems and discuss the main photochemical paths. We finish with an outlook on the current challenges in the field and future directions of research.

Solvent Induced Transformations of n-π* Absorption in Formaldehyde, Acetaldehyde, and Acetone

Absorption spectra of formaldehyde (FA), acetaldehyde (AA), and acetone are compared in the vapor phase, nonpolar, and polar solutions at 295 K. The vibronic n-π* transition of carbonyl chromophore is mainly composed of the overtones of >C═O stretching vibration. A new phenomenon is observed in liquid solutions, consisting of a relative increase of Franck-Condon factors for the second and third harmonics in FA, and the second to fourth replica in AA, with respect to the gas phase. In AA and acetone with poorly resolved vibronic structure, the redistribution of intensities produces a false "solvent shift" of the band maximum between the vapor and nonpolar liquid phase by -250 ± 50 cm(-1). Modification in vibronic coupling can also explain unusual narrowing of the band contour in the solution, reported earlier for acetone (Renge , I. J. Phys. Chem. A 2009, 113, 10678). No detectable shift occurs as a function of solvent polarizability (refractive index function (n(2) - 1)/(n(2) + 2)) in n-alkanes for FA, AA, and acetone, as well as for cyclopentanone and camphor. Incidentally, the bathochromic dispersive shift is almost exactly compensated by a hypsochromic induction shift. The latter is due to the diminishing dipole moment in the excited state of the carbonyl chromophore. Differences in polarizability α and dipole moments μ were estimated for FA (Δα = 0.33 ± 0.1 Å(3)), AA (Δμ = -1.05 ± 0.2 D, Δα = 0.5 ± 0.2 Å(3)), and acetone (Δμ = -1.3 ± 0.2 D, Δα = 0.65 ± 0.2 Å(3)). The increase of α by ∼10% upon excitation is plausible for a weak n-π* transition. By contrast, near doubling of α in the upper state has been reported recently for several ketones, with Δα reaching 10 Å(3) (Catalán, J.; Catalán, J. P. Phys. Chem. Chem. Phys. 2011, 13, 4072). Empirical partitioning of solvent shifts into repulsive-dispersive, induction, dipole-dipole, and hydrogen bonding contributions was proposed to serve as a benchmark in computer chemical calculations.

Solvation effects on the ground and excited states of p–n diblock-conjugated polymers

A semiempirical study was carried out to give theoretical insights of the solvation effects on the absorption and emission wavelength of oxadiazole-containing p–n diblock-conjugated polymers at both ground and excited states on the basis of experimental results. The Conductor-like Screening Model was used to calculate the absorption and emission wavelengths of the carefully designed model molecules in solvents with different dielectric constants (approximately ∊ = 1–35) from the optimized ground states and the lowest excited states by coupling AM1 formalism to the new configuration interaction scheme. The calculated results agree well with the experimental observations of the corresponding polymers. The geometrical modifications, charge changes upon photoexcitation, frontier orbital distributions, and dipoles of the two model molecules were investigated and discussed to explain the different solvation behavior of the corresponding diblock copolymers. The different charge-transfer behavior was suggested to be the major factor that causes the varied solvation effects of the p–n diblock molecular systems at excited states upon photoluminescence. These findings may provide important insights at the molecular level for selectively tuning the solvatochromic properties of p–n diblock-conjugated polymers for practical applications.

Theoretical calculation of reorganization energy for electron self-exchange reaction by constrained density functional theory and constrained equilibrium thermodynamics

Within the framework of constrained density functional theory (CDFT), the diabatic or charge localized states of electron transfer (ET) have been constructed. Based on the diabatic states, inner reorganization energy λin has been directly calculated. For solvent reorganization energy λs, a novel and reasonable nonequilibrium solvation model is established by introducing a constrained equilibrium manipulation, and a new expression of λs has been formulated. It is found that λs is actually the cost of maintaining the residual polarization, which equilibrates with the extra electric field. On the basis of diabatic states constructed by CDFT, a numerical algorithm using the new formulations with the dielectric polarizable continuum model (D-PCM) has been implemented. As typical test cases, self-exchange ET reactions between tetracyanoethylene (TCNE) and tetrathiafulvalene (TTF) and their corresponding ionic radicals in acetonitrile are investigated. The calculated reorganization energies λ are 7293 cm(-1) for TCNE/TCNE(-) and 5939 cm(-1) for TTF/TTF(+) reactions, agreeing well with available experimental results of 7250 cm(-1) and 5810 cm(-1), respectively.