Matrix site effects on vibrational frequencies of HXeCCH, HXeBr, and HXeI: a hybrid quantum-classical simulation

The matrix shifts of the H–Xe stretching frequency of noble-gas hydrides, HXeCCH, HXeBr, and HXeI in various noble-gas matrices (in Ne, Ar, Kr, and Xe matrices) are investigated via the hybrid quantum-classical simulations. The order of the H–Xe stretching frequencies is found to be n (gas) o n (Ne) o n (Xe) o n (Kr) o n (Ar) for HXeCCH and HXeBr, while it is n (gas) o n (Ne) o n (Xe) o n (Ar) o n (Kr) for HXeI. This order is anomalous with respect to the matrix dielectric constants, and the calculated results reproduce the experimentally observed shifts quite successfully. We also find that the matrix shifts from the gas-phase values are D n (HXeCCH) E D n (HXeCl) o D n (HXeBr) o D n (HXeI) in the same noble-gas matrix environments, which implies that the weakly bound molecules exhibit large matrix shifts. The local trapping site is analyzed in detail, and it is shown that a realistic modeling of the surrounding matrix environments is essential to describe the unusual matrix shifts accurately.


I. Introduction
The last two decades have witnessed a rebirth of interest in the field of noble-gas chemistry mainly due to the discovery of noblegas hydrides. [1][2][3][4][5][6][7] These compounds have the structural formula of HNgY, where Ng is a noble gas atom and Y is an electronegative fragment. The UV photolysis of the HY precursors in the noble-gas matrix produces these HNgY compounds, and a unique electronic character and also the formation process of the HNgY compounds has been the focus of intensive research. The preparation of the noble-gas compounds in large Xe clusters by UV photolysis [8][9][10][11] and also by fast electrons and X-ray irradiation in the Xe matrix 12,13 have been reported in the literature as well.
Since these noble-gas hydrides have a strong ion-pair character of the form (HNg) + Y À , the resulting large dipole moment and also the weak bonding of these compounds sometimes induce unique solvation effects. 7,[14][15][16] For example, the vibrational spectroscopy measurements of HXeCl in various noble-gas matrices exhibited an unusual vibrational shift of the H-Xe stretching frequency. The order of the vibrational frequency was observed to be n(Ne) o n(Xe) o n(Kr) o n(Ar), and it is a non-monotonous function of the matrix dielectric constant. 1,[17][18][19] These results are in contrast to the ''normal'' order observed for HCl: n(Xe) o n(Kr) o n(Ar) o n(Ne) o n(gas). 20,21 The vibrational spectroscopy measurements of HXeCCH and HXeBr in various noble-gas matrices also revealed the unusual vibrational shifts of the H-Xe stretching frequency in the order n(Ne) o n(Xe) o n(Kr) o n(Ar), 1,17,[22][23][24][25] which is the same as that observed for HXeCl.
There have been several attempts to explain these unique features of the matrix shifts. First, it has been reported that simple calculations of neither HXeY with the polarizable continuum model (PCM) nor single-atom complex HXeY-Ng can reproduce the order of the vibrational shifts. [24][25][26] Very recently, Kalinowski et al. performed quantum chemical calculations by MP4(SDQ) for HXeCl embedded in a single noble-gas layer cage and successfully reproduced the order of the experimentally observed vibrational shifts. 19 They discussed these calculated shifts in terms of the size available for HXeCl in the cages and its implications for the stresses acting on the embedded HXeCl compound. Curiously, the vibrational frequencies of single-layer calculations by Ne atoms were strongly red-shifted from that of the isolated HXeCl; the shifts were À50 cm À1 and À192 cm À1 for harmonic and anharmonic calculations, respectively, which are the largest among the single-layer cages of four Ng atoms (Ng = Ne, Ar, Kr, and Xe). This is in contrast to the general observations that the interaction of an embedded compound with Ne atoms is the weakest and that the vibrational frequency in a Ne matrix should be closest to the situation in the gas phase. 27 At almost the same time as the work of Kalinowski et al., 19 we performed the hybrid quantum-classical simulations on the vibrational shift of HXeCl, where the interaction energy between HXeCl and surrounding Ng atoms was modeled by the pairwise-additive form and the pair-interaction energies were determined by high-level quantum chemical calculations. 28 We demonstrated that the experimentally observed vibrational shifts are well reproduced by the simulations and that the local matrix morphology is an important factor affecting the vibrational shifts. The largest shifts were predicted for the Ar matrix by our simulations and at that time we were not aware of the experimental results of the Ar matrix, but the very recent experimental results are in accord with our prediction. 19 It is noteworthy that the two different approaches of Kalinowski et al. 19 and ours predict the same order of vibrational frequency for HXeCl embedded in four noble-gas matrices.
This paper presents an application of our simulation method to other interesting noble-gas hydrides, HXeCCH, HXeBr, and HXeI in order to provide further information about the matrix shifts. The H-Xe stretching frequencies of HXeCCH 22-25 and HXeBr 1,17,25 have recently been reported in four different Ng matrices, and, therefore, a comparison can be made with our simulations. The spectroscopic measurements of HXeI and its complexes with various molecules (i.e. HXeIÁ Á ÁHI, HXeIÁ Á ÁH 2 O) have been reported in the Xe matrix before, 1,29,30 and very recently the vibrational frequencies of HXeI have been measured in Ar and Kr matrices. 31 It has been shown that the vibrational frequency of HXeI exhibits the order n(Xe) o n(Ar) o n(Kr) and it is slightly different from those observed for HXeCCH, HXeCl, and HXeBr (the experimental results are summarized in Table 1). In that paper, we performed the simulations for HXeI in Ar, Kr, and Xe matrices, and the calculated order is also in accord with the experimental data. 31

II. Computational details
The vibrational spectra of HXeCCH, HXeBr, and HXeI in noblegas (Ng) matrices are computed using the hybrid quantumclassical simulations. The computational methodology is essentially the same as that employed in our previous studies; see ref. 28 and 32 for details. In its implementation, the total potential energy is represented by the pairwise-additive form, where Y = (CCH, Br, or I) and q are the normal coordinates of HXeY. R (i) is the coordinate of i-th Ng atoms (Ng = Ne, Ar, Kr, and Xe). V HXeY-Ng is the interaction potential energy between HXeY and the Ng atom, and V Ng-Ng is that of the Ng dimer. In the simulations, the normal coordinate corresponding to the H-Xe stretching motion, q(H-Xe), is only exploited to represent V HXeY , and the other normal coordinates are fixed to zero (denoted as DVR-1D hereafter). The inclusion of the doubly degenerate bending vibrations is only considered for the vibrational analysis of HXeY and HXeY-Ng (denoted as DVR-3D hereafter). In our previous work on HXeCl, 28 it has been shown that the inclusion of Xe-Cl stretching vibration shows a negligible effect on the H-Xe stretching frequency (B1 cm À1 ). Also, it has been shown that the inclusion of bending motion by DVR-3D reduces the H-Xe stretching frequency by 20-30 cm À1 in comparison to DVR-1D, but the frequency shifts from the gas-phase values are close to each other. In a recent paper, 31 we reported the simulation results of HXeI in Ar, Kr, and Xe matrices and the vibrational frequencies have been calculated using the DVR-3D method. In the present paper, we provide the simulation results using the DVR-1D method and also add the results of HXeI in the Ne matrix, and comparative investigation of matrix shifts of HXeCCH, HXeCl, HXeBr, and HXeI in four different Ng matrices is carried out. In obtaining potential energy surfaces of V HXeY and also pair interaction energies of V HXeY-Ng and V Ng-Ng , the electronic structure calculations employing the coupled-cluster singles and doubles including the perturbative contributions of connected triple excitations [CCSD(T)] method are performed. For Kr, Xe, Br, and I atoms, the relativistic pseudopotentials 33 are employed. The cc-pVQZ basis sets are used for H, C, Br, and I atoms, while the aug-cc-pVQZ basis sets are employed for Ng atoms. It is known that the equilibrium distances and binding energies of Ng dimers with these basis sets yield quite reliable values (see ref. 28 for details) and these basis sets were employed in our previous study on HXeCl. All ab initio calculations are carried out using the MOLPRO2010.1 package. 34 The interaction potential energies between HXeY and an Ng atom (V HXeY-Ng ) are obtained at each point of the normal coordinate of q(H-Xe). The position of the Ng atom is represented by polar coordinates (r, y), where r is the distance of Ng from the center-of-mass of HXeY and y is the polar angle measured from the molecular axis of HXeY. The numbers of grid points are 25 and 21 for r and y, respectively. These potential energies are tabulated for use in the subsequent Monte Carlo (MC) simulations, and in the MC runs the potential energies are a The shifts, n(X)-n(Ne) (X = Ar, Kr, and Xe), were shown in the parentheses. evaluated by the interpolation. The number of grid points for the normal coordinate of the H-Xe stretching vibration corresponds to the number of optimized basis set for representing the vibrational wavefunction in the framework of the potentialoptimized discrete variable representation (PO-DVR) method. 35 In this work, the optimized basis set is constructed from the solution of the sinc-DVR method 36 for HXeY in the gas phase, and the number of optimized basis sets is determined to be 8, where the convergence of the H-Xe fundamental frequency with respect to the number of optimized basis sets is confirmed from the preliminary calculations. The MC simulations are carried out under the periodic boundary conditions, and 500 Ng atoms are first placed in a cubic box with the fcc lattice structure. Then, the HXeY compound is inserted into this fcc lattice by removing a couple of Ng atoms. There are several possible sites at which Ng atom(s) can be substituted by HXeY. Three possible orientations in the fcc structures are considered for the insertion of HXeY: h110i (nearest neighbor), h100i (four-atomic window), and h111i (three-atomic windows). The number of substituted Ng atoms (N r ) is determined by a sensible choice based on the distance between the adjacent Ng atoms along each orientation. The simulation parameters are tabulated in Table 2, and the size of the unit cell has been adjusted to yield the density of experimental values under the saturated vapor conditions. 37 A cutoff distance of the half the length of the unit cell is applied to all interactions. After equilibrating the system, a total of 10 7 MC steps are taken to obtain statistical averages.
The stability of different sites is compared using the stabilization energy, which is defined as where E 0 is the vibrational ground-state energy of H-Xe stretching motion in the presence of surrounding Ng atoms that is obtained using the PO-DVR method and E 0 gas is the corresponding value in the gas phase. hi indicates the ensemble average over MC configurations, and m Ng is the chemical potential of Ng atoms. In order to obtain m Ng , the MC simulations of pure solid Ng atoms (500 atoms) have been performed and m Ng is obtained by dividing the ensemble average of the total potential energy by the number of Ng atoms in the system.
The vibrational spectra are obtained from the following simple form: where E 1 is the vibrational first excited-state energy of H-Xe stretching motion in the matrices, and r(R) represents the distribution function of the surrounding Ng atoms, which is generated by the MC simulation.

III. Results and discussion
A. Equilibrium structures and vibrational analysis of HXeCCH, HXeBr, and HXeI The geometry optimizations were performed for HXeCCH, HXeBr, and HXeI using the CCSD(T) method, and the vibrational frequency analysis followed it. The structural parameters at equilibrium geometries and the harmonic and fundamental frequencies of H-Xe stretching vibration are summarized in Tables 3 and 4, respectively (see also Fig. 1 for equilibrium structures). The results of HXeCl are also shown in the tables for reference. 28 As seen in the table, the inclusion of anharmonicity induces large red shifts of 85, 137, and 179 cm À1 for HXeCCH, HXeBr, and HXeI, respectively. The inclusion of the bending motions by DVR-3D further reduces the frequencies by 16, 28, and 44 cm À1 for HXeCCH, HXeBr, and HXeI, respectively, with respect to those by DVR-1D. The natural charge and dipole moments of these compounds are given in Table 5, where one can clearly see that the compound exhibits a strong ion-pair character of the form (HXe) + Y À , in particular, for HXeCl, HXeBr, and HXeI.

B. Equilibrium structures and vibrational analysis of
HXeCCH-Ng, HXeBr-Ng, and HXeI-Ng (Ng = Ne, Ar, Kr, Xe) In this subsection, the equilibrium structures, binding energies, and vibrational frequencies of HXeY-Ng (Y = CCH, Br, I) complexes are examined. Here, the internal coordinates of HXeY are fixed at the equilibrium values of the monomer (given in Table 3). In fact, the intramolecular structures of HXeY are negligibly perturbed upon the formation of HXeY-Ng complexes, and the changes of bond lengths of HXeY between the monomer and the complexes in the fully relaxed configurations are less than 0.003 Å for all complexes. As expected from the previous studies, 28,38 three minima have been located for all complexes. Two of them are on both sides of HXeY in collinear geometry (H-side and Y-side structures) and the third one is located between and almost equidistant from Xe and Y (bent structure); see Fig. 2 for equilibrium structures of HXeCCH-Xe. The structural parameters and binding energies of an Ng atom at these minima are given in Tables 6 and 7, respectively. The strong  binding site is the bent structure for all complexes, and in collinear geometry the Y-side structure exhibits a slightly stronger binding energy than the X-side structure, except for HXeCCH-Kr and HXeCCH-Xe. Note that these binding energies are several times larger than those of the Ng dimer, where they are 0.095, 0.274, 0.374, and 0.489 kcal mol À1 for Ne 2 , Ar 2 , Kr 2 , and Xe 2 , respectively. 28 Fig. 3 shows the two-dimensional contour plot of the potential energy surface of HXeCCH-Xe, where the internal coordinates of HXeCCH are fixed at the equilibrium values of the monomer. One can see that the potential energies at the two minima in the collinear geometry are very shallow with respect to the displacement of Xe in the direction perpendicular to the molecular axis. In Table 7, it is seen that for each Ng the binding energies to HXeY in the bent structure are very similar to each other. Also we find that for each Ng the binding energies in the H-side and Y-side structures are very close to each other for HXeCl-Ng, HXeBr-Ng, and HXeI-Ng, and it decreases slightly as the halogen atom becomes larger. The binding energies of HXeCl-Ng, HXeBr-Ng, and HXeI-Ng are generally larger than those of HXeCCH-Ng, and this can be explained by the differences in the dipole moment of monomers (see Table 5).
The frequency analysis is made at these minimum energy structures, and the results are shown in Table 8. Note that harmonic frequencies are calculated at fully relaxed configurations and that the fundamental frequencies are obtained by fixing the position of Ng atom at the respective minima. As seen in the table, the frequencies of the H-Xe stretching vibration are most strongly affected at the H-side structures, whereas the Y-side and bent structures have a negligible effect on the H-Xe stretching frequency. For the H-side structures, the blue shifts are observed and this shift increases with the Ng size. These blue shifts were previously explained by the enhanced charge separation (HXe) + Y À upon complex formation. 39 The inclusion of anharmonicity induces large red shifts with respect    to the harmonic frequencies, but it is also noted that the frequency shifts from the gas-phase values become larger in the anharmonic calculations. The inclusion of the bending motion by DVR-3D further reduces the frequencies, but the frequency shifts from the gas-phase values are very similar in DVR-3D and DVR-1D calculations with the difference of B2 cm À1 .
The order of the H-Xe stretching frequency is n(Ne) o n(Ar) o n(Kr) o n(Xe) for HXeY-Ng complexes, which is a monotonous function of the matrix dielectric constants, and it is in disagreement with the experimental results in the matrices. The PCM calculations 25 for HXeCCH and HXeBr also predicted the same order of n(gas) o n(Ne) o n(Ar) o n(Kr) o n(Xe), simply indicating that the H-Xe stretching frequency becomes higher in matrices with larger dielectric constants, but again, this PCM model cannot explain the experimental observations. As already demonstrated in our previous work, 28 these results indicate that a more realistic modeling of the matrix environments is also necessary for HXeCCH, HXeBr, and HXeI. Table 9 shows a summary of the MC simulations, where the peak positions of the vibrational spectra of H-Xe stretching motion are listed. As noted above, the DVR-1D method is employed for solving the vibrational Schrödinger equation in the MC runs. For each orientation of the substitution site, the MC simulations were initiated with different N r , and we show only the results for the case of lower stabilization energy. As seen in the table, the h110i site exhibits a lower stabilization energy in all cases, and therefore the h110i site would be energetically favored. For the simulations in the Ne matrix, we found that the configurations of Ne atoms around HXeY become rather structureless after equilibration, and therefore the classification by the orientation in the fcc structure is somewhat vague. On the other hand, the configurations of Ar, Kr, and Xe atoms in the matrices are highly structured. For the MC simulations starting from the h110i or h100i site, the HXeY compound does not change its orientation, and the lattice structure of Ar, Kr, and Xe atoms surrounding HXeY is almost maintained. However, when the simulations are initiated from the h111i site, in most cases the compound changes its orientation and the configurations of Ar, Kr, and Xe atoms become almost identical to those of the h100i site after equilibration.   In all cases shown in Table 9, blue shifts from the isolated gas-phase values are observed. The observed blue shifts are about 2-3 times larger than those estimated from the calculations of isolated HXeY-Ng complexes (see the fundamental frequencies by DVR-1D given in Table 8). The order of H-Xe stretching frequencies is n(Ne) o n(Xe) o n(Kr) o n(Ar) for HXeCCH and HXeBr in the h110i site, which is in accord with the experimental observations. 1,17,[22][23][24][25] The same order has been seen for HXeCl in our previous work. 28 The observed order is slightly different for HXeI in the h110i site, where it is n(Ne) o n(Xe) o n(Ar) o n(Kr), and this order is also in accord with the experiment. 1,31 It is noted that the order of frequencies is different between the h110i and h100i sites.

C. Hybrid quantum-classical simulations in the solid matrices
The contour plots of two-dimensional distribution functions of Ng (= Ar, Kr, Xe) atoms around HXeCCH and HXeBr are shown in Fig. 4 in the h110i site. Here the x-axis is taken as the molecular axis of HXeY and r is the radial distance from the x-axis, and the origin is set to the center-of-mass of HXeY. Due to the highly structured configurations of solid matrices, the positions of the Ng atoms are quite localized. In the same figure, the positions of Ng atoms in the three minimum energy structures of the isolated HXeY-Ng complexes (H-side, Y-side, and bent structures) are also plotted.
As already seen in the frequency analysis of the isolated HXeY-Ng complex, the blue shifts are observed in all minimum energy structures of HXeY-Ng, and since there are several Ng atoms surrounding HXeY, the calculated larger blue shifts in the matrix environments are naturally expected. Larger blue shifts are observed in Ar and Kr matrices compared with the Xe matrix, and this could partly be explained by the observations (Fig. 4) that the Ar or Kr atom in collinear geometry on the H-side is found to be slightly closer to the hydrogen atom than that found in the minimum energy structures of HXeY-Ng, in comparison to the situation in the Xe matrix. Similar observations have been made for HXeCl and HXeI in our previous studies. 28,31 As shown above in the frequency analysis of HXeY-Ng, the H-Xe stretching frequency is strongly affected by Ng atoms that are close to the hydrogen atom. In this situation, a displacement of the Ng atom toward the hydrogen atom results in a blue shift because the potential energy curve is quite repulsive when the Ng atom is located closer to the hydrogen atom.
In matrix environments, we find that four Ng atoms are located in the position of the bent structure of the HXeY-Ng    complex, forming a ring around the molecular axis. Since the bent structure exhibits the strongest binding energy to HXeY (see Fig. 3 and Table 7), it is plausible that HXeY adjusts its coordinate along the molecular axis to have these surrounding four Ng atoms at the position of the bent structure. As a result of this, an Ng atom in the collinear geometry on the H-side can be close to the hydrogen atom in the case of the Ar or Kr matrix, causing the large blue shifts. Of course, the Ng atom cannot be too close to the hydrogen atom since there is quite a repulsive potential between the Ng and hydrogen atoms. Therefore, the location of HXeY in matrix environments depends on the delicate balance among stabilization and destabilization arising from the interactions between HXeY and the surrounding Ng atoms. As another explanation of larger blue shifts in Ar or Kr matrix, which is somewhat related to the above arguments, the Ar cage would be too tight for HXeY, while the Xe cage would be rather relaxed due to the difference in the lattice constants. The Ne matrix exhibits the shortest lattice constant, but the Ne matrix is ''soft'' due to the weak interaction between Ne atoms, and therefore HXeY would be mildly surrounded by Ne atoms. This is reflected in the small stabilization energies of ca. À1 to 2 kcal mol À1 in the Ne matrix in comparison to those in the other matrices (from À9 to À11 kcal mol À1 ) and also in the small frequency shifts from the gas-phase values (from +10 to +30 cm À1 ). The stress acting on HXeY would be an important factor affecting the vibrational shifts, which was also suggested by Kalinowski et al. 19 We note in passing that the small stabilization energies in the Ne matrix (sometimes positive values) can be related to the rather structureless configurations of Ne atoms around HXeY. Also, as pointed out in the previous work, 31 it may be the reason for a failure to prepare HXeI in the Ne matrix experimentally. It is also interesting to note that the frequency shifts from the gas-phase values exhibit a trend Dn(HXeCCH) E Dn(HXeCl) o Dn(HXeBr) o Dn(HXeI) in the same Ng matrix environments. The binding energies of an Ng atom to HXeCl, HXeBr and HXeI are very similar to each other in the HXeY-Ng complex and also the stabilization energies are quite close to each other for these compounds in the same Ng matrix, but there are appreciable differences in the frequency shifts. This trend is in agreement with the experiments (see Table 1), and it is consistent with the chemical intuition that the matrix effects are stronger for more weakly bound molecules.
Finally, we comment on the reliability of our calculations. First, the potential energy surfaces of embedded compounds and also the pair interaction potentials are determined using the CCSD(T) method with the (aug-)cc-pVQZ basis sets, which is known to be quite accurate, but of course further improvement would be possible (e.g., taking the complete basis set limit). We speculate that discrepancies in the absolute values of frequencies between experiments and our calculations (50-150 cm À1 ) stem mainly from the quality of potentials. The large discrepancies were found for HXeI, and this might be related to the weakness of the chemical bond of this compound. Second, the higher-body interactions between Ng atoms should be examined carefully because it has been reported that the threebody contributions are not negligible in bulk solids ranging from 3% (in Ne) to 7% (in Xe) in cohesive energy. 40 The inclusion of the three-body contributions also increases the lattice constants by B1%, which might shift slightly the vibrational frequencies to the red. Lastly, the present MC simulations have been performed using the DVR-1D method for practical consideration of the computational cost. The vibrational frequencies of isolated HXeY-Ng complexes show that the H-Xe stretching frequencies decrease by B20 cm À1 with the inclusion of the bending motions by DVR-3D, but the shifts from the gas-phase values are very similar in DVR-1D and DVR-3D calculations with the difference of at most 4 cm À1 (see Table 8). In our previous study of HXeCl, 28 the MC simulations employing DVR-3D have predicted the frequencies of 1738, 1735, and 1710 cm À1 in Ar, Kr, and Xe matrices, respectively, and therefore the shifts from the gas-phase values are +119, +116, and +91 cm À1 , which is close to the results obtained by DVR-1D shown above. For HXeI, the shifts from the gas-phase values in DVR-3D calculations were +226, +230, and +196 cm À1 in Ar, Kr, and Xe matrices, 31 respectively, and these shifts are slightly larger than those obtained from DVR-1D calculations (see Table 9). We expect that the differences in relative shifts between DVR-1D and DVR-3D are smaller for HXeCCH and HXeBr than HXeI considering the trend in the vibrational frequencies of these compounds, and the order of the vibrational frequency will almost certainly be unchanged.

IV. Conclusions and outlook
The hybrid quantum-classical simulations are performed to investigate the unusual vibrational spectral shifts of noble-gas hydrides, HXeCCH, HXeBr, and HXeI in matrix environments (in Ne, Ar, Kr, and Xe matrices). It is demonstrated that the H-Xe stretching frequencies are blue-shifted from the isolated gas-phase values in all matrix environments and that the amounts of blue shifts are a non-monotonous function of the matrix dielectric constants. The calculated results reproduce the experimentally observed shifts quite successfully.
Our calculations have clarified that the local matrix morphology is an important factor affecting the vibrational shifts and that a realistic modeling of the surrounding matrix environments is essential to reproduce the unusual vibrational shifts observed in the experiments. The orders of the H-Xe stretching frequencies are n(Ne) o n(Xe) o n(Kr) o n(Ar) for HXeCCH and HXeBr and n(Ne) o n(Xe) o n(Ar) o n(Kr) for HXeI, and therefore the order of the matrix shifts depends on the embedded compound. Also our calculations have shown that the matrix shifts from the gas-phase values exhibit the trend Dn(HXeCCH) E Dn(HXeCl) o Dn(HXeBr) o Dn(HXeI) in the same Ng matrix environments. This is consistent with the experimental results and also with the speculation that the stronger matrix effects are observed for more weakly bound molecules.
The theoretical investigations of the vibrational spectra in other matrices such as N 2 or CO 2 will further enhance our understanding of the matrix shifts. In N 2 or CO 2 matrices, the large matrix shifts were observed experimentally due to the dipole-quadrupole interactions between HXeY and matrix molecules. 41 The elucidation of these shifts is also an exciting future challenge.

Conflicts of interest
The authors declare no competing financial interest.