1
|
Tajouo Tela H, Quintas-Sánchez E, Dubernet ML, Scribano Y, Dawes R, Gatti F, Ndengué S. Rovibrational states calculations of the H 2O-HCN heterodimer with the multiconfiguration time dependent Hartree method. Phys Chem Chem Phys 2023; 25:31813-31824. [PMID: 37966067 DOI: 10.1039/d3cp03225f] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2023]
Abstract
Water and hydrogen cyanide are two of the most common species in space and the atmosphere with the ability of binding to form dimers such as H2O-HCN. In the literature, while calculations characterizing various properties of the H2O-HCN cluster (equilibrium distance, vibrational frequencies and rotational constants) have been done in the past, extensive calculations of the rovibrational states of this system using a reliable quantum dynamical approach have yet to be reported. In this work, we intend to mend that by performing the first calculation of the rovibrational states of the H2O-HCN van der Waals complex on a recently developed potential energy surface. We use the block improved relaxation procedure implemented in the Heidelberg MultiConfiguration Time-Dependent Hartree (MCTDH) package to compute the states of the H2O-HCN isomer, from which we extract the transition frequencies and rotational constants of the complex. We further adapt an approach first suggested by Wang and Carrington-and supported here by analysis routines of the Heidelberg MCTDH package-to properly characterize the computed rovibrational states. The subsequent assignment of rovibrational states was done by theoretical analysis and visual inspection of the wavefunctions. Our simulations provide a Zero Point Energy (ZPE) and intermolecular vibrational frequencies in good agreement with past ab initio calculations. The transition frequencies and rotational constants obtained from our simulations match well with the available experimental data. This work has the broad aim to propose the MCTDH approach as a reliable option to compute and characterize rovibrational states of van der Waals complexes such as the current one.
Collapse
Affiliation(s)
- Hervé Tajouo Tela
- ICTP-East African Institute for Fundamental Research, University of Rwanda, Kigali, Rwanda.
| | - Ernesto Quintas-Sánchez
- Department of Chemistry, Missouri University of Science and Technology, 65409 Rolla, Missouri, USA
| | - Marie-Lise Dubernet
- LERMA, Observatoire de Paris, PSL Research University, CNRS, Sorbonne University, UPMC Univ Paris 06, 75014 Paris, France
| | - Yohann Scribano
- Laboratoire Univers et Particules de Montpellier, UMR-CNRS 5299, Université de Montpellier, Place Eugène Bataillon, 34095 Montpellier, France
| | - Richard Dawes
- Department of Chemistry, Missouri University of Science and Technology, 65409 Rolla, Missouri, USA
| | - Fabien Gatti
- Institut de Sciences Moleculaires d'Orsay, UMR 8214, Université Paris-Sud - Université Paris-Saclay, 91405 Orsay, France
| | - Steve Ndengué
- ICTP-East African Institute for Fundamental Research, University of Rwanda, Kigali, Rwanda.
- Department of Physics, Trinity College, 06106 Hartford, Connecticut, USA
| |
Collapse
|
2
|
Dupuy L, Talotta F, Agostini F, Lauvergnat D, Poirier B, Scribano Y. Adiabatic and Nonadiabatic Dynamics with Interacting Quantum Trajectories. J Chem Theory Comput 2022; 18:6447-6462. [DOI: 10.1021/acs.jctc.2c00744] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022]
Affiliation(s)
- Lucien Dupuy
- Laboratoire Univers et Particules de Montpellier, UMR-CNRS 5299, Université de Montpellier, Place Eugène Bataillon, 34095Montpellier, France
| | - Francesco Talotta
- Université Paris-Saclay, CNRS, Institut de Chimie Physique, UMR-CNRS 8000, 91405Orsay, France
| | - Federica Agostini
- Université Paris-Saclay, CNRS, Institut de Chimie Physique, UMR-CNRS 8000, 91405Orsay, France
| | - David Lauvergnat
- Université Paris-Saclay, CNRS, Institut de Chimie Physique, UMR-CNRS 8000, 91405Orsay, France
| | - Bill Poirier
- Department of Chemistry and Biochemistry, and Department of Physics, Texas Tech University, Box 41061, 79409-1061Lubbock, Texas, United States
| | - Yohann Scribano
- Laboratoire Univers et Particules de Montpellier, UMR-CNRS 5299, Université de Montpellier, Place Eugène Bataillon, 34095Montpellier, France
| |
Collapse
|
3
|
Chen A, Benoit DM, Scribano Y, Nauts A, Lauvergnat D. Smolyak Algorithm Adapted to a System-Bath Separation: Application to an Encapsulated Molecule with Large-Amplitude Motions. J Chem Theory Comput 2022; 18:4366-4372. [PMID: 35584357 DOI: 10.1021/acs.jctc.2c00108] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
A Smolyak algorithm adapted to system-bath separation is proposed for rigorous quantum simulations. This technique combines a sparse grid method with the system-bath concept in a specific configuration without limitations on the form of the Hamiltonian, thus achieving a highly efficient convergence of the excitation transitions for the "system" part. Our approach provides a general way to overcome the perennial convergence problem for the standard Smolyak algorithm and enables the simulation of floppy molecules with more than a hundred degrees of freedom. The efficiency of the present method is illustrated on the simulation of H2 caged in an sII clathrate hydrate including two kinds of cage modes. The transition energies are converged by increasing the number of normal modes of water molecules. Our results confirm the triplet splittings of both translational and rotational (j = 1) transitions of the H2 molecule. Furthermore, they show a slight increase of the translational transitions with respect to the ones in a rigid cage.
Collapse
Affiliation(s)
- Ahai Chen
- Maison de la Simulation, UVSQ, CNRS, CEA, Université Paris-Saclay, 91191 Gif-sur-Yvette, France.,Institut de Chimie Physique, UMR-CNRS 8000, Université Paris-Saclay, CNRS, 91405 Orsay, France
| | - David M Benoit
- E.A. Milne Centre for Astrophysics, The University of Hull, Cottingham Road, Kingston upon Hull HU6 7RX, U.K
| | - Yohann Scribano
- Laboratoire Univers et Particules de Montpellier, UMR-CNRS 5299, Université de Montpellier, 34095 Montpellier Cedex, France
| | - André Nauts
- Institut de Chimie Physique, UMR-CNRS 8000, Université Paris-Saclay, CNRS, 91405 Orsay, France.,Institute of Condensed Matter and Nanosciences (NAPS), Université Catholique de Louvain, 1348 Louvain-la-Neuve, Belgium
| | - David Lauvergnat
- Institut de Chimie Physique, UMR-CNRS 8000, Université Paris-Saclay, CNRS, 91405 Orsay, France
| |
Collapse
|
4
|
Bossion D, Ndengué S, Meyer HD, Gatti F, Scribano Y. Theoretical investigation of the H + HD → D + H 2 chemical reaction for astrophysical applications: A state-to-state quasi-classical study. J Chem Phys 2020; 153:081102. [PMID: 32872883 DOI: 10.1063/5.0017697] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We report a large set of state-to-state rate constants for the H + HD reactive collision, using Quasi-Classical Trajectory (QCT) simulations on the accurate H3 global potential energy surface of Mielke et al. [J. Chem. Phys. 116, 4142 (2002)]. High relative collision energies (up to ≈56 000 K) and high rovibrational levels of HD (up to ≈50 000 K), relevant to various non thermal equilibrium astrophysical media, are considered. We have validated the accuracy of our QCT calculations with a new efficient adaptation of the Multi Configuration Time Dependent Hartree (MCTDH) method to compute the reaction probability of a specific reactive channel. Our study has revealed that the high temperature regime favors the production of H2 in its highly rovibrationnally excited states, which can de-excite radiatively (cooling the gas) or collisionally (heating the gas). Those new state-to-state QCT reaction rate constants represent a significant improvement in our understanding of the possible mechanisms leading to the destruction of HD by its collision with a H atom.
Collapse
Affiliation(s)
- Duncan Bossion
- Laboratoire Univers et Particules de Montpellier, Université de Montpellier, UMR-CNRS 5299, 34095 Montpellier, France
| | - Steve Ndengué
- ICTP-East African Institute for Fundamental Research, University of Rwanda, Kigali, Rwanda
| | - Hans-Dieter Meyer
- Physikalisch-Chemisches Institut, Universität Heidelberg, Im Neuenheimer Feld 229, D-69120 Heidelberg, Germany
| | - Fabien Gatti
- Institut de Sciences Moleculaires d'Orsay, UMR-CNRS 8214, Université Paris-Saclay, 91405 Orsay, France
| | - Yohann Scribano
- Laboratoire Univers et Particules de Montpellier, Université de Montpellier, UMR-CNRS 5299, 34095 Montpellier, France
| |
Collapse
|
5
|
González-Lezana T, Bossion D, Scribano Y, Bhowmick S, Suleimanov YV. Dynamics of H + HeH +( v = 0, j = 0) → H 2+ + He: Insight on the Possible Complex-Forming Behavior of the Reaction. J Phys Chem A 2019; 123:10480-10489. [PMID: 31725286 DOI: 10.1021/acs.jpca.9b06122] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The H + HeH+→ He + H2+ reaction has been studied by means of a combination of theoretical approaches: a statistical quantum method (SQM), ring polymer molecular dynamics (RPMD), and the quasiclassical trajectory (QCT) method. Cross sections and rate constants have been calculated in an attempt to investigate the dynamics of the process. The comparison with previous calculations and experimental results reveals that despite the fact that statistical predictions seem to reproduce some of the overall observed features, the analysis at a more detailed state-to-state level shows noticeable deviations from a complex-forming dynamics. We find some differences in cross sections and rate constants obtained in the QCT calculation with a Gaussian binning procedure with respect to previous works in which the standard histogram binning was employed.
Collapse
Affiliation(s)
| | - Duncan Bossion
- Laboratoire Univers et Particules de Montpellier, UMR-CNRS 5299 , Université de Montpellier , 34095 Montpellier Cedex , France
| | - Yohann Scribano
- Laboratoire Univers et Particules de Montpellier, UMR-CNRS 5299 , Université de Montpellier , 34095 Montpellier Cedex , France
| | - Somnath Bhowmick
- Computation-based Science and Technology Research Center , The Cyprus Institute , 20 Konstantinou Kavafi Street , Nicosia 2121 , Cyprus
| | - Yury V Suleimanov
- Computation-based Science and Technology Research Center , The Cyprus Institute , 20 Konstantinou Kavafi Street , Nicosia 2121 , Cyprus
| |
Collapse
|
6
|
Ndengué S, Scribano Y, Gatti F, Dawes R. State-to-state inelastic rotational cross sections in five-atom systems with the multiconfiguration time dependent Hartree method. J Chem Phys 2019; 151:134301. [PMID: 31594314 DOI: 10.1063/1.5119381] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We present a MultiConfiguration Time Dependent Hartree (MCTDH) method as an attractive alternative approach to the usual quantum close-coupling method that approaches some computational limits in the calculation of rotational excitation (and de-excitation) between polyatomic molecules (here collisions between triatomic and diatomic rigid molecules). We have performed a computational investigation of the rotational (de-)excitation of the benchmark rigid rotor H2O-H2 system on a recently developed Potential Energy Surface of the complex using the MCTDH method. We focus here on excitations and de-excitations from the 000, 111, and 110 states of H2O with H2 in its ground rotational state, looking at all the potential transitions in the energy range 1-200 cm-1. This work follows a recently completed study on the H2O-H2 cluster where we characterized its spectroscopy and more generally serves a broader goal to describe inelastic collision processes of high dimensional systems using the MCTDH method. We find that the cross sections obtained from the MCTDH calculations are in excellent agreement with time independent calculations from previous studies but does become challenging for the lower kinetic energy range of the de-excitation process: that is, below approximately 20 cm-1 of collision energy, calculations with a relative modest basis become unreliable. The MCTDH method therefore appears to be a useful complement to standard approaches to study inelastic collision for various collision partners, even at low energy, though performing better for rotational excitation than for de-excitation.
Collapse
Affiliation(s)
- Steve Ndengué
- ICTP-East African Institute for Fundamental Research, University of Rwanda, Kigali, Rwanda
| | - Yohann Scribano
- Laboratoire Univers et Particule de Montpellier, Université de Montpellier, UMR-CNRS 5299, 34095 Montpellier Cedex, France
| | - Fabien Gatti
- Institut de Sciences Moléculaires d'Orsay, UMR 8214, Université Paris-Sud - Université Paris-Saclay, 91405 Orsay, France
| | - Richard Dawes
- Department of Chemistry, Missouri University of Science and Technology, Rolla, Missouri 65409, USA
| |
Collapse
|
7
|
Felker PM, Lauvergnat D, Scribano Y, Benoit DM, Bačić Z. Intramolecular stretching vibrational states and frequency shifts of (H2)2 confined inside the large cage of clathrate hydrate from an eight-dimensional quantum treatment using small basis sets. J Chem Phys 2019; 151:124311. [DOI: 10.1063/1.5124051] [Citation(s) in RCA: 20] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Peter M. Felker
- Department of Chemistry and Biochemistry, University of California, Los Angeles, California 90095-1569, USA
| | - David Lauvergnat
- Laboratoire de Chimie Physique, UMR-CNRS 8000, Université de Paris-Sud, Orsay F-91405, France
| | - Yohann Scribano
- Laboratoire Univers et Particule de Montpellier, Université de Montpellier, UMR-CNRS 5299, 34095 Montpellier Cedex, France
| | - David M. Benoit
- Department of Physics and Mathematics, E. A. Milne Centre for Astrophysics and G. W. Gray Centre for Advanced Materials, The University of Hull, Cottingham Road, Kingston upon Hull HU6 7RX, United Kingdom
| | - Zlatko Bačić
- Department of Chemistry, New York University, New York, New York 10003, USA
- NYU-ECNU Center for Computational Chemistry at NYU Shanghai, 3663 Zhongshan Road North, Shanghai 200062, China
| |
Collapse
|
8
|
Stoecklin T, Denis-Alpizar O, Clergerie A, Halvick P, Faure A, Scribano Y. Rigid-Bender Close-Coupling Treatment of the Inelastic Collisions of H2O with para-H2. J Phys Chem A 2019; 123:5704-5712. [DOI: 10.1021/acs.jpca.9b04052] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Thierry Stoecklin
- UMR5255-CNRS, Université de Bordeaux, 351 cours de la libération, F-33405 Talence, France
| | - Otoniel Denis-Alpizar
- Instituto de Ciencias Químicas Aplicadas, Facultad de Ingeniería, Universidad Autónoma de Chile, El Llano Subercaseaux 2801, San Miguel, Santiago, Chile
| | - Alexandre Clergerie
- UMR5255-CNRS, Université de Bordeaux, 351 cours de la libération, F-33405 Talence, France
| | - Philippe Halvick
- UMR5255-CNRS, Université de Bordeaux, 351 cours de la libération, F-33405 Talence, France
| | - Alexandre Faure
- CNRS, IPAG, Université Grenoble Alpes, F-38000 Grenoble, France
| | - Yohann Scribano
- Laboratoire Univers et Particules de Montpellier, UMR-CNRS 5299, Université de Montpellier, F-34095 Montpellier Cedex, France
| |
Collapse
|
9
|
Lauvergnat D, Felker P, Scribano Y, Benoit DM, Bačić Z. H2, HD, and D2in the small cage of structure II clathrate hydrate: Vibrational frequency shifts from fully coupled quantum six-dimensional calculations of the vibration-translation-rotation eigenstates. J Chem Phys 2019; 150:154303. [DOI: 10.1063/1.5090573] [Citation(s) in RCA: 21] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- David Lauvergnat
- Laboratoire de Chimie Physique, UMR-CNRS 8000, Université de Paris-Sud, Orsay F-91405, France
| | - Peter Felker
- Department of Chemistry and Biochemistry, University of California, Los Angeles, California 90095-1569, USA
| | - Yohann Scribano
- Laboratoire Univers et Particules de Montpellier, Université de Montpellier, UMR-CNRS 5299, 34095 Montpellier Cedex, France
| | - David M. Benoit
- Department of Physics and Mathematics, E.A. Milne Centre for Astrophysics and G. W. Gray Centre for Advanced Materials, The University of Hull, Cottingham Road, Kingston upon Hull HU6 7RX, United Kingdom
| | - Zlatko Bačić
- Department of Chemistry, New York University, New York, New York 10003, USA
- NYU-ECNU Center for Computational Chemistry at NYU Shanghai, 3663 Zhongshan Road North, Shanghai, 200062, China
| |
Collapse
|
10
|
Bossion D, Scribano Y, Parlant G. State-to-state quasi-classical trajectory study of the D + H2 collision for high temperature astrophysical applications. J Chem Phys 2019; 150:084301. [DOI: 10.1063/1.5082158] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Duncan Bossion
- Laboratoire Univers et Particules de Montpellier, Université de Montpellier, UMR-CNRS 5299, 34095 Montpellier Cedex, France
| | - Yohann Scribano
- Laboratoire Univers et Particules de Montpellier, Université de Montpellier, UMR-CNRS 5299, 34095 Montpellier Cedex, France
| | - Gérard Parlant
- Institut Charles Gerhardt, Université de Montpellier, UMR-CNRS 5253, 34095 Montpellier, France
| |
Collapse
|
11
|
Bhowmick S, Bossion D, Scribano Y, Suleimanov YV. The low temperature D + + H 2→ HD + H + reaction rate coefficient: a ring polymer molecular dynamics and quasi-classical trajectory study. Phys Chem Chem Phys 2018; 20:26752-26763. [PMID: 30324962 DOI: 10.1039/c8cp05398g] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The reaction between D+ and H2 plays an important role in astrochemistry at low temperatures and also serves as a prototype for a simple ion-molecule reaction. Its ground X[combining tilde]1A' state has a very small thermodynamic barrier (up to 1.8 × 10-2 eV) and the reaction proceeds through the formation of an intermediate complex lying within the potential well with a depth of at least 0.2 eV, thus representing a challenge for dynamical studies. In the present work, we analyze the title reaction within the temperature range of 20-100 K by means of ring polymer molecular dynamics (RPMD) and quasi-classical trajectory (QCT) methods over the full-dimensional global potential energy surface developed by Aguado et al. [A. Aguado, O. Roncero, C. Tablero, C. Sanz and M. Paniagua, J. Chem. Phys., 2000, 112, 1240]. The computed thermal RPMD and QCT rate coefficients are found to be almost independent of temperature and fall within the range of 1.34-2.01 × 10-9 cm3 s-1. They are also in very good agreement with previous time-independent quantum mechanical and statistical quantum method calculations. Furthermore, we observe that the choice of asymptotic separation distance between the reactants can markedly alter the rate coefficient in the low temperature regime (20-50 K). Therefore it is of utmost importance to correctly assign the value of this parameter for dynamical studies, particularly at very low temperatures of astrochemical importance. We finally conclude that the experimental rate measurements for the title reaction are highly desirable in future.
Collapse
Affiliation(s)
- Somnath Bhowmick
- Computation-based Science and Technology Research Center, The Cyprus Institute, 20 Konstantinou Kavafi Street, Nicosia 2121, Cyprus.
| | - Duncan Bossion
- Laboratoire Univers et Particules de Montpellier, UMR-CNRS 5299, Université de Montpellier, Place Eugène Bataillon, 34095 Montpellier, France
| | - Yohann Scribano
- Laboratoire Univers et Particules de Montpellier, UMR-CNRS 5299, Université de Montpellier, Place Eugène Bataillon, 34095 Montpellier, France
| | - Yury V Suleimanov
- Computation-based Science and Technology Research Center, The Cyprus Institute, 20 Konstantinou Kavafi Street, Nicosia 2121, Cyprus.
| |
Collapse
|
12
|
Scribano Y, Parlant G, Poirier B. Communication: Adiabatic quantum trajectory capture for cold and ultra-cold chemical reactions. J Chem Phys 2018; 149:021101. [PMID: 30007377 DOI: 10.1063/1.5041091] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The Langevin capture model is often used to describe barrierless reactive collisions. At very low temperatures, quantum effects may alter this simple capture image and dramatically affect the reaction probability. In this paper, we use the trajectory-ensemble reformulation of quantum mechanics, as recently proposed by one of the authors (Poirier) to compute adiabatic-channel capture probabilities and cross-sections for the highly exothermic reaction Li + CaH(v = 0, j = 0) → LiH + Ca, at low and ultra-low temperatures. Each captured quantum trajectory takes full account of tunneling and quantum reflection along the radial collision coordinate. Our approach is found to be very fast and accurate, down to extremely low temperatures. Moreover, it provides an intuitive and practical procedure for determining the capture distance (i.e., where the capture probability is evaluated), which would otherwise be arbitrary.
Collapse
Affiliation(s)
- Yohann Scribano
- Laboratoire Univers et Particules de Montpellier, UMR-CNRS 5299, Université de Montpellier, Place Eugène Bataillon, 34095 Montpellier, France
| | - Gérard Parlant
- Institut Charles Gerhardt, UMR-CNRS 5253, Université de Montpellier, Place Eugène Bataillon, 34095 Montpellier, France
| | - Bill Poirier
- Max-Planck-Institut für Physik komplexer Systeme, Nöthnitzer Strasse 38, 01187 Dresden, Germany and Department of Chemistry and Biochemistry, and Department of Physics, Texas Tech University, P.O. Box 41061, Lubbock, Texas 79409-1061, USA
| |
Collapse
|
13
|
Powers A, Scribano Y, Lauvergnat D, Mebe E, Benoit DM, Bačić Z. The effect of the condensed-phase environment on the vibrational frequency shift of a hydrogen molecule inside clathrate hydrates. J Chem Phys 2018; 148:144304. [DOI: 10.1063/1.5024884] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Anna Powers
- Department of Chemistry, New York University, New York, New York 10003, USA
| | - Yohann Scribano
- Laboratoire Univers et Particules de Montpellier, Université de Montpellier, LUPM-UMR CNRS 5299, 34095 Montpellier Cedex, France
| | - David Lauvergnat
- Laboratoire de Chimie Physique UMR CNRS 8000-Université de Paris-Sud, Orsay F-91405, France
| | - Elsy Mebe
- Laboratoire de Chimie Physique UMR CNRS 8000-Université de Paris-Sud, Orsay F-91405, France
| | - David M. Benoit
- E.A. Milne Centre for Astrophysics & G.W. Gray Centre for Advanced Materials, Chemistry, The University of Hull, Cottingham Road, Kingston upon Hull HU6 7RX, United Kingdom
| | - Zlatko Bačić
- Department of Chemistry, New York University, New York, New York 10003, USA
- NYU-ECNU Center for Computational Chemistry at NYU Shanghai, 3663 Zhongshan Road North, Shanghai 200062, China
| |
Collapse
|
14
|
Benoit DM, Lauvergnat D, Scribano Y. Does cage quantum delocalisation influence the translation–rotational bound states of molecular hydrogen in clathrate hydrate? Faraday Discuss 2018; 212:533-546. [DOI: 10.1039/c8fd00087e] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
In this study, we examine the effect of a flexible description of the clathrate hydrate framework on the translation–rotation (TR) eigenstates of guest molecules such as molecular hydrogen.
Collapse
Affiliation(s)
- David M. Benoit
- E.A. Milne Centre for Astrophysics
- G.W. Gray Centre for Advanced Materials, Chemistry
- The University of Hull
- UK
| | - David Lauvergnat
- Laboratoire de Chimie Physique
- CNRS
- Univ. Paris-Sud
- Université Paris-Saclay
- Orsay
| | - Yohann Scribano
- Laboratoire Univers et Particules de Montpellier
- UMR-CNRS 5299
- Université de Montpellier
- 34095 Montpellier Cedex
- France
| |
Collapse
|
15
|
Scribano Y, Faure A. Note: On the inclusion of a diagonal Born-Oppenheimer correction in the reduced dimensional treatment of the H2O–para-H2 complex. J Chem Phys 2017; 146:226102. [DOI: 10.1063/1.4986074] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Yohann Scribano
- Laboratoire Univers et Particules de Montpellier, Université de Montpellier, LUPM - UMR CNRS 5299, 34095 Montpellier Cedex, France
| | - Alexandre Faure
- Université Grenoble Alpes, CNRS, IPAG, F-38000 Grenoble, France
| |
Collapse
|
16
|
Spielfiedel A, Senent ML, Kalugina Y, Scribano Y, Balança C, Lique F, Feautrier N. A new ab initio potential energy surface for the collisional excitation of N2H(+) by H2. J Chem Phys 2015; 143:024301. [PMID: 26178099 DOI: 10.1063/1.4923440] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We compute a new potential energy surface (PES) for the study of the inelastic collisions between N2H(+) and H2 molecules. A preliminary study of the reactivity of N2H(+) with H2 shows that neglecting reactive channels in collisional excitation studies is certainly valid at low temperatures. The four dimensional (4D) N2H(+)-H2 PES is obtained from electronic structure calculations using the coupled cluster with single, double, and perturbative triple excitation level of theory. The atoms are described by the augmented correlation consistent triple zeta basis set. Both molecules were treated as rigid rotors. The potential energy surface exhibits a well depth of ≃2530 cm(-1). Considering this very deep well, it appears that converged scattering calculations that take into account the rotational structure of both N2H(+) and H2 should be very difficult to carry out. To overcome this difficulty, the "adiabatic-hindered-rotor" treatment, which allows para-H2(j = 0) to be treated as if it were spherical, was used in order to reduce the scattering calculations to a 2D problem. The validity of this approach is checked and we find that cross sections and rate coefficients computed from the adiabatic reduced surface are in very good agreement with the full 4D calculations.
Collapse
Affiliation(s)
- Annie Spielfiedel
- LERMA, Observatoire de Paris, Sorbonne Université, UPMC Univ Paris 06, CNRS-UMR 8112, F-92195 Meudon, France
| | - Maria Luisa Senent
- Departamento de Química y Física Teóricas, IEM-CSIC, Serrano 121, Madrid 28006, Spain
| | - Yulia Kalugina
- LOMC-UMR 6294, CNRS-Université du Havre, 25 rue Philippe Lebon, BP 1123, 76063 Le Havre, France
| | - Yohann Scribano
- LUPM-UMR 5299, CNRS-Université de Montpellier, Place Eugene Bataillon, 34095 Montpellier Cedex, France
| | - Christian Balança
- LERMA, Observatoire de Paris, Sorbonne Université, UPMC Univ Paris 06, CNRS-UMR 8112, F-92195 Meudon, France
| | - François Lique
- LERMA, Observatoire de Paris, Sorbonne Université, UPMC Univ Paris 06, CNRS-UMR 8112, F-92195 Meudon, France
| | - Nicole Feautrier
- LERMA, Observatoire de Paris, Sorbonne Université, UPMC Univ Paris 06, CNRS-UMR 8112, F-92195 Meudon, France
| |
Collapse
|
17
|
Loreau J, Liévin J, Scribano Y, van der Avoird A. Potential energy surface and bound states of the NH3-Ar and ND3-Ar complexes. J Chem Phys 2015; 141:224303. [PMID: 25494745 DOI: 10.1063/1.4903047] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022] Open
Abstract
A new, four-dimensional potential energy surface for the interaction of NH3 and ND3 with Ar is computed using the coupled-cluster method with single, double, and perturbative triple excitations and large basis sets. The umbrella motion of the ammonia molecule is explicitly taken into account. The bound states of both NH3-Ar and ND3-Ar are calculated on this potential for total angular momentum values from J = 0 to 10, with the inclusion of Coriolis interactions. The energies and splittings of the rovibrational levels are in excellent agreement with the extensive high-resolution spectroscopic data accumulated over the years in the infrared and microwave regions for both complexes, which demonstrates the quality of the potential energy surface.
Collapse
Affiliation(s)
- J Loreau
- Service de Chimie Quantique et Photophysique, Université Libre de Bruxelles (ULB) CP 160/09, 50 av. F.D. Roosevelt, 1050 Brussels, Belgium
| | - J Liévin
- Service de Chimie Quantique et Photophysique, Université Libre de Bruxelles (ULB) CP 160/09, 50 av. F.D. Roosevelt, 1050 Brussels, Belgium
| | - Y Scribano
- Laboratoire Univers et Particules de Montpellier, Université de Montpellier II, LUPM - UMR CNRS 5299, 34095 Montpellier Cedex, France
| | - A van der Avoird
- Theoretical Chemistry, Institute for Molecules and Materials, Radboud University Nijmegen, Heyendaalseweg 135, 6525 AJ Nijmegen, The Netherlands
| |
Collapse
|
18
|
González-Lezana T, Scribano Y, Honvault P. The D(+) + H2 reaction: differential and integral cross sections at low energy and rate constants at low temperature. J Phys Chem A 2014; 118:6416-24. [PMID: 24802076 DOI: 10.1021/jp501446y] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Abstract
The D(+) + H2 reaction is investigated by means of a time independent quantum mechanical (TIQM) and statistical quantum mechanical (SQM) methods. Differential cross sections and product rotational distributions obtained with these two theoretical approaches for collision energies between 1 meV and 0.1 eV are compared to analyze the dynamics of the process. The agreement observed between the TIQM differential cross sections and the SQM predictions as the energy increases revealed the role played by the complex-forming mechanism. The importance of a good description of the asymptotic regions is also investigated by calculating rate constants for the title reaction at low temperature.
Collapse
|
19
|
González-Lezana T, Honvault P, Scribano Y. Dynamics of the D(+) + H2 → HD + H(+) reaction at the low energy regime by means of a statistical quantum method. J Chem Phys 2014; 139:054301. [PMID: 23927256 DOI: 10.1063/1.4816638] [Citation(s) in RCA: 23] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The D(+) +H2(v = 0, j = 0, 1) → HD+H(+) reaction has been investigated at the low energy regime by means of a statistical quantum mechanical (SQM) method. Reaction probabilities and integral cross sections (ICSs) between a collisional energy of 10(-4) eV and 0.1 eV have been calculated and compared with previously reported results of a time independent quantum mechanical (TIQM) approach. The TIQM results exhibit a dense profile with numerous narrow resonances down to Ec ~ 10(-2) eV and for the case of H2(v = 0, j = 0) a prominent peak is found at ~2.5 × 10(-4) eV. The analysis at the state-to-state level reveals that this feature is originated in those processes which yield the formation of rotationally excited HD(v' = 0, j' > 0). The statistical predictions reproduce reasonably well the overall behaviour of the TIQM ICSs at the larger energy range (Ec ≥ 10(-3) eV). Thermal rate constants are in qualitative agreement for the whole range of temperatures investigated in this work, 10-100 K, although the SQM values remain above the TIQM results for both initial H2 rotational states, j = 0 and 1. The enlargement of the asymptotic region for the statistical approach is crucial for a proper description at low energies. In particular, we find that the SQM method leads to rate coefficients in terms of the energy in perfect agreement with previously reported measurements if the maximum distance at which the calculation is performed increases noticeably with respect to the value employed to reproduce the TIQM results.
Collapse
|
20
|
Honvault P, Scribano Y. Correction to “State-to-State Quantum Mechanical Calculations of Rate Coefficients for the D+ + H2 → HD + H+ Reaction at Low Temperature”. J Phys Chem A 2013. [DOI: 10.1021/jp4113377] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
|
21
|
Honvault P, Scribano Y. State-to-state quantum mechanical calculations of rate coefficients for the D+ + H2 → HD + H+ reaction at low temperature. J Phys Chem A 2013; 117:9778-84. [PMID: 23452294 DOI: 10.1021/jp3124549] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The dynamics of the D(+) + H2 → HD + H(+) reaction on a recent ab initio potential energy surface (Velilla, L.; Lepetit, B.; Aguado, A.; Beswick, J. A.; Paniagua, M. J. Chem. Phys. 2008, 129, 084307) has been investigated by means of a time-independent quantum mechanical approach. Cross-sections and rate coefficients are calculated, respectively, for collision energies below 0.1 eV and temperatures up to 100 K for astrophysical application. An excellent accord is found for collision energy above 5 meV, while a disagreement between theory and experiment is observed below this energy. We show that the rate coefficients reveal a slightly temperature-dependent behavior in the upper part of the temperature range considered here. This is in agreement with the experimental data above 80 K, which give a temperature independent value. However, a significant decrease is found at temperatures below 20 K. This decrease can be related to quantum effects and the decay back to the reactant channel, which are not considered by simple statistical approaches, such as the Langevin model. Our results have been fitted to appropriate analytical expressions in order to be used in astrochemical and cosmological models.
Collapse
Affiliation(s)
- P Honvault
- Laboratoire Interdisciplinaire Carnot de Bourgogne, UMR CNRS 6303, Université de Bourgogne , 21078 Dijon Cedex, and UFR Sciences et Techniques, Université de Franche-Comté , 25030 Besançon cedex, France
| | | |
Collapse
|
22
|
Ouk CM, Zvereva-Loëte N, Scribano Y, Bussery-Honvault B. Transition state theory thermal rate constants and RRKM-based branching ratios for the N((2)D) + CH(4) reaction based on multi-state and multi-reference ab initio calculations of interest for the Titan's chemistry. J Comput Chem 2012; 33:2211-24. [PMID: 22782670 DOI: 10.1002/jcc.23054] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/23/2011] [Revised: 04/26/2012] [Accepted: 05/30/2012] [Indexed: 11/11/2022]
Abstract
Multireference single and double configuration interaction (MRCI) calculations including Davidson (+Q) or Pople (+P) corrections have been conducted in this work for the reactants, products, and extrema of the doublet ground state potential energy surface involved in the N((2)D) + CH(4) reaction. Such highly correlated ab initio calculations are then compared with previous PMP4, CCSD(T), W1, and DFT/B3LYP studies. Large relative differences are observed in particular for the transition state in the entrance channel resolving the disagreement between previous ab initio calculations. We confirm the existence of a small but positive potential barrier (3.86 ± 0.84 kJ mol(-1) (MR-AQCC) and 3.89 kJ mol(-1) (MRCI+P)) in the entrance channel of the title reaction. The correlation is seen to change significantly the energetic position of the two minima and five saddle points of this system together with the dissociation channels but not their relative order. The influence of the electronic correlation into the energetic of the system is clearly demonstrated by the thermal rate constant evaluation and it temperature dependance by means of the transition state theory. Indeed, only MRCI values are able to reproduce the experimental rate constant of the title reaction and its behavior with temperature. Similarly, product branching ratios, evaluated by means of unimolecular RRKM theory, confirm the NH production of Umemoto et al., whereas previous works based on less accurate ab initio calculations failed. We confirm the previous findings that the N((2)D) + CH(4) reaction proceeds via an insertion-dissociation mechanism and that the dominant product channels are CH(2)NH + H and CH(3) + NH.
Collapse
Affiliation(s)
- Chanda-Malis Ouk
- Laboratoire ICB, UMR, CNRS-Université de Bourgogne, Dijon, France
| | | | | | | |
Collapse
|
23
|
Ziemkiewicz MP, Pluetzer C, Nesbitt DJ, Scribano Y, Faure A, van der Avoird A. Overtone vibrational spectroscopy in H2-H2O complexes: A combined high level theoretical ab initio, dynamical and experimental study. J Chem Phys 2012; 137:084301. [DOI: 10.1063/1.4732581] [Citation(s) in RCA: 22] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
24
|
Scribano Y, Faure A, Lauvergnat D. Rotational excitation of H2O by para-H2 from an adiabatically reduced dimensional potential. J Chem Phys 2012; 136:094109. [DOI: 10.1063/1.3690881] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022] Open
|
25
|
Scribano Y, Akin-Ojo O, Faure A. Note: Second virial coefficient of the water-hydrogen complex from an explicitly correlated potential energy surface. J Chem Phys 2011; 135:116101. [PMID: 21950886 DOI: 10.1063/1.3639191] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
26
|
Benoit DM, Madebene B, Ulusoy I, Mancera L, Scribano Y, Chulkov S. Towards a scalable and accurate quantum approach for describing vibrations of molecule-metal interfaces. Beilstein J Nanotechnol 2011; 2:427-447. [PMID: 22003450 PMCID: PMC3190614 DOI: 10.3762/bjnano.2.48] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 04/19/2011] [Accepted: 07/14/2011] [Indexed: 05/31/2023]
Abstract
We present a theoretical framework for the computation of anharmonic vibrational frequencies for large systems, with a particular focus on determining adsorbate frequencies from first principles. We give a detailed account of our local implementation of the vibrational self-consistent field approach and its correlation corrections. We show that our approach is both robust, accurate and can be easily deployed on computational grids in order to provide an efficient computational tool. We also present results on the vibrational spectrum of hydrogen fluoride on pyrene, on the thiophene molecule in the gas phase, and on small neutral gold clusters.
Collapse
Affiliation(s)
- David M Benoit
- Nachwuchsgruppe Theorie – SFB 569, Albert-Einstein-Allee 11, University of Ulm, D-89081 Ulm, Germany
- Department of Chemistry, The University of Hull, Cottingham Road, Kingston upon Hull HU6 7RX, United Kingdom
| | - Bruno Madebene
- L.A.D.I.R Boite 49, Bat F74, Université Pierre et Marie Curie, 4 Place Jussieu, F-75252 Paris, France
| | - Inga Ulusoy
- Technische Universität München, Theoretical Chemistry, Lichtenbergstr. 4, D-85747 Garching, Germany
| | - Luis Mancera
- Theoretical Chemistry, Albert-Einstein-Allee 11, University of Ulm, D-89081 Ulm, Germany
| | - Yohann Scribano
- Laboratoire Interdisciplinaire Carnot de Bourgogne,UMR 5209 CNRS – Université de Bourgogne, 9 av. A. Savary, BP 47870, F-21078 Dijon Cedex, France
| | - Sergey Chulkov
- Nachwuchsgruppe Theorie – SFB 569, Albert-Einstein-Allee 11, University of Ulm, D-89081 Ulm, Germany
| |
Collapse
|
27
|
Ulusoy IS, Scribano Y, Benoit DM, Tschetschetkin A, Maurer N, Koslowski B, Ziemann P. Vibrations of a single adsorbed organic molecule: anharmonicity matters! Phys Chem Chem Phys 2011; 13:612-8. [PMID: 21031209 DOI: 10.1039/c0cp01289k] [Citation(s) in RCA: 21] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Affiliation(s)
- I S Ulusoy
- Nachwuchsgruppe Theorie - SFB 569, Albert-Einstein-Allee 11, Ulm University, D-89081 Ulm, Germany
| | | | | | | | | | | | | |
Collapse
|
28
|
|
29
|
Scribano Y, Faure A, Wiesenfeld L. Communication: Rotational excitation of interstellar heavy water by hydrogen molecules. J Chem Phys 2010; 133:231105. [PMID: 21186851 DOI: 10.1063/1.3507877] [Citation(s) in RCA: 19] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Yohann Scribano
- Laboratoire Interdisciplinaire Carnot de Bourgogne-UMR 5209, CNRS-Université de Bourgogne, 9 Av. Alain Savary, B.P. 47870, F-21078 Dijon Cedex, France.
| | | | | |
Collapse
|
30
|
Yang CH, Sarma G, ter Meulen JJ, Parker DH, McBane GC, Wiesenfeld L, Faure A, Scribano Y, Feautrier N. Communication: Mapping water collisions for interstellar space conditions. J Chem Phys 2010; 133:131103. [DOI: 10.1063/1.3475517] [Citation(s) in RCA: 26] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
|
31
|
Scribano Y, Lauvergnat DM, Benoit DM. Fast vibrational configuration interaction using generalized curvilinear coordinates and self-consistent basis. J Chem Phys 2010; 133:094103. [DOI: 10.1063/1.3476468] [Citation(s) in RCA: 56] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
32
|
Belpassi L, Reca ML, Tarantelli F, Roncaratti LF, Pirani F, Cappelletti D, Faure A, Scribano Y. Charge-Transfer Energy in the Water−Hydrogen Molecular Aggregate Revealed by Molecular-Beam Scattering Experiments, Charge Displacement Analysis, and ab Initio Calculations. J Am Chem Soc 2010; 132:13046-58. [DOI: 10.1021/ja1056642] [Citation(s) in RCA: 77] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Leonardo Belpassi
- Dipartimento di Chimica, Università di Perugia, and CNR−Istituto di Scienze e Tecnologie Molecolari, Via Elce di Sotto 8, 06123 Perugia, Italy, Dipartimento di Ingegneria Civile e Ambientale, Università di Perugia, Via G. Duranti 93, 06125 Perugia, Italy, Laboratoire d’Astrophysique, Observatoire de Grenoble, Université Joseph Fourier, CNRS UMR5571, B.P. 53, 38041 Grenoble Cedex 09, France, and Laboratoire Interdisciplinaire Carnot de Bourgogne, UMR 5209 CNRS−Université de Bourgogne, 9 Avenue Alain
| | - Michael L. Reca
- Dipartimento di Chimica, Università di Perugia, and CNR−Istituto di Scienze e Tecnologie Molecolari, Via Elce di Sotto 8, 06123 Perugia, Italy, Dipartimento di Ingegneria Civile e Ambientale, Università di Perugia, Via G. Duranti 93, 06125 Perugia, Italy, Laboratoire d’Astrophysique, Observatoire de Grenoble, Université Joseph Fourier, CNRS UMR5571, B.P. 53, 38041 Grenoble Cedex 09, France, and Laboratoire Interdisciplinaire Carnot de Bourgogne, UMR 5209 CNRS−Université de Bourgogne, 9 Avenue Alain
| | - Francesco Tarantelli
- Dipartimento di Chimica, Università di Perugia, and CNR−Istituto di Scienze e Tecnologie Molecolari, Via Elce di Sotto 8, 06123 Perugia, Italy, Dipartimento di Ingegneria Civile e Ambientale, Università di Perugia, Via G. Duranti 93, 06125 Perugia, Italy, Laboratoire d’Astrophysique, Observatoire de Grenoble, Université Joseph Fourier, CNRS UMR5571, B.P. 53, 38041 Grenoble Cedex 09, France, and Laboratoire Interdisciplinaire Carnot de Bourgogne, UMR 5209 CNRS−Université de Bourgogne, 9 Avenue Alain
| | - Luiz F. Roncaratti
- Dipartimento di Chimica, Università di Perugia, and CNR−Istituto di Scienze e Tecnologie Molecolari, Via Elce di Sotto 8, 06123 Perugia, Italy, Dipartimento di Ingegneria Civile e Ambientale, Università di Perugia, Via G. Duranti 93, 06125 Perugia, Italy, Laboratoire d’Astrophysique, Observatoire de Grenoble, Université Joseph Fourier, CNRS UMR5571, B.P. 53, 38041 Grenoble Cedex 09, France, and Laboratoire Interdisciplinaire Carnot de Bourgogne, UMR 5209 CNRS−Université de Bourgogne, 9 Avenue Alain
| | - Fernando Pirani
- Dipartimento di Chimica, Università di Perugia, and CNR−Istituto di Scienze e Tecnologie Molecolari, Via Elce di Sotto 8, 06123 Perugia, Italy, Dipartimento di Ingegneria Civile e Ambientale, Università di Perugia, Via G. Duranti 93, 06125 Perugia, Italy, Laboratoire d’Astrophysique, Observatoire de Grenoble, Université Joseph Fourier, CNRS UMR5571, B.P. 53, 38041 Grenoble Cedex 09, France, and Laboratoire Interdisciplinaire Carnot de Bourgogne, UMR 5209 CNRS−Université de Bourgogne, 9 Avenue Alain
| | - David Cappelletti
- Dipartimento di Chimica, Università di Perugia, and CNR−Istituto di Scienze e Tecnologie Molecolari, Via Elce di Sotto 8, 06123 Perugia, Italy, Dipartimento di Ingegneria Civile e Ambientale, Università di Perugia, Via G. Duranti 93, 06125 Perugia, Italy, Laboratoire d’Astrophysique, Observatoire de Grenoble, Université Joseph Fourier, CNRS UMR5571, B.P. 53, 38041 Grenoble Cedex 09, France, and Laboratoire Interdisciplinaire Carnot de Bourgogne, UMR 5209 CNRS−Université de Bourgogne, 9 Avenue Alain
| | - Alexandre Faure
- Dipartimento di Chimica, Università di Perugia, and CNR−Istituto di Scienze e Tecnologie Molecolari, Via Elce di Sotto 8, 06123 Perugia, Italy, Dipartimento di Ingegneria Civile e Ambientale, Università di Perugia, Via G. Duranti 93, 06125 Perugia, Italy, Laboratoire d’Astrophysique, Observatoire de Grenoble, Université Joseph Fourier, CNRS UMR5571, B.P. 53, 38041 Grenoble Cedex 09, France, and Laboratoire Interdisciplinaire Carnot de Bourgogne, UMR 5209 CNRS−Université de Bourgogne, 9 Avenue Alain
| | - Yohann Scribano
- Dipartimento di Chimica, Università di Perugia, and CNR−Istituto di Scienze e Tecnologie Molecolari, Via Elce di Sotto 8, 06123 Perugia, Italy, Dipartimento di Ingegneria Civile e Ambientale, Università di Perugia, Via G. Duranti 93, 06125 Perugia, Italy, Laboratoire d’Astrophysique, Observatoire de Grenoble, Université Joseph Fourier, CNRS UMR5571, B.P. 53, 38041 Grenoble Cedex 09, France, and Laboratoire Interdisciplinaire Carnot de Bourgogne, UMR 5209 CNRS−Université de Bourgogne, 9 Avenue Alain
| |
Collapse
|
33
|
Scribano Y, Benoit DM. Calculation of vibrational frequencies through a variational reduced-coupling approach. J Chem Phys 2007; 127:164118. [DOI: 10.1063/1.2798104] [Citation(s) in RCA: 36] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
34
|
Scribano Y, Leforestier C. Contribution of water dimer absorption to the millimeter and far infrared atmospheric water continuum. J Chem Phys 2007; 126:234301. [PMID: 17600414 DOI: 10.1063/1.2746038] [Citation(s) in RCA: 58] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022] Open
Abstract
We present a rigorous calculation of the contribution of water dimers to the absorption coefficient alpha(nu,T) in the millimeter and far infrared domains, over a wide range (276-310 K) of temperatures. This calculation relies on the explicit consideration of all possible transitions within the entire rovibrational bound state manifold of the dimer. The water dimer is described by the flexible 12-dimensional potential energy surface previously fitted to far IR transitions [C. Leforestier et al., J. Chem. Phys. 117, 8710 (2002)], and which was recently further validated by the good agreement obtained for the calculated equilibrium constant Kp(T) with experimental data [Y. Scribano et al., J. Phys. Chem. A. 110, 5411 (2006)]. Transition dipole matrix elements were computed between all rovibrational states up to an excitation energy of 750 cm(-1), and J=K=5 rotational quantum numbers. It was shown by explicit calculations that these matrix elements could be extrapolated to much higher J values (J=30). Transitions to vibrational states located higher in energy were obtained from interpolation of computed matrix elements between a set of initial states spanning the 0-750 cm(-1) range and all vibrational states up to the dissociation limit (approximately 1200 cm(-1)). We compare our calculations with available experimental measurements of the water continuum absorption in the considered range. It appears that water dimers account for an important fraction of the observed continuum absorption in the millimeter region (0-10 cm(-1)). As frequency increases, their relative contribution decreases, becoming small (approximately 3%) at the highest frequency considered nu=944 cm(-1).
Collapse
Affiliation(s)
- Yohann Scribano
- Institut Charles Gerhardt (CTMM), CC 1501, Université Montpellier II-CNRS, 34095 Montpellier, Cedex 05, France
| | | |
Collapse
|
35
|
Harker HA, Keutsch FN, Leforestier C, Scribano Y, Han JX, Saykally RJ. Refinements in the description of excited VRT states of the water dimer. Mol Phys 2007. [DOI: 10.1080/00268970601153373] [Citation(s) in RCA: 12] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/23/2022]
|
36
|
Harker HA, Keutsch FN, Leforestier C, Scribano Y, Han JX, Saykally RJ. Toward a precise determination of the acceptor switching splitting in the water dimer. Mol Phys 2007. [DOI: 10.1080/00268970601063812] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/23/2022]
|
37
|
Scribano Y, Goldman N, Saykally RJ, Leforestier C. Water Dimers in the Atmosphere III: Equilibrium Constant from a Flexible Potential. J Phys Chem A 2006; 110:5411-9. [PMID: 16623469 DOI: 10.1021/jp056759k] [Citation(s) in RCA: 115] [Impact Index Per Article: 6.4] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/03/2023]
Abstract
We present new results for the water dimer equilibrium constant K(p)(T) in the range 190-390 K, using a flexible potential energy surface fitted to spectroscopical data. The increased numerical complexity due to explicit consideration of the monomer vibrations is handled via an adiabatic (6 + 6)d decoupling between intra- and intermolecular modes. The convergence of the canonical partition function of the dimer is ensured by computing all energy levels up to dissociation for total angular momentum values J = 0-5 and using an extrapolation scheme to higher values. The newly calculated values for K(p)(T) are in very good agreement with available experimental data at room temperature. At higher temperatures, an analysis of the convergence of the partition function reveals that quasi-bound states are likely to contribute to the equilibrium constant. Additional thermodynamical quantities (deltaG, deltaH, deltaS, and C(p)) have also been determined and fit to quadratic expressions a + bT + cT2.
Collapse
Affiliation(s)
- Yohann Scribano
- UMR-CNRS 5636, CC 014, Université Montpellier II, 34095 Montpellier, France
| | | | | | | |
Collapse
|