1
|
García-Alfonso E, Barranco M, Bonhommeau DA, Halberstadt N, Pi M, Calvo F. Clustering, collision, and relaxation dynamics in pure and doped helium nanoclusters: Density- vs particle-based approaches. J Chem Phys 2022; 157:014106. [DOI: 10.1063/5.0091942] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The clustering, collision, and relaxation dynamics of pristine and doped helium nanodroplets is theoretically investigated in cases of pickup and clustering of heliophilic argon, collision of heliophobic cesium atoms, and coalescence of two droplets brought into contact by their mutual long-range van der Waals interaction. Three approaches are used and compared with each other. The He time-dependent density functional theory method considers the droplet as a continuous medium and accounts for its superfluid character. The ring-polymer molecular dynamics method uses a path-integral description of nuclear motion and incorporates zero-point delocalization while bosonic exchange effects are ignored. Finally, the zero-point averaged dynamics approach is a mixed quantum–classical method in which quantum delocalization is described by attaching a frozen wavefunction to each He atom, equivalent to classical dynamics with effective interaction potentials. All three methods predict that the growth of argon clusters is significantly hindered by the helium host droplet due to the impeding shell structure around the dopants and kinematic effects freezing the growing cluster in metastable configurations. The effects of superfluidity are qualitatively manifested by different collision dynamics of the heliophilic atom at high velocities, as well as quadrupole oscillations that are not seen with particle-based methods, for droplets experiencing a collision with cesium atoms or merging with each other.
Collapse
Affiliation(s)
- Ernesto García-Alfonso
- Laboratoire Collisions, Agrégats, Réactivité (LCAR), Université de Toulouse, CNRS, 31062 Toulouse, France
| | - Manuel Barranco
- Laboratoire Collisions, Agrégats, Réactivité (LCAR), Université de Toulouse, CNRS, 31062 Toulouse, France
- Department FQA, Facultat de Física, Universitat de Barcelona, Diagonal 645, 08028 Barcelona, Spain
- Institute of Nanoscience and Nanotechnology (IN2UB), Universitat de Barcelona, Barcelona, Spain
| | - David A. Bonhommeau
- Université de Reims Champagne Ardenne, CNRS, GSMA UMR 7331, 51100 Reims, France
| | - Nadine Halberstadt
- Laboratoire Collisions, Agrégats, Réactivité (LCAR), Université de Toulouse, CNRS, 31062 Toulouse, France
| | - Martí Pi
- Department FQA, Facultat de Física, Universitat de Barcelona, Diagonal 645, 08028 Barcelona, Spain
- Institute of Nanoscience and Nanotechnology (IN2UB), Universitat de Barcelona, Barcelona, Spain
| | - Florent Calvo
- Université Grenoble Alpes, CNRS, LIPHY, F38000 Grenoble, France
| |
Collapse
|
2
|
Ellis AM, Davies JA, Yurtsever E, Calvo F. Dimerization dynamics of carboxylic acids in helium nanodroplets. J Chem Phys 2022; 156:174304. [PMID: 35525638 DOI: 10.1063/5.0087957] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The dimerization of molecules in helium nanodroplets is known to preferentially yield structures of higher energy than the global energy minimum structure for a number of quite different monomers. Here, we explore dimerization in this environment using an atomistic model within statistically converged molecular dynamics (MD) trajectories, treating the solvent implicitly through the use of a thermostat, or more explicitly by embedding one monomer in a He100 cluster. The focus is on the two simplest carboxylic acids, formic and acetic, both of which have been studied experimentally. While the global minimum structure, which comprises two CO⋯HO hydrogen bonds, is predicted to be the most abundant dimer in the absence of the helium solvent, this is no longer the case once helium atoms are included. The simulations confirm the importance of kinetic trapping effects and also shed light on the occurrence of specific dynamical effects, leading to the occasional formation of high-energy structures away from minima, such as saddle configurations. Theoretically predicted infrared spectra, based on the MD statistics, are in good agreement with the experimental spectra.
Collapse
Affiliation(s)
- Andrew M Ellis
- School of Chemistry, University of Leicester, University Road, Leicester LE1 7RH, United Kingdom
| | - Julia A Davies
- School of Chemistry, University of Leicester, University Road, Leicester LE1 7RH, United Kingdom
| | - Ersin Yurtsever
- Department of Chemistry, Koç University, Rumelifeneri Yolu, Sariyer, 34450 Istanbul, Turkey
| | - Florent Calvo
- Université Grenoble Alpes, CNRS, LiPhy, F38000 Grenoble, France
| |
Collapse
|
3
|
|
4
|
|
5
|
Dontot L, Spiegelman F, Rapacioli M. Structures and Energetics of Neutral and Cationic Pyrene Clusters. J Phys Chem A 2019; 123:9531-9543. [PMID: 31589446 PMCID: PMC6917508 DOI: 10.1021/acs.jpca.9b07007] [Citation(s) in RCA: 17] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/31/2023]
Abstract
The low energy structures of neutral and cationic pyrene clusters containing up to seven molecules are searched through a global exploration scheme combining parallel tempering Monte Carlo algorithm and local quenches. The potential energies are computed at the density functional based tight binding level for neutrals and configuration interaction density functional based tight binding for cations in order to treat properly the charge resonance. New simplified versions of these schemes are also presented and used during the global exploration. Neutral clusters are shown to be made of compact assemblies of sub-blocs containing up to three units whereas cations present a charged dimer or trimer core surrounded by neutral units. The structural features of the clusters are analyzed and correlated for the cation with the charge distribution. The stability of clusters is also discussed in terms of cohesive and evaporation energies. Adiabatic and vertical ionization potentials are also discussed.
Collapse
Affiliation(s)
- Léo Dontot
- Laboratoire de Chimie et Physique Quantiques LCPQ/IRSAMC, UMR5626, Université de Toulouse (UPS) and CNRS, 118 Route de Narbonne, F-31062 Toulouse, France
| | - Fernand Spiegelman
- Laboratoire de Chimie et Physique Quantiques LCPQ/IRSAMC, UMR5626, Université de Toulouse (UPS) and CNRS, 118 Route de Narbonne, F-31062 Toulouse, France
| | - Mathias Rapacioli
- Laboratoire de Chimie et Physique Quantiques LCPQ/IRSAMC, UMR5626, Université de Toulouse (UPS) and CNRS, 118 Route de Narbonne, F-31062 Toulouse, France
| |
Collapse
|
6
|
Schouder C, Chatterley AS, Calvo F, Christiansen L, Stapelfeldt H. Structure determination of the tetracene dimer in helium nanodroplets using femtosecond strong-field ionization. STRUCTURAL DYNAMICS (MELVILLE, N.Y.) 2019; 6:044301. [PMID: 31463336 PMCID: PMC6711753 DOI: 10.1063/1.5118005] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/02/2019] [Accepted: 08/09/2019] [Indexed: 05/29/2023]
Abstract
Dimers of tetracene molecules are formed inside helium nanodroplets and identified through covariance analysis of the emission directions of kinetic tetracene cations stemming from femtosecond laser-induced Coulomb explosion. Next, the dimers are aligned in either one or three dimensions under field-free conditions by a nonresonant, moderately intense laser pulse. The experimental angular covariance maps of the tetracene ions are compared to calculated covariance maps for seven different dimer conformations and found to be consistent with four of these. Additional measurements of the alignment-dependent strong-field ionization yield of the dimer narrow the possible conformations down to either a slipped-parallel or parallel-slightly rotated structure. According to our quantum chemistry calculations, these are the two most stable gas-phase conformations of the dimer and one of them is favorable for singlet fission.
Collapse
Affiliation(s)
- Constant Schouder
- Department of Physics and Astronomy, Aarhus University, Ny Munkegade 120, 8000 Aarhus C, Denmark
| | - Adam S Chatterley
- Department of Chemistry, Aarhus University, Langelandsgade 140, 8000 Aarhus C, Denmark
| | - Florent Calvo
- Université Grenoble Alpes CNRS, LIPHY, F-38000 Grenoble, France
| | - Lars Christiansen
- Department of Chemistry, Aarhus University, Langelandsgade 140, 8000 Aarhus C, Denmark
| | - Henrik Stapelfeldt
- Department of Chemistry, Aarhus University, Langelandsgade 140, 8000 Aarhus C, Denmark
| |
Collapse
|
7
|
Goulart M, Kuhn M, Rasul B, Postler J, Gatchell M, Zettergren H, Scheier P, Echt O. The structure of coronene cluster ions inferred from H 2 uptake in the gas phase. Phys Chem Chem Phys 2018; 19:27968-27973. [PMID: 29022968 DOI: 10.1039/c7cp04999d] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Mass spectra of helium nanodroplets doped with H2 and coronene feature anomalies in the ion abundance that reveal anomalies in the energetics of adsorption sites. The coronene monomer ion strongly adsorbs up to n = 38 H2 molecules indicating a commensurate solvation shell that preserves the D6h symmetry of the substrate. No such feature is seen in the abundance of the coronene dimer through tetramer complexed with H2; this observation rules out a vertical columnar structure. Instead we see evidence for a columnar structure in which adjacent coronenes are displaced in parallel, forming terraces that offer additional strong adsorption sites. The experimental value for the number of adsorption sites per terrace, approximately six, barely depends on the number of coronene molecules. The displacement estimated from this number exceeds the value reported in several theoretical studies of the bare, neutral coronene dimer.
Collapse
Affiliation(s)
- Marcelo Goulart
- Institut für Ionenphysik und Angewandte Physik, Universität Innsbruck, Technikerstr. 25, A-6020 Innsbruck, Austria.
| | | | | | | | | | | | | | | |
Collapse
|
8
|
Seki Y, Takayanagi T, Shiga M. Photoexcited Ag ejection from a low-temperature He cluster: a simulation study by nonadiabatic Ehrenfest ring-polymer molecular dynamics. Phys Chem Chem Phys 2017; 19:13798-13806. [DOI: 10.1039/c7cp00888k] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Nonadiabatic ring-polymer molecular dynamics simulations were performed to understand the photoexcitation dynamics of a low-temperature Ag·He500 cluster.
Collapse
Affiliation(s)
- Yusuke Seki
- Department of Chemistry
- Saitama University
- Saitama City
- Japan
| | | | - Motoyuki Shiga
- Center for Computational Science and E-Systems
- Japan Atomic Energy Agency
- Kashiwa
- Japan
| |
Collapse
|
9
|
Silva AF, Duarte LJ, Bruns RE. QTAIM-Based Characteristic Group Infrared Intensities of Amino Acids and Their Transference to Peptides. J Phys Chem A 2016; 120:8387-8399. [DOI: 10.1021/acs.jpca.6b07690] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Arnaldo F. Silva
- Institute
of Chemistry, University of Campinas, UNICAMP, Cidade Universitária Campinas, PO Box 6154, Campinas, São Paulo 13083970, Brazil
| | - Leonardo J. Duarte
- Institute
of Chemistry, University of Campinas, UNICAMP, Cidade Universitária Campinas, PO Box 6154, Campinas, São Paulo 13083970, Brazil
| | - Roy E. Bruns
- Institute
of Chemistry, University of Campinas, UNICAMP, Cidade Universitária Campinas, PO Box 6154, Campinas, São Paulo 13083970, Brazil
| |
Collapse
|
10
|
Calvo F, Yurtsever E. Solvation of carbonaceous molecules by para-H2 and ortho-D2 clusters. II. Fullerenes. J Chem Phys 2016; 145:084304. [PMID: 27586919 DOI: 10.1063/1.4961159] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The coating of various fullerenes by para-hydrogen and ortho-deuterium molecules has been computationally studied as a function of the solvent amount. Rotationally averaged interaction potentials for structureless hydrogen molecules are employed to model their interaction with neutral or charged carbonaceous dopants containing between 20 and 240 atoms, occasionally comparing different fullerenes having the same size but different shapes. The solvation energy and the size of the first solvation shell obtained from path-integral molecular dynamics simulations at 2 K show only minor influence on the dopant charge and on the possible deuteration of the solvent, although the shell size is largest for ortho-D2 coating cationic fullerenes. Nontrivial finite size effects have been found with the shell size varying non-monotonically close to its completion limit. For fullerenes embedded in large hydrogen clusters, the shell size and solvation energy both follow linear scaling with the fullerene size. The shell sizes obtained for C60 (+) and C70 (+) are close to 49 and 51, respectively, and agree with mass spectrometry experiments.
Collapse
Affiliation(s)
- F Calvo
- Université Grenoble Alpes, LIPHY, F-38000 Grenoble, France and CNRS, LIPHy, F-38000 Grenoble, France
| | - E Yurtsever
- Koç University, Rumelifeneriyolu, Sariyer, Istanbul 34450, Turkey
| |
Collapse
|
11
|
Calvo F, Yurtsever E. Solvation of carbonaceous molecules by para-H2 and ortho-D2 clusters. I. Polycyclic aromatic hydrocarbons. J Chem Phys 2016; 144:224302. [DOI: 10.1063/1.4952957] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022] Open
Affiliation(s)
- F. Calvo
- Univ. Grenoble Alpes, LIPHY, F-38000 Grenoble, France and CNRS, LIPHY, F-38000 Grenoble, France
| | - E. Yurtsever
- Koç University, Rumelifeneriyolu, Sariyer, Istanbul 34450, Turkey
| |
Collapse
|