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
|
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
|
3
|
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
|
4
|
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
|
5
|
Abstract
This paper summarizes the scientific trends associated with the rapid development of the technique of high-energy X-ray diffraction over the past decade pertaining to the field of liquids, glasses, and amorphous materials. The measurement of high-quality X-ray structure factors out to large momentum transfers leads to high-resolution pair distribution functions which can be directly compared to theory or combined with data from other experimental techniques. The advantages of combining highly penetrating radiation with low angle scattering are outlined together with the data analysis procedure and formalism. Also included are advances in high-energy synchrotron beamline instrumentation, sample environment equipment, and an overview of the role of simulation and modeling for interpreting data from disordered materials. Several examples of recent trends in glass and liquid research are described. Finally, directions for future research are considered within the context of past and current developments in the field.
Collapse
Affiliation(s)
- C. J. Benmore
- Department of Physics, Arizona State University, Tempe, AZ 85287-1604, USA
| |
Collapse
|
6
|
McGrath MJ, Kuo IFW, Siepmann JI. Liquid structures of water, methanol, and hydrogen fluoride at ambient conditions from first principles molecular dynamics simulations with a dispersion corrected density functional. Phys Chem Chem Phys 2011; 13:19943-50. [DOI: 10.1039/c1cp21890e] [Citation(s) in RCA: 59] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
7
|
Giura P, Angelini R, Datchi F, Ruocco G, Sette F. High frequency dynamics and structural relaxation process in liquid ammonia. J Chem Phys 2007; 127:084508. [PMID: 17764270 DOI: 10.1063/1.2753161] [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
The dynamic structure factor S(Q,omega) of liquid ammonia has been measured by inelastic x-ray scattering in the terahertz frequency region as a function of the temperature in the range of 220-298 K at a pressure P=85 bars. The data have been analyzed using the generalized hydrodynamic formalism with a three term memory function to take into account the thermal, the structural, (alpha) and the microscopic (mu) relaxation processes affecting the dynamics of the liquid. This allows to extract the temperature dependence of the structural relaxation time (tau(alpha)) and strength (Delta(alpha)). The former quantity follows an Arrhenius behavior with an activation energy E(a)=2.6+/-0.2 kcal/mol, while the latter is temperature independent suggesting that there are no changes in the interparticle potential and arrangement with T. The obtained results, compared with those already existing in liquid water and liquid hydrogen fluoride, suggest the strong influence of the connectivity of the molecular network on the structural relaxation.
Collapse
Affiliation(s)
- P Giura
- IMPMC, Université Pierre et Marie Curie, Paris 6, UMR7590, 140 Rue de Lourmel, Paris F-75015, France
| | | | | | | | | |
Collapse
|
8
|
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
|
9
|
Fernandez-Alonso F, Bermejo FJ, McLain SE, Turner JFC, Molaison JJ, Herwig KW. Observation of fractional Stokes-Einstein behavior in the simplest hydrogen-bonded liquid. PHYSICAL REVIEW LETTERS 2007; 98:077801. [PMID: 17359062 DOI: 10.1103/physrevlett.98.077801] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/24/2006] [Indexed: 05/14/2023]
Abstract
Quasielastic neutron scattering has been used to investigate the single-particle dynamics of hydrogen fluoride across its entire liquid range at ambient pressure. For T>230 K, translational diffusion obeys the celebrated Stokes-Einstein relation, in agreement with nuclear magnetic resonance studies. At lower temperatures, we find significant deviations from the above behavior in the form of a power law with exponent xi=-0.71+/-0.05. More striking than the above is a complete breakdown of the Debye-Stokes-Einstein relation for rotational diffusion. Our findings provide the first experimental verification of fractional Stokes-Einstein behavior in a hydrogen-bonded liquid, in agreement with recent computer simulations [S. R. Becker, Phys. Rev. Lett. 97, 055901 (2006)10.1103/PhysRevLett.97.055901].
Collapse
Affiliation(s)
- F Fernandez-Alonso
- ISIS Facility, Rutherford Appleton Laboratory, Chilton, Didcot, Oxfordshire OX11 0QX, United Kingdom.
| | | | | | | | | | | |
Collapse
|
10
|
Bermejo FJ, Taylor JW, McLain SE, Bustinduy I, Turner JFC, Ruiz-Martin MD, Cabrillo C, Fernandez-Perea R. Evidence of the presence of opticlike collective modes in a liquid from neutron scattering experiments. PHYSICAL REVIEW LETTERS 2006; 96:235501. [PMID: 16803381 DOI: 10.1103/physrevlett.96.235501] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/24/2005] [Indexed: 05/10/2023]
Abstract
Inelastic neutron scattering data from liquid DF close to the melting point show, in addition to spectra comprising quasielastic and heavily damped acoustic motions, an intense, nondispersive band centered at about 27 meV along with a broader higher energy feature. Observation of the former band provides the first direct verification of the existence within the liquid state of collective opticlike excitations as predicted by molecular dynamics simulations. The latter corresponds to mainly reorientational motions assigned from mode eigenvector analysis carried out by computer simulations.
Collapse
Affiliation(s)
- F J Bermejo
- CSIC, Department of Electricity and Electronics, University Basque Country, P.O. Box 644, Bilbao 48080, Spain
| | | | | | | | | | | | | | | |
Collapse
|
11
|
Pártay L, Jedlovszky P, Vallauri R. Development of a new polarizable potential model of hydrogen fluoride and comparison with other effective models in liquid and supercritical states. J Chem Phys 2006; 124:184504. [PMID: 16709123 DOI: 10.1063/1.2192771] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
Development of a new polarizable potential of hydrogen fluoride through the reparametrization of the JV-P model is presented: The length of the H-F bond has been shortened and the other parameters of the model have been readjusted accordingly. The structural, thermodynamic, and liquid-vapor equilibrium properties of the new model are compared with those of other effective potential models of HF as well as with experimental data in a broad range of thermodynamic states, from near-freezing to supercritical conditions. It is found that although the reparametrization does not change the structural properties of the HF model noticeably at the level of the pair correlations, it improves the reproduction of the thermodynamic properties of hydrogen fluoride over the entire range of existence of a thermodynamically stable liquid phase and also that of the vapor-liquid coexistence curve. However, the new model, which still overestimates the close-contact separation of the HF molecules, underestimates the density of the coexisting liquid phase and overestimates the saturation pressure, probably due to the too steep repulsion of the potential function.
Collapse
Affiliation(s)
- Lívia Pártay
- Institute of Chemistry, Eötvös Loránd University, Pázmány Péter sétány 1/a, H-1117 Budapest, Hungary
| | | | | |
Collapse
|
12
|
Wierzchowski SJ, Kofke DA. Trimer Based Polarization as a Multibody Molecular Model. Application to Hydrogen Fluoride. J Am Chem Soc 2004; 127:690-8. [PMID: 15643894 DOI: 10.1021/ja031877b] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
A molecular modeling approach is introduced as a way to treat multibody (more than two molecules) contributions to the intermolecular potential. There are two key features to the method. First, it employs polarizable electrostatics on the molecules, but converges the charges and fields for only three molecules at a time, taken separately for all trimers (three molecules falling within a cutoff distance) in the system. This feature introduces significant computational savings when applied in Monte Carlo simulation (in comparison to a full N-body polarization treatment), as movement of a single molecule does not require re-converging of the polarization of all molecules, and it achieves this without approximations that cause the value of the energy to depend on the history of the simulation. Second, the approach defines the polarization energy in excess of the pairwise contribution, meaning that the trimer energy has subtracted from it the sum of the energies obtained by converging the polarization of each molecule pair in the trimer. This feature is advantageous because it removes the need (often found in polarizable models) to stiffen inappropriately the repulsive part of the pair potential. The polarization contribution is thus a purely three-body potential. The approach is applied to model hydrogen fluoride, which in experiments exhibits unusual properties that have proven difficult to capture well by molecular models. The new HF model is shown to be much more successful than previous modeling efforts in obtaining agreement with a broad range of experimental data (volumetric properties, heat effects, molecular structure, and vapor-liquid equilibria).
Collapse
Affiliation(s)
- Scott J Wierzchowski
- Department of Chemical and Biological Engineering, University at Buffalo, The State University of New York, Buffalo, New York 14260-4200, USA
| | | |
Collapse
|