The sulfur dioxide-water complex: CCSD(T)/CBS anharmonic vibrational spectroscopy of stacked and hydrogen-bonded dimers.

This study examines the structures, energies, and IR vibrational spectra of the sulfur dioxide-water SO2(H2O) complexes by employing coupled cluster theory CCSD(T) with Dunning style correlation consistent type basis sets aug-cc-pV(n+d)Z (n = D, T, Q, 5). Complete basis set (CBS) extrapolations have been carried out to predict binding energies for two isomers of the SO2(H2O) complex: a stacked global minimum (1A) structure and a hydrogen-bonded local minimum (1B) structure. The CCSD(T)/CBS extrapolation predicts an intermolecular S-O distance rS⋯O = 2.827 Å for the stacked isomer, which is in excellent agreement with an experimental measurement of 2.824 Å [K. Matsumura et al., J. Chem. Phys., 91, 5887 (1989)]. The CCSD(T)/CBS binding energy for the stacked dimer 1A and hydrogen-bonded form 1B is De = -4.37 kcal/mol and De = -2.40 kcal/mol, respectively. This study also employs anharmonic VPT2 MP2/aug-cc-pV(n+d)Z level corrections to CCSD(T)/aug-cc-pV(n+d)Z vibrational frequencies in both forms of SO2(H2O). The anharmonic CCSD(T)/aug-cc-pV(Q+d)Z OH stretching frequencies in the stacked structure 1A are 3743 cm-1 (ν3) and 3647 cm-1 (ν1), and these align well with the recorded IR spectroscopic values of 3745 and 3643 cm-1, respectively [C. Wang et al., J. Phys. Chem. Lett., 13, 5654 (2022)]. If we combine CCSD(T)/aug-cc-pV(n+d)Z De values with VPT2 vibrational frequencies, we obtain a new CCSD(T)/aug-cc-pV(Q+d)Z anharmonic dissociation energy D0 = -3.48 kcal/mol for 1A and D0 = -1.74 kcal/mol for 1B. In summary, the results presented here demonstrate that the application of CCSD(T) calculations with aug-cc-pV(n+d)Z basis sets and CBS extrapolations is critical in probing the structure and IR spectroscopic properties of the sulfur dioxide-water complex.

Structure and solvation dynamics of the hydroxide ion in ice-like water clusters: a CCSD(T) and car-parrinello molecular dynamics study.

Using MP2, CCSD(T) electronic structure theory and ab initio molecular dynamics simulations, we explore the structure, solvation dynamics and vibrational spectra of OH-(H2O)n clusters. Our study reports new cubic and fused cubic global minima structures of OH-(H2O)n for n = 8-26 with surface and interior solvation arrangements. In the case of OH-(H2O)26, we show that MP2 and CCSD(T) calculations predict global minima structures with the hydroxide ion occupying the interior region of a densely packed cubic cluster that is secured by ionic hydrogen bonds. More importantly, results from ab initio molecular dynamics simulations of OH-(H2O)26 demonstrate that the hydroxide ion remains in the cluster interior and hexa-coordinated, irrespective of the temperature, up to around 175 K, then incrementally transitions from a surface-exposed penta- (170-200 K), to a tetra- (225 K) to a tri-coordinated OH-(H2O)3 structure at 300 K. Building on our temperature-dependent vibrational power spectra, we are also able to disentangle structure and temperature effects on individual spectral contributions arising from water molecules located in the inner and outer shell of OH-(H2O)26. Some of these theoretical results provide valuable guidance for the interpretation of IRMPD spectra of small hydroxide-water clusters, but there are also several intriguing implications of these results, in particular, for the solvation of the OH- ion at the surface of water nanodroplets and aqueous interfaces.

The ozone-water complex: CCSD(T)/CBS structures and anharmonic vibrational spectroscopy of O3(H2O)n, (n = 1 - 2).

Ozone-water complexes O3(H2O)n (n = 1-2) have been studied using coupled cluster theory with triple excitations CCSD(T) with correlation consistent basis sets aug-cc-pVnZ (n = D, T, Q) and complete basis set (CBS) extrapolation techniques. We identified seven dimer (n = 1) and nine trimer species (n = 2) with open C2v and cyclic D3h ozone. Calculations at the CCSD(T)/CBS level of theory for C2v O3(H2O) on the counterpoise (CP)-corrected potential energy surface yield a dissociation energy of De = 2.31 kcal/mol and an O3 central-oxygen (Oc) H2O oxygen (Ow) distance r[Oc⋯Ow] of 3.097 Å, which is in good agreement with an experimental value of 2.957 Å [J. Z. Gillies et al., J. Mol. Spectrosc. 146, 493 (1991)]. Combining our CCSD(T)/CBS value of De for C2v O3(H2O) with our best estimate anharmonic CCSD(T)/aVTZ ΔZPE yields a Do value of 1.82 kcal/mol; the CCSD(T)/CBS value of De for D3h O3(H2O) is 1.51 kcal/mol and yields an anharmonic CCSD(T)/aVTZ Do = 0.99 kcal/mol. CCSD(T)/aVTZ dissociation energies and structures for C2v O3(H2O)2 are De = 4.15 kcal/mol, (Do = 3.08 kcal/mol) and r[Oc⋯Ow] = 2.973 Å, and De = 2.64 kcal/mol (Do = 1.68 kcal/mol) with r[Oc⋯Ow] = 2.828 Å for D3h O3(H2O)2. The results from ab initio molecular dynamics simulations, which consider dynamic and thermal effects in O3(H2O), show that the O3(H2O) complex remains stable at 50 K and dynamically interconverts between two hydrogen-bonded conformers with short Oc⋯Ow contacts (3.85 Å). Carr-Parrinello molecular dynamic (CPMD) simulations for O3(H2O) and O3(H2O)2 at 100 K demonstrate that O3(H2O)2 remains structurally intact, whereas O3(H2O) dissociates to free ozone and water, a feature consistent with the larger average binding energy in O3(H2O)2 (2.2 kcal/mol) vs that in O3(H2O) (1.8 kcal/mol). Finally, the results from CCSD(T)/CBS and CPMD simulations demonstrate that the large inter-trimer binding energies in O3(H2O)2 would give rise to an elevated trimer/dimer population ratio, making O3(H2O)2 a particularly stable and spectroscopically detectable complex.

Structure and binding energy of the H2S dimer at the CCSD(T) complete basis set limit.

This study presents results for the binding energy and geometry of the H2S dimer which have been computed using Møller-Plesset perturbation theory (MP2, MP4) and coupled cluster (CCSD, CCSD(T)) calculations with basis sets up to aug-cc-pV5Z. Estimates of De, EZPE, Do, and dimer geometry have been obtained at each level of theory by taking advantage of the systematic convergence behavior toward the complete basis set (CBS) limit. The CBS limit binding energy values of De are 1.91 (MP2), 1.75 (MP4), 1.41 (CCSD), and 1.69 kcal/mol (CCSD[T]). The most accurate values for the equilibrium S-S distance rSS (without counterpoise correction) are 4.080 (MP2/aug-cc-pV5Z), 4.131 (MP4/aug-cc-pVQZ), 4.225 (CCSD/aug-cc-pVQZ), and 4.146 Å (CCSD(T)/aug-cc-pVQZ). This study also evaluates the effect of counterpoise correction on the H2S dimer geometry and binding energy. As regards the structure of (H2S)2, MPn, CCSD, and CCSD(T) level values of rSS, obtained by performing geometry optimizations on the counterpoise-corrected potential energy surface, converge systematically to CBS limit values of 4.099 (MP2), 4.146 (MP4), 4.233 (CCSD), and 4.167 Å (CCSD(T)). The corresponding CBS limit values of the equilibrium binding energy De are 1.88 (MP2), 1.76 (MP4), 1.41 (CCSD), and 1.69 kcal/mol (CCSD(T)), the latter in excellent agreement with the measured binding energy value of 1.68 ± 0.02 kcal/mol reported by Ciaffoni et al. [Appl. Phys. B 92, 627 (2008)]. Combining CBS electronic binding energies De with EZPE predicted by CCSD(T) vibrational second-order perturbation theory calculations yields Do = 1.08 kcal/mol, which is around 0.6 kcal/mol smaller than the measured value of 1.7 ± 0.3 kcal/mol. Overall, the results presented here demonstrate that the application of high level calculations, in particular CCSD(T), in combination with augmented correlation consistent basis sets provides valuable insight into the structure and energetics of the hydrogen sulfide dimer.

Gold chloride clusters with Au(III) and Au(I) probed by FT-ICR mass spectrometry and MP2 theory.

Microsolvated clusters of gold chloride are probed by electrospray ionization mass spectrometry (ESI-MS) and scalar relativistic electronic structure calculations. Electrospray ionization of aqueous AuCl3 leads to mononuclear clusters of types [AuCl2](+)(H2O)n (n = 0-4), [AuOHCl](+)(H2O)n (n = 0-1) and [AuCl2](+)(HCl)2(H2O)n (n = 0-4). In addition, strong ion signals due to dinuclear [Au2Cl5-xOHx](+)(H2O)n (x = 0-1) are present in ESI mass spectra of aqueous AuCl3, with the abundance of individual dinuclear species controlled by the concentration-dependent variation of the precursor complexes [AuCl2-xOHx](+)(H2O)n and AuCl3. Equilibrium structures, energies and thermodynamic properties of mono- and dinuclear gold clusters have been predicted using MP2 and CCSD(T) theory, and these data have been applied to examine the influence of microsolvation on cluster stability. Specifically, results from CCSD(T) calculations indicate that non-covalently bound ion-neutral complexes Au(+)(Cl2)(H2O)n, with formal Au(I), are the dominant forms of mononuclear gold with n = 0-2, while higher hydrates (n > 2) are covalently bound [AuCl2](+)(H2O)n complexes in which gold exists as Au(III). MP2 calculations show that the lowest energy structure of dinuclear gold is an ion-molecule cluster [Au2Cl(Cl2)2](+) consisting of a single-bridged digold-chloronium ion bound end-on to two dichlorine ligands, with two higher energy isomers, single-bridged [Au2Cl3(Cl2)](+) and double-bridged [Au2Cl5](+) clusters. Finally, AuAu interactions in the singly-bridged clusters [Au2Cl(Cl2)2](+)(H2O)n and [Au2Cl3(Cl2)](+)(H2O)n are examined employing a wide range of computational tools, including natural bond order (NBO) analysis and localized orbital locator (LOL) profiles.

Relativistic-consistent electron densities of the coinage metal clusters M2, M4, M4(2-), and M4Na2 (M = Cu, Ag, Au): a QTAIM study.

We employ second-order Møller-Plesset perturbation theory level in combination with recently developed pseudopotential-based correlation consistent basis sets to obtain accurate relativistic-consistent electron densities for small coinage metal clusters. Using calculated electron densities, we employ Bader's quantum theory of atoms in molecules (QTAIM) to gain insights into the nature of metal-metal bonding in the clusters M(2), M(4), M(4)(2-), and M(4)Na(2) (M = Cu, Ag, Au). For the simplest case of the metal dimer, M(2), we correlate the strength of the metal-metal bond with the value of the electron density at the bond critical point, the total energy density at the bond critical point, the sharing (delocalization) index, and the values of the two principle negative curvatures. We then consider changes to the metal-metal bonding and charge density distribution upon the addition of two metal atoms to form the metal tetramer, M(4), and then followed by the addition of an electron pair to form M(4)(2-) and finally followed by the addition of two alkali metal (sodium) ions to form M(4)Na(2). Using topological properties of the electron density, we present evidence for the existence of σ-aromaticity in Au(4)(2-). We also report the existence of two non-nuclear attractors in the molecular graph of Cu(4)(2-) and large negative charge accumulation in the nonbonded Cu basins of this cluster.

Peptide synthesis in early Earth hydrothermal systems.

We report here results from experiments and thermodynamic calculations that demonstrate a rapid, temperature-enhanced synthesis of oligopeptides from the condensation of aqueous glycine. Experiments were conducted in custom-made hydrothermal reactors, and organic compounds were characterized with ultraviolet-visible procedures. A comparison of peptide yields at 260 degrees C with those obtained at more moderate temperatures (160 degrees C) gives evidence of a significant (13 kJ . mol(-1)) exergonic shift. In contrast to previous hydrothermal studies, we demonstrate that peptide synthesis is favored in hydrothermal fluids and that rates of peptide hydrolysis are controlled by the stability of the parent amino acid, with a critical dependence on reactor surface composition. From our study, we predict that rapid recycling of product peptides from cool into near-supercritical fluids in mid-ocean ridge hydrothermal systems will enhance peptide chain elongation. It is anticipated that the abundant hydrothermal systems on early Earth could have provided a substantial source of biomolecules required for the origin of life.

Extreme accumulation of nucleotides in simulated hydrothermal pore systems.

We simulate molecular transport in elongated hydrothermal pore systems influenced by a thermal gradient. We find extreme accumulation of molecules in a wide variety of plugged pores. The mechanism is able to provide highly concentrated single nucleotides, suitable for operations of an RNA world at the origin of life. It is driven solely by the thermal gradient across a pore. On the one hand, the fluid is shuttled by thermal convection along the pore, whereas on the other hand, the molecules drift across the pore, driven by thermodiffusion. As a result, millimeter-sized pores accumulate even single nucleotides more than 10(8)-fold into micrometer-sized regions. The enhanced concentration of molecules is found in the bulk water near the closed bottom end of the pore. Because the accumulation depends exponentially on the pore length and temperature difference, it is considerably robust with respect to changes in the cleft geometry and the molecular dimensions. Whereas thin pores can concentrate only long polynucleotides, thicker pores accumulate short and long polynucleotides equally well and allow various molecular compositions. This setting also provides a temperature oscillation, shown previously to exponentially replicate DNA in the protein-assisted PCR. Our results indicate that, for life to evolve, complicated active membrane transport is not required for the initial steps. We find that interlinked mineral pores in a thermal gradient provide a compelling high-concentration starting point for the molecular evolution of life.