-
Targeting spectroscopic accuracy for dispersion bound systems from ab initio techniques: translational eigenstates of Ne@C$_{70}$
Authors:
K. Panchagnula,
D. Graf,
E. R. Johnson,
A. J. W. Thom
Abstract:
We investigate the endofullerene system Ne@C$_{70}$, by constructing a three-dimensional Potential Energy Surface (PES) describing the translational motion of the Ne atom. We compare a plethora of electronic structure methods including: MP2, SCS-MP2, SOS-MP2, RPA@PBE, C(HF)-RPA, which were previously used for He@C$_{60}$ in J. Chem. Phys. 160, 104303 (2024), alongside B86bPBE-25X-XDM and B86bPBE-5…
▽ More
We investigate the endofullerene system Ne@C$_{70}$, by constructing a three-dimensional Potential Energy Surface (PES) describing the translational motion of the Ne atom. We compare a plethora of electronic structure methods including: MP2, SCS-MP2, SOS-MP2, RPA@PBE, C(HF)-RPA, which were previously used for He@C$_{60}$ in J. Chem. Phys. 160, 104303 (2024), alongside B86bPBE-25X-XDM and B86bPBE-50X-XDM. The reduction in symmetry moving from C$_{60}$ to C$_{70}$ introduces a double well potential along the anisotropic direction, which forms a test of the sensitivity and effectiveness of the methods. Due to the large cost of these calculations, the PES is interpolated using Gaussian Process Regression due to its effectiveness with sparse training data. The nuclear Hamiltonian is diagonalised using a symmetrised double minimum basis set outlined in J. Chem. Phys. 159, 164308 (2023), with translational energies having error bars $\pm 1$ and $\pm 2$ cm$^{-1}$. We quantify the shape of the ground state wavefunction by considering its prolateness and kurtosis, and compare the eigenfunctions between electronic structure methods from their Hellinger distance. 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 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@C$_{60}$, due to the lack of concordance between all electronic structure methods we require more experimental data in order to properly validate the choice.
△ Less
Submitted 11 June, 2024;
originally announced June 2024.
-
Translational eigenstates of He@C$_{60}$ from four-dimensional \textit{ab initio} Potential Energy Surfaces interpolated using Gaussian Process Regression
Authors:
K. Panchagnula,
D. Graf,
F. E. A. Albertani,
A. J. W. Thom
Abstract:
We investigate the endofullerene system $^3$He@C$_{60}$ with a four-dimensional Potential Energy Surface (PES) to include the three He translational degrees of freedom and C$_{60}$ cage radius. We compare MP2, SCS-MP2, SOS-MP2, RPA@PBE and C(HF)-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 usi…
▽ More
We investigate the endofullerene system $^3$He@C$_{60}$ with a four-dimensional Potential Energy Surface (PES) to include the three He translational degrees of freedom and C$_{60}$ cage radius. We compare MP2, SCS-MP2, SOS-MP2, RPA@PBE and C(HF)-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 diagonalised 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 96cm$^{-1}$ and 110cm$^{-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 approximately $\pm$ 1.5cm$^{-1}$. Wavefunctions are also compared by considering their overlap and Hellinger distance to the one-dimensional empirical potential. As with the energies, the two \textit{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.
△ Less
Submitted 18 March, 2024; v1 submitted 9 January, 2024;
originally announced January 2024.
-
Exploring the parameter space of an endohedral atom in a cylindrical cavity
Authors:
K. Panchagnula,
A. J. W. Thom
Abstract:
Endohedral fullerenes, or endofullerenes, are chemical systems of fullerene cages encapsulating single atoms or small molecules. These species provide an interesting challenge of Potential Energy Surface (PES) determination as examples of non-covalently bonded, bound systems. While the majority of studies focus on C$_{60}$ as the encapsulating cage, introducing some anisotropy by using a different…
▽ More
Endohedral fullerenes, or endofullerenes, are chemical systems of fullerene cages encapsulating single atoms or small molecules. These species provide an interesting challenge of Potential Energy Surface (PES) determination as examples of non-covalently bonded, bound systems. While the majority of studies focus on C$_{60}$ as the encapsulating cage, introducing some anisotropy by using a different fullerene, e.g., C$_{70}$ can unveil a double well potential along the unique axis. By approximating the potential as a pairwise Lennard-Jones (LJ) summation over the fixed C cage atoms, the parameter space of the Hamiltonian includes three tunable variables: $(M,\varepsilon,σ)$ representing the mass of the trapped species, the LJ energy, and length scales respectively. Fixing the mass and allowing the others to vary can imitate the potentials of endohedral species trapped in more elongated fullerenes. We choose to explore the LJ parameter space of an endohedral atom in C$_{70}$ with $\varepsilon\in$ [20cm$^{-1}$, 150cm$^{-1}$], and $σ\in$ [2.85Å , 3.05Å].
As the barrier height and positions of these wells vary between [1cm$^{-1}$, 264cm$^{-1}$] and [0.35Å, 0.85Å] respectively, using a 3D direct product basis of 1D harmonic oscillator (HO) wavefunctions centred at the origin where there is a local maximum is unphysical. Instead we propose the use of a non-orthogonal basis set, using 1D HO wavefunctions centred in each minimum and compare this to other choices. The ground state energy of the X@C$_{70}$ is tracked across the LJ parameter space, along with its corresponding nuclear translational wavefunctions. A classification of the wavefunction characteristics, namely the prolateness and ``peanut-likeness'' based on its statistical moments is also proposed.
△ Less
Submitted 20 November, 2023; v1 submitted 14 June, 2023;
originally announced June 2023.