Carlos M. Estévez

Ultrafast Ring-Opening/Closing and Deactivation Channels for a Model Spiropyran–Merocyanine System

The photochemistry of a model merocyanine-spiropyran system was analyzed theoretically at the MS-CASPT2//SA-CASSCF(14,12) level. Several excited singlet states were studied in both the closed spiropyran and open merocyanine forms, and the paths to the different S(1)/S(0) conical intersections found were analyzed. After absorption of UV light from the spiropyran form, there are two possible ultrafast routes to efficient conversion to the ground state; one involves the rupture of the C(spiro)-O bond leading to the open form and the other involves the lengthening of the C(spiro)-N bond with no photoreaction. From the merocyanine side the excited state can reach a very broad S(1)/S(0) conical intersection region that leads the system to the closed form after rotation of the central methine bond. Alternatively, rotation of the other methine bonds connects the system through different S(1)/S(0) conical intersections to several merocyanine isomers. The present work provides a theoretical framework for the recent experimental results (Buback , J. J. Am. Chem. Soc. 2010, 132, 1610-1619) and sheds light on the complex photochemistry of these kinds of compounds.

Product energy distributions for the four-center HF elimination from 1,1-difluoroethylene. A direct dynamics study

Product energy distributions (PEDs) were computed on the four-center HF elimination from 1,1-difluoroethylene by using direct trajectory calculations. The vibrational and rotational populations of HF obtained with a quasi-classical normal mode/rigid rotor excitation model compare very well with the experimental results. Also, the translational energy distributions obtained with an efficient microcanonical sampling (EMS) at the barrier are in excellent accord with experiment and do not substantially change as the excitation energy increases.

Anion-π Aromatic Neutral Tweezers Complexes: Are They Stable in Polar Solvents?

The impact of the solvent environment on the stabilization of the complexes formed by fluorine (T-F) and cyanide (T-CN) substituted tweezers with halide anions has been investigated theoretically. The study was carried out using computational methodologies based on density functional theory (DFT) and symmetry adapted perturbation theory (SAPT). Interaction energies were obtained at the M05-2X/6-31+G* level. The obtained results show a large stability of the complexes in solvents with large dielectric constant and prove the suitability of these molecular tweezers as potential hosts for anion recognition in solution. A detailed analysis of the effects of the solvent on the electron withdrawing ability of the substituents and its influence on the complex stability has been performed. In particular, the interaction energy in solution was split up into intermonomer and solvent-complex terms. In turn, the intermonomer interaction energy was partitioned into electrostatic, exchange, and polarization terms. Polar resonance structures in T-CN complexes are favored by polar solvents, giving rise to a stabilization of the intermonomer interaction, the opposite is found for T-F complexes. The solvent-complex energy increases with the polarity of the solvent in T-CN complexes, nonetheless the energy reaches a maximum and then decreases slowly in T-F complexes. An electron density analysis was also performed before and after complexation, providing an explanation to the trends followed by the interaction energies and their different components in solution.

A computational study of the protonation of simple amines in water clusters.

The microsolvation study of a group of amines with a variable number of water molecules was performed by conducting a theoretical analysis of the properties of the clusters formed by the amines with up to seven molecules of water. We describe the microsolvation of several amines focusing on the dissociation of a water molecule that transfers a proton to the amine and forms a hydroxide ion. Ab initio calculations were performed on these clusters employing the DFT/B3LYP and MP2 methods with the 6-311++G(2d,p) basis set. Several stationary points for each cluster were thus located and characterized as minima from frequency calculations. Intermolecular BSSE corrected interaction energies were obtained. The protonation mechanism of the amines was examined in terms of some parameters that include the lengths of the bonds involved in the process of proton transfer and the frequencies associated with certain O-H and N-H stretching modes. On the basis of the calculations, all studied amines present similar behavior but trimethylamine, whose limitations to be integrated in the water hydrogen bond network cause the instability of some of their complexes. The cyclic configurations are the most stable structures up to five water molecules due to the presence of cooperative effects associated with the hydrogen bonds of water molecules. However, when the number of water molecules increases the spatial forms become the most stable configurations. The dissociated forms were not found to have the most stable configuration in any of the studied systems but energetic differences between the dissociated and non-dissociated forms decrease with the number of water molecules.

A theoretical study of the low-lying electronic states of SC3

Abstract A theoretical study of the excited electronic states of the SC3 system has been carried out by means of density functional (B3LYP) and ab inito methods (QCI, MRCI, RHF-CC) in combination with relatively large basis sets. The structures of several singlet and triplet species are reported, together with their vibrational frequencies and infrared intensities. The electronic structures are analyzed by means of the natural bond orbital (NBO) method. The vertical electronic excitation energies and oscillator strengths of the ground state SCCC ( 1 Σ + ) have also been computed. The lowest-energy allowed transitions are X 1 Σ + → C 1 Π and X 1 Σ + → D 1 Σ + , and should have comparable intensities.

Influence of Multiple Conformations and Paths on Rate Constants and Product Branching Ratios. Thermal Decomposition of 1-Propanol Radicals

The potential energy surface involved in the thermal decomposition of 1-propanol radicals was investigated in detail using automated codes (tsscds2018 and Q2DTor). From the predicted elementary reactions, a relevant reaction network was constructed to study the decomposition at temperatures in the range 1000-2000 K. Specifically, this relevant network comprises 18 conformational reaction channels (CRCs), which in general exhibit a large wealth of conformers of reactants and transition states. Rate constants for all the CRCs were calculated using two approaches within the formulation of variational transition-state theory (VTST), as incorporated in the TheRa program. The simplest, one-well (1W) approach considers only the most stable conformer of the reactant and that of the transition state. In the second, more accurate approach, contributions from all the reactant and transition-state conformers are taken into account using the multipath (MP) formulation of VTST. In addition, kinetic Monte Carlo (KMC) simulations were performed to compute product branching ratios. The results show significant differences between the values of the rate constants calculated with the two VTST approaches. In addition, the KMC simulations carried out with the two sets of rate constants indicate that, depending on the radical considered as reactant, the 1W and the MP approaches may display different qualitative pictures of the whole decomposition process.