1
|
Prosmiti R, González-Lezana T. Computational Modeling: Up-to-Date Approaches and Cutting-Edge Applications from Clusters, Nanostructures to Bulk Systems. Chemphyschem 2024:e202400207. [PMID: 38837591 DOI: 10.1002/cphc.202400207] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/26/2024] [Indexed: 06/07/2024]
Abstract
The contributions in this special theme collection, in honor to Prof. P. Villarreal, cover a broad variety of computational methodologies and experimental techniques, containing studies on gas phase, clusters and condensed phase systems.
Collapse
Affiliation(s)
- Rita Prosmiti
- Instituto de Fìsica Fundamental, IFF-CSIC, Serrano 123, 28006, Madrid, Spain
| | | |
Collapse
|
2
|
Montes de Oca-Estévez MJ, Valdés Á, Prosmiti R. A kernel-based machine learning potential and quantum vibrational state analysis of the cationic Ar hydride (Ar 2H +). Phys Chem Chem Phys 2024; 26:7060-7071. [PMID: 38345626 DOI: 10.1039/d3cp05865d] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/23/2024]
Abstract
One of the most fascinating discoveries in recent years, in the cold and low pressure regions of the universe, was the detection of ArH+ and HeH+ species. The identification of such noble gas-containing molecules in space is the key to understanding noble gas chemistry. In the present work, we discuss the possibility of [Ar2H]+ existence as a potentially detectable molecule in the interstellar medium, providing new data on possible astronomical pathways and energetics of this compound. As a first step, a data-driven approach is proposed to construct a full 3D machine-learning potential energy surface (ML-PES) via the reproducing kernel Hilbert space (RKHS) method. The training and testing data sets are generated from CCSD(T)/CBS[56] computations, while a validation protocol is introduced to ensure the quality of the potential. In turn, the resulting ML-PES is employed to compute vibrational levels and molecular spectroscopic constants for the cation. In this way, the most common isotopologue in ISM, [36Ar2H]+, was characterized for the first time, while simultaneously, comparisons with previously reported values available for [40Ar2H]+ are discussed. Our present data could serve as a benchmark for future studies on this system, as well as on higher-order cationic Ar-hydrides of astrophysical interest.
Collapse
Affiliation(s)
- María Judit Montes de Oca-Estévez
- Institute of Fundamental Physics (IFF-CSIC), CSIC, Serrano 123, 28006 Madrid, Spain.
- Atelgraphics S.L., Mota de Cuervo 42, 28043, Madrid, Spain
| | - Álvaro Valdés
- Escuela de Física, Universidad Nacional de Colombia, Sede Medellín, A. A., 3840, Medellín, Colombia
| | - Rita Prosmiti
- Institute of Fundamental Physics (IFF-CSIC), CSIC, Serrano 123, 28006 Madrid, Spain.
| |
Collapse
|
3
|
Mabrouk N, Dhiflaoui J, Saidi S, Bejaoui M, Alharzali N, Berriche H. Potential Energy Surface and Bound States of Ne-Li 2+( X2Σ g+) van der Waals Complex Based on Ab Initio Calculations. J Phys Chem A 2023; 127:9167-9177. [PMID: 37890154 PMCID: PMC10641847 DOI: 10.1021/acs.jpca.3c03811] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/06/2023] [Revised: 10/07/2023] [Accepted: 10/10/2023] [Indexed: 10/29/2023]
Abstract
Theoretical studies of the potential energy surface and vibrational bound states calculations were performed for the ground state of the Ne-Li2+(X2Σg+) van der Waals (vdW) complex. The intermolecular interactions were investigated by using an accurate monoconfigurational RCCSD(T) method and large basis sets (aug-cc-pVnZ, n = T, Q, 5), extrapolated to the complete basis set (CBS) limit. In turn, the obtained raw data from RCCSD(T)/CBS(Q5) calculations were numerically interpolated using the Morse + vdW model and the Reproducing Kernel Hilbert Space (RKHS) polynomial method to generate analytic expressions for the 2D-PES. The RKHS interpolated PES was then used to assess the bound states of the Ne-Li2+(X2Σg+) system through nuclear quantum calculations. By studying the aspect of the potential energy surface, the analysis sheds light on the behavior of the Ne-Li2+(X2Σg+) complex and its interactions between repulsive and attractive forces with other particles. By examining the vibrational states and wave functions of the system, the researchers were able to gain a better understanding of the behavior of the Ne-Li2+(X2Σg+) complex. The calculated radial and angular distributions for all even and odd symmetries are discussed in detail. We observe that the radial distributions exhibit a more complicated nodal structure, representing stretching vibrational behavior in the neon atom along its radial coordinate. For the highest bound states, the situation is very different, and the energies surpass the angular barrier.
Collapse
Affiliation(s)
- Nesrine Mabrouk
- Laboratory
of Interfaces and Advanced Materials, Physics Department, Faculty
of Sciences of Monastir, University of Monastir, Monastir 5019, Tunisia
| | - Jamila Dhiflaoui
- Laboratory
of Interfaces and Advanced Materials, Physics Department, Faculty
of Sciences of Monastir, University of Monastir, Monastir 5019, Tunisia
| | - Samah Saidi
- Laboratory
of Interfaces and Advanced Materials, Physics Department, Faculty
of Sciences of Monastir, University of Monastir, Monastir 5019, Tunisia
- Department
of Physics, College of Science and Humanities in Al-Kharj, Prince Sattam bin Abdulaziz University, Al-Kharj 16273, Saudi Arabia
| | - Mohamed Bejaoui
- Laboratory
of Interfaces and Advanced Materials, Physics Department, Faculty
of Sciences of Monastir, University of Monastir, Monastir 5019, Tunisia
| | - Nissrin Alharzali
- Laboratory
of Interfaces and Advanced Materials, Physics Department, Faculty
of Sciences of Monastir, University of Monastir, Monastir 5019, Tunisia
- Department
of Physical and Theoretical Chemistry, Faculty of Natural Sciences, Comenius University in Bratislava, Bratislava 814 99, Slovakia
| | - Hamid Berriche
- Laboratory
of Interfaces and Advanced Materials, Physics Department, Faculty
of Sciences of Monastir, University of Monastir, Monastir 5019, Tunisia
- Mathematics
and Physics Department, School of Arts and Sciences, American University of Ras Al Khaimah, Ras Al-Khaimah 10021, United Arab Emirates
| |
Collapse
|
4
|
Saidi S, Mabrouk N, Dhiflaoui J, Berriche H. Structural, Spectroscopic, and Dynamic Properties of Li2+(X2∑g+) in Interaction with Krypton Atom. Molecules 2023; 28:5512. [PMID: 37513385 PMCID: PMC10385072 DOI: 10.3390/molecules28145512] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/23/2023] [Revised: 07/10/2023] [Accepted: 07/17/2023] [Indexed: 07/30/2023] Open
Abstract
We report a computational study of the potential energy surface (PES) and vibrational bound states for the ground electronic state of Li2+Kr. The PES was calculated in Jacobi coordinates at the Restricted Coupled Cluster method RCCSD(T) level of calculation and using aug-cc-pVnZ (n = 4 and 5) basis sets. Afterward, this PES is extrapolated to the complete basis set (CBS) limit for correction. The obtained interaction energies were, then, interpolated numerically using the reproducing kernel Hilbert space polynomial (RKHS) approach to produce analytic expressions for the 2D-PES. The analytical PES is used to solve the nuclear Schrodinger equation to determine the bound states' eigenvalues of Li2+Kr for a J = 0 total angular momentum configuration and to understand the effects of orientational anisotropy of the forces and the interplay between the repulsive and attractive interaction within the potential surface. In addition, the radial and angular distributions of some selected bound state levels, which lie below, around, and above the T-shaped 90° barrier well, are calculated and discussed. We note that the radial distributions clearly acquire a more complicated nodal structure and correspond to bending and stretching vibrational motions "mode" of the Kr atom along the radial coordinate, and the situation becomes very different at the highest bound states levels with energies higher than the T-shaped 90° barrier well. The shape of the distributions becomes even more complicated, with extended angular distributions and prominent differences between even and odd states.
Collapse
Affiliation(s)
- Samah Saidi
- Laboratory of Interfaces and Advanced Materials, Physics Department, Faculty of Sciences of Monastir, Avenue de l'Environnement, Monastir 5019, Tunisia
- Department of Physics, College of Science and Humanities in Al-Kharj, Prince Sattam Bin Abdulaziz University, Al-Kharj 11942, Saudi Arabia
| | - Nesrine Mabrouk
- Laboratory of Interfaces and Advanced Materials, Physics Department, Faculty of Sciences of Monastir, Avenue de l'Environnement, Monastir 5019, Tunisia
| | - Jamila Dhiflaoui
- Laboratory of Interfaces and Advanced Materials, Physics Department, Faculty of Sciences of Monastir, Avenue de l'Environnement, Monastir 5019, Tunisia
| | - Hamid Berriche
- Laboratory of Interfaces and Advanced Materials, Physics Department, Faculty of Sciences of Monastir, Avenue de l'Environnement, Monastir 5019, Tunisia
- Mathematics and Natural Sciences Department, School of Arts and Sciences, American University of Ras Al Khaimah, Ras Al-Khaimah 10021, United Arab Emirates
| |
Collapse
|
5
|
Ghanmi C, Nakbi H, Al-Qarni HJ, Alharzali N, Berriche H. Structure, energetics, and spectroscopy of the K 2+(X 2Σ +g) interacting with the noble gas atoms Ar, Kr and Xe. J Mol Graph Model 2023; 120:108413. [PMID: 36758327 DOI: 10.1016/j.jmgm.2023.108413] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/07/2022] [Revised: 12/31/2022] [Accepted: 01/12/2023] [Indexed: 02/05/2023]
Abstract
The structure, energetic, and spectroscopy properties of the ionic system K2+(X2Σ+g) interacting with the noble gas atoms Argon, Krypton and Xenon are studied. The computations are done by an accurate ab initio approach based on the pseudo-potential technique, Gaussian basis sets, parameterized l-dependent polarization potentials and an analytic potential form for the K+Ar, K+Kr and K+Xe interactions. These interactions are added via the CCSD(T) potential taken from literature and fitted applying the analytical expression of Tang and Toennies. The application of the pseudo-potential approach reduces the number of active electrons of each complex to only one electron. The potential energy surfaces are analyzed on a large range of the Jacobi coordinates, R and θ. By the general interpolation outline based on the RKHS (Reproducing Kernel Hilbert Space) procedure, we have reproduced for each complex from our ab initio results the two-dimensional contour plots of an analytical potential. To evaluate the stability of each complex, we have determined from the potential energy surfaces the equilibrium distance (Re), the well depth (De), the quantum energy (D0), the zero-point-energy (ZPE) and the ZPE%. The results showed that the linear configurations, where the noble gas atom connected to the K2+(X2Σ+g) system, are the more stable.
Collapse
Affiliation(s)
- Chedli Ghanmi
- Laboratory of Interfaces and Advanced Materials, Faculty of Science, University of Monastir, 5019, Monastir, Tunisia.
| | - Haifa Nakbi
- Laboratory of Interfaces and Advanced Materials, Faculty of Science, University of Monastir, 5019, Monastir, Tunisia
| | - Hind Jahman Al-Qarni
- Physics Department, College of Science and Arts, Balqarn, Bisha University, Saudi Arabia
| | - Nisrin Alharzali
- Laboratory of Interfaces and Advanced Materials, Faculty of Science, University of Monastir, 5019, Monastir, Tunisia
| | - Hamid Berriche
- Laboratory of Interfaces and Advanced Materials, Faculty of Science, University of Monastir, 5019, Monastir, Tunisia; Mathematics and Natural Sciences Department, School of Arts and Sciences, American University of Ras Al Khaimah, Ras Al Khaimah, United Arab Emirates.
| |
Collapse
|
6
|
Kollotzek S, Campos-Martínez J, Bartolomei M, Pirani F, Tiefenthaler L, Hernández MI, Lázaro T, Zunzunegui-Bru E, González-Lezana T, Bretón J, Hernández-Rojas J, Echt O, Scheier P. Helium nanodroplets as an efficient tool to investigate hydrogen attachment to alkali cations. Phys Chem Chem Phys 2022; 25:462-470. [PMID: 36477158 PMCID: PMC9768848 DOI: 10.1039/d2cp03841b] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/19/2022] [Accepted: 11/16/2022] [Indexed: 11/18/2022]
Abstract
We report a novel method to reversibly attach and detach hydrogen molecules to positively charged sodium clusters formed inside a helium nanodroplet host matrix. It is based on the controlled production of multiply charged helium droplets which, after picking up sodium atoms and exposure to H2 vapor, lead to the formation of Nam+(H2)n clusters, whose population was accurately measured using a time-of-flight mass spectrometer. The mass spectra reveal particularly favorable Na+(H2)n and Na2+(H2)n clusters for specific "magic" numbers of attached hydrogen molecules. The energies and structures of these clusters have been investigated by means of quantum-mechanical calculations employing analytical interaction potentials based on ab initio electronic structure calculations. A good agreement is found between the experimental and the theoretical magic numbers.
Collapse
Affiliation(s)
- Siegfried Kollotzek
- University of Innsbruck, Institute for Ion Physics and Applied Physics, Innsbruck, Austria.
| | | | | | - Fernando Pirani
- Dipartimento di Chimica, Biologia e Biotecnologie, Università di Perugia, Perugia, Italy
| | - Lukas Tiefenthaler
- University of Innsbruck, Institute for Ion Physics and Applied Physics, Innsbruck, Austria.
| | | | - Teresa Lázaro
- Instituto de Física Fundamental, C.S.I.C., Madrid, Spain.
| | | | | | - José Bretón
- Departamento de Física and IUdEA, Universidad de La Laguna, La Laguna, Tenerife, Spain
| | | | - Olof Echt
- University of Innsbruck, Institute for Ion Physics and Applied Physics, Innsbruck, Austria.
- Department of Physics, University of New Hampshire, Durham, NH 03824, USA
| | - Paul Scheier
- University of Innsbruck, Institute for Ion Physics and Applied Physics, Innsbruck, Austria.
| |
Collapse
|
7
|
Rock CA, Arradondo SN, Tschumper GS. Solvation of Isoelectronic Halide and Alkali Metal Ions by Argon Atoms. J Phys Chem A 2021; 125:10524-10531. [PMID: 34851634 DOI: 10.1021/acs.jpca.1c08069] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/18/2023]
Abstract
This work systematically examines the interactions of alkali metal cations and their isoelectronic halide counterparts with up to six solvating Ar atoms (M+Arn and X-Arn, where M = Li, Na, K, and Rb; X = H, F, Cl, and Br; and n = 1-6) via full geometry optimizations with the MP2 method and robust, correlation-consistent quadruple-ζ (QZ) basis sets. 116 unique M+Arn and X-Arn stationary points have been characterized on the MP2/QZ potential energy surface. To the best of our knowledge, approximately two dozen of these stationary points have been reported here for the first time. Some of these new structures are either the lowest-energy stationary point for a particular cluster or energetically competitive with it. The CCSD(T) method was employed to perform additional single-point energy computations upon all MP2/QZ-optimized structures using the same basis set. CCSD(T)/QZ results indicate that internally solvated structures with the ion at/near the geometric center of the cluster have appreciably higher energies than those placing the ion on the periphery. While this study extends the prior investigations of M+Arn clusters found within the literature, it notably provides one of the first thorough characterizations of and comparisons to the corresponding negatively charged X-Arn clusters.
Collapse
Affiliation(s)
- Carly A Rock
- Department of Chemistry and Biochemistry, University of Mississippi, University, Mississippi 38677-1848, United States
| | - Sarah N Arradondo
- Department of Chemistry, Washington College, Chestertown, Maryland 21620-1438, United States
| | - Gregory S Tschumper
- Department of Chemistry and Biochemistry, University of Mississippi, University, Mississippi 38677-1848, United States
| |
Collapse
|
8
|
Alharzali N, Rodríguez-Segundo R, Prosmiti R. Modelling interactions of cationic dimers in He droplets: microsolvation trends in He nK 2+ clusters. Phys Chem Chem Phys 2021; 23:7849-7859. [PMID: 33220666 DOI: 10.1039/d0cp05406b] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
We report the results of a detailed theoretical investigation of small K2+-doped He clusters. The structural characteristics and stabilities of such cations are determined from ab initio electronic structure calculations at the MRCI+Q level of theory. The underlying interactions show a multireference character and such effects are analyzed. The interaction potentials are constructed employing an interpolation technique within the inverse problem theory method, while the nuclear quantum effects are computed for the trimers, their spatial arrangements are discussed, and information was extracted on the orientational anisotropy of the forces. We found that energetically the most stable conformer corresponds to linear arrangements that are taking place under large amplitude vibrations, with high zero-point energy. We have further looked into the behavior of higher-order species with various He atoms surrounding the cationic dopant. By using a sum of potentials approach and an evolutionary programming method, we analyzed the structural stability of clusters with up to six He atoms in comparison with interactions energies obtained from MRCI+Q quantum chemistry computations. Structures containing Hen motifs that characterize pure rare gas clusters, appear for the larger K2+-doped He clusters, showing selective growth during the microsolvation process of the alkali-dimer cation surrounded by He atoms. Such results indicate the existence of local solvation microstructures in these aggregates, where the cationic impurity could get trapped for a short time, contributing to the slow ionic mobility observed experimentally in ultra-cold He-droplets.
Collapse
Affiliation(s)
- Nissrin Alharzali
- Institute of Fundamental Physics (IFF-CSIC), CSIC, Serrano 123, 28006 Madrid, Spain.
| | | | | |
Collapse
|