Communication: The H2@C60 inelastic neutron scattering selection rule: Expanded and explained

H2O inside the fullerene C60: Inelastic neutron scattering spectrum from rigorous quantum calculations

We present a methodology that, for the first time, allows rigorous quantum calculation of the inelastic neutron scattering (INS) spectra of a triatomic molecule in a nanoscale cavity, in this case, H₂O inside the fullerene C₆₀. Both moieties are taken to be rigid. Our treatment incorporates the quantum six-dimensional translation-rotation (TR) wave functions of the encapsulated H₂O, which serve as the spatial parts of the initial and final states of the INS transitions. As a result, the simulated INS spectra reflect the coupled TR dynamics of the nanoconfined guest molecule. They also exhibit the features arising from symmetry breaking observed for solid H₂O@C₆₀ at low temperatures. Utilizing this methodology, we compute the INS spectra of H₂O@C₆₀ for two incident neutron wavelengths and compare them with the corresponding experimental spectra. Good overall agreement is found, and the calculated spectra provide valuable additional insights.

The Endofullerene HF@C60: Inelastic Neutron Scattering Spectra from Quantum Simulations and Experiment, Validity of the Selection Rule, and Symmetry Breaking

Accurate quantum simulations of the low-temperature inelastic neutron scattering (INS) spectra of HF@C₆₀ are reported for two incident neutron wavelengths. They are distinguished by the rigorous inclusion of symmetry-breaking effects in the treatment and having the spectra computed with HF as the guest, rather than H₂ or HD, as in the past work. The results demonstrate that the precedent-setting INS selection rule, originally derived for H₂ and HD in near-spherical nanocavities, applies also to HF@C₆₀, despite the large mass asymmetry of HF and the strongly mixed character of its translation-rotation eigenstates. This lends crucial support to the theoretical prediction made earlier that the INS selection rule is valid for any diatomic molecule in near-spherical nanoconfinement. The selection rule remains valid in the presence of symmetry breaking but is modified slightly in an interesting way. Comparison is made with the recently published experimental INS spectrum of HF@C₆₀. The agreement is very good, apart from one peak for which our calculations suggest a reassignment. This reassignment is consistent with the measured INS spectrum presented in this work, which covers an extended energy range.

Interactions between a water molecule and C60 in the endohedral fullerene H2O@C60

A water molecule encapsulated inside a C60 fullerene cage behaves almost like an asymmetric top rotor, as would be expected of an isolated water molecule. However, inelastic neutron scattering (INS) experiments show evidence of interactions between the water molecule and its environment [Goh et al., Phys. Chem. Chem. Phys., 2014, 16, 21330]. In particular, a resolved splitting of the 101 rotational level into a singlet and a doublet indicates that the water molecule experiences an environment of lower symmetry than the icosahedral symmetry of a C60 cage. Recent calculations have shown that the splitting can be explained in terms of electrostatic quadrupolar interactions between the water molecule and the electron clouds of nearest-neighbour C60 molecules, which results in an effective environment of S6 symmetry [Felker et al., Phys. Chem. Chem. Phys., 2017, 19, 31274 and Bačić et al., Faraday Discussions, 2018, 212, 547-567]. We use symmetry arguments to obtain a simple algebraic expression, expressed in terms of a linear combination of products of translational and rotational basis functions, that describes the effect on a water molecule of any potential of S6 symmetry. We show that we can reproduce the results of the electrostatic interaction model up to ≈12 meV in terms of two unknown parameters only. The resulting potential is in a form that can readily be used in future calculations, without needing to use density functional theory (DFT) for example. Adjusting parameters in our potential would help identify whether other symmetry-lowering interactions are also present if experimental results that resolve splittings in higher-energy rotational levels are obtained in the future. As another application of our model, we show that the results of DFT calculations of the variation in energy as a water molecule moves inside the cage of an isolated C60 molecule, where the water molecule experiences an environment of icosahedral symmetry, can also be reproduced using our model.

Effects of symmetry breaking on the translation-rotation eigenstates of H2, HF, and H2O inside the fullerene C60

Splittings of the translation-rotation (TR) eigenstates of the solid light-molecule endofullerenes M@C60 (M = H2, H2O, HF) attributed to the symmetry breaking have been observed in the infrared (IR) and inelastic neutron scattering spectra of these species in the past couple of years. In a recent paper [Felker et al., Phys. Chem. Chem. Phys., 2017, 19, 31274], we established that the electrostatic, quadrupolar interaction between the guest molecule M and the twelve nearest-neighbor C60 cages of the solid is the main source of the symmetry breaking. The splittings of the three-fold degenerate ground states of the endohedral ortho-H2, ortho-H2O and the j = 1 level of HF calculated using this model were found to be in excellent agreement with the experimental results. Utilizing the same electrostatic model, this theoretical study investigates the effects of the symmetry breaking on the excited TR eigenstates of the three species, and how they manifest in their simulated low-temperature (5-6 K) near-IR (NIR) and far-IR (FIR) spectra. The TR eigenstates are calculated variationally for both the major P and minor H crystal orientations. For the H orientation, the calculated splittings of all of the TR levels of these species are less than 0.1 cm-1. For the dominant P orientation, the splittings vary strongly depending on the character of the excitations involved. In all of the species, the splittings of the higher rotationally excited levels are comparable in magnitude to those for the j = 1 levels. For the levels corresponding to purely translational excitations, the calculated splittings are about an order of magnitude smaller than those of the purely rotational eigenstates. Based on the computed TR eigenstates, the low-temperature NIR (for M = H2) and FIR (for M = HF and H2O) spectra are simulated for both the P and H orientations, and also combined as their weighted sum (0.15H + 0.85P). The weighted sum spectra computed for M = H2 and HF match quantitatively the corresponding measured spectra, while for M = H2O, the weighted sum FIR spectrum predicts features that can potentially be observed experimentally.

Perspective: Accurate treatment of the quantum dynamics of light molecules inside fullerene cages: Translation-rotation states, spectroscopy, and symmetry breaking

In this perspective, I review the current status of the theoretical investigations of the quantum translation-rotation (TR) dynamics and spectroscopy of light molecules encapsulated inside fullerenes, mostly C₆₀ and C₇₀. The methodologies developed in the past decade allow accurate quantum calculations of the TR eigenstates of one and two nanoconfined molecules and have led to deep insights into the nature of the underlying dynamics. Combining these bound-state methodologies with the formalism of inelastic neutron scattering (INS) has resulted in the novel and powerful approach for the quantum calculation of the INS spectra of a diatomic molecule in a nanocavity with an arbitrary geometry. These simulations have not only become indispensable for the interpretation and assignment of the experimental spectra but are also behind the surprising discovery of the INS selection rule for diatomics in near-spherical nanocavities. Promising directions for future research are discussed.

Explaining the symmetry breaking observed in the endofullerenes H2@C60, HF@C60, and H2O@C60

Symmetry breaking has been recently observed in the endofullerenes M@C₆₀ (M = H₂, HF, H₂O), manifesting in the splittings of the three-fold degenerate ground states of the endohedral ortho-H₂, ortho-H₂O and the j = 1 level of HF.

Experimental, theoretical and computational investigation of the inelastic neutron scattering spectrum of a homonuclear diatomic molecule in a nearly spherical trap: H2@C60

In this paper we report a methodology for calculating the inelastic neutron scattering spectrum of homonuclear diatomic molecules confined within nano-cavities of spherical symmetry.

Symmetry-breaking in the H2@C60 endofullerene revealed by inelastic neutron scattering at low temperature

The fine structure of the rotational ground state of molecular ortho-hydrogen confined inside the fullerene cage C₆₀ is investigated by inelastic neutron scattering (INS).