Reactivity of Hydrated Hydroxide Anion Clusters with H and Rb: An ab Initio Study.

We present a theoretical investigation of the hydrated hydroxide anion clusters, OH(H2O)n-, and of the collisional complexes, H-OH(H2O)n- and Rb-OH(H2O)n- (with n = 1-4). The MP2 and CCSD(T) methods are used to calculate interaction energies, optimized geometries, and vertical detachment energies. Parts of the potential energy surfaces are explored with a focus on the autodetachment region. We point out the importance of diffuse functions to correctly describe the latter. We use our results to discuss the different water loss and electronic detachment channels, which are the main reaction routes at room temperature and below. In particular, we have considered a direct and an indirect process for the electronic detachment, depending on whether water loss follows or precedes the detachment of the excess electron. We use our results to discuss the implications for astrochemistry and hybrid trap experiments in the context of cold chemistry.

Vibronic structure of the cyanobutadiyne cation. II. Theoretical exploration of the complex energy landscape of HC5N.

The results of an extensive ab initio study of the cyanobutadiyne cation, initially motivated by threshold-photoelectron spectroscopy experiments [see the study by Gans et al., J. Chem. Phys. 150, 244304 (2019)], are reported in the present paper. Calculations at the internally contracted multireference configuration interaction level of theory have been performed to derive the rovibronic properties of the seven lowest electronic states of HC5N+. Equilibrium geometries, rotational constants, vibrational frequencies, electric dipole moments, and spin-orbit constants have been calculated and compared with experimental data when available. Adiabatic and vertical ionization energies from the neutral ground state as well as transition energies within the cation electronic manifold are predicted, using the convergence to the complete basis set limit. The accurate description of the complex energy landscape up to 32 000 cm-1 above the ionization potential allows us to perform Franck-Condon simulations of the photoionization spectrum to the X+ 2Π, A+ 2Π, B+ 2Σ+, and C+ 2Π states and allows us to simulate the A+ 2Π â†’ X+ 2Π emission spectrum. The vibronic perturbations occurring on the excited potential energy surfaces are revealed and discussed, in particular, for the 3 2Π surface, which presents a double-well topography.

Vibronic structure of the cyanobutadiyne cation. I. VUV photoionization study of HC5N.

We report the vacuum-ultraviolet threshold-photoelectron spectrum of HC5N recorded over a wide spectral range, from 84 000 to 120 000 cm-1, with a 120 cm-1 spectral resolution, better than what was achieved in previous photoelectron studies, and with mass selectivity. The adiabatic ionization potential of cyanobutadiyne is measured at 85 366 (±40) cm-1. Assignment of the vibrational bands of the four lowest electronic states X+2Π, A+2Π, B+2Σ+, and C+2Π are performed, supported by high level ab initio calculations which are fully detailed in Paper II [B. Gans et al., J. Chem. Phys. 150, 244303 (2019)] and by Franck-Condon simulations. Only vibrational stretching modes are observed in the threshold-photoelectron spectra. The ground state of HC5N+ exhibits a vibrational progression in the ν2 stretching mode involving mainly the elongation of the C≡C triple bonds, whereas the A+ and C+ excited electronic states show a progression in the stretching mode mainly associated with the elongation of the C≡N bond, i.e., ν4 and ν3, respectively. The B+ state appears almost as a vibrationless structure in close vicinity to the A+ state.

Ab initio study of the neutral and anionic alkali and alkaline earth hydroxides: Electronic structure and prospects for sympathetic cooling of OH.

We have performed a systematic ab initio study on alkali and alkaline earth hydroxide neutral (MOH) and anionic (MOH-) species where M = Li, Na, K, Rb, Cs or Be, Mg, Ca, Sr, Ba. The CCSD(T) method with extended basis sets and Dirac-Fock relativistic effective core potentials for the heavier atoms has been used to study their equilibrium geometries, interaction energies, electron affinities, electric dipole moment, and potential energy surfaces. All neutral and anionic species exhibit a linear shape with the exception of BeOH, BeOH-, and MgOH-, for which the equilibrium structure is found to be bent. Our analysis shows that the alkaline earth hydroxide anions are valence-bound whereas the alkali hydroxide anions are dipole bound. In the context of sympathetic cooling of OH- by collision with ultracold alkali and alkaline earth atoms, we investigate the 2D MOH- potential energy surfaces and the associative detachment reaction M + OH→- MOH + e-, which is the only energetically allowed reactive channel in the cold regime. We discuss the implication for the sympathetic cooling of OH- and conclude that Li and K are the best candidates for an ultracold buffer gas.

Experimental and theoretical investigations of H2O-Ar.

We have used continuous-wave cavity ring-down spectroscopy to record the spectrum of H2O-Ar in the 2OH excitation range of H2O. 24 sub-bands have been observed. Their rotational structure (Trot = 12 K) is analyzed and the lines are fitted separately for ortho and para species together with microwave and far infrared data from the literature, with a unitless standard deviation σ=0.98 and 1.31, respectively. Their vibrational analysis is supported by a theoretical input based on an intramolecular potential energy surface obtained through ab initio calculations and computation of the rotational energy of sub-states of the complex with the water monomer in excited vibrational states up to the first hexad. For the ground and (010) vibrational states, the theoretical results agree well with experimental energies and rotational constants in the literature. For the excited vibrational states of the first hexad, they guided the assignment of the observed sub-bands. The upper state vibrational predissociation lifetime is estimated to be 3 ns from observed spectral linewidths.

Ab initio study of reactive collisions between Rb((2)S) or Rb((2)P) and OH(-)((1)Σ(+)).

A theoretical rate constant for the associative detachment reaction Rb((2)S) + OH(-)((1)Σ(+)) → RbOH((1)Σ(+)) + e(-) of 4 × 10(-10) cm(3) s(-1) at 300 K has been calculated. This result agrees with the experimental rate constant of 2-1 (+2)×10(-10)cm(3)s(-1) obtained by Deiglmayr et al. [Phys. Rev. A 86, 043438 (2012)] for a temperature between 200 K and 600 K. A Langevin-based dynamics which depends on the crossing point between the anion (RbOH(-)) and neutral (RbOH) potential energy surfaces has been used. The calculations were performed using the ECP28MDF effective core potential to describe the rubidium atom at the CCSD(T) level of theory and extended basis sets. The effect of ECPs and basis set on the height of the crossing point, and hence the rate constant, has been investigated. The temperature dependence of the latter is also discussed. Preliminary work on the potential energy surface for the excited reaction channel Rb((2)P) + OH(-)((1)Σ(+)) calculated at the CASSCF-icMRCI level of theory is presented. We qualitatively discuss the charge transfer and associative detachment reactions arising from this excited entrance channel.

Experimental and ab initio characterization of HC3N+ vibronic structure. II. High-resolution VUV PFI-ZEKE spectroscopy.

Vacuum-ultraviolet pulsed-field-ionization zero-kinetic-energy photoelectron spectra of X+Π2←XΣ+1 and B+Π2←XΣ+1 transitions of the HC314N and HC315N isotopologues of cyanoacetylene have been recorded. The resolution of the photoelectron spectra allowed us to resolve the vibrational structures and the spin-orbit splittings in the cation. Accurate values of the adiabatic ionization potentials of the two isotopologues (EI/hc(HC314N)=93 909(2) cm-1 and EI/hc(HC315N)=93 912(2) cm-1), the vibrational frequencies of the ν2, ν6, and ν7 vibrational modes, and the spin-orbit coupling constant (ASO = -44(2) cm-1) of the X+Π2 cationic ground state have been derived from the measurements. Using ab initio calculations, the unexpected structure of the B+Π2←XΣ+1 transition is tentatively attributed to a conical intersection between the A+ and B+ electronic states of the cation.

Experimental and ab initio characterization of HC3N+ vibronic structure. I. Synchrotron-based threshold photo-electron spectroscopy.

Threshold-photoionization spectroscopy of cyanoacetylene (HC3N) and its 15N isotopologue has been investigated in the vacuum-ultraviolet range with a synchrotron-based experiment allowing to record threshold-photoelectron spectrum and photoion yield over a large energy range (from 88 500 to 177 500 cm-1, i.e., from 11 to 22 eV). Adiabatic ionization energies towards the three lowest electronic states X+2Π, A+ Σ+2, and B+ Π2 are derived from the threshold-photoelectron spectrum. A detailed description of the vibrational structure of these states is proposed leading to the determination of the vibrational frequencies for most modes. The vibrational assignments and the discussion about the electronic structure are supported by multireference ab initio calculations (CASPT2, MRCI). Unprecedented structures are resolved and tentatively assigned in the region of the B+← X transition. Exploratory calculations highlight the complexity of the electronic landscape of the cation up to approximately 10 eV above its ground state.

Theoretical description of electronically excited vinylidene up to 10 eV: first high level ab initio study of singlet valence and Rydberg states.

The first quantitative description of the Rydberg and valence singlet electronic states of vinylidene lying in the 0-10 eV region is performed by using large scale ab initio calculations. A deep analysis of Rydberg-valence interactions has been achieved thanks to the comprehensive information contained in the accurate Multi-Reference Configuration Interaction wavefunctions and an original population analysis highlighting the respective role played by orbital and state mixing in such interactions. The present theoretical approach is thus adequate for dealing with larger than diatomic Rydberg systems. The nine lowest singlet valence states have been optimized. Among them, some are involved in strong Rydberg-valence interactions in the region of the Rydberg state equilibrium geometry. The Rydberg states of vinylidene present a great similarity with the acetylene isomer, concerning their quantum defects and Rydberg molecular orbital character. As in acetylene, strong s-d mixing is revealed in the n = 3 s-d supercomplex. Nevertheless, unlike in acetylene, the close-energy of the two vinylidene ionic cores (2)A1 and (2)B1 results into two overlapped Rydberg series. These Rydberg series exhibit local perturbations when an accidental degeneracy occurs between them and results in avoided crossings. In addition, some Δl = 1 (s-p and p-d) mixings arise for some Rydberg states and are rationalized in term of electrostatic interaction from the electric dipole moment of the ionic core. The strongest dipole moment of the (2)B1 cationic state also stabilizes the lowest members of the n = 3 Rydberg series converging to this excited state, as compared to the adjacent series converging toward the (2)A1 ionic ground state. The overall energies of vinylidene Rydberg states lie above their acetylene counterpart. Finally, predictions for optical transitions in singlet vinylidene are suggested for further experimental spectroscopic characterization of vinylidene.

Ab initio intermolecular potential of Ar-C2H2 refined using high-resolution spectroscopic data.

The high-resolution infrared spectra of the ν1 + ν3 (2CH) band of the Ar-C2H2 complex has been recorded from 6544 to 6566 cm(-1). The previously reported K(a) = 1 ← 0, 2 ← 1, and 0 ← 1 subbands were observed and the K(a) = 1 ← 2, 2 ← 3, and 3 ← 2 subbands were assigned for the first time. The intermolecular potential energy surface of this complex has been calculated ab initio and optimized by fitting the new high-resolution data. Refined intermolecular potential energy surfaces have been obtained for the ground vibrational state and for the excited v1 = v3 = 1 stretching state. For the former state, the results of the analysis are satisfactory and the microwave transitions of the complex are reproduced with a root-mean-square deviation of 5 MHz. For the latter state, systematic discrepancies arise in the analysis.

Determination of the ground electronic state in transition metal halides: ZrF.

The spectroscopy of the ZrF radical, produced by a laser ablation-molecular beam experimental setup, has been investigated for the first time using a two-color two-photon (1 + 1') REMPI scheme and time-of-flight (TOF) mass spectrometry detection. The region of intense bands 400-470 nm has been studied, based upon the first spectroscopic observations of the isovalent ZrCl radical by Carroll and Daly. The overall spectrum observed is complex. However, simultaneous and individual ion detection of the five naturally occurring isotopologues of ZrF has provided a crucial means of identifying band origins and characterization via the isotopic shift, δ(iso), of the numerous vibronic transitions recorded. Hence, five (0-0) transitions, of which only two were free of overlap with other transitions, have been identified. The most intense (0-0) transition at 23113 cm(-1) presented an unambiguously characteristic RQP rotational structure. From rotational contour simulations of the observed spectra, the nature of the ground electronic state is found to be unambiguously of (2)Δ symmetry, leading to the assignment of this band as (0-0) (2)Δ â† X(2)Δ at 23113 cm(-1). A set of transitions (1-0) (2)Δ â† X(2)Δ at 22105 cm(-1) and (2-0) (2)Φ â† X(2)Δ at 22944 cm(-1) involving the X(2)Δ state has also been identified and analyzed. Furthermore, a second series of transitions with lesser intensity has also been related to the long-lived metastable (4)Σ(-) state: (3-0) (4)Π(-1/2) ← (4)Σ(-) at 21801 cm(-1), (2-0) (4)Π(-1/2) ← (4)Σ(-) at 21285 cm(-1) and (2-0) (4)Σ(-) ← (4)Σ(-) at 23568 cm(-1). These spectroscopic assignments are supported by MRCI ab initio calculations, performed using the MOLPRO quantum chemistry package, and show that the low-lying excited states of the ZrF radical are the (4)Σ(-) and (4)Φ states lying at 2383 and 4179 cm(-1) respectively above the ground X(2)Δ state. The difference in the nature of ground state and ordering of the first electronic states with TiF (X(4)Φ)(2-4) and ZrCl,(5) respectively, is examined in terms of the ligand field theory (LFT)(7) applied to diatomic molecules. These results give a precise description of the electronic structure of the low lying electronic states of the ZrF transition metal radical.

The search for a deterministic origin for the presence of nonracemic amino-acids in meteorites: a computational approach.

Amino-acid enantiomeric excesses (ee's) have been detected in different types of carbonaceous chondrites, all in favor of the L enantiomer. In this article, we discuss possible deterministic causes to the presence of these amino-acid ee's in meteorites and evaluate in particular enantioselective photolysis by circularly polarized light (CPL). The electronic circular dichroism spectra of a set of amino- and hydroxy-acids, all detected in chondritic matter but some with ee's and others without ee's, were calculated and compared. The spectra were calculated for the most stable conformation(s) of the considered molecules using quantum mechanical methods (density functional theory). Our results suggest that CPL photolysis in the gas phase was perhaps not at the origin of the presence of ee's in meteorites and that the search for another, but still unknown, deterministic cause must be seriously undertaken.

Ab initio study of the electron transfer in an ionized stacked complex of guanines.

The charge transfer process in an ionized stacking of two consecutive guanines (G(5')G(3'))(+) has been studied by means of state-averaged CASSCF/MRCI and RASSCF/RASPT2 calculations. The ground and two first excited states of the radical cation have been characterized, and the topology of the corresponding potential energy surfaces (PESs) has been studied as a function of all intermolecular geometrical parameters. The results demonstrate that the charge transfer process in (G(5')G(3'))(+) is governed by the avoiding crossing between the ground and first excited states of the complex. Relative translation motions of both guanines in their molecular planes are shown to lead to the charge migration between G(5') and G(3'). Five stationary points (three minima and two saddle points) have been characterized along the reaction path describing the passage of the positive charge from G(5') to G(3'). The global minimum on the PES is found to correspond to the charge configuration G(5')(+)G(3'). The existence of an intermediate minimum along the reaction path has been established, characterizing a structure where the positive charge is equally distributed between the two guanines. The calculated energy profile allowed us to determine the height of the potential energy barrier (7.33 kcal/mol) and to evaluate the electronic coupling at a geometry close to the avoiding crossing (3.6 kcal/mol). Test calculations showed that the topology of the ground state PES of the complex GG(+) is qualitatively conserved upon optimization of the intramolecular geometrical parameters of the stationary points.

Revisiting Mulliken's concepts about Rydberg states and Rydberg-valence interactions from large-scale Ab initio calculations on the acetylene molecule.

A quantitative characterization of the Rydberg and valence singlet electronic states of acetylene lying in the 5-10.7 eV region is performed by using large-scale ab initio calculations. A special attention is paid on the comparison between the present calculations and Mulliken's concepts for Rydberg states, based on single-electron and single-configuration description. Most of the properties of the Rydberg states have been qualitatively understood via this comparison, mainly shown by the shape and size of the outer Rydberg molecular orbital. More quantitatively, Rydberg-valence mixing has been evaluated in several excited energy regions, as for instance, the interaction between the ' (1pig)2 1Ag doubly excited valence state and the manifold of electronic components of the np series, or the interaction between the 1pig 1Bu valence state and the 3dpig 1Sigma(u)+ Rydberg state. The rapid predissociation of the lowest 3s(sigma) 1Piu Rydberg state has been interpreted as a case of Rydbergization, earlier predicted by Mulliken.

Equilibrium structure and torsional barrier of BH3NH3.

Born-Oppenheimer equilibrium structures, r(e)(BO), of the electronic ground state of the borazane (BH3NH3) molecule of C3v point-group symmetry are computed ab initio using the CCSD(T) method with basis sets up to quintuple-zeta quality. Inclusion of the counterpoise correction and extrapolation of the structural parameters to the complete basis set limit yield a best estimate of r(e)(BO) of BH3NH3. The anharmonic force field of BH3NH3, computed at the CCSD(T) level of theory with a basis set of triple-zeta quality, allows the determination of semi-experimental equilibrium rotational constants, which in turn result in a semi-experimental equilibrium structure, r(e)(SE). The r(e)(BO) and r(e)(SE) structures are in excellent agreement, indicating the validity of the methods used for their determination. The empirical mass-dependent structure, r(m)(1), of BH3NH3 is also determined. Although it is inferior in quality to the previous two structures, it is much more accurate than the standard empirical r0 and r(s) structures reported earlier for BH3NH3. The semi-experimental r(e)(SE) as well as the empirical r(m)(1) structures determined are based on experimental ground-state rotational constants available from the literature for nine isotopologues of borazane. The effective barrier to the internal rotation of BH3NH3, a molecule isoelectronic with CH3CH3, has been computed ab initio, employing the focal-point analysis (FPA) approach, to be 699 +/- 11 cm(-1). This compares favorably with an empirical redetermination of the effective barrier based on the above r(e)(SE) structure, V3 = 718(17) cm(-1).

Sequence dependence of electron-induced DNA strand breakage revealed by DNA nanoarrays.

The electronic structure of DNA is determined by its nucleotide sequence, which is for instance exploited in molecular electronics. Here we demonstrate that also the DNA strand breakage induced by low-energy electrons (18 eV) depends on the nucleotide sequence. To determine the absolute cross sections for electron induced single strand breaks in specific 13 mer oligonucleotides we used atomic force microscopy analysis of DNA origami based DNA nanoarrays. We investigated the DNA sequences 5'-TT(XYX)3TT with X = A, G, C and Y = T, BrU 5-bromouracil and found absolute strand break cross sections between 2.66 · 10(-14) cm(2) and 7.06 · 10(-14) cm(2). The highest cross section was found for 5'-TT(ATA)3TT and 5'-TT(ABrUA)3TT, respectively. BrU is a radiosensitizer, which was discussed to be used in cancer radiation therapy. The replacement of T by BrU into the investigated DNA sequences leads to a slight increase of the absolute strand break cross sections resulting in sequence-dependent enhancement factors between 1.14 and 1.66. Nevertheless, the variation of strand break cross sections due to the specific nucleotide sequence is considerably higher. Thus, the present results suggest the development of targeted radiosensitizers for cancer radiation therapy.

Electron-attachment-induced DNA damage: instantaneous strand breaks.

Low energy electron-attachment-induced damage in DNA, where dissociation channels may involve multiple bonds including complex bond rearrangements and significant nuclear motions, is analyzed here. Quantum mechanics/molecular mechanics (QM/MM) calculations reveal how rearrangements of electron density after vertical electron attachment modulate the position and dynamics of the atomic nuclei in DNA. The nuclear motions involve the elongation of the P-O (P-O(3') and P-O(5')) and C-C (C(3')-C(4') and C(4')-C(5')) bonds for which the acquired kinetic energy becomes high enough so that the neighboring C(3')-O(3') or C(5')-O(5') phosphodiester bond may break almost immediately. Such dynamic behavior should happen on a very short time scale, within 15-30 fs, which is of the same order of magnitude as the time scale predicted for the excess electron to localize around the nucleobases. This result indicates that the C-O phosphodiester bonds can break before electron transfer from the backbone to the base.

Conformations consistent with charge migration observed in DNA and RNA X-ray structures.

Electron holes are known to migrate along the DNA or RNA duplexes and to localize preferentially on successive guanines. The stationary point conformations of Gua pairs that can trap or let pass these holes have been characterized by quantum chemistry calculations. Here we show their recurrent occurrence in DNA and RNA X-ray structures, often in quadruplex conformations or in interaction with proteins, ligands or metal ions. These findings give support to the biological, possibly regulatory, roles of charge migration in cell functioning.

A comparison of two methods for selecting vibrational configuration interaction spaces on a heptatomic system: ethylene oxide.

Two recently developed methods for solving the molecular vibrational Schrodinger equation, namely, the parallel vibrational multiple window configuration interaction and the vibrational mean field configuration interaction, are presented and compared on the same potential energy surface of ethylene oxide, c-C(2)H(4)O. It is demonstrated on this heptatomic system with strong resonances that both approaches converge towards the same fundamental frequencies. This confirms their ability to tackle the vibrational problem of large molecules for which full configuration interaction calculations are not tractable.

Ab initio study of the ionization of the DNA bases: ionization potentials and excited states of the cations.

The ionization of the four DNA bases is investigated by means of ab initio calculations. Accurate values of the gas-phase vertical and adiabatic ionization potentials (IP) are obtained at the MP2/6-31G(2d(0.8,alpha(d)),p) level of theory. The need of introducing extra polarization to the standard 6-31G(d,p) basis set is demonstrated by test calculations and an optimal value of alpha(d) = 0.1 is obtained. Ionization to electronically excited radical cations is also considered. The low-lying excited states of the cations are characterized for the first time. The topology of the corresponding potential energy surfaces is qualitatively described in terms of the stationary points (minima and saddle points) located on these surfaces. A conical intersection is characterized for the first time on the ground-state potential energy surface of all cations. It arises from the crossing of the adiabatic surfaces of the ground and first excited state at planar geometries. A nonplanar minimum is observed for the cytosine cation only. The geometry and electronic changes occurring along these surfaces are analyzed, leading to a comparison between the different nucleobase cations. The study of larger ionized systems related to DNA is rendered possible thanks to the optimized medium size basis set proposed in this work, as exemplified by the calculation of the IP of a stacked dimer of guanines.