Perfluorination of Aromatic Compounds Reinforce Their van der Waals Interactions with Rare Gases: The Rotational Spectrum of Pentafluoropyridine-Ne

The rotational spectrum of the pentafluoropyridine-Ne complex, generated in a supersonic jet, has been investigated using chirped-pulse microwave Fourier transform spectroscopy in the 2–8 GHz range. The spectra of the 20Ne and 22Ne species have been observed, and the rotational constants have been used to determine the structure of the complex. This structure, and those of the previously experimentally studied complexes benzene-Ne and pyridine-Ne, are an excellent benchmark for the theoretical calculations on these adducts. These complexes and hexafluorobenzene-Ne have been investigated at the CCSD/6-311++G(2d,p) level. The calculations reproduce the experimental structures well and show how the van der Waals complexes are stronger for the perfluorinated compound.


Introduction
Supersonic jet expansions have been widely used to generate intermolecular complexes of diverse nature, which are studied by spectroscopic methods to characterize the noncovalent forces responsible for their formation. Fourier transform microwave spectroscopy (FTMW) combined with supersonic jets with reinforced high resolution and sensitivity [1][2][3][4] has contributed to the study of challenging molecular systems, including weakly bound molecular complexes [5]. This technique makes possible the unambiguous discrimination between different isomers, including tautomers, conformers, or isotopomers. The spatial mass distribution of such species has unique spectroscopic constants and distinct individual rotational spectra, a key feature that makes FTMW techniques powerful tools [6] in cases where other techniques may struggle. Complexes generated in supersonic expansions are isolated species that can be characterized in absence of intermolecular forces existing in condensed phases. This provides an ideal environment to investigate their intrinsic properties, structure, and interactions, which can be benchmarked with the results of theoretical calculations [7].
Aggregates formed by a molecule and a rare gas (RG) atom, generally that used to drive the expansion, have been studied to characterize the van der Waals interactions due to dispersion forces [8]. These forces have an important role in the stabilization of biological systems and materials [9], and have been also studied in complexes or rare gases present in the upper atmosphere with oxygen molecules in excited electronic states [10]. Aromatic rings offer a relatively large surface for dispersion interaction forces, so these kinds of molecules have often been used to model the interactions with rare gas atoms. Complexes of the prototype rings, benzene (BZ), with He [11], Ne [12,13], Ar [13,14], Kr [15], and Xe [13], or pyridine (PY), with He [16], Ne [17][18][19][20], Ar [19,21,22], and Kr [21], have been studied. Complexes of RG atoms with other six-membered ring aromatic molecules such as pyrimidine [23,24], pyridazine [25,26], and fluorobenzene [27][28][29] or fluoropyridine [30] derivatives have also been characterized.
The structure of molecule-RG atom complexes has been investigated through the observation of the rotational spectra of different isotopologues, including those of RG atoms. Although those parameters are often affected by intermolecular vibrations, it has been possible to investigate the geometry of the heterodimers, which in the case of the benzene-RG complex corresponds to a C 6v symmetry where the RG atom lies along the C 6 benzene symmetry axis [13,14]. For pyridine-RG complexes, the RG atom is above the ring plane but is considerably shifted towards the N atom [16][17][18][19][20][21][22]. The rotation of the RG atom to explore both sides of the ring would give rise to two equivalent forms; however, with the exception of He [11,16], the barriers hindering this flipping motion are high enough to prevent the observation of tunneling effects. In many cases, under the assumption that the RG stretching motion can be isolated from bending vibrations, the dissociation energies of these complexes have been estimated from the ∆ J centrifugal distortion constants based on the pseudodiatomic approximation [17,20,[22][23][24][25]. Additionally, in some cases, information about the bending motions has been obtained [24,31].
The complexes of perfluorinated aromatic rings such as hexafluorobenzene (F 6 BZ) or pentafluoropyridine (F 5 PY) with water [32,33] ammonia [34], and formaldehyde [35] have led to important conclusions about the changes induced by perfluorination in the aromatic π-cloud. The aromatic π-cloud acts, in general, as a Lewis base in noncovalent interactions. Fluorine atoms withdraw the electron density of the π-cloud region, creating a π-hole [36] which acts as a Lewis acid [37][38][39][40] that can interact with lone-pair-bearing atoms and enhance anion-π interactions [39,41,42]. In fact, in the complexes of F 6 BZ or F 5 PY with water, ammonia, or formaldehyde, the O or N lone pairs of these small molecules point directly toward the center of the ring, showing the existence of lone-pair···π-hole interactions. These interactions show almost equal contributions of electrostatic and dispersion forces [34,35], but the changes in the structures from O-H···π, N-H···π in the BZ-H 2 O or BZ-NH 3 or N···C=O (n→ π*) in PY-CH 2 O to the H 2 O···π-hole, H 3 N···π-hole or H 2 C=O···π-hole in F 6 BZ and F 5 PY complexes can probably be attributed mainly to electrostatic changes. An investigation of the effects of perfluorination on the π-cloud and their implications on the dispersion forces contributing to the formation of complexes has not been conducted. In this paper we characterized the structure of the complex F 5 PY-Ne from the rotational spectrum. We compared the structure of this complex with previous experimental results on PY-Ne. Furthermore, a comparative study on the theoretically calculated energies and binding properties of BZ-Ne, PY-Ne, F 6 BZ-Ne and F 5 PY-Ne was carried out.

Rotational Spectrum
We recently studied the broadband rotational spectrum of the complex F 5 PY-CH 2 O [35], generated in a supersonic jet driven by an expansion of Ne. In the spectrum, the most intense lines were those of the monomers, followed by those of the complex F 5 PY-H 2 O. Several sets of rotational lines having a similar intensity to those of F 5 PY-CH 2 O lines remained unassigned (see Figure S1 [35]). Those lines show the quadrupole coupling hyperfine structure (hfs) of a species carrying a 14 N nucleus (see Figure 1). This hfs results from the interaction of the nuclear quadrupole moment of the 14 N nucleus (I = 1), eQ, with the electric field gradient, q, at the site of the N nucleus. The corresponding spectroscopic constants, χ αβ (α, β = a, b, c), are the elements of the nuclear quadrupole coupling tensor, set up in the inertial principal axis system representation, and are a direct measure of the electric field gradient at the nitrogen nucleus (χ αβ = eQq αβ ). FTMW spectroscopy has enough resolution to detect and resolve the individual quadrupole hyperfine components (see Figure 1). Once the lines of all known species were dropped from the spectrum, a search of possible species was conducted to identify the unknown spectra. Among the possible species, those corresponding to complexes of the molecules in the mixture with the carrier gas should be considered; in this case, the complex F 5 PY-Ne. Using the calculated rotational parameters for this complex, a prediction of the spectrum in the 2-8 GHz region was made. According to the predictions, this van der Waals complex has only the µ b type spectrum. The results of the theoretical calculations were valuable in the identification of rotational lines based on the frequency patterns and quadrupole coupling hyperfine structure. The quadrupole coupling constants are related to the electronic environment of this nucleus, and in this case, strongly depend on the orientation of the nitrogen atom in the principal inertial axis system. Figure 1 shows the comparison of the predicted spectrum for F 5 PY-Ne (upper blue trace), with the observed one (black lower trace) in the 5.0-5.6 GHz range. In this figure, it is illustrated how the R-branch transitions J + 1 0,J + 1 ← J 1,J and J + 1 1,J + 1 ← J 0,J or J + 1 J + 1,1 ← J J,0 and J + 1 J + 1,0 ← J J,1 formed well-defined doublets which were identified from frequency and intensity patterns for the parent species, and for the 22 Ne isotopologue with a natural abundance ratio 22 Ne/ 20 Ne = 1/3. In the excerpt, it can be seen how the comparison of the predicted and observed quadrupole coupling hyperfine structure patterns confirms the assignment. constants, χαβ (α, β = a, b, c), are the elements of the nuclear quadrupole coupling tensor, set up in the inertial principal axis system representation, and are a direct measure of the electric field gradient at the nitrogen nucleus (χαβ = eQqαβ). FTMW spectroscopy has enough resolution to detect and resolve the individual quadrupole hyperfine components (see Figure 1). Once the lines of all known species were dropped from the spectrum, a search of possible species was conducted to identify the unknown spectra. Among the possible species, those corresponding to complexes of the molecules in the mixture with the carrier gas should be considered; in this case, the complex F5PY-Ne. Using the calculated rotational parameters for this complex, a prediction of the spectrum in the 2-8 GHz region was made. According to the predictions, this van der Waals complex has only the μb type spectrum. The results of the theoretical calculations were valuable in the identification of rotational lines based on the frequency patterns and quadrupole coupling hyperfine structure. The quadrupole coupling constants are related to the electronic environment of this nucleus, and in this case, strongly depend on the orientation of the nitrogen atom in the principal inertial axis system. Figure 1 shows the comparison of the predicted spectrum for F5PY-Ne (upper blue trace), with the observed one (black lower trace) in the 5.0-5.6 GHz range. In this figure, it is illustrated how the R-branch transitions J + 10,J + 1 ← J1,J and J + 11,J + 1 ← J0,J or J + 1J + 1,1 ← JJ,0 and J + 1J + 1,0 ← JJ,1 formed well-defined doublets which were identified from frequency and intensity patterns for the parent species, and for the 22 Ne isotopologue with a natural abundance ratio 22 Ne/ 20 Ne = 1/3. In the excerpt, it can be seen how the comparison of the predicted and observed quadrupole coupling hyperfine structure patterns confirms the assignment. The assignments of the complete μb-type spectra were followed from successive fitting-prediction-measurement cycles. The analysis of the spectrum was conducted [43] using a Hamiltonian H = HR (A) + HQ, where HR (A) represents Watson's A-reduced semirigid rotor Hamiltonian in the I r representation [44] and HQ is the nuclear quadrupole coupling Figure 1. A section of the CP-FTMW spectrum in the 5.0-5.6 GHz region recorded from a supersonic jet expansion of F 5 PY diluted in Ne. The comparison with the prediction of the spectrum carried out from the rotational parameters predicted at CCSD/6-311++G(2d,p) (upper blue inverted trace) for F 5 PY-Ne complex helps to identify the patterns of the 4 0,4 ←3 1,3 , 4 1,4 ←3 0,3 and 3 3,1 ←2 2,0 , 3 3,0 ←2 2,1 doublets for both the parent (blue arrow) and 22 Ne (red arrow) isotopologues. The excerpt compares the observed and predicted 14 N quadrupole coupling hyperfine structure for the 3 3,0 ←2 2,1 transition where the components are identified by the quantum number F, taking values from (J + I) to (J − I) reflecting the coupling between the overall rotation moment (J) and 14 N nuclear spin moments (I = 1).
The assignments of the complete µ b -type spectra were followed from successive fittingprediction-measurement cycles. The analysis of the spectrum was conducted [43] using a Hamiltonian H = H R (A) + H Q , where H R (A) represents Watson's A-reduced semirigid rotor Hamiltonian in the I r representation [44] and H Q is the nuclear quadrupole coupling interaction term [45]. The determinable spectroscopic parameters are the rotational constants, the centrifugal distortion constants, and the elements of the nuclear quadrupole coupling tensor χ. Usually, for 14 N only the diagonal elements of the tensor are determined. The experimental constants determined for the F 5 PY-20 Ne and F 5 PY- 22 Ne isotopologues of the complex are collected in Table 1, where they are compared to those calculated at the CCSC/6-311++G(2d,p) level. The parameters calculated at other levels are listed in Table S1. The observed frequencies are collected in Tables S9 and S10.
are the planar moments of inertia, derived from the inertial moments P cc = (

Structure of the Complex
The comparison of the rotational constants, the corresponding inertial Iα and planar moments Pαα (α = a, b, or c), and the quadrupole coupling constants χαα of F5PY-Ne with those of the monomer, F5PY, provides experimental insights into the structure of the complex (see Table 2 and Figure 2). F5PY is a planar molecule with C2v symmetry. The C2 axis is coincident with the inertial axis b and the plane of the molecule is coincident with the

Structure of the Complex
The comparison of the rotational constants, the corresponding inertial I α and planar moments P αα (α = a, b, or c), and the quadrupole coupling constants χ αα of F 5 PY-Ne with those of the monomer, F 5 PY, provides experimental insights into the structure of the complex (see Table 2 and Figure 2). F 5 PY is a planar molecule with C 2v symmetry. The C 2 axis is coincident with the inertial axis b and the plane of the molecule is coincident with the σ ab inertial plane, while the symmetry plane perpendicular to the molecule can be identified as the σ bc . In correspondence, the planar moment P cc = (I a + I b − I c )/2 = ∑ i m i c 2 i , giving the mass extension out of the ab inertial plane, is nearly zero for F 5 PY (see Table 2 and Figure 2, where the principal axes of F 5 PY are labeled with a prime). The changes in the rotational parameters upon complexation give interesting clues about the structure and symmetry of the F 5 PY-Ne adduct. P cc is rather large for the cluster (see Table 2), indicating that Ne lies above (or below) the F 5 PY ring plane. The planar moment P aa , giving the mass extension out of the bc plane, is nearly the same for the monomer and the cluster (see Table 2). This proves that the Ne lies in the σ bc molecular plane, which is preserved as a symmetry plane upon formation of the complex, thus giving a C s symmetry of F 5 PY-Ne. Furthermore, the rotational constant C changes just about 1% indicating that the Ne atom lies close to the c inertial axis. It is worth noting that the P aa values show a slight decrease when going from F 5 PY to F 5 PY-20 Ne or F 5 PY- 22 Ne, revealing the contribution of large-amplitude intermolecular vibrations in the complex to this planar moment.  A, B, C), the moments of inertia (I a , I b , I c ), and the planar moments (P aa , P bb , P cc ) of F 5 PY and F 5 PY-Ne.  (22) a Taken from reference [35]. b Standard errors in parenthesis in units of the last digit.

Parameters
The χ aa , χ bb , and χ cc quadrupole coupling constants shown in Table 2 are the diagonal elements of the quadrupole coupling tensor set up in the principal inertial axis system. For F 5 PY, due to the C 2v symmetry, the principal quadrupole coupling axis and inertial axis systems are coincident (i.e., a→y, b→z, c→x) [35,46], so the quadrupole coupling tensor is diagonal. The conservation of the σ bc plane as a symmetry plane in the complex implies that the coincidence of the a inertial and y quadrupole coupling axes is also preserved in the complex. This can be deduced from the near equality of the χ aa quadrupole coupling constant in F 5 PY and F 5 PY-Ne (see Table 2). The orientation of the b and c inertial axes of the aggregate results from the rotation of the monomer axes by an angle θ bb' as defined in Figure 2. Using a value of θ bb' = 11.0 • to transform the experimental quadrupole coupling tensor constants of the monomer into those expected for the adduct, it gives χ bb = −3.737 MHz and χ cc = 1.771 MHz, values, very close to the experimental values of F 5 PY-20 Ne (see Table 2). An angle of 13.1 • gives χ bb = −3.648 MHz and χ cc = 1.682 MHz, close to the values of F 5 PY- 22 Ne. This evidences that the electric field gradient, and thus the electronic environment around the N atom, is not altered upon formation of the complex.
The degree of agreement between the experimental and calculated inertial data is better appreciated if we take the values of the planar moments of inertia P αα . The three translational motions of the isolated Ne are replaced by three low-energy large-amplitude vibrational modes upon formation of the complex. Due to these van der Waals motions, which may add vibration-rotation interaction contributions to the moments of inertia, the usual methods for structure determination may supply poor results. One of these vibrations is a stretching ν s , and the other two, ν ra and ν rb , can be associated to rotations of Ne around F 5 PY. For small displacements around the equilibrium structure these are described as nearly degenerate displacements of Ne in trajectories parallel to a (ν ra ) and b (ν rb ) axes. The frequencies calculated at B3LYP-D3BJ/6-311++G(2d,p) are ν s = 53.1 cm −1 , ν ra = 22.2 cm −1 , and ν rb = 24.6 cm −1 . The frequencies calculated with the MP2 method and the same basis set are ν s = 48.1 cm −1 , ν ra = 18.4 cm −1 and ν rb = 18.5 cm −1 . It can be seen in Table 1, Table 2, and Table S1 that the planar moments P aa and P bb are well-described theoretically in correspondence with the good description of F 5 PY rotational parameters at the levels of calculation used (see Table S2) [35]. The main discrepancies come from the planar moment P cc , which is highly dependent on the position of the Ne atom relative to the ring. In this case, the best description is made at the CCSD level. Part of these discrepancies may arise from the intermolecular vibrational contributions, in particular stretching vibrations, which may result in larger P cc experimental values corresponding with average ground vibrational state Ne-ring distances larger than those corresponding to equilibrium. As observed previously [35], the best agreement for the quadrupole coupling constants is obtained at MP2/6-311++G(2d,p) level.
The structure of F 5 PY-Ne has been further investigated from the analysis of the rotational parameters of the two observed isotopologues, by using the substitution (r s [47,48]), effective (r 0 [49]), and mass dependence (r m [49,50]) methods, as illustrated in Figure 2. The substitution, r s , method compares the inertial moments of the parent and monosubstituted isotopologues through the Kraitchman equations [47] to give the absolute values of the substituted atom coordinates-in this case, Ne-in the principal inertial axis system of the parent molecule. The signs of the coordinates can be assigned from any reasonable geometry as effective (r 0 ) or theoretical (r e ) structures. This experimental approach has some limitations in locating atoms close to the principal inertial axes, or for light atoms such as hydrogen. The r s coordinates of Ne in the complex are given in Table 3, where these are compared with those from the r 0 method and theoretical CCSD calculations. A total or partial effective r0 structure can be obtained when bond distances and angles are obtained from a least-squares fit of the inertial moments of all observed isotopologues. For F5PY-Ne, we have assumed that the structure of F5PY does not change upon complexation, so the previously determined r0 structure [35] was kept fixed to determine Figure 3. A sketch where the CCSD/6-311++G(2d,p) van der Waals predicted structure of F 5 PY-Ne (blue) is compared with that of PY-Ne (red) [20]. Table 3. Principal inertial axis coordinates for the neon atom of pentafluoropyridine···Ne complex. The table compares the substitution (r s ) and effective (r 0 ) coordinates with those (r e ) calculated at the CCSD/6-311++G(2d,p) (CCSD) computational methods. The table also compares the effective bonding parameters determined from a least-squares fit of the observed rotational constants, and the angle θ of rotation between the principal inertial axis systems of F 5 [48]. c Distance from Ne to the ring-centroid (rc, see Figure 3). d Distance from Ne to the F 5 PY center of mass (rcm, see Figure 3). e angle between the line from Ne to the ring center of mass (rcm) and a line perpendicular to the ring. f see Figure 2 for definition. g Differences between the experimental constants and those calculated from the determined r 0 structure.
A total or partial effective r 0 structure can be obtained when bond distances and angles are obtained from a least-squares fit of the inertial moments of all observed isotopologues. For F 5 PY-Ne, we have assumed that the structure of F 5 PY does not change upon complexation, so the previously determined r 0 structure [35] was kept fixed to determine the r(N-Ne) distance and ∠C 4 -N-Ne angle. The results are given in Table 3 and summarized in Figure 2, where it can be seen that the r 0 value of the angle θ is in good agreement with that calculated from the quadrupole coupling parameters, if we take into account that these may be affected by the effects of intermolecular vibrations. The comparison of the experimental structures with those derived theoretically at the CCSD/6-311++G(2d,p) level shown in Table 3, or with other methods (see Table S3) reflect the same degree of agreement pointed out when comparing the experimental values of the planar moments (see above); the main discrepancies come from the distance between the Ne atom and the ring plane. The CCSD method gives the closest results and nicely reproduces the r s values for F 5 PY-Ne. The CCSD predictions could then be considered reasonable, taking into account that the experimental r 0 data might be more affected by intermolecular vibration contributions, giving rise to larger distances.

Comparison with Related Complexes
The experimental and CCSD geometry of F 5 PY-Ne is compared in Table 4 with those of the three related complexes, PY-Ne, F 6 BZ-Ne, and BZ-Ne. The geometrical parameters given are the distance R from Ne to the ring center of mass, and the angle ϕ between the line from Ne to the ring center of mass and a line perpendicular to the ring plane ( Figure 3). These parameters allow a direct comparison with the geometries of two of these complexes, PY-Ne and BZ-Ne, already described experimentally. In all cases, the computational results reproduce well the r s values for F 5 PY-Ne and the equilibrium r e values derived from fitting the rotational parameters to Lennard-Jones type potential energy functions for the intramolecular stretching motions for PY-Ne [17] and BZ-Ne [13]. In the case of the r 0 values, the experimental distances are always larger than the predicted equilibrium distances, and this fact may reflect the contribution of large-amplitude intermolecular vibrations, especially the stretching vibrations to the ground vibrational state parameters. The longest predicted intermolecular R distances are obtained for the complexes of the hydrogenated aromatic rings (PY-Ne and BZ-Ne). The perfluorination of the aromatic systems decreases the intermolecular distance around 0.1 Å. The changes in the structures of PY-Ne and F 5 PY-Ne are sketched in Figure 3. In both, the RG atom is shifted towards the N atom, but in F 5 PY-Ne it is closer to the ring plane. Table 4. Intermolecular distance between the center of mass of the aromatic system and the Neon atom (Å), R, and the deviation of the perpendicular axes ( • ), ϕ. The values from experimental data sources are provided in parenthesis.

Binding Energies
The calculated binding energies at the CCSD(T) computational level and the DFT-SAPT contributions have been gathered in Table 5. As in the case of the intermolecular distances, the binding energies for the neon complexes of the hydrogenated aromatic systems are smaller than the perfluorinated ones by approx. 2 kJ mol −1 . The DFT-SAPT partition shows that in all cases, dispersion is the most important contribution, with values between −3.1 and −3.9 kJ mol −1 . The values of the dispersion term increase with the molecular weight of the aromatic system. A comparison of these values with those obtained with the empirical D3(BJ) method provides an excellent linear correlation (R 2 = 0.99, see Figure S2). The second most important term is the repulsive exchange that also increases with the molecular weight of the aromatic system from 2.2 kJ mol −1 in PY-Ne to 3.3 kJ mol −1 in Bz-Ne. The other three attractive terms have small contributions with the following trend in absolute value in each complex: electrostatic > δHF > Induction. Table 5. Bonding energies (BE in kJmol −1 ) calculated at the CCSD(T) level and DFT-SAPT contributions for the complexes PY-Ne, BZ-Ne, F 5 PY-Ne, and F 6 BZ-Ne.

BE[CCSD(T)]
Elect. The electron density analysis of the complexes shows the presence of three intermolecular BCPs in the F 5 PY-Ne complex (see Figure 4), one in the PY-Ne and six in the BZ-Ne and F 6 BZ-Ne ones (see the molecular graphs in Tables S5-S8). The values of the electron density of all the intermolecular BCPs are very close to what is considered to be the van der Waals limit in QTAIM theory (0.002 au). The largest value is obtained in the F 6 BZ-Ne complex (0.0027 au) and the smallest is obtained in the PY-Ne one (0.0022 au). Thus, the electron density values also are in concordance of very weak interactions between the two molecules. The NCIPlot analysis shows a region of weak interaction between the neon atom and the aromatic system (see Figure 4 for the F 5 PY-Ne complex).

Experimental Methods
A commercial sample of F5PY was used. The spectrum of the F5PY-Ne complex was created in a supersonic jet of Ne and investigated using a chirped-pulse Fourier transform microwave spectrometer (CP-FTMW) [66], described elsewhere [67], and a narrowband Fabry-Perot Fourier transform microwave spectrometer (FP-FTMW) [68,69]. F5PY was kept at room temperature in a deposit inserted in the gas line close to the solenoid valve in both spectrometers. In the CP-FTMW instrument, which covers the 2-8 GHz frequency range, the spectra were recorded in steps of 2 GHz. The carrier gas was Ne at backing pressures of about 2 bar, expanding through a 0.8 mm nozzle in pulses of 700 μs duration. Chirp pulses of 4 μs were created by an arbitrary waveform generator and amplified to 20 W. The polarization signal was radiated from a horn antenna in a direction perpendicular to that of the expanding gas. A molecular transient emission spanning 40 μs was then detected through a second horn, recorded with a digital oscilloscope, and Fourier-transformed to the frequency domain. The accuracy of frequency measurements was higher than 10 kHz. In the FP-FTMW instrument, operated in the 5-13 GHz frequency range, Ne was also used at stagnation pressures ranging up to 2 bar, expanding in pulses of about 800 μs through a 0.8 mm nozzle. Short (typ. 0.3 μs, 10-300 mW) microwave pulses were used for polarization purposes. Typically, a ca. 400 μs-length time-domain spectrum was recorded at 40-100 ns intervals and converted to the frequency domain by a fast Fourier

Experimental Methods
A commercial sample of F 5 PY was used. The spectrum of the F 5 PY-Ne complex was created in a supersonic jet of Ne and investigated using a chirped-pulse Fourier transform microwave spectrometer (CP-FTMW) [66], described elsewhere [67], and a narrowband Fabry-Perot Fourier transform microwave spectrometer (FP-FTMW) [68,69]. F 5 PY was kept at room temperature in a deposit inserted in the gas line close to the solenoid valve in both spectrometers. In the CP-FTMW instrument, which covers the 2-8 GHz frequency range, the spectra were recorded in steps of 2 GHz. The carrier gas was Ne at backing pressures of about 2 bar, expanding through a 0.8 mm nozzle in pulses of 700 µs duration. Chirp pulses of 4 µs were created by an arbitrary waveform generator and amplified to 20 W. The polarization signal was radiated from a horn antenna in a direction perpendicular to that of the expanding gas. A molecular transient emission spanning 40 µs was then detected through a second horn, recorded with a digital oscilloscope, and Fourier-transformed to the frequency domain. The accuracy of frequency measurements was higher than 10 kHz. In the FP-FTMW instrument, operated in the 5-13 GHz frequency range, Ne was also used at stagnation pressures ranging up to 2 bar, expanding in pulses of about 800 µs through a 0.8 mm nozzle. Short (typ. 0.3 µs, 10-300 mW) microwave pulses were used for polarization purposes. Typically, a ca. 400 µs-length time-domain spectrum was recorded at 40-100 ns intervals and converted to the frequency domain by a fast Fourier transformation. Due to the collinear arrangement of the jet and resonator axis, each rotational transition split into two Doppler components, so the resonant frequencies were taken as the arithmetic mean of both components. Frequency accuracy was lower than 3 kHz. The rotational spectra of the 22 Ne isotopologue was observed in its natural abundance.

Conclusions
In this work, we studied the rotational spectrum of the complex F 5 PY-Ne to characterize its structure and properties. These data, together with those from the rotational spectra of the complexes of BZ [12,13] and PY [17][18][19][20] with Ne, were used as a benchmark for the theoretical data obtained on these complexes using CCSD/6-311++G(2d,p) level with excellent results. The calculations were extended to the complex F 6 BZ-Ne to investigate the effects of perfluorination of BZ and PY in the interaction of these compounds with Ne. The experimental and theoretical results show that the van der Waals distances are shorter by ca. 0.1 Å in the perfluorinated compounds, indicating the latter are strongly bound. The binding energies calculated for the complexes corroborate this conclusion, since the energies of BZ-Ne and PY-Ne (−2.9 kcal mol −1 ) are smaller than the values (−4.9 kcal mol −1 ) of perfluorinated compounds.