Leak-out spectroscopy as alternative method to rare-gas tagging for the Renner–Teller perturbed HCCH+ and DCCD+ ions

Infrared messenger-tagging predissociation action spectroscopy (IRPD) is a well-established technique to record vibrational spectra of reactive molecular ions. One of its major drawbacks is that the spectrum of the messenger-ion complex is taken instead of that of the bare ion. In particular for small open-shell species, such as the Renner–Teller (RT) affected HCCH+ and DCCD+, the attachment of the tag may have a significant impact on the spectral features. Here we present the application of the novel leak-out spectroscopy (LOS) as a tag-free method to record the cis-bending of the HCCH+ (∼700 cm−1) and DCCD+ cations (∼520 cm−1), using a cryogenic ion trap end user station at the FELIX laboratory. We demonstrate that the obtained LOS spectrum is equivalent to a previously recorded laser-induced reactions (LIR) spectrum of HCCH+. The bending modes are the energetically lowest-lying vibrational modes targeted with LOS so far, showing its potential as a universal broadband spectroscopic technique. Furthermore, we have investigated the effect of the rare gas attachment by recording the vibrational spectra of Ne- and Ar-tagged HCCH+. We found that the Ne-attachment led to a shift in band positions and change in relative intensities, while the Ar-attachment even led to a complete quenching of the RT splitting, showing the importance of using a tag-free method for RT affected systems. The results are interpreted with the help of high-level ab initio calculations in combination with an effective Hamiltonian approach.


Introduction
3][4] It is proposed to be a key building block in the formation of larger organic molecules such as polycyclic aromatic hydrocarbons, e.g., via the hydrogen abstraction-acetylene addition (HACA) mechanism in circumstellar shells, 5,6 or via ion-molecule reactions in the cold interstellar medium. 7Acetylene itself has been detected in space, 8,9 but its cationic counterpart, HCCH + , likely formed upon UV or cosmic-ray ionization in the ISM and expected to play an important part in the formation of larger hydrocarbon cations, [10][11][12] has so far remained undetected.
This interesting cation, exhibiting a 2 P 3/2 electronic ground state, is also relevant from a fundamental spectroscopic point of view due to its open-shell linear nature.It is subjected to both spin-orbit and Renner-Teller (RT) coupling effects complicating the rovibrational spectra of its trans-(n 4 ) and cis-bending (n 5 ) modes of which only the latter is IR active.][16][17][18][19] The spectroscopic studies on the RT affected bending modes of this ion include low-and high-resolution pulsed-field induced zero electron kinetic energy photoelectron spectroscopy (PFI-ZEKE PES), [20][21][22] laser induced reactions (LIR), 23,24 and Ar-tagging infrared predissociation (IRPD) spectroscopy. 25Investigations into its doubly deuterated form, DCCD + , are limited to two theoretical studies by Peric ´and Radic ´-Peric ´26 and Jutier and Le ´onard 16 and a low-resolution photoelectron spectroscopic study by Reutt et al., 27 revealing one feature of the RT split trans-bending, but none of the cis-bending.
This abundance of spectroscopic and theoretical data makes HCCH + an excellent system to test the potential of the newly developed action spectroscopic technique leak-out spectroscopy (LOS). 28With this technique, the mass-selected ions, stored in a low-temperature ion trap, are vibrationally excited after which they collide with a neutral buffer gas, in turn transferring their internal energy to kinetic energy.As a result of this extra kinetic energy, the ions ''leak-out'' of the trap and are then detected, giving rise to a nearly background-free vibrational spectrum.The LOS method is thus very universal, in particular compared to another tag-free method, laser-induced reactions (LIR 29,30 ), which has also been applied previously to the cis-bending of HCCH + , 23,24 but requires a suitable slightly endothermic reaction to be applicable.Leak-out spectroscopy has already been applied for the high-resolution ro-vibrational spectroscopy of C 3 H + , 28 ro-vibrational and double-resonance rotational spectroscopy of c-C 3 H 2 D + , H 2 CCCH + , and HCCCO + , 31 and to the RT inactive antisymmetric C-H stretch (n 3 ) of HCCH + . 32n these studies, the infrared excitation leading to collisional energy transfer was done in the 3 micron region, using narrowbandwidth continuous-wave, or in some cases also ns-pulsed, radiation.Our work presented here on the low-energetic bending modes of HCCH + (B700 cm À1 ) and DCCD + (B500 cm À1 ) using a ms-pulsed infrared free-electron laser may therefore act as proofof-principle for LOS as a broadband spectroscopic technique.
Another spectroscopic technique to record broadband vibrational spectra of reactive species under isolated conditions is IRPD action spectroscopy.8][39][40][41] However, for small open-shell species, such as HCCH + and DCCD + , the interaction between the RG tag and the ion may have a significant impact on the bending potential energy curve and thus its spectral features, as seen previously for Netagged HC 3 N + . 42So far, to our knowledge, no systematic study exists on the effect of the rare-gas attachment on Renner-Teller perturbed systems.
In the present work we have recorded the spectrum of the IR active cis-bending of HCCH + with LOS and Ne-and Ar-tagging IRPD, adding to the existing LIR work of Schlemmer et al. 23 and Asvany et al. 24 Not only do these measurements act as proof-ofprinciple for the application of LOS on low-energetic vibrational modes, but they also allow for a systematic investigation of the RG attachment to RT perturbed systems.Furthermore, we have recorded the first vibrational spectrum of the DCCD + ion by means of LOS and Ne-tagging IRPD.Using an effective Hamiltonian approach the observed structures were analyzed and the resulting spectroscopic constants were compared to those obtained in the earlier PFI-ZEKE PES and LIR studies.

Experimental methods
The experiments were performed in the cryogenic 22-pole ion trap instrument FELion located at the Free Electron Lasers for Infrared eXperiments (FELIX 43 ) laboratory.A detailed account of the instrument and the employed action spectroscopic technique IRPD has been given previously. 44For the LOS experiments we follow the method described by Schmid et al. 28 Here we give a brief account on the details specifically related to the IRPD and LOS measurements of the HCCH + and DCCD + ions.
The HCCH + ions are produced by electron impact (EI) ionization of acetylene (C 2 H 2 Z 98%, Linde Gas) and the DCCD + ions from acetylene-d2 (C 2 D 2 , 99%-d2, CDN Isotopes).Both gases were diluted with He in a 1 : 3 (precursor : He) mixing ratio to aid ionization and the pressure in the source was kept at B10 À5 mbar.Fig. 1 shows the EI mass spectra (30 eV electron energy) of the HCCH and DCCD precursors overlayed with the background signal from residual gas.We dominantly form the radical cations with m/z 26 (HCCH + ) and m/z 28 (DCCD + ), respectively, where for the latter Fig. 1 Mass spectra of ions produced from electron impact ionization (B30 eV) of (a) HCCH and (b) DCCD overlayed with residual gas mass spectra.
Here the first quadrupole functions as an ion guide from the source to the trap and the second quadrupole is used for scanning over the m/z range.
This journal is © the Owner Societies 2024 a small admixture of CO + and N 2 + is present in the same mass channel.
The ions are extracted from the source in a B15-50 ms long pulse and are introduced into a quadrupole mass filter to mass select the ions of interest.After mass selection the ions are guided into the 22-pole ion trap, 45 where either an IRPD or LOS spectroscopic scheme is applied.
For IRPD action spectroscopic experiments complexes of the target ions with Ne or Ar are formed inside the ion trap via three-body collisions by storing them at 6-7 K for Ne and 13-18 K for Ar and exposing them to a short (few 10 ms) He : Ne(Ar) pulse with high number density (B10 15 cm À3 ) with mixing ratio 3 : 1 for He : Ne and 5 : 1 for He : Ar.Under these conditions B17% of m/z 26 (i.e.HCCH + ) forms a weakly bound complex with Ne and B5% with Ar (see Fig. 2).For m/z 28 only B10% binds with Ne, which is likely due to the presence of both N 2 + and CO + in the m/z 28 mass channel.These ions may have lower Ne binding energies than DCCD + , resulting in an overall lower tagging efficiency.In the trap the formed complexes are then exposed to FELIX infrared light for 1-3 s, inducing predissociation upon on-resonant irradiation.Then, the content of the trap is analyzed with a second quadrupole mass selector and a single-ion counting detector, monitoring the depletion of the ion-complex signal as a function of the FELIX wavelength.
For leak-out spectroscopy the target ions are stored and collisionally cooled in the ion trap at 15-17 K by exposing them to a short (few 10 ms) He pulse at high number density (B10 15 cm À3 ).The trap exit electrode potential is set so that the ions are just stored and that only a small kinetic energy gain is needed for them to leak out.These settings are continuously adjusted to account for changes in the effective voltages due to condensation.While the ions are stored, neon gas is continuously flowing into the ion trap (B10 12 cm À3 ) to act as a collision partner, which is needed to convert the internal energy of the excited ions to kinetic energy.The number density and trap temperature were optimized for maximum LOS signal strength and avoidance of formation of Ne-complexes.In the trap, the bare ions are exposed to FELIX infrared radiation for 1-3 s and leak-out upon on-resonant irradiation.After each cycle the trap content is discarded and the process repeats.To record the vibrational spectrum one may either monitor the depletion of the bare-ion signal as a function of the FELIX wavelength after each trapping cycle, or directly monitor the ions leaking out during the trapping period.In this work the former method was used.
The 440-600 cm À1 and 650-840 cm À1 wavenumber ranges were covered in this study using the free-electron IR laser FEL-2 of the FELIX Laboratory ‡ with a macropulse repetition rate of 10 Hz, maximum pulse energy in the trap region of o20 mJ (at 700 cm À1 ), and a full width at half maximum (FWHM) of around 0.7% of the center wavenumber.To account for varying laser pulse energy E, pulse number n, and for saturation effects, the signal is normalized prior to averaging using , giving the intensity I in units of crosssection per joule as a function of the frequency (n), where N ON (n) gives the on-and N OFF the off-resonance ion count.The individual measurements are averaged using statistical binning with a bin size of 1 cm À1 to obtain the final spectrum.The final spectra were normalized with maximum intensity equal to one, to facilitate comparison between different methods.
2.2 Theoretical methods 2.2.1 Ab initio calculations.Ab initio quantum chemical calculations were used to optimize the geometric structure and calculate the harmonic vibrational frequencies of the bare HCCH + and DCCD + species and to compute the potential energy curve of the Ne-and Ar-tag around the HCCH + ion.All calculations were performed at the partially spin-restricted coupled cluster level of theory, with single, double, and perturbative triple excitations [RCCSD(T)]. 46For geometry optimization and subsequent harmonic wavenumber calculations on Here the first quadrupole is used to mass select the ion of interest before entering the trap and the second quadrupole is used for scanning over the m/z range after trapping.‡ https://www.ru.nl/en/hfml-felix.
the HCCH + and DCCD + ions Dunning's aug-cc-pV5Z basis set [47][48][49] was used.To investigate the interaction of HCCH + with the Ne and Ar atom a one-dimensional cut of the potential energy curve (aug-cc-pVTZ) [47][48][49] was made in Jacobi coordinates by attaching the Ne/Ar atom to the center of mass of the HCCH + .In these calculations the attachment angle was kept fixed while all other geometry parameters were allowed to relax.The computed interaction energies were counterpoise corrected to account for the basis set superposition error.All quantum chemical calculations were performed using the MOLPRO suite, version 2015.1. 50.2.2 Effective Hamiltonian.To simulate the spin-vibronic levels of HCCH + and DCCD + an effective Hamiltonian approach was used similar to that of Yang and Mo, 20 Tang et al.21 and Steenbakkers et al. 42 The chosen basis is given by: |nPJM;pi (1)   where J is the total angular momentum quantum number, M and P are quantum numbers corresponding to projection of the total angular momentum on the space-fixed and molecule-fixed axes, respectively, p = AE1 is the parity, and Here, S = 1/2 is the total electron spin quantum number and v 5 = 0, 1,. . . is the vibrational quantum number of the cisbending mode n 5 .We only include the cis-bending mode since the trans-bending mode (n 4 ) is IR inactive.Coupling with the symmetric (n 1 ) and anti-symmetric C-H (n 3 ) stretch as well as with the C-C stretch (n 2 ) were ignored.Furthermore, L = AE1, S = AE1/2, and l 5 = Àv 5 , Àv 5 + 2,. ..,v 5 are the quantum numbers corresponding to projection on the molecular axis of the electronic orbital, total electron spin, and vibrational angular momenta, respectively.Quantum number K = L + l 5 represents the projection of the total angular momentum excluding electron spin and P = K + S. Quantum numbers J, M, and p are good quantum numbers and the Hamiltonian is independent of M.
In our model (excluding rotation) K and P are also good, since we only include diagonal spin-orbit coupling and first order RT and neglect Coriolis coupling.Because of this approximation states with parity p = +1 and p = À1 are degenerate, and we can compute the Hamiltonian by restricting the basis to P 4 0, without parity adapting the basis.For a basis truncated at v 5 = 10 we find that energy levels are converged up to v 5 = 4.We approximate the total effective Hamiltonian by: where H ˆvib represents the vibrational energy, with o 5 the harmonic frequency of the cis-bending modes n 5 and x 5 its anharmonic correction.The spin-orbit Hamiltonian is given by 51 where A SO represents the spin orbit constant and L ˆz and S ˆz the molecule fixed components of the electronic orbital and spin angular momenta operators, respectively.The effective RT Hamiltonian is: 52 where the constant g 5 is related to the dimensionless RT constants e 5 as In this work o 5 and e 5 are either fitted or approximated from ab initio harmonic frequency calculations as and where o 5 (A 0 ) is the ab initio harmonic frequency of the in-plane and o 5 (A 00 ) of the out-of-plane cis-bending.Here eqn ( 8) and ( 9) are derived for a first order Renner-Teller problem described in ref. 53.
The complex normal mode operators Q ˆ5,AE are related to the Cartesian normal modes Q ˆ5,x and Q ˆ5,y by and L ˆAE2 represent the orbital angular momentum ladder operators coupling states L = 1 and L = À1.Finally, the g K 5 term, first introduced by Brown, 52 is a correction to the RT Hamiltonian as a result of the K-dependent mixing of the vibrational levels of a P electronic state with the vibrational levels of S and D states.This term is dependent on the molecule fixed components of the electronic orbital (L ˆz) and vibrational angular (G ˆ5,z ) momenta.Finally, the rotational Hamiltonian is given by where B is the rotational constant.In this equation the final two terms take P-state mixing into account rendering P to not be a good quantum number.For the S states this interaction is strongest, while for the P and D states it is moderate for the lower J-states.From this point on we will therefore label the P and D states with their largest P-state contribution (as approximated by the expectation value of J z ).The wave function is expanded in the basis The expansion coefficients u (Jp) n,i are determined variationally by solving the matrix eigenvalue problem given in eqn (13) below, using the free and open source numerical software SCILAB version 6.1.1. 54 i = E i u i , (13)   where E i are the eigenvalues and H the Hamiltonian matrix of which the matrix elements are given in ref. 55.The population of the computed energy levels is then approximated with the Boltzmann distribution.
The line strengths for the DK = AE1 transitions are computed by A derivation of this expression is given in the ESI, † where we also define the ''time-reversal notation'' for % P 0 in the 3Àj symbol.This derivation of the line intensity is based on literature and text-book results.It is provided for the convenience of the reader, and to define exactly how the calculation was done.The dipole operator mk,5 is approximated by The perpendicular dipole moment, m > 5 , is set equal to 1, since all computed intensities are normalized and only one vibrational mode is considered.

Leak-out spectroscopy (LOS) of HCCH +
Fig. 3 shows the spectrum of the cis-bending mode of the HCCH + ion recorded with LOS overlayed with the previously recorded LIR spectrum. 23The positions of the most intense features match within the experimental uncertainty of the band center (B2 cm À1 , taking into account the wavenumber calibration), but a clear difference is observed in the band profiles.To explain these differences the ro-vibrational spectra were simulated with the effective Hamiltonian approach described in Section 2.2.2.The spectroscopic constants were kept fixed to experimentally obtained values from Tang et al. 21and the temperature and line width (FWHM) of individual rovibrational lines were adjusted (''fitted'') by hand to match best the observed band profile (see Table 1).The resulting simulations (shown in Fig. 4) are in good agreement with the experimental spectra, though the intensities of the RT affected S bands are either over-predicted (for the mS or under-predicted (for the kS)).This may be a consequence of (i) neglecting coupling to other vibrational modes in our model, (ii) a potential discrepancy between the collisional cross-sections of the mS and kS modes or (iii) the sensitivity of the leak-out efficiency possibly leading to effects in the intensity of the vibrational energy (see Section 3.2.1 for further discussion).High-level quantum dynamical computations are thus likely needed to interpret the observed intensities.Using the simulated spectra the observed band profiles can be explained and the most intense bands in both spectra can be assigned to belong dominantly to the mS, D 5/2 , and kS transitions from the P 3/2 vibrational ground state at approximately 700, 712, and 750 cm À1 , respectively.Furthermore, the transition to the D 3/2 state from the excited P 1/2 spin-orbit state is about 2 cm À1 shifted from the P 3/2 -D 5/2 band and will contribute to the 712 cm À1 feature.These transitions obey the DK AE 1 and DP = AE1 selection rules (see panel in Fig. 4), of which the latter is somewhat relaxed due to the P-state mixing.Clearly, the rotational structures of these bands overlap, complicating the assignment of the convoluted spectra.At higher T, this convolution becomes more problematic and can be attributed to several causes.First, more J levels are occupied resulting in an effective ''smearing'' of the substructure.Second, with the occupation of higher J-levels the P-state mixing becomes more prevalent, loosening the DP = AE1 selection rule.Finally, the contribution of transitions from the P 1/2 excited state (B30 cm À1 above the ground state) starts to increase: at 38 K (simulation temperature for the LOS experiment) this contribution, determined assuming a Boltzmann distribution, is 26.5%, while at 90 K (LIR experiment) it is 40.2%.A full list of calculated transitions with their respective intensities can be found in ESI: † Table S1 for the S and ESI: † Table S2 for the D states.
The good match of the experimental infrared spectra with the simulations using the spectroscopic parameters obtained by Tang et al. 21provides some clarification regarding the discrepancies between the Renner-Teller parameters for the cis-binding mode reported in the three experimental works: Yang and Mo 20 reported o 5 = 704.05cm À1 and e 5 = 0.0219, Tang et al. 21o 5 = 719.39cm À1 and e 5 = À0.0274 and Schlemmer et al. 23 o 5 = 710 cm À1 and e 5 = 0.032.Most striking is the significant difference between spectroscopic parameters obtained by Yang and Mo 20 and Tang et al., 21 since both works use the same, high-resolution, experimental technique (zero-kinetic-energy photo-electron spectroscopy).The effective Hamiltonians used for the fitting were also very similar with the largest difference being the inclusion of an anharmonicity term and the Brown g K 5 Renner-Teller correction by Tang et al. 21We have simulated the vibrational spectra using both Hamiltonians with their corresponding parameters and found a far better agreement with the parameters of Tang et al., 21 where the inclusion of the g K 5 term turned out to be critical to fit our experimental data.
It is noteworthy that the fitted rotational temperature of the LOS spectrum (38 K) is significantly higher than the nominal trap temperature (16 K), while for the LIR spectrum this is not the case (both recorded and fitted temperatures are 90 K).This effect of the LOS technique has been previously observed by Asvany et al., 31 where they discuss that collisions of the ions with heaver species (i.e.Ne) in an RF field leads to heating of the stored ions.The line width is expected to be dominated by the resolution of the free-electron laser whose FWHM may vary between 0.6-0.8% of the center wavenumber, corresponding to Do = 4.2-5.6 cm À1 at 700 cm À1 and Do = 3-4 cm À1 at 500 cm À1 .The fitted FWHM are in agreement with these values.The fact that the structure of both the LIR and LOS spectra can be reproduced by the simulation at different temperatures, indicates that the discrepancy between the obtained spectra is mainly a result of the temperature difference.Furthermore, the simulations show that recording a vibrational spectrum at different temperatures may be of great help in assigning the observed structure and validating the used model, especially for a relatively low-resolution technique.

Leak-out spectroscopy of DCCD +
Expanding on this proof-of-principle experiment, LOS was used to record the first vibrational spectrum of the RT-perturbed cisbending mode of DCCD + , after which it was simulated with the model described above (see Fig. 5).Regarding the spectral shape a similar band profile is observed compared to that of the HCCH + ion (an overlay of the structures is given in the ESI: † Fig. S1).Since isotopic substitution should not alter the electronic structure of the ion the A SO , e 5 (=g 5 /o 5 ) and x 5 parameters were kept fixed to those of HCCH + .The rotational constant B = 0.80 cm À1 was calculated the ab initio at RCCSD(T)/aug-cc-pV5Z level of theory and the parameters that are dependent on the ion mass (o 5 and g K 5 ), temperature and line width were fitted to the band contour.The best fit spectroscopic and experimental parameters are given in Table 1.Fig. 4 Infrared spectra simulated with the effective Hamiltonian model described in Section 2.2.2.The spectroscopic parameters were taken from Tang et al. 21and the temperature and FWHM were fitted to the spectra obtained with LOS (a) and LIR (b), 23 recorded at a nominal trap temperature of 16 K and 90 K, respectively.The best fit internal temperature is 38 K for the LOS and 90 K for the LIR measurement.The corresponding spectroscopic and fit parameters are given in Table 1.The inlay in (b) shows an illustrated energy level diagram (not to scale).
This journal is © the Owner Societies 2024 From ab initio calculations the harmonic vibrational frequency and Renner-Teller constant can be extracted for both species.The g K 5 parameter, which captures the mixing of the excited 2 S in the ground 2 P electronic state upon the bending of the ion, decreases with deuterium substitution.This is in line with previous findings of Smith et al. 56 , who investigated deuterium substitution in SiCH and found that . This equation only holds when the coupling and energy gap between the excited 2 S in the ground 2 P is isotopically invariant.From the fit we get 719:4 2 528:5 2 ¼ 1:85, which is in reasonable agreement considering that the DCCD + parameters were fitted by hand.This agreement indicates that the isotope substitution of HCCH + indeed only has a minimal impact on the electronic landscape and shows that isotope substitution can be extremely useful in understanding the complex ro-vibrational structure of Renner-Teller affected species and can aid in justifying the used model and parameter space.
Since this is the first vibrational spectrum of DCCD + a comparison with earlier theoretical works is imperative.Table 2 shows the rotational band origins of the four vibrational bands of n 5 (m 2 S 1/2 , 2 D 5/2 , 2 D 3/2 , and k 2 S 1/2 ) obtained from the fit together with the ab initio values of earlier theoretical work by Peric ´and Radic ´-Peric ´26 and Jutier and Le ´onard. 16Both works are based on a variational approach, where Peric ´and Radic ´-Peric ´used reduced degrees of freedom and Jutier and Le ´onard a full set of valence coordinates.Though both works report fairly different values, the divergence from the experimental work, here captured by the root-means-squared (RMS) of the difference between the calculated and observed values, is almost equal.The largest deviations are found for the RT affected 2 S 1/2 states, where Peric ´and Radic ´-Peric ´26 underestimate and Jutier and Le ´onard 16 overestimate the Renner-Teller splitting.Furthermore, both works overestimate most of the predicted band positions.Therefore, with lower E vib also DE kin will decrease.The cisbending of HCCH + has a vibrational energy of B700 cm À1 (B0.087 eV) and will consequentially result in a maximum kinetic energy gain of only B300 cm À1 or B0.037 eV.For the case of DCCD + the vibrational energy is even lower (B520 cm À1 ) resulting in a maximum gain of B216 cm À1 (B0.027 eV).
To date, these bending modes are by far the lowest-energy vibrational modes recorded with LOS and act as proof-ofprinciple for obtaining broadband vibrational spectra, showing the capabilities of this technique in combination with the FELIX free electron lasers.
To determine the effect of the vibrational energy on the observed leak-out intensity one needs to take the height of the trapping barrier into account.Ideally the height of the trapping barrier is set so that (almost) all ions are trapped but can still leak-out upon resonant irradiation.The setting of the trap barrier therefore affects the leak-out efficiency as shown by Bast et al. 57 As a consequence, a lower gain of kinetic energy would reduce the leak-out intensity.In order to find out to what extent intensity variations within a band are influenced by this effect the intensities of the cis-bending vibrations for HCCH + at around 700 cm À1 and DCCD + at around 520 cm À1 should be Fig. 5 Vibrational infrared spectrum of the cis-bending of DCCD + recorded by means of LOS at a nominal trap temperature of 16 K.The experimental spectrum is overlayed with a spectrum simulated with the effective Hamiltonian model described in Section 2.2.2 using a best fit internal temperature of 27 K.Used spectroscopic parameters are given in Table 1.carefully analyzed.Only then one could judge whether a significant intensity variation within one of the abovementioned bands could be associated with a changing efficiency of the leak-out process.Since LOS is a considerably new method many such systematic studies are needed to get a better account on its limits.Notably, the discrepancies between the measured and calculated intensities within the cis-bending vibrations of HCCH + and DCCD + only occur in the Sigma bands.This is an indicator that the deviation might be associated with an incomplete theoretical treatment, e.g., by ignoring couplings to other modes.Also here more systematic studies are needed to explain the current findings.

Infrared predissociation spectroscopy of HCCH +
The well-established infrared predissociation action spectroscopy (IRPD) provides an alternative method to record the spectra of small reactive ions.With this technique a messenger atom is weakly bound to the ion and its vibrational spectrum is taken as a proxy for that of the bare ion.We recorded vibrational spectra of HCCH + with the IRPD method using Ne and Ar as tags.Earlier work on the Ar-tagged IRPD spectrum of HCCH + by Relph et al. 25 did not allow for an analysis on the Renner-Teller structure due to the low signal-to-noise in the 600-1000 cm À1 region.The resulting spectra are shown in Fig. 6 together with the LOS spectrum of HCCH + , previously discussed in Section 3.1.
Very clear differences can be observed between the three spectra: first, for both IRPD spectra no rotational substructure is observed due to the decrease of the rotational constant upon rare-gas attachment.Second, the Ne-tag induces a shift in band position together with a change in relative intensity altering the Renner-Teller pattern.The Ar-tag results in a collapse of this structure, indicating a complete quenching of the Renner-Teller splitting.To explain this behavior, we performed calculations on the geometrical structure of the two RG-HCCH + complexes.A scan of the potential energy curve was computed at RCCSD(T)/aug-cc-pVTZ level of theory, where the Ar or Ne atom were displaced around the center of mass of HCCH + .The counterpoise corrected binding energy was then calculated as a function of the attachment angle, which was kept frozen, while all other geometric parameters were allowed to relax.For both Ar-and Ne-HCCH + the overall structure remained planar so that a symmetry plane can be distinguished (here xz).The label A 0 corresponds to a symmetric and A 00 to an antisymmetric wavefunction with regard to this plane, corresponding to configurations (p x ) 1 (p y ) 2 and (p x ) 2 (p y ) 1 , respectively.The resulting potential energy curves are shown in Fig. 7.
For the Ar-attachment the global minimum is found at a bent position (B701), where the Ar atom attaches with a binding energy of B1200 cm À1 , which is in agreement with earlier computational work of Dopfer et al. 58 and experimental work of Relph et al. 25 , where the n 2 C-C stretch was observed, indicating symmetry breaking effects as a result of strong Ar attachment.Furthermore, the here recorded spectrum clearly shows that the bent attachment results in a quenching of the Renner-Teller splitting, meaning that the complexation with the argon atom has lifted the degeneracy of the P orbitals enough so that the Born-Oppenheimer approximation holds.
It is striking that the estimated binding energy (B1200 cm À1 ) is much larger than the observed mode energy (B720 cm À1 ).This could mean that the binding energy of the Ar is greatly overestimated by the calculations, but this seems rather unlikely due to the high level of theory used together with the Counter-Poise correction applied here.Another explanation is that the fragmentation of the Ar-complex occurs through a two-photon absorption.In order for this to happen, the difference between the v = 0-1 and v = 1-2 transitions should fall within the FELIX bandwidth.If the system would still be RT affected one would not expect this to occur, since the splitting pattern in v = 2 is different than in v = 1.However, since for the bent Ar-HCCH + complex the RT splitting is quenched this difference is only determined by the anharmonicity and can be calculated with 2o 5 x 5 = 5.37 cm À1 , which falls well within the FELIX bandwidth (FWHM = 5.5 cm À1 ).
In contrast, the Ne-HCCH + PES reveals two minima of similar binding energy (D 0 B 250 cm À1 ) at linear and bent (B701) geometry with a B40 cm À1 barrier.However, after correcting for the zero-point energy (see Table 3) the bent geometry is clearly favored (D 0 B 211 cm À1 ) compared to the linear one (D 0 B 96 cm À1 ).Seeing that the Renner-Teller splitting remains visible one might argue that the Ne must be attached in a linear fashion.However, the Ne atom is much more weakly bound than the Ar-atom, so that the lifting of the degeneracy of the P orbitals in the bend minimum might not be enough to fully quench the Renner-Teller effect.To test this hypothesis we have recorded the (originally IR inactive) C-C stretch region of the Ne-HCCH + complex (see ESI: † Fig. S2), but only a very weak signal could be observed (within the noise level), indeed pointing towards a weak attachment of the Ne atom.Furthermore, with a barrier of only 40 cm À1 between the bent and linear configuration, one might expect a large This journal is © the Owner Societies 2024 amplitude motion of the Ne sampling all geometries.It is beyond the scope of this work to calculate the behavior of the Ne-tag and its impact on the vibrational landscape ab initio, since it requires a full diabatic treatment including the messenger atom.However, if we assume that the Ne-atom is bound in a linear fashion we can calculate its rotational constant ab initio, fit the Renner-Teller parameters and make a comparison with parameters of the bare ion.To verify the used model this procedure was repeated for the Ne-tagged DCCD + , which shows similar behavior upon Ne tagging as HCCH + (ESI: † Fig. S3).To reduce the number of fitting parameters the spin-orbit constant A SO and anharmonic constant x 5 were kept fixed.Furthermore, the temperature T was set to the recorded trap temperature (7 K), since no rotational substructure was visible.The resulting fits are shown in Fig. 8 and the corresponding spectroscopic parameters in Table 4.
For both the Ne-HCCH + and Ne-DCCD + all fitted parameters seem to change in a similar manner: the harmonic Table 3 Frequencies of linear and bent Ne-HCCH + complexes and the bare HCCH + ion calculated at RCCSD(T)/aug-cc-pVTZ level of theory.All values are given in cm À1 except the dimensionless Renner-Teller constant e 5 , which is given between brackets for the bending modes as (o 5 , e 5 ) Fig. 8 Infrared spectra of (a) Ne-HCCH + and (b) Ne-DCCD + simulated with the effective Hamiltonian model described in Section 2.2.2.The corresponding spectroscopic parameters are given in Table 4.
frequencies (o 5 ) are blue-shifted, the RT constants (e 5 ) decrease, the Brown parameters (g K 5 ) increase and the FWHMs increase with respect to the bare ions.The blue-shift of the harmonic frequency is in line with the ab initio calculations for the linear Ne-HCCH + complex (see Table 3), but is much less pronounced than predicted.Furthermore, the RT parameter e 5 decreases, while an increase is predicted from the harmonic frequency calculations and the g K 5 increases with about one third, which is much more than would be expected if the electronic structure would be invariant to the neon attachment (then equation would hold).These findings indicate that the Ne-attachment has a significant impact on the electronic landscape of Ne-HCCH + and Ne-DCCD + .Finally, the fitted FWHM are about a factor 2 larger than those of the bare-ion LOS and the linewidth of FELIX (3.3-6.6 cm À1 ).This cannot be fully attributed to a saturation effect, since the maximum depletion of the DCCD + was only B30% and we correct for this in the data analysis.One could argue that the temperature of the Ne-HCCH + and Ne-DCCD + complexes might be significantly higher than the set trap temperature of 7 K, which may result in a broader rotational profile.To test this we have simulated the Ne-HCCH + spectrum at 30 K and 50 K (see ESI: † Fig. S4).These figures clearly show that the broadness of the simulated spectra is dominated by the FWHM of the Gaussian convolution at all temperatures.The observed linewidth can therefore not be related to the rotational temperature of the complexes.Another cause of this broadening can be linked to life-time or pre-dissociation broadening effects, where a linewidth of B8 cm À1 would correspond to a very short lifetime of about 0.66 ps, which seems unlikely looking at the stability of the complex.Furthermore, it seems that the FWHM of the lines in the D band (highest intensity band) is larger than that of the S bands, since their contour is not well matched when using a universal FWHM for all lines.This could indicate that the vibronic angular momentum K may have an impact on the stability of the Ne-ion complex, but more likely is that the D band is overlayed with the rotational structure of the bent Ne-HCCH + complex.The latter hypothesis is strengthened by the spectrum of the Ar-HCCH + complex, which is much stronger bound than the Ne-HCCH + (and should therefore have a longer lifetime), but is of similar linewidth.Of course higher-level calculations are needed to fully describe the effect of the Ne-attachment on the electronic structure, but these calculations and measurements are in line with a large amplitude motion of the Ne atom around the HCCH + , sampling both geometric minima.

Conclusions
In this work we have recorded the vibrational spectrum of the low-lying (B700 cm À1 ) cis-bending of HCCH + with the novel spectroscopic method leak-out spectroscopy (LOS), 28 acting as a proof-of-principle for the application of this method to lowenergy vibrational modes.The effective Hamiltonian fit allowed to directly compare our work to the spectrum previously recorded at 90 K by means of LIR 23,24 , allowing to attribute the differences between the spectra to the ion temperature.Finally, this work provided some much-needed clarification regarding the large discrepancies between the spectroscopic parameters of the RT-perturbed cis-bending mode of HCCH + reported in different experimental works.The excellent agreement between our work and the parameters derived from PFI-ZEKE PES by Tang et al. 21leads us to believe that their data-set provides the most accurate description of the HCCH + ion.Extending on these measurements we have recorded the first vibrational spectrum of the cis-bending mode of DCCD + and the corresponding spectroscopic parameters were extracted with the same model as used for HCCH + .The obtained parameters may be used for future benchmarking studies related to isotope substitution in RT affected species.
In addition, we have recorded the vibrational spectra of HCCH + and DCCD + by means of the commonly used infrared predissociation (IRPD) action spectroscopic method, using Ne and Ar as messenger atoms.We found that the Ne-attachment leads to a distortion of the RT splitting and Ar to a complete quenching, which is supported by both ab initio calculations and an effective Hamiltonian analysis.The IRPD method is often used under the assumption that the spectrum of the weakly-bound complex can be taken as proxy for that of the bare ion.Seeing that this is clearly not the case for HCCH + and DCCD + and that similar effects were reported for HC 3 N + , 42 where it was found that the Ne attachment affected the RT splitting, we argue that spectra obtained with the IRPD method should not be taken as a proxy for the bare ion spectra for RT affected species.In these cases LOS provides an excellent universal alternative.
In contrast to most closed-shell species, the state of ab initio calculations for the simulation of the ro-vibrational spectra of Renner-Teller perturbed molecules is relatively poor, as shown here for DCCD + .It is therefore particularly important to obtain accurate experimental spectroscopic constants that are used to simulate the radio-or infrared spectra through which these species are detected.The newly launched James-Webb Space Telescope (JWST) provides mid-resolution observational infrared data at unprecedented sensitivity, and a means of detection for low-abundance symmetric species with no rotational spectrum like HCCH + .The spectroscopic parameters now derived and confirmed for the highly characteristic cis-bendings of DCCD + and HCCH + may be used for follow-up high-resolution studies or may even be sufficient for direct comparison to JWST data, as has recently been demonstrated for the hydrocarbon ion CH 3 + . 59,60This paper also highlights the complementary nature of PFI-ZEKE PES and infrared LOS, where PFI-ZEKE may provide rotationally resolved spectra and accurate spectroscopic constants and LOS information on expected line intensities and selection rules.The demonstration of LOS as a universal broadband infrared action spectroscopic method as shown here, opens the door for future studies of other astroand fundamentally relevant RT affected cations, such as linear HCCCH + , CCH + and HCCN + .

Fig. 2
Fig. 2 Mass selected HCCH + (a) and DCCD + (b) together with the mass spectra of the trap content after 600 ms trapping under Ne-tagging, Ar-tagging and leak-out conditions.Mass channels m/z +2 (i.e.m/z 28 for HCCH + and m/z 30 for DCCD + ) are a consequence of H 2 trace contamination in the Ne.Here the first quadrupole is used to mass select the ion of interest before entering the trap and the second quadrupole is used for scanning over the m/z range after trapping.

3 . 2 . 1
Discussion on measured and calculated intensities.Finally, we wish to draw attention to the low energy of the recorded cis-bendings of HCCH + and DCCD + .As described by Schmid et al., 28 the LOS signal should critically depend on the energy of the vibrational transition and may depend on the ion temperature as well.In principle, the kinetic energy release is divided among the two collision partners, here Ne (m = 20 u) and HCCH + (m = 26 u) or DCCD + (m = 28 u), resulting in a maximum energy gain of DE kin ¼ m Ne m Ne þ m ion E vib for the ion.

Fig. 6
Fig. 6 Comparison of the infrared spectra of the bare HCCH + ion taken with LOS and of Ne-HCCH + and Ar-HCCH + taken with IRPD.

Fig. 7
Fig. 7 Calculated counterpoise corrected interaction energies as a function of (a) Ne angle and (b) Ar angle with respect to the molecular axis.The symmetric A 0 state corresponds to configuration (p x ) 1 (p y ) 2 and the antisymmetric A 00 state to (p x ) 2 (p y ) 1 .

Table 1
Spectroscopic and experimental parameters used for the simulated LOS and LIR spectra of HCCH + shown in Fig.4.Unmarked parameters are taken from ref.21 a Parameters are fitted by hand.b Parameters are calculated ab initio on RCCSD(T)/aug-cc-pV5Z level of theory.

Table 2
Comparison of rotational band origins of cis-bending mode n 5 of DCCD + (in cm À1 )

Table 4
Spectroscopic parameters used for the simulated spectra of Ne-HCCH + and Ne-DCCD + shown in Fig.8.Unmarked parameters are taken from ref.21 aThis journal is © the Societies 2024