Lutidyl Radical Photoelectron Spectra Reveal Additive Substituent Effects on Benzyl Derivatives’ Ionization Energy

Understanding how isomerism influences photoelectron spectra helps in the assignment and analysis of reactive mixtures, especially for heterocycles with numerous isomers. Threshold photoelectron spectra of lutidyl radical isomers, i. e., benzyl derivatives with a nitrogen heteroatom and a methyl substituent, are recorded using vacuum ultraviolet synchrotron radiation. The radicals are produced by flash pyrolysis from aminomethyl methylpyridine precursors. Experimental ionization energies are determined to be 7.54, 7.50, and 7.45 eV for 2,4-, 2,6- and 3,5-lutidyl, respectively, in excellent agreement with composite method calculations. Franck–Condon simulations aid the TPES assignment but are also shown to exhibit artifacts if large-amplitude motions, notably the methyl internal rotation are assumed to be active in the double harmonic approximation. Based on calculated adiabatic ionization energies (AIE) of benzyl, picolyl, and xylyl radicals, the N and CH₃ substituent effects are found to be additive, position-dependent and decrease in the para>ortho≳meta order in magnitude with the nitrogen heteroatom increasing and the methyl substituent decreasing the AIE. These effects are discussed in light of the charge distribution upon ionization. The additivity of the substituent effects also helps predict the influence of substituents on the binding energy of the unpaired electron in analogous radicals.     

Calculation of the Vibrationally Resolved,Circularly Polarized Luminescence of d‐Camphorquinone and (S,S)‐trans‐β‐Hydrindanone

Circularly polarized luminescence (CPL), the differential emission of left‐ and right‐handed circularly polarized light from a molecule, is modeled by using time‐dependent density functional theory. Calculations of the CPL spectra for the first electronic excited states of d‐camphorquinone and (S,S)‐trans‐β‐hydrindanone under the Franck–Condon approximation and using various functionals are presented, as well as calculations of absorption, emission, and circular dichroism spectra. The functionals B3LYP, BHLYP, and CAM‐B3LYP are employed, along with the TZVP and aug‐cc‐pVDZ Gaussian‐type basis sets. For the lowest‐energy transitions, all functionals and basis sets perform comparably, with the long‐range‐corrected CAM‐B3LYP better reproducing the excitation energy of camphorquinone but leading to a blue shift with respect to experiment for hydrindanone. The vibrationally resolved spectra of camphorquinone are very well reproduced in terms of peak location, widths, shapes, and intensities. The spectra of hydrindanone are well reproduced in terms of overall envelope shape and width, as well as the lack of prominent vibrational structure in the emission and CPL spectra. Overall the simulated spectra compare well with experiment, and reproduce the band shapes, emission red shifts, and presence or absence of visible vibrational fine structure.     

Difluorocarbene Studied with Threshold Photoelectron Spectroscopy (TPES): Measurement of the First Adiabatic Ionization Energy (AIE) of CF2

The first photoelectron band of difluorocarbene CF₂, has been studied by threshold photoelectron (TPE) spectroscopy. CF₂ was prepared by microwave discharge of a flowing mixture of hexafluoropropene, C₃F₆, and argon. A vibrationally resolved band was observed in which at least twenty‐two components were observed. In the first PE band of CF₂, the adiabatic ionization energy differs significantly from the vertical ionization energy because, for the ionization CF₂⁺ (X?²A₁)+e^? ← CF₂ (X?¹A₁), there is an increase in the FCF bond angle (by ≈20°) and a decrease in the C? F bond length (by ≈0.7 Å). The adiabatic component was not observed in the experimental TPE spectrum. However, on comparing this spectrum with an ab initio/Franck–Condon simulation of this band, using results from high‐level ab initio calculations, the structure associated with the vibrational components could be assigned. This led to alignment of the experimental TPE spectrum and the computed Franck–Condon envelope, and a determination of the first adiabatic ionization energy of CF₂ as (11.362±0.005) eV. From the assignment of the vibrational structure, values were obtained for the harmonic and fundamental frequencies of the symmetric stretching mode (ν₁′) and symmetric bending mode (ν₂′) in CF₂⁺ (X?²A₁).     

Valence Photoionization of Thymine: Ionization Energies,Vibrational Structure,and Fragmentation Pathways from the Slow to the Ultrafast

The photoionization of thymine has been studied by using vacuum ultraviolet radiation and imaging photoelectron photoion coincidence spectroscopy after aerosol flash vaporization and bulk evaporation. The two evaporation techniques have been evaluated by comparison of the photoelectron spectra and breakdown diagrams. The adiabatic ionization energies for the first four electronic states were determined to be 8.922±0.008, 9.851±0.008, 10.30±0.02, and 10.82±0.01 eV. Vibrational features have been assigned for the first three electronic states with the help of Franck–Condon factor calculations based on density functional theory and wave function theory vibrational analysis within the harmonic approximation. The breakdown diagram of thymine, as supported by composite method ab initio calculations, suggests that the main fragment ions are formed in sequential HNCO-, CO-, and H-loss dissociation steps from the thymine parent ion, with the first step corresponding to a retro-Diels–Alder reaction. The dissociation rate constants were extracted from the photoion time-of-flight distributions and used together with the breakdown curves to construct a statistical model to determine 0 K appearance energies of 11.15±0.16 and 11.95±0.09 eV for the m/z 83 and 55 fragment ions, respectively. These results have allowed us to revise previously proposed fragmentation mechanisms and to propose a model for the final, nonstatistical H-loss step in the breakdown diagram, yielding the m/z 54 fragment ion at an appearance energy of 13.24 eV.     

IR Low‐Temperature Matrix and ab Initio Study on β‐Alanine Conformers

For the first time the argon‐matrix low‐temperature FTIR spectra of β‐alanine are recorded. They reveal a quite complicated spectral pattern which suggests the presence of several β‐alanine conformers in the matrix. To interpret the spectra, the eighteen β‐alanine conformers, stable in the gas phase, are estimated at the B3LYP and MP2 levels combined with the aug‐cc‐pVDZ. Ten low‐energy structures are reoptimized at the QCISD/aug‐cc‐pVDZ and B3LYP and MP2 levels by using the aug‐cc‐pVTZ basis sets. Assignment of the experimental spectra is undertaken on the basis of the calculated B3LYP/aug‐cc‐pVDZ anharmonic IR frequencies as well as careful estimation of the conformer population. The presence of at least three β‐alanine conformers is demonstrated. The detailed analysis of IR spectra points to the possible presence of five additional β‐alanine conformers.     

Improving on equilibrium isotropic nuclear shielding constants

Isotropic nuclear shielding constants at the equilibrium molecular structure σ_eq and zero‐point vibrational corrections (ZPVCs) to σ_eq are evaluated using the B3LYP/aug‐cc‐pVTZ level of theory, as well as the KT2/aug‐cc‐pVTZ level of theory. Various scaling factors and systematic corrections are obtained by linear regression to experimental shielding constants. Comparisons of the scaled and systematically corrected equilibrium and vibrationally averaged shielding constants reveal that, at the 99% confidence level, the ZPVCs via second‐order perturbation theory do not improve the agreement of B3LYP/aug‐cc‐pVTZ and KT2/aug‐cc‐pVTZ calculated shielding constants with experiment. This holds true when the same analysis is applied to CCSD(T)/aug‐cc‐pCV[TQ]Z calculated σ_eq of Teale et al. [Journal of Chemical Physics 2013, 138, 024111]. In addition, at the 99% confidence level, B3LYP/aug‐cc‐pVTZ and KT2/aug‐cc‐pVTZ scaled and systematically corrected shielding constants are found to be statistically no different from CCSD(T)/aug‐cc‐pCV[TQ]Z calculated σ_eq. The use of scaling factors and systematic corrections could thus provide a cheaper but yet reasonably accurate alternative for the study of nuclear shielding constants of larger systems.     

Improved Ionization Energies for the Two Isomers of Phenylpropargyl Radical

The ionization of two resonantly stabilized radicals, namely 1‐phenylpropargyl (1PPR) and 3‐phenylpropargyl (3PPR) are reinvestigated applying vacuum ultraviolet synchrotron radiation and threshold photoelectron spectroscopy. Ionization energies of 7.24±0.02 and 7.25±0.01 eV are obtained for 1 and 3PPR respectively, which compare well with ab initio calculations. The quality of the spectra are significantly improved mostly due to the application of a new high‐photon‐flux grating available at the VUV‐beamline of Swiss Light Source. Resolved vibrational features are assigned according to a Franck–Condon approach.     

A theoretical study of the photoelectron spectra of dichloroketene with accurate computation of ionization energies via complete basis set limit extrapolation

We calculated the equilibrium geometries and harmonic vibrational frequencies of the ground state and five cationic states of dichloroketene using (TD-)B3LYP, PBE0, and M06/M06-2X approaches. The photoelectron spectra of dichloroketene were simulated by computing Franck-Condon factors. The ionization energies were computed using the CCSD(T) approach with extrapolation to the complete basis set (CBS) limit. We propose two new CBS energy formulas (E = E_CBS + Aexp(-x) + B/(x−1) ⁿ, n = 2 or 3) and compare the performance of different CBS approaches. A new ionic state of dichloroketene belonging to the C_s point group is reported. This state is identified as the first excited state of Cl₂CCO⁺ having a double-well potential-energy curve along the CCO bending mode with a barrier height of 1.335 eV. The simulated photoelectron spectra are in agreement with the experiment. The vertical ionization energies calculated via spectral simulation are more accurate compared with those obtained at the ground-state structure. Among the CBS formulas used, the proposed ansatz with n = 2 performs best, with a mean absolute error of 0.021 and 0.012 eV for the adiabatic and vertical ionization energies, respectively.     

Structural and vibrational analysis of the OH torsional motion in difluorohydroxyborane

In this work, we present a complete structural and vibrational analysis of the OH torsional motion in difluorohydroxyborane (BF₂OH) at the HF/aug‐cc‐pVTZ, MP2(full)/aug‐cc‐pVTZ, and CCSD/aug‐cc‐pVTZ theory levels. After full relaxation of the geometry, the equilibrium structure is found in a planar conformation of C_s symmetry. The difference in the two BF distances suggests the existence of a nonbonded interaction between the fluorine and oxygen atoms. The structural and energetic variation of BF₂OH as a function of the OH torsional angle is considered. The torsional barrier, at the CCSD/aug‐cc‐pVTZ level, and including the effect of the zero‐point energy of the remaining vibrations, is found 2,728 cm^?1. In addition, an anharmonic Hamiltonian for the OH torsional mode is presented and variationally solved. To simplify the treatment and to classify the energy levels, BF₂OH is classified under a G₄ nonrigid group accounting for the inversion symmetry of the molecule and the interchange of the fluorine atoms. The computed torsional energy levels exhibit a very small inversion splitting. The torsional spectrum is simulated considering the dipole moment components along the principal axes of inertia as an explicit function of the torsional coordinate. We observe three dominant bands in the spectrum formed by doublets corresponding to ν₉ = 0 → 1, 2 transitions. The fundamental is an a‐type, Franck–Condon, transition. This is the strongest and appears at 466.80 cm^?1 with relative intensity 0.4312. The ν₉ = 0 → 2 bands correspond to doublets of b‐ and c‐type, i.e., Herzberg–Teller transitions. These are two overlapping bands found at 890.92 and 890.94 cm^?1 with intensity 0.2207 for the b‐type band and 0.2193 for the c‐type band. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

A Comparative Study of Hydrogen Bonding Using Density Functional Theory

The hydrogen‐bond energies of water dimer and water‐formaldehyde complexes have been studied using density functional theory (DFT). Basis sets up to aug‐cc‐pVXZ (X=D, T, Q) were used. It was found that counterpoise corrected binding energies using the aug‐cc‐pVDZ basis set are very close to those predicted with the aug‐cc‐pVQZ set. Comparative studies using various DFT functionals on these two systems show that results from B3LYP, mPW1PW91 and PW91PW91 functionals are in better agreements with those predicted using high‐level ab initio methods. These functionals were applied to the study of hydrogen bonding between guanine (G) and cytosine (C), and between adenine (A) and thy mine (T) base pairs. With the aug‐cc‐pVDZ basis set, the predicted binding energies of base pairs are in good agreement with the most elaborate ab initio results.     

Absorption,resonance, and near‐resonance Raman studies of the tetracyanoquinodimethane neutral and its monoanion in terms of density functional theory and complete active space self‐consistent field methods

The electronic structure of the 1¹B_1u and 1²B_3u excited electronic states of the tetracyanoquinodimethane (TCNQ) neutral and its charged derivative are studied within the framework of complete active space self‐consistent field (CASSCF) and Becke's three‐parameter hybrid method with Lee–Yang–Parr correlation functional (B3LYP) methods applied to the level aug‐cc‐p‐VDZ basis set. Both CASSCF/aug‐cc‐p‐VDZ and B3LYP/aug‐cc‐p‐VDZ treatments provide the ground‐state and the excited state geometries; these are then used to assess the Franck–Condon (FC) parameters in the 1¹B_1u state of the neutral TCNQ and in the 1²B_3u state of the TCNQ monoanion. The quality of numerical results is then tested on the base of available experimental near‐resonance and resonance Raman data. The studies are performed in terms of the vibronic model, which takes both FC and mode‐mixing (Dushinsky) effects into account. This somewhat simplified vibronic model leads to very good agreement between the theory and the Raman experiments concerning both neutral TCNQ and its monoanion. In particular, the calculated excitation profiles of the ν₂ = 2215 cm^?1, ν₄ = 1389 cm^?1, ν₅ = 1195 cm^?1, and ν₉ = 336 cm^?1 fundamentals are shown to be in excellent agreement with those for the TCNQ monoanion. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2006     

Characterization of the potential energy surfaces of two small but challenging noncovalent dimers: (P2)2 and (PCCP)2

This work characterizes eight stationary points of the P₂ dimer and six stationary points of the PCCP dimer, including a newly identified minimum on both potential energy surfaces. Full geometry optimizations and corresponding harmonic vibrational frequencies were computed with the second‐order Møller–Plesset (MP2) electronic structure method and six different basis sets: aug‐cc‐pVXZ, aug‐cc‐pV(X+d)Z, and aug‐cc‐pCVXZ where X = T, Q. A new L‐shaped structure with C₂ symmetry is the only minimum for the P₂ dimer at the MP2 level of theory with these basis sets. The previously reported parallel‐slipped structure with C_2h symmetry and a newly identified cross configuration with D₂ symmetry are the only minima for the PCCP dimer. Single point energies were also computed using the canonical MP2 and CCSD(T) methods as well as the explicitly correlated MP2‐F12 and CCSD(T)‐F12 methods and the aug‐cc‐pVXZ (X = D, T, Q, 5) basis sets. The energetics obtained with the explicitly correlated methods were very similar to the canonical results for the larger basis sets. Extrapolations were performed to estimate the complete basis set (CBS) limit MP2 and CCSD(T) binding energies. MP2 and MP2‐F12 significantly overbind the P₂ and PCCP dimers relative to the CCSD(T) and CCSD(T)‐F12 binding energies by as much as 1.5 kcal mol^?1 for the former and 5.0 kcal mol^?1 for the latter at the CBS limit. The dominant attractive component of the interaction energy for each dimer configuration was dispersion according to several symmetry‐adapted perturbation theory analyses. © 2014 Wiley Periodicals, Inc.     

A DFT investigation of dissociation of nitrous acid,transformation of its monomeric isomers,their dimeric isomers,and dimerization

The structure optimization of the monomeric isomers of the nitrous acid HONO, their cyclic dimers (HONO)₂, and transition states of their transformations has been carried out using the density functional theory calculations at the B3LYP/6‐311++G(d,p) and B3LYP/aug‐cc‐pVQZ levels of theory. Three monomeric isomers (c, t, and z isomers) and six dimeric isomers (cc, ct, cz, tt, zz, and tz isomers) were found. Dissociation of the most stable of the HONO isomer (t isomer) was investigated and its rate constant was obtained at the UB3LYP/aug‐cc‐pVQZ level. Energetics, thermodynamic properties, rate, and equilibrium constants of the monomeric and dimeric isomerizations of the nitrous acid were obtained. Dimerization energies, enthalpies, and free energies of the dimeric isomers, (HONO)₂ were also obtained. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2010     

Franck–Condon breakdown in core-level photoelectron spectroscopy of chemisorbed CO

The photon energy dependence of the vibrational fine structure in the C1s and O1s X-ray photoelectron main lines of chemisorbed CO on Ni(100) and Ru(0001) has been measured from 6 to 150 eV above the core-level thresholds. Significant deviations from the behavior in gas-phase CO are found. A strong dominance of the adiabatic peak towards threshold is found for the C1s, but not the O1s, lines. In the C1s lines, we observe a broad maximum of vibrational excitation 5 eV above the shape resonance. At high photon energies, Franck–Condon behavior is observed in both the C1s and O1s lines. This behavior is discussed in terms of the adsorbate electronic structure and the dynamic metallic screening upon core ionization.     

Many‐body interactions of carbon monoxide cyclic oligomers: A computational study

Structural properties and energetics for the optimized carbon monoxide cyclic oligomers are analyzed at the correlated ab initio second‐order Møller–Plesset (MP2) and density functional methods (B3LYP and mPW1PW), using Dunning's cc‐pVXZ (X = D, T, Q) basis set, augmented with diffuse functions. Many‐body interactions of the stable carbon monoxide cyclic oligomers, (CO)₄ and (CO)₅ are computed at the MP2/aug‐cc‐pVTZ level. Contributions of two‐ to five‐body terms to each of these oligomers for their interaction energies, including corrections for basis set superposition error (BSSE) are investigated by using function counterpoise and its generalized version. It has been found that three‐body terms are attractive in nature and essential in order to describe the cooperative effects in the stable cyclic CO oligomers. © 2005 Wiley Periodicals, Inc. Int J Quantum Chem, 2005     

Theoretical computation of low‐lying electronic states of HCNS: A CASPT2 study

Complete active space self‐consistent field (CASSCF) and complete active space second‐order perturbation theory (CASPT2) calculations in conjunction with the aug‐cc‐pVTZ basis set have been used to investigate the low‐lying electronic states of thiofulminic acid (HCNS), HCNS⁺, and HCNS^?. The result of geometry optimization using CASPT2/aug‐cc‐pVTZ shows that theoretically determined geometric parameters and harmonic vibrational frequencies for the HCNS ground state X¹Σ⁺(X¹A′) are in agreement with previous studies. The ionization energies, the electron affinity energies, the adiabatic excitation energies, and vertical excitation energies have been calculated and the corresponding cation and anion states are identified. By calculating adiabatic electron affinity, the states of HCNS^? have been identified to contain both π orbital states (X²A′ and 1²A″) and dipole‐bond states (1⁴A′ and 1⁴A″). © 2012 Wiley Periodicals, Inc.     

Computational thermochemistry of glycolaldehyde

We use a variant of the focal point analysis to refine estimates of the relative energies of the four low‐energy torsional conformers of glycolaldehyde. The most stable form is the cis‐cis structure which enjoys a degree of H‐bonding from hydroxyl H to carbonyl O; here dihedral angles τ₁ (O?C? C? O) and τ₂ (C? C? O? H) both are zero. We optimized structures in both CCSD(T)/aug‐cc‐pVDZ and aug‐cc‐pVTZ; the structures agree within 0.01 Å for bond lengths and 1.0 degrees for valence angles, but the larger basis brings the rotational constants closer to experimental values. According to our extrapolation of CCSD(T) energies evaluated in basis sets ranging to aug‐cc‐pVQZ the trans‐trans form (180°, 180°) has a relative energy of 12.6 kJ/mol. The trans‐gauche conformer (160°, ±75°) is situated at 13.9 kJ/mol and the cis‐trans form (0°, 180°) at 18.9 kJ/mol. Values are corrected for zero point vibrational energy by MP2/aug‐cc‐pVTZ frequencies. Modeling the vibrational spectra is best accomplished by MP2/aug‐cc‐pVTZ with anharmonic corrections. We compute the Watsonian parameters that define the theoretical vibrational‐rotational spectra for the four stable conformers, to assist the search for these species in the interstellar medium. Six transition states are located by G4 and CBS‐QB3 methods as well as extrapolation using energies for structures optimized in CCSD(T)/aug‐cc‐pVDZ structures. We use two isodesmic reactions with two well‐established thermochemical computational schemes G4 and CBS‐QB3 to estimate energy enthalpy and Gibbs energy of formation as well as the entropy of the gas phase system. Our extrapolated electronic energies of species appearing in the isodesmic reactions produce independent values of thermodynamic quantities consistent with G4 and CBS‐QB3. © 2013 Wiley Periodicals, Inc.     

Potential‐Energy Surfaces,the Born–Oppenheimer Approximations,and the Franck–Condon Principle: Back to the Roots

The concept of a potential‐energy surface (PES) is central to our understanding of spectroscopy, photochemistry, and chemical kinetics. However, the terminology used in connection with the basic approximations is variously, and somewhat confusingly, represented with such phrases as “adiabatic”, “Born–Oppenheimer”, or “Born–Oppenheimer adiabatic” approximation. Concerning the closely relevant and important Franck–Condon principle (FCP), the IUPAC definition differentiates between a classical and quantum mechanical formulation. Consequently, in many publications we find terms such as “Franck–Condon (excited) state”, or a vertical transition to the “Franck–Condon point” with the “Franck–Condon geometry” that relaxes to the excited‐state equilibrium geometry. The Born–Oppenheimer approximation and the “classical” model of the Franck–Condon principle are typical examples of misused terms and lax interpretations of the original theories. In this essay, we revisit the original publications of pioneers of the PES concept and the FCP to help stimulate a lively discussion and clearer thinking around these important concepts.     

Compatibility of correlation‐consistent basis sets with a hybrid Hartree–Fock/density functional method

The present study examines the feasibility of combining the correlation‐consistent basis sets developed by Dunning and coworkers with the hybrid Hartree–Fock/density functional method B3LYP. Furthermore, extrapolation to the complete basis set (CBS) limit minimizes errors due to the presence of an incomplete basis set and can act as a rigorous test of the limitations of the B3LYP method. Equilibrium geometries, energies, and harmonic vibrational frequencies were determined for a series of well‐studied, yet computationally challenging, small inorganics and their respective ions. The results were then extrapolated to the CBS limit, where applicable, and compared to experiment. It was found that a union between the hybrid Hartree–Fock/density functional B3LYP method and Dunning's augmented correlation‐consistent basis sets gave results that were comparable to molecular orbital methods that explicitly account for electron correlation. Furthermore, the minimum basis set necessary to attain reasonable results for the systems studied was aug‐cc‐pVTZ. Upgrading to the aug‐cc‐pVQZ level and subsequent extrapolation to the CBS limit further improved the overall agreement with the experiment. ©1999 John Wiley & Sons, Inc. J Comput Chem 20: 207–216, 1999     

Multipolar electrostatics for proteins: Atom–atom electrostatic energies in crambin

Accurate electrostatics necessitates the use of multipole moments centered on nuclei or extra point charges centered away from the nuclei. Here, we follow the former alternative and investigate the convergence behavior of atom‐atom electrostatic interactions in the pilot protein crambin. Amino acids are cut out from a Protein Data Bank structure of crambin, as single amino acids, di, or tripeptides, and are then capped with a peptide bond at each side. The atoms in the amino acids are defined through Quantum Chemical Topology (QCT) as finite volume electron density fragments. Atom‐atom electrostatic energies are computed by means of a multipole expansion with regular spherical harmonics, up to a total interaction rank of L = ?_A+ ?_B + 1 = 10. The minimum internuclear distance in the convergent region of all the 15 possible types of atom‐atom interactions in crambin that were calculated based on single amino acids are close to the values calculated from di and tripeptides. Values obtained at B3LYP/aug‐cc‐pVTZ and MP2/aug‐cc‐pVTZ levels are only slightly larger than those calculated at HF/6‐31G(d,p) level. This convergence behavior is transferable to the well‐known amyloid beta polypeptide Aβ_1–42. Moreover, for a selected central atom, the influence of its neighbors on its multipole moments is investigated, and how far away this influence can be ignored is also determined. Finally, the convergence behavior of AMBER becomes closer to that of QCT with increasing internuclear distance. © 2013 Wiley Periodicals, Inc.