1
|
Slama M, Habli H, Laajimi M, Ghalla H, Ben El Hadj Rhouma M. Microsolvation of lithium cation in xenon clusters: An octahedral growth pattern. J Mol Graph Model 2022; 116:108229. [PMID: 35671571 DOI: 10.1016/j.jmgm.2022.108229] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/19/2022] [Revised: 05/09/2022] [Accepted: 05/19/2022] [Indexed: 01/18/2023]
Abstract
The structural and energetic proprieties for the Li + Xen (n = 1-18) clusters are investigated using both Basin-Hopping combined with Potential Model description (BH-PM) and DFT methods. A structural transition from tetrahedral (4 coordination) form to octahedral (6 coordination) one is observed for n = 6. Above this size, all structures have an octahedral core. The cubic-face-centered arrangement for xenon atoms is detected for Li + Xe14. To the best of our knowledge, the Li + Xen (n = 1-18) clusters are studied in the present work for the first time using the DFT theoretical approach. The M062X functional combined with aug-cc-pVDZ (for Li) and def2-TZVP (for Xe) basis sets reproduces accurately the CCSD(T) potential energy curve of Li + Xe system. Atom-Centered Density Matrix Propagation (ADMP) molecular dynamic calculations have been carried. Moreover, we investigate the larger sizes n = 31-35, 44, and 55 for the first time using the BH-PM theoretical approach. The closing of the first and second octahedron shells are proved for the n = 6 and 34 sizes, respectively. The relative stabilities of the Li + Xen molecules are also studied by computing the total energy, the binding energy per atoms for each size n. Then, the second energy difference between the size n and its two near neighbors allows identifying the magic number series. Our present data are analyzed, discussed and compared with the available theoretical and experimental data.
Collapse
Affiliation(s)
- Marwa Slama
- Université de Monastir, Institut Préparatoire aux Études des Ingénieurs de Monastir, Laboratoire d'Études des Milieux Ionisés et Réactifs (EMIR), 5000 Monastir, Tunisia.
| | - Hela Habli
- Université de Monastir, Faculté des Sciences de Monastir, Laboratoire de Physique Quantique et Statistique, Avenue de l'Environnement 5019 Monastir, Tunisia; Université de Sousse, Institut Supérieur des Sciences Appliquées et de Technologie de Sousse, Rue ibn Khaldun, Cité Taffala, 4003 Sousse, Tunisia
| | - Maha Laajimi
- Université de Monastir, Faculté des Sciences de Monastir, Laboratoire de Physique Quantique et Statistique, Avenue de l'Environnement 5019 Monastir, Tunisia
| | - Houcine Ghalla
- Université de Monastir, Faculté des Sciences de Monastir, Laboratoire de Physique Quantique et Statistique, Avenue de l'Environnement 5019 Monastir, Tunisia
| | - Mounir Ben El Hadj Rhouma
- Université de Monastir, Institut Préparatoire aux Études des Ingénieurs de Monastir, Laboratoire d'Études des Milieux Ionisés et Réactifs (EMIR), 5000 Monastir, Tunisia
| |
Collapse
|
2
|
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
|
3
|
Andrade MDD, Jesus WS, Prudente FV, Marques JMC. On the stabilization of the Li$$^+$$-Li$$^+$$ interaction by microsolvation with rare-gas atoms. Theor Chem Acc 2021. [DOI: 10.1007/s00214-021-02763-8] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/21/2022]
|
4
|
Jesus WS, Prudente FV, Marques JMC, Pereira FB. Modeling microsolvation clusters with electronic-structure calculations guided by analytical potentials and predictive machine learning techniques. Phys Chem Chem Phys 2021; 23:1738-1749. [PMID: 33427847 DOI: 10.1039/d0cp05200k] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/16/2023]
Abstract
We propose a new methodology to study, at the density functional theory (DFT) level, the clusters resulting from the microsolvation of alkali-metal ions with rare-gas atoms. The workflow begins with a global optimization search to generate a pool of low-energy minimum structures for different cluster sizes. This is achieved by employing an analytical potential energy surface (PES) and an evolutionary algorithm (EA). The next main stage of the methodology is devoted to establish an adequate DFT approach to treat the microsolvation system, through a systematic benchmark study involving several combinations of functionals and basis sets, in order to characterize the global minimum structures of the smaller clusters. In the next stage, we apply machine learning (ML) classification algorithms to predict how the low-energy minima of the analytical PES map to the DFT ones. An early and accurate detection of likely DFT local minima is extremely important to guide the choice of the most promising low-energy minima of large clusters to be re-optimized at the DFT level of theory. In this work, the methodology was applied to the Li+Krn (n = 2-14 and 16) microsolvation clusters for which the most competitive DFT approach was found to be the B3LYP-D3/aug-pcseg-1. Additionally, the ML classifier was able to accurately predict most of the solutions to be re-optimized at the DFT level of theory, thereby greatly enhancing the efficiency of the process and allowing its applicability to larger clusters.
Collapse
Affiliation(s)
- W S Jesus
- Instituto de Física, Universidade Federal da Bahia, 40170-115 Salvador, BA, Brazil.
| | - F V Prudente
- Instituto de Física, Universidade Federal da Bahia, 40170-115 Salvador, BA, Brazil.
| | - J M C Marques
- CQC, Department of Chemistry, University of Coimbra, 3004-535 Coimbra, Portugal.
| | - F B Pereira
- Coimbra Polytechnic - ISEC, Coimbra, Portugal and Centro de Informática e Sistemas da Universidade de Coimbra (CISUC), Coimbra, Portugal.
| |
Collapse
|
5
|
Guimarães MN, M. de Almeida M, Marques JMC, Prudente FV. A thermodynamic view on the microsolvation of ions by rare gas: application to Li+ with argon. Phys Chem Chem Phys 2020; 22:10882-10892. [DOI: 10.1039/d0cp01283a] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Parallel tempering Monte Carlo calculations on the Li+Arn microsolvation clusters have shown that the two peaks appearing in the heat capacity curve as a function of temperature correspond to the melting of the second and first solvation shells.
Collapse
Affiliation(s)
- M. N. Guimarães
- Instituto de Física
- Universidade Federal da Bahia
- 40170-115 Salvador
- Brazil
| | - M. M. de Almeida
- Instituto de Física
- Universidade Federal da Bahia
- 40170-115 Salvador
- Brazil
| | - J. M. C. Marques
- CQC
- Department of Chemistry
- University of Coimbra
- 3004-535 Coimbra
- Portugal
| | - F. V. Prudente
- Instituto de Física
- Universidade Federal da Bahia
- 40170-115 Salvador
- Brazil
| |
Collapse
|
6
|
Ortiz de Zárate J, Bartolomei M, González-Lezana T, Campos-Martínez J, Hernández MI, Pérez de Tudela R, Hernández-Rojas J, Bretón J, Pirani F, Kranabetter L, Martini P, Kuhn M, Laimer F, Scheier P. Snowball formation for Cs + solvation in molecular hydrogen and deuterium. Phys Chem Chem Phys 2019; 21:15662-15668. [PMID: 31271179 DOI: 10.1039/c9cp02017a] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/18/2023]
Abstract
Interactions of atomic cations with molecular hydrogen are of interest for a wide range of applications in hydrogen technologies. These interactions are fairly strong despite being non-covalent, hence one can ask whether hydrogen molecules would form dense, solid-like, solvation shells around the ion (snowballs) or rather a more weakly bound compound. In this work, the interactions between Cs+ and H2 are studied both experimentally and computationally. Isotopic substitution of H2 by D2 is also investigated. On the one hand, helium nanodroplets doped with cesium and hydrogen or deuterium are ionized by electron impact and the (H2/D2)nCs+ (up to n = 30) clusters formed are identified via mass spectrometry. On the other hand, a new analytical potential energy surface, based on ab initio calculations, is developed and used to study cluster energies and structures by means of classical and quantum-mechanical Monte Carlo methods. The most salient features of the measured ion abundances are remarkably mimicked by the computed evaporation energies, particularly for the clusters composed of deuterium. This result supports the reliability of the present potential energy surface and allows us to recommend its use in related systems. Clusters with either twelve H2 or D2 molecules stand out for their stability and quasi-rigid icosahedral structures. However, the first solvation shell involves thirteen or fourteen molecules for hydrogenated or deuterated clusters, respectively. This shell retains its internal structure when extra molecules are added to the second shell and is nearly solid-like, especially for the deuterated clusters. The role played by three-body induction interactions as well as the rotational degrees of freedom is analyzed and they are found to be significant (up to 15% and 18%, respectively) for the molecules belonging to the first solvation shell.
Collapse
Affiliation(s)
- Josu Ortiz de Zárate
- Instituto de Física Fundamental, Consejo Superior de Investigaciones Científicas (IFF-CSIC), Serrano 123, 28006 Madrid, Spain.
| | - Massimiliano Bartolomei
- Instituto de Física Fundamental, Consejo Superior de Investigaciones Científicas (IFF-CSIC), Serrano 123, 28006 Madrid, Spain.
| | - Tomás González-Lezana
- Instituto de Física Fundamental, Consejo Superior de Investigaciones Científicas (IFF-CSIC), Serrano 123, 28006 Madrid, Spain.
| | - José Campos-Martínez
- Instituto de Física Fundamental, Consejo Superior de Investigaciones Científicas (IFF-CSIC), Serrano 123, 28006 Madrid, Spain.
| | - Marta I Hernández
- Instituto de Física Fundamental, Consejo Superior de Investigaciones Científicas (IFF-CSIC), Serrano 123, 28006 Madrid, Spain.
| | | | - Javier Hernández-Rojas
- Departamento de Fsica and IUdEA, Universidad de La Laguna, 38205, La Laguna, Tenerife, Spain
| | - José Bretón
- Departamento de Fsica and IUdEA, Universidad de La Laguna, 38205, La Laguna, Tenerife, Spain
| | - Fernando Pirani
- Dipartimento di Chimica, Biologia e Biotecnologie, Università di Perugia, Perugia, Italy
| | - Lorenz Kranabetter
- Universität Innsbruck, Institut für Ionenphyisk und Angewandte Physik, Technikerstraße 25, 6020 Innsbruck, Austria
| | - Paul Martini
- Universität Innsbruck, Institut für Ionenphyisk und Angewandte Physik, Technikerstraße 25, 6020 Innsbruck, Austria
| | - Martin Kuhn
- Universität Innsbruck, Institut für Ionenphyisk und Angewandte Physik, Technikerstraße 25, 6020 Innsbruck, Austria
| | - Felix Laimer
- Universität Innsbruck, Institut für Ionenphyisk und Angewandte Physik, Technikerstraße 25, 6020 Innsbruck, Austria
| | - Paul Scheier
- Universität Innsbruck, Institut für Ionenphyisk und Angewandte Physik, Technikerstraße 25, 6020 Innsbruck, Austria
| |
Collapse
|