Complex between Formic Acid and Nitrous Oxide: A Matrix-Isolation

Oct 25, 2017 - (5, 6) In addition, FA and N2O have been detected in interstellar molecular clouds and are relevant to studies of prebiotic species in ...
0 downloads 6 Views 735KB Size
Subscriber access provided by University of Florida | Smathers Libraries

Article

Complex Between Formic Acid and Nitrous Oxide: A Matrix-Isolation and Computational Study Luis Duarte, Iiris Rekola, and Leonid Khriachtchev J. Phys. Chem. A, Just Accepted Manuscript • DOI: 10.1021/acs.jpca.7b09586 • Publication Date (Web): 25 Oct 2017 Downloaded from http://pubs.acs.org on October 29, 2017

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

The Journal of Physical Chemistry A is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 38

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Complex between Formic Acid and Nitrous Oxide: A Matrix-Isolation and Computational Study Luís Duarte, Iiris Rekola, and Leonid Khriachtchev* Department of Chemistry, University of Helsinki, P.O. Box 55, FI-00014 Helsinki, Finland

ABSTRACT: The complex of formic acid (FA, HCOOH) with nitrous oxide (N2O) was studied experimentally and computationally. Eight structures of the trans-FA···N2O complex and nine structures of the cis-FA···N2O complex were found at the DFT (M06-2X and wB97XD), MP2(full), and CCSD(T)-F12a levels of theory. Two structures of the trans-FA···N2O (1t and 3t) complex and two structures of the cis-FA···N2O (1c and 3c) complex were identified by infrared spectroscopy in an argon matrix. Structure 1t with the bonded OH group appears in the matrices after deposition. Structures 3t and 3c with the free OH groups were prepared by vibrational excitation of the trans-FA conformer in structures 1t and 3t, respectively. Structure 3t is thermally unstable and relaxes to structure 1t at ~15 K. Structure 1c with the bonded OH group is made by vibrational excitation of trans-FA monomer combined with thermal annealing. The lifetimes of the cis-FA···N2O complex structures 1c (81.5 min) and 3c (18.7 min) in an argon matrix at 4.3 K are longer than that of cis-FA monomer (6.3 min). This difference agrees with the calculated stabilization barriers of these species.

ACS Paragon Plus Environment

1

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 2 of 38

1. INTRODUCTION Formic acid (HCOOH, FA) and nitrous oxide (N2O) participate in atmospheric chemistry. Formic acid is one of the major contributors to rainwater acidity in polluted and remote regions.12

Nitrous oxide has a global warming potential nearly 300 times that of carbon dioxide and is one

of the largest contributors to ozone depletion in the stratosphere.3-4 Recent studies indicate that the FA and N2O annual emissions from known sources may have been significantly underestimated and suggest that these species might be more important to atmospheric chemistry than previously believed.5-6 In addition, FA and N2O have been detected in interstellar molecular clouds and are relevant to studies of prebiotic species in astrochemistry.7-9 The interest in noncovalent interactions of FA and N2O is considerable.10-21 The formation of weakly bound complexes can cause changes in the stability and reactivity of the constituent monomers.22-24 Quantum mechanical tunneling observed in FA is an exciting phenomenon that can affect the dynamics of the reactions and open new chemical reaction channels.25-26 Matrix-isolation infrared spectroscopy is a powerful method to study weakly bonded species,27-29 and a fair number of the FA and N2O complexes have been prepared in cryogenic matrices.30-45 The conformational change is another successful direction of research at low temperatures. The higher-energy conformer of formic acid, cis-FA, can be prepared in cryogenic matrices by selective vibrational excitation of the ground-state conformer, trans-FA (see Scheme 1 for the structures of these conformers).46-50 The cis-FA molecules isolated in cryogenic matrices are not stable and they convert back to trans-FA by hydrogen-atom tunneling.49-51 The cis-FA complexes can be made by using selective vibrational excitation of trans-FA often combined with annealing of the matrix. In this way, the complexes of cis-FA with H2O, N2, CO2, Xe, and several trans-cis and cis-cis (the latter of cis-HCOOD) dimers have been prepared.33-35,

ACS Paragon Plus Environment

2

Page 3 of 38

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

38, 43, 52-54

This complexation can stabilize the higher-energy conformer. Remarkably, the

conformation-dependent FA + O reaction was found in a krypton matrix.55

Scheme 1. Structures of trans- and cis-FA.

In the present work, we study the complexes of the trans-FA and cis-FA conformers with N2O. These complexes are identified by using infrared spectroscopy in an argon matrix supported by extensive quantum chemical calculations. Selective vibrational excitation of transFA is used to prepare higher-energy structures of the FA···N2O complex. Hydrogen-atom tunneling decay of the cis-FA···N2O complex is measured and discussed based on the calculated stabilization barriers.

2. MATERIALS AND METHODS 2.1. Computational details. The DFT (M06-2X56-57 and wB97XD58 functionals) and MP2(full)59 quantum chemical calculations were carried out with the Gaussian 16 program.60 The calculations employed a very tight optimization convergence criteria and the standard augcc-pVTZ basis set.61 The DFT calculations used the Gaussian 16 default ultrafine integration grid. The explicitly correlated CCSD(T)-F12a calculations were performed with the MOLPRO program.62-63 The cc-pVDZ-F12 (VDZ-F12) correlation consistent basis set was used together with the recommended value of the geminal Slater exponent (0.9)64 and with the Molpro default

ACS Paragon Plus Environment

3

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 4 of 38

auxiliary basis sets. The geometry optimizations were followed by harmonic frequency calculations at the same levels of theory, which also gave the zero-point vibrational energies (ZPVE) and verified the nature of the obtained minima. The interaction energies were calculated as the difference of the total energies of the complex and isolated species (with the geometries in the complex) and corrected for the ZPVE and the basis set superposition error (BSSE).65

2.2. Experimental details. The FA/N2O/Ar mixtures were prepared with concentration ratios of ~1/1/1000 and 1/3/1000. FA (≥98%, Merck) was degassed by several freeze-pump-thaw cycles. N2O (99.9999%, AGA) and argon (≥99.9999%, AGA) were used as supplied. Prior to the mixture preparation, the glass bulb was passivated with FA vapors by several fill–keep–evacuate cycles. The gas mixtures were deposited onto a CsI window held at 15 K in a closed-cycle helium cryostat (RDK-408D2, Sumitomo Heavy Industries, Ltd.). The matrix thickness was ~100 µm. The FTIR spectra in the 4000−500 cm−1 range were measured at 4.3 K with a Bruker Vertex 80 spectrometer by co-adding 500 scans at a spectral resolution of 1 cm−1. In the kinetic measurements, a long-pass filter that transmitted below 1850 cm−1 was inserted between the Globar source and the cryostat in order to eliminate the light-induced cis-to-trans conversion.51 The Globar light was blocked between the measurements. Conformational changes were promoted by an optical parametric oscillator (LaserVision) that produced infrared light with a pulse duration of ~5 ns, line width of ~0.1 cm−1, and repetition rate of 10 Hz.

3. RESULTS 3.1. Computational results. To probe the potential energy surface of the complexes, several initial geometries were analyzed at the DFT (M06-2X and wB97XD) and MP2 (full)

ACS Paragon Plus Environment

4

Page 5 of 38

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

levels of theory. Each of the stable minima was then optimized at the CCSD(T)-F12a level. As a result, eight trans-FA···N2O and nine cis-FA···N2O structures were found at different levels of theory. The calculated complex structures with the interaction energies are presented in Figures 1 and 2 (see also Tables S1 and S2 in the Supporting Information). Structures 1t/8t, 2t/7t, 3t/5t, and 4t/6t of the trans-FA···N2O complex are stabilized by O−H···O, O−H···N, C−H···O, and C−H···N interactions, respectively. With the exception of 7t (C1 symmetry), all other transFA···N2O structures are planar (Cs symmetry). The interaction energies of the three most stable structures 1t, 2t, and 3t, are −15.9, −10.6, and −9.6 kJ mol−1 (after ZPVE and BSSE corrections) at the CCSD(T)-F12a/VDZ-F12 level of theory. The O−H···O, O−H···N, and C−H···O bond lengths for structures 1t, 2t, and 3t are 2.02, 2.20, and 2.64 Å (Table 1). Two trans-FA···N2O structures, similar to 1t and 2t, were reported by Solimannejad et al.16 Structures 1c, 2c, 3c/5c, and 6c of the cis-FA···N2O complex are stabilized by O−H···O, O−H···N, C−H···O, and C−H···N interactions, respectively. The interaction energies of the most stable structures 1c, 2c, and 3c are −9.5, −8.6, and −10.4 kJ mol−1 (after ZPVE and BSSE corrections) at the CCSD(T)-F12a/VDZ-F12 level of theory. The O−H···O, O−H···N, C−H···O bond lengths for 1c, 2c, and 3c are 2.04, 2.16, and 2.60 Å. Structure 7c is an energy minimum only at the M06-2X/aug-cc-pVTZ level and structures 8c and 9c are energy minima only at the MP2(full)/aug-cc-pVTZ level.

ACS Paragon Plus Environment

5

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 6 of 38

Figure 1. Structures of the trans-FA···N2O complex optimized at the CCSD(T)-F12a/VDZ-F12 level of theory. The interaction energies (in kJ mol−1; after ZPVE and BSSE corrections) are given in parentheses. The geometrical parameters r, θ, and β represent the shortest distance and the angles between the complex units and are given in Table 1.

ACS Paragon Plus Environment

6

Page 7 of 38

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Figure 2. Structures of the cis-FA···N2O complex optimized at the CCSD(T)-F12a/VDZ-F12 level of theory. The interaction energies (in kJ mol−1; after ZPVE and BSSE corrections) are given in parentheses. Structures 7c and 8c/9c are not energy minima at this level and are given at the M06-2X/aug-cc-pVTZ and MP2(full)/aug-cc-pVTZ levels (denoted as n.m.). The geometrical parameters r, θ, and β represent the shortest distance and the angles between the complex units and are given in Table 1.

ACS Paragon Plus Environment

7

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 8 of 38

Table 1. Selected geometrical parameters of the optimized structures of the trans-FA···N2O and cis-FA···N2O complexes.a

a

trans-FA···N2O

r

θ

β

cis-FA···N2O

r

θ

β

1t

2.02

166.2

166.9

1c

2.04

171.0

123.2

2t

2.20

161.0

121.1

2c

2.16

178.8

170.6

3t

2.64

108.9

113.8

3c

2.60

110.4

113.5

4t

2.79

106.9

111.4

4c

3.02

94.0

84.7

5t

2.69

117.7

113.6

5c

2.96

94.1

85.7

6t

2.89

121.3

89.8

6c

2.74

108.8

111.8

7t

2.87

110.9

89.6

7c

2.10

135.4

106.4

8t

2.09

153.1

119.1

8c

2.94

98.7

90.1

9c

2.95

88.5

90.2

At the CCSD(T)-F12a/VDZ-F12 level except structures 7c and 8c/9c that are at the M06-

2X/aug-cc-pVTZ and MP2(full)/aug-cc-pVTZ levels. The shortest distance between the monomers (r) is in angstroms and the angles (θ and β) are in degrees.

The complexation-induced spectral shifts of the trans- and cis-FA complexes at CCSD(T)-F12a level of theory are given in Tables 2 and 3. The most significant shifts are predicted for the OH stretching and COH torsional modes of structures 1t, 2t, 8t, 1c, and 2c. The red shift of the OH stretching mode and the blue shift of the COH torsional mode are caused by the O−H···O and O−H···N interactions in structures 1t/1c/8t, and 2t/2c, respectively. For the structures with the C−H···O interactions (3t, 5t, 3c, and 5c), the largest shifts are obtained for the CH stretching mode. For the structures with the C−H···N interactions (4t, 6t, and 6c), the shifts are rather small. The vibrational frequencies for the monomers and the trans-FA···N2O and cisFA···N2O complexes at the M06-2X, wB97XD, MP2(full), and CCSD(T)-F12a levels of theory

ACS Paragon Plus Environment

8

Page 9 of 38

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

are given in Tables S3, S4, and S5 in the Supporting Information. It is seen that the different levels of theory predict similar spectral changes in the complexes. An exception occurs for the NN and NO stretching modes of N2O at the MP2(full) level of theory, which gives contradictory results compared to the DFT (M06-2X and wB97XD) and CCSD(T)-F12a levels (Tables S4 and S5).

Table 2. Calculated and experimental vibrational spectra (cm−1) of trans-FA and N2O monomers and the trans-FA···N2O complex.a monomers exptlb

calcd

tB

1t

2t

3t

shift exptl

4t

5t

6t

7t

8t

shift calcd

assignment

3550.5

3756.7

−82.3

−5.9

−77.3

−44.8

−1.7

0.3

−1.2

−0.1

−12.6

−43.5

OH stretch.

2953.0

3089.1

−84.7 −8.8

−6.2

−0.7

−1.1

14.3

6.3

13.4

0.6

8.0

−7.4

CH stretch.

2218.5

2281.9

14.5

10.7

14.6

14.4

9.6

7.4

8.4

7.0

3.1

2.8

NN stretch.

1767.4

1808.9

−11.5

−6.7

−14.9

−10.8

−10.1

−6.5

0.5

2.0

−3.3

−4.4

C=O stretch.

1381.0

1410.1

33.8

−2.6

9.6

7.2

−1.5

9.9

−1.3

−0.2

3.4

6.0

CH rock.

1282.9

1298.1

−20.6

−3.3

−18.9

22.2

−3.9

11.7

−0.4

6.2

3.9

−13.8

NO stretch.

1305.6c

1315.5

29.5e

6.1

34.3

18.6

4.4

10.6

−7.7

−6.4

2.3

14.9

CO−COH def.

1215.4c

a

tU

9.0

1103.6

1134.6

28.1

7.5

31.1

18.6

5.7

9.7

−8.2

−6.8

1.6

10.1

COH−CO def.

1038.2

1052.2

1.8

2.8

7.0

5.3

8.7

13.2

2.5

0.5

1.2

4.3

CH wagg.

635.4

670.5

70.1

6.1

80.5

49.2

6.3

14.9

−2.8

−2.6

12.4

56.7

COH tors.

629.2

630.7

18.4

2.3

20.8

13.8

3.6

10.0

−0.6

−0.8

5.6

12.9

OCO bend.

589.0

594.6 d

−9.5

−5.4

−8.6

7.6

−1.1

7.2

−0.7

3.0

5.1

−0.9

NNO bend.



594.6 d





−9.6

6.6

−5.9

2.8

−2.0

1.6

−1.3

−7.1

NNO bend.

At the CCSD(T)-F12a/VDZ-F12 level of theory. The shifts are the difference between the

complex and monomer frequencies. The shifts of the strongest absorptions are given in bold.

b

ACS Paragon Plus Environment

9

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 10 of 38

trans-FA in site 2 according to Ref. 48. c Fermi resonance doublet between the CO−COH def. fundamental mode and the first overtone of OH torsional mode. dDoubly degenerated mode.

e

Calculated relatively to the average of the Fermi resonance values. Table 3. Calculated and experimental vibrational spectra (cm−1) of cis-FA and N2O monomers and the cis-FA···N2O complex.a monomers exptlb

calcd

3615.8

3819.9

cB

cU

1c

2c

shift exptl −20.0

3c

4c

5c

6c

shift calcd 1.9

−44.4

−26.3

−0.9

assignment 0.2

0.0

−1.2 OH stretch.

−30.4 2896.4

3004.3

−1.2

−6.3

−0.9

−2.5

19.6

13.3

16.0

7.7 CH stretch.

−17.5 2218.5

2281.9

9.8

10.8

7.8

22.2

10.0

7.6

7.4

8.3 NN stretch.

7.4 1807.0

1850.6

−6.4

−5.1

−5.0

−4.5

−9.8

−4.1

−7.7

−7.9

C=O stretch.

1391.9

1422.3



−2.2

2.6

2.3

0.7

−0.3

6.3

2.9

CH rock.

1282.9

1298.1

−6.2

−4.7

−16.9

22.6

−6.9

−2.2

−7.7

13.8

NO stretch.

1248.9

1293.0

16.1

9.3

38.7

34.9

3.8

−0.5

5.2

2.1

4.6

−1.7

CO−COH def.

−6.0

a

1107.8

1118.0



6.6

22.9

21.4

7.5

2.5

9.3

7.5

COH−CO def.



1030.7





6.8

6.1

9.5

3.8

4.1

7.3

CH wagg.

505.1

520.6

43.1

8.0

97.8

96.3

11.4

6.0

9.7

8.8

COH tors.

662.3

660.3

5.7

2.9

13.8

10.2

3.9

4.2

4.3

3.5

OCO bend.

589.0

594.6 c

−18.7

−5.6

−4.4

3.7

−1.9

−1.1

−3.0

5.8

NNO bend.



594.6 c

−11.3

1.2

−7.1

−3.6

−5.0

1.1

NNO bend.





At the CCSD(T)-F12a/VDZ-F12 level of theory. Complex 7c, 8c, and 9c are not energy minima

at this level. The shifts are the difference between the complex and monomer frequencies. The

ACS Paragon Plus Environment

10

Page 11 of 38

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

shifts of the strongest absorptions are given in bold.

b

cis-FA in site 2 according to Ref. 48.

c

Doubly degenerated mode.

The cis-to-trans torsional energy scans of the cis-FA monomer and cis-FA···N2O complex (structures 1c and 3c) are shown in Figure 3. The energies were calculated at the CCSD(T)-F12a/VDZ-F12 level by changing the torsional angle from the initial torsional angle of 180 degrees (cis-FA) to 0 degrees (trans-FA) with an increment of 10 degrees. In the torsional scan, the coordinates of all other atoms are frozen, which corresponds to the adiabatic approximation for fast tunneling.66 The energies are not corrected for ZPVE or BSSE. The performed simple torsional scans is known to slightly overestimate the barrier height.67 The cisto-trans barriers of 3109, 4115, 4085, 3257, and 3232 cm−1 were obtained for the cis-FA monomer and structures 1c, 2c, 3c, and 6c, respectively. For structures 4c and 5c, two different barriers exist for the torsional movement (i) toward the N2O unit and (ii) away from the N2O unit. The cis-to-trans barriers of structures 4c and 5c are 3088/3152 cm−1 and 3095/3292 cm−1 for the toward/away torsional movements. It is worth comparing the calculated properties of the FA···N2O complex and the isoelectronic FA···CO2 complex.33 The most stable structures of the trans-FA···CO2 and cisFA···CO2 complexes are formed by interaction of the OH group of FA with an oxygen atom of CO2. We have calculated these structures at the CCSD(T)-F12a/VDZ-F12 level used in the present work. The interaction energies of these structures (−15.7 and –10.4 kJ mol−1; ZPVE and BSSE corrected) are similar to those of structures 1t (trans-FA···N2O, −15.9 kJ mol−1) and 1c (cis-FA···N2O, −9.5 kJ mol−1). The O−H···O bond lengths are 2.04 and 2.07 Å for the strongest trans-FA···CO2 and cis-FA···CO2 complexes. These lengths are slightly larger than those of

ACS Paragon Plus Environment

11

The Journal of Physical Chemistry

structures 1t (trans-FA···N2O; 2.02 Å) and 1c (cis-FA···N2O; 2.04 Å). The shifts predicted for the OH stretching/COH torsional modes of the trans-FA···CO2 and cis-FA···CO2 complexes are −56/+75 and −19/+81 cm−1, respectively. The shifts obtained for structures 1t (trans-FA···N2O) and 1c (cis-FA···N2O) are similar (−77/+81 and −44/+98 cm−1, respectively), although some difference occurs for the OH stretching mode of 1c.

5000

Relative energy (cm−1)

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 12 of 38

cis-FA 1c 3c

4000 3000 2000 1000 0 0

30

60

90

120

150

180

Torsional angle (degree) Figure 3. CCSD(T)-F12a/VDZ-F12 energy scans for the cis-to-trans reorganization of FA monomer (red) and structures 1c (green) and 3c (black). The torsional angles of 0° and 180° correspond to the trans- and cis-FA conformations.

3.2 Experimental results. After deposition of an FA/N2O/Ar matrix, a set of absorptions is observed that do not belong to FA and N2O monomers or their multimers. These bands are consistent with the formation of the 1:1 complex between these two species. Figure 4a presents

ACS Paragon Plus Environment

12

Page 13 of 38

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

selected regions of the IR spectrum of an as-deposited FA/N2O/Ar matrix. For the O−H and C=O stretching modes, the trans-FA···N2O complex bands (marked as tB) are observed at 3468.2/3465.8 and 1755.9 cm−1, together with the known absorptions of the trans-FA monomer (3550 and 1767 cm−1)48 and trans-FA dimers (tt1: 1728 cm−1; tt2: 1748 cm−1 following Ref. 54). Other bands of tB are observed for the NN stretching (2218.5 cm−1), CO–COH deformation (1290.0 cm−1), NO stretching (1262.3 cm−1), COH–CO deformation (1131.7 cm−1), COH torsional (705.5 cm−1), and OCO bending (647.6 cm−1) modes. After annealing at 35 K, the monomer bands of trans-FA and N2O decrease in intensity and the tB bands increase (Figure 4b). The complexation-induced shifts of the OH stretching and torsional modes of tB from those of trans-FA monomer are significant. Based on this fact, we experimentally assign the tB bands to the trans-FA···N2O complex, in which N2O interacts with the OH group of trans-FA. The experimental shifts of the tB bands with respect to the monomer are given in Table 2. Vibrational excitation at 3466 cm−1 (the strongest component of the OH stretching mode of tB after annealing at 35 K) at 4.3 K consumes the tB bands and produces two new sets of bands, denoted here as tU (major channel) and cU (minor channel) (Figure 4c). The most characteristic bands of tU belong to the OH stretching (3544.6 cm−1), NN stretching (2229.2 cm−1), C=O stretching (1760.6 cm−1), NO stretching (1279.6 cm−1), COH−CO deformation (1111.1 cm−1), and COH torsional (641.5 cm−1) modes. The complexation-induced shifts of OH stretching and COH torsional modes of tU from those of trans-FA monomer are rather small. Based on this fact, we experimentally assign the tU bands to the trans-FA···N2O complex without an interaction between N2O and the OH group of trans-FA. The experimental shifts of the tU species are given in Table 2. The cU species is a minor product of the 3466 cm−1 pumping of tB and it is discussed below.

ACS Paragon Plus Environment

13

The Journal of Physical Chemistry

Excitation at 3544 cm−1 (the OH stretching frequency of tU) at 4.3 K consumes tU and produces new bands (denoted earlier as cU) at 3617.8, 2229.4, 1801.9, 1278.2, 1258.2/1247.2, and 513.2 cm−1 (Figure 4d). Based on the characteristic shifts, we experimentally assign the cU bands to the cis-FA···N2O complex without an interaction between N2O and the OH group of cis-FA. Excitation of the OH stretching mode of tU also produces a small amount of tB. The experimental shifts of the cU species are given in Table 3.

trans-FA 0.6

Absorbance

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 14 of 38

trans-FA

0.4

tB

(a) tB (b) tU

0.2 cU (c)

cU

cU

tU

tB

tB cU

(d) x3 tU

0.0

tU 3600

3550

3500

1800

1750

−1

Wavenumber (cm )

Figure 4. FA···N2O complex in an argon matrix. FTIR spectra of (a) as deposited FA/N2O/Ar (1/1/1000) matrix; (b) after annealing at 35 K; (c) result of excitation at 3466 cm−1 at 4.3 K (difference spectrum); (d) result of excitation at 3544 cm−1 at 4.3 K (difference spectrum). The spectra were measured at 4.3 K.

ACS Paragon Plus Environment

14

Page 15 of 38

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 ACS Paragon Plus Environment

15

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 16 of 38

As pointed out above, vibrational excitation of tB at 4.3 K efficiently produces structure tU (Figure 5a). Upon annealing of the matrix at ~15 K, tU (prepared by vibrational excitation of tB at 4.3 K) converts back to tB at a subminute scale (Figure 5b). This annealing-induced conversion confirms that the tB and tU species are two structures of the trans-FA···N2O complex, tB (with the bonded OH group) being the more stable one. Figure 5c shows the result of vibrational excitation at 3466 cm−1 of a FA/N2O/Ar matrix performed at 15 K. Under these conditions, cU is mainly produced from tB. After stopping irradiation, cU decays mostly to tU (Figure 5d). As previously reported, selective vibrational excitation of the OH stretching mode of trans-FA at 3550 cm−1 efficiently produces cis-FA molecules (Figure 6a).46-50 If the matrix is annealed at 35 K during this excitation, a new set of bands denoted as cB in Figure 6b is produced (3595.9/3585.0, 2228.3/2225.9, 1800.5, 1276.7, 1265.0/1253.5/1242.9, and 548.2). In time, cB slowly decays to tU (Figure 6c). Based on the characteristic shifts, we assign the cB bands to the cis-FA···N2O complex, in which N2O interacts with the OH group of cis-FA. In the same way, the H-bonded complexes of cis-FA with CO2, N2, Xe, and H2O have been prepared.3335, 38

The experimental shifts of the cB species are given in Table 3.

ACS Paragon Plus Environment

16

Page 17 of 38

0.8

tU

tU cU

cU (a) 0.6

tB Absorbance

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

tB

tB 0.4

(b) cU

tU cU 0.2

(c) tB (d)

0.0

tU

tU cU

cU 3600

tB

3550

3500

1800

1750

−1

Wavenumber (cm )

Figure 5. FA···N2O complex in an argon matrix. Difference FTIR spectra showing the results of (a) excitation at 3466 cm−1 at 4.3 K of a FA/N2O/Ar (1/1/1000) matrix annealed at 35 K (similar to spectrum c in Figure 4); (b) annealing at 15 K of the previous matrix; (c) excitation at 3466 cm−1 at 15 K of a FA/N2O/Ar (1/3/1000) matrix annealed at 35 K after deposition; (d) waiting for ~50 min under Globar irradiation at 4.3 K for the previous matrix. The spectra were measured at 4.3 K.

ACS Paragon Plus Environment

17

The Journal of Physical Chemistry

cis-FA 0.8

cis-FA

trans-FA

1.0

trans-FA 0.6

Absorbance

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 18 of 38

(a)

*

tB (b) x5

0.4

tB

(b) x5 0.2

.

(b) tU tU

0.0

(c) x5 cB

cB -0.2 3600

3550

3500

1800

1750

Wavenumber (cm−1)

Figure 6. FA···N2O complex in an argon matrix. FTIR spectra (a) after excitation at 3550 cm−1 at 4.3 K of a FA/N2O/Ar (1/3/1000) matrix; (b) after annealing at 35 K of a FA/N2O matrix upon excitation at 3550 cm−1; (c) the result of waiting for ~70 min under Globar irradiation at 4.3 K for the previous matrix (difference spectrum). The spectra were measured at 4.3 K. The bands marked with an asterisk and a dot are from the C=O stretching overtone of cis-FA monomer and the C=O stretching mode of a trans-cis-FA dimer (tc3 according to the notation of Ref. 54), respectively.

ACS Paragon Plus Environment

18

Page 19 of 38

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

The “dark” decay curves of two structures of the cis-FA···N2O complex and of cis-FA monomer are shown in Figure 7a. These data were obtained by measuring through a long-pass filter (1850 cm−1), which suppressed the effect of the spectrometer beam on the decay rate. The C=O stretching band of each species were integrated and the results were fitted with a single exponential function. The lifetimes of the cU and cB species and cis-FA monomer at 4.3 K are 18.7, 81.5, and 6.3 min, respectively. Our data for the decay of cis-FA monomer agrees with the literature.51 The decays of cU and cis-FA monomer were also measured at elevated temperatures (Figure 7b). The data show the low-temperature limit of the reaction, which is typical for H-atom tunneling.68 To recall, the activation of the reaction at higher temperatures (at ~20 K) is not an over-barrier process but explained in terms of matrix-reorganization energy.66 The temperature dependence of the cB was not measured because its band intensities are relatively weak and every data point required the preparation of a new matrix.

ACS Paragon Plus Environment

19

The Journal of Physical Chemistry

Relative concentration

1.0

(a) cB

0.5

cU

cis-FA 0.0 0

10

20

30

40

50

Time (min) (b)

-5

cis-FA ln k

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 20 of 38

-6

cU

-7

0.05

0.10

0.15

0.20

0.25

1/T (K−1) Figure 7. (a) “Dark” decay of cis-FA, cU, and cB in an argon matrix at 4.3 K. The following bands were integrated: 1807.0 cm−1 (cis-FA, squares); 1801.9 cm−1 (cU, down-triangles); 1800.5 cm−1 (cB, circles). The lines are single exponential fits. (b) Arrhenius plots for the decay rate of the cis-FA and cU species. The lines guide the eye.

ACS Paragon Plus Environment

20

Page 21 of 38

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

4. DISCUSSION 4.1. Spectral assignment. The set of bands observed after deposition and increasing upon annealing (marked as tB in Figure 4) is assigned to structure 1t of the trans-FA···N2O complex (Table 2). The observed shifts of the OH stretching and COH torsional modes are consistent with the presence of the O−H···O and O−H···N interactions in 1t and 2t. However, structure 2t can be ruled out, because the calculated shift for the NO stretching mode (+22.2 cm−1) is of the opposite sign than the experimental shift (−20.6 cm−1) observed for tB. On the other hand, all calculated values obtained for 1t are in good agreement with the experimental shifts of tB (Table 2). The OH stretching, CH stretching, NN stretching, C=O stretching, NO stretching, COH−CO deformation, and COH torsional bands of this complex are shifted from the monomer bands, respectively, by −82.3/−84.7, –8.8, +14.5, −11.5, −20.6, +28.1 and +70.1 cm−1 in the experiment (tB) and by −77.3, –0.7, +14.6, −14.9, −18.9, +31.1 and +80.5 cm−1 in the calculations (1t). The CO−COH deformation mode of trans-FA monomer is split by Fermi resonance with the first overtone of the COH torsional mode whereas this mode of the tB complex is free of Fermi resonance. Therefore, an accurate comparison of the experimental and calculated shifts is not possible for this mode. However, the averaged shift from the Fermi components (+29.5 cm−1) agrees well with the calculated value (+34.3 cm−1). The experimental species tU is assigned to structure 3t of the trans-FA···N2O complex. The small shifts of the OH stretching and COH torsional modes show weak perturbation of the O−H bond as in structures 3t and 4t. Structure 4t can be ruled out based on the NO stretching shifts (theory: +11.7 cm−1; experiment: −3.3 cm−1). The OH stretching, NN stretching, C=O stretching, NO stretching, COH−CO deformation, and COH torsional modes of tU are shifted, respectively, by −5.9, +10.7, −6.7, −3.3, +7.5 and +6.1 cm−1, which is in a good agreement with

ACS Paragon Plus Environment

21

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 22 of 38

the calculations for structure 3t (−1.7, +9.6, −10.1, −3.9, +5.7 and +6.3 cm−1) (Table 2). In the CO−COH deformation region, the observation of two bands shifted from the Fermi resonance components of the monomer by +6.1 and +9.0 cm−1 suggest that CO−COH mode in the complex is also split. One can notice that the CH stretching mode is shifted in the experiment and in the calculations by −6.2 and +14.3 cm−1, respectively. This mismatch can be due to changes in interaction with the matrix (solvation) for this mode involved in the C–H···O interaction. It has been suggested that in this case, it is better to compare the complex in a matrix with the monomer in a vacuum.69-70 The CH stretching band of trans–FA in the gas phase is at 2942.8 cm−1,71 indicating a shift of +4 cm−1 which is closer to the calculated value. Furthermore, the geometry of the weakly bound complex can be modified by interactions with the surrounding argon atoms. Indeed, the calculations are performed for the species in a vacuum, whereas the experiment is made in a polarizable matrix. For the experimental species cU, the OH stretching, NN stretching, C=O stretching, NO stretching, COH–CO deformation, and COH torsional modes are shifted by +1.9, +10.8, −5.1, −4.7, +9.3/−1.7 and +8.0 cm−1, respectively. Based on the calculated shifts, structures 3c, 4c, and 5c are possible for this species (Table 3). However, the IR-induced tU → cU and tunnelinginduced cU → tU processes are dominant, the tU complex having structure 3t. From the geometrical considerations, structure 3c is the most suitable candidate for the experimental complex cU; therefore, we assign it accordingly. The mismatch between the experimental and calculated shifts of the CH stretching mode can be explained by the same rationale used above for tU. It is interesting to notice that the OH stretching mode of cU is slightly blue-shifted. To our knowledge, it is the first example of an experimental blue shift of this mode of FA upon complexation.

ACS Paragon Plus Environment

22

Page 23 of 38

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

The cB experimental species is assigned to structure 1c of the cis-FA···N2O complex. The shifts observed for the OH stretching and COH torsional modes suggest the O−H···O or O−H···N interactions as in 1c or 2c. However, structure 2c can be ruled out based on the NO stretching mode (Table 3). The experimental shifts of the OH stretching, NN stretching, C=O stretching, NO stretching, COH−CO stretching, and COH torsional modes of cB (−20.0/−30.4, +9.8/+7.4, −6.4, −6.2, +16.1/+4.6/−6.0 and +43.1 cm−1) are in a reasonable agreement with the calculations for 1c (−44.4, +7.8, −5.0, −16.9, +38.7 and +97.8 cm−1) (Table 3). The splitting of the OH stretching, CH stretching, NN stretching, and CO−COH deformation bands is presumably due to matrix site effects. As previously, the agreement in less satisfactory for the CH stretching mode. 4.2. Tunneling reaction. Figure 7a shows that the “dark” decay of species cB (lifetime of 81.5 min) and cU (lifetime of 18.7 min) is substantially slower than that of cis-FA monomer (lifetime of 6.3 min). The decay of cis-FA monomer as well as of a number of other complexes of cis-FA occurs via H-atom tunneling through the torsional barrier, as discussed elsewhere.49-50, 67

We suggest that the decay of the cB and cU species also occurs by H-atom tunneling. Several

factors can affect the tunneling rates in low-temperature rare-gas matrices. The most important factor controlling the decay of cis-carboxylic acids is the cis-to-trans barrier.49-52 The cis-to-trans torsional scans of cis-FA monomer and structures 1c and 3c of the cis-FA···N2O complex show barrier heights of 3109, 4115, 3257 cm−1 for these species, respectively (Figure 3). Thus, the calculated barriers are in full agreement with the order of the experimental decay rates. One experimental fact is difficult to fully understand. The cis-to-trans barrier for the cisFA···CO2 complex (4089 cm−1 at the CCSD(T)F12a/VDZ-F12 level) is close to that of structure 1c of the cis-FA···N2O complex with a similar arrangement of the complex units (4115 cm−1).

ACS Paragon Plus Environment

23

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 24 of 38

Despite the similar stabilization barriers, the lifetime of cis-FA···CO2 complex was reported to be 13.3 h at 10 K,33 which is an order of magnitude longer than that of structure 1c of the cisFA···N2O complex. We can suggest three possible reasons for this mismatch. First, it has been found that the vibrational levels behind the torsional barriers are important for the tunneling process due to energy conservation.51, 72 The energy mismatch between the energy of the initial configuration should not differ much from the energy of the accepting level (in the scale of the Debye frequency of the matrix). The interaction energies in these two complexes somewhat differ from each other, which changes the energy mismatches. Second, the local matrix morphology can be different for these two complexes. As we have recently discussed, the true reaction coordinate is not exclusively the hydrogen motion but rather the collective coordinate involving the reorganization of matrix molecules.73 Finally, as we pointed out above, the calculations are performed for the species in a vacuum, whereas the experiment is made in a polarizable matrix. The geometry of the weakly bound complexes can be modified by interactions with the surrounding argon atoms and the matrix solvation effect can differ between these two complexes. Despite these possibilities, the relatively low stability of the 1c structure of the cis-FA···N2O complex is surprising, in our opinion. Explanation of this observation is beyond the scope of the present work.

ACS Paragon Plus Environment

24

Page 25 of 38

infrared 3c

*

annealing

3t infrared

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

tunneling

* 1c

1t

Figure 8. Main processes observed in the present work. For the infrared excitation of 1t and 3t, the processes marked with an asterisk have a relatively low efficiency compared to the second channels.

5. CONCLUSIONS In the present work, we have studied the interaction of formic acid (FA, HCOOH) with nitrous oxide (N2O). FA and N2O participate in atmospheric chemistry, where non-covalent interactions play an important role. These two molecules appear in the list of interstellar species. This fact is also a motivation because the detection of weak complexes in interstellar space is a challenge and the laboratory experiments on relevant species are needed.74 The studied complex has a number of interesting fundamental features discussed below. The geometries, interaction energies, and vibrational spectra of the FA···N2O complex were calculated at the DFT (M06-2X

ACS Paragon Plus Environment

25

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 26 of 38

and wB97XD), MP2(full), and CCSD(T)-F12a levels of theory. Eight structures of the transFA···N2O complex and nine structures of the cis-FA···N2O complex were obtained at these levels (Figures 1 and 2). Two structures of the trans-FA···N2O (1t and 3t) complex and two structures of the cis-FA···N2O (1c and 3c) complex were identified by infrared spectroscopy in an argon matrix. The calculated and experimental spectra of these species are in good agreement (Tables 2 and 3). The observed processes are summarized in Figure 8. Structure 1t of the trans-FA···N2O complex is observed after deposition of a FA/N2O/Ar matrix and its amount increases upon annealing. Excitation of the OH stretching mode of 1t at 4.3 K mainly leads to structure 3t of the trans-FA···N2O complex, i.e. the O−H···O interaction in 1t reorganizes to the C−H···O interaction in 3t. This process probably occurs by a flip of the trans-FA unit. A minor formation of structure 3c with the C−H···O interaction is observed at this pumping, which means the transto-cis conversion with a subsequent reorganization of the complex. Excitation of the OH stretching mode of 3t mainly leads to structure 3c of the cis-FA···N2O complex, i.e. the C−H···O interaction retains and the free OH group reorganizes to the cis configuration. Structure 3c is spectroscopically unusual because it shows a blue shift of the free OH stretching mode; however, we have no strong theoretical support to this experimental fact. Structure 3t of the trans-FA···N2O complex reorganizes to the more stable structure 1t upon annealing at 15 K at sub-minute scale. This decay indicates a very low stabilization barrier of 3t, which is unusual for these kinds of complexes. If excitation of the OH stretching mode of 1t is performed at 15 K, mainly structure 3c of the cis-FA···N2O complex is formed. This temperature-induced switching of the IR-pumping product is a remarkable observation. Structure 3c decays at ≤10 K exclusively to structure 3t of the trans-FA···N2O complex even in the dark at

ACS Paragon Plus Environment

26

Page 27 of 38

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

a scale of about 20 min presumably by H-atom tunneling. When the decay of 3c is monitored at 15 K and above, exclusively the formation of structure 1t is observed. Of course, it is not a true temperature-induced switching of the tunneling decay product, because the 3c → 3t → 1t cascade process presumably occurs: the first step operates by H-atom tunneling and the second one by thermal relaxation. Annealing of the matrices containing both cis-FA and N2O produces the cis-FA···N2O complex having structure 1c with the O−H···O interaction. This complex decays at 4.3 K in the dark to structure 3t of the trans-FA···N2O complex by H-atom tunneling. In this process, the O−H···O interaction in 1c breaks and the C−H···O interaction in 3t forms. The remaining question is why the lifetime of 1c is much shorter than that of the cis-FA···CO2 complex with the similar stabilization barrier and the structural motif. No 3c → 1c process was observed upon annealing at 20 K. This observation shows that the 3c → 1c stabilization barrier is higher that the 3t → 1t stabilization barrier.

ASSOCIATED CONTENT Supporting Information. Calculated interaction energies of the trans-FA···N2O and cisFA···N2O complexes (Table S1), selected geometrical parameters of the optimized transFA···N2O and cis-FA···N2O complexes (Table S2), calculated vibrational frequencies and infrared intensities of trans-FA, cis-FA, and N2O monomers (Table S3), calculated vibrational frequencies and infrared intensities of the trans-FA···N2O complex (Table S4), calculated vibrational frequencies and infrared intensities of the cis-FA···N2O complex (Table S5), and full references.

ACS Paragon Plus Environment

27

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 28 of 38

AUTHOR INFORMATION Corresponding Author *E-mail: [email protected] Notes The authors declare no competing financial interest. ACKNOWLEDGMENT The work was supported by the Academy of Finland through the Project KUMURA (No. 1277993). The CSC-IT Center for Science is thanked for computational resources.

ACS Paragon Plus Environment

28

Page 29 of 38

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

REFERENCES (1)

Galloway, J. N.; Likens, G. E.; Keene, W. C.; Miller, J. M. The Composition of

Precipitation in Remote Areas of the World. J. Geophys. Res. 1982, 87, 8771-8786. (2)

Chameides, W. L.; Davis D. D. Aqueous-Phase Source of Formic Acid in Clouds. Nature

1983, 304, 427-429. (3)

Ravishankara, A. R.; Daniel, J. S.; Portmann, R. W. Nitrous Oxide (N2O): The Dominant

Ozone-Depleting Substance Emitted in the 21st Century. Science 2009, 326, 123. (4)

Myhre, G.; Shindell, D.; Bréon, F.-M.; Collins, W.; Fuglestvedt, J.; Huang, J.; Koch, D.;

Lamarque, J.-F.; Lee, D.; Mendoza, B.; Nakajima, T.; et al. In Climate Change 2013: The Physical Science Basis. Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change; Stocker, T. F., Qin, D., Plattner, G.-K., Tignor, M., Allen, S. K., Boschung, J., Nauels, A., Xia, Y., Bex, V., Midgley, P. M., Eds.; Cambridge University Press: Cambridge, United Kingdom and New York, NY, USA, 2013; pp 659–740. (5)

Stavrakou, T.; Muller, J. F.; Peeters, J.; Razavi, A.; Clarisse, L.; Clerbaux, C.; Coheur, P.

F.; Hurtmans, D.; De Maziere, M.; Vigouroux, C.; et al. Satellite Evidence for a Large Source of Formic Acid from Boreal and Tropical Forests. Nat. Geosci. 2012, 5, 26-30. (6)

Voigt, C.; Marushchak, M. E.; Lamprecht, R. E.; Jackowicz-Korczyński, M.; Lindgren,

A.; Mastepanov, M.; Granlund, L.; Christensen, T. R.; Tahvanainen, T.; Martikainen, P. J.; et al. Increased Nitrous Oxide Emissions from Arctic Peatlands after Permafrost Thaw. Proc. Natl. Acad. Sci. U.S.A. 2017, 114, 6238-6243. (7)

Zuckerman, B.; Ball, J. A.; Gottlieb, C. A. Microwave Detection of Interstellar Formic

Acid. Astrophys. J. 1971, 163, L41-L45.

ACS Paragon Plus Environment

29

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(8)

Page 30 of 38

Winnewisser, G.; Churchwell, E. Detection of Formic Acid in Sagittarius B2 by Its 211-

212 Transition. Astrophys. J. 1975, 200, L33-L36. (9)

Ziurys, L. M.; Apponi, A. J.; Hollis, J. M.; Snyder, L. E. Detection of Interstellar N2O: A

New Molecule Containing an N−O Bond. Astrophys. J. 1994, 436, L181-L184. (10)

Zolandz, D.; Yaron, D.; Peterson, K. I.; Klemperer, W. Water in Weak Interactions: The

Structure of the Water–Nitrous Oxide Complex. J. Chem. Phys. 1992, 97, 2861-2868. (11)

Lundell, J.; Räsänen, M.; Latajka, Z. Complexes between Formic Acid and Carbon

Monoxide: An Ab Initio Investigation. J. Phys. Chem. 1993, 97, 1152-1157. (12)

Peebles, R. A.; Peebles, S. A.; Kuczkowski, R. L.; Leung, H. O. Isotopic Studies,

Structure and Modeling of the Nitrous Oxide−Acetylene Complex. J. Phys. Chem. A 1999, 103, 10813-10818. (13)

Priem, D.; Ha, T.-K.; Bauder, A. Rotational Spectra and Structures of Three Hydrogen-

Bonded Complexes between Formic Acid and Water. J. Chem. Phys. 2000, 113, 169-175. (14)

Zhou, Z.; Shi, Y.; Zhou, X. Theoretical Studies on the Hydrogen Bonding Interaction of

Complexes of Formic Acid with Water. J. Phys. Chem. A 2004, 108, 813-822. (15)

Panek, J. J.; Wawrzyniak, P. K.; Latajka, Z.; Lundell, J. Interaction Energy

Decomposition Analysis for Formic Acid Complexes with Argon and Krypton Atoms. Chem. Phys. Lett. 2006, 417, 100-104. (16)

Solimannejad, M.; Alkorta, I.; Elguero, J. Weakly Bound Complexes of N2O:  An Ab

Initio Theoretical Analysis toward the Design of N2O Receptors. J. Phys. Chem. A 2007, 111, 2077-2083. (17)

Salmon, S. R.; de Lange, K. M.; Lane, J. R. Structure and Abundance of Nitrous Oxide

Complexes in Earth’s Atmosphere. J. Phys. Chem. A 2016, 120, 2096-2105.

ACS Paragon Plus Environment

30

Page 31 of 38

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

(18)

Murdachaew, G.; Nathanson, G. M.; Gerber, R. B.; Halonen, L. Deprotonation of Formic

Acid in Collisions with a Liquid Water Surface Studied by Molecular Dynamics and Metadynamics Simulations. Phys. Chem. Chem. Phys. 2016, 18, 29756-29770. (19)

Farfán, P.; Echeverri, A.; Diaz, E.; Tapia, J. D.; Gómez, S.; Restrepo, A. Dimers of

Formic Acid: Structures, Stability, and Double Proton Transfer. J. Chem. Phys. 2017, 147, 044312. (20)

Barclay, A. J.; Lauzin, C.; Sheybani-Deloui, S.; Michaelian, K. H.; Moazzen-Ahmadi, N.

Detection of a Higher Energy Isomer of the CO-N2O van der Waals Complex and Determination of Two of Its Intermolecular Frequencies. Phys. Chem. Chem. Phys. 2017, 19, 1610-1613. (21)

Rezaei, M.; Michaelian, K. H.; McKellar, A. R. W.; Moazzen-Ahmadi, N. Infrared

Spectra of the Ne2–N2O, Ar2–N2O Trimers. J. Mol. Spectrosc. 2012, 278, 17-22. (22)

Vohringer-Martinez, E.; Hansmann, B.; Hernandez, H.; Francisco, J. S.; Troe, J.; Abel,

B. Water Catalysis of a Radical-Molecule Gas-Phase Reaction. Science 2007, 315, 497-501. (23)

Tachikawa, H.; Kawabata, H. Effects of a Single Water Molecule on the Reaction Barrier

of Interstellar CO2 Formation Reaction. J. Phys. Chem. A 2016, 120, 6596-6603. (24)

Zhang, L.; Yang, L.; Zhao, Y.; Zhang, J.; Feng, D.; Sun, S. Effects of Water Molecule on

CO Oxidation by OH: Reaction Pathways, Kinetic Barriers, and Rate Constants. J. Phys. Chem. A 2017, 121, 4868-4880. (25)

Meisner, J.; Kästner, J. Atom Tunneling in Chemistry. Angew. Chem. Int. Ed. 2016, 55,

5400-5413. (26)

Atom Tunneling Phenomena in Physics, Chemistry and Biology; Miyazaki, T., Ed.;

Springer Berlin Heidelberg: Berlin, Heidelberg, 2004.

ACS Paragon Plus Environment

31

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(27)

Page 32 of 38

Gong, Y.; Zhou, M.; Andrews, L. Spectroscopic and Theoretical Studies of Transition

Metal Oxides and Dioxygen Complexes. Chem. Rev. 2009, 109, 6765-6808. (28)

Young, N. A. Main Group Coordination Chemistry at Low Temperatures: A Review of

Matrix Isolated Group 12 to Group 18 Complexes. Coord. Chem. Rev. 2013, 257, 956-1010. (29)

Khriachtchev, L. Matrix-Isolation Studies of Noncovalent Interactions: More

Sophisticated Approaches. J. Phys. Chem. A 2015, 119, 2735-2746. (30)

Cao, Q.; Gor, G. Y.; Krogh-Jespersen, K.; Khriachtchev, L. Non-Covalent Interactions of

Nitrous Oxide with Aromatic Compounds: Spectroscopic and Computational Evidence for the Formation of 1:1 Complexes. J. Chem. Phys. 2014, 140, 144304. (31)

George, L.; Sanchez-García, E.; Sander, W. Matrix Isolation Infrared and Ab Initio Study

of Formic Acid−Acetylene Interaction:  Example of H···π and C−H···O Interaction. J. Phys. Chem. A 2003, 107, 6850-6858. (32)

Sánchez-García, E.; George, L.; Montero, L. A.; Sander, W. 1:2 Formic Acid/Acetylene

Complexes:  Ab Initio and Matrix Isolation Studies of Weakly Interacting Systems. J. Phys. Chem. A 2004, 108, 11846-11854. (33)

Tsuge, M.; Marushkevich, K.; Räsänen, M.; Khriachtchev, L. Infrared Characterization

of the HCOOH···CO2 Complexes in Solid Argon: Stabilization of the Higher-Energy Conformer of Formic Acid. J. Phys. Chem. A 2012, 116, 5305-5311. (34)

Marushkevich, K.; Räsänen, M.; Khriachtchev, L. Interaction of Formic Acid with

Nitrogen: Stabilization of the Higher-Energy Conformer. J. Phys. Chem. A 2010, 114, 1058410589.

ACS Paragon Plus Environment

32

Page 33 of 38

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

(35)

Cao, Q.; Melavuori, M.; Lundell, J.; Räsänen, M.; Khriachtchev, L. Matrix-Isolation and

Ab Initio Study of the Complex between Formic Acid and Xenon. J. Mol. Struct. 2012, 1025, 132-139. (36)

Banerjee, P.; Bhattacharya, I.; Chakraborty, T. Matrix Isolation Infrared Spectroscopy of

an O–H···π Hydrogen-Bonded Complex between Formic Acid and Benzene. J. Phys. Chem. A 2016, 120, 3731-3739. (37)

Jin, X.; Wang, G.; Zhou, M. Matrix Isolation Infrared Spectroscopic and Theoretical

Study of Nickel, Palladium, and Platinum Nitrous Oxide Complexes. J. Phys. Chem. A 2006, 110, 8017-8022. (38)

Marushkevich, K.; Khriachtchev, L.; Räsänen, M. Hydrogen Bonding between Formic

Acid and Water:  Complete Stabilization of the Intrinsically Unstable Conformer. J. Phys. Chem. A 2007, 111, 2040-2042. (39)

Zins, E. L.; Krim, L. Photochemistry of a 1 : 1 Hydrogen-Bonded CH3CN : HCOOH

Complex under Astrochemically-Relevant Conditions. Phys. Chem. Chem. Phys. 2014, 16, 33883398. (40)

Sánchez-García, E.; Studentkowski, M.; Montero, L. A.; Sander, W. Noncovalent

Complexes between Dimethyl Ether and Formic Acid—an Ab Initio and Matrix Isolation Study. ChemPhysChem 2005, 6, 618-624. (41)

George, L.; Sander, W. Matrix Isolation Infrared and Ab Initio Study of the Hydrogen

Bonding between Formic Acid and Water. Spectrochim. Acta A Mol. Biomol. Spectrosc. 2004, 60, 3225-3232.

ACS Paragon Plus Environment

33

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(42)

Page 34 of 38

Rozenberg, M.; Loewenschuss, A.; Nielsen, C. J. H-Bonding of Formic Acid with Its

Decomposition Products: A Matrix Isolation and Computational Study of the HCOOH/CO and HCOOH/CO2 Complexes. J. Phys. Chem. A 2015, 119, 8497-8502. (43)

Marushkevich, K.; Siltanen, M.; Räsänen, M.; Halonen, L.; Khriachtchev, L.

Identification of New Dimers of Formic Acid: The Use of a Continuous-Wave Optical Parametric Oscillator in Matrix Isolation Experiments. J. Phys. Chem. Lett. 2011, 2, 695-699. (44)

Sánchez-García, E.; Mardyukov, A.; Studentkowski, M.; Montero, L. A.; Sander, W.

Furan−Formic Acid Dimers:  An Ab Initio and Matrix Isolation Study. J. Phys. Chem. A 2006, 110, 13775-13785. (45)

Bahou, M.; Schriver-Mazzuoli, L.; Camy-Peyret, C.; Schriver, A.; Chiavassa, T.; Aycard,

J. P. Infrared Matrix Spectra of the N2O···O2 Complex in Solid Nitrogen. The N2O···O2+O Thermal Diffusion Limited Reaction. Chem. Phys. Lett. 1997, 265, 145-153. (46)

Pettersson, M.; Lundell, J.; Khriachtchev, L.; Räsänen, M. IR Spectrum of the Other

Rotamer of Formic Acid, Cis-HCOOH. J. Am. Chem. Soc. 1997, 119, 11715-11716. (47)

Pettersson, M.; Maçôas, E. M. S.; Khriachtchev, L.; Fausto, R.; Räsänen, M.

Conformational Isomerization of Formic Acid by Vibrational Excitation at Energies Below the Torsional Barrier. J. Am. Chem. Soc. 2003, 125, 4058-4059. (48)

Maçôas, E. M. S.; Lundell, J.; Pettersson, M.; Khriachtchev, L.; Fausto, R.; Räsänen, M.

Vibrational Spectroscopy of Cis- and Trans-Formic Acid in Solid Argon. J. Mol. Spectrosc. 2003, 219, 70-80. (49)

Maçôas, E. M. S.; Khriachtchev, L.; Pettersson, M.; Fausto, R.; Räsänen, M. Rotational

Isomerization of Small Carboxylic Acids Isolated in Argon Matrices: Tunnelling and Quantum Yields for the Photoinduced Processes. Phys. Chem. Chem. Phys. 2005, 7, 743-749.

ACS Paragon Plus Environment

34

Page 35 of 38

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

(50)

Khriachtchev, L. Rotational Isomers of Small Molecules in Noble-Gas Solids: From

Monomers to Hydrogen-Bonded Complexes. J. Mol. Struct. 2008, 880, 14-22. (51)

Pettersson, M.; Maçôas, E. M. S.; Khriachtchev, L.; Lundell, J.; Fausto, R.; Räsänen, M.

Cis→Trans Conversion of Formic Acid by Dissipative Tunneling in Solid Rare Gases: Influence of Environment on the Tunneling Rate. J. Chem. Phys. 2002, 117, 9095-9098. (52)

Marushkevich, K.; Khriachtchev, L.; Lundell, J.; Räsänen, M. Cis−Trans Formic Acid

Dimer:  Experimental Observation and Improved Stability against Proton Tunneling. J. Am. Chem. Soc. 2006, 128, 12060-12061. (53)

Marushkevich, K.; Khriachtchev, L.; Räsänen, M.; Melavuori, M.; Lundell, J. Dimers of

the Higher-Energy Conformer of Formic Acid: Experimental Observation. J. Phys. Chem. A 2012, 116, 2101-2108. (54)

Marushkevich, K.; Khriachtchev, L.; Lundell, J.; Domanskaya, A.; Räsänen, M. Matrix

Isolation and Ab Initio Study of Trans−Trans and Trans−Cis Dimers of Formic Acid. J. Phys. Chem. A 2010, 114, 3495-3502. (55)

Khriachtchev, L.; Domanskaya, A.; Marushkevich, K.; Räsänen, M.; Grigorenko, B.;

Ermilov, A.; Andrijchenko, N.; Nemukhin, A. Conformation-Dependent Chemical Reaction of Formic Acid with an Oxygen Atom. J. Phys. Chem. A 2009, 113, 8143-8146. (56)

Zhao, Y.; Truhlar, D. G. The M06 Suite of Density Functionals for Main Group

Thermochemistry, Thermochemical Kinetics, Noncovalent Interactions, Excited States, and Transition Elements: Two New Functionals and Systematic Testing of Four M06-Class Functionals and 12 Other Functionals. Theor. Chem. Acc. 2008, 120, 215-241. (57)

Zhao, Y.; Truhlar, D. G. Density Functionals with Broad Applicability in Chemistry. Acc.

Chem. Res. 2008, 41, 157-167.

ACS Paragon Plus Environment

35

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(58)

Page 36 of 38

Chai, J.-D.; Head-Gordon, M. Long-Range Corrected Hybrid Density Functionals with

Damped Atom-Atom Dispersion Corrections. Phys. Chem. Chem. Phys. 2008, 10, 6615-6620. (59)

Møller, C.; Plesset, M. S. Note on an Approximation Treatment for Many-Electron

Systems. Phys. Rev. 1934, 46, 618-622. (60)

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 16, Revison A.03; Gaussian, Inc.: Wallingford CT, 2016. (61)

Dunning, T. H. Gaussian Basis Sets for Use in Correlated Molecular Calculations. I. The

Atoms Boron through Neon and Hydrogen. J. Chem. Phys. 1989, 90, 1007-1023. (62)

Werner, H. J.; Knowles, P. J.; Knizia, G.; Manby, F. R.; Schütz, M. Molpro: A General‐

Purpose Quantum Chemistry Program Package. WIREs Comput Mol Sci 2012, 2, 242-253. (63)

Werner, H. J.; Knowles, P. J.; Knizia, G.; Manby, F. R.; Schütz, M.; Celani, P.; Györffy,

W.; Kats, D.; Korona, T.; Lindh, R.; et al. Molpro, Version 2015.1, a Package of Ab Initio Programs, 2015. (64)

Peterson, K. A.; Adler, T. B.; Werner, H.-J. Systematically Convergent Basis Sets for

Explicitly Correlated Wavefunctions: The Atoms H, He, B–Ne, and Al–Ar. J. Chem. Phys. 2008, 128, 084102. (65)

Boys, S. F.; Bernardi, F. The Calculation of Small Molecular Interactions by the

Differences of Separate Total Energies. Some Procedures with Reduced Errors. Mol. Phys. 1970, 19, 553-566. (66)

Trakhtenberg, L. I. In Atom Tunneling Phenomena in Physics, Chemistry and Biology;

Miyazaki, T., Ed.; Springer Berlin Heidelberg: Berlin, Heidelberg, 2004; pp 33-58.

ACS Paragon Plus Environment

36

Page 37 of 38

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

(67)

Tsuge, M.; Khriachtchev, L. Tunneling Isomerization of Small Carboxylic Acids and

Their Complexes in Solid Matrixes: A Computational Insight. J. Phys. Chem. A 2015, 119, 26282635. (68)

Goldanskii, V. I. Chemical Reactions at Very Low Temperatures. Annu. Rev. Phys.

Chem. 1976, 27, 85-126. (69)

Lundell, J.; Jolkkonen, S.; Khriachtchev, L.; Pettersson, M.; Räsänen, M. Matrix

Isolation and Ab Initio Study of the Hydrogen-Bonded H2O2–CO Complex. Chem. Eur. J. 2001, 7, 1670-1678. (70)

Khriachtchev, L.; Tapio, S.; Räsänen, M.; Domanskaya, A.; Lignell, A. HY⋯N2 and

HXeY⋯N2 Complexes in Solid Xenon (Y = Cl and Br): Unexpected Suppression of the Complex Formation for Deposition at Higher Temperature. J. Chem. Phys. 2010, 133, 084309. (71)

Shimanouchi, T. Tables of Molecular Vibrational Frequencies Consolidated Volume I;

National Bureau of Standards, 1972. (72)

Domanskaya, A.; Marushkevich, K.; Khriachtchev, L.; Räsänen, M. Spectroscopic Study

of Cis-to-Trans Tunneling Reaction of HCOOD in Rare Gas Matrices. J. Chem. Phys. 2009, 130, 154509. (73)

Ryazantsev, S. V.; Feldman, V. I.; Khriachtchev, L. Conformational Switching of HOCO

Radical: Selective Vibrational Excitation and Hydrogen-Atom Tunneling. J. Am. Chem. Soc. 2017, 139, 9551-9557. (74)

Potapov, A. Weakly Bound Molecular Complexes in the Laboratory and in the

Interstellar Medium: A Lost Interest? Mol. Astrophys. 2017, 6, 16-21.

ACS Paragon Plus Environment

37

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 38 of 38

TOC Graphic

ACS Paragon Plus Environment

38