1
|
Cosinschi M, Preda AT, Pantis-Simut CA, Filipoiu N, Ghitiu I, Dulea MA, Ion L, Manolescu A, Nemnes GA. Collective dynamics of Ca atoms encapsulated in C 60 endohedral fullerenes. Phys Chem Chem Phys 2024; 26:22090-22098. [PMID: 39118483 DOI: 10.1039/d4cp01048e] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/10/2024]
Abstract
Endohedral C60 fullerenes with up to four encapsulated Ca atoms were investigated by ab initio molecular dynamics simulations (AIMD). The relatively long runs allow us to describe the correlated movement of the Ca atoms inside the fullerene cage. For the systems with one or two Ca atoms a relatively unimpeded rotation was conjectured by earlier nuclear magnetic resonance experiments and supported by previous ab initio calculations used to sample the potential energy landscape. Here, by AIMD calculations, we confirm not only the circular motion, but also the correlated movement of the two Ca atoms, which is due to electric dipole interactions on the inner surface of the C60 molecule. Furthermore, systems with three and four Ca atoms present highly symmetric configurations of the embedded atoms, which are shown to rotate consistently within the fullerene cage, while more complex charge density patterns emerge. Employing artificial neural network models we perform a force-field mapping, which enables us to reproduce the main characteristics of the actual dynamics, such as the circular motion and the correlated movement of the Ca atoms.
Collapse
Affiliation(s)
- Mihaela Cosinschi
- University of Bucharest, Faculty of Physics, 077125 Magurele-Ilfov, Romania.
- Horia Hulubei National Institute for Physics and Nuclear Engineering, 077126 Magurele-Ilfov, Romania
| | - Amanda T Preda
- University of Bucharest, Faculty of Physics, 077125 Magurele-Ilfov, Romania.
- Horia Hulubei National Institute for Physics and Nuclear Engineering, 077126 Magurele-Ilfov, Romania
| | - C-A Pantis-Simut
- University of Bucharest, Faculty of Physics, 077125 Magurele-Ilfov, Romania.
- Horia Hulubei National Institute for Physics and Nuclear Engineering, 077126 Magurele-Ilfov, Romania
| | - N Filipoiu
- University of Bucharest, Faculty of Physics, 077125 Magurele-Ilfov, Romania.
- Horia Hulubei National Institute for Physics and Nuclear Engineering, 077126 Magurele-Ilfov, Romania
| | - I Ghitiu
- University of Bucharest, Faculty of Physics, 077125 Magurele-Ilfov, Romania.
- National Institute for Laser, Plasma and Radiation Physics, 077125 Magurele-Ilfov, Romania
| | - M A Dulea
- Horia Hulubei National Institute for Physics and Nuclear Engineering, 077126 Magurele-Ilfov, Romania
| | - L Ion
- University of Bucharest, Faculty of Physics, 077125 Magurele-Ilfov, Romania.
| | - A Manolescu
- Department of Engineering, School of Technology, Reykjavik University, Menntavegur 1, IS-102 Reykjavik, Iceland
| | - G A Nemnes
- University of Bucharest, Faculty of Physics, 077125 Magurele-Ilfov, Romania.
- Research Institute of the University of Bucharest (ICUB), 90 Panduri Street, 050663 Bucharest, Romania
- Horia Hulubei National Institute for Physics and Nuclear Engineering, 077126 Magurele-Ilfov, Romania
| |
Collapse
|
2
|
Panchagnula K, Graf D, Johnson ER, Thom AJW. Targeting spectroscopic accuracy for dispersion bound systems from ab initio techniques: Translational eigenstates of Ne@C70. J Chem Phys 2024; 161:054308. [PMID: 39092939 DOI: 10.1063/5.0223298] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/13/2024] [Accepted: 07/17/2024] [Indexed: 08/04/2024] Open
Abstract
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.
Collapse
Affiliation(s)
- K Panchagnula
- Yusuf Hamied Department of Chemistry, University of Cambridge, Cambridge, United Kingdom
| | - D Graf
- Yusuf Hamied Department of Chemistry, University of Cambridge, Cambridge, United Kingdom
- Department of Chemistry, University of Munich (LMU), Munich, Germany
| | - E R Johnson
- Yusuf Hamied Department of Chemistry, University of Cambridge, Cambridge, United Kingdom
- Department of Chemistry, Dalhousie University, 6243 Alumni Crescent, Halifax, Nova Scotia B3H 4R2, Canada
| | - A J W Thom
- Yusuf Hamied Department of Chemistry, University of Cambridge, Cambridge, United Kingdom
| |
Collapse
|
3
|
Panchagnula K, Graf D, Albertani FEA, Thom AJW. Translational eigenstates of He@C60 from four-dimensional ab initio potential energy surfaces interpolated using Gaussian process regression. J Chem Phys 2024; 160:104303. [PMID: 38465682 DOI: 10.1063/5.0197903] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/15/2024] [Accepted: 02/22/2024] [Indexed: 03/12/2024] Open
Abstract
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.
Collapse
Affiliation(s)
- K Panchagnula
- Yusuf Hamied Department of Chemistry, University of Cambridge, Cambridge, United Kingdom
| | - D Graf
- Yusuf Hamied Department of Chemistry, University of Cambridge, Cambridge, United Kingdom
| | - F E A Albertani
- Yusuf Hamied Department of Chemistry, University of Cambridge, Cambridge, United Kingdom
| | - A J W Thom
- Yusuf Hamied Department of Chemistry, University of Cambridge, Cambridge, United Kingdom
| |
Collapse
|
4
|
Bloodworth S, Whitby RJ. Synthesis of endohedral fullerenes by molecular surgery. Commun Chem 2022; 5:121. [PMID: 36697689 PMCID: PMC9814919 DOI: 10.1038/s42004-022-00738-9] [Citation(s) in RCA: 16] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/30/2022] [Accepted: 09/21/2022] [Indexed: 01/28/2023] Open
Abstract
Encapsulation of atoms or small molecules inside fullerenes provides a unique opportunity for study of the confined species in the isolated cavity, and the synthesis of closed C60 or C70 fullerenes with enclosed atoms or molecules has recently developed using the method of 'molecular surgery'; in which an open-cage intermediate fullerene is the host for encapsulation of a guest species, before repair of the cage opening. In this work we review the main methods for cage-opening and closure, and the achievements of molecular surgery to date.
Collapse
Affiliation(s)
- Sally Bloodworth
- grid.5491.90000 0004 1936 9297Chemistry, Faculty of Engineering and Physical Sciences, University of Southampton, Southampton, SO17 1BJ UK
| | - Richard J. Whitby
- grid.5491.90000 0004 1936 9297Chemistry, Faculty of Engineering and Physical Sciences, University of Southampton, Southampton, SO17 1BJ UK
| |
Collapse
|
5
|
Morcillo-Arencibia MF, Alcaraz-Pelegrina JM, Sarsa AJ, Randazzo JM. An off-center endohedrally confined hydrogen molecule. Phys Chem Chem Phys 2022; 24:22971-22977. [PMID: 36125249 DOI: 10.1039/d2cp03456e] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
In this study, we address the problem of a C60 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 C60 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.
Collapse
Affiliation(s)
- Milagros F Morcillo-Arencibia
- Departamento de Física, Campus de Rabanales, Edif. C2. Universidad de Córdoba, E-14071 Córdoba, Spain. .,Centro Atómico Bariloche, CNEA and CONICET, S. C. de Bariloche, Río Negro, Argentina
| | | | - Antonio J Sarsa
- Departamento de Física, Campus de Rabanales, Edif. C2. Universidad de Córdoba, E-14071 Córdoba, Spain.
| | - Juan M Randazzo
- Centro Atómico Bariloche, CNEA and CONICET, S. C. de Bariloche, Río Negro, Argentina
| |
Collapse
|