High resolution spectroscopy of carboxylic acid in the gas phase: Observation of proton transfer in (DCOOH)2

Hydrated Formic Acid Clusters and their Interaction with Electrons

We investigate ion formation in hydrated formic acid (FA) clusters upon collision with electrons of variable energy, focusing on electron ionization at 70 eV (EI) and low-energy (1.5-15 eV) electron attachment (EA). To uncover details about the composition of neutral clusters, we aim to elucidate the ion formation processes in FAM ⋅ WN clusters initiated by interaction with electrons and determine the extent of cluster fragmentation. EI predominantly produces protonated [FAm+H]+ ions, and in FA-rich clusters, the stable ring structures surrounding H3O+ ions are formed. In contrast, EA leads to a competition between the formation of intact [FAm ⋅ Wn]- and dissociated [FAm ⋅ Wn-H]- fragment ions, influenced by the cluster size, level of hydration, and electron energy. Our findings reveal a predisposition of low-energy EA towards forming [FAm ⋅ Wn]-, while higher electron energies tend to favor the formation of [FAm ⋅ Wn-H]- due to intracluster ion-molecule reactions. The comparison of positive and negative ion spectra suggests that the mass spectra of FA-rich clusters may indicate their actual size and composition. On the other hand, the more weakly bound water evaporation from the clusters depends strongly on the ionization. Thus, for the hydrated clusters, the neutral cluster size can hardly be estimated from the mass spectra.

Ultrafast Proton Transfer and Contact Ion-Pair Formation in Formic Acid Clusters

The ultrafast proton transfer dynamics of homogeneous formic acid clusters (FA)n, n < 10, are investigated with femtosecond time-resolved mass spectrometry. We monitor the proton transfer pathway following Rydberg state electronic relaxation and find that successful ion pair formation increases logarithmically with cluster size. Ab initio calculations demonstrate similar excitation/relaxation behavior for each cluster, revealing a contact ion pair forms between two molecules composing the cluster before finally a formate anion (HCOO-) is dissociated by the probe pulse. The sub-ps time scale for rearrangement and proton transfer increases almost linearly with cluster size, requiring ∼67 fs per additional formic acid molecule and ranging from 213 ± 51 fs for the trimer to 667 ± 116 fs for FA9. The near-linear trends measured for both rearrangement lifetime and ion pair formation suggest that proton transfer is unlikely in the formic acid dimer but becomes prominent in small clusters.

An Ab Initio Neural Network Potential Energy Surface for the Dimer of Formic Acid and Further Quantum Tunneling Dynamics

We construct a full-dimensional ab initio neural network potential energy surface (PES) for the isomerization system of the formic acid dimer (FAD). This is based upon ab initio calculations using the DLPNO-CCSD(T) approach with the aug-cc-pVTZ basis set, performed at over 14000 symmetry-unique geometries. An accurate fit to the obtained energies is generated using a general neural network fitting procedure combined with the fundamental invariant method, and the overall energy-weighted root-mean-square fitting error is about 6.4 cm^-1. Using this PES, we present a multidimensional quantum dynamics study on tunneling splittings with an efficient theoretical scheme developed by our group. The ground-state tunneling splitting of FAD calculated with a four-mode coupled method is in good agreement with the most recent experimental measurements. The PES can be applied for further dynamics studies. The effectiveness of the present scheme for constructing a high-dimensional PES is demonstrated, and this scheme is expected to be feasible for larger molecular systems.

Wavepacket dynamical study of H-atom tunneling in catecholate monoanion: the role of intermode couplings and energy flow

We present a study of H-atom tunneling in catecholate monoanion through wavepacket dynamical simulations. In our earlier study of this symmetrical double-well system [Phys. Chem. Chem. Phys., 2022, 24, 10887], a limited number of transition state modes were identified as being important for the tunneling process. These include the imaginary frequency mode Q₁, the CO scissor mode Q₁₀, and the OHO bending mode Q₂₉. In this work, starting from non-stationary initial states prepared with excitations in these modes, we have carried out wavepacket dynamics in two and three dimensional spaces. We analyse the dynamical effects of the intermode couplings, in particular the role of energy flow between the studied modes on H-atom tunneling. We find that while Q₁₀ strongly modulates the donor-acceptor distance, it does not exchange energy with Q₁. However, excitation in Q₂₉ or Q₁ does lead to rapid energy exchange between these modes, which modifies the tunneling rate at early times.

Multidimensional H-atom tunneling in the catecholate monoanion

We present the catecholate monoanion as a new model system for the study of multidimensional tunneling. It has a symmetrical O-H double-well structure, and the H atom motion between the two wells is coupled to both low and high frequency modes with different strengths. With a view to studying mode-specific tunneling in the catecholate monoanion, we have developed a full (33) dimensional potential energy surface in transition state (TS) normal modes using a Distributed Gaussian Empirical Valence Bond (DGEVB) based approach. We have computed eigenstates in different subspaces using both unrelaxed and relaxed potentials based on the DGEVB model. With unrelaxed potentials, we present results up to 7D subspaces that include the imaginary frequency mode and six modes coupled to it. With relaxed potentials, we focus on the two most important coupling modes. The structures of the ground and vibrationally excited eigenstates are discussed for both approaches and mode-specific tunneling splitting and their trends are presented.

Vibrational Stark fields in carboxylic acid dimers

Carboxylic acids form exceptionally stable dimers and have been used to model proton and double proton transfer processes. The stabilization energies of the carboxylic acid dimers are very weakly dependent on the nature of substitution. However, the electric field experienced by the OH group of a particular carboxylic acid is dependent more on the nature of the substitution on the dimer partner. In general, the electric field was higher when the partner was substituted with an electron-donating group and lower with an electron-withdrawing substituent on the partner. The Stark tuning rate (Δ) of the O-H stretching vibrations calculated at the MP2/aug-cc-pVDZ level was found to be weakly dependent on the nature of substitution on the carboxylic acid. The average Stark tuning rate of O-H stretching vibrations of a particular carboxylic acid when paired with other acids was 5.7 cm^-1 (MV cm^-1)^-1, while the corresponding average Stark tuning rate of the partner acids due to a particular carboxylic acid was 21.9 cm^-1 (MV cm^-1)^-1. The difference in the Stark tuning rate is attributed to the primary and secondary effects of substitution on the carboxylic acid. The average Stark tuning rate for the anharmonic O-D frequency shifts is about 40-50% higher than the corresponding harmonic O-D frequency shifts calculated at the B3LYP/aug-cc-pVDZ level, much greater than the typical scaling factors used, indicating the strong anharmonicity of O-H/O-D oscillators in carboxylic acid dimers. Finally, the linear correlation observed between pK_a and the electric field was used to estimate the pK_a of fluoroformic acid to be around 0.9.

CC-stretched formic acid: isomerisation, dimerisation, and carboxylic acid complexation

The cis-trans-isomerism of the propiolic acid monomer (HC[triple bond, length as m-dash]C-COOH) is examined with linear Raman jet spectroscopy, yielding the first environment-free vibrational band centres of a higher-energy cis-rotamer beyond formic acid (HCOOH) in addition to all fundamentals and a large number of hot and combination/overtone bands of the trans-conformer. Two near-isoenergetic trans-fundamentals of different symmetry (CC[double bond, length as m-dash]O bend and OH torsion) prove to be a sensitive benchmarking target, as their energetic order is susceptible to the choice of electronic structure method, basis set size, and inclusion of vibrational anharmonicity. For the infrared- and Raman-active C[double bond, length as m-dash]O stretching fundamentals of the cyclic (C_2h) trans-propiolic acid dimer, resonance couplings are found that in part extend to the C_s-symmetric heterodimer of trans-propiolic and trans-formic acid. Exploratory vibrational perturbation theory (VPT2) calculations show that all perturbing states involve displacements of the OH moieties located on the doubly hydrogen bonded ring. The comparison of the infrared spectra of the propiolic acid dimer and its heterodimer with formic acid to that of several other carboxylic acid dimers from the literature reveals a notable similarity regarding a non-fundamental dimer band around 1800 cm^-1, which in most cases is so far unassigned. VPT2 calculations and a simple harmonic model suggest an assignment to a combination vibration of the symmetric and antisymmetric OH torsion.

Fingerprint region of the formic acid dimer: variational vibrational computations in curvilinear coordinates

Curvilinear kinetic energy models are developed for variational nuclear motion computations including the inter- and the low-frequency intra-molecular degrees of freedom of the formic acid dimer. The coupling of the inter- and intra-molecular modes is studied by solving the vibrational Schrödinger equation for a series of vibrational models, from two up to ten active vibrational degrees of freedom by selecting various combinations of active modes and constrained coordinate values. Vibrational states, nodal assignment, and infrared vibrational intensity information is computed using the full-dimensional potential energy surface (PES) and electric dipole moment surface developed by Qu and Bowman [Phys. Chem. Chem. Phys., 2016, 18, 24835; J. Chem. Phys., 2018, 148, 241713]. Good results are obtained for several fundamental and combination bands in comparison with jet-cooled vibrational spectroscopy experiments, but the description of the ν₈ and ν₉ fundamental vibrations, which are close in energy and have the same symmetry, appears to be problematic. For further progress in comparison with experiment, the potential energy surface, and in particular, its multi-dimensional couplings representation, requires further improvement.

Efficient Quantum Mechanical Calculations of Mode-Specific Tunneling Splittings upon Fundamental Excitation in the Dimer of Formic Acid

The formic acid dimer (FAD) is an important benchmark system for understanding the double hydrogen transfer process. Most recently, Zhang et al. measured a few tunneling splittings upon fundamental excitation of FAD precisely (Zhang, Y. et al. J. Chem. Phys. 2017, 146, 244306); however, relevant theoretical studies are very limited. Here, we present a multidimensional quantum dynamics study on mode-specific tunneling splittings upon fundamental excitation in FAD with an efficient theoretical scheme developed by our group in which the process-oriented basis function customization strategy is combined with the preconditioned inexact spectral transform method. Various mode-specific tunneling splittings upon fundamental excitation are systematically calculated, and interesting mode-specific excitation effects on tunneling rate are identified. In particular, the calculated tunneling splittings for the ν₂₂ and ν₂₁ states are in good agreement with experiment, and the remarkable mode-specific suppression effects upon excitation should result from that the antisymmetric vibrational modes hinder the concerted double H-transfer. The present work is helpful to acquire a better understanding of the mode-specific excitation effects on tunneling processes.

Concerted Pair Motion Due to Double Hydrogen Bonding: The Formic Acid Dimer Case

Formic acid dimer as the prototypical doubly hydrogen-bonded gas-phase species is discussed from the perspective of the three translational and the three rotational degrees of freedom which are lost when two formic acid molecules form a stable complex. The experimental characterisation of these strongly hindered translations and rotations is reviewed, as are attempts to describe the associated fundamental vibrations, their combinations, and their thermal shifts by different electronic structure calculations and vibrational models. A remarkable match is confirmed for the combination of a CCSD(T)-level harmonic treatment and an MP2-level anharmonic VPT2 correction. Qualitatively correct thermal shifts of the vibrational spectra can be obtained from classical molecular dynamics in CCSD(T)-quality force fields. A detailed analysis suggests that this agreement between experiment and composite theoretical treatment is not strongly affected by fortuitous error cancellation but fully converged variational treatments of the six pair or intermolecular modes and their overtones and combinations in this model system would be welcome.

Double Proton Transfer in the Dimer of Formic Acid: An Efficient Quantum Mechanical Scheme

Double proton transfer plays an important role in biology and chemistry, such as with DNA base pairs, proteins and molecular clusters, and direct information about these processes can be obtained from tunneling splittings. Carboxylic acid dimers are prototypes for multiple proton transfer, of which the formic acid dimer is the simplest one. Here, we present efficient quantum dynamics calculations of ground-state and fundamental excitation tunneling splittings in the formic acid dimer and its deuterium isotopologues. These are achieved with a multidimensional scheme developed by us, in which the saddle-point normal coordinates are chosen, the basis functions are customized for the proton transfer process, and the preconditioned inexact spectral transform method is used to solve the resultant eigenvalue problem. Our computational results are in excellent agreement with the most recent experiments (Zhang et al., 2017; Li et al., 2019).

S–H rotamerizationviatunneling in a thiol form of thioacetamide

Rotamerization of the S–H groupviahydrogen tunneling is reported for the first time.

The collective and quantum nature of proton transfer in the cyclic water tetramer on NaCl(001)

Proton tunneling is an elementary process in the dynamics of hydrogen-bonded systems. Collective tunneling is known to exist for a long time. Atomistic investigations of this mechanism in realistic systems, however, are scarce. Using a combination of ab initio theoretical and high-resolution experimental methods, we investigate the role played by the protons on the chirality switching of a water tetramer on NaCl(001). Our scanning tunneling spectroscopies show that partial deuteration of the H₂O tetramer with only one D₂O leads to a significant suppression of the chirality switching rate at a cryogenic temperature (T), indicating that the chirality switches by tunneling in a concerted manner. Theoretical simulations, in the meantime, support this picture by presenting a much smaller free-energy barrier for the translational collective proton tunneling mode than other chirality switching modes at low T. During this analysis, the virial energy provides a reasonable estimator for the description of the nuclear quantum effects when a traditional thermodynamic integration method cannot be used, which could be employed in future studies of similar problems. Given the high-dimensional nature of realistic systems and the topology of the hydrogen-bonded network, collective proton tunneling may exist more ubiquitously than expected. Systems of this kind can serve as ideal platforms for studies of this mechanism, easily accessible to high-resolution experimental measurements.

High-Resolution Absorption and Electronic Circular Dichroism Spectra of (R)-(+)-1-Phenylethanol. Confident Interpretation Based on the Synergy between Experiments and Computations

Using density functional theory and its time-dependent extension for excited states, the S₀ →S₁ high-resolution vibronic absorption and electronic circular dichroism spectra of (R)-(+)-1-phenylethanol are computed and compared to experimental spectra measured in jet-cooled conditions in the region within 1000 cm^-1 of the 0-0 transition. The agreement between theory and computation is satisfactory and allows a confident assignment of several experimental bands in terms of fundamentals of different modes. Cases are documented for which the analysis of optical anisotropy factors, owing to their signed nature, remarkably enhances the possibility of a robust assignment of the experimental absorption bands. Computational analysis shows that the experimental spectra are dominated by Herzberg-Teller contributions and that the electronic circular dichroism spectrum and the anisotropy factors are also strongly modulated by the effect of Duschinsky mixings.

Collective proton transfer in ordinary ice: local environments, temperature dependence and deuteration effects

Ordinary ice at low temperature: what about collective nuclear quantum effects in its chiral six rings?

Full- and reduced-dimensionality instanton calculations of the tunnelling splitting in the formic acid dimer

Nearly all degrees of freedom need to be included for accurate theoretical predictions of quantum dynamics.

Predicting collision-induced dissociation mass spectra: understanding the role of the mobile proton in small molecule fragmentation

RATIONALE

Intramolecular proton migration has been reported to be required for fragmentation by collision-induced dissociation (CID). If the collision energy is required to provide energy for proton movement to a ‘dissociative’ site, it may be possible to predict the optimal collision energy for fragmentation using quantum computational chemistry software. A greater understanding of the mechanism(s) of proton migration is necessary.

METHODS

The product ion spectra of seven compounds were obtained at collision energies stepped in the range from 5 to 50 eV, with precursor ions being generated in positive ion mode by both atmospheric pressure chemical ionization (APCI) and electrospray ionisation (ESI) (using an ESCi ionisation source with or without corona discharge, respectively). The products ions observed at each collision energy were assessed in terms of structure to ascertain if they were formed as a result of protonation at the initial ionisation site or if the proton had migrated to a dissociative site.

RESULTS

Proton migration was shown to be independent of collision energy, stability of the protonated molecule and the distance that the proton moved. Therefore, proton migration is not a barrier to fragmentation as the proton appears to be fully mobile at 5 eV. As proton migration is independent of collision energy for these compounds, whereas fragmentation is energy dependent, protonation at the dissociative site alone is not sufficient to cause bond cleavage.

CONCLUSIONS

The role of collision energy in bond cleavage may be to increase the vibrational energy of the bond and/or increase the rate of bond cleavage such that it occurs within the residence time of the ion within the collision cell rather than to supply the energy for proton migration. Therefore, quantum chemistry alone cannot predict the collision energies appropriate for fragmentation on the basis of modelling proton movements.

An ab initio potential energy surface for the formic acid dimer: zero-point energy, selected anharmonic fundamental energies, and ground-state tunneling splitting calculated in relaxed 1–4-mode subspaces

We report a full-dimensional, permutationally invariant potential energy surface (PES) for the cyclic formic acid dimer.

Exceptional isotopic-substitution effect: breakdown of collective proton tunneling in hexagonal ice due to partial deuteration

Multiple proton transfer controls many chemical reactions in hydrogen-bonded networks. However, in contrast to well-understood single proton transfer, the mechanisms of correlated proton transfer and of correlated proton tunneling in particular have remained largely elusive. Herein, fully quantized ab initio simulations are used to investigate H/D isotopic-substitution effects on the mechanism of the collective tunneling of six protons within proton-ordered cyclic water hexamers that are contained in proton-disordered ice, a prototypical hydrogen-bonded network. At the transition state, isotopic substitution leads to a Zundel-like complex, [HO⋅⋅⋅D⋅⋅⋅OH], which localizes ionic defects and thus inhibits perfectly correlated proton tunneling. These insights into fundamental aspects of collective proton tunneling not only rationalize recent neutron-scattering experiments, but also stimulate investigations into multiple proton transfer in hydrogen-bonded networks much beyond ice.

Understanding collision-induced dissociation of dofetilide: a case study in the application of density functional theory as an aid to mass spectral interpretation

Fragmentation of molecules under collision-induced dissociation (CID) conditions is not well-understood. This may make interpretation of MSMS spectra difficult and limit the effectiveness of software tools intended to aid mass spectral interpretation. Density Functional Theory (DFT) has been successfully applied to explain the thermodynamics of fragmentation in the gas phase by the modelling the effect that protonation has on the bond lengths (and hence bond strengths). In this study, dofetilide and four methylated analogues were used to investigate further the potential for using DFT to understand and predict the CID fragmentation routes. The products ions present in the CID spectra of all five compounds were consistent with charge-directed fragmentation, with protonation adjacent to the cleavage site being required to initiate fragmentation. Protonation at the dissociative site may have occurred either directly or via proton migration. A correlation was observed between protonation-induced bond lengthening and the bonds which were observed to break in the CID spectra. This correlation was quantitative in that the bonds calculated to elongate to the greatest extent gave rise to the most abundant of the major product ions. Thus such quantum calculations may offer the potential for contributing to a predictive tool for aiding the accuracy and speed mass spectral interpretation by generating numerical data in the form of bond length increases to act as descriptors flagging potential bond cleavages.

Molecular structure, spectroscopic (FT-IR, FT-Raman), NBO and HOMO-LUMO analyses, computation of thermodynamic functions for various temperatures of 2, 6-dichloro-3-nitrobenzoic acid

The FT-IR and Raman spectra of 2, 6-dichloro-3-nitrobenzoic acid (DCNBA) have been recorded and analyzed. The equilibrium geometry, various bonding and harmonic vibrational wavenumbers have been calculated with the help of density functional theory (DFT/B3LYP/cc-pvdz/6-311++G(d,p)) method. The optimized geometrical parameters obtained by B3LYP/6-311++G(d,p) method show good agreement with experimental X-ray data. Most of the vibrational modes are observed in the expected range. The Mulliken population analysis shows the interactions of C-N-O···H-C and C-O···H-C. The most possible interaction is explained using natural bond orbital (NBO) analysis. The effects of molecular association through O-H···O hydrogen bonding have been described by the single dimer structure. The strengthening and polarization of the C=O bond increases due to the degree of conjugation. HOMO-LUMO energy and the thermodynamic parameters are also evaluated. The thermodynamic functions (heat capacity, internal heat energy, Gibbs energy and entropy) from spectroscopic data by statistical methods were obtained for the range of temperature 100-1000 K.

Proton Tunneling in Heterodimers of Carboxylic Acids: A Rotational Study of the Benzoic Acid-Formic Acid Bimolecule

Tunneling effects have been measured in the pulsed jet Fourier transform microwave spectra of two isotopologues of the benzoic acid-formic acid bimolecule. The tunneling splittings are originated by the concerted proton transfer of the two carboxylic hydrogens. From the values of these splittings for the OH-OH and OD-OD species, it has been possible to model/size the barrier to the concerted double proton transfer.

Hydrated alizarin complexes: hydrogen bonding and proton transfer

We investigated the hydrogen bonding structures and proton transfer for the hydration complexes of alizarin (Az) produced in a supersonic jet using fluorescence excitation (FE), dispersed laser induced fluorescence (LIF), visible-visible hole burning (HB), and fluorescence detected infrared (FDIR) spectroscopy. The FDIR spectrum of bare Az with two O-H groups exhibits two vibrational bands at 3092 and 3579 cm(-1), which, respectively, correspond to the stretching vibration of O1-H1 that forms a strong intramolecular hydrogen bond with the C9=O9 carbonyl group and the stretching vibration of O2-H2 that is weakly hydrogen-bonded to O1-H1. For the 1:1 hydration complex Az(H(2)O)(1), we identified three conformers. In the most stable conformer, the water molecule forms hydrogen bonds with the O1-H1 and O2-H2 groups of Az as a proton donor and proton acceptor, respectively. In the other conformers, the water binds to the C10=O10 group in two nearly isoenergetic configurations. In contrast to the sharp vibronic peaks in the FE spectra of Az and Az(H(2)O)(1), only broad, structureless absorption was observed for Az(H(2)O)(n) (n≥ 2), indicating a facile decay process, possibly due to proton transfer in the electronic excited state. The FDIR spectrum with the wavelength of the probe laser fixed at the broad band exhibited a broad vibrational band near the O2-H2 stretching vibration frequency of the most stable conformer of Az(H(2)O)(1). With the help of theoretical calculations, we suggest that the broad vibrational band may represent the occurrence of proton transfer by tunnelling in the electronic ground state of Az(H(2)O)(n) (n≥ 2) upon excitation of the O2-H2 vibration.