1
|
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
|
2
|
Rodríguez-Segundo R, Gijón A, Prosmiti R. Quantum molecular simulations of micro-hydrated halogen anions. Phys Chem Chem Phys 2022; 24:14964-14974. [PMID: 35686995 DOI: 10.1039/d2cp01396g] [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
We report the results of a detailed and accurate investigation focused on structures and energetics of poly-hydrated halides employing first-principles polarizable halide-water potentials to describe the underlying forces. Following a bottom-up data-driven potential approach, we initially looked into the classical behavior of higher-order X-(H2O)N clusters. We have located several low-lying energies, such as global and local minima, structures for each cluster, with various water molecules (up to N = 8) surrounding the halide anion (X- = F-, Cl-, Br-, I-), employing an evolutionary programming method. It is found that the F--water clusters exhibit different structural configurations than the heavier halides, however independently of the halide anion, all clusters show in general a selective growth with the anion preferring to be connected to the outer shell of the water molecule arrangements. In turn, path-integral molecular dynamics simulations are performed to incorporate explicitly nuclear quantum and thermal effects in describing the nature of halide ion microsolvation in such prototypical model systems. Our data reveal that at low finite temperatures, nuclear quantum effects affect certain structural properties, such as weakening hydrogen bonding between the halide anion and water molecules, with minor distortions in the water network beyond the first hydration shell, indicating local structure rearrangements. Such structural characteristics and the promising cluster size trends observed in the single-ion solvation energies motivated us to draw connections of small size cluster data to the limits of continuum bulk values, toward the investigation of the challenging computational modeling of bulk single ion hydration.
Collapse
Affiliation(s)
- Raúl Rodríguez-Segundo
- Institute of Fundamental Physics (IFF-CSIC), CSIC, Serrano 123, 28006 Madrid, Spain. .,Atelgraphics S.L., Mota de Cuervo 42, 28043, Madrid, Spain
| | - Alfonso Gijón
- Materials Science Institute of Madrid (ICMM-CSIC), CSIC, Sor Juana Inés de la Cruz 3, 28049 Madrid, Spain
| | - Rita Prosmiti
- Institute of Fundamental Physics (IFF-CSIC), CSIC, Serrano 123, 28006 Madrid, Spain.
| |
Collapse
|
3
|
A Benchmark Protocol for DFT Approaches and Data-Driven Models for Halide-Water Clusters. MOLECULES (BASEL, SWITZERLAND) 2022; 27:molecules27051654. [PMID: 35268757 PMCID: PMC8924895 DOI: 10.3390/molecules27051654] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 02/04/2022] [Revised: 02/18/2022] [Accepted: 02/26/2022] [Indexed: 11/17/2022]
Abstract
Dissolved ions in aqueous media are ubiquitous in many physicochemical processes, with a direct impact on research fields, such as chemistry, climate, biology, and industry. Ions play a crucial role in the structure of the surrounding network of water molecules as they can either weaken or strengthen it. Gaining a thorough understanding of the underlying forces from small clusters to bulk solutions is still challenging, which motivates further investigations. Through a systematic analysis of the interaction energies obtained from high-level electronic structure methodologies, we assessed various dispersion-corrected density functional approaches, as well as ab initio-based data-driven potential models for halide ion-water clusters. We introduced an active learning scheme to automate the generation of optimally weighted datasets, required for the development of efficient bottom-up anion-water models. Using an evolutionary programming procedure, we determined optimized and reference configurations for such polarizable and first-principles-based representation of the potentials, and we analyzed their structural characteristics and energetics in comparison with estimates from DF-MP2 and DFT+D quantum chemistry computations. Moreover, we presented new benchmark datasets, considering both equilibrium and non-equilibrium configurations of higher-order species with an increasing number of water molecules up to 54 for each F, Cl, Br, and I anions, and we proposed a validation protocol to cross-check methods and approaches. In this way, we aim to improve the predictive ability of future molecular computer simulations for determining the ongoing conflicting distribution of different ions in aqueous environments, as well as the transition from nanoscale clusters to macroscopic condensed phases.
Collapse
|
4
|
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
|
5
|
Andreev AS, Baturo VV, Lukashov SS, Poretsky SA, Pravilov AM, Zhironkin AI. Experimental and theoretical investigations of HeNeI 2 trimer. J Chem Phys 2020; 152:234307. [DOI: 10.1063/5.0008760] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Aleksei S. Andreev
- Department of Physics, Saint Petersburg State University, SPbSU, 7/9 Universitetskaya nab., Saint Petersburg 199034, Russia
| | - Vera V. Baturo
- Department of Physics, Saint Petersburg State University, SPbSU, 7/9 Universitetskaya nab., Saint Petersburg 199034, Russia
| | - Sergey S. Lukashov
- Department of Physics, Saint Petersburg State University, SPbSU, 7/9 Universitetskaya nab., Saint Petersburg 199034, Russia
| | - Sergey A. Poretsky
- Department of Physics, Saint Petersburg State University, SPbSU, 7/9 Universitetskaya nab., Saint Petersburg 199034, Russia
| | - Anatoly M. Pravilov
- Department of Physics, Saint Petersburg State University, SPbSU, 7/9 Universitetskaya nab., Saint Petersburg 199034, Russia
| | - Anatoly I. Zhironkin
- Department of Physics, Saint Petersburg State University, SPbSU, 7/9 Universitetskaya nab., Saint Petersburg 199034, Russia
| |
Collapse
|
6
|
Valdés Á, Prosmiti R. Quantum effects on the stability of the He 5 I 2 van der Waals conformers. J Comput Chem 2019; 40:2200-2206. [PMID: 31148224 DOI: 10.1002/jcc.25870] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/11/2019] [Revised: 04/10/2019] [Accepted: 05/16/2019] [Indexed: 11/05/2022]
Abstract
We present 15-dimensional quantum multiconfiguration time-dependent Hartree calculations of the vibrational levels of the He5 I2 van der Waals (vdW) complex employing an ab initio-based potential energy surface (PES). The energies and spatial features of such bound structures are analyzed, providing predictions on the structures and relative stabilities of its three lowest isomers. We found that the most stable isomer corresponds to all five He atoms encircling the I2 molecule, indicating that in this case the anharmonic quantum effects do not stabilize the isomers involving a He atom in a linear configuration as reported previously for the smaller HeN I2 systems. Such finding provides information on the overall structuring of the finite-size-solvent systems, highlighting the intriguing interplay between weak intermolecular interactions and quantum effects. © 2019 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Álvaro Valdés
- Departamento de Física, Universidad Nacional de Colombia, Calle 26, Cra 39, Edificio 404, Bogotá, Colombia
| | - Rita Prosmiti
- Departamento PAMS, Instituto de Física Fundamental (IFF-CSIC), CSIC, Serrano 123, 28006, Madrid, Spain
| |
Collapse
|
7
|
Alharzali N, Berriche H, Villarreal P, Prosmiti R. Theoretical Study of Cationic Alkali Dimers Interacting with He: Li 2+-He and Na 2+-He van der Waals Complexes. J Phys Chem A 2019; 123:7814-7821. [PMID: 31442041 DOI: 10.1021/acs.jpca.9b05551] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
We present a theoretical study on the potential energy surface and bound states of He-A2+ complexes, where A is one of the alkali Li or Na atoms. The intermolecular interactions were systematically investigated by high-level ab initio electronic structure computations, and the corresponding raw data were then employed to reproduce accurate analytical expressions of the potential surfaces. In turn, we used these potentials to evaluate bound configurations of the trimers from nuclear quantum calculations and to extract information on the effect of orientational anisotropy of the forces and the interplay between repulsive and attractive interaction within the potential surfaces. The spatial features of the bound states are analyzed and discussed in detail. We found that both systems are going under large amplitude stretching and bending motions with high zero-point energies. Despite the large differences in the potential well-depths, the correct treatment of nuclear quantum effects provides insights on the effect of different strength of the ionic interaction on the spectral structure of such cationic alkali van der Waals complexes, related to the mobility of ions and the formation of cold-molecules in He-controlled environments.
Collapse
Affiliation(s)
- Nissrin 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.,Department of Mathematics and Natural Sciences, School of Arts and Sciences , American University of Ras Al Khaimah , RAK , P.O. Box 10021 , UAE
| | - Pablo Villarreal
- Institute of Fundamental Physics (IFF-CSIC) , CSIC , Serrano 123 , 28006 Madrid , Spain
| | - Rita Prosmiti
- Institute of Fundamental Physics (IFF-CSIC) , CSIC , Serrano 123 , 28006 Madrid , Spain
| |
Collapse
|
8
|
Valdés Á, Carrillo-Bohórquez O, Prosmiti R. Fully Coupled Quantum Treatment of Nanoconfined Systems: A Water Molecule inside a Fullerene C 60. J Chem Theory Comput 2018; 14:6521-6531. [PMID: 30419169 DOI: 10.1021/acs.jctc.8b00801] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023]
Abstract
We implemented a systematic procedure for treating the quantal rotations by including all translational and vibrational degrees of freedom for any triatomic bent molecule in any embedded or confined environment, within the MCTDH framework. Fully coupled quantum treatments were employed to investigate unconventional properties in nanoconfined molecular systems. In this way, we facilitate a complete theoretical analysis of the underlying dynamics that enables us to compute the energy levels and the nuclear spin isomers of a single water molecule trapped in a C60 fullerene cage. The key point lies in the full 9D description of both nuclear and electronic degrees of freedom, as well as a reliable representation of the guest-host interaction. The presence of occluded impurities or inhomogeneities due to noncovalent interactions in the interfullerene environment could modify aspects of the potential, causing significant coupling between otherwise uncoupled modes. Using specific n-mode model potentials, we obtained splitting patterns that confirm the effects of symmetry breaking observed by experiments in the ground ortho-H2O state. Further, our investigation reveals that the first rotationally excited states of the encapsulated ortho- and para-H2O have also raised their 3-fold degeneracy. In view of the complexity of the problem, our results highlight the importance of accurate and computational demanding approaches for building up predictive models for such nanoconfined molecules.
Collapse
Affiliation(s)
- Álvaro Valdés
- Departamento de Física , Universidad Nacional de Colombia , Calle 26, Cra 39, Edicio 404 , Bogotá , Colombia
| | - Orlando Carrillo-Bohórquez
- Departamento de Física , Universidad Nacional de Colombia , Calle 26, Cra 39, Edicio 404 , Bogotá , Colombia
| | - Rita Prosmiti
- Institute of Fundamental Physics (IFF-CSIC), CSIC , Serrano 123 , 28006 Madrid , Spain
| |
Collapse
|
9
|
Valdés Á, Prosmiti R. Preferential stabilization of HeI2 van der Waals isomers: the effect of energetics and temperature. RSC Adv 2017. [DOI: 10.1039/c7ra01378g] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022] Open
Abstract
Energetics and temperature dependence on the preferential stabilization of the linear/T-shaped HeI2 isomers.
Collapse
Affiliation(s)
- Álvaro Valdés
- Departamento de Física
- Universidad Nacional de Colombia
- Bogotá
- Colombia
| | - Rita Prosmiti
- Instituto de Física Fundamental (IFF-CSIC)
- CSIC
- 28006 Madrid
- Spain
| |
Collapse
|
10
|
Carrillo-Bohórquez O, Valdés Á, Prosmiti R. Temperature Dependence of HeBr 2 Isomers' Stability through Rovibrational Multiconfiguration Time-Dependent Hartree Calculations. J Phys Chem A 2016; 120:9458-9464. [PMID: 27802030 DOI: 10.1021/acs.jpca.6b09107] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/27/2022]
Abstract
The multiconfiguration time-dependent Hartree (MCTDH) method using a six-dimensional Hamiltonian that includes all rotational and vibrational degrees of freedom and an ab initio potential energy surface was employed to calculate the rovibronic states of the HeBr2 van der Waals complex. All rotational states of energies within 7 cm -1 with respect to the energy of the linear ground state were calculated without restriction of the total angular momentum. In total, we obtained 500 and 320 rotationally excited states of the ground vibrational T-shaped and linear isomers of the HeBr2, respectively, and compared them with those predicted by the rigid rotor model. A thermodynamic model was then introduced to determine the relative stability of the two conformers as a function of the temperature. On the basis of the present results, the linear conformers were found to be energetically more stable than the T-shaped ones by 1.14 cm-1 at T = 0 K, whereas conversion from linear to T-shaped complexes was observed at temperatures above 2.87 K.
Collapse
Affiliation(s)
- Orlando Carrillo-Bohórquez
- Departamento de Física, Universidad Nacional de Colombia , Calle 26, Cra 39, Edificio 404, Bogotá, Colombia
| | - Álvaro Valdés
- Departamento de Física, Universidad Nacional de Colombia , Calle 26, Cra 39, Edificio 404, Bogotá, Colombia
| | - Rita Prosmiti
- Institute of Fundamental Physics (IFF-CSIC), CSIC , Serrano 123, 28006 Madrid, Spain
| |
Collapse
|
11
|
Rodríguez-Cantano R, González-Lezana T, Villarreal P. Path integral Monte Carlo investigations on doped helium clusters. INT REV PHYS CHEM 2016. [DOI: 10.1080/0144235x.2015.1132595] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/22/2022]
|
12
|
Valdés Á, Prosmiti R. Vibrational Calculations of Higher-Order Weakly Bound Complexes: The He3,4I2 Cases. J Phys Chem A 2015; 119:12736-41. [PMID: 26634405 DOI: 10.1021/acs.jpca.5b10398] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
The structure and relative stability of higher-order He3,4I2 clusters are investigated by carrying out full-dimensional quantum calculations within the multiconfiguration time-dependent Hartree framework. The full interaction between the I2 molecule and the He atoms is based on analytical three-body ab initio He-I2 potentials obtained from high level ab initio calculations plus the He-He interaction. The low-lying minima on the potential surfaces are found to be very close in energy with the He atoms in a ring encircling the dopant for the global minimum structure, while for the local minima one or two of the He atoms prefer the linear arrangements along the I2-axis. Such classical description on the basis of the potential energy landscape is corrected by including anharmonic quantum effects, present in highly floppy systems, by carrying out full dimensional quantum calculations. The potential energy operator was constructed by natural potential fits, while a mode combination scheme was employed to optimize the computational cost of the improved relaxation calculations. The obtained results predict the relative stability of the He3,4I2 isomers at zero temperature and provide benchmark data on binding energies and structural properties of these van der Waals systems. The (2,1) and (2,2), involving two He atoms in the T-shape and one or two He atoms in the linear configurations, respectively, are found to be the most stable isomers, although extremely close in energy with the (3,0) and (4,0) ones as predicted by classical optimizations. Comparison with experimental data on similar systems at low temperatures is also discussed. This analysis indicates once more the importance of quantum delocalization and the need of accurate quantum-mechanical treatments to characterize such doped helium nanosystems.
Collapse
Affiliation(s)
- Álvaro Valdés
- Departamento de Física, Universidad Nacional de Colombia , Calle 26, Cra 39, Edificio 404, Bogotá, Colombia
| | - Rita Prosmiti
- Instituto de Física Fundamental (IFF-CSIC), CSIC , Serrano 123, 28006 Madrid, Spain
| |
Collapse
|
13
|
Rodríguez-Cantano R, González-Lezana T, Villarreal P, Gianturco FA. A configurational study of helium clusters doped with He∗− and He2∗−. J Chem Phys 2015; 142:104303. [DOI: 10.1063/1.4913958] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2023] Open
Affiliation(s)
| | | | - Pablo Villarreal
- Instituto de Física Fundamental, CSIC, Serrano 123, 28006 Madrid, Spain
| | - Franco A. Gianturco
- Institute of Ion Physics, The University of Innsbruck, Technikerstr. 25, 6020 Innsbruck, Austriaand Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy
| |
Collapse
|