Explicitly correlated multireference configuration interaction with multiple reference functions: Avoided crossings and conical intersections

Accelerated basis-set convergence of coupled-cluster excitation energies using the density-based basis-set correction method

We present the first application to real molecular systems of the recently proposed linear-response theory for the density-based basis-set correction method [J. Chem. Phys., 158, 234107 (2023)]. We apply this approach to accelerate the basis-set convergence of excitation energies in the equation-of-motion coupled-cluster singles and doubles (EOM-CCSD) method. We use an approximate linear-response framework that neglects the second-order derivative of the basis-set correction density functional and consists in simply adding to the usual Hamiltonian the one-electron potential generated by the first-order derivative of the functional. This additional basis-set correction potential is evaluated at the Hartree-Fock density, leading to a very computationally cheap basis-set correction. We tested this approach over a set of about 30 excitation energies computed for five small molecular systems and found that the excitation energies from the ground state to Rydberg states are the main source of basis-set error. These excitation energies systematically increase when the size of the basis set is increased, suggesting a biased description in favour of the excited state. Despite the simplicity of the present approach, the results obtained with the basis-set-corrected EOM-CCSD method are encouraging as they yield a mean absolute deviation of 0.02 eV for the aug-cc-pVTZ basis set, while it is 0.04 eV using the standard EOM-CCSD method. This might open the path to an alternative to explicitly correlated approaches to accelerate the basis-set convergence of excitation energies.

Photoelectron spectrum and breakdown diagram of ethanolamine: conformers, excited states, and thermochemistry

Advanced theoretical methodologies and photoelectron photoion coincidence spectroscopy were used to investigate the photoionization of ethanolamine in the 8-18 eV energy range. We identified the low-lying cation conformers and the excited cation electronic states after vertical excitation from the most stable neutral conformer computationally. The TPES is composed of broad, structureless bands because of unfavorable Franck-Condon factors for origin transitions upon ionization, populating the D0-D7 cationic states from the most stable neutral conformer, g'Gg'. The adiabatic ionization energy of ethanolamine is calculated at 8.940 ± 0.010 eV, and the 0 K appearance energies of aminomethylium, NH2CH2+ (+CH2OH), and methyleneammonium, NH3CH2+ (+H2CO), are determined experimentally to be 9.708 ± 0.010 eV and 9.73 ± 0.03 eV, respectively. The former is used to re-evaluate the ethanolamine enthalpy of formation in the gas and liquid phases as ΔfH⊖298K[NH2(CH2)2OH, g] = -208.2 ± 1.2 kJ mol-1 and ΔfH⊖298K[NH2(CH2)2OH, l] = -267.8 ± 1.2 kJ mol-1. This represents a substantial correction of the previous experimental determination and is supported by ab initio calculations.

Fully quantal description of combined internal conversion and intersystem crossing processes in the smallest Criegee intermediate CH2OO

A quantal description of nuclear motion using coupled fifteen-state potential energy and spin-orbit coupling surfaces for studying the photodissociation of CH2OO to H2CO(X1A1) + O1D and H2CO(X1A1) + O3P channels is presented. For the evaluation of surfaces, multireference electronic wave functions are employed. For the fully quantal description of the nuclear motion, we diabatize the PESs of the two and four lowest excited singlet and triplet states, respectively, within the three sets of vibronically coupled states, i.e. (B1A', C1A'), (a3A', b3A') and (a3A'', b3A''), employing the diabatization by ansatz method. This yields three different adiabatic-to-diabatic mixing angles, which are used for the diabatization of the spin-orbit coupling surfaces and allow to investigate simultaneously the internal conversion and intersystem crossing processes in CH2OO using a nuclear quantum dynamical approach for the first time. Our calculation predicts the presence of a weak spectral band with irregular and discrete structures, which originates from the role of spin-orbit couplings. This band of the spectrum is mainly located between the minimum energy of the a3A' state and the onset of the B1A' ← X1A' spectral band. Furthermore, sizable SOC between the B1A' and a3A'' states mixes them via intersystem crossing. Due to the vibronic interaction between the a3A'' and a3A' states, the molecule finally relaxes to the a3A' state and is dissociated to H2CO(X1A1) and O3P products.

New Potential Energy Surface for the H + Cl2 Reaction and Quantum Dynamics Studies

The reaction of H + Cl2 → HCl + Cl plays a crucial role in various fields. However, no previous study has investigated this reaction using accurate quantum mechanical methods. In this paper, we construct a global potential energy surface (PES) using the neural network method with more than 20,000 ab initio energies obtained by the MRCI-F12+Q method with the aug-cc-pV5Z basis and extrapolated to the complete basis set limit. The spin-orbit coupling of the Cl atom has been considered in the PES. With this new PES, product state-resolved quantum dynamics calculations for the H + Cl2 (v0 = 0, j0 = 0-2) → HCl + Cl reaction was carried out. Numerical results show that the initial rotational excitation of the Cl2 has negligible effects on the reactivity. Product state-resolved integral cross sections (ICS) and rate constants reveal that the HCl is most favorably produced in its v' = 2 vibrational state. The calculated product vibrational state-resolved and total reaction rate constants suggest that the new global PES is accurate enough, as compared with the available experimental measurements.

High-accuracy DMBE potential energy surface for CNO(A''4) and the rate coefficients for the C + NO reaction in the A'2, A''2, and A''4 states

An accurate potential energy surface (PES) for the lowest lying A''4 state of the CNO system is presented based on explicitly correlated multi-reference configuration interaction calculations with quadruple zeta basis set (MRCI-F12/cc-pVQZ-F12). The ab initio energies are fitted using the double many-body expansion method, thus incorporating long-range energy terms that can accurately describe the electrostatic and dispersion interactions with physically motivated decaying functions. Together with the previously fitted lowest A'2 and A''2 states using the same theoretical framework, this constitutes a new set of PESs that are suitable to predict rate coefficients for all atom-diatom reactions of the CNO system. We use this set of PESs to calculate thermal rate coefficients for the C(P3) + NO(Π2) reaction and compare the temperature dependence and product branching ratios with experimental results. The comparison between theory and experiment is shown to be improved over previous theoretical studies. We highlight the importance of the long-range interactions for low-temperature rate coefficients.

Potential energy surfaces for singlet and triplet states of the LiH2+ system and quasi-classical trajectory cross sections for H + LiH+ and H+ + LiH

A new set of six accurate ab initio potential energy surfaces (PESs) is presented for the first three singlet and triplet states of LiH2+ (1,21A', 11A'', 1,23A', and 13A'' states, where four of them are investigated for the first time), which have allowed new detailed studies gaining a global view on this interesting system. These states are relevant for the study of the most important reactions of lithium chemistry in the early universe. More than 45 000 energy points were calculated using the multi-reference configuration interaction level of theory using explicitly correlated methods (ic-MRCI-F12), and the results obtained for each individual electronic state were fitted to an analytical function. Using quasiclassical trajectories and considering the initial diatomic fragment in the ground rovibrational state, we have determined the integral cross sections for the H + LiH+(X2Σ+, C2Π) and H+ + LiH(X1Σ+, B1Π) reactions. In these calculations all available reaction channels were considered: the chemically most important H or H+ transfer/abstraction as well as atom exchange and collision induced dissociation for up to 1.0 eV of collision energy.

New accurate diabatic potential energy surfaces for the two lowest 1A'' states of H2S and photodissociation dynamics in its first absorption band

In this work, state-to-state photodissociation dynamics of H2S in its first absorption band has been studied quantum mechanically with a new set of coupled potential energy surfaces (PESs) for the first two 1A'' excited states, which were developed at the explicitly correlated internally contracted multi-reference configuration interaction level with the cc-pVQZ-F12 basis set and a large active space. The calculated absorption spectrum, product state distributions, and angular distributions are in excellent agreement with available experimental data, validating the accuracy of the PESs and the non-adiabatic couplings. Detailed analysis of the dynamics reveals that there are strong non-adiabatic couplings between the bound 11B1 and dissociative 11A2 states around the Franck-Condon region, leading to very fast predissociation to ro-vibrationally cold SH(X̃) fragments, during which marginal angular anisotropy of the PESs is involved. This study provides quantitatively accurate characterization of the electronic structure and detailed fragmentation dynamics of this prototypical photodissociation system, which is desirable for improving astrochemical modelling.

Correlated quantum treatment of the photodissociation dynamics of formaldehyde oxide CH2OO

An extended theoretical analysis of the photodissociation of the smallest Criegee intermediate CH₂OO following excitation to the B state is presented. It relies on explicitly correlated multireference electronic wavefunctions combined with a quantum dynamical treatment for two interacting (B-C) electronic states and three coupled nuclear degrees of freedom. The 3D model relies on PESs along the O-O and C-O stretching as well as C-O-O bending modes for the two lowest excited states with A' symmetry, and is sufficiently accurate to reproduce the experimental B¹A'-X¹A' absorption spectrum, especially at the low-energy range to unprecedented accuracy. The existence of a deep well (∼0.4 eV) on the (diabatic) B state causes a considerable amount of the wavepacket to be reflected by the B state well, which can explain the vibronic structures appearing in the long wavelength range of 360-470 nm of the spectrum. The main progression appearing in the energy range from 360 to 470 nm is assigned to the O-O stretching mode while finer details are affected by couplings to the C-O stretching and C-O-O bending modes. The weakly avoided crossing between the B-state and C-state potential energy surfaces appearing near 3.1 eV excitation energy (for RS2-F12 method) causes considerable disturbance in the vibronic fine structure of the bands. The description of the latter is quite strongly affected by the type of electron correlation treatment adopted, either fully variational (MRCI type) or perturbation theoretic (RS2 type). The results give novel insight into the complex interactions governing that intriguing process.

Spectroscopic Properties and Spin-Orbit Coupling of Electronic Excited States of GeCl. J Phys Chem A 2022;126:653-662. [PMID: 35084857 DOI: 10.1021/acs.jpca.1c08042] [Citation(s) in RCA: 0] [Impact Index Per Article: 0]

The electronic structure and spectroscopic properties of GeCl⁺ are studied by high-level ab initio calculations by considering spin-orbit coupling (SOC). The potential energy curves (PECs) and spectroscopic constants of 12 Λ-S states and 23 Ω states are calculated using the multi-reference configuration interaction plus Davidson correction method (MRCI + Q), which are in good agreement with the experiment. Based on the calculated SO matrix and the PECs of the Ω states, the interaction between the d³Π state and other states caused by the SOC and the double-potential well structure caused by the avoided crossing rule and the properties of transitions d³Π₀₊-X¹Σ⁺₀₊, d³Π₁-X¹Σ⁺₀₊, and a³Σ⁺₁-X¹Σ⁺₀₊ are studied. Our results indicate that the previously observed spectra of GeCl⁺ in the 290-325 nm range should be assigned as the a³Σ⁺₁-X¹Σ⁺₀₊ transition. Moreover, the predissociation behavior of the d³Π state between the vibrational levels v' = 1 and v' = 10 is discussed, and the radiative lifetimes of transitions d³Π₀₊-X¹Σ⁺₀₊ and d³Π₁-X¹Σ⁺₀₊ are evaluated on the order of microseconds, while a³Σ⁺₁-X¹Σ⁺₀₊ is on the order of milliseconds. We estimate that the strongest bands of a³Σ⁺₁-X¹Σ⁺₀₊ are the 0-16, 0-17, and 0-18 bands. This study will promote our understanding of the detailed electronic structure and spectra of the GeCl⁺ radical cation.

Cryogenic IR and UV spectroscopy of isomer-selected cytosine radical cation

The UV photodissociation of cryogenic-cooled isomer-selected cytosine–silver complex leads to the production of cytosine radical cation without isomerization.

Probing the dynamics of the photo-induced decarboxylation of neutral and ionic pyruvic acid

The dynamics of the electronically excited pyruvic acid (PA) and of its unimolecular decomposition upon single photon ionization are investigated by means of a table top fs laser and VUV...

Correlation consistent basis sets for explicitly correlated wavefunctions: Pseudopotential-based basis sets for the group 11 (Cu, Ag, Au) and 12 (Zn, Cd, Hg) elements

New correlation consistent basis sets for the group 11 (Cu, Ag, Au) and 12 (Zn, Cd, Hg) elements have been developed specifically for use in explicitly correlated F12 calculations. This includes orbital basis sets for valence only (cc-pVnZ-PP-F12, n = D, T, Q) and outer core-valence (cc-pCVnZ-PP-F12) correlation, along with both of these augmented with additional high angular momentum diffuse functions. Matching auxiliary basis sets required for density fitting and resolution-of-the-identity approaches to conventional and F12 integrals have also been optimized. All of the basis sets are to be used in conjunction with small-core relativistic pseudopotentials [Figgen et al., Chem. Phys. 311, 227 (2005)]. The accuracy of the basis sets is determined through benchmark calculation at the explicitly correlated coupled-cluster level of theory for various properties of atoms and diatomic molecules. The convergence of the properties with respect to the basis set is dramatically improved compared to conventional coupled-cluster calculations, with cc-pVTZ-PP-F12 results close to conventional estimates of the complete basis set limit. The patterns of convergence are also greatly improved compared to those observed from the use of conventional correlation consistent basis sets in F12 calculations.

Threshold photoelectron spectroscopy of 9-methyladenine: theory and experiment

We present a combined experimental and theoretical study of single-photon ionization of 9-methyladenine (9MA) in the gas phase. In addition to tautomerism, several rotamers due to the rotation of the methyl group may exist. Computations show, however, that solely one rotamer contributes because of low population in the molecular beam and/or unfavorable Franck-Condon factors upon ionization. Experimentally, we used VUV radiation available at the DESIRS beamline of the synchrotron radiation facility SOLEIL to record the threshold photoelectron spectrum of this molecule between 8 and 11 eV. This spectrum consists of a well-resolved band assigned mainly to vibronic levels of the D₀ cationic state, plus a contribution from the D₁ state, and two large bands corresponding to the D₁, D₂ and D₃ electronically excited states. The adiabatic ionization energy of 9MA is measured at 8.097 ± 0.005 eV in close agreement with the computed value using the explicitly correlated coupled cluster approach including core valence, scalar relativistic and zero-point vibrational energy corrections. This work sheds light on the complex pattern of the lowest doublet electronic states of 9MA⁺. The comparison to canonical adenine reveals that methylation induces further electronic structure complication that may be important to understand the effects of ionizing radiation and the charge distribution in these biological entities at different time scales.

Study of the CN(X2Σ+) + N(4S) Reaction at High Temperatures: Potential Energy Surface and Thermal Rate Coefficients

Reactions involving C and N play an essential role in the chemistry around the surface of a hypersonic spacecraft during its atmospheric re-entry. The collision of CN with other molecules and atoms has particular interest in aerothermodynamic modeling. This work focuses on the study of the CN + N → N₂ + C reaction in the triplet manifold ³A″ of CN₂. A high-level full-dimensional potential energy surface for this system is developed from ab initio calculations at the MRCI-F12 + Q level of theory. This surface is employed in quasiclassical trajectory calculations, and thermal rate coefficients from 100 to 20,000 K are computed. The rates for the formation of N₂ are compared with the available experimental data, and good agreement is found. At low and intermediate temperatures, the N₂ formation is more efficient than the N-exchange process, while at high temperatures, the rates for both processes are comparable. Finally, analytically modified Arrhenius expressions for the reaction rates of N₂ formation and N-exchange are reported.

Nonadiabatic reaction dynamics to silicon monosulfide (SiS): A key molecular building block to sulfur-rich interstellar grains

Sulfur- and silicon-containing molecules are omnipresent in interstellar and circumstellar environments, but their elementary formation mechanisms have been obscure. These routes are of vital significance in starting a chain of chemical reactions ultimately forming (organo) sulfur molecules-among them precursors to sulfur-bearing amino acids and grains. Here, we expose via laboratory experiments, computations, and astrochemical modeling that the silicon-sulfur chemistry can be initiated through the gas-phase reaction of atomic silicon with hydrogen sulfide leading to silicon monosulfide (SiS) via nonadiabatic reaction dynamics. The facile pathway to the simplest silicon and sulfur diatomic provides compelling evidence for the origin of silicon monosulfide in star-forming regions and aids our understanding of the nonadiabatic reaction dynamics, which control the outcome of the gas-phase formation in deep space, thus expanding our view about the life cycle of sulfur in the galaxy.

Accurate DMBE potential-energy surface for CNO(2A″) and rate coefficients in C(3P)+NO collisions

A realistic double many-body expansion potential energy surface (PES) is developed for the ²A″ state of the carbon-nitrogen-oxygen (CNO) system based on MRCI-F12/cc-pVQZ-F12 ab initio energies. The new PES reproduces the fitted points with chemical accuracy (root mean square deviation up to 0.043 eV) and explicitly incorporates long range energy terms that can accurately describe the electrostatic and dispersion interactions. Thermal rate coefficients were computed for the C(³P) + NO(²Π) reaction for temperatures ranging from 15 K to 10 000 K, and the values are compared to previously reported results. The differences are rationalized, and the major importance of long range forces in predicting the rate coefficients for barrierless reactions is emphasized.

An accurate full-dimensional potential energy surface for the reaction OH + SO → H + SO2

An accurate full-dimensional PES for the OH + SO ↔ H + SO₂ reaction is developed by the permutation invariant polynomial-neural network approach.

The Molpro quantum chemistry package

Molpro is a general purpose quantum chemistry software package with a long development history. It was originally focused on accurate wavefunction calculations for small molecules but now has many additional distinctive capabilities that include, inter alia, local correlation approximations combined with explicit correlation, highly efficient implementations of single-reference correlation methods, robust and efficient multireference methods for large molecules, projection embedding, and anharmonic vibrational spectra. In addition to conventional input-file specification of calculations, Molpro calculations can now be specified and analyzed via a new graphical user interface and through a Python framework.

Spectroscopic characterization of the first excited state and photochemistry of the HO3 radical

We report the one-dimensional cuts of the six-dimensional potential energy surfaces (PESs) of the ground and lowest doublet and quartet electronic states of trans-HO₃ at the MRCI-F12/aug-cc-pVTZ level of theory. Theoretical calculations predict that the first excited state (A²A) presents a real minimum on its PES and possesses a nonplanar structure. The adiabatic excitation energy at the MRCI+Q and MRCI-F12 levels shows that the A²A state lies in the near-infrared region. Both the transition dipole moment and the oscillator strength were predicted to be weak, which suggests that photodissociation of HO₃ to produce OH and O₂ after UV-Vis absorption is not a plausible mechanism. The harmonic vibrational frequencies and rotational constants of the weakly bound complex OH-O₂ in the two electronic states were predicted to help in its detection. Our PES shows that the reactions of H + O₃ or HO₂ + O in their ground states do not lead to trans-HO₃ in its ground electronic state if one of the component fragments, i.e., HO₂(A²A') + O(³P) or H(²S) + O₃(³B₂), is excited.

Quasiclassical Study of the C(3P) + NO(X2Π) and O(3P) + CN(X2Σ+) Collisional Processes on an Accurate DMBE Potential Energy Surface

The predicted rate constants for C + NO and O + CN collisions in three potential energy surfaces (PESs) for the ²A' state of the CNO molecule are compared using quasiclassical trajectories. Different temperature dependencies are obtained for the C + NO reaction, which are explained in terms of the long-range properties of the PESs. Recommended values and mechanistic details are also reported. For O + CN collisions, a better agreement between the theoretical results is found, except for temperatures below 100 K.

Total Absorption Cross Section for UV Excitation of Sulfur Monoxide

In a new study on sulfur monoxide, new ab initio potential energy curves are developed for excited states A³Π, B³Σ^-, C³Π, and C'³Π and for three states that are unassigned in the literature, which we numerically name 3³Σ^-, 4³Π, and 5³Π. All these excited states have allowed transitions from ground state, X³Σ^-. The ab initio calculations were performed using the MRCI-F12+Q/aug-cc-pV(5+d)Z level of theory implemented in MOLPRO2015. On the basis of close-coupling R-matrix theory, fine structure absorption cross sections of isotopically substituted sulfur monoxide are calculated for wavelengths of 190-300 nm. The spectra are shown at the highest possible resolution (FWHH ≈ 0.15-0.18 cm^-1) for reference and future studies. The effects of self-shielding and possible mutual shielding are discussed.

Quantum Behavior of Spin-Orbit Inelastic Scattering of C-Atoms by D2 at Low Energy

Fine-structure populations and collision–induced energy transfer in atoms are of interest for many fields, from combustion to astrophysics. In particular, neutral carbon atoms are known to play a role in interstellar media, either as probes of physical conditions (ground state ³P_j spin-orbit populations), or as cooling agent (collisional excitation followed by radiative decay). This work aims at investigating the spin-orbit excitation of atomic carbon in its ground electronic state due to collisions with molecular deuterium, an isotopic variant of H₂, the most abundant molecule in the interstellar medium. Spin-orbit excitations of C(³P_j) by H₂ or D₂ are governed by non-adiabatic and spin-orbit couplings, which make the theoretical treatment challenging, since the Born-Oppenheimer approximation no longer holds. Inelastic collisional cross-sections were determined for the C(³P₀) + D₂ → C(³P_j) + D₂ (with j = 1 and 2) excitation process. Experimental data were acquired in a crossed beam experiment at low collision energies, down to the excitation thresholds (at 16.42 and 43.41 cm⁻¹, respectively). C-atoms were produced mainly in their ground spin-orbit state, ³P₀, by dissociation of CO in a dielectric discharge through an Even-Lavie pulsed valve. The C-atom beam was crossed with a D₂ beam from a second valve. The state-to-state cross-sections were derived from the C(³P_j) (j = 1 or 2) signal measured as a function of the beam crossing angle, i.e., as a function of the collision energy. The results show different quantum behaviors for excitation to C(³P₁) or C(³P₂) when C(³P₀) collides with ortho-D₂ or normal-D₂. These experimental results are analyzed and discussed in the light of highly accurate quantum calculations. A good agreement between experimental and theoretical results is found. The present data are compared with those obtained for the C-He and C-H₂ collisional systems to get new insights into the dynamics of collision induced spin-orbit excitation/relaxation of atomic carbon.

Two-State, Three-Mode Parametrization of the Force Field of a Retinal Chromophore Model

In recent years, the potential energy surfaces of the penta-2,4-dieniminium cation have been investigated using several electronic structure methods. The resulting pool of geometrical, electronic, and energy data provides a suitable basis for the construction of a topographically correct analytical model of the molecule force field and, therefore, for a better understanding of this class of molecules, which includes the chromophore of visual pigments. In the present contribution, we report the construction of such a model for regions of the force field that drive the photochemical and thermal isomerization of the central double bound of the cation. While previous models included only two modes, it is here shown that the proposed three-mode model and corresponding set of parameters are able to reproduce the complex topographical and electronic structure features seen in electronically correlated data obtained at the XMCQDPT2//CASSCF/6-31G* level of theory.

EVB and polarizable MM study of energy relaxation in fluorine–acetonitrile reactions

Many-body effects can impact on rates of energy transfer from a ‘hot’ DF solute to acetonitrile solvent.

The effect of nonadiabaticity on the C+ + HF reaction

The chemistry of fluorine in the interstellar medium is particularly simple, with only a few key species and important reactions. Of the latter, the rate of the reaction of C⁺ ions with HF is not well established but is one of the key reactions that sets the relative abundance of HF and the CF⁺ ion, the two fluorine-bearing species that have been observed in interstellar clouds. The C⁺ + HF → CF⁺ + H reaction proceeds through a deeply bound HCF⁺ well. In this work, statistical methods, namely, the statistical adiabatic channel method originally developed by Quack and Troe and the quantum statistical method of Manolopoulos and co-workers, are applied to compute the total cross section as a function of energy for this reaction. This reaction proceeds on the ground 1² A' potential energy surface (PES), and there are also two non-reactive PES's, 1² A″ and 2² A', correlating with the C⁺(² P _1/2,3/2) + HF reactants. Two sets of scattering calculations were carried out, namely, a single-surface calculation on the 1² A' PES and the one in which all three PES's and the spin-orbit splitting of C⁺ are included in the description of the entrance channel. In the latter, reactivity of the spin-orbit excited ² P _3/2 level can be computed, and not just assumed to be zero, as in the single-state adiabatic approximation.

Unveiling the complex vibronic structure of the canonical adenine cation

Adenine, a DNA base, exists as several tautomers and isomers that are closely lying in energy and that may form a mixture upon vaporization of solid adenine. Indeed, it is challenging to bring adenine into the gas phase, especially as a unique tautomer. The experimental conditions were tuned to prepare a jet-cooled canonical adenine (9H-adenine). This isolated DNA base was ionized by single VUV photons from a synchrotron beamline and the corresponding slow photoelectron spectrum was compared to ab initio computations of the neutral and ionic species. We report the vibronic structure of the X+ 2A'' (D0), A+ 2A' (D1) and B+ 2A'' (D2) electronic states of the 9H adenine cation, from the adiabatic ionization energy (AIE) up to AIE + 1.8 eV. Accurate AIEs are derived for the 9H-adenine (X[combining tilde] 1A') + hν → 9H-adenine+ (X+ 2A'', A+ 2A', B+ 2A'') + e- transitions. Close to the AIE, we fully assign the rich vibronic structure solely to the 9H-adenine (X 1A') + hν → 9H-adenine+ (X+ 2A'') transition. Importantly, we show that the lowest cationic electronic states of canonical adenine are coupled vibronically. The present findings are important for understanding the effects of ionizing radiation and the charge distribution on this elementary building block of life, at ultrafast, short, and long timescales.

Crystal Field Splittings in Lanthanide Complexes: Inclusion of Correlation Effects beyond Second Order Perturbation Theory

State-averaged complete active space self-consistent field (CASSCF) calculations and a subsequent spin-orbit calculation mixing the CASSCF wave functions (CASSCF/state-interaction with spin-orbit coupling) is the conventional approach used for ab initio calculations of crystal-field splittings and magnetic properties of lanthanide complexes. However, this approach neglects dynamical correlation. Complete active space second-order perturbation theory (CASPT2) can be used to account for dynamical correlation but suffers from the well-known problems of multireference perturbation theory, e.g., intruder state problems. Variational multireference configuration interaction (MRCI) calculations do not show these problems but are usually not feasible due to the large size of real lanthanide complexes. Here, we present a quasi-local projected internally contracted MRCI approach which makes MRCI calculations of lanthanide complexes feasible and allows assessing the influence of dynamical correlation beyond second-order perturbation theory. We apply the method to two well-studied molecules, namely, [Er{N(SiMe₃)₂}₃] and {C(NH₂)₃}₅[Er(CO₃)₄]·11H₂O.

UV photodissociation dynamics of CHI2Cl and its role as a photolytic precursor for a chlorinated Criegee intermediate

Photolysis of geminal diiodoalkanes in the presence of molecular oxygen has become an established route to the laboratory production of several Criegee intermediates, and such compounds also have marine sources. Here, we explore the role that the trihaloalkane, chlorodiiodomethane (CHI₂Cl), may play as a photolytic precursor for the chlorinated Criegee intermediate ClCHOO. CHI₂Cl has been synthesized and its UV absorption spectrum measured; relative to that of CH₂I₂ the spectrum is shifted to longer wavelength and the photolysis lifetime is calculated to be less than two minutes. The photodissociation dynamics have been investigated using DC slice imaging, probing ground state I and spin-orbit excited I* atoms with 2 + 1 REMPI and single-photon VUV ionization. Total translational energy distributions are bimodal for I atoms and unimodal for I*, with around 72% of the available energy partitioned in to the internal degrees of freedom of the CHICl radical product, independent of photolysis wavelength. A bond dissociation energy of D₀ = 1.73 ± 0.11 eV is inferred from the wavelength dependence of the translational energy release, which is slightly weaker than typical C-I bonds. Analysis of the photofragment angular distributions indicate dissociation is prompt and occurs primarily via transitions to states of A'' symmetry. Complementary high-level MRCI calculations, including spin-orbit coupling, have been performed to characterize the excited states and confirm that states of A'' symmetry with highly mixed singlet and triplet character are predominantly responsible for the absorption spectrum. Transient absorption spectroscopy has been used to measure the absorption spectrum of ClCHOO produced from the reaction of CHICl with O₂ over the range 345-440 nm. The absorption spectrum, tentatively assigned to the syn conformer, is at shorter wavelengths relative to that of CH₂OO and shows far weaker vibrational structure.

Isomer Identification in Flames with Double-Imaging Photoelectron/Photoion Coincidence Spectroscopy (i2PEPICO) using Measured and Calculated Reference Photoelectron Spectra

Double-imaging photoelectron/photoion coincidence (i2PEPICO) spectroscopy using a multiplexing, time-efficient, fixed-photon-energy approach offers important opportunities of gas-phase analysis. Building on successful applications in combustion systems that have demonstrated the discriminative power of this technique, we attempt here to push the limits of its application further to more chemically complex combustion examples. The present investigation is devoted to identifying and potentially quantifying compounds featuring five heavy atoms in laminar, premixed low-pressure flames of hydrocarbon and oxygenated fuels and their mixtures. In these combustion examples from flames of cyclopentene, iso-pentane, iso-pentane blended with dimethyl ether (DME), and diethyl ether (DEE), we focus on the unambiguous assignment and quantitative detection of species with the sum formulae C5H6, C5H7, C5H8, C5H10, and C4H8O in the respective isomer mixtures, attempting to provide answers to specific chemical questions for each of these examples. To analyze the obtained i2PEPICO results from these combustion situations, photoelectron spectra (PES) from pure reference compounds, including several examples previously unavailable in the literature, were recorded with the same experimental setup as used in the flame measurements. In addition, PES of two species where reference spectra have not been obtained, namely 2-methyl-1-butene (C5H10) and the 2-cyclopentenyl radical (C5H7), were calculated on the basis of high-level ab initio calculations and Franck-Condon (FC) simulations. These reference measurements and quantum chemical calculations support the early fuel decomposition scheme in the cyclopentene flame towards 2-cyclopentenyl as the dominant fuel radical as well as the prevalence of branched intermediates in the early fuel destruction reactions in the iso-pentane flame, with only minor influences from DME addition. Furthermore, the presence of ethyl vinyl ether (EVE) in DEE flames that was predicted by a recent DEE combustion mechanism could be confirmed unambiguously. While combustion measurements using i2PEPICO can be readily obtained in isomer-rich situations, we wish to highlight the crucial need for high-quality reference information to assign and evaluate the obtained spectra.

Theoretical investigation of rotationally inelastic collisions of CH(X2Π) with hydrogen atoms

We report calculations of state-to-state cross sections for collision-induced rotational transitions of CH(X²Π) with atomic hydrogen. These calculations employed the four adiabatic potential energy surfaces correlating CH(X²Π) + H(²S), computed in this work through the multi-reference configuration interaction method [MRCISD + Q(Davidson)]. Because of the presence of deep wells on three of the potential energy surfaces, the scattering calculations were carried out using the quantum statistical method of Manolopoulos and co-workers [Chem. Phys. Lett. 343, 356 (2001)]. The computed cross sections included contributions from only direct scattering since the CH₂ collision complex is expected to decay predominantly to C + H₂. Rotationally energy transfer rate constants were computed for this system since these are required for astrophysical modeling.

Advances in spectroscopy and dynamics of small and medium sized molecules and clusters

Investigations of the spectroscopy and dynamics of small- and medium-sized molecules and clusters represent a hot topic in atmospheric chemistry, biology, physics, atto- and femto-chemistry and astrophysics.

Stereoisomers of hydroxymethanes: Probing structural and spectroscopic features upon substitution

Ab initio studies on CH_x(OH)_4-x (x = 0-3) polyols are carried out to derive their structural and spectroscopic features. Several stereoisomers (both equilibrium structures and transition states) are found. Some are predicted here for the first time. We determined hence their geometrical parameters, vibrational frequencies, electronic excitation energies for the singlet manifold, and IR spectra. While the IR spectra for all polyols present similar shapes, their UV spectra exhibit however distinct band origin that are specific to each polyol and more interestingly to each diasteroisomer. Stereoelectronic effects are also noticed and discussed. It is suggested that UV spectroscopy is an efficient probe to experimentally identify polyols in mixtures involving polyols.

Induction of unidirectional π-electron rotations in low-symmetry aromatic ring molecules using two linearly polarized stationary lasers

A new laser-control scenario of unidirectional π-electron rotations in a low-symmetry aromatic ring molecule having no degenerate excited states is proposed. This scenario is based on dynamic Stark shifts of two relevant excited states using two linearly polarized stationary lasers. Each laser is set to selectively interact with one of the two electronic states, the lower and higher excited states are shifted up and down with the same rate, respectively, and the two excited states become degenerate at their midpoint. One of the four control parameters of the two lasers, i.e. two frequencies and two intensities, determines the values of all the other parameters. The direction of π-electron rotations, clockwise or counter-clockwise rotation, depends on the sign of the relative phase of the two lasers at the initial time. An analytical expression for the time-dependent expectation value of the rotational angular momentum operator is derived using the rotating wave approximation (RWA). The control scenario depends on the initial condition of the electronic states. The control scenario with the ground state as the initial condition was applied to toluene molecules. The derived time-dependent angular momentum consists of a train of unidirectional angular momentum pulses. The validity of the RWA was checked by numerically solving the time-dependent Schrödinger equation. The simulation results suggest an experimental realization of the induction of unidirectional π-electron rotations in low-symmetry aromatic ring molecules without using any intricate quantum-optimal control procedure. This may open up an effective generation method of ring currents and current-induced magnetic fields in biomolecules such as amino acids having aromatic ring molecules for searching their interactions.