Infrared Spectroscopy of the Astrochemically Relevant Protonated

Dec 13, 2017 - Special Issue. Published as part of The Journal of Physical Chemistry virtual special issue “W. Lester S. Andrews Festschrift”...
0 downloads 0 Views 899KB Size
Article pubs.acs.org/JPCA

Cite This: J. Phys. Chem. A 2018, 122, 192−198

Infrared Spectroscopy of the Astrochemically Relevant Protonated Formaldehyde Dimer Published as part of The Journal of Physical Chemistry virtual special issue “W. Lester S. Andrews Festschrift”. J. Philipp Wagner, David C. McDonald, II, and Michael A. Duncan* Department of Chemistry, University of Georgia, 140 Cedar Street, Athens, Georgia 30602, United States

Downloaded via UNIV OF SUSSEX on June 27, 2018 at 21:29:40 (UTC). See https://pubs.acs.org/sharingguidelines for options on how to legitimately share published articles.

S Supporting Information *

ABSTRACT: The protonated formaldehyde dimer (H2CO)2H+ was generated in an electrical discharge and supersonic expansion of argon saturated with formalin solution vapor. Its infrared spectrum was measured in the region from 900 to 4000 cm−1 employing infrared laser photodissociation and messenger atom tagging. Comparison of the experiment to quantum chemical computations at the CCSD(T)/cc-pVQZ//MP2/cc-pVTZ level reveals that the experimentally observed structure is the head-to-tail dimer and not the more stable proton-bound dimer. This is consistent with the usually observed C−O bond formation upon formaldehyde oligomerization under acidic conditions in solution and resembles the structure of the neutral (H2CO)2 dimer in the gas phase. There is no evidence for the formation of other isomers, most notably protonated glycolaldehyde, that could result from covalent bond formation. These findings may be relevant to a proposed carbohydrate formation mechanism in the interstellar medium starting from protonated formaldehyde dimer. stoichiometry C2H5O2+ in the gas phase starting from formaldehyde vapor and probe its structure with infrared laser photodissociation spectroscopy. In solution, formaldehyde undergoes two kinds of oligomerization reactions. While C−O bond formation is preferred under acidic conditions, yielding linear polyoxymethylenes or cyclic oxanes, basic conditions favor reaction over the carbon atoms producing glycolaldehyde and monosaccharides.21 The latter reaction was described by Butlerow as early as 1861 and is generally referred to as the formose reaction.22 This reaction is now generally regarded to be important for the formation of prebiotic molecules.23,24 Breslow and others have rationalized its mechanism as a sequence of aldol reactions.25,26 However, the first step of formaldehyde dimerization cannot be explained this way, and the participation of “activated formaldehyde” has been invoked.26 Most interestingly, glycolaldehyde was detected toward Sagittarius B2(N) by millimeter-wave spectroscopy27,28 causing speculation about a gas-phase version of the formose reaction in the interstellar medium. Given the abundance of protonating agents like H3+ in space,29 a Nazarov-type cyclization of proton-bound formaldehyde dimer with subsequent [1,2]H-shift (Scheme 1) has been proposed to be a source of protonated glycolaldehyde.9,10 Although this

1. INTRODUCTION Formaldehyde (H2CO) is a well-known, abundant polyatomic organic in the interstellar medium (ISM).1 Although it is relatively unreactive in the gas phase, formaldehyde’s interaction with its interstellar protonated analogue (H2COH+)2 provides access to a rich chemistry due to the general importance of barrierless ion−molecule reactions in space.3−8 It has been proposed that a reaction beginning with the protonated dimer (Ia) could potentially lead to protonated glycolaldehyde (IIa) in the ISM, which formally corresponds to a two-carbon aldose (Scheme 1).9,10 Otherwise, in mass spectrometry and computational chemistry the structure of the ion is usually attributed to the proton-bound dimer (Ia).11−14 With new techniques in ion infrared spectroscopy, it has become possible to determine the structures of more complex organic ions in the gas phase.15−20 Therefore, to investigate this structural issue, we set out to prepare the ion of the Scheme 1. Proposed Gas-Phase Dimerization Reactiona of Proton-Bound Formaldehyde Dimer with Subsequent Rearrangement to Protonated Glycolaldehyde

Received: October 25, 2017 Revised: December 9, 2017 Published: December 13, 2017

a

Drawn in analogy to Figure 5 in ref 9. © 2017 American Chemical Society

192

DOI: 10.1021/acs.jpca.7b10573 J. Phys. Chem. A 2018, 122, 192−198

Article

The Journal of Physical Chemistry A

fragment ion yield from argon loss as a function of the energy of the laser’s photon energy. Computational. Conceivable noncovalently bound isomers of protonated formaldehyde dimer were explored and optimized at the MP2/cc-pVTZ level of theory utilizing the Gaussian09 program package.55,56 Additionally, geometries of other molecular ions of the same stoichiometry were localized at the same level.57 We chose this level of theory, because a recent theoretical study on the neutral formaldehyde dimer revealed that MP2 binding energies differ by no more than ∼0.2 kcal mol−1 from those obtained with explicitly correlated CCSD(T)-F12 theory.58 MP2 also seems a reasonable choice, since it has been found that it performs better for hydrogenbonded systems than the more modern SCS-MP2,59 and DFT often has difficulties when it comes to processes related to proton accommodation and transfer.60 Although density functional theory (DFT) methods with dispersion-corrected functionals are often employed for improved energetics of noncovalent systems, we found that such methods are no better for the vibrational patterns of small ions.47 However, to explore this issue further, we also conducted DFT computations at the B3LYP-D3(BJ)/cc-pVTZ level on these various structures, as shown in the Supporting Information. There was no significant difference in the vibrational spectra predicted by these different methods. For each of the methods employed, harmonic vibrational frequencies were computed to ascertain the nature of the stationary points on the potential energy surfaces (PES). For comparison to experimental infrared spectra, an empirical scaling factor of 0.961 was applied to the MP2 frequencies, derived from a comparison of computed frequencies for neutral formaldehyde to its known vibrations.61 Zero-point vibrational energies (ZPVE) are included in all the energies reported here (ΔH0). For further refinement of energies, CCSD(T)/cc-pVQZ single points were computed on top of the optimized structures.62,63

idea is tempting, the proposed reaction seems unlikely based on simple polarity arguments. On the one hand, a flow drift tube experiment in which H3O+ was reacted with formaldehyde gas yielded an ion with the m/z ratio of protonated glycolaldehyde, but no structural evidence was provided.10 On the other hand, it was demonstrated that in its unimolecular decomposition protonated glycolaldehyde exclusively loses water and not formaldehyde.30 In a similar vein, protonated glyceraldehyde and dihydroxyacetone (C3 aldose and ketose, respectively) lost formaldehyde in unimolecular decomposition, but the back reactions were computed to have nonnegligible barriers.31 A more likely pathway to interstellar glycolaldehyde might be grain-surface hydrogenation of carbon monoxide in combination with C−C bond formation by radical recombination,32 similar to the formation of interstellar formaldehyde itself.33,34 The starting point for the proposed reaction sequence (Scheme 1) is the proton-bound dimer of formaldehyde, in which a proton links the two oxygens in an O−H+−O structural motif. Such proton-bound dimers have been studied for many years with gas-phase infrared spectroscopy to develop a microscopic understanding of proton accommodation and proton transport.35−49 Infrared bands associated with the shared proton (νsp) can occur in a broad spectral range, and their frequencies have been found to correlate with the difference in proton affinity of the two involved bases. In symmetrical systems, νsp is found near or below 1000 cm−1 and is often split into multiplets due to coupling with other vibrational modes.40 In this regard, the aggregation of aldehydes upon protonation was investigated with mass spectrometry. Whereas cluster ions with up to 11 monomer units have been observed for acetaldehyde,50 in the formaldehyde case only the dimer was found.11 Quantum chemical studies of the (H2CO)2H+ complex assigned the lowest energy structure to be the proton-bound dimer displaying a doublewell minimum at low levels of theory.12,13 When correlation is included into the computations, the proton transfer barrier collapses to a single minimum with equal proton sharing.10 Spectroscopically, only the protonated formaldehyde monomer has been studied so far.51−53 Our group investigated the chemically related protonated acetone dimer with infrared photodissociation spectroscopy and assigned its structure to a proton-bound dimer.44 The proton stretch in that symmetric system was found at 900 cm−1, with a series of intense combination bands extending to higher frequencies just above that.44 We therefore aim to further elucidate the gas-phase chemistry of the protonated formaldehyde dimer with a combined infrared and computational study.

3. RESULTS AND DISCUSSION The experimentally recorded photodissociation spectrum of the (H2CO)2H+Ar ion is presented in the upper black trace in Figure 1. The spectrum is dominated by a strong band at 3401 cm−1 in the O−H stretching region accompanied by two weaker features at 3344 and 3560 cm−1. A band at 3145 cm−1 indicates a C−H stretch at an unsaturated carbon center. Three bands at 1618, 1683, and 1710 cm−1 suggest the presence of one or more carbonyl groups. In the fingerprint region, various transitions of medium (1489, 1248, and 983 cm−1) to weak (1399 and 1162 cm−1) intensity are found. There are no strong infrared transitions near or below 1000 cm−1 like that expected for a proton-bound dimer, although such a band might still be found outside the experimental range. On the contrary, the dominant absorption at 3401 cm−1 indicates that there is a free or weakly coordinated hydroxy group. To further understand this spectrum, we explored possible isomers of protonated formaldehyde dimer in which the monomeric subunits are still intact. The optimized structures are displayed as complexes Ia−d in Figure 2. Recognizing the possibility of reactions in the plasma, molecular ions isobaric to protonated formaldehyde dimer were optimized as well and are shown as structures IIa−e (Figure 2). For the noncovalent ions it turns out that the C2h symmetric proton-bound dimer Ia is the most stable isomer in agreement with a previous theoretical study.13 We compute a D0 binding energy of

2. METHODS Experimental Section. Ions of the composition C2H5O2+ were generated by a pulsed electrical discharge in a supersonic expansion of 5% hydrogen in argon saturated with formalin solution vapor at ambient temperature, as used in a recent study of the formaldehyde radical cation.54 The argon-tagged complex of this ion was mass-selected in a reflectron time-offlight spectrometer, and its infrared spectrum was recorded with laser photodissociation spectroscopy.20 Infrared transitions were explored in the range of 900−4000 cm−1 with a Nd:YAG-pumped infrared OPO/OPA laser system (LaserVision). For the low-energy region (900−2050 cm−1), a AgGaSe2 crystal was used in a second stage of difference frequency generation. The action spectra were collected as 193

DOI: 10.1021/acs.jpca.7b10573 J. Phys. Chem. A 2018, 122, 192−198

Article

The Journal of Physical Chemistry A

providing added stability to this configuration. The C−O bond length of the protonated aldehyde unit is substantially elongated to 1.33 Å, indicating that the double-bond character of the attacked carbonyl group is greatly reduced. This is also expected to lower the C−O vibrational frequency significantly compared to a normal carbonyl (e.g., 1628 cm−1 in the protonated formaldehyde monomer53). Dimer Ic would not lead to a sugar-like structure upon reaction but rather to an oxane (IIb) due to the preformed connectivity. Finally, a less important CH−O bound isomer Id was calculated at an energy of 17.7 kcal mol−1. The two most important covalently bound isomers result from C−C and C−O bond formation, respectively: protonated glycolaldehyde IIa and protonated oxane IIb (formally 1,3dioxethane). While formation of protonated glycolaldehyde is favorable by 7.1 kcal mol−1 relative to Ia, the reaction to IIb is highly endothermic (+29.2 kcal mol−1), probably due to a high strain in the four-membered ring.65 Other chemically meaningful, isobaric ions are less likely to form and would require substantial molecular rearrangement. The most stable of these is protonated acetic acid IIc that is 35.7 kcal mol−1 lower in energy than Ia. The second most-stable ion is protonated methyl formate (IId, −17.3 kcal mol−1). Both acetic acid and methyl formate are known molecules in the ISM.66,67 The least energetically favorable ion is the peroxide IIe (+82.7 kcal mol−1), which would formally derive from a [2 + 2] cycloaddition. Figure 1 shows a comparison of our experimental spectrum to the scaled harmonic infrared spectra predicted for these various isomeric structures. This comparison suggests that the ion detected in our experiment is not the expected protonbound formaldehyde dimer Ia. Dimer Ia lacks an appropriate high-energy mode to explain the strongest band at 3401 cm−1 (green trace, Figure 1). The computed spectral pattern of protonated glycolaldehyde also does not match the observed spectrum (purple trace, Figure 1). Moreover, ion IIa displays an intramolecular ionic hydrogen bond, which often appears as a broadened band in the spectrum;68,69 we do not observe such a feature. The other covalent isomers IIb−e can also be excluded on the basis of similar spectral comparisons. Figure S14 in the Supporting Information shows the computed infrared spectra of all the covalent structures IIa−e compared to the experiment. As the formalin solution used is stabilized with methanol, protonated methanol-carbon monoxide com-

Figure 1. Comparison of the experimental infrared photodissociation spectrum of the C2H5O2+Ar ion and the predicted scaled harmonic spectra of the most relevant isomers employing MP2 theory together with a cc-pVTZ basis set.

formaldehyde to protonated formaldehyde of 31.6 kcal mol−1. In contrast to earlier lower level computations, this ion does not display a symmetric double-well potential but exhibits a single minimum with equal O−H distances (C2h).13,14 The cisrotamer Ib (C2v), which could lead to the proposed Nazarov cyclization, turns out to be a transition state 2.9 kcal mol−1 higher in energy. The head-to-tail dimer Ic featuring a C−O connectivity is only 3.9 kcal mol−1 less stable than Ia. In this structure, the nucleophilic, nonprotonated aldehyde oxygen attacks the carbonyl group carbon of the protonated monomer at an O−C−O angle of 110.8° close to the “Bürgi-Dunitz approach trajectory” known for carbonyl reactions.64 Also, the structure exhibits a nonclassical CH−O hydrogen bond

Figure 2. Structures and energies (ΔH0) of ions of the stoichiometry C2H5O2+ at CCSD(T)/cc-pVQZ//MP2/cc-pVTZ level including ZPVE. 194

DOI: 10.1021/acs.jpca.7b10573 J. Phys. Chem. A 2018, 122, 192−198

Article

The Journal of Physical Chemistry A

at 3344 cm−1 is not explained easily. One might consider the possibility of a Fermi resonance resulting from mixing with overtones of the vibrations around 1700 cm−1, but they are probably not intense enough to cause such an effect. Other isomers in which the argon binds to different sites of the molecule would be blue-shifted from the main O−H stretching feature. Also, the presence of further rotamers with respect to the intermolecular C−O bond of Ic that might cause this absorption can be considered unlikely based on a relaxed potential energy scan (Figure S18, Supporting Information). The band at 3145 cm−1 can be assigned to the asymmetric C− H stretch of the nonprotonated formaldehyde. The carbonyl stretch of the attacking formaldehyde unit is assigned to the broad band at 1618 cm−1. Unfortunately, the two accompanying bands at 1683 and 1710 cm−1 are not easy to understand. The carbonyl stretch of the protonated formaldehyde is expected to be at a much lower frequency because of its reduced double bonding, as noted earlier. The modes below 1500 cm−1 are difficult to describe, as their character is often a mixture of H−C−H scissoring, wagging and rocking, H−O−C deformation, and C−O bond stretching (of the protonated formaldehyde unit). All these modes are similar in energy, which can give rise to complex coupling and calls for a better theoretical treatment in the future. Even though there are unassigned bands, the main features of the experimental infrared spectrum are understandable. We therefore conclude that the C2H5O2+ in our experiment is unlikely to be the proton-bound dimer of formaldehyde or the protonated glycolaldehyde ion but instead corresponds to the previously unanticipated head-to-tail dimer Ic. This is surprising, as this complex is computed to be less stable than its competing structural isomers. Therefore, one might speculate that this structure forms for kinetic reasons in the ionized and rapidly cooled supersonic expansion. Because of its large excess in the expansion, it is conceivable that argon attaches to the protonated formaldehyde monomer first, blocking the O−H group for further interactions. Because of the low temperature, the second formaldehyde would then attach to the carbonyl group in an antiparallel fashion. If this is the case, tagging with a more weakly bound atom such as neon may make it possible for the second formaldehyde molecule to displace the rare gas and find its more stable binding site. We tried to make such a neon-tagged species but were unfortunately not able to produce enough of this ion for study. Another possibility is that the neutral formaldehyde dimer (H2CO)2 forms first in the expansion and then is subsequently protonated. The structure seen here is essentially the same as that of the neutral dimer according to a microwave study,70 which is in agreement with computations at high levels of theory.71 The efficiency of formation for the neutral dimer should vary with the formaldehyde concentration and beam conditions, also perhaps providing some control over the ion cluster growth. Again, our attempts to produce other spectra by varying these conditions were unsuccessful. The most stable proton-bound structure would likely be produced by a reaction of the protonated monomer with a neutral monomer. A bias against this chemistry is conceivable if the monomer concentration in the beam is low compared to larger formaldehyde aggregates. Polymers of formaldehyde are of course well-known, but these species would have extremely low vapor pressures compared to monomers under our roomtemperature conditions. Additionally, we are confident that formaldehyde monomers are present, as they are intense in the

plexes could conceivably form and contribute to the experimental IR spectrum. However, the spectral comparison in Figure S15 in the Supporting Information also makes it possible to exclude this heterogeneous complex. The only ion whose predicted infrared spectrum is in convincing agreement with the observed vibrational transitions is the head-to-tail dimer Ic, whose argon-tagged isomer we also present in Figure 1 (red and blue traces, respectively). Band assignments for the argon-tagged ion based on scaled harmonic frequencies are given in Table 1. Table 1. Comparison of Observed and Scaled (0.961) Harmonic Infrared Transitions (in cm−1) at the MP2/ccpVTZ Level Together with the Approximate Description of the Mode experimental

ωa (I)

3560 3401 3145

3597 (391) 3143 (17)

1710 1683 1618 1489

1639 (59) 1439 (52)

1452 1399 1248

1366 (23) 1328 (31) 1225 (97)

1162

1192 (45)

983

1026 (25)

approximate description O−H stretch + Ar−H combination O−H stretch asym. C−H stretch, nonprotonated formaldehyde

CO stretch, nonprotonated formaldehyde H−C−H scissoring, nonprotonated formaldehyde H−O−C bending/H−C−H rocking H−C−H wagging, protonated formaldehyde H−C−H rocking, nonprotonated formaldehyde H−O−C deformation/C−O stretch, protonated formaldehyde H−O−C bending/H−C−H rocking

a

Intensities in kilometers per mole.

The strongest infrared transition for isomer Ic at 3401 cm−1 is assigned to the O−H stretch of the protonated molecule. The experimentally observed frequency shows a large deviation from the scaled harmonic MP2 computation. This suggests that this mode either displays a very large anharmonicity or that the electronic structure treatment of the noncovalent interaction is still insufficient. A more pronounced nucleophilic interaction would increase the single-bond character of the attacked protonated formaldehyde converting it into an alcohol. This would be accompanied by an expected spectral blue shift of the O−H stretch to the region typically observed for these functional groups (e.g., in methanol at 3681 cm−1).61 If the interaction is less evolved, the absorption should be closer to that of a protonated carbonyl group. In the case of the argon-tagged protonated formaldehyde monomer, this absorption was seen at 3182 cm−1.53 These considerations suggest that the intermolecular C−O bond length would elongate when higher levels of theory are applied or that there is a bondstretch isomer with a greater intermolecular distance. Unfortunately, all attempts to localize such a structure failed. More sophisticated electronic structure treatments and consideration of anharmonic effects could improve the agreement between experiment and theory in future studies. The band at 3560 cm−1 probably corresponds to a combination of the O−H with the Ar−H stretching. An analogous band 159 cm−1 to the blue of the O−H stretch was observed for the argon-tagged protonated formaldehyde monomer.53 This is in excellent agreement with the interval of 158 cm−1 that we observe here. The relatively intense band 195

DOI: 10.1021/acs.jpca.7b10573 J. Phys. Chem. A 2018, 122, 192−198

Article

The Journal of Physical Chemistry A

lations would be challenging. It would also be helpful to explore different ion−molecule reaction conditions to investigate the formation mechanism for the unexpected structure Ic and the possible formation of the more stable Ia structure.

mass spectrum produced (Figure S1), and we used these same conditions previously to measure spectra for the protonated formaldehyde monomer and the formaldehyde cation.53,54 Even allowing for these potential kinetic effects, it is remarkable that there is no signal indicating formation of the stable proton-bound dimer structure. In other ion systems that we studied previously, we found several examples of coexisting isomers separated by activation barriers that precluded the equilibration to the most stable structures.20 However, this is the first system in which we find no evidence for the most stable isomer. Note that the strong shared proton stretch vibration characteristic of the proton-bound dimer structure is notoriously difficult to handle computationally.35−48 It is conceivable that this vibration occurs at a frequency significantly lower than that predicted by theory, and this would then be lower than the range covered by the experiment. If the predicted spectrum is accurate in other respects, there would not be any strong vibrations for this isomer in the higher-frequency region. It could therefore be possible that the proton-bound dimer structure is present but just not detected. An argument against this possibility is that most proton-bound dimer spectra studied in the past,35−48 especially including that for the corresponding acetone species,44 have extensive vibrational progressions extending to higher frequencies from the proton stretch fundamental. These additional strong bands correspond to the strong couplings of the proton stretch with other low-frequency modes. Such progressions would fall in the range of the present spectrum, but we also do not detect such bands. The simplest interpretation of the overall experiment is therefore that we produced primarily isomer Ic. It is instructive to consider the implication of the present work for astrochemistry. The formation of interstellar protonated glycolaldehyde as in Scheme 1 requires the existence of the proton-bound dimer form of (CH2O)2H+, which was not observed in this experiment. However, the present experiment includes the complicating factors of argon tagging and the presence of high concentrations of neutral formaldehyde, either one of which could conceivably account for the inefficient formation of the proton-bound structure. Neither of these factors would apply in the conditions of interstellar space. Therefore, it is not clear from the present results what the likelihood is for the formation of protonbound dimer structures in space. This would presumeably involve the reaction of a protonated monomer with a neutral monomer. Such reactions could be studied with standard ion− molecule methods in different kinds of mass spectrometers or flowing afterglow instruments. However, the reaction would need subsequent spectroscopy to verify the structure. Our present instrumentation does not allow such an experiment, but it may be possible in other laboratories in the future.



ASSOCIATED CONTENT

S Supporting Information *

The Supporting Information is available free of charge on the ACS Publications website at DOI: 10.1021/acs.jpca.7b10573. The complete citation for ref 55, the mass spectrum, computed geometric structures and electronic energies for different isomers, comparison of the experimental spectrum to those of other isomers, and additional comparisons to spectra predicted by DFT (PDF)



AUTHOR INFORMATION

Corresponding Author

*E-mail: [email protected]. ORCID

Michael A. Duncan: 0000-0003-4836-106X Notes

The authors declare no competing financial interest.



ACKNOWLEDGMENTS We gratefully acknowledge support for this work by the National Science Foundation (MAD Grant No. CHE1464708). J.P.W. acknowledges the Alexander von Humboldt Foundation for a Feodor Lynen Postdoctoral Fellowship.



REFERENCES

(1) Snyder, L. E.; Buhl, D.; Zuckerman, B.; Palmer, P. Microwave Detection of Interstellar Formaldehyde. Phys. Rev. Lett. 1969, 22, 679−681. (2) Ohishi, M.; Ishikawa, S.-i.; Amano, T.; Oka, H.; Irvine, W. M.; Dickens, J. E.; Ziurys, L. M.; Apponi, A. J. Detection of A New Interstellar Molecular Ion, H2COH+ (Protonated Formaldehyde). Astrophys. J. 1996, 471, L61−L64. (3) Smith, D. The Ion Chemistry of Interstellar Clouds. Chem. Rev. 1992, 92, 1473−1485. (4) Petrie, S.; Bohme, D. K. Ions in Space. Mass Spectrom. Rev. 2007, 26, 258−280. (5) Snow, T. P.; Bierbaum, V. M. Ion Chemistry in the Interstellar Medium. Annu. Rev. Anal. Chem. 2008, 1, 229−259. (6) Klemperer, W. Astronomical Chemistry. Annu. Rev. Phys. Chem. 2011, 62, 173−184. (7) Geppert, W. D.; Larsson, M. Experimental Investigations into Astrophysically Relevant Ionic Reactions. Chem. Rev. 2013, 113, 8872−8905. (8) Herbst, E. The Synthesis of Large Interstellar Molecules. Int. Rev. Phys. Chem. 2017, 36, 287−331. (9) Halfen, D. T.; Apponi, A. J.; Woolf, N.; Polt, R.; Ziurys, L. M. A Systematic Study of Glycolaldehyde in Sagittarius B2(N) at 2 and 3 mm: Criteria for Detecting Large Interstellar Molecules. Astrophys. J. 2006, 639, 237−245. (10) Jalbout, A. F.; Abrell, L.; Adamowicz, L.; Polt, R.; Apponi, A. J.; Ziurys, L. M. Sugar Synthesis from a Gas-Phase Formose Reaction. Astrobiology 2007, 7, 433−442. (11) Karpas, Z. E.; Klein, F. S. Ion−Molecule Chemistry in Formaldehyde by Ion Cyclotron Resonance Mass Spectrometry. Int. J. Mass Spectrom. Ion Phys. 1975, 16, 289−306. (12) Hagler, A. T.; Karpas, Z.; Klein, F. S. Ion Cyclotron Resonance Mass Spectrometric and Ab Initio Studies of the Structure and Mechanism of Formation of Protonated Dimers in Simple Carbonyl

4. CONCLUSION AND OUTLOOK The protonated formaldehyde dimer has been studied with infrared photodissociation spectroscopy and computational chemistry. Comparison of the spectrum for the argon-tagged ion with scaled harmonic frequency computations at the MP2/ cc-pVTZ level reveals that the structure most likely corresponds to a noncovalently bound head-to-tail dimer Ic and not to the expected proton-bound dimer Ia or its proposed reaction product protonated glycolaldehyde IIa. An anharmonic vibrational as well as a higher-level electronic structure treatment of ions Ia and Ic would be desirable to better understand the experimental spectrum, although such calcu196

DOI: 10.1021/acs.jpca.7b10573 J. Phys. Chem. A 2018, 122, 192−198

Article

The Journal of Physical Chemistry A Compounds XYCO (X or Y = Hydrogen or Fluorine). J. Am. Chem. Soc. 1979, 101, 2191−2196. (13) Aviyente, V.; Varnali, T. Protonated Aldehyde Clusters: A Semiempirical Approach. J. Mol. Struct. 1993, 299, 191−195. (14) Chu, C.-H.; Ho, J.-J. Calculated Effects of Formaldehyde Substituents on Proton Transfer in (H2CO-H-OCX2)+. J. Am. Chem. Soc. 1995, 117, 1076−1082. (15) Ebata, T.; Fujii, A.; Mikami, N. Vibrational Spectroscopy of Small-Sized Hydrogen-Bonded Clusters and Their Ions. Int. Rev. Phys. Chem. 1998, 17, 331−361. (16) Bieske, E. J.; Dopfer, O. High-Resolution Spectroscopy of Cluster Ions. Chem. Rev. 2000, 100, 3963−3998. (17) Duncan, M. A. Frontiers in the Spectroscopy of Mass-Selected Molecular Ions. Int. J. Mass Spectrom. 2000, 200, 545−569. (18) Fridgen, T. D. Infrared Consequence Spectroscopy of Gaseous Protonated and Metal Ion Cationized Complexes. Mass Spectrom. Rev. 2009, 28, 586−607. (19) Baer, T.; Dunbar, R. C. Ion Spectroscopy: Where Did It Come From; Where Is It Now; and Where Is It Going? J. Am. Soc. Mass Spectrom. 2010, 21, 681−693. (20) Duncan, M. A. Infrared Laser Spectroscopy of Mass-Selected Carbocations. J. Phys. Chem. A 2012, 116, 11477−11491. (21) Kua, J.; Avila, J. E.; Lee, C. G.; Smith, W. D. Mapping the Kinetic and Thermodynamic Landscape of Formaldehyde Oligomerization under Neutral Conditions. J. Phys. Chem. A 2013, 117, 12658−12667. (22) Butlerow, A. Bildung einer Zuckerartigen Substanz durch Synthese. Justus Liebigs Ann. Chem. 1861, 120, 295−298. (23) Gabel, N. W.; Ponnamperuma, C. Model for Origin of Monosaccharides. Nature 1967, 216, 453−455. (24) Orgel, L. E. Prebiotic Chemistry and the Origin of the RNA World. Crit. Rev. Biochem. Mol. Biol. 2004, 39, 99−123. (25) Langenbeck, W. Ü ber Organische Katalysatoren-L: Entwicklungslinien der Organischen Katalysatoren. Tetrahedron 1958, 3, 185−196. (26) Breslow, R. On the Mechanism of the Formose Reaction. Tetrahedron Lett. 1959, 1, 22−26. (27) Hollis, J. M.; Lovas, F. J.; Jewell, P. R. Interstellar Glycolaldehyde: The First Sugar. Astrophys. J. 2000, 540, L107−L110. (28) Hollis, J. M.; Vogel, S. N.; Snyder, L. E.; Jewell, P. R.; Lovas, F. J. The Spatial Scale of Glycolaldehyde in the Galactic Center. Astrophys. J. 2001, 554, L81−L85. (29) Geballe, T. R.; Oka, T. Detection of H3+ in Interstellar Space. Nature 1996, 384, 334−335. (30) Bouchoux, G.; Penaud-Berruyer, F.; Bertrand, W. Structure, Thermochemistry and Reactivity of Protonated Glycolaldehyde. Eur. Mass Spectrom. 2001, 7, 351−357. (31) Simakov, A.; Sekiguchi, O.; Bunkan, A. J. C.; Uggerud, E. Energetics and Mechanisms for the Unimolecular Dissociation of Protonated Trioses and Relationship to Proton-mediated Formaldehyde Polymerization to Carbohydrates in Interstellar Environments. J. Am. Chem. Soc. 2011, 133, 20816−20822. (32) Fedoseev, G.; Cuppen, H. M.; Ioppolo, S.; Lamberts, T.; Linnartz, H. Experimental Evidence for Glycolaldehyde and Ethylene Glycol Formation by Surface Hydrogenation of CO Molecules under Dense Molecular Cloud Conditions. Mon. Not. R. Astron. Soc. 2015, 448, 1288−1297. (33) Watanabe, N.; Kouchi, A. Efficient Formation of Formaldehyde and Methanol by the Addition of Hydrogen Atoms to CO in H2OCO Ice at 10 K. Astrophys. J. 2002, 571, L173−L176. (34) Chuang, K. J.; Fedoseev, G.; Ioppolo, S.; van Dishoeck, E. F.; Linnartz, H. H-Atom Addition and Abstraction Reactions in Mixed CO, H2CO and CH3OH Ices - An Extended View on Complex Organic Molecule Formation. Mon. Not. R. Astron. Soc. 2016, 455, 1702−1712. (35) Asmis, K. R.; Pivonka, N. L.; Santambrogio, G.; Brümmer, M.; Kaposta, C.; Neumark, D. M.; Wöste, L. Gas-Phase Infrared Spectrum of the Protonated Water Dimer. Science 2003, 299, 1375−1377.

(36) Headrick, J. M.; Bopp, J. C.; Johnson, M. A. Predissociation Spectroscopy of the Argon-Solvated H5O2+ ″Zundel″ Cation in the 1000−1900 cm−1 Region. J. Chem. Phys. 2004, 121, 11523−11526. (37) Headrick, J. M.; Diken, E. G.; Walters, R. S.; Hammer, N. I.; Christie, R. A.; Cui, J.; Myshakin, E. M.; Duncan, M. A.; Johnson, M. A.; Jordan, K. D. Spectral Signatures of Hydrated Proton Vibrations in Water Clusters. Science 2005, 308, 1765−1769. (38) Fridgen, T. D.; MacAleese, L.; Maître, P.; McMahon, T. B.; Boissel, P.; Lemaire, J. Infrared Spectra of Homogeneous and Heterogeneous Proton-Bound Dimers in the Gas Phase. Phys. Chem. Chem. Phys. 2005, 7, 2747−2755. (39) Fridgen, T. D.; MacAleese, L.; McMahon, T. B.; Lemaire, J.; Maître, P. Gas Phase Infrared Multiple-Photon Dissociation Spectra of Methanol, Ethanol and Propanol Proton-Bound Dimers, Protonated Propanol and the Propanol/Water Proton-Bound Dimer. Phys. Chem. Chem. Phys. 2006, 8, 955−966. (40) Roscioli, J. R.; McCunn, L. R.; Johnson, M. A. Quantum Structure of the Intermolecular Proton Bond. Science 2007, 316, 249− 254. (41) Asmis, K. R.; Yang, Y.; Santambrogio, G.; Bruemmer, M.; Roscioli, J. R.; McCunn, L. R.; Johnson, M. A.; Kuehn, O. Gas-Phase Infrared Spectroscopy and Multidimensional Quantum Calculations of the Protonated Ammonia Dimer N2H7+. Angew. Chem., Int. Ed. 2007, 46, 8691−8694. (42) Gardenier, G. H.; Roscioli, J. R.; Johnson, M. A. Intermolecular Proton Binding in the Presence of a Large Electric Dipole: Ar-Tagged Vibrational Predissociation Spectroscopy of the CH3CN·H+·OH2 and CH3CN·D+·OD2 Complexes. J. Phys. Chem. A 2008, 112, 12022− 12026. (43) Douberly, G. E.; Ricks, A. M.; Ticknor, B. W.; Duncan, M. A. Structure of Protonated Carbon Dioxide Clusters: Infrared Photodissociation Spectroscopy and ab Initio Calculations. J. Phys. Chem. A 2008, 112, 950−959. (44) Douberly, G. E.; Ricks, A. M.; Ticknor, B. W.; Duncan, M. A. The Structure of Protonated Acetone and its Dimer: Infrared Photodissociation Spectroscopy from 800 to 4000 cm−1. Phys. Chem. Chem. Phys. 2008, 10, 77−79. (45) Ricks, A. M.; Douberly, G. E.; Duncan, M. A. Infrared Spectroscopy of the Protonated Nitrogen Dimer: The Complexity of Shared Proton Vibrations. J. Chem. Phys. 2009, 131, 104312. (46) Cheng, T. C.; Jiang, L.; Asmis, K. R.; Wang, Y.; Bowman, J. M.; Ricks, A. M.; Duncan, M. A. Mid- and Far-IR Spectra of H5+ and D5+ Compared to the Predictions of Anharmonic Theory. J. Phys. Chem. Lett. 2012, 3, 3160−3166. (47) Cheng, T. C.; Bandyopadhyay, B.; Duncan, M. A.; Mosley, J. D. IR Spectroscopy of Protonation in Benzene-Water Nanoclusters: Hydronium, Zundel and Eigen at a Hydrophobic Interface. J. Am. Chem. Soc. 2012, 134, 13046−13055. (48) McDonald, D. C., II; Mauney, D. T.; Leicht, D.; Marks, J. H.; Tan, J. A.; Kuo, J.-L.; Duncan, M. A. Communication: Trapping a Proton in Argon: Spectroscopy and Theory of the Proton-Bound Argon Dimer and its Solvation. J. Chem. Phys. 2016, 145, 231101. (49) Wolke, C. T.; Fournier, J. A.; Dzugan, L. C.; Fagiani, M. R.; Odbadrakh, T. T.; Knorke, H.; Jordan, K. D.; McCoy, A. B.; Asmis, K. R.; Johnson, M. A. Spectroscopic Snapshots of the Proton-Transfer Mechanism in Water. Science 2016, 354, 1131−1135. (50) Tzeng, W. B.; Wei, S.; Castleman, A. W. Protonated Acetaldehyde Clusters: Stability, Structure and Metastable Unimolecular Decomposition. Chem. Phys. Lett. 1990, 168, 30−36. (51) Amano, T.; Warner, H. E. Laboratory Detection of Protonated Formaldehyde (H2COH+). Astrophys. J. 1989, 342, L99−L101. (52) Dore, L.; Cazzoli, G.; Civiš, S.; Scappini, F. Extended Measurements of the Millimeter Wave Spectrum of H2COH+. Chem. Phys. Lett. 1995, 244, 145−148. (53) Mosley, J. D.; Cheng, T. C.; McCoy, A. B.; Duncan, M. A. Infrared Spectroscopy of the Mass 31 Cation: Protonated Formaldehyde vs Methoxy. J. Phys. Chem. A 2012, 116, 9287−9294. (54) Mauney, D. T.; Mosley, J. D.; Madison, L. R.; McCoy, A. B.; Duncan, M. A. Infrared Spectroscopy and Theory of Formaldehyde 197

DOI: 10.1021/acs.jpca.7b10573 J. Phys. Chem. A 2018, 122, 192−198

Article

The Journal of Physical Chemistry A Cation and its Hydroxymethylene Isomer. J. Chem. Phys. 2016, 145, 174303. (55) Frisch, M. J., Trucks, G. W.; Schlegel, H. B.; Scuseria, G. E.; Robb, M. A.; Cheeseman, J. R.; Scalmani, G.; Barone, V.; Mennucci, B.; Petersson, G. A.; et al. Gaussian 09, Revision D.01; Gaussian, Inc.: Wallingford, CT, 2009. (56) Dunning, T. H., Jr. Gaussian Basis Sets for Use in Correlated Molecular Calculations. I. The Atoms Boron Through Neon and Hydrogen. J. Chem. Phys. 1989, 90, 1007−1023. (57) Holmes, J. L.; Aubrey, C.; Mayer, P. M. Assigning Structures to Ions in Mass Spectrometry; CRC Press: Boca Raton, FL, 2007. (58) Van Dornshuld, E.; Holy, C. M.; Tschumper, G. S. Homogeneous and Heterogeneous Noncovalent Dimers of Formaldehyde and Thioformaldehyde: Structures, Energetics, and Vibrational Frequencies. J. Phys. Chem. A 2014, 118, 3376−3385. (59) Antony, J.; Grimme, S. Is Spin-Component Scaled SecondOrder Møller-Plesset Perturbation Theory an Appropriate Method for the Study of Noncovalent Interactions in Molecules? J. Phys. Chem. A 2007, 111, 4862−4868. (60) Mangiatordi, G. F.; Brémond, E.; Adamo, C. DFT and Proton Transfer Reactions: A Benchmark Study on Structure and Kinetics. J. Chem. Theory Comput. 2012, 8, 3082−3088. (61) Shimanouchi, T. Molecular Vibrational Frequencies. In NIST Chemistry WebBook, NIST Standard Reference Database No. 69; Linstrom, P. J., Mallard, W. G., Eds.; National Institute of Standards and Technology: Gaithersburg, MD, http://webbook.nist.gov. (62) Purvis, G. D., III; Bartlett, R. J. A Full Coupled-Cluster Singles and Doubles Model: The Inclusion of Disconnected Triples. J. Chem. Phys. 1982, 76, 1910−1918. (63) Pople, J. A.; Head-Gordon, M.; Raghavachari, K. Quadratic Configuration Interaction. A General Technique for Determining Electron Correlation Energies. J. Chem. Phys. 1987, 87, 5968−5975. (64) Bü rgi, H. B.; Dunitz, J. D.; Lehn, J. M.; Wipff, G. Stereochemistry of Reaction Paths at Carbonyl Centres. Tetrahedron 1974, 30, 1563−1572. (65) We note that the formation of isomer IIb would still be slightly exothermic with respect to the isolated monomers formaldehyde and protonated formaldehyde. (66) Mehringer, D. M.; Snyder, L. E.; Miao, Y.; Lovas, F. J. Detection and Confirmation of Interstellar Acetic Acid. Astrophys. J. 1997, 480, L71−L74. (67) Brown, R. D.; Crofts, J. G.; Gardner, F. F.; Godfrey, P. D.; Robinson, B. J.; Whiteoak, J. B. Discovery of Interstellar Methyl Formate. Astrophys. J. 1975, 197, L29−L31. (68) DeBlase, A. F.; Bloom, S.; Lectka, T.; Jordan, K. D.; McCoy, A. B.; Johnson, M. A. Origin of the Diffuse Vibrational Signature of a Cyclic Intramolecular Proton Bond: Anharmonic Analysis of Protonated 1,8-Disubstituted Naphthalene Ions. J. Chem. Phys. 2013, 139, 024301. (69) Mauney, D. T.; Maner, J. A.; Duncan, M. A. IR Spectroscopy of Protonated Acetylacetone and Its Water Clusters: Enol-Keto Tautomers and Ion→Solvent Proton Transfer. J. Phys. Chem. A 2017, 121, 7059−7069. (70) Lovas, F. J.; Suenram, R. D.; Coudert, L. H.; Blake, T. A.; Grant, K. J.; Novick, S. E. The Torsional-Rotational Spectrum and Structure of the Formaldehyde Dimer. J. Chem. Phys. 1990, 92, 891− 898. (71) Dolgonos, G. A. Which Isomeric Form of Formaldehyde Dimer is the Most Stable - a High-level Coupled-Cluster Study. Chem. Phys. Lett. 2013, 585, 37−41.

198

DOI: 10.1021/acs.jpca.7b10573 J. Phys. Chem. A 2018, 122, 192−198