1
|
Feigl E, Jedlovszky P, Sega M. Percolation transition and bimodal density distribution in hydrogen fluoride. J Chem Phys 2024; 160:204503. [PMID: 38785286 DOI: 10.1063/5.0207202] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/06/2024] [Accepted: 05/06/2024] [Indexed: 05/25/2024] Open
Abstract
Hydrogen-bond networks in associating fluids can be extremely robust and characterize the topological properties of the liquid phase, as in the case of water, over its whole domain of stability and beyond. Here, we report on molecular dynamics simulations of hydrogen fluoride (HF), one of the strongest hydrogen-bonding molecules. HF has more limited connectivity than water but can still create long, dynamic chains, setting it apart from most other small molecular liquids. Our simulation results provide robust evidence of a second-order percolation transition of HF's hydrogen bond network occurring below the critical point. This behavior is remarkable as it underlines the presence of two different cohesive mechanisms in liquid HF, one at low temperatures characterized by a spanning network of long, entangled hydrogen-bonded polymers, as opposed to short oligomers bound by the dispersion interaction above the percolation threshold. This second-order phase transition underlines the presence of marked structural heterogeneity in the fluid, which we found in the form of two liquid populations with distinct local densities.
Collapse
Affiliation(s)
- Elija Feigl
- Faculty of Physics, University of Vienna, Boltzmanngasse 5, Wien A-1090, Austria
| | - Pál Jedlovszky
- Department of Chemistry, Eszterházy Károly Catholic University, Leányka utca 12, H-3300 Eger, Hungary
| | - Marcello Sega
- Department of Chemical Engineering, University College London, WC1E 7JE London, United Kingdom
| |
Collapse
|
2
|
Ovsyannikov RI, Makhnev VY, Zobov NF, Koput J, Tennyson J, Polyansky OL. Highly accurate HF dimer ab initio potential energy surface. J Chem Phys 2022; 156:164305. [PMID: 35490002 DOI: 10.1063/5.0083563] [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
A highly accurate, (HF)2 potential energy surface (PES) is constructed based on ab initio calculations performed at the coupled-cluster single double triple level of theory with an aug-cc-pVQZ-F12 basis set at about 152 000 points. A higher correlation correction is computed at coupled-cluster single double triple quadruple level for 2000 points and is considered alongside other more minor corrections due to relativity, core-valence correlation, and Born-Oppenheimer failure. The analytical surface constructed uses 500 constants to reproduce the ab initio points with a standard deviation of 0.3 cm-1. Vibration-rotation-inversion energy levels of the HF dimer are computed for this PES by variational solution of the nuclear-motion Schrödinger equation using the program WAVR4. Calculations over an extended range of rotationally excited states show very good agreement with the experimental data. In particular, the known empirical rotational constants B for the ground vibrational states are predicted to better than about 2 MHz. B constants for excited vibrational states are reproduced several times more accurately than by previous calculations. This level of accuracy is shown to extend to higher excited inter-molecular vibrational states v and higher excited rotational quantum numbers (J, Ka).
Collapse
Affiliation(s)
- Roman I Ovsyannikov
- Institute of Applied Physics, Russian Academy of Sciences, 46 Ulyanov Street, Nizhny Novgorod 603950, Russia
| | - Vladimir Yu Makhnev
- Institute of Applied Physics, Russian Academy of Sciences, 46 Ulyanov Street, Nizhny Novgorod 603950, Russia
| | - Nikolai F Zobov
- Institute of Applied Physics, Russian Academy of Sciences, 46 Ulyanov Street, Nizhny Novgorod 603950, Russia
| | - Jacek Koput
- Faculty of Chemistry, Adam Mickiewicz University, Poznan, Poland
| | - Jonathan Tennyson
- Department of Physics and Astronomy, University College London, Gower Street, London WC1E 6BT, United Kingdom
| | - Oleg L Polyansky
- Institute of Applied Physics, Russian Academy of Sciences, 46 Ulyanov Street, Nizhny Novgorod 603950, Russia
| |
Collapse
|
3
|
Saito K, Torii H. Hidden Halogen-Bonding Ability of Fluorine Manifesting in the Hydrogen-Bond Configurations of Hydrogen Fluoride. J Phys Chem B 2021; 125:11742-11750. [PMID: 34662140 DOI: 10.1021/acs.jpcb.1c07211] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Abstract
Elucidating how the intermolecular interactions of a covalently bonded fluorine atom are similar to and different from those of the other halogen atoms will be helpful for a better unified understanding of them. In the present study, the case of hydrogen fluoride is theoretically studied from this viewpoint by using the techniques of electron density analysis, molecular dynamics of liquid, and others. It is shown that the extra-point model, which locates an additional charge site on the line extended from (not within) the covalent bond and has been adopted for halogen-bonding systems as a key to the generation of proper stability and directionality, works well also in this case. A significantly bent hydrogen-bond configuration, which is characteristic of the intermolecular interactions of hydrogen fluoride, is reasonably well reproduced, meaning that it is a manifestation of the latent halogen-bonding ability, which is hidden by the strongly electronegative nature.
Collapse
|
4
|
Anderson DT, Fajardo ME, Lindsay CM. High resolution infrared spectroscopy of (HCl) 2 and (DCl) 2 isolated in solid parahydrogen: Interchange-tunneling in a quantum solid. J Chem Phys 2021; 154:164309. [PMID: 33940830 DOI: 10.1063/5.0049599] [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
Infrared spectroscopic studies of weakly bound clusters isolated in solid parahydrogen (pH2) that exhibit large-amplitude tunneling motions are needed to probe how quantum solvation perturbs these types of coherent dynamics. We report high resolution Fourier transform infrared absorption spectra of (HCl)2, HCl-DCl, and (DCl)2 isolated in solid pH2 in the 2.4-4.8 K temperature range. The (HCl)2 spectra show a remarkable amount of fine structures that can be rigorously assigned to vibration-rotation-tunneling transitions of (HCl)2 trapped in double substitution sites in the pH2 matrix where end-over-end rotation of the cluster is quenched. The spectra are assigned using a combination of isotopically (H/D and 35Cl/37Cl) enriched samples, polarized IR absorption measurements, and four-line combination differences. The interchange-tunneling (IT) splitting in the ground vibrational state for in-plane and out-of-plane H35Cl-H37Cl dimers is 6.026(1) and 6.950(1) cm-1, respectively, which are factors of 2.565 and 2.224 smaller than in the gas phase dimer. In contrast, the (DCl)2 results show larger perturbations where the ground vibrational state IT splitting in D35Cl-D37Cl is 1.141(1) cm-1, which is a factor of 5.223 smaller than in the gas phase, and the tunneling motion is quenched in excited intramolecular vibrational states. The results are compared to similar measurements on (HCl)2 made in liquid helium nanodroplets to illustrate the similarities and differences in how both these quantum solvents interact with large amplitude tunneling motions of an embedded chromophore.
Collapse
Affiliation(s)
- David T Anderson
- Department of Chemistry, University of Wyoming, Laramie, Wyoming 82071, USA
| | - Mario E Fajardo
- Air Force Research Laboratory, Munitions Directorate, 2306 Perimeter Rd., Eglin AFB, Florida 32542-5910, USA
| | - C Michael Lindsay
- Air Force Research Laboratory, Munitions Directorate, 2306 Perimeter Rd., Eglin AFB, Florida 32542-5910, USA
| |
Collapse
|
5
|
Orabi EA, Faraldo-Gómez JD. New Molecular-Mechanics Model for Simulations of Hydrogen Fluoride in Chemistry and Biology. J Chem Theory Comput 2020; 16:5105-5126. [PMID: 32615034 DOI: 10.1021/acs.jctc.0c00247] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Abstract
Hydrogen fluoride (HF) is the most polar diatomic molecule and one of the simplest molecules capable of hydrogen-bonding. HF deviates from ideality both in the gas phase and in solution and is thus of great interest from a fundamental standpoint. Pure and aqueous HF solutions are broadly used in chemical and industrial processes, despite their high toxicity. HF is a stable species also in some biological conditions, because it does not readily dissociate in water unlike other hydrogen halides; yet, little is known about how HF interacts with biomolecules. Here, we set out to develop a molecular-mechanics model to enable computer simulations of HF in chemical and biological applications. This model is based on a comprehensive high-level ab initio quantum chemical investigation of the structure and energetics of the HF monomer and dimer; (HF)n clusters, for n = 3-7; various clusters of HF and H2O; and complexes of HF with analogs of all 20 amino acids and of several commonly occurring lipids, both neutral and ionized. This systematic analysis explains the unique properties of this molecule: for example, that interacting HF molecules favor nonlinear geometries despite being diatomic and that HF is a strong H-bond donor but a poor acceptor. The ab initio data also enables us to calibrate a three-site molecular-mechanics model, with which we investigate the structure and thermodynamic properties of gaseous, liquid, and supercritical HF in a wide range of temperatures and pressures; the solvation structure of HF in water and of H2O in liquid HF; and the free diffusion of HF across a lipid bilayer, a key process underlying the high cytotoxicity of HF. Despite its inherent simplifications, the model presented significantly improves upon previous efforts to capture the properties of pure and aqueous HF fluids by molecular-mechanics methods and to our knowledge constitutes the first parameter set calibrated for biomolecular simulations.
Collapse
Affiliation(s)
- Esam A Orabi
- Theoretical Molecular Biophysics Laboratory, National Heart, Lung and Blood Institute, National Institutes of Health, 10 Center Drive, Bethesda, Maryland 20814, United States
| | - José D Faraldo-Gómez
- Theoretical Molecular Biophysics Laboratory, National Heart, Lung and Blood Institute, National Institutes of Health, 10 Center Drive, Bethesda, Maryland 20814, United States
| |
Collapse
|
6
|
Remya PR, Mishra BK, Ramachandran C, Sathyamurthy N. Effect of confinement on structure, energy and vibrational spectra of (HF) , n = 1–4. Chem Phys Lett 2019. [DOI: 10.1016/j.cplett.2019.136670] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/26/2022]
|
7
|
Huang J, Yang D, Zhou Y, Xie D. A new full-dimensional ab initio intermolecular potential energy surface and vibrational states for (HF)2 and (DF)2. J Chem Phys 2019; 150:154302. [DOI: 10.1063/1.5090225] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Jing Huang
- Institute of Theoretical and Computational Chemistry, Key Laboratory of Mesoscopic Chemistry, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, China
| | - Dongzheng Yang
- Institute of Theoretical and Computational Chemistry, Key Laboratory of Mesoscopic Chemistry, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, China
| | - Yanzi Zhou
- Institute of Theoretical and Computational Chemistry, Key Laboratory of Mesoscopic Chemistry, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, China
| | - Daiqian Xie
- Institute of Theoretical and Computational Chemistry, Key Laboratory of Mesoscopic Chemistry, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, China
| |
Collapse
|
8
|
Ab initio studies of the structures and vibrational spectra of the hydrogen halide and lithium halide homo- and heterodimers and some mixed hydrogen halide/lithium halide heterodimers. Theor Chem Acc 2016. [DOI: 10.1007/s00214-016-1949-y] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/21/2022]
|
9
|
Samanta AK, Wang Y, Mancini JS, Bowman JM, Reisler H. Energetics and Predissociation Dynamics of Small Water, HCl, and Mixed HCl–Water Clusters. Chem Rev 2016; 116:4913-36. [DOI: 10.1021/acs.chemrev.5b00506] [Citation(s) in RCA: 41] [Impact Index Per Article: 5.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Amit K. Samanta
- Department
of Chemistry, University of Southern California, Los Angeles, California 90089-0482, United States
| | - Yimin Wang
- Department
of Chemistry and Cherry L. Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, United States
| | - John S. Mancini
- Department
of Chemistry and Cherry L. Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, United States
| | - Joel M. Bowman
- Department
of Chemistry and Cherry L. Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, United States
| | - Hanna Reisler
- Department
of Chemistry, University of Southern California, Los Angeles, California 90089-0482, United States
| |
Collapse
|
10
|
Jelil M, Abaydulla A. Graph theoretical enumeration of topology-distinct structures for hydrogen fluoride clusters (HF)n (n ≤ 6). J Chem Phys 2015; 143:044301. [PMID: 26233123 DOI: 10.1063/1.4926939] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
A graph theoretical procedure to generate all the possible topology-distinct structures for hydrogen fluoride (HF) clusters is presented in this work. The hydrogen bond matrix is defined and used to enumerate the topology-distinct structures of hydrogen fluoride (HF)n (n = 2-8) clusters. From close investigation of the structural patterns obtained, several restrictions that should be satisfied for a structure of the HF clusters to be stable are found. The corresponding digraphs of generated hydrogen bond matrices are used as the theoretical framework to obtain all the topology-distinct local minima for (HF)n (n ≤ 6), at the level of MP2/6-31G**(d, p) of ab initio MO method and B3LYP/6-31G**(d, p) of density functional theory method. For HF clusters up to tetramers, the local minimum structures that we generated are same as those in the literature. For HF pentamers and hexamers, we found some new local minima structures which had not been obtained previously.
Collapse
Affiliation(s)
- Mahmutjan Jelil
- Department of Chemistry and Environmental Science, Kashgar University, 029, Xueyuan Road, Kashgar, Xinjiang 844008, China
| | - Alimjan Abaydulla
- Department of Chemistry and Environmental Science, Kashgar University, 029, Xueyuan Road, Kashgar, Xinjiang 844008, China
| |
Collapse
|
11
|
Alkorta I, Elguero J, Del Bene JE. Exploring the PX3:NCH and PX3:NH3 potential surfaces, with X = F, Cl, and Br. Chem Phys Lett 2015. [DOI: 10.1016/j.cplett.2015.10.050] [Citation(s) in RCA: 21] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/22/2022]
|
12
|
Mazack MJM, Gao J. Quantum mechanical force field for hydrogen fluoride with explicit electronic polarization. J Chem Phys 2015; 140:204501. [PMID: 24880295 DOI: 10.1063/1.4875922] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The explicit polarization (X-Pol) theory is a fragment-based quantum chemical method that explicitly models the internal electronic polarization and intermolecular interactions of a chemical system. X-Pol theory provides a framework to construct a quantum mechanical force field, which we have extended to liquid hydrogen fluoride (HF) in this work. The parameterization, called XPHF, is built upon the same formalism introduced for the XP3P model of liquid water, which is based on the polarized molecular orbital (PMO) semiempirical quantum chemistry method and the dipole-preserving polarization consistent point charge model. We introduce a fluorine parameter set for PMO, and find good agreement for various gas-phase results of small HF clusters compared to experiments and ab initio calculations at the M06-2X/MG3S level of theory. In addition, the XPHF model shows reasonable agreement with experiments for a variety of structural and thermodynamic properties in the liquid state, including radial distribution functions, interaction energies, diffusion coefficients, and densities at various state points.
Collapse
Affiliation(s)
- Michael J M Mazack
- Department of Chemistry and Supercomputing Institute, University of Minnesota, 207 Pleasant Street, SE, Minneapolis, Minnesota 55455-0431, USA
| | - Jiali Gao
- Department of Chemistry and Supercomputing Institute, University of Minnesota, 207 Pleasant Street, SE, Minneapolis, Minnesota 55455-0431, USA
| |
Collapse
|
13
|
Silva AF, Richter WE, Terrabuio LA, Haiduke RLA, Bruns RE. Quantum theory of atoms in molecules/charge-charge flux-dipole flux models for fundamental vibrational intensity changes on H-bond formation of water and hydrogen fluoride. J Chem Phys 2014; 140:084306. [DOI: 10.1063/1.4865938] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
14
|
Asselin P, Soulard P, Madebène B, Goubet M, Huet TR, Georges R, Pirali O, Roy P. The cyclic ground state structure of the HF trimer revealed by far infrared jet-cooled Fourier transform spectroscopy. Phys Chem Chem Phys 2014; 16:4797-806. [DOI: 10.1039/c3cp55047h] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
|
15
|
Terrabuio LA, Richter WE, Silva AF, Bruns RE, Haiduke RLA. An atom in molecules study of infrared intensity enhancements in fundamental donor stretching bands in hydrogen bond formation. Phys Chem Chem Phys 2014; 16:24920-8. [DOI: 10.1039/c4cp03489a] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
A semi-quantitative explanation for infrared intensity enhancements in hydrogen bonding is provided by a charge–charge flux interaction contribution.
Collapse
Affiliation(s)
- Luiz A. Terrabuio
- Departamento de Química e Física Molecular
- Instituto de Química de São Carlos
- Universidade de São Paulo
- São Carlos, Brazil
| | - Wagner E. Richter
- Instituto de Química
- Universidade Estadual de Campinas
- Campinas, Brazil
| | - Arnaldo F. Silva
- Instituto de Química
- Universidade Estadual de Campinas
- Campinas, Brazil
| | - Roy E. Bruns
- Instituto de Química
- Universidade Estadual de Campinas
- Campinas, Brazil
| | - Roberto L. A. Haiduke
- Departamento de Química e Física Molecular
- Instituto de Química de São Carlos
- Universidade de São Paulo
- São Carlos, Brazil
| |
Collapse
|
16
|
The Interplay of Nuclear and Electron Wavepacket Motion in the Control of Molecular Processes: A Theoretical Perspective. ACTA ACUST UNITED AC 2013. [DOI: 10.1007/978-3-642-45290-1_8] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/21/2023]
|
17
|
Suhm MA, Kollipost F. Femtisecond single-mole infrared spectroscopy of molecular clusters. Phys Chem Chem Phys 2013; 15:10702-21. [DOI: 10.1039/c3cp51515j] [Citation(s) in RCA: 69] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
18
|
Ji-Min J. M. Yan Y, Williams DE. Quantum-chemical investigation on hydrogen bonding interaction of hydrogen fluoride dimer at various mutual orientations. CHINESE J CHEM 2011. [DOI: 10.1002/cjoc.19900080103] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
|
19
|
Friedrich J, Perlt E, Roatsch M, Spickermann C, Kirchner B. Coupled Cluster in Condensed Phase. Part I: Static Quantum Chemical Calculations of Hydrogen Fluoride Clusters. J Chem Theory Comput 2011; 7:843-51. [DOI: 10.1021/ct100131c] [Citation(s) in RCA: 36] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Joachim Friedrich
- Institute for Chemistry, Chemnitz University of Technology, Strasse der Nationen 62, 09111 Chemnitz, Germany
| | - Eva Perlt
- Wilhelm-Ostwald-Institut für Physikalische und Theoretische Chemie, Universität Leipzig, Linnéstrasse 2, D-04103 Leipzig, Germany
| | - Martin Roatsch
- Wilhelm-Ostwald-Institut für Physikalische und Theoretische Chemie, Universität Leipzig, Linnéstrasse 2, D-04103 Leipzig, Germany
| | - Christian Spickermann
- Wilhelm-Ostwald-Institut für Physikalische und Theoretische Chemie, Universität Leipzig, Linnéstrasse 2, D-04103 Leipzig, Germany
| | - Barbara Kirchner
- Wilhelm-Ostwald-Institut für Physikalische und Theoretische Chemie, Universität Leipzig, Linnéstrasse 2, D-04103 Leipzig, Germany
| |
Collapse
|
20
|
VALLAURI PALJEDLOVSZKYRENZO. Structural properties of liquid HF: a computer simulation investigation. Mol Phys 2010. [DOI: 10.1080/002689798169393] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/16/2022]
|
21
|
Abstract
Spectroscopic characterizations of the stereochemistry of complexes of ammonia (NH(3)) have strongly confirmed some long-held ideas about the weak interactions of NH(3) while casting doubt on others. As expected, NH(3) is observed to be a nearly universal proton acceptor, accepting hydrogen bonds from even some of the weakest proton donors. Surprisingly, no evidence has been found to support the view that NH(3) acts as a proton donor through hydrogen bonding. A critical evaluation of the work that has been done to gather such evidence, as well as of earlier work involving condensed-phase observations, suggests that NH(3) might well be best described as a powerful hydrogen-bond acceptor with little propensity to donate hydrogen bonds.
Collapse
|
22
|
McGrath MJ, Ghogomu JN, Mundy CJ, Kuo IFW, Siepmann JI. First principles Monte Carlo simulations of aggregation in the vapor phase of hydrogen fluoride. Phys Chem Chem Phys 2010; 12:7678-87. [DOI: 10.1039/b924506e] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
23
|
He Y, Müller HB, Quack M, Suhm MA. High Resolution FTIR and Diode Laser Supersonic Jet Spectroscopy of the N = 2 HF Stretching Polyad in (HF)2 and (HFDF): Hydrogen Bond Switching and Predissociation Dynamics. ACTA ACUST UNITED AC 2009. [DOI: 10.1524/zpch.2007.221.11-12.1581] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/24/2022]
Abstract
We report Fourier transform infrared (FTIR) and high resolution diode laser spectra (∼ 1MHz instrumental bandwidth) obtained in cooled absorption cells as well as in a supersonic jet expansion for the N = 2 polyad region of the HF-stretching vibrations of (HF)2, HFDF and DFHF. Three vibrational transitions have been observed for (HF)2 and two for both monodeuterated isotopomers. For (HF)2 we have identified and analysed the observed transitions of the polyad member 22 of the type Δ K
a = 0 and Δ K
a = ± 1 up to rotational sublevel Δ K
a = 3. Band centers as well as rotational constants of all four K
a states have been determined. The tunneling splittings due to hydrogen bond switching for these four K
a states have been investigated, with the Δ K
a = 0 up to Δ K
a = 2 sublevels having tunneling symmetry Γ
vt = A
+ for the lower tunneling states, and switching periods ranging from 158ps for K
a = 0 to 1.35ns for K
a = 2. A tunneling level inversion is found at Δ K
a = 3, leading to a symmetry Γ
vt = B
+ for the lower tunneling state of this K
a-sublevel. The vibrational assignment of the measured spectra of (HF)2 was established by comparison with the monodeuterated isotopomers HFDF and DFHF. For HFDF we have identified and analysed five subbands between 7600cm-1 and 7730cm-1. We have determined the spectroscopic constants of the rotational levels Δ K
a = 0 and Δ K
a = 1 for the vibrationally excited state and of the levels of Δ K
a = 1 and Δ K
a = 2 of the ground state, the latter from combination differences. From the measurements in a supersonic jet expansion we determined the predissociation line width of the N = 22, K
a = 1 to be about 120MHz for the Γ
vt = A
+ tunneling state of (HF)2 and about 90MHz for Γ
vt = B
+. For the Δ K
a = 0 level of N = 22 we obtained predissociation line widths ranging around 100MHz, similar to those of the Δ K
a = 1 level. In the case of HFDF, the predissociation line width of Δ K
a = 1 is about 80MHz. Predissociation lifetimes for these levels with the unbonded HF stretching excited thus are in the range of about 1 to 2ns. The predissociation width in the N = 21 level is uncertain by about a factor three with lg(Δν/MHz) = (3 ± 0.5) and in N = 23 it is about 600MHz corresponding to rounded lifetimes of 0.1ns and 0.3ns when the bonded HF stretching is excited thereby demonstrating strongly mode selective predissociation rates in the N = 2 polyad. Under thermal equilibrium conditions we derived the pressure broadening coefficient for (HF)2 (γ = (6 ± 1) × 10-4cm-1/mbar in the wavenumber range between 7713cm-1 and 7721cm-1 for total gas pressures between 10 and 60mbar, all values as full widths half maximum). For absolute frequency calibrations we have remeasured the first overtone transitions of the monomer HF with much improved precision between P(5) (7515.80151cm-1) and R(7) (7966.22188cm-1).
Collapse
|
24
|
Zhang XP, Lee WB, Zhao DF, Hsiao MK, Chen YL, Lin KC. Photodissociation of (ICN)2 van der Waals dimer using velocity imaging technique. J Chem Phys 2009; 130:214305. [DOI: 10.1063/1.3148376] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
25
|
Thar J, Kirchner B. Relativistic all-electron molecular dynamics simulations. J Chem Phys 2009; 130:124103. [DOI: 10.1063/1.3086655] [Citation(s) in RCA: 28] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
26
|
Lii JH, Allinger NL. The important role of lone-pairs in force field (MM4) calculations on hydrogen bonding in alcohols. J Phys Chem A 2008; 112:11903-13. [PMID: 18942820 DOI: 10.1021/jp804581h] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
An expanded treatment of hydrogen bonding has been developed for MM4 force field calculations, which is an extension from the traditional van der Waals-electrostatic model. It adds explicit hydrogen-bond angularity by the inclusion of lone-pair directionality. The vectors that account for this directionality are placed along the hydrogen acceptor and its chemically intuitive electron pairs. No physical lone-pairs are used in the calculations. Instead, an H-bond angularity function, and a lone-pair directionality function, are incorporated into the hydrogen-bond term. The inclusion of the lone-pair directionality results in improved accuracy in hydrogen-bonded geometries and interaction energies. In this work is described hydrogen bonding in alcohols, and also in water and hydrogen fluoride dimer. The extension to other compounds such as aldehydes, ketones, amides, and so on is straightforward and will be discussed in future work. The conformational energies of ethylene glycol are discussed.
Collapse
Affiliation(s)
- Jenn-Huei Lii
- Department of Chemistry, University of Georgia, Athens, Georgia 30602-2526, USA
| | | |
Collapse
|
27
|
Avilés MW, McCandless ML, Curotto E. Stereographic projection path integral simulations of (HCl)n clusters (n=2–5): Evidence of quantum induced melting in small hydrogen bonded networks. J Chem Phys 2008; 128:124517. [DOI: 10.1063/1.2837802] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
28
|
Tsuchida N, Yamazaki S, Yamabe S. A new mechanism for the Favorskii rearrangement. Org Biomol Chem 2008; 6:3109-17. [DOI: 10.1039/b806577b] [Citation(s) in RCA: 24] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
29
|
Hippler M, Oeltjen L, Quack M. High-Resolution Continuous-Wave-Diode Laser Cavity Ring-Down Spectroscopy of the Hydrogen Fluoride Dimer in a Pulsed Slit Jet Expansion: Two Components of the N = 2 Triad near 1.3 μm. J Phys Chem A 2007; 111:12659-68. [DOI: 10.1021/jp076894s] [Citation(s) in RCA: 24] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Michael Hippler
- Physical Chemistry, ETH Zürich, CH-8093 Zürich, Switzerland, and Department of Chemistry, University of Sheffield, Sheffield S3 7HF, England
| | - Lars Oeltjen
- Physical Chemistry, ETH Zürich, CH-8093 Zürich, Switzerland, and Department of Chemistry, University of Sheffield, Sheffield S3 7HF, England
| | - Martin Quack
- Physical Chemistry, ETH Zürich, CH-8093 Zürich, Switzerland, and Department of Chemistry, University of Sheffield, Sheffield S3 7HF, England
| |
Collapse
|
30
|
Baburao B, Visco DP, Albu TV. Association patterns in (HF)(m)(H2O)(n) (m + n = 2-8) clusters. J Phys Chem A 2007; 111:7940-56. [PMID: 17636970 DOI: 10.1021/jp0727648] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022]
Abstract
In an attempt to understand the phase behavior of aqueous hydrogen fluoride, the clustering in the mixture is investigated at the molecular level. The study is performed at the mPW1B95/6-31+G(d,p) level of theory. Several previous studies attempted to describe the dissociation of HF in water, but in this investigation, the focus is only on the association patterns that are present in this binary mixture. A total of 214 optimized geometries of (HF)n(H2O)m clusters, with m + n as high as 8, were investigated. For each cluster combination, several different conformations are investigated, and the preferred conformations are presented. Using multiple linear regressions, the average strengths of the four possible H-bonding interactions are obtained. The strongest H-bond interaction is reported to be the H2O...H-F interaction. The most probable distributions of mixed clusters as a function of composition are also deduced. It is found that the larger (HF)n(H2O)m clusters are favored both energetically and entropically compared to the ones that are of size m + n < or = 3. Also, the clusters with equimolar contributions of HF and H2O are found to have the strongest interactions.
Collapse
Affiliation(s)
- Barath Baburao
- Department of Chemical Engineering, Tennessee Technological University, Box 5013, Cookeville, Tennessee 38505, USA
| | | | | |
Collapse
|
31
|
Izvekov S, Voth GA. Effective force field for liquid hydrogen fluoride from ab initio molecular dynamics simulation using the force-matching method. J Phys Chem B 2007; 109:6573-86. [PMID: 16851738 DOI: 10.1021/jp0456685] [Citation(s) in RCA: 62] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
A recently developed force-matching method for obtaining effective force fields for condensed matter systems from ab initio molecular dynamics (MD) simulations has been applied to fit a simple nonpolarizable two-site pairwise force field for liquid hydrogen fluoride. The ab initio MD in this case was a Car-Parrinello (CP) MD simulation of 64 HF molecules at nearly ambient conditions within the Becke-Lee-Yang-Parr approximation to the electronic density functional theory. The force-matching procedure included a fit of short-ranged nonbonded forces, bonded forces, and atomic partial charges. The performance of the force-match potential was examined for the gas-phase dimer and for the liquid phase at various temperatures. The model was able to reproduce correctly the bent structure and energetics of the gas-phase dimer, while the results for the structural properties, self-diffusion, vibrational spectra, density, and thermodynamic properties of liquid HF were compared to both experiment and the CP MD simulation. The force-matching model performs well in reproducing nearly all of the liquid properties as well as the aggregation behavior at different temperatures. The model is computationally cheap and compares favorably to many more computationally expensive potential energy functions for liquid HF.
Collapse
Affiliation(s)
- Sergei Izvekov
- Department of Chemistry and Center for Biophysical Modeling and Simulation, University of Utah, 315 South 1400 East Room 2020, Salt Lake City, Utah 84112-0850, USA
| | | |
Collapse
|
32
|
Wang FF, Hou JH, Li ZR, Wu D, Li Y, Lu ZY, Cao WL. Unusual halogen-bonded complex FBrδ+⋯Brδ+F and hydrogen-bonded complex FBrδ+⋯Hδ+F formed by interactions between two positively charged atoms of different polar molecules. J Chem Phys 2007; 126:144301. [PMID: 17444706 DOI: 10.1063/1.2715559] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Using ab initio calculations, the authors' predicted for the first time that the halogen-bonded complex FBrdelta+...delta+BrF and hydrogen-bonded complex FBrdelta+...delta+HF formed by the interactions between two positively charged atoms of different polar molecules can be stable in gas phase. It shows that halogen bond or hydrogen bond not only exists between oppositely charged atoms but also between like-charged atoms. That the attraction arising from the special halogen bond or hydrogen bond can exceed the electrostatic repulsion between two contact positively charged atoms stabilizes the complex. Of course, from the point of view of physics they can consider the interactions in FBrdelta+...delta+BrF and FBrdelta+...delta+HF as mainly the sum of the long range molecular interactions, namely, electrostatic, induction, and dispersion with some short-range repulsion. They found that the intermolecular electron correlation contribution representing dispersion interaction plays a crucial role in the stabilities of seemingly repulsive complexes FBrdelta+...delta+BrF and FBrdelta+...delta+HF.
Collapse
Affiliation(s)
- Fang-Fang Wang
- State Key Laboratory of Theoretical and Computational Chemistry, Institute of Theoretical Chemistry, Jilin University, Changchun 130023, China
| | | | | | | | | | | | | |
Collapse
|
33
|
Küpper J, Merritt JM. Spectroscopy of free radicals and radical containing entrance-channel complexes in superfluid helium nanodroplets. INT REV PHYS CHEM 2007. [DOI: 10.1080/01442350601087664] [Citation(s) in RCA: 32] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/23/2022]
|
34
|
Vibrational Predissociation Dynamics of Van Der Waals Complexes: Product Rotational State Distributions. ADVANCES IN CHEMICAL PHYSICS 2007. [DOI: 10.1002/9780470141557.ch2] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register]
|
35
|
Isayev O, Gorb L, Leszczynski J. Theoretical calculations: Can Gibbs free energy for intermolecular complexes be predicted efficiently and accurately? J Comput Chem 2007; 28:1598-1609. [PMID: 17340602 DOI: 10.1002/jcc.20696] [Citation(s) in RCA: 31] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Abstract
The theoretical study has been performed to refine the procedure for calculations of Gibbs free energy with a relative accuracy of less than 1 kcal/mol. Three benchmark intermolecular complexes are examined via several quantum-chemical methods, including the second-order Moller-Plesset perturbation (MP2), coupled cluster (CCSD(T)), and density functional (BLYP, B3LYP) theories augmented by Dunnings correlation-consistent basis sets. The effects of electron correlation, basis set size, and anharmonicity are systematically analyzed, and the results are compared with available experimental data. The results of the calculations suggest that experimental accuracy can be reached only by extrapolation of MP2 and CCSD(T) total energies to the complete basis set. The contribution of anharmonicity to the zero point energy and TDeltaSint values is fairly small. The new, economic way to reach chemical accuracy in the calculations of the thermodynamic parameters of intermolecular interactions is proposed. In addition, interaction energy (De) and free energy change (DeltaA) for considered species have been evaluated by Carr-Parrinello molecular dynamics (CPMD) simulations and static BLYP-plane wave calculations. The free energy change along the reaction paths were determined by the thermodynamic integration/"Blue Moon Ensemble" technique. Comparison between obtained values, and available experimental and conventional ab initio results has been made. We found that the accuracy of CPMD simulations is affected by several factors, including statistical uncertainty and convergence of constrained forces (TD integration), and the nature of DFT (density functional theory) functional. The results show that CPMD technique is capable of reproducing interaction and free energy with an accuracy of 1 kcal/mol and 2-3 kcal/mol respectively.
Collapse
Affiliation(s)
- Olexandr Isayev
- Computational Center for Molecular Structure and Interactions, Jackson State University, Jackson, Mississippi 39217
| | - Leonid Gorb
- Computational Center for Molecular Structure and Interactions, Jackson State University, Jackson, Mississippi 39217
| | - Jerzy Leszczynski
- Computational Center for Molecular Structure and Interactions, Jackson State University, Jackson, Mississippi 39217
| |
Collapse
|
36
|
Swalina C, Wang Q, Chakraborty A, Hammes-Schiffer S. Analysis of Nuclear Quantum Effects on Hydrogen Bonding. J Phys Chem A 2007; 111:2206-12. [PMID: 17388289 DOI: 10.1021/jp0682661] [Citation(s) in RCA: 35] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
Abstract
The impact of nuclear quantum effects on hydrogen bonding is investigated for a series of hydrogen fluoride (HF)n clusters and a partially solvated fluoride anion, F-(H2O). The nuclear quantum effects are included using the path integral formalism in conjunction with the Car-Parrinello molecular dynamics (PICPMD) method and using the second-order vibrational perturbation theory (VPT2) approach. For the HF clusters, a directional change in the impact of nuclear quantum effects on the hydrogen-bonding strength is observed as the clusters evolve toward the condensed phase. Specifically, the inclusion of nuclear quantum effects increases the F-F distances for the (HF)n=2-4 clusters and decreases the F-F distances for the (HF)n>4 clusters. This directional change occurs because the enhanced electrostatic interactions between the HF monomers become more dominant than the zero point energy effects of librational modes as the size of the HF clusters increases. For the F-(H2O) system, the inclusion of nuclear quantum effects decreases the F-O distance and strengthens the hydrogen bonding interaction between the fluoride anion and the water molecule because of enhanced electrostatic interactions. The vibrationally averaged 19F shielding constant for F-(H2O) is significantly lower than the value for the equilibrium geometry, indicating that the electronic density on the fluorine decreases as a result of the quantum delocalization of the shared hydrogen. Deuteration of this system leads to an increase in the vibrationally averaged F-O distance and nuclear magnetic shielding constant because of the smaller degree of quantum delocalization for deuterium.
Collapse
Affiliation(s)
- Chet Swalina
- Department of Chemistry, 104 Chemistry Building, Pennsylvania State University, University Park, Pennsylvania 16802, USA
| | | | | | | |
Collapse
|
37
|
|
38
|
Inada Y, Orita H. Efficiency of numerical basis sets for predicting the binding energies of hydrogen bonded complexes: Evidence of small basis set superposition error compared to Gaussian basis sets. J Comput Chem 2007; 29:225-32. [PMID: 17565500 DOI: 10.1002/jcc.20782] [Citation(s) in RCA: 257] [Impact Index Per Article: 15.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Abstract
Binding energies of selected hydrogen bonded complexes have been calculated within the framework of density functional theory (DFT) method to discuss the efficiency of numerical basis sets implemented in the DFT code DMol3 in comparison with Gaussian basis sets. The corrections of basis set superposition error (BSSE) are evaluated by means of counterpoise method. Two kinds of different numerical basis sets in size are examined; the size of the one is comparable to Gaussian double zeta plus polarization function basis set (DNP), and that of the other is comparable to triple zeta plus double polarization functions basis set (TNDP). We have confirmed that the magnitudes of BSSE in these numerical basis sets are comparative to or smaller than those in Gaussian basis sets whose sizes are much larger than the corresponding numerical basis sets; the BSSE corrections in DNP are less than those in the Gaussian 6-311+G(3df,2pd) basis set, and those in TNDP are comparable to those in the substantially large scale Gaussian basis set aug-cc-pVTZ. The differences in counterpoise corrected binding energies between calculated using DNP and calculated using aug-cc-pVTZ are less than 9 kJ/mol for all of the complexes studied in the present work. The present results have shown that the cost effectiveness in the numerical basis sets in DMol3 is superior to that in Gaussian basis sets in terms of accuracy per computational cost.
Collapse
Affiliation(s)
- Yasuji Inada
- Accelrys K K, Nishishinbashi TS Building, 3-3-1 Nishishinbashi, Minato-ku, Tokyo 105-0003, Japan.
| | | |
Collapse
|
39
|
Giribet CG, Ruiz de Azúa MC. The Sign and Magnitude of 2hJ(F,F) and 1hJ(F,H) in FH···FH. A CLOPPA Analysis of Their Distance Dependence. J Phys Chem A 2006; 110:11575-83. [PMID: 17020272 DOI: 10.1021/jp060095c] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The sign change of the intermolecular (2h)J(F,F) coupling in the (HF)2 dimer as a function of the F-F distance is discussed by means of the CLOPPA method. It is found that it is due to the competition of positive and negative contributions involving the interaction of the sigma lone pair of the acceptor nucleus with vacant molecular orbitals localized in the F-H...F moiety and with other molecular orbitals localized in the donor molecule. The origin of the sign of each contribution is fully determined by analyzing the response of the electronic system to the magnetic perturbation at the acceptor F nucleus. (2h)J(F,F) coupling in the FH...F-, which is positive for all F-F distances, is also analyzed in order to look for the differences with the former case.
Collapse
Affiliation(s)
- Claudia G Giribet
- Department of Physics, Facultad de Ciencias Exactas y Naturales, University of Buenos Aires, Ciudad Universitaria, Pab. I, (1428) Buenos Aires, Argentina.
| | | |
Collapse
|
40
|
Avilés MW, Gray PT, Curotto E. Stereographic projection path-integral simulations of (HF)n clusters. J Chem Phys 2006; 124:174305. [PMID: 16689569 DOI: 10.1063/1.2192773] [Citation(s) in RCA: 16] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We perform several quantum canonical ensemble simulations of (HF)(n) clusters. The HF stretches are rigid, and the stereographic projection path-integral method is employed for the simulation in the resulting curved configuration space. We make use of the reweighted random series techniques to accelerate the convergence of the path-integral simulation with respect to the number of path coefficients. We develop and test estimators for the total energy and heat capacity based on a finite difference approach for non-Euclidean spaces. The quantum effects at temperatures below 400 K are substantial for all sizes. We observe interesting thermodynamic behaviors in the quantum simulations of the octamer and the heptamer.
Collapse
Affiliation(s)
- Michael W Avilés
- Department of Chemistry and Physics, Arcadia University, Glenside, Pennsylvania 19038, USA
| | | | | |
Collapse
|
41
|
Hobza P, Zahradník R, Müller-Dethlefs K. The World of Non-Covalent Interactions: 2006. ACTA ACUST UNITED AC 2006. [DOI: 10.1135/cccc20060443] [Citation(s) in RCA: 166] [Impact Index Per Article: 9.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/05/2022]
Abstract
The review focusses on the fundamental importance of non-covalent interactions in nature by illustrating specific examples from chemistry, physics and the biosciences. Laser spectroscopic methods and both ab initio and molecular modelling procedures used for the study of non-covalent interactions in molecular clusters are briefly outlined. The role of structure and geometry, stabilization energy, potential and free energy surfaces for molecular clusters is extensively discussed in the light of the most advanced ab initio computational results for the CCSD(T) method, extrapolated to the CBS limit. The most important types of non-covalent complexes are classified and several small and medium size non-covalent systems, including H-bonded and improper H-bonded complexes, nucleic acid base pairs, and peptides and proteins are discussed with some detail. Finally, we evaluate the interpretation of experimental results in comparison with state of the art theoretical models: this is illustrated for phenol...Ar, the benzene dimer and nucleic acid base pairs. A review with 270 references.
Collapse
|
42
|
Kim KH, Kim Y. The basis set superposition error in multilevel methods: a test on the H2O and HF dimer. Theor Chem Acc 2005. [DOI: 10.1007/s00214-005-0069-x] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
|
43
|
Swalina C, Hammes-Schiffer S. Impact of Nuclear Quantum Effects on the Molecular Structure of Bihalides and the Hydrogen Fluoride Dimer. J Phys Chem A 2005; 109:10410-7. [PMID: 16833338 DOI: 10.1021/jp053552i] [Citation(s) in RCA: 45] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022]
Abstract
The structural impact of nuclear quantum effects is investigated for a set of bihalides, [XHX](-), X = F, Cl, and Br, and the hydrogen fluoride dimer. Structures are calculated with the vibrational self-consistent-field (VSCF) method, the second-order vibrational perturbation theory method (VPT2), and the nuclear-electronic orbital (NEO) approach. In the VSCF and VPT2 methods, the vibrationally averaged geometries are calculated for the Born-Oppenheimer electronic potential energy surface. In the NEO approach, the hydrogen nuclei are treated quantum mechanically on the same level as the electrons, and mixed nuclear-electronic wave functions are calculated variationally with molecular orbital methods. Electron-electron and electron-proton dynamical correlation effects are included in the NEO approach using second-order perturbation theory (NEO-MP2). The nuclear quantum effects are found to alter the distances between the heavy atoms by 0.02-0.05 A for the systems studied. These effects are of similar magnitude as the electron correlation effects. For the bihalides, inclusion of the nuclear quantum effects with the NEO-MP2 or the VSCF method increases the X-X distance. The bihalide X-X distances are similar for both methods and are consistent with two-dimensional grid calculations and experimental values, thereby validating the use of the computationally efficient NEO-MP2 method for these types of systems. For the hydrogen fluoride dimer, inclusion of nuclear quantum effects decreases the F-F distance with the NEO-MP2 method and increases the F-F distance with the VSCF and VPT2 methods. The VPT2 F-F distances for the hydrogen fluoride dimer and the deuterated form are consistent with the experimentally determined values. The NEO-MP2 F-F distance is in excellent agreement with the distance obtained experimentally for a model that removes the large amplitude bending motions. The analysis of these calculations provides insight into the significance of electron-electron and electron-proton correlation, anharmonicity of the vibrational modes, and nonadiabatic effects for hydrogen-bonded systems.
Collapse
Affiliation(s)
- Chet Swalina
- Department of Chemistry, 104 Chemistry Building, The Pennsylvania State University, University Park, Pennsylvania 16802, USA
| | | |
Collapse
|
44
|
Raynaud C, Maron L, Jolibois F, Daudey JP, Esteves PM, Ramírez-Solís A. Ab initio molecular dynamics: Plane waves vs. local basis. Chem Phys Lett 2005. [DOI: 10.1016/j.cplett.2005.08.034] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/25/2022]
|
45
|
|
46
|
Reyes A, Pak MV, Hammes-Schiffer S. Investigation of isotope effects with the nuclear-electronic orbital approach. J Chem Phys 2005; 123:64104. [PMID: 16122297 DOI: 10.1063/1.1990116] [Citation(s) in RCA: 48] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
This paper addresses fundamental issues that arise in the application of the nuclear-electronic orbital (NEO) approach to systems with equivalent quantum nuclei. Our analysis illustrates that Hartree-Fock nuclear wave functions do not provide physically reasonable descriptions of systems comprised of equivalent low-spin fermions or equivalent bosons. The physical basis for this breakdown is that the ionic terms dominate due to the localized nature of the nuclear orbitals. Multi-configurational wave functions that include only covalent terms provide physically reasonable descriptions of these types of systems. The application of the NEO approach to a variety of chemical systems is presented to elucidate the isotope effects on the geometries and electronic wave functions. Deuteration of hydrogen halides, water, ammonia, and hydronium ion decreases the bond length and the magnitude of negative partial atomic charge on the heavy atom. These results are consistent with experimental spectroscopic data. Deuteration at the beta position for formate anion and a series of amines increases the magnitude of negative partial atomic charge on the protonation site for the unprotonated species. This observation is consistent with the experimentally observed increase in basicity upon deuteration at the beta position for carboxylic acids and amines.
Collapse
Affiliation(s)
- Andrés Reyes
- Department of Chemistry, 104 Chemistry Building, Pennsylvania State University, University Park, Pennsylvania 16802, USA
| | | | | |
Collapse
|
47
|
|
48
|
Abstract
The deltaK = 0 and 1 subbands of the (4,0) <-- (0,0) transition of (HF)2, near 14,700 cm(-1), have been measured by molecular-beam intracavity laser-induced fluorescence. The hydrogen interchange tunneling is basically quenched in (4, 0) for both K = 0 and 1 levels, consistent with the early suggestion from a phenomenological model [H.-C. Chang and W. Klemperer, J. Chem. Phys. 104, 7830 (1996)]. The band origin upsilon0 = 14,700.458(7) cm(-1) and rotational constant (B + C)/2 = 0.22278(31) cm(-1) are determined for K = 0 of the (4, 0) mode. From the observed deltaK = 1 <-- 0 spectrum, we determined that A = 24.3 cm(-1), (B + C)/2 = 0.22296(20) cm(-1), and (B-C) = 4.5(2) x 10(-3) cm(-1). The predissociation linewidths of both K = 0 and 1 levels are 470(30) MHz with no apparent rotational dependence.
Collapse
Affiliation(s)
- Zhenhong Yu
- Department of Chemistry and Chemical Biology Harvard University, Cambridge, Massachusetts 02138, USA
| | | | | |
Collapse
|
49
|
Kreitmeir M, Heusel G, Bertagnolli H, Tödheide K, Mundy CJ, Cuello GJ. Structure of dense hydrogen fluoride gas from neutron diffraction and molecular dynamics simulations. J Chem Phys 2005; 122:154511. [PMID: 15945649 DOI: 10.1063/1.1877232] [Citation(s) in RCA: 12] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The gas phase of hydrogen fluoride has been investigated by neutron diffraction experiments at three different particle densities. All investigated states are within the liquid-gas coexistence region of hydrogen fluoride. From the obtained diffraction data we deduced information about the local structure of the gas phase, which consists of small agglomerates. This has been expected as liquid hydrogen fluoride forms the strongest hydrogen bonds known. Molecular dynamics simulations with a modified potential have been carried out for all experimentally investigated states. The results confirmed that the size of the formed agglomerates in the gas phase is growing with increasing density of the gas phase.
Collapse
Affiliation(s)
- Markus Kreitmeir
- Institut für Physikalische Chemie, Universität Stuttgart, Pfaffenwaldring 55, D-70569 Stuttgart, Germany.
| | | | | | | | | | | |
Collapse
|
50
|
Benedek NA, Snook IK, Latham K, Yarovsky I. Application of numerical basis sets to hydrogen bonded systems: A density functional theory study. J Chem Phys 2005; 122:144102. [PMID: 15847511 DOI: 10.1063/1.1876152] [Citation(s) in RCA: 77] [Impact Index Per Article: 4.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We have investigated and compared the ability of numerical and Gaussian-type basis sets to accurately describe the geometries and binding energies of a selection of hydrogen bonded systems that are well studied theoretically and experimentally. The numerical basis sets produced accurate results for geometric parameters but tended to overestimate binding energies. However, a comparison of the time taken to optimize phosphinic acid dimer, the largest complex considered in this study, shows that calculations using numerical basis sets offer a definitive advantage where geometry optimization of large systems is required.
Collapse
Affiliation(s)
- N A Benedek
- School of Applied Sciences, Applied Physics and Applied Chemistry, RMIT University, GPO Box 2476V, Melbourne 3001, Australia
| | | | | | | |
Collapse
|