Inelastic neutron scattering spectrum of H2@C60 and its temperature dependence decoded using rigorous quantum calculations and a new selection rule

Targeting spectroscopic accuracy for dispersion bound systems from ab initio techniques: Translational eigenstates of Ne@C70

We investigate the endofullerene system Ne@C70 by constructing a three-dimensional Potential Energy Surface (PES) describing the translational motion of the Ne atom. This is constructed from electronic structure calculations from a plethora of methods, including MP2, SCS-MP2, SOS-MP2, RPA@PBE, and C(HF)-RPA, which were previously used for He@C60 in Panchagnula et al. [J. Chem. Phys. 160, 104303 (2024)], alongside B86bPBE-25X-XDM and B86bPBE-50X-XDM. The reduction in symmetry moving from C60 to C70 introduces a double well potential along the anisotropic direction, which forms a test of the sensitivity and effectiveness of the electronic structure methods. The nuclear Hamiltonian is diagonalized using a symmetrized double minimum basis set outlined in Panchagnula and Thom [J. Chem. Phys. 159, 164308 (2023)], with translational energies having error bars ±1 and ±2 cm-1. We find no consistency between electronic structure methods as they find a range of barrier heights and minima positions of the double well and different translational eigenspectra, which also differ from the Lennard-Jones (LJ) PES given in Mandziuk and Bačić [J. Chem. Phys. 101, 2126-2140 (1994)]. We find that generating effective LJ parameters for each electronic structure method cannot reproduce the full PES nor recreate the eigenstates, and this suggests that the LJ form of the PES, while simple, may not be best suited to describe these systems. Even though MP2 and RPA@PBE performed best for He@C60, due to the lack of concordance between all electronic structure methods, we require more experimental data in order to properly validate the choice.

Translational eigenstates of He@C60 from four-dimensional ab initio potential energy surfaces interpolated using Gaussian process regression

We investigate the endofullerene system 3He@C60 with a four-dimensional potential energy surface (PES) to include the three He translational degrees of freedom and C60 cage radius. We compare second order Møller-Plesset perturbation theory (MP2), spin component scaled-MP2, scaled opposite spin-MP2, random phase approximation (RPA)@Perdew, Burke, and Ernzerhof (PBE), and corrected Hartree-Fock-RPA to calibrate and gain confidence in the choice of electronic structure method. Due to the high cost of these calculations, the PES is interpolated using Gaussian Process Regression (GPR), owing to its effectiveness with sparse training data. The PES is split into a two-dimensional radial surface, to which corrections are applied to achieve an overall four-dimensional surface. The nuclear Hamiltonian is diagonalized to generate the in-cage translational/vibrational eigenstates. The degeneracy of the three-dimensional harmonic oscillator energies with principal quantum number n is lifted due to the anharmonicity in the radial potential. The (2l + 1)-fold degeneracy of the angular momentum states is also weakly lifted, due to the angular dependence in the potential. We calculate the fundamental frequency to range between 96 and 110 cm-1 depending on the electronic structure method used. Error bars of the eigenstate energies were calculated from the GPR and are on the order of ∼±1.5 cm-1. Wavefunctions are also compared by considering their overlap and Hellinger distance to the one-dimensional empirical potential. As with the energies, the two ab initio methods MP2 and RPA@PBE show the best agreement. While MP2 has better agreement than RPA@PBE, due to its higher computational efficiency and comparable performance, we recommend RPA as an alternative electronic structure method of choice to MP2 for these systems.

A nuclear configuration interaction approach to study nuclear spin effects: an application to ortho- and para-3 He2 @C60

We introduce a non-orthogonal configuration interaction approach to investigate nuclear quantum effects on energies and densities of confined fermionic nuclei. The Hamiltonian employed draws parallels between confined systems and many-electron atoms, where effective non-Coulombic potentials represent the interactions of the trapped particles. One advantage of this method is its generality, as it offers the potential to study the nuclear quantum effects of various confined species affected by effective isotropic or anisotropic potentials. As a first application, we analyze the quantum states of two 3 He atoms encapsulated in C60 . At the Hartree-Fock level, we observe the breaking of spin and spatial symmetries. To ensure wavefunctions with the correct symmetries, we mix the broken-symmetry Hartree-Fock states within the non-orthogonal configuration interaction expansion. Our proposed approach predicts singly and triply degenerate ground states for the singlet (para-3 He2 @C60 ) and triplet (ortho-3 He2 @C60 ) nuclear spin configurations, respectively. The ortho-3 He2 @C60 ground state is 5.69 cm-1 higher in energy than the para-3 He2 @C60 ground state. The nuclear densities obtained for these states exhibit the icosahedral symmetry of the C60 embedding potential. Importantly, our calculated energies for the lowest 85 states are in close agreement with perturbation theory results based on a harmonic oscillator plus rigid rotor model of 3 He2 @C60 .

Real-time hydrogen molecular dynamics satisfying the nuclear spin statistics of a quantum rotor

Apparent presence of the nuclear-spin species of a hydrogen molecule, para-hydrogen and ortho-hydrogen, associated with the quantum rotation is a manifestation of the nuclear quantum nature of hydrogen, governing not only molecular structures but also physical and chemical properties of hydrogen molecules. It has been a great challenge to observe and calculate real-time dynamics of such molecularized fermions. Here, we developed the non-empirical quantum molecular dynamics method that enables real-time molecular dynamics simulations of hydrogen molecules satisfying the nuclear spin statistics of the quantum rotor. While reproducing the species-dependent quantum rotational energy, population ratio, specific heat, and H-H bond length and frequency, we found that their translational, orientational and vibrational dynamics becomes accelerated with the higher rotational excitation, concluding that the nuclear quantum rotation stemmed from the nuclear spin statistics can induce various kinds of dynamics and reactions intrinsic to each hydrogen species.

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.

Experimental determination of the interaction potential between a helium atom and the interior surface of a C60 fullerene molecule

The interactions between atoms and molecules may be described by a potential energy function of the nuclear coordinates. Nonbonded interactions between neutral atoms or molecules are dominated by repulsive forces at a short range and attractive dispersion forces at a medium range. Experimental data on the detailed interaction potentials for nonbonded interatomic and intermolecular forces are scarce. Here, we use terahertz spectroscopy and inelastic neutron scattering to determine the potential energy function for the nonbonded interaction between single He atoms and encapsulating C₆₀ fullerene cages in the helium endofullerenes ³He@C₆₀ and ⁴He@C₆₀, synthesized by molecular surgery techniques. The experimentally derived potential is compared to estimates from quantum chemistry calculations and from sums of empirical two-body potentials.

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.

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).

General Selection Rule in the Inelastic Neutron Scattering Spectroscopy of a Diatomic Molecule Confined Inside a Near-Spherical Nanocavity

Knowledge of the relevant selection rules is crucial for the accurate interpretation of experimental spectra in general. There has been a consensus for a long time that the incoherent inelastic neutron scattering (INS) spectroscopy of the vibrations of discrete molecular compunds is free from any selection rules. We contradict this widely held view by presenting an analytical derivation of the general selection rule for the INS spectroscopy of a hydrogen molecule inside a near-spherical nanocavity. It defines all forbidden transitions, originating in a range of initial translation-rotation (TR) states, ground and excited, of the caged para- and ortho-H2, as well as HD, that are unobservable in the INS spectra. These predictions are amenable to experimental verification. In addition, we demonstrate that the general selection rule applies to the INS spectroscopy of any diatomic molecule in a nanocavity with near-spherical symmetry, which exhibits strong TR coupling. Its existence strongly suggests that similar selection rules apply to the INS spectra of other molecular and supramolecular systems, and need to be identified.

Confirming a predicted selection rule in inelastic neutron scattering spectroscopy: the quantum translator-rotator H2 entrapped inside C60

We report an inelastic neutron scattering (INS) study of a H2 molecule encapsulated inside the fullerene C60 which confirms the recently predicted selection rule, the first to be established for the INS spectroscopy of aperiodic, discrete molecular compounds. Several transitions from the ground state of para-H2 to certain excited translation-rotation states, forbidden according to the selection rule, are systematically absent from the INS spectra, thus validating the selection rule with a high degree of confidence. Its confirmation sets a precedent, as it runs counter to the widely held view that the INS spectroscopy of molecular compounds is not subject to any selection rules.

Symmetry-breaking in the endofullerene H2O@C60 revealed in the quantum dynamics of ortho and para-water: a neutron scattering investigation

The splitting of the ortho-H₂O ground state is clearly revealed by inelastic neutron scattering.