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

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 .

About Perfluoropolyhedranes, Their Electron-Accepting Ability and Questionable Supramolecular Hosting Capacity

Polyhedral molecules are appealing for their eye-catching architecture and distinctive chemistry. Perfluorination of such, often greatly strained, compounds is a momentous challenge. It drastically changes the electron distribution, structure and properties. Notably, small high-symmetry perfluoropolyhedranes feature a centrally located, star-shaped low-energy unoccupied molecular orbital that can host an extra electron within the polyhedral frame, thus producing a radical anion, without loss of symmetry. This predicted electron-hosting capacity was definitively established for perfluorocubane, the first perfluorinated Platonic polyhedrane to be isolated pure. Hosting atoms, molecules, or ions in such "cage" structures is, however, all but forthright, if not illusionary, offering no easy access to supramolecular constructs. While adamantane and cubane have fostered numerous applications in materials science, medicine, and biology, specific uses for their perfluorinated counterparts remain to be established. Some aspects of highly fluorinated carbon allotropes, such as fullerenes and graphite, are briefly mentioned for context.

On the nature of the Schottky anomaly in endohedral water

In this work, we study the heat capacity contribution of a rigid water molecule encapsulated in C₆₀ by performing six-dimensional eigenstate calculations with the inclusion of its quantized rotational and translational degrees of freedom. Two confinement model potentials are considered: in the first, confinement is described using distributed pairwise Lennard-Jones interactions, while in the second, the water molecule is trapped within an eccentric but isotropic 3D harmonic effective confinement potential [Wespiser et al., J. Chem. Phys. 156, 074304 (2022)]. Contributions to the heat capacity from both the ortho and para nuclear spin isomers of water are considered to enable the effects of their interconversion to be assessed. By including a symmetry-breaking quadrupolar potential energy term in the Hamiltonian, we can reproduce the experimentally observed Schottky anomaly at ∼2 K [Suzuki et al., J. Phys. Chem. Lett. 10, 1306 (2019)]. Furthermore, our calculations predict a second Schottky anomaly at ∼0.1 K resulting from the H configuration, a different orientational arrangement of the fullerene cages in crystalline solid C₆₀. Contributions from the H configuration to C_V also explain the second peak observed at ∼7 K in the experimentally measured heat capacity.

H2 + H2O → H4O: Synthesizing Hyperhydrogenated Water in Small-Sized Fullerenes?

Nanoscale confinement provides an ideal platform to rouse some exceptional reactions which cannot happen in the open space. Intuitively, H₂ and H₂O cannot react. Herein, through utilizing small-sized fullerenes (C₂₄, C₂₆, C₂₈, and C₃₀) as nanoreactors, we demonstrate that a hyperhydrogenated water species, H₄O, can be easily formed using H₂ and H₂O under ambient conditions by ab initio molecular dynamics simulations. The H₄O molecule rotates freely in the cavity of the cages and maintains its structure during the simulations. Further theoretical analysis indicates that H₄O in the fullerene possesses high stability thermodynamically and chemically, which can be rationalized by the electron transfer between H₄O and the fullerene. This work highlights the possibility of utilizing fullerene as a nanoreactor to provide confinement constraints for unexpected chemistry.

Smolyak Algorithm Adapted to a System-Bath Separation: Application to an Encapsulated Molecule with Large-Amplitude Motions

A Smolyak algorithm adapted to system-bath separation is proposed for rigorous quantum simulations. This technique combines a sparse grid method with the system-bath concept in a specific configuration without limitations on the form of the Hamiltonian, thus achieving a highly efficient convergence of the excitation transitions for the "system" part. Our approach provides a general way to overcome the perennial convergence problem for the standard Smolyak algorithm and enables the simulation of floppy molecules with more than a hundred degrees of freedom. The efficiency of the present method is illustrated on the simulation of H₂ caged in an sII clathrate hydrate including two kinds of cage modes. The transition energies are converged by increasing the number of normal modes of water molecules. Our results confirm the triplet splittings of both translational and rotational (j = 1) transitions of the H₂ molecule. Furthermore, they show a slight increase of the translational transitions with respect to the ones in a rigid cage.

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

We explore the origin of the anomalous splitting of the 101 levels reported experimentally for the H2 O@C60 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 C60 -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 C60 -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.

First principles molecular dynamics calculations of the mechanical properties of endofullerenes containing noble gas atoms or small molecules

The mechanical properties of endofullerenes have been investigated by performing compression tests using finite temperature first principles molecular dynamics calculations. We considered various X@C₆₀ systems, with X a single noble gas atom (He, Ne, Ar, Kr, or Xe), small molecules (H₂O, CH₄), or small helium clusters. In the absence of compression, it is observed that there is no or at best a negligible effect of X on the properties of C₆₀. The compression simulations revealed several original findings. First, the influence of X on the stiffness of X@C₆₀ can be quantified, although it is at most 12% for the studied cases. Next, both energy and contact force variations as a function of strain are demonstrated to depend on X. However, this is not the case for the yield strain and for the failure mechanism of the C₆₀ shell. Finally, it is shown that the X@C₆₀ compression could bring X to be in a high stress state. In the specific cases of H₂O and CH₄ molecules, a mechanism of stress assisted dissociation is observed.

Vibronic Coupling in Spherically Encapsulated, Diatomic Molecules: Prediction of a Renner-Teller-like Effect for Endofullerenes

In the year 1933, Herzberg and Teller realized that the potential energy surface of a triatomic, linear molecule splits into two as soon as the molecule is bent. The phenomenon, later dubbed the Renner–Teller effect due to the detailed follow-up work of Renner on the subject, describes the coupling of a symmetry-reducing molecular vibration with degenerate electronic states. In this article, we show that a very similar type of nonadiabatic coupling can occur for certain translational degrees of freedom of diatomic, electronically degenerate molecules when trapped in a nearly spherical or cylindrical quantum confinement, e.g., realized through electromagnetic fields or molecular encapsulation. We illustrate this on the example of fullerene-encapsulated nitric oxide, and provide a prediction of its interesting, perturbed vibronic spectrum.