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

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.

Fullerenes containing water molecules: a study of reactive molecular dynamics simulations

A different technique was used to investigate fullerenes encapsulating a polar guest species. By reactive molecular dynamics simulations, three types of fullerenes were investigated on a gold surface: an empty C60, a single H2O molecule inside C60 (H2O@C60), and two water molecules inside C60 ((H2O)2@C60). Our findings revealed that despite the free movement of all fullerenes on gold surfaces, confined H2O molecules within the fullerenes result in a distinct pattern of motion in these systems. The (H2O)2@C60 complex had the highest displacement and average velocity, while C60 had the lowest displacement and average velocity. The symmetry of molecules and the polarity of water seem to be crucial in these cases. ReaxFF simulations showed that water molecules in an H2O molecule, H2O@C60, and (H2O)2@C60 have dipole moments of 1.76, 0.42, and 0.47 D, respectively. A combination of the non-polar C60 and polar water demonstrated a significant reduction in the dipole moment of H2O molecules due to encapsulation. The dipole moments of water molecules agreed with those in other studies, which can be useful in the development of biocompatible and high-efficiency nanocars.

An off-center endohedrally confined hydrogen molecule

In this study, we address the problem of a C₆₀ endohedrally confined hydrogen molecule through a configuration-interaction approach to electronic dynamics. Modeling the confinement by means of a combination of two Woods-Saxon potentials, we analyze the stability of the system as a function of the nuclei position through the behavior of the electronic spectrum. After studying the convergence of two different partial wave expansions, one related to the molecular Coulomb centers and the other related to the off-centering of the C₆₀ well, we found that the second approach provides a more accurate description of the system. Furthermore, we observed that the inter-atomic distance changes based on the position of the atoms inside the cavity. Thus, to obtain the most favourable energetic configuration for the molecule, it should be positioned inside the cavity next to the structure, where its size decreases.

Unraveling the Origin of Symmetry Breaking in H2 O@C60 Endofullerene Through Quantum Computations

We explore the origin of the anomalous splitting of the 1₀₁ levels reported experimentally for the H₂O@C₆₀ endofullerene, in order to give some insight about the physical interpretations of the symmetry breaking observed. We performed fully‐coupled quantum computations within the multiconfiguration time‐dependent Hartree approach employing a rigorous procedure to handle such computationally challenging problems. We introduce two competing physical models, and discuss the observed unconventional quantum patterns in terms of anisotropy in the interfullerene interactions, caused by the change in the off‐center position of the encapsulated water molecules inside the cage or the uniaxial C₆₀‐cage distortion, arising from noncovalent bonding upon water's encapsulation, or exohedral fullerene perturbations. Our results show that both scenarios could reproduce the experimentally observed rotational degeneracy pattern, although quantitative agreement with the available experimental rotational levels splitting value has been achieved by the model that considers an uniaxial elongation of the C₆₀‐cage. Such finding supports that the observed symmetry breaking could be mainly caused by the distortion of the fullerene cage. However, as nuclear quantum treatments rely on the underlying interactions, a decisive conclusion hinges on the availability of their improved description, taken into account both endofullerene and exohedral environments, from forthcoming highly demanding electronic structure many‐body interaction studies.

Noncovalently bound molecular complexes beyond diatom–diatom systems: full-dimensional, fully coupled quantum calculations of rovibrational states

The methodological advances made in recent years have significantly extended the range and dimensionality of noncovalently bound molecular complexes for which full-dimensional quantum calculations of their rovibrational states are feasible.

Encapsulation of a Water Molecule inside C60 Fullerene: The Impact of Confinement on Quantum Features

We introduce an efficient quantum fully coupled computational scheme within the multiconfiguration time-dependent Hartree (MCTDH) approach to handle the otherwise extremely costly computations of translational–rotational–vibrational states and energies of light-molecule endofullenes. Quantum calculations on energy levels are reported for a water molecule inside C₆₀ fullerene by means of such a systematic approach that includes all nine degrees of freedom of H₂O@C₆₀ and does not consider restrictions above them. The potential energy operator is represented as a sum of natural potentials employing the n-mode expansion, along with the exact kinetic energy operator, by introducing a set of Radau internal coordinates for the H₂O molecule. On the basis of the present rigorous computations, various aspects of the quantized intermolecular dynamics upon confinement of H₂O@C₆₀ are discussed, such as the rotational energy level splitting and the significant frequency shifts of the encapsulated water molecule vibrations. The impact of water encapsulation on quantum features is explored, and insights into the nature of the underlying forces are provided, highlighting the importance of a reliable first-principles description of the guest–host interactions.

Exclusion principle repulsion effects on the covalent bond beyond the Born-Oppenheimer approximation

The changes in the covalent bond of the hydrogen molecule limited in space by a spherical hard boundary are studied. The sphere is moved along an axis parallel or orthogonal to the molecular axis. The diffusion Monte Carlo approach is used to solve the Schrödinger equation with the relevant boundary conditions and to evaluate the changes in the bond energy versus the location of the sphere. The vertical and lateral quantum forces exerted on the sphere are evaluated by calculating the energy derivative versus the distances to the sphere. The results show that the quantum forces present an important dependence on the distance and vanish rapidly as the separation between the sphere and the molecule increases. In the limiting case the molecular bond breaks due to the electronic depletion induced in the covalent bond. An application of this study is the modelisation of the forces exerted on the passivated cantilever of an atomic force microscope probing the electron cloud in the contact mode in the Pauli exclusion regime.

Noble gas dimers confined inside C70

The potential energy surfaces for the interior rotation of a series of pairs of noble gas atoms encapsulated in the C₇₀ cavity have been explored.