Katharina A. E.
Meyer
* and
Arman
Nejad
Institute of Physical Chemistry, University of Göttingen, Tammannstr. 6, 37077 Göttingen, Germany. E-mail: katharina.meyer@chemie.uni-goettingen.de; anejad@gwdg.de
First published on 5th August 2021
The cis–trans-isomerism of the propiolic acid monomer (HCC–COOH) is examined with linear Raman jet spectroscopy, yielding the first environment-free vibrational band centres of a higher-energy cis-rotamer beyond formic acid (HCOOH) in addition to all fundamentals and a large number of hot and combination/overtone bands of the trans-conformer. Two near-isoenergetic trans-fundamentals of different symmetry (CC
O bend and OH torsion) prove to be a sensitive benchmarking target, as their energetic order is susceptible to the choice of electronic structure method, basis set size, and inclusion of vibrational anharmonicity. For the infrared- and Raman-active C
O stretching fundamentals of the cyclic (C2h) trans-propiolic acid dimer, resonance couplings are found that in part extend to the Cs-symmetric heterodimer of trans-propiolic and trans-formic acid. Exploratory vibrational perturbation theory (VPT2) calculations show that all perturbing states involve displacements of the OH moieties located on the doubly hydrogen bonded ring. The comparison of the infrared spectra of the propiolic acid dimer and its heterodimer with formic acid to that of several other carboxylic acid dimers from the literature reveals a notable similarity regarding a non-fundamental dimer band around 1800 cm−1, which in most cases is so far unassigned. VPT2 calculations and a simple harmonic model suggest an assignment to a combination vibration of the symmetric and antisymmetric OH torsion.
![]() | ||
Fig. 1 The cis- and trans-rotamers of carboxylic acid monomers as well as the global minimum structure of the dimer. |
In contrast to the chaotic, resonance broadened O–H stretching spectrum spanning several hundreds of cm−1 that many carboxylic acid dimers seem to share,7,9,19,20 the vibrational spectrum of the soft modes up to the carbonyl stretching vibration around 1700–1800 cm−1 is less convoluted and more suitable for a systematic vibrational characterisation. Mediated by the hydrogen bonds, these cyclic dimers are centrosymmetric (C2h) and pairs of identical monomer vibrations are split into so-called Davydov pairs,21–23 with the symmetric, in-phase combinations of monomer vibrations being exclusively Raman-active and the antisymmetric, out-of-phase combinations being infrared-active. The systematic comparison of the infrared and Raman spectra of several different acids reveals interesting phenomena which to the best of our knowledge have so far not been systematically analysed. In the infrared spectra, for example, a non-fundamental dimer band is observed at a similar band position for formic (1810.6 cm−1),24 acetic (1834.9 cm−1),25 pivalic (1834/1838 cm−1),26 and benzoic acid (1824 cm−1)27,28 which raises the question of a possible similar origin. For formic acid, it was ascribed to a combination vibration of the two (hindered) OH torsions by Kollipost,24 but for the other acids, it has yet to be assigned. The Raman spectra of several carboxylic acids also exhibit a similarity in this spectral region, namely strong resonance perturbations of the symmetric CO stretches (cf.ref. 29 and 30). While a combination band was predicted in proximity to the symmetric C
O stretch of the formic acid dimer nearly forty years ago,8 conclusive evidence for a potential resonance partner has yet to be provided.
Slightly larger aliphatic, saturated carboxylic acids such as acetic (CH3COOH) or pivalic acid ((CH3)3COOH), in addition to the challenge constituted by their sheer size, significantly complicate routine anharmonic treatments31 due to a manifold of internal rotations introduced by the methyl groups.32 Propiolic acid (HCC–COOH), as the smallest unsaturated, aliphatic carboxylic acid, constitutes with only seven atoms and fifteen relatively stiff internal degrees of freedom an excellent intermediate between the well studied formic and the larger saturated carboxylic acids. The vibrational spectrum of its dimer has been studied in solid rare gas matrices,33,34 and in the gas phase using low-resolution infrared19 and jet-cooled Raman spectroscopy.29,35 Building on a previous study,29 we will present and scrutinise the vibrational spectrum of the propiolic acid homodimer in the C
O stretching region by means of FTIR and Raman jet spectroscopy and characterise the vibrational spectrum of the formic–propiolic acid heterodimer for the first time, which has been thoroughly investigated with microwave spectroscopy.36–39 Our results will show how certain features in the C
O stretching spectra of the formic and propiolic acid homo- and heterodimers can qualitatively be understood as perturbations due to vibrational motion of the cyclic ring shared by all three dimers.
One very intriguing feature of the monomers of carboxylic acids is their torsional isomerism that connects the higher-energy cis- and global minimum trans-rotamers (Fig. 1). The first vibrational detection of a cis-carboxylic acid was achieved for formic acid in a cryogenic argon matrix by Pettersson et al. in 1997,40 nearly sixty years after the first vibrational characterisation of its trans-form.41 Ever since, these matrix isolation studies have been extended to a multitude of other acids34,42–45 and matrix hosts.44,46–48 One advantage of this technique are the long storing times achievable, which enable the formation of higher-energy conformers via light-induced isomerisation.49 In the gas phase, where a comparison with theory is more easily achievable due to the absence of external perturbations, vibrational band positions of the higher-energy cis-carboxylic acids are sparse – a consequence of their low abundance at room temperature (RT). The only higher-energy cis-carboxylic acid where perturbation-free reference data are available to this date is formic acid (0.1% cis-abundance at RT)50 with a total of sixteen fundamentals for all four H/D isotopologues.30,51–53
The propiolic acid monomer has a significantly lower cis–trans-energy difference compared to formic acid, resulting in a room temperature abundance that is an order of magnitude higher (about 1%), yet vibrational data on the cis-rotamer are so far restricted to matrix isolation studies.34,44 The global minimum trans-rotamer has been examined with low-resolution gas phase infrared spectroscopy,54 in the liquid phase with infrared and Raman spectroscopy,54 with matrix isolation spectroscopy,33,34,44 and Raman jet spectroscopy29,35 focussing mostly on fundamentals. Following the approach that allowed for the detection of cis-formic acid using a heatable Raman jet set-up,30,52,53 we will present the first perturbation-free vibrational data for cis-propiolic acid in addition to a manifold of vibrational information beyond fundamentals on trans-propiolic acid, highlighting that the torsional isomerism of the propiolic acid monomer as well as its homo- and heterodimer with formic acid can prove to be meaningful reference systems for theory.
For the Raman jet expansions, propiolic acid (Acros Organics, 98%) was seeded in helium in a temperature controlled saturator (temperature range 15–20 °C) and expanded continuously at reservoir pressures of 400–600 mbar through a temperature-controlled slit nozzle (4 × 0.15 mm2) and feed-line (RT – 190 °C) into an evacuated jet-chamber with background pressures below 1.5 mbar during the expansion. For measurements of the formic–propiolic acid heterodimer, formic acid (Acros Organics, 98+%) was also seeded in helium (saturator temperature of −30 °C) and both acid-in-helium mixtures were premixed in the reservoir prior to the expansion (for further details, see Section S1 in the ESI†). The acid-in-helium expansion is probed at a nozzle distance of 1 mm with a Spectra Physics Millennia eV 532 nm laser operated at 20–24 W. The scattered radiation is collected perpendicular to both the expanding flow and the laser, and is focussed onto the entry slit of a 1 m monochromator, which disperses the radiation onto 1340 pixel columns of a liquid nitrogen cooled CCD camera. Vertical binning is performed along 400 pixel rows. For each spectrum, 6–24 spectra were co-added with exposure times of 300 s. Spectral calibration is carried out with atomic neon and krypton transitions in vacuum.59 The experimental uncertainty of the band positions resulting from the spectral calibration and resolving power of the monochromator depends on the spectral range and amounts to about 1–2 cm−1. We therefore assign Raman band position errors of ±2 cm−1 throughout this work.
The polarisation of the incident laser radiation of the Raman set-up is perpendicular ⊥ with respect to the 90° scattering plane.55 In some cases, additional measurements with parallel laser polarisation ‖ were carried out to compare these to the ⊥-spectra, which can be used to distinguish totally symmetric from non-totally symmetric modes.60 Throughout this work, the higher detection sensitivity of perpendicularly compared to parallelly polarised laser light of the set-up is neglected.61
Infrared spectra were recorded with a Bruker IFS 66v spectrometer equipped with a globar, a potassium bromide (KBr) beam splitter as well as KBr optics. The IR beam is focussed onto a pulsed supersonic expansion (reservoir pressures of 500–700 mbar) from a 600 mm long and 0.2 mm wide slit nozzle. In contrast to the Raman set-up, the nozzle is not heatable. The distance of IR beam and nozzle amounts to 10 mm in the focal point. For detection of the infrared radiation, a liquid nitrogen cooled mercury cadmium telluride detector is used.
To assist in the interpretation of the recorded vibrational spectra, harmonic Raman scattering cross-sections and anharmonic (VPT2,62–64 vibrational perturbation theory of second-order) vibrational energies have been computed using Gaussian 16 version A.03.65 To solve the electronic Schrödinger equation, several ab initio and DFT methods including MP2,66 B3LYP,67,68 PBE0,69 and B2PLYP70 were used. For the latter three, Grimme's two-body dispersion correction71 with Becke–Johnson damping has been employed.72
In the VPT2 calculations, default thresholds for resonance detection were used, but Darling–Dennison resonances were deactivated as low-frequency fundamental energies seem to be very sensitive to the resonance thresholds as implemented in Gaussian 16 (version A.03).56 Additional calculations upon inclusion of Darling–Dennison resonances are shown in Fig. S1 in the ESI.† If not further specified, symmetry and very tight optimisation thresholds are used, and in case of DFT calculations the integration grid is set to superfine. To enhance the quality and remove some deficiencies of the harmonic MP2 force field,73 the Hessian is replaced by the corresponding harmonic CCSD(T)-F12a74 force field which has been computed using Molpro version 2018.1.75,76 For the DFT and MP2 calculations, Dunning's correlation consistent aug-cc-pVXZ (aVXZ)77 and for CCSD(T)-F12a, explicitly correlated VXZ-F12 basis sets by Peterson et al.78 are used. A comprehensive summary of the computational details and used keywords can be found in the ESI† (Tables S1 and S2).
![]() | ||
Fig. 2 Schematic representation of the fifteen normal modes of trans-propiolic acid alongside their VPT2 band positions (in cm−1, see Table 1). Modes are labelled according to the Herzberg nomenclature. ν1–ν11 are totally symmetric in-plane (A′) and ν12–ν15 non-totally symmetric out-of-plane vibrations (A′′). |
Fundamental | trans-Propiolic acid | cis-Propiolic acid | ||||||||
---|---|---|---|---|---|---|---|---|---|---|
VPT2a | σ Ra | Raman jetb | FTIR jetb | N2 matrixc | VPT2a | σ Ra | Raman jetb | N2 matrixc | ||
a Calculated at the MP2/aVQZ level where the Hessian is replaced with the corresponding CCSD(T)-F12a/VTZ-F12 force field. b This work. FTIR spectra are only recorded between 2000–700 cm−1. c Infrared matrix isolation data from ref. 34 (matrix site 1). d Strong Fermi resonance predicted between ν5 and ν12 + ν14 of trans-propiolic acid. Note that while the assignment inverts when changing the anharmonic MP2 basis set from aVQZ to aVTZ, both calculations predict higher infrared activity for the lower-wavenumber band (cf. Table S4 in the ESI). e Band maxima are listed, see band contours in Fig. S4 in the ESI. The band centre likely corresponds to the dip between both peaks. | ||||||||||
A′ | ||||||||||
ν 1 | ν(OH) | 3578 | 99 | 3578 | 3534.5 | 3622 | 59 | 3621 | 3575.5 | |
ν 2 | ν(CH) | 3330 | 47 | 3329 | 3404.5 | 3324 | 47 | 3298.0 | ||
ν 3 |
ν(C![]() |
2160 | 340 | 2143 | 2147.0, 2142.5 | 2136 | 324 | 2121.0 | ||
ν 4 |
ν(C![]() |
1772 | 114 | 1764 | 1765 | 1750.0, 1748.0 | 1805 | 115 | 1800 | 1781.0, 1776.5 |
ν 5 | δ(OH) | 1361 | 8 | 1362 | 1363 | 1326.5 | 1289 | 15 | 1248.0 | |
1299 | 1302 | 1303 | ||||||||
ν 6 | ν(C–O) | 1154 | 15 | 1154 | 1154 | 1152.5 | 1152 | 41 | n.o. | |
ν 7 | ν(C–C) | 814 | 91 | 816 | 816 | 822.0 | 820 | 61 | 832.0 | |
ν 8 | δ(CH) | 652 | 13 | 654 | 665.5 | 667 | 14 | n.o. | ||
ν 9 |
δ(CC![]() |
585 | 10 | 584 | 586.7 | 584 | 9 | 591.0 | ||
ν 10 | δ(CC–O) | 527 | 19 | 527 | 534.5 | 537 | 22 | n.o. | ||
ν 11 | δ(CCC) | 178 | 191 | 180 | n.i. | 179 | 216 | n.i. | ||
A′′ | ||||||||||
ν 12 | γ(OCO) | 754 | 1 | n.o. | 755 | 760.0 | 741 | 3 | 749.5 | |
ν 13 | γ(CH) | 682 | 17 | 696/688e | 760.0 | 689 | 18 | 703.5 | ||
ν 14 | γ(OH) | 574 | 21 | 578/569e | 606.0 | 476 | 8 | 508.0 | ||
ν 15 | γ(CCC) | 229 | 49 | 233/225e | n.i. | 230 | 67 | n.i. |
![]() | ||
Fig. 3 Raman jet spectra of propiolic acid in helium recorded at nozzle temperatures between 100–190 °C. The Raman spectra are intensity-scaled to a trans-fundamental (tPr) in the spectrum recorded at 190 °C (red) in each spectral region. cis-Propiolic acid monomer bands are labelled ‘cPr’ and bands corresponding to propiolic acid decomposition products are marked with an asterisk. Additional measurements with parallel laser polarisation (nozzle temperature 100 °C) are shown in green and have been scaled by a factor of 7/6 (see text and ref. 60 for details). Below the Raman spectra, anharmonic VPT2 fundamentals of cPr and tPr and two-quantum states (tPr) calculated at the MP2/aV(T,Q)Z level are shown using harmonic Raman scattering cross-sections. For the VPT2 calculations, the Hessian is replaced with the corresponding CCSD(T)-F12a/VTZ-F12 force field. The symmetry of the vibrations (A′ or A′′) is encoded in the colour shade of the symbols (see legend in panel (B)). |
Raman jet spectra of the six stretching vibrations (ν1–ν4 and ν6–ν7) and other selected spectral regions promising for cis-detection are shown in Fig. 3 for nozzle temperatures between 100–190 °C. The Raman jet spectra of the A′′ fundamentals are shown in the ESI† (Fig. S4). By scaling the spectra recorded at different nozzle temperatures to a known trans-propiolic acid (tPr) fundamental in each spectral region, cluster and hot bands can be distinguished from fundamentals, combination and overtone bands of tPr, as cluster bands decrease in intensity with temperature whereas hot bands increase.52 Two types of hot bands are visible in the spectra – isomeric hot bands, i.e., those corresponding to the higher-energy cis-rotamer of propiolic acid (cPr), and non-isomeric hot bands originating from thermally populated vibrational states νj in the trans-rotamer. In the context of perturbation theory,62 the shift of a non-isomeric hot band νi + νj − νj with respect to the fundamental νi amounts to the anharmonic coupling constant xi,j (2xi,i for 2νi − νi). Assuming that the scattering cross-sections of the fundamental νi and the hot band νi + νj − νj are similar, the relative intensity of νi + νj − νj with respect to the fundamental reflects the thermal population in νj.52 Hence, a comparison of band position shifts and intensities with respect to the fundamental in each spectral region can help to distinguish isomeric from non-isomeric hot bands. Below the Raman jet spectra in Fig. 3, CCSD(T)-edited MP2 VPT2 energy levels are plotted, where the MP2 Hessian is replaced with the corresponding CCSD(T)-F12a/VTZ-F12 force field. To visualise convergence of the anharmonic MP2 force field, results using an aVTZ (dashed) and aVQZ basis set (solid lines) are shown. For brevity, we will refer to the CCSD(T)-edited MP2/aVQZ VPT2 calculation as ‘the VPT2’ calculation throughout this section.
The O–H stretching spectrum of propiolic acid is shown in Fig. 3A. It exhibits one fairly strong band of a combination or overtone vibration at 3571(2) cm−1, 7 cm−1 downshifted from the trans-fundamental (3578(2) cm−1). The cis-band is observed upshifted from both at 3621(2) cm−1. Due to the large cis–trans-shift of +43 cm−1 together with a moderate Raman scattering cross-section (cf.Table 1), ν1 of cis-propiolic acid is one of the most evident cis-assignments. The corresponding formic acid cis–trans-shift is with +67 cm−1 about a third larger.52 Close to the trans-propiolic acid band at 3571(2) cm−1, a combination of the C–H stretching vibration ν2 with the CCC out-of-plane bending vibration ν15 is predicted by the VPT2 calculation (3560 cm−1, A′′ symmetry). However, in line with the high intensity likely caused by resonance interaction with the fundamental, depolarisation measurements reveal A′ symmetry for the band (Fig. S3 in the ESI†). The energetically closest A′-symmetric two-quantum state is predicted at 3528 cm−1 (2ν4). Thus, an assignment of the band at 3571(2) cm−1 based on the VPT2 calculation is not possible. In the Raman jet spectrum of the formic acid OH stretching vibration, two combination bands were observed downshifted from the trans-OH fundamental. One was attributed to a combination vibration of the C–H stretch with the OCO in-plane bend,79 a mode that is distinctly higher in energy in case of propiolic acid due to the stiffer CH stretch. The second combination band was attributed to a state composed of more than two vibrational quanta,79 which might also be the case for the trans-propiolic acid band at 3571(2) cm−1. In that context, it is worth pointing to the numerical similarities between the formic and propiolic acid vibrational energy levels with regard to vibrations of the carboxylic –COOH moiety that indicate a possible similar origin for some perturbers in the OH stretching spectra (see note in ref. 80 for further details).
A similarly straightforward cis-assignment as for the OH-stretch of propiolic acid is also feasible for the CO stretch ν4, which is shown in Fig. 3D. The cis-band is observed at 1800(2) cm−1, 36 cm−1 upshifted from the trans-fundamental at 1764(2) cm−1. Expectedly, the cis–trans-shift resembles that of formic acid (+41 cm−1).52 For the other stretching fundamentals, distinctly smaller cis–trans-shifts are predicted (Table 1), therefore significantly complicating the allocation of isomeric cis-bands.
One such example is the C–H stretch (ν2) with a predicted cis–trans-shift of solely −6 cm−1. For formic acid, a significantly larger shift of −69 cm−1 was observed52 – likely a result of the closer proximity of the OH and CH groups.81 In the propiolic acid C–H stretching spectrum shown in Fig. 3B, a hot band is observed on the lower-wavenumber shoulder of the trans-fundamental (3329(2) cm−1). Since non-isomeric hot bands originating from thermally populated low-lying vibrational energy levels are also predicted to overlap with ν2 of cPr (for xi,j values, see Table S5 in the ESI†), an unambiguous cis-assignment is currently not possible. In addition to this shoulder of the ν2-fundamental, there is another hot band at 3310(2) cm−1, downshifted by 19 cm−1 from ν2 (tPr). The largest x2,j values of the VPT2 calculation for energy levels with significant thermal population are found for the CH in-plane ν8 (−14.7 cm−1) and out-of-plane bending vibration ν13 (−13.6 cm−1). At the B3LYP-D3(BJ)/aVQZ level, which was shown to perform well for formic acid monomer hot band assignments,52 slightly larger values of −19.3 (x2,8) and −19.0 cm−1 (x2,13) are obtained, which are compatible with the observed shift. This interpretation of a relatively strong coupling (large xi,j value) between the CH stretch and corresponding bends matches the findings for formic acid, where a similar experimental shift of −17 cm−1 is observed for the ν2 hot band originating from the C–H out-of-plane bending vibration.79 Therefore, it seems likely that the propiolic acid hot band at 3310(2) cm−1 originates from either one or both of these states.
In the same spectral region (Fig. 3B), three additional bands are observed that exhibit a temperature scaling that differs from that of other (isomeric and non-isomeric) hot bands. In the spectra recorded below 160 °C, these are barely visible but increase significantly in intensity at the highest two nozzle temperatures. This indicates that these bands are due to thermal decomposition of the acid. Due to the steeper intensity-increase of such decomposition bands with temperature, they can be clearly separated from hot bands and therefore do not impair cis-assignments. The most pronounced decomposition band in Fig. 3B at 3372(2) cm−1 corresponds to the CH stretching vibration of acetylene (3372.845(1)82); the other two decomposition bands correspond to its hot bands.82 The observation of bands of the second decomposition product CO2 at 1408(2), 1387(2), 1284(2), and 1264(2) cm−1 (Fig. 3G) further confirms this interpretation and indicates that significant decarboxylation of the acid takes place at temperatures above 130 °C. A full list of all observed acetylene and CO2 decomposition bands and their assignments can be found in the ESI† (Table S7). Using computed harmonic Raman scattering cross-sections for propiolic acid and acetylene at the MP2/aVQZ level, we can estimate the amount of decarboxylation to about 14% at 190 °C. Under similar conditions and the same temperature range, no decomposition bands are observed in the Raman jet spectrum of formic acid (see ref. 79).
The CC stretching vibration ν3 of propiolic acid is by far the most Raman-active mode of all fifteen fundamentals and the predicted cis–trans-shift is also fairly large (−24 cm−1, cf.Table 1). Therefore one might expect it to be one of the most suitable vibrations for cis-detection. The trans-C
C stretch band shown in Fig. 3C (2143(2) cm−1), however, is framed by cluster bands on the higher-wavenumber and an extensive hot band structure on the lower-wavenumber side that spans approximately 25 cm−1. In contrast to most other VPT2 fundamentals, which agree with the Raman jet values within the experimental uncertainty (±2 cm−1), the trans-ν3-band position deviation is particularly large (+17 cm−1). This might be a result of an insufficient description of the electron correlation of the triple bond. As such, the cis-ν3 band position error might also be significant. In addition, there are a vast number of x3,j values between 0 and −9 cm−1 (Table S5 in the ESI†), further complicating an allocation of the cis-propiolic acid ν3-band. Therefore, an assignment is currently not possible. This richness of non-isomeric hot bands afforded by thermal excitation might prove to be very useful in cross-checking assignments of combination and overtone bands, as they are direct probes of anharmonicity.79 To fully profit from the large abundance of spectral information contained in these C
C-stretch Raman jet spectra, further guidance from theory would be particularly useful.
In addition to the broad hot band structure seen in Fig. 3C, there are three A′ bands at 2177(2), 2164(2), and 2110(2) cm−1 that correspond to cold monomer transitions of tPr, i.e., combination or overtone bands. Between 2250 and 2050 cm−1, the VPT2 calculation predicts only one totally symmetric two-quantum state at 2131 cm−1 (ν5 + ν7). It therefore seems likely that some of these trans-monomer bands correspond to energy levels involving more than two vibrational quanta.
For the C–O stretch ν6, an even smaller cis–trans-shift (−2 cm−1) than for the C–H stretch (−6 cm−1) is predicted. Next to the predicted cis-band position, a weak hot band is observed in the spectra (Fig. 3E), which exhibits an intensity that roughly matches the expected cis-intensity with respect to the trans-fundamental at 1154(2) cm−1, considering the Raman scattering cross-sections advantage of the cis- compared to the trans-rotamer (σcPr/σtPr = 3). As the spectral intensity of the band is very low and there are non-isomeric hot bands with similar predicted shifts x6,j (Table S5 in the ESI†) and thermal populations, a cis-assignment is not possible. In addition to the hot band, there is another band downshifted from the trans-fundamental at 1118(2) cm−1 that corresponds to a trans-combination or overtone band with a relative intensity of 25% with respect to ν6 of tPr. The depolarised spectrum confirms that it is of A′ symmetry (Fig. S3 in the ESI†). While the VPT2 calculation predicts two A′ states, 2ν14 and ν9 + ν10, in close proximity to the band at 1118(2) cm−1, and 2ν14 is even predicted to interact with ν6via Fermi resonance, no clear assignment can be made, as Fig. 3E shows that the absolute energies of 2ν14 and ν9 + ν10 are very sensitive to the basis set size of the anharmonic MP2 force field, both seemingly converging towards the band at 1118(2) cm−1.
A substantial basis set sensitivity is also observed for vibrational states in the C–C stretching region (cf.Fig. 3F), yet in contrast to the C–O stretching region, the assignment of anharmonic features in the spectra seems more evident. The trans-fundamental ν7 is observed at 816(2) cm−1 in the Raman jet spectrum and the cis-band is predicted to be upshifted by 6 cm−1 from it. The prediction coincides with a hot band observed at 822(2) cm−1. However, there are a few non-isomeric hot bands for which a positive spectral shift of up to +4 cm−1 with respect to the fundamental is predicted (cf. Table S5 in the ESI†). While a conclusive assignment of ν7 of cis-propiolic acid is not possible, the band at 822(2) cm−1 appears to be the most likely candidate considering the good agreement of the other two cis–trans-shifts. Near the trans-fundamental, there are three additional tPr bands at 832(2), 800(2), and 762(2) cm−1 that all correspond to A′-symmetric transitions (Fig. S3 in the ESI†). The best energy matches for the former two are ν8 + ν11 and ν14 + ν15, respectively. The third band at 762(2) cm−1 is close to the prediction of the OCO out-of-plane bending vibration ν12 (A′′ symmetry) and was assigned to it in a previous Raman jet study by Xue,29 yet ν12 can be ruled out via symmetry and its near-zero Raman scattering cross-section (Table 1). Thus, a more plausible assignment for this band is ν9 + ν11. As ν12 is predicted to have significant infrared activity (Table S4 in the ESI†) but was not reported in the previous gas phase infrared study by Katon and McDevitt,54 we have measured additional FTIR jet spectra between 2000–700 cm−1 (Fig. S2 in the ESI†), where we observe ν12 of trans-propiolic acid at 755 cm−1.
In addition to the stretching vibrations, the VPT2 calculation predicts some A′-symmetric bending vibrations between 1450–500 cm−1 such as the OH (ν5) and the CC–O bend (ν10) to be promising candidates for cis-detection. The other three A′ fundamentals either have too small cis–trans shifts (ν9 and ν11) and/or are predicted in spectrally congested regions in the vicinity of broad A′′ fundamentals of trans-propiolic acid (ν8 and ν9, cf. Fig. S4, ESI†).
The Raman jet spectrum between 1450–1250 cm−1 is shown in Fig. 3G. In this spectral region, only one propiolic acid fundamental is predicted, namely ν5. In the peripheries around 1440 and 1270 cm−1 are contributions from clusters that decrease in intensity with temperature and CO2 bands that increase due to thermal decomposition of the acid. Overall, there are only a few propiolic acid monomer contributions in this spectral region – the band at 1298(2) cm−1, a weak hot band at 1292(2) cm−1, and two broad features centred at 1302(2) and 1362(2) cm−1. ν5 of trans-propiolic acid is predicted to be near-completely depolarised (ρ = 0.73, cf. Table S4 in the ESI†) and should therefore persist in the spectrum measured with parallel laser polarisation (green in Fig. 3G, multiplied by a factor of 7/660). In the parallel spectrum, the band at 1298(2) cm−1 almost completely vanishes, whereas the broad features at 1362(2) and 1302(2) cm−1 largely persist. The former must therefore correspond to a A′-symmetric combination or overtone band, whereas the fundamental can be assigned to one of the latter two bands. In line with the two potential ν5 features observed, the VPT2 calculation predicts a near 1:
1 Fermi resonance between ν5 (OH in-plane bend) and the ν12 + ν14 (OCO out-of-plane bend and OH torsion) combination band at 1361 and 1299 cm−1, in excellent agreement with experiment. As ν5 is predicted to have significant infrared activity (cf. Table S4 in the ESI†), we have measured additional infrared spectra which corroborate the assignment (cf. inset in Fig. 3G). An assignment of the non-fundamental A′ band at 1298(2) cm−1 is beyond the reach of this study, as the density of states is quite large in this spectral region. Based on the good agreement of the VPT2 calculation with experiment, one might be inclined to assign the hot band at 1292(2) cm−1 in the Raman spectrum to ν5 of cPr. However, the overall weak intensity and possible hot band contributions from tPr (x5,j values ranging from 2 to −17 cm−1) prevent any firm assignment. We note that the predicted infrared intensity of ν5 of cPr is the highest amongst all cPr fundamentals (cf. Table S4 in the ESI†), so it might be worthwhile to measure the (jet cooled) infrared spectrum at lower concentrations to prevent overlap with cluster contributions (C in the inset of Fig. 3G).
In search of cPr bands, we lastly turn to the spectral region between 610–410 cm−1 which exhibits three broad features (Fig. 3H). Their assignment is straightforward when comparing it with the VPT2 calculation. The band at 527(2) cm−1 can be assigned to the CC–O in-plane bend ν10. As the predicted depolarisation ratio of the A′ vibration (0.73) is close to that of an A′′ band (0.75), it is not surprising that the experimental band mostly persists in the spectrum measured with parallel laser polarisation (green trace in Fig. 3H). On the higher-wavenumber side of ν10 of tPr, there is a hot band at 538(2) cm−1, which matches the predicted band position of ν10 of cis-propiolic acid within the experimental uncertainty. However, the hot band originating from the lowest-energy A′′ fundamental ν10 + ν15 − ν15 is predicted to be upshifted by 6 cm−1 by the VPT2 calculation. With B3LYP-D3(BJ)/aVQZ, an even larger anharmonicity constant x10,15 of +10 cm−1 is obtained that matches the band centre within the experimental uncertainty. As the expected thermal population of ν15 at the highest nozzle temperature is significant, the hot band at 538(2) cm−1 could also correspond to a trans-propiolic acid hot band.
Between 500–250 cm−1, no fundamental transition of trans-propiolic acid is expected (cf.Table 1). Thus, the tPr band at 455(2) cm−1 in Fig. 3H is likely due to a combination or overtone vibration of the two lowest-energy fundamentals ν11 and ν15. In the spectrum measured with parallel laser polarisation (green in Fig. 3H), only the broad socket of the band persists. Therefore, the band is assigned to the A′-symmetric 2ν15 overtone of trans-propiolic acid, which is predicted at 458 cm−1 by the VPT2 calculation.
Another interesting feature in this spectral region is the broad feature between 600–550 cm−1 with maxima at 584(2), 578(2), and 569(2) cm−1 as well as additional cluster contributions in the spectra recorded at the lowest two nozzle temperatures. From the VPT2 calculation, we expect two tPr fundamentals, ν9 (A′, ρ = 0.45) and ν14 (A′′, ρ = 0.75), to contribute to it. The depolarised spectrum facilitates the assignment of ν9 to 584(2) cm−1, as this band nearly vanishes in the parallel spectrum (green trace in Fig. 3H), and an assignment of the two maxima at 578(2) and 569(2) cm−1 to ν14, as these persist. A similar double maximum band shape is observed for the other A′′-symmetric fundamentals ν13 and ν15 (cf. Fig. S4 in the ESI†). It is interesting to note that the energetic order of ν9 and ν14 seems to be sensitive to subtle perturbations of the environment. The energetic order reported in argon, xenon, and krypton matrices44 is the same as in the gas phase, whereas an inverse order is observed in cryogenic nitrogen matrices (cf.Table 1) with ν14 shifted to higher wavenumbers.34 This sensitivity of the energetic order of ν9 and ν14 also extends to quantum chemical calculations, rendering the ν9–ν14 energy sequence excellent for benchmarking, which will be explored in the following section.
In addition to these computational means, the experimental uncertainty of the band positions has to be considered, which is a result of the resolving power of the monochromator and the wavenumber calibration. It amounts to ±2 cm−1 for all experimental band positions and is displayed as a grey area for ν9 in Fig. 4. As the band position of the OH torsion ν14 cannot be clearly ascertained from the spectrum, agreement with experiment is achieved when the shift is positive. An upper bound for the shift can be estimated by taking the lower-wavenumber maximum of the band at 569 cm−1 as the ν14 band position, which results in a (ν9–ν14)-shift of +15 cm−1, indicated by a dashed line in Fig. 4, yielding a ‘hard’ (>0) and a ‘soft’ criterium (≤15 cm−1) for agreement with experiment.
As the ν9–ν14 energetic shift is rather small and DFT VPT2 energies have previously been found to be very sensitive to the calculation settings (see for example Fig. 5 in ref. 52),84–86 we have averaged the anharmonic vibrational frequency calculations over different settings for each method shown in Fig. 4. The anharmonic force fields were computed with and without symmetry, tight and very tight geometry optimisation thresholds, and for the DFT methods using two different integration grids (ultrafine and superfine). The error bars shown in Fig. 4 correspond to the standard deviation of these eight calculations (four for MP2 and CCSD(T)). For all calculations, a basis set progression from aVDZ to aVQZ has been carried out, with the size encoded in the filling of the symbols (aVDZ □, aVTZ , aVQZ ■).
Fig. 4 shows that harmonically (circles), most methods predict the energetic order of the two fundamentals incorrectly. The two exceptions are the PBE0-D3(BJ) and B3LYP-D3(BJ) calculations with the larger aVTZ and aVQZ basis sets. Upon inclusion of vibrational anharmonicity (squares), most methods yield the correct sign for the shift. One exception is MP2 (red), which predicts the correct energetic order only for the aVQZ basis set with the very tight optimisation threshold and Cs symmetry. Generally, the standard deviation of the calculations with different settings (symmetry, optimisation criterion, and DFT integration grid) increases significantly with basis set size and is distinctly larger for the OH torsion than for the CCO in-plane bend in case of the DFT methods. Similarly large deviations were observed for a benchmark study of formic acid monomer where the OH torsion was also shown to be particularly sensitive to these settings.52 Unsurprisingly, the band position of ν9 is mostly overestimated by the harmonic calculations, which also extends to most VPT2 calculations using DFT. Only MP2 (with very tight optimisation and Cs symmetry) and CCSD(T)-edited VPT2 with large basis sets yield quantitative agreement with experiment for both experimental observables. With only a double-zeta quality basis set, electron correlation of the triple bond is severely underestimated, leading, amongst other things, to a soft force constant for this bond. Overall, it can be concluded that CCSD(T)-edited VPT2 yields the best result for the right reason. Out of the DFT methods, B2PLYP-D3(BJ) performs best with inclusion of vibrational anharmonicity. Similarly good agreement was previously observed for the formic acid monomer.52 Harmonically, B3LYP-D3(BJ) with an aVTZ basis set is the best choice. This good performance of harmonic vibrational frequency calculations at the B3LYP-D3(BJ) level for relative quantities has been observed for other carboxylic acids30,52,87 and various other hydrogen bonded clusters.88
In the following section, dimerisation of trans-propiolic acid shall be explored. As seen in Fig. 3 and 4, increasing the basis set size for the harmonic CCSD(T) and the anharmonic MP2 VPT2 force field can yield significant improvement of the anharmonic energy levels. CCSD(T)-F12a/VTZ-F12//MP2/aVQZ VPT2, which was shown to perform particularly well for the propiolic acid monomer, however, is computationally too demanding for the dimers. Therefore, a compromise for the two basis sets (harmonic and anharmonic) needs to be found. This will be explored for the fifteen fundamentals of trans-propiolic acid. For the harmonic force field, the performance of the smaller VDZ-F12 basis is compared to that of VTZ-F12. In case of the anharmonic MP2 force field, aVTZ as well as a split aVTZ + VQZ basis are tested. For the split basis, we use quadruple-zeta basis sets (VQZ) on the carbon atoms and aVTZ on the hydrogen and oxygen atoms. In Fig. 5, the performance of different basis set combinations is compared to the CCSD(T)-F12a/VTZ-F12//MP2/aVQZ reference, which illustrates two important findings: firstly, the difference between the VDZ-F12 and VTZ-F12 basis set for the harmonic CCSD(T)-F12a force field (yellow bars) is reasonably small with differences below 10 cm−1 and a mean absolute deviation of 2.5 cm−1. Notably, the VDZ-F12 basis yields smaller energies compared to the VTZ-F12 results. The second finding is that decreasing the basis set size of the anharmonic MP2 force field (aVTZ) while using CCSD(T)-F12a/VTZ-F12 harmonic wavenumbers yields on average an overestimation of the band positions compared to the MP2/aVQZ VPT2 reference. Fortunately, most of the missing electron correlation (compared to an all-aVQZ basis) seems to be recovered when the basis set size is increased only on the carbon atoms (split basis, green) which significantly reduces the shortcomings of the all-aVTZ basis (mean absolute deviation of 2.8 cm−1 instead of 4.8 cm−1). If both reductions of the basis set size (harmonic and anharmonic) are implemented at the same time (red in Fig. 5), the introduced errors tend to cancel fortuitously, which is reflected in a mean absolute deviation of 2.5 cm−1. We will therefore use CCSD(T)-F12a/VDZ-F12//MP2/aVTZ + VQZ throughout Section 4.
As the OH stretching region of carboxylic acid dimers is notoriously congested with a broad, complicated resonance structure, we focus on the CO stretching region in the following. Mixed cis–trans- or cis–cis-dimers are not expected to a considerable amount in the jet spectra due to the relatively low abundance of the cis-rotamers in the expansion at the nozzle temperatures explored for the dimer investigation (≤60 °C, resulting in 1–2% of cis-propiolic acid monomer). We will start with the infrared spectrum of the propiolic acid dimer (tPr–tPr) before it is compared to the infrared spectrum of the formic–propiolic acid heterodimer (tF–tPr), followed by an analysis of the Raman spectra of both dimers in Section 4.2. Heterodimers such as (tF–tPr) are no longer centrosymmetric and therefore, the C
O oscillators are not strictly equivalent, yet due to the local inversion symmetry of the carboxylic acid dimer ring, the infrared-activity of the symmetric and the Raman-activity of the antisymmetric C
O stretch are fairly small. As in the monomer section, the vibrational modes are labelled according to the Herzberg nomenclature of the respective dimer. To refer to the two carbonyl stretching vibrations, we introduce the short-hand notation ‘a’ (antisymmetric) and ‘s’ (symmetric) as indices to refer to the infrared- and Raman-active C
O stretching fundamental, respectively, e.g., (tPr–tPr)a for the infrared-active antisymmetric C
O stretch of the propiolic acid homodimer.
For a closer investigation of potential resonance partners of (tPr–tPr)a, anharmonic (VPT2) vibrational frequency calculations have been carried out using the same electronic structure methods as for the monomer benchmark (Fig. 4), namely B3LYP-D3(BJ), PBE0-D3(BJ), B2PLYP-D3(BJ), and MP2 all using an aVTZ basis set. Additional CCSD(T)-edited VPT2 results are shown, illustrating the convergence of the anharmonic MP2 force field for an all-aVTZ and the aVTZ + VQZ split basis. The results of the VPT2 calculations are plotted below the FTIR spectra in Fig. 6.
Fig. 6 shows that all methods consistently predict two Bu-symmetric combination vibrations with significant intensity between 2000–1600 cm−1, one in proximity to the antisymmetric CO stretching fundamental of the dimer (ν7 + ν31) and one upshifted from the monomer band (ν14 + ν20). The former is a combination of the symmetric (ν7) and antisymmetric OCO in-plane bending vibration (ν31) and the latter a combination of the symmetric (ν20) and antisymmetric OH torsion (ν14) (see insets of Fig. 6 for schematic representation of the vibrational motion). CCSD(T)-edited VPT2 predicts a third combination state with significant infrared intensity in that spectral region, namely a combination of the symmetric OH in-plane bend and a low-frequency mode that modulates the hydrogen bond length (ν5 + ν35). The ν7 + ν31 band position relative to the monomer and cyclic dimer varies across the methods – for all DFT methods, it is predicted closer to the monomer than to the dimer fundamental, with B3LYP-D3(BJ) showing the overall best band position agreement. In case of MP2, the relative distance between ν7 + ν31 and the antisymmetric C
O stretch (21 cm−1) roughly matches the observed splitting of 15 cm−1 between both bands, but the relative distance with respect to the monomer band as well as the intensity ratio between both bands is too small. An assignment of the 1712/1727 cm−1 doublet to the antisymmetric C
O stretch and ν7 + ν31 is consistent with the CCSD(T)-edited VPT2 results that predict an energy difference of 13 cm−1 (exp. 15 cm−1, intensity inversely). A firm assignment, however, is currently not possible as the alternative assignment to ν5 + ν35 cannot be ruled out with the available computational results. Note that except for MP2/aVTZ, which predicts a resonance to ν7 + ν31, neither DFT nor CCSD(T)-edited VPT2 predicts Fermi resonances involving the antisymmetric C
O stretch when default resonance thresholds are used.
The predicted band position of the second combination vibration with significant intensity (ν14 + ν20) matches the third, smaller dimer feature at 1784 cm−1 within the experimental uncertainty for CCSD(T)-edited VPT2 with the smaller aVTZ basis for the anharmonic force field. However, its band position is not yet converged, as it shifts down significantly (25 cm−1) when changing to the split basis for the anharmonic force field. While MP2 and B2PLYP-D3(BJ) predict ν14 + ν20 slightly higher in energy than CCSD(T)-edited VPT2, particularly large deviations are observed for B3LYP-D3(BJ) and PBE0-D3(BJ). Interestingly, the predicted infrared intensity of ν14 + ν20 is of similar order of magnitude across all methods (87–250 km mol−1, compare to (tPr–tPr)a in Fig. 6), though the band position difference to the fundamental varies from 268 cm−1 (B3LYP-D3(BJ)) to 37 cm−1 (coupled-cluster), with the smaller band position differences coinciding, perhaps somewhat accidentally, with smaller predicted infrared intensities. Altogether, this relatively high infrared intensity despite fairly large energy differences to the fundamental indicates that ν14 + ν20 may have intrinsic oscillator strength. For formic acid, such a low-intense dimer band upshifted from the monomeric CO stretching fundamental is also observed in the infrared spectra. In line with the results presented here for the propiolic acid dimer, it was assigned to a combination vibration of the symmetric and antisymmetric OH torsion by Kollipost,24 which is also predicted to exhibit significant oscillator strength.26
To cross-check this assignment, we look at an intermediate system, i.e., the formic–propiolic acid heterodimer (tF–tPr). FTIR spectra of the mixture of formic and propiolic acid are shown in Fig. 7 alongside predicted band positions of the heterodimer using CCSD(T)-edited VPT2. The FTIR difference spectrum obtained by subtracting monomer-scaled single substance spectra from that of the mixture (orange in Fig. 7) shows one distinct band at 1732 cm−1, which can be clearly attributed to the antisymmetric CO stretching vibration of the heterodimer (tF–tPr)a. In addition to the fundamental of the mixed dimer, which agrees within 2 cm−1 with experiment, the VPT2 calculation predicts two combination vibrations of (tF–tPr) with significant infrared intensity between 1820–1600 cm−1, namely ν23 + ν24 and ν8 + ν18 (cf. insets in Fig. 7 for vibrational motion). ν23 + ν24 corresponds to the aforementioned OH torsion combination vibration and ν8 + ν18 is predicted to the lower-wavenumber side of the mixed dimer fundamental. In the spectrum of the mixture (black in Fig. 7), there seems to be some intensity in-between the two dimer features of the propiolic acid homodimer, which could correspond to ν8 + ν18, yet due to the overall low-intensity, no conclusions can be drawn. Close to the predicted band position of the combination vibration of the symmetric and antisymmetric OH torsion (ν23 + ν24), there is a dimer band in the mixed spectrum at 1798 cm−1 (black in Fig. 7, see also Fig. S6 in the ESI†) similar to the dimer band in the propiolic acid homodimer spectrum (1784 cm−1). Interestingly, increasing the basis set size of the anharmonic MP2 force field (all-aVTZ to aVTZ + VQZ) shifts the prediction from 1826 cm−1 to 1788 cm−1 towards experiment (1798 cm−1) whereas the exact opposite is the case for the propiolic acid homodimer (1784 to 1758 cm−1; exp. 1784 cm−1), demonstrating that the anharmonic couplings of this OH torsion combination vibration are not yet fully converged.
Such an infrared-active, non-fundamental dimer band upshifted from the monomer is also observed for other acids such as acetic (1834.9 cm−1),25 pivalic (1838/1834 cm−1),26 and benzoic acid (1824 cm−1),27,28 which for all three acids so far remains unassigned. Considering that the vibrational motion of this OH torsion combination vibration only involves the hydrogen-bonded ring, one would expect it at a similar band position independent of the carboxylic acid substituent, given that it has sufficient oscillator strength or gains it via resonance. As seen for the propiolic acid and formic–propiolic acid dimer, the anharmonic VPT2 band positions of this combination state are not sufficiently converged and the quantitative agreement with experiment varies greatly depending on the basis set chosen for the anharmonic MP2 force field. To further investigate if this combination vibration can plausibly be ascribed to the so far unassigned infrared bands around 1800 cm−1 in the spectra of other carboxylic acid dimers, we turn to a simple harmonic model using B3LYP-D3(BJ) and alternatively explicitly correlated coupled-cluster, assuming similar anharmonic corrections for this vibration across all dimers. In Fig. 8, the sum of the harmonic wavenumbers of the symmetric and antisymmetric OH torsion of the formic, propiolic, acetic, pivalic, and benzoic acid homodimers and selected heterodimers26,30 is plotted against the experimental band position (for individual results, see Table S8 in the ESI†). One should note that for the formic-acetic and formic-pivalic acid heterodimers (daggers † in Fig. 8), the corresponding normal modes are no longer concerted pair vibrations, but mostly localised OH torsions on each monomer and instead of one, there are two normal vibrations in each case that can be labelled ‘Tora’ (cf. Fig. S8 in the ESI†). We therefore average the harmonic wavenumbers of the two respective normal modes. In Fig. 8, it can be clearly seen that for the homodimers (squares) and selected heterodimers (crosses) all points lie on straight lines with different slopes for B3LYP-D3(BJ) and coupled-cluster, strengthening this assignment. For formic acid, deuteration experiments have been carried out by Kollipost,24 showing that such a dimer band is present at a similar band position in the infrared spectrum of DCOOH (1821 cm−1versus 1811 cm−1 for HCOOH), but no dimer band is observed for HCOOD between 1780–1840 cm−1. A systematic extension to the dimers of oxygen-deuterated carboxylic acids therefore seems most promising for verification.
![]() | ||
Fig. 8 Sum of harmonic B3LYP-D3(BJ)/aVTZ and CCSD(T)-F12a/VDZ-F12 band positions of the symmetric (Tors) and antisymmetric OH torsion (Tora) for selected homo- and heterodimers plotted against the experimental band position of a non-fundamental dimer feature upshifted from the antisymmetric C![]() |
Below the spectra in Fig. 9, CCSD(T)-edited VPT2 results are shown using default Fermi resonance thresholds and inclusion of two-quantum states. The VPT2 calculations do predict resonances for the homo- (tPr–tPr)s and heterodimer (tF–tPr)s CO stretch. In the absence of anharmonic Raman scattering cross-sections, we inspect the wavefunction overlap with the deperturbed fundamental to check the assignments for plausibility.89,90 In Fig. 9, the squared wavefunction contribution of the (deperturbed) C
O stretch to the respective eigenstate is plotted against the anharmonic energy after diagonalisation of the effective resonance Hamiltonian for the propiolic acid homo- (Fig. 9A) and the formic–propiolic acid heterodimer (Fig. 9B). The absolute wavenumbers of the three main peaks in the propiolic acid homodimer spectrum (Fig. 9A) at 1667/1651/1645 cm−1 are surprisingly well reproduced by the VPT2 calculations. The assignments to the respective peaks vary between both basis sets used for the anharmonic MP2 force field, but given the experimental near 1
:
1
:
1 intensity, harmonic oscillator labels become insignificant in this context. With the smaller all-aVTZ basis, the triad is predicted at 1667/1655/1646 cm−1, with the larger split basis at 1672/1659/1644 cm−1. While the squared wavefunction coefficients do not quantitatively reproduce the observed intensity pattern, they demonstrate significant delocalisation of the symmetric C
O stretch, as is observed in the Raman jet spectrum. Future fine-tuning should include higher excited states in a variational treatment to account for the fine structure.
For the formic–propiolic acid heterodimer, two main peaks are seen at 1670 and 1657 cm−1 in the Raman jet difference spectrum (orange in Fig. 9B). The VPT2 calculations predict a Fermi resonance between the fundamental (tF–tPr)s and the combination vibration of the symmetric OH (ν24) and OH/CH out-of-plane bending vibration (ν25) at 1668/1652 cm−1 (aVTZ + VQZ) and 1671/1662 cm−1 (aVTZ) in good agreement with the experimental band positions. Depending on the basis set size of the anharmonic MP2 force field, the resonance between (tF–tPr)s and ν24 + ν25 is predicted to be either week (aVTZ + VQZ) or strong (all-aVTZ). A third, much weaker perturber is predicted further downshifted from the experimental bands at 1628 cm−1 (aVTZ + VQZ; 1624 cm−1 aVTZ), though is not observed in the spectra. Due to the low acid concentrations used for the measurement of the heterodimer spectra, however, additional weaker mixed dimer bands cannot be excluded.
Altogether, it can be concluded that the VPT2 calculations using default Fermi resonance treatment and including only binary states can predict the experimental band patterns surprisingly well. It is worth noting that the two predicted resonance partners for the homo- and heterodimer are in fact very similar normal modes (cf. normal modes in Fig. 9). Interestingly, all normal modes involved in these resonances almost exclusively comprise displacements of the cyclic ring, particularly (in- and out-of-plane) bending of the hydrogen bonded H atoms. Future work will have to show whether this also extends to the resonance perturbed CO stretching fundamentals of other carboxylic acid dimers.
The systematic comparison of the infrared and Raman jet spectra of the propiolic acid dimer, the formic–propiolic acid dimer as well as other carboxylic acid homo- and heterodimers in the CO stretching region revealed interesting similarities in their vibrational signature. A weak non-fundamental dimer band that is observed in the propiolic acid (1784 cm−1) and formic–propiolic acid infrared spectrum (1798 cm−1) can be assigned to a combination vibration of the symmetric and antisymmetric (hindered) OH torsion, which has also been observed for the formic acid dimer (1810.6 cm−1).24 A simple harmonic model (Fig. 8) indicates that this combination vibration might be a plausible explanation for the so far unassigned features in the infrared spectra of many other carboxylic acid dimers, such as the benzoic acid,27,28 the acetic acid,25 and the pivalic acid dimer,26,30 and the mixed dimers of formic, acetic, and pivalic acid.26,30 For validation, a systematic experimental extension to the dimers of oxygen-deuterated carboxylic acids is necessary.
As observed for other carboxylic acid dimers,29,30 the Raman-active symmetric CO stretching fundamentals of the propiolic acid dimer and its heterodimer with formic acid are strongly perturbed by resonances (see Fig. S5 in the ESI†). In the Raman spectrum of propiolic acid, a broad feature is observed spanning about 40 cm−1 with three distinct peaks, whereas for the corresponding formic–propiolic acid heterodimer, a resonance doublet is observed. CCSD(T)-edited MP2 VPT2 calculations used to guide spectral assignments throughout this work were shown to be able to reproduce the main features in the Raman spectra. In both cases, the calculations clearly illustrate that while the resulting resonance pattern varies, all perturbing states share a notable similarity: all involved normal modes have significant OH in- and out-of-plane bending character. This raises the question of whether the resonance perturbations of the Raman-active C
O stretching fundamentals of other carboxylic acid dimers can also be explained by only a few couplings between the C
O stretch and the (in- and out-of-plane) OH bend in the internal coordinate picture which map onto several vibrations in the normal mode representation (cf. related work by Florio et al.13 concerning the OH stretching spectrum of the formic and benzoic acid dimer). The present results call for a similar reduced-dimensional model calculation in the C
O stretching range of the propiolic acid and other carboxylic acid dimers.
Footnote |
† Electronic supplementary information (ESI) available: Depolarised Raman jet spectra of the propiolic acid monomer (including A′′ fundamentals), FTIR jet spectra of the propiolic acid monomer, list of all observed experimental monomer band centres, band positions and assignments of decomposition products, Raman jet spectra of symmetric C![]() |
This journal is © the Owner Societies 2021 |