HD in C₆₀: Theoretical prediction of the inelastic neutron scattering spectrum and its temperature dependence

We report rigorous quantum calculations of the inelastic neutron scattering (INS) spectra of HD@C₆₀, over a range of temperatures from 0 to 240K and for two incident neutron wavelengths used in recent experimental investigations. The computations were performed using our newly developed methodology, which incorporates the coupled five-dimensional translation-rotation (T-R) eigenstates of the guest molecule as the initial and final states of the INS transitions, and yields highly detailed spectra. Depending on the incident neutron wavelength, the number of computed INS transitions varies from almost 500 to over 2000. The low-temperature INS spectra display the fingerprints of the coupling between the translational and rotational motions of the entrapped HD molecule, which is responsible for the characteristic splitting patterns of the T-R energy levels. INS transitions from the ground T-R state of HD to certain sublevels of excited T-R multiplets have zero intensity and are absent from the spectra. This surprising finding is explained by the new INS selection rule introduced here. The calculated spectra exhibit strong temperature dependence. As the temperature increases, numerous new peaks appear, arising from the transitions originating in excited T-R states which become populated. Our calculations show that the higher temperature features typically comprise two or more transitions close in energy and with similar intensities, interspersed with numerous other transitions whose intensities are negligible. This implies that accurately calculated energies and intensities of INS transitions which our methodology provides will be indispensable for reliable interpretation and assignment of the experimental spectra of HD@C₆₀ and related systems at higher temperatures.