Probing the Effects of Size and Charge on the Monohydration and Dihydration of SiF5- and SiF62- via Comparisons with BF4- and PF6.

This study systematically examines the interactions of the trigonal bipyramidal silicon pentafluoride and octahedral silicon hexafluoride anions with either one or two water molecules, (SiF5-(H2O)n and SiF62-(H2O)n, respectively, where n = 1, 2). Full geometry optimizations and subsequent harmonic vibrational frequency computations are performed using the CCSD(T) ab initio method with a triple-ζ correlation consistent basis set augmented with diffuse functions on all non-hydrogen atoms (cc-pVTZ for H and aug-cc-pVTZ for Si, O, and F; denoted as haTZ). Two monohydrate and six dihydrate minima have been identified for the SiF5-(H2O)n systems, whereas one monohydrate and five dihydrate minima have been identified for the SiF62-(H2O)n systems. Both monohydrated anions have a minimum in which the water molecule adopts a symmetric double ionic hydrogen bond (DIHB) motif with C2v symmetry. However, a second unique monohydrate minimum has been identified for SiF5- in which the water molecule adopts an asymmetric DIHB motif along the edge of the trigonal bipyramidal anion between one axial and one equatorial F atom. This Cs structure is more than 2 kcal mol-1 lower in energy than the C2v local minimum at the CCSD(T)/haTZ level of theory. While the interactions between the solvent and ionic solute are quite strong for the monohydrated anions (electronic dissociation energies of ≈12 and ≈24 kcal mol-1 for the SiF5-(H2O)1 and SiF62-(H2O)1 global minima, respectively), these values are nearly perfectly doubled for the dihydrates, with the lowest-energy SiF5-(H2O)2 and SiF62-(H2O)2 minima exhibiting dissociation energies of ≈24 and ≈47 kcal mol-1, respectively. Structures that form hydrogen bonds between the solvating water molecules also exhibit the largest shifts in the harmonic OH stretching frequencies for the waters of hydration. These shifts can exceed -100 cm-1 for the SiF5-(H2O)2 minimum and -300 cm-1 for the SiF62-(H2O)2 minimum relative to an isolated H2O molecule at the CCSD(T)/haTZ level of theory. This work also corrects the OH stretching frequency shifts for two dihydrate minima of PF6- that were previously erroneously reported ( J. Phys. Chem. A 2020, 124, 8744-8752, DOI: 10.1021/acs.jpca.0c06466).

Phospha-RosIndolizine Dye with Shortwave Infrared (SWIR) Absorption and Emission.

Shortwave infrared (SWIR, 1000-1700 nm) absorbing and emitting dyes are needed for infrared diodes and sensors used in a wide variety of industrial and medical applications. Herein, an electron-withdrawing phosphine oxide (P═O) substituted xanthene is coupled with strong indolizine donors to produce a SWIR absorbing (λabs = 1294 nm in DCM) and emitting (λemis = 1450 nm in DCM) dye called PRos1450. The unique properties of this dye are characterized via photophysical, electrochemical, and computational analyses.

Systematic analysis of electronic barrier heights and widths for concerted proton transfer in cyclic hydrogen bonded clusters: (HF)n, (HCl)n and (H2O)n where n = 3, 4, 5.

The MP2 and CCSD(T) methods are paired with correlation consistent basis sets as large as aug-cc-pVQZ to optimize the structures of the cyclic minima for (HF)n, (HCl)n and (H2O)n where n = 3-5, as well as the corresponding transition states (TSs) for concerted proton transfer (CPT). MP2 and CCSD(T) harmonic vibrational frequencies confirm the nature of each minimum and TS. Both conventional and explicitly correlated CCSD(T) computations are employed to assess the electronic dissociation energies and barrier heights for CPT near the complete basis (CBS) limit for all 9 clusters. Results for (HF)n are consistent with prior studies identifying Cnh and Dnh point group symmetry for the minima and TSs, respectively. Our computations also confirm that CPT proceeds through Cs TS structures for the C1 minima of (H2O)3 and (H2O)5, whereas the process goes through a TS with D2d symmetry for the S4 global minimum of (H2O)4. This work corroborates earlier findings that the minima for (HCl)3, (HCl)4 and (HCl)5 have C3h, S4 and C1 point group symmetry, respectively, and that the Cnh structures are not minima for n = 4 and 5. Moreover, our computations show the TSs for CPT in (HCl)3, (HCl)4 and (HCl)5 have D3h, D2d, and C2 point group symmetry, respectively. At the CCSD(T) CBS limit, (HF)4 and (HF)5 have the smallest electronic barrier heights for CPT (≈15 kcal mol-1 for both), followed by the HF trimer (≈21 kcal mol-1). The barriers are appreciably higher for the other clusters (around 27 kcal mol-1 for (H2O)4 and (HCl)3; roughly 30 kcal mol-1 for (H2O)3, (H2O)5 and (HCl)4; up to 38 kcal mol-1 for (HCl)5). At the CBS limit, MP2 significantly underestimates the CCSD(T) barrier heights (e.g., by ca. 2, 4 and 7 kcal mol-1 for the pentamers of HF, H2O and HCl, respectively), whereas CCSD overestimates these barriers by roughly the same magnitude. Scaling the barrier heights and dissociation energies by the number of fragments in the cluster reveals strong linear relationships between the two quantities and with the magnitudes of the imaginary vibrational frequency for the TSs.

Insight into the Binding of Argon to Cyclic Water Clusters from Symmetry-Adapted Perturbation Theory.

This work systematically examines the interactions between a single argon atom and the edges and faces of cyclic H2O clusters containing three-five water molecules (Ar(H2O)n=3-5). Full geometry optimizations and subsequent harmonic vibrational frequency computations were performed using MP2 with a triple-ζ correlation consistent basis set augmented with diffuse functions on the heavy atoms (cc-pVTZ for H and aug-cc-pVTZ for O and Ar; denoted as haTZ). Optimized structures and harmonic vibrational frequencies were also obtained with the two-body-many-body (2b:Mb) and three-body-many-body (3b:Mb) techniques; here, high-level CCSD(T) computations capture up through the two-body or three-body contributions from the many-body expansion, respectively, while less demanding MP2 computations recover all higher-order contributions. Five unique stationary points have been identified in which Ar binds to the cyclic water trimer, along with four for (H2O)4 and three for (H2O)5. To the best of our knowledge, eleven of these twelve structures have been characterized here for the first time. Ar consistently binds more strongly to the faces than the edges of the cyclic (H2O)n clusters, by as much as a factor of two. The 3b:Mb electronic energies computed with the haTZ basis set indicate that Ar binds to the faces of the water clusters by at least 3 kJ mol-1 and by nearly 6 kJ mol-1 for one Ar(H2O)5 complex. An analysis of the interaction energies for the different binding motifs based on symmetry-adapted perturbation theory (SAPT) indicates that dispersion interactions are primarily responsible for the observed trends. The binding of a single Ar atom to a face of these cyclic water clusters can induce perturbations to the harmonic vibrational frequencies on the order of 5 cm-1 for some hydrogen-bonded OH stretching frequencies.

Competition between Solvent···Solvent and Solvent···Solute Interactions in the Microhydration of the Tetrafluoroborate Anion, BF4-(H2O)n=1,2,3,4.

This study systematically examines the interactions of the tetrafluoroborate anion (BF4-) with up to four water molecules (BF4-(H2O)n=1,2,3,4). Full geometry optimizations and subsequent harmonic vibrational frequency computations are performed using a variety of density functional theory (DFT) methods (B3LYP, B3LYP-D3BJ, and M06-2X) and the MP2 ab initio method with a triple-ζ correlation consistent basis set augmented with diffuse functions on all non-hydrogen atoms (cc-pVTZ for H and aug-cc-pVTZ for B, O, and F; denoted as haTZ). Optimized structures and harmonic vibrational frequencies were also obtained with the CCSD(T) ab initio method and the haTZ basis set for the mono- and dihydrate (n = 1, 2) structures. The 2-body:Many-body (2b:Mb) technique, in which CCSD(T) computations capture the 1- and 2-body contributions to the interactions and MP2 computations recover all higher-order contributions, was used to extend these demanding computations to the tri- and tetrahydrate (n = 3, 4) systems. Four, five, and eight new stationary points have been identified for the di-, tri-, and tetrahydrate systems, respectively. The global minimum of the monohydrate adopts a symmetric double ionic hydrogen bond motif with C2v symmetry and an electronic dissociation energy of 13.17 kcal mol-1 at the CCSD(T)/haTZ level of theory. This strong solvent···solute interaction, however, competes with solute···solute interactions in the lowest-energy BF4-(H2O)n=2,3,4 minima that are not seen in the other di-, tri-, or tetrahydrate minima. The latter interactions help increase the 2b:Mb dissociation energies to more than 26, 41, and 51 kcal mol-1 for n = 2, 3, and 4, respectively. Structures that form hydrogen bonds between the solvating water molecules also exhibit the largest shifts in the harmonic OH stretching frequencies for the waters of hydration. These shifts can exceed -280 cm-1 relative to an isolated H2O molecule at the 2b:Mb/haTZ level of theory.

Characterization of Competing Halogen- and Hydrogen-Bonding Motifs in Simple Mixed Dimers of HCN and HX (X = F, Cl, Br, and I).

This work performs the first systematic comparison of hydrogen- and halogen-bonded configurations of the HCN/HX mixed dimer, where X = F, Cl, Br, and I. Eleven different minima have been characterized for these four heterogeneous dimers near the CCSD(T) complete basis set (CBS) limit. For each complex, two different hydrogen-bonded minima were identified: the global minimum where HX acts as the hydrogen bond donor and a local minimum where HX acts as the hydrogen bond acceptor. A halogen-bonded local minimum was also identified for all but the fluorine mixed dimer. To the best of our knowledge, three of the minima are identified here for the first time. The hydrogen- and halogen-bonded local minima of each complex become more energetically competitive with the global minimum as the atomic radius of the halogen atom increases. CCSD(T) relative energies of the hydrogen-bonded local minima computed near the CBS limit decrease from 4.5 kcal mol-1 for HCN/HF to 2.9, 2.4, and 1.2 kcal mol-1 for X = Cl, Br, and I, respectively. Corresponding relative energies for the halogen-bonded local minima range from 4.0 kcal mol-1 for X = Cl to 2.7 kcal mol-1 for X = Br and to as little as 0.5 kcal mol-1 X = I. Harmonic vibrational frequency shifts reported here suggest that it may be feasible to differentiate between the various minima for X = Cl, Br, and I via spectroscopic analysis, as was the case for the HCN/HF dimer.

Solvation of Isoelectronic Halide and Alkali Metal Ions by Argon Atoms.

This work systematically examines the interactions of alkali metal cations and their isoelectronic halide counterparts with up to six solvating Ar atoms (M+Arn and X-Arn, where M = Li, Na, K, and Rb; X = H, F, Cl, and Br; and n = 1-6) via full geometry optimizations with the MP2 method and robust, correlation-consistent quadruple-ζ (QZ) basis sets. 116 unique M+Arn and X-Arn stationary points have been characterized on the MP2/QZ potential energy surface. To the best of our knowledge, approximately two dozen of these stationary points have been reported here for the first time. Some of these new structures are either the lowest-energy stationary point for a particular cluster or energetically competitive with it. The CCSD(T) method was employed to perform additional single-point energy computations upon all MP2/QZ-optimized structures using the same basis set. CCSD(T)/QZ results indicate that internally solvated structures with the ion at/near the geometric center of the cluster have appreciably higher energies than those placing the ion on the periphery. While this study extends the prior investigations of M+Arn clusters found within the literature, it notably provides one of the first thorough characterizations of and comparisons to the corresponding negatively charged X-Arn clusters.

Relative energetics of CH3CH2O, CH3CHOH, and CH2CH2OH radical products from ethanol dehydrogenation.

This study has examined the relative energetics of nine stationary points associated with the three different radical isomers generated by removing a H atom from ethanol at the O atom (ethoxy, CH3CH2O), the α C atom (CH3CHOH), and the ß C atom (CH2CH2OH). For the first time, CCSD(T) geometry optimizations and harmonic vibrational frequency computations with the cc-pVTZ and aug-cc-pVTZ basis sets have been carried out to characterize two unique minima for each isomer along with three transition state structures with Cs symmetry. Explicitly correlated CCSD(T) computations were also performed to estimate the relative energetics of these nine stationary points near the complete basis set limit. These benchmark results were used to assess the performance of various density functional theory (DFT) and wave function theory methods, and they will help guide method selection for future studies of alcohols and their radicals. The structures generated by abstracting H from the α C atom have significantly lower electronic energies (by at least 7 kcal mol-1) than the CH3CH2O and CH2CH2OH radicals. Although previously reported as a minimum on the ground-state surface, the 2A″ Cs structure of the ethoxy radical was found to be a transition state in this study with MP2, CCSD(T), and a number of DFT methods. An implicit solvation model used in conjunction with DFT and MP2 methods did not qualitatively change the relative energies of the isomers, but the results suggest that the local minima for the CH3CHOH and CH2CH2OH radicals could become more energetically competitive in condensed phase environments, such as liquid water and ethanol.

Torsional Profiles of Thiophene and Furan Oligomers: Probing the Effects of Heterogeneity and Chain Length.

A systematic analysis of the torsional profiles of 55 unique oligomers composed of two to four thiophene and/or furan rings (n = 2 to 4) has been conducted using three density functional theory (DFT) methods along with MP2 and three different coupled-cluster methods. Two planar or quasi-planar minima were identified for each n = 2 oligomer system. In every case, the torsional angle (τ) between the heteroatoms about the carbon-carbon bond connecting the two rings is at or near 180° for the global minimum and 0° for the local minimum, referred to as anti and syn conformations, respectively. These oligomers have rotational barrier heights ranging from ca. 2 kcal mol-1 for 2,2'-bithiophene to 4 kcal mol-1 for 2,2'-bifuran, based on electronic energies computed near the CCSD(T) complete basis set (CBS) limit. The corresponding rotational barrier for the heterogeneous 2-(2-thienyl)furan counterpart falls approximately halfway between those values. The energy differences between the minima are approximately 2 and 0.4 kcal mol-1 for the homogeneous 2,2'-bifuran and 2,2'-bithiophene, respectively, whereas the energy difference between the planar local and global minima (at τ = 0 and 180°, respectively) is only 0.3 kcal mol-1 for 2-(2-thienyl)furan. Extending these three oligomers by adding one or two additional thiophene and/or furan rings resulted in only minor changes to the torsional profiles when rotating around the same carbon-carbon bond as the two-ring profiles. Relative energy differences between the syn and anti conformations were changed by no more than 0.4 kcal mol-1 for the corresponding n = 3 and 4 oligomers, while the rotational barrier height increased by no more than 0.8 kcal mol-1.

Anharmonic vibrational frequencies of ammonia borane (BH3NH3).

The fundamental vibrational frequency of the B-N stretch in BH3NH3 has eluded gas-phase experimental observation for decades. This work offers a theoretical anharmonic prediction of this mode to be 644 cm-1, using a Cartesian quartic force field at the CCSD(T)-F12/cc-pVTZ-F12 level of theory. The other fundamental frequencies reported herein have a mean absolute error of only 5 cm-1 from the seven available gas-phase experimental frequencies, making the anharmonic vibrational frequencies and rotational constants the most accurate computational data available for BH3NH3 to date. The inclusion of Fermi, Coriolis, and Darling-Dennison resonances is a major source of this accuracy, with the non-resonance-corrected frequencies having a mean absolute error of 10 cm-1. In particular, the inclusion of the 2ν6 = ν5 type 1 Fermi resonance increases the B-N stretching frequency by 14 cm-1 compared to previous work. Ammonia borane also represents one of the largest molecules ever studied by quartic force fields, making this work an important step in extending the breadth of application for these theoretical rovibrational techniques.

Competition between Solvent-Solvent and Solvent-Solute Interactions in the Microhydration of the Hexafluorophosphate Anion, PF6-(H2O)n=1,2.

This study systematically examines the interactions of the hexafluorophosphate anion (PF6-) with one or two solvent water molecules (PF6-(H2O)n where n = 1, 2). Full geometry optimizations and subsequent harmonic vibrational frequency computations are performed on each stationary point using a variety of common density functional theory methods (B3LYP, B3LYP-D3, M06-2X, and ωB97XD) and the MP2 and CCSD(T) ab initio methods with a triple-ζ correlation consistent basis set augmented with diffuse functions on all non-hydrogen atoms (cc-pVTZ for H and aug-cc-pVTZ for P, O, and F; denoted as haTZ). Five new stationary points of PF6-(H2O)2 have been identified, one of which has an electronic energy of approximately 2 kcal mol-1 lower than the only other dihydrate structure reported for this system. The CCSD(T) computations also reveal that the detailed interactions between PF6- and H2O can be quite difficult to model reliably, with some methods struggling to correctly characterize stationary points for n = 1 or accurately reproduce the vibrational frequency shifts induced by the formation of the hydrated complex. Although the interactions between the solvent and ionic solute are quite strong (CCSD(T) electronic dissociation energy ≈10 kcal mol-1 for the monohydrate minimum), the solvent-solvent interactions in the lowest-energy PF6-(H2O)2 minimum give rise to appreciable cooperative effects not observed in the other dihydrate minima. In addition, this newly identified structure exhibits the largest frequency shifts in the OH stretching vibrations for the waters of hydration (with Δω exceeding -100 cm-1 relative to the values for an isolated H2O molecule).

Anchoring the hydrogen sulfide dimer potential energy surface to juxtapose (H2S)2 with (H2O)2.

Twelve stationary points have been characterized on the (H2S)2 potential energy surface using the MP2 and CCSD(T) methods with large, correlation consistent basis sets. To the best of our knowledge, five of the structures have not been identified elsewhere and are presented here for the first time. A similar analysis was performed on the ten, well-known structures of the water dimer in order to facilitate direct comparisons between the corresponding (H2O)2 and (H2S)2 configurations. Harmonic vibrational frequency computations identify three (H2S)2 configurations as minima, four as transition states, and five as higher-order saddle points (ni = 0, ni = 1, and ni ≥ 2, respectively, where ni is the number of imaginary frequencies). The two local minima and four transition state structures identified have electronic energies within 0.73 kJ mol-1 of the global minimum near the CCSD(T) complete basis set (CBS) limit, and the five higher-order saddle points range from 1.90 kJ mol-1 to 4.31 kJ mol-1 above the global minimum at the same level of theory. One of the more substantial differences observed between the H2S and H2O systems is that (H2O)2 has only a single minimum, while the other nine stationary points are significantly higher in energy ranging from 2.15 kJ mol-1 to 14.89 kJ mol-1 above the global minimum near the CCSD(T) CBS limit. For (H2S)2, the electronic dissociation energy of the global minimum is only 7.02 kJ mol-1 at the CCSD(T) CBS limit, approximately three times smaller than the dissociation energy of (H2O)2.

Observation of the Low-Frequency Spectrum of the Water Trimer as a Sensitive Test of the Water-Trimer Potential and the Dipole-Moment Surface.

Intermolecular interactions in bulk water are dominated by pairwise and non-pairwise cooperative interactions. While accurate descriptions of the pairwise interactions are available and can be tested by precise low-frequency spectra of the water dimer up to 550 cm-1 , the same does not hold for the three-body interactions. Here, we report the first comprehensive spectrum of the water trimer in the frequency region from 70 to 620 cm-1 using helium-nanodroplet isolation and free-electron lasers. By comparison to accompanying high-level quantum calculations, the experimentally observed intermolecular bands can be assigned. The transition frequencies of the degenerate translation, the degenerate in-plane and the non-degenerate out-of-plane libration, as well as additional bands of the out-of-plane librational mode are reported for the first time. These provide a benchmark for state-of-the-art water potentials and dipole-moment surfaces, especially with respect to three-body interactions.

Co-Localization of DNA i-Motif-Forming Sequences and 5-Hydroxymethyl-cytosines in Human Embryonic Stem Cells.

G-quadruplexes (G4s) and i-motifs (iMs) are tetraplex DNA structures. Sequences capable of forming G4/iMs are abundant near the transcription start sites (TSS) of several genes. G4/iMs affect gene expression in vitro. Depending on the gene, the presence of G4/iMs can enhance or suppress expression, making it challenging to discern the underlying mechanism by which they operate. Factors affecting G4/iM structures can provide additional insight into their mechanism of regulation. One such factor is epigenetic modification. The 5-hydroxymethylated cytosines (5hmCs) are epigenetic modifications that occur abundantly in human embryonic stem cells (hESC). The 5hmCs, like G4/iMs, are known to participate in gene regulation and are also enriched near the TSS. We investigated genomic co-localization to assess the possibility that these two elements may play an interdependent role in regulating genes in hESC. Our results indicate that amongst 15,760 G4/iM-forming locations, only 15% have 5hmCs associated with them. A detailed analysis of G4/iM-forming locations enriched in 5hmC indicates that most of these locations are in genes that are associated with cell differentiation, proliferation, apoptosis and embryogenesis. The library generated from our analysis is an important resource for investigators exploring the interdependence of these DNA features in regulating expression of selected genes in hESC.

Cis/Trans Energetics in Epoxide, Thiirane, Aziridine and Phosphirane Containing Cyclopentanols: Effects of Intramolecular OH⋯O, S, N and P Contacts.

A recent computational analysis of the stabilizing intramolecular OH⋯O contact in 1,2-dialkyl-2,3-epoxycyclopentanol diastereomers has been extended to thiiriane, aziridine and phosphirane analogues. Density functional theory (DFT), second-order Møller-Plesset perturbation theory (MP2) and CCSD(T) coupled-cluster computations with simple methyl and ethyl substituents indicate that electronic energies of the c i s isomers are lowered by roughly 3 to 4 kcal mol-1 when the OH group of these cyclopentanol systems forms an intramolecular contact with the O, S, N or P atom on the adjacent carbon. The results also suggest that S and P can participate in these stabilizing intramolecular interactions as effectively as O and N in constrained molecular environments. The stabilizing intramolecular OH⋯O, OH⋯S, OH⋯N and OH⋯P contacts also increase the covalent OH bond length and significantly decrease the OH stretching vibrational frequency in every system with shifts typically on the order of -41 cm-1.

Computational Investigation on Electronic Structures and Properties of 4,6-Bis(nitroimino)-1,3,5-triazinan-2-one: An Insensitive Munition Compound.

In order to minimize unintentional detonation, munitions researchers have focused on the development of chemical compounds that are insensitive to external stimuli while maintaining their effectiveness. Although these compounds, known as high-performance insensitive munition compounds, are promising in terms of potency and stability, their environmental impacts have either not been fully understood or are yet to be investigated. In the present research, we have performed a quantum chemical investigation on electronic structures and properties of an insensitive munition compound 4,6-bis(nitroimino)-1,3,5-triazinan-2-one (DNAM). The density functional theory using the B3LYP and M06-2X functionals and MP2 methodology were used for geometry optimization of various tautomeric forms of DNAM. The effect of bulk water solution was evaluated using the conductor-like polarizable continuum model and the density-based solvation model. Ionization potentials, electron affinities, redox properties, and p Ka values were also computed and compared with the available experimental data. These physical and chemical properties of DNAM have been discussed with regard to the varying tautomeric forms in which DNAM can exist.

Examination of the structures, energetics, and vibrational frequencies of small sulfur-containing prototypical dimers, (H2 S)2 and H2 O/H2 S.

The optimized geometries, vibrational frequencies, and dissociation energies from MP2 and CCSD(T) computations with large correlation consistent basis sets are reported for (H2 S)2 and H2 O/H2 S. Anharmonic vibrational frequencies have also been computed with second-order vibrational perturbation theory (VPT2). As such, the fundamental frequencies, overtones, and combination bands reported in this study should also provide a useful road map for future spectroscopic studies of the simple but important heterogeneous H2 O/H2 S dimer in which the hydrogen bond donor and acceptor can interchange, leading to two unique minima with very similar energies. Near the CCSD(T) complete basis set limit, the HOH⋯SH2 configuration (H2 O donor) lies only 0.2 kcal mol-1 below the HSH⋯OH2 structure (H2 S donor). When the zero-point vibrational energy is included, however, the latter configuration becomes slightly lower in energy than the former by <0.1 kcal mol-1 . © 2018 Wiley Periodicals, Inc.

Communication: Gas phase vibrational spectroscopy of the azide-water complex.

The vibrational spectra of the azide-water complex, N3 -(H2O), and its fully deuterated isotopologue are studied using infrared photodissociation (IRPD) spectroscopy (800-3800 cm-1) and high-level ab initio computations. The IRPD spectrum of the H2-tagged complex exhibits four fundamental transitions at 3705, 3084, 2003, and 1660 cm-1, which are assigned to the free OH stretching, the hydrogen-bonded O-H stretching, the antisymmetric N3 stretching, and the water bending mode, respectively. The IRPD spectrum is consistent with a planar, singly hydrogen-bonded structure according to an MP2 and CCSD(T) anharmonic analysis via generalized second-order vibrational perturbation theory. The red-shift of the hydrogen-bonded OH stretching fundamental of 623 cm-1 associated with this structure is computed within 6 cm-1 (or 1%) and is used to estimate the proton affinity of azide (1410 kJ mol-1). Born-Oppenheimer molecular dynamics simulations show that large amplitude motions are responsible for the observed band broadening at cryogenic temperature. Temperature-dependent (6-300 K) IR multiphoton dissociation spectra of the untagged complex are also presented and discussed in the context of spectral diffusion observed in the condensed phase.

Binding of the atomic cations hydrogen through argon to water and hydrogen sulfide.

Water and hydrogen sulfide will bind with every atomic cation from the first three rows of the periodic table. While some atoms bind more tightly than others, explicitly correlated coupled cluster theory computations show that energy is required to be put into the system in order to dissociate these bonds even for noble gas atoms. The most promising systems have shallow entrance potential energy surfaces (PESs) that lie above deeper wells of a different spin. These wells are shown explicitly for H2OO+, H2SS+, and H2OS+ where relaxed PESs of the heavy atom bond lengths indicate that quartet states will cross more deeply-bound doublet states allowing for relatively easy association but much more difficult dissociation. In astrophysical regions that are cold and diffuse, such associations could lead to the formation of novel molecules utilizing water (or H2S) as the building blocks of more rich subsequent chemistry. Recent work has hypothesized that oxywater (H2OO) may be an intermediate in the formation of molecular oxygen in comets, and this work supports such a conclusion at least from a molecular cation perspective.