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; 25: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
|
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
|
3
|
A Path Integral Molecular Dynamics Simulation of a Harpoon-Type Redox Reaction in a Helium Nanodroplet. Molecules 2021; 26:molecules26195783. [PMID: 34641327 PMCID: PMC8510490 DOI: 10.3390/molecules26195783] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/02/2021] [Revised: 09/14/2021] [Accepted: 09/21/2021] [Indexed: 11/17/2022] Open
Abstract
We present path integral molecular dynamics (PIMD) calculations of an electron transfer from a heliophobic Cs2 dimer in its (3Σu) state, located on the surface of a He droplet, to a heliophilic, fully immersed C60 molecule. Supported by electron ionization mass spectroscopy measurements (Renzler et al., J. Chem. Phys.2016, 145, 181101), this spatially quenched reaction was characterized as a harpoon-type or long-range electron transfer in a previous high-level ab initio study (de Lara-Castells et al., J. Phys. Chem. Lett.2017, 8, 4284). To go beyond the static approach, classical and quantum PIMD simulations are performed at 2 K, slightly below the critical temperature for helium superfluidity (2.172 K). Calculations are executed in the NVT ensemble as well as the NVE ensemble to provide insights into real-time dynamics. A droplet size of 2090 atoms is assumed to study the impact of spatial hindrance on reactivity. By changing the number of beads in the PIMD simulations, the impact of quantization can be studied in greater detail and without an implicit assumption of superfluidity. We find that the reaction probability increases with higher levels of quantization. Our findings confirm earlier, static predictions of a rotational motion of the Cs2 dimer upon reacting with the fullerene, involving a substantial displacement of helium. However, it also raises the new question of whether the interacting species are driven out-of-equilibrium after impurity uptake, since reactivity is strongly quenched if a full thermal equilibration is assumed. More generally, our work points towards a novel mechanism for long-range electron transfer through an interplay between nuclear quantum delocalization within the confining medium and delocalized electronic dispersion forces acting on the two reactants.
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: 5] [Impact Index Per Article: 1.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
|
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
|
6
|
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: 9] [Impact Index Per Article: 1.8] [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
|
7
|
Hauser AW, de Lara-Castells MP. Spatial quenching of a molecular charge-transfer process in a quantum fluid: the Cs x-C 60 reaction in superfluid helium nanodroplets. Phys Chem Chem Phys 2017; 19:1342-1351. [PMID: 27975088 DOI: 10.1039/c6cp06858h] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
A recent experimental study [Renzler et al., J. Chem. Phys., 2016, 145, 181101] on superfluid helium nanodroplets reported different reactivities for Cs atoms and Cs2 dimers with C60 fullerenes inside helium droplets. Alkali metal atoms and clusters are heliophobic, therefore typically residing on the droplet surface, while fullerenes are fully immersed into the droplet. In this theoretical study, which combines standard methods of computational chemistry with orbital-free helium density functional theory, we show that the experimental findings can be interpreted in the light of a quenched electron-transfer reaction between the fullerene and the alkali dopant, which is additionally hindered by a reaction barrier stemming from the necessary extrusion of helium upon approach of the two reactants.
Collapse
Affiliation(s)
- Andreas W Hauser
- Graz University of Technology, Institute of Experimental Physics, Petersgasse 16, 8010 Graz, Austria.
| | | |
Collapse
|
8
|
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]
|
9
|
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
|
10
|
Villarreal P, Rodríguez-Cantano R, González-Lezana T, Prosmiti R, Delgado-Barrio G, Gianturco FA. Quantum Features of Anionic Species He*⁻ and He₂*⁻ in Small He(N) Clusters. J Phys Chem A 2015; 119:11574-82. [PMID: 26517581 DOI: 10.1021/acs.jpca.5b08470] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
We present variational calculations on systems containing a few boson helium atoms attached to electronically excited atomic and molecular helium anions He*⁻ and He₂*⁻ and characterize their structures and energetics. Previously reported high-level ab initio results [Huber, S. E.; Mauracher, A. Mol. Phys. 2014, 112, 794] to describe the interactions between excited (metastable) anions and a neutral He atom have been employed. For the case of the atomic species He*⁻, the corresponding interaction with He suggests large anharmonicity effects due to the presence of a deep well of ∼17,500 cm⁻¹ at short distances, together with a more external shallow secondary well of ∼4 cm⁻¹, both supporting bound levels. Moreover, when a sum of pairwise interactions is assumed to describe the full PES corresponding to the presence of several neutral He atoms, geometrical constraints already predict the complete solvation of the anionic impurity by six helium atoms, giving rise to a bipyramidal structure. In turn, for the anisotropic weak interaction He-He₂*⁻, where the anionic dimer is considered as a rigid rotor, the obtained structures show the tendency of the helium atoms to pack themselves together and largely far away from the dopant, thereby confirming the heliophobic character of He₂*⁻.
Collapse
Affiliation(s)
- Pablo Villarreal
- Instituto de Fı́sica Fundamental, IFF-CSIC , Serrano 123, 28006 Madrid, Spain
| | | | | | - Rita Prosmiti
- Instituto de Fı́sica Fundamental, IFF-CSIC , Serrano 123, 28006 Madrid, Spain
| | | | - Franco A Gianturco
- Institute of Ion Physics, The University of Innsbruck , Technikerstrasse 25, 6020 Innsbruck, Austria.,Scuola Normale Superiore ,Piazza dei Cavalieri 7, 56126 Pisa, Italy
| |
Collapse
|
11
|
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
|
12
|
Albertí M, Aguilar A, Pirani F. Propensities in the solvation of M+–Benzene systems (M=Na, K, Rb) investigated by cluster dynamics. Chem Phys 2012. [DOI: 10.1016/j.chemphys.2011.07.030] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/27/2022]
|
13
|
Rodríguez-Cantano R, López-Durán D, González-Lezana T, Delgado-Barrio G, Villarreal P, Yurtsever E, Gianturco FA. Spin-polarized Rb2 interacting with bosonic He atoms: potential energy surface and quantum structures of small clusters. J Phys Chem A 2012; 116:2394-404. [PMID: 22320704 DOI: 10.1021/jp3004932] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Abstract
A new full-dimension potential energy surface of the three-body He-Rb₂(³Σ(u)(+)) complex and a quantum study of small (⁴He)(N)-Rb₂(³Σ(u)(+)) clusters, 1 ≤ N ≤ 4, are presented. We have accurately fitted the ab initio points of the interaction to an analytical form and addressed the dopant's vibration, which is found to be negligible. A Variational approach and a Diffusion Monte Carlo technique have been applied to yield energy and geometric properties of the selected species. Our quantum structure calculations show a transition in the arrangements of the helium atoms from N = 2, where they tend to be separated across the diatomic bond, to N = 4, in which a closer packing of the rare gas particles is reached, guided by the dominance of the He-He potential over the weaker interaction of the latter adatoms with the doping dimer. The deepest well of the He-Rb₂ interaction is placed at the T-shape configuration, a feature which causes the dopant to be located as parallel to the helium "minidroplet". Our results are shown to agree with previous findings on this and on similar systems.
Collapse
|
14
|
Pérez de Tudela R, López-Durán D, González-Lezana T, Delgado-Barrio G, Villarreal P, Gianturco FA, Yurtsever E. Quantum features of a barely bound molecular dopant: Cs2(3Σu) in bosonic helium droplets of variable size. J Phys Chem A 2011; 115:6892-902. [PMID: 21585200 DOI: 10.1021/jp111825n] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
We present in this work the study of small (4)He(N)-Cs(2)((3)Σ(u)) aggregates (2 ≤ N ≤ 30) through combined variational, diffusion Monte Carlo (DMC), and path integral Monte Carlo (PIMC) calculations. The full surface is modeled as an addition of He-Cs(2) interactions and He-He potentials. Given the negligible strength and large range of the He-Cs(2) interaction as compared with the one for He-He, a propensity of the helium atoms to pack themselves together, leaving outside the molecular dopant is to be expected. DMC calculations determine the onset of helium gathering at N = 3. To analyze energetic and structural properties as a function of N, PIMC calculations with no bosonic exchange, i.e., Boltzmann statistics, at low temperatures are carried out. At T = 0.1 K, although acceptable one-particle He-Cs(2) distributions are obtained, two-particle He-He distributions are not well described, indicating that the proper symmetry should be taken into account. PIMC distributions at T = 1 K already compare well with DMC ones and show minor exchange effects, although binding energies are still far from having converged in terms of the number of quantum beads. As N increases, the He-He PIMC pair correlation function shows a clear tendency to coincide with the experimental boson-liquid helium one at that temperature. It supports the picture of a helium droplet which carries the molecular impurity on its surface, as found earlier for other triplet dimers.
Collapse
|
15
|
Guillon G, Zanchet A, Leino M, Viel A, Zillich RE. Theoretical Study of Rb2 in HeN: Potential Energy Surface and Monte Carlo Simulations. J Phys Chem A 2011; 115:6918-26. [DOI: 10.1021/jp112053b] [Citation(s) in RCA: 24] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Grégroire Guillon
- Institut de Physique de Rennes, UMR 6251, CNRS & Université de Rennes I, F-35042 Rennes, France
| | | | | | | | | |
Collapse
|