skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Formic Acid–Ammonia Heterodimer: A New Δ-Machine Learning CCSD(T)-Level Potential Energy Surface Allows Investigation of the Double Proton Transfer
The formic acid-ammonia dimer is an important example of a hydrogen-bonded complex in which a double proton transfer can occur. Its microwave spectrum has recently been reported and rotational constants and quadrupole coupling constants were determined. Calculated estimates of the double-well barrier and the internal barriers to rotation were also reported. Here we report a full-dimensional potential energy surface (PES) for this complex, using two closely related Δ-machine learning methods to bring it to the CCSD(T) level of accuracy. The PES dissociates smoothly and accurately. Using a 2d quantum model the ground vibrational-state tunneling splitting is estimted to be less than 10−4 cm−1. The dipole moment along the intrinsic reaction coordinate is calculated along with a Mullikan charge analysis and supports mildly ionic character of the minimum and strongly ionic character at the double-well barrier.  more » « less
Award ID(s):
1952289
PAR ID:
10518570
Author(s) / Creator(s):
; ; ; ; ; ; ;
Editor(s):
Gagliardi, Laura
Publisher / Repository:
American Chemical Society
Date Published:
Journal Name:
Journal of Chemical Theory and Computation
Edition / Version:
1
Volume:
20
Issue:
5
ISSN:
1549-9618
Page Range / eLocation ID:
1821 to 1828
Subject(s) / Keyword(s):
Energy Mathematical methods' Polarity' Potential energy Tunneling
Format(s):
Medium: X Size: 2.5Mb Other: pdf
Size(s):
2.5Mb
Sponsoring Org:
National Science Foundation
More Like this
  1. Heaven, Michael (Ed.)
    New analysis and spectra are reported for the gas-phase ammonia-formic acid complex. Calculations to deter- mine the theoretical barrier to internal rotation were conducted and led to the new internal rotation analysis of the dimer. Using the new analysis and calculations, 12 new lines were measured and assigned and included in the present analysis. This is the !rst internal rotation analysis for this complex. The measurements were made in the 7–22 GHz range using two Flygare-Balle type pulsed beam Fourier transform microwave (FTMW) spectrometers. The complex was analyzed as a hindered rotor and 20 A and 16 E state transitions were !t with the XIAM5 program. The rotational constants were determined to have the following values: A = 11970.19(9) MHz, B = 4331.479(4) MHz, and C = 3227.144(4) MHz. Rotational constants, quadrupole coupling constants, and internal rotor parameters were !t to the spectrum. Double resonance was used to verify line assignments and access higher frequencies. The barrier to internal rotation was found to be 195.18(7) cm 1. High level calculations are in good agreement with experimental values. The calculated V3 barrier values range from 168.3 to 212.8 cm 1. 
    more » « less
  2. In this work the H2O–HCN complex is quantitatively characterized in two ways. First, we report a new rigid-monomer 5D intermolecular potential energy surface (PES) for this complex, calculated using the symmetry-adapted perturbation theory based on density functional theory method. The PES is based on 2833 ab initio points computed employing the aug-cc-pVQZ basis set, utilizing the autoPES code, which provides a site-site analytical fit with the long-range region given by perturbation theory. Next, we present the results of the quantum 5D calculations of the fully coupled intermolecular rovibrational states of the H2O–HCN complex for the total angular momentum J values of 0, 1, and 2, performed on the new PES. These calculations rely on the quantum bound-state methodology developed by us recently and applied to a variety of noncovalently bound binary molecular complexes. The vibrationally averaged ground-state geometry of H2O–HCN determined from the quantum 5D calculations agrees very well with that from the microwave spectroscopic measurements. In addition, the computed ground-state rotational transition frequencies, as well as the B and C rotational constants calculated for the ground state of the complex, are in excellent agreement with the experimental values. The assignment of the calculated intermolecular vibrational states of the H2O–HCN complex is surprisingly challenging. It turns out that only the excitations of the intermolecular stretch mode can be assigned with confidence. The coupling among the angular degrees of freedom (DOFs) of the complex is unusually strong, and as a result most of the excited intermolecular states are unassigned. On the other hand, the coupling of the radial, intermolecular stretch mode and the angular DOFs is weak, allowing straightforward assignment of the excitation of the former. 
    more » « less
  3. Abstract Potential energy surface (PES) analyses at the SMD[MP2/6–311++G(d,p)] level and higher‐level energies up to MP4(fc,SDTQ) are reported for the fluorinated tertiary carbamateN‐ethyl‐N‐(2,2,2‐trifluoroethyl) methyl carbamate (VII) and its parent systemN,N‐dimethyl methyl carbamate (VI). Emphasis is placed on the analysis of the rotational barrier about the CN carbamate bond and its interplay with the hybridization of theN‐lone pair (NLP). All rotational transition state (TS) structures were found by computation of 1D relaxed rotational profiles but only 2D PES scans revealed the rotation‐inversion paths in a compelling fashion. We found four unique chiral minima ofVII, one pair each ofE‐andZ‐rotamers, and we determined theeightunique rotational TS structures associated with every possibleE/Z‐isomerization path. It is a significant finding that all TS structures featureN‐pyramidalization whereas the minima essentially contain sp2‐hybridized nitrogen. We will show that the TS stabilities are affected by the synergetic interplay between NLP/CO2repulsion minimization, NLP→σ*(CO) negative hyperconjugation, and two modes of intramolecular through‐space electrostatic stabilization. We demonstrate how Boltzmann statistics must be applied to determine the predicted experimental rotational barrier based on the energetics of all eight rotamerization pathways. The computed barrier forVIIis in complete agreement with the experimentally measured barrier of the very similar fluorinated carbamateN‐Boc‐N‐(2,2,2‐trifluoroethyl)‐4‐aminobutan‐1‐olII. NMR properties ofVIIwere calculated with a variety of density functional/basis set combinations and Boltzmann averaging over theE‐andZ‐rotamers at our best theoretical level results in good agreement with experimental chemical shifts δ(13C) andJ(13C,19F) coupling constants ofII(within 6 %). 
    more » « less
  4. Abstract A typical task for classical and quantum computing in chemistry is finding a potential energy surface (PES) along a reaction coordinate, which involves solving the quantum chemistry problem for many points along the reaction path. Developing algorithms to accomplish this task on quantum computers has been an active area of development, yet finding all the relevant eigenstates along the reaction coordinate remains a difficult problem, and determining PESs is thus a costly proposal. In this paper, we demonstrate the use of a eigenvector continuation—a subspace expansion that uses a few eigenstates as a basis—as a tool for rapidly exploring PESs. We apply this to determining the binding PES or torsion PES for several molecules of varying complexity. In all cases, we show that the PES can be captured using relatively few basis states; suggesting that a significant amount of (quantum) computational effort can be saved by making use of already calculated ground states in this manner. 
    more » « less
  5. Chiral transition-metal complexes are of interest in many fields ranging from asymmetric catalysis and molecular materials science to optoelectronic applications or fundamental physics including parity violation effects. We present here a combined theoretical and experimental investigation of gas-phase valence-shell photoelectron circular dichroism (PECD) on the challenging open-shell ruthenium( iii )-tris-(acetylacetonato) complex, Ru(acac) 3 . Enantiomerically pure Δ- or Λ-Ru(acac) 3 , characterized by electronic circular dichroism (ECD), were vaporized and adiabatically expanded to produce a supersonic beam and photoionized by circularly-polarized VUV light from the DESIRS beamline at Synchrotron SOLEIL. Photoelectron spectroscopy (PES) and PECD experiments were conducted using a double imaging electron/ion coincidence spectrometer, and compared to density functional theory (DFT) and time-dependent DFT (TDDFT) calculations. The open-shell character of Ru(acac) 3 , which is not taken into account in our DFT approach, is expected to give rise to a wide multiplet structure, which is not resolved in our PES signals but whose presence might be inferred from the additional striking features observed in the PECD curves. Nevertheless, the DFT-based assignment of the electronic bands leads to the characterisation of the ionized orbitals. In line with other recent works, the results confirm that PECD persists independently on the localization and/or on the achiral or chiral nature of the initial orbital, but is rather a probe of the molecular potential as a whole. Overall, the measured PECD signals on Ru(acac) 3 , a system exhibiting D 3 propeller-type chirality, are of similar magnitude compared to those on asymmetric-carbon-based chiral organic molecules which constitute the vast majority of species investigated so far, thus suggesting that PECD is a universal mechanism, inherent to any type of chirality. 
    more » « less