1
|
Kříž K, van der Spoel D. Quantification of Anisotropy in Exchange and Dispersion Interactions: A Simple Model for Physics-Based Force Fields. J Phys Chem Lett 2024; 15:9974-9978. [PMID: 39314113 DOI: 10.1021/acs.jpclett.4c02034] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/25/2024]
Abstract
In some compounds, exchange repulsion is orientation dependent. However, in contrast to quantum chemical methods that treat exchange explicitly, empirical models assume exchange to be spherically symmetric, yielding an average description only. Here we quantify the anisotropy of exchange and dispersion energy for hydrogen halides and water by probing these compounds with a helium atom using the symmetry-adapted perturbation theory (SAPT). The exchange interaction is reduced by up to 33% due to the σ-hole in hydrogen iodide, depending on the location of the probe. We demonstrate how this anisotropy can be modeled in empirical force fields either using an angle-dependent potential or by introducing virtual sites, reducing the error in the empirical model by a factor of 5 compared to isotropic atoms. Lone-pairs on water, positioned close to perpendicular to the plane of the molecule, on a line with the oxygen atom, and, surprisingly, σ-holes on water both modulate the exchange interaction strongly. Both lone-pairs and σ-holes can be modeled by virtual sites, leading to an 80% reduced error.
Collapse
Affiliation(s)
- Kristian Kříž
- Science for Life Laboratory, Department of Cell and Molecular Biology, Uppsala University, Husargatan 3, Box 596, SE-75124 Uppsala, Sweden
| | - David van der Spoel
- Science for Life Laboratory, Department of Cell and Molecular Biology, Uppsala University, Husargatan 3, Box 596, SE-75124 Uppsala, Sweden
| |
Collapse
|
2
|
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
|
3
|
Savin AV, Kivshar YS. Stabilization of hydrogen-bonded molecular chains by carbon nanotubes. CHAOS (WOODBURY, N.Y.) 2024; 34:043111. [PMID: 38572948 DOI: 10.1063/5.0197401] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/13/2024] [Accepted: 03/14/2024] [Indexed: 04/05/2024]
Abstract
We study numerically nonlinear dynamics of several types of molecular systems composed of hydrogen-bonded chains placed inside carbon nanotubes with open edges. We demonstrate that carbon nanotubes provide a stabilization mechanism for quasi-one-dimensional molecular chains via the formation of their secondary structures. In particular, a polypeptide chain (Gly)N placed inside a carbon nanotube can form a stable helical chain (310-, α-, π-, and β-helix) with parallel chains of hydrogen-bonded peptide groups. A chain of hydrogen fluoride molecules ⋯FH⋯FH⋯FH can form a hydrogen-bonded zigzag chain. Remarkably, we demonstrate that for molecular complexes (Gly)N∈CNT and (FH)N∈CNT, the hydrogen-bonded chains will remain stable even at T=500 K. Thus, our results suggest that the use of carbon nanotubes with encapsulated hydrogen fluoride molecules may be important for the realization of high proton conductivity at high temperatures.
Collapse
Affiliation(s)
- Alexander V Savin
- Semenov Institute of Chemical Physics, Russian Academy of Sciences, Moscow 119991, Russia
- Academic Department of Innovational Materials and Technologies Chemistry, Plekhanov Russian University of Economics, Moscow 117997, Russia
| | - Yuri S Kivshar
- Nonlinear Physics Center, Research School of Physics, Australian National University, Canberra, ACT 2601, Australia
| |
Collapse
|
4
|
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
|
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
|
Pethes I, Bakó I, Pusztai L. Chloride ions as integral parts of hydrogen bonded networks in aqueous salt solutions: the appearance of solvent separated anion pairs. Phys Chem Chem Phys 2020; 22:11038-11044. [DOI: 10.1039/d0cp01806f] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023]
Abstract
Characteristic hydrogen bonded motifs, including solvent separated anion pairs, in concentrated aqueous LiCl solutions.
Collapse
Affiliation(s)
- Ildikó Pethes
- Wigner Research Centre for Physics
- Konkoly Thege út 29-33
- H-1121 Budapest
- Hungary
| | - Imre Bakó
- Research Centre for Natural Sciences
- Magyar tudósok körútja 2
- H-1117 Budapest
- Hungary
| | - László Pusztai
- Wigner Research Centre for Physics
- Konkoly Thege út 29-33
- H-1121 Budapest
- Hungary
- International Research Organisation for Advanced Science and Technology (IROAST)
| |
Collapse
|
7
|
Experimental and theoretical assessment of the interactions of ionic liquids (ILs) with fluoridated compounds (HF, R-F) in organic medium. J Mol Liq 2019. [DOI: 10.1016/j.molliq.2018.12.040] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
|
8
|
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
|
9
|
Christiansen PL, Savin AV, Zolotaryuk AV. Lattice stretching bistability and dynamic heterogeneity. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2012; 85:046601. [PMID: 22680592 DOI: 10.1103/physreve.85.046601] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/03/2011] [Revised: 02/20/2012] [Indexed: 06/01/2023]
Abstract
A simple one-dimensional lattice model is suggested to describe the experimentally observed plateau in force-stretching diagrams for some macromolecules. This chain model involves the nearest-neighbor interaction of a Morse-like potential (required to have a saturation branch) and a harmonic second-neighbor coupling. Under an external stretching applied to the chain ends, the intersite Morse-like potential results in the appearance of a double-well potential within each chain monomer, whereas the interaction between the second neighbors provides a homogeneous bistable (degenerate) ground state, at least within a certain part of the chain. As a result, different conformational changes occur in the chain under the external forcing. The transition regions between these conformations are described as topological solitons. With a strong second-neighbor interaction, the solitons describe the transition between the bistable ground states. However, the key point of the model is the appearance of a heterogenous structure, when the second-neighbor coupling is sufficiently weak. In this case, a part of the chain has short bonds with a single-well potential, whereas the complementary part admits strongly stretched bonds with a double-well potential. This case allows us to explain the existence of a plateau in the force-extension diagram for DNA and α-helix protein. Finally, the soliton dynamics are studied in detail.
Collapse
Affiliation(s)
- P L Christiansen
- Department of Informatics, Technical University of Denmark, Kgs. Lyngby, Denmark
| | | | | |
Collapse
|
10
|
KOVALENKO ANDRIY, HIRATA FUMIO. TOWARDS A MOLECULAR THEORY FOR THE VAN DER WAALS–MAXWELL DESCRIPTION OF FLUID PHASE TRANSITIONS. JOURNAL OF THEORETICAL & COMPUTATIONAL CHEMISTRY 2012. [DOI: 10.1142/s0219633602000282] [Citation(s) in RCA: 38] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
Abstract
We briefly review developments of theories for phase transitions of molecular fluids and mixtures, from semi-phenomenological approaches providing equations of state with adjustable parameters to first-principles microscopic methods qualitatively correct for a variety of molecular models with realistic interaction potentials. We further present the generalization of the van der Waals–Maxwell description of fluid phase diagrams to account for chemical specificities of polar molecular fluids, such as hydrogen bonding. Our theory uses the reference interaction site model (RISM) integral equation approach complemented with the new closure we have proposed (KH approximation), successful over a wide range of density from gas to liquid. The RISM/KH theory is applied to the known three-site models of water, methanol, and hydrogen fluoride. It qualitatively reproduces their vapor-liquid phase diagrams and the structure in the gas as well as liquid phases, including hydrogen bonding. Furthermore, phase transitions of water and methanol sorbed in nanoporous carbon aerogel are described by means of the replica generalization of the RISM approach we have developed. The changes as compared to the bulk fluids are in agreement with simulations and experiment. The RISM/KH theory is promising for description of phase transitions in various associating fluids, in particular, electrolyte as well as non-electrolyte solutions and ionic liquids.
Collapse
Affiliation(s)
- ANDRIY KOVALENKO
- Institute for Molecular Science, Myodaiji, Okazaki 444-8585, Japan
| | - FUMIO HIRATA
- Institute for Molecular Science, Myodaiji, Okazaki 444-8585, Japan
| |
Collapse
|
11
|
Assessing group-based cutoffs and the Ewald method for electrostatic interactions in clusters and in saturated, superheated, and supersaturated vapor phases of dipolar molecules. Theor Chem Acc 2011. [DOI: 10.1007/s00214-011-0973-1] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
|
12
|
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
|
13
|
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]
|
14
|
JEDLOVSZKY BPAL, VALLAURI RENZO. Computer simulation study of liquid HF with a new effective pair potential model. Mol Phys 2010. [DOI: 10.1080/002689797170536] [Citation(s) in RCA: 34] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/17/2022]
|
15
|
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]
|
16
|
Bhargava BL, Balasubramanian S. Ab Initio Molecular Dynamics Simulation of a 1-Ethyl-3-methylimidazolium Fluoride−Hydrogen Fluoride Mixture. J Phys Chem B 2008; 112:7566-73. [DOI: 10.1021/jp801323g] [Citation(s) in RCA: 26] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- B. L. Bhargava
- Chemistry and Physics of Materials Unit, Jawaharlal Nehru Centre for Advanced Scientific Research, Jakkur, Bangalore 560 064, India
| | - S. Balasubramanian
- Chemistry and Physics of Materials Unit, Jawaharlal Nehru Centre for Advanced Scientific Research, Jakkur, Bangalore 560 064, India
| |
Collapse
|
17
|
Investigation of the saturated adsorption layer of 5-cyano-biphenyl and 5-cyano-terphenyl at the free water surface by Monte Carlo simulation. J Mol Liq 2007. [DOI: 10.1016/j.molliq.2007.08.013] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
|
18
|
Woittequand S, Duflot D, Monnerville M, Pouilly B, Toubin C, Briquez S, Meyer HD. Classical and quantum studies of the photodissociation of a HX (X=Cl,F) molecule adsorbed on ice. J Chem Phys 2007; 127:164717. [DOI: 10.1063/1.2799519] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
|
19
|
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
|
20
|
Jedlovszky P, Pártay LB. Adsorption of Octyl Cyanide at the Free Water Surface as Studied by Monte Carlo Simulation. J Phys Chem B 2007; 111:5885-95. [PMID: 17487997 DOI: 10.1021/jp068566i] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Monte Carlo simulations of the adsorption layer of octyl cyanide have been performed on the canonical (N, V, T) ensemble at 300 K. The systems simulated cover the range of octyl cyanide surface densities from 0.27 to 7.83 mumol/m2. The surface density value at which the saturation of the adsorption layer occurs is estimated to be 1.7 mumol/m2. At low surface densities, the main driving force of the adsorption is found to be the formation of hydrogen bonds between the water and octyl cyanide molecules, whereas at higher surface concentrations, the dipole-dipole attraction between the neighboring adsorbed octyl cyanide molecules becomes more important. At low surface concentrations, the water-octyl cyanide hydrogen bonds prefer tilted alignments relative to the interface; however, in the case of the saturated adsorption layer, the number of such hydrogen bonds is maximized, leading to the preference of these bonds for the orientation perpendicular to the interface. Contrary to nonionic surfactants of multiple hydrogen bonding abilities (e.g., 1-octanol, C8E3), the increasing surface concentration of octyl cyanide was not found to lead to considerable competition of the molecules for positions of optimal arrangement. As a consequence, the energy and geometry of the water-octyl cyanide hydrogen bonds are found to be insensitive to the octyl cyanide surface concentration.
Collapse
Affiliation(s)
- Pal Jedlovszky
- Laboratory of Interfaces and Nanosize Systems, Institute of Chemistry, Eötvös Loránd University, Pázmány Péter stny. 1/a, Budapest, Hungary.
| | | |
Collapse
|
21
|
Dubois V, Bourasseau E, Maillet JB. New potential model for molecular dynamic simulation of liquid HF. II – Parameter optimization for repulsion-dispersion potential. Mol Phys 2007. [DOI: 10.1080/00268970601148258] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/23/2022]
|
22
|
McLain SE, Benmore CJ, Siewenie JE, Molaison JJ, Turner JFC. On the variation of the structure of liquid deuterium fluoride with temperature. J Chem Phys 2007; 121:6448-55. [PMID: 15446944 DOI: 10.1063/1.1790432] [Citation(s) in RCA: 11] [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 structure of liquid deuterium fluoride has been measured using pulsed neutron diffraction and high energy x-ray diffraction techniques as a function of temperature. The neutron experiments were performed at T=296+/-2 K, 246+/-2 K, and 193+/-2 K and the x-ray measurements carried out at 296+/-2 K and 195+/-2 K. The x-ray pair correlation functions, which are dominated by fluorine-fluorine interactions, show the first peak at approximately 2.53+/-0.05 A remains very nearly invariant with decreasing temperature. Peaks around 4.5 and 5.0 A also appear at both temperatures in the x-ray data. In contrast, the intermolecular peaks in the total neutron pair correlation function show that significant systematic local structural changes occur as the temperature is lowered. The first intermolecular peak position shortens from 1.64+/-0.05 A at 296 K to 1.56+/-0.05 A at 195 K. Although there are overlapping contributions from the intermolecular hydrogen-fluorine and hydrogen-hydrogen correlations, it is clear that the temperature dependent structural changes are largely due to a rearrangement of the deuterium atom positions in the fluid. By comparison with partial structure factor data the hydrogen bonds appear to become more linear at lower temperatures.
Collapse
Affiliation(s)
- S E McLain
- Department of Chemistry, University of Tennessee, Knoxville, Tennessee 37996-1600, USA
| | | | | | | | | |
Collapse
|
23
|
Deraman M, Dore J, Powles J, Holloway J, Chieux P. Structural studies of liquid hydrogen fluoride by neutron diffraction. Mol Phys 2006. [DOI: 10.1080/00268978500102061] [Citation(s) in RCA: 42] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/24/2022]
|
24
|
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
|
25
|
Wierzchowski SJ, Fang ZH, Kofke DA, Tilson JL. Three-body effects in hydrogen fluoride: survey of potential energy surfaces. Mol Phys 2006. [DOI: 10.1080/00268970500424321] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/25/2022]
|
26
|
Bourasseau E, Maillet JB, Mondelain L, Anglade PM. New potential model for molecular dynamic simulation of liquid HF. I—Parameter optimization for charge equilibration method. MOLECULAR SIMULATION 2005. [DOI: 10.1080/08927020500231817] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/23/2022]
|
27
|
Vlček L, Nezbeda * I. From realistic to simple models of fluids. III. Primitive models of carbon dioxide, hydrogen sulphide and acetone, and their properties. Mol Phys 2005. [DOI: 10.1080/00268970500083630] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/25/2022]
|
28
|
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
|
29
|
Coarse-grained approach to nonequilibrium molecular dynamics: Application to relaxation process of vibrationally excited hydrogen fluoride in aqueous solution. Chem Phys Lett 2005. [DOI: 10.1016/j.cplett.2005.03.129] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/23/2022]
|
30
|
Chialvo AA, Kettler M, Nezbeda I. Effect of the Range of Interactions on the Properties of Fluids. 2. Structure and Phase Behavior of Acetonitrile, Hydrogen Fluoride, and Formic Acid. J Phys Chem B 2005; 109:9736-50. [PMID: 16852173 DOI: 10.1021/jp050922u] [Citation(s) in RCA: 17] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
To complete the study on the effect of the long-range part of Coulombic interactions on properties of complex polar and associating fluids, we have investigated in detail three compounds with extreme features: acetonitrile for its unusually large dipole moment, hydrogen fluoride with very strong hydrogen bonding, and formic acid for its potential formation of different n-mers in liquid and gaseous phases. The effect of the long-range Coulombic interactions on both the structure and thermodynamics of the homogeneous phase, and on the vapor-liquid equilibria has been examined using the same decomposition of realistic potential models into a short-range part and a residual part as in the previous paper [Kettler, M.; et al. J. Phys. Chem. B 2002, 106, 7537-7546]. The present results fully confirm the previous findings that the properties of polar and associating systems are determined primarily by the short-range interactions regardless of their nature, i.e., contributions arising from the long-range interactions constitute only a small portion of the total properties, and thus that the short-range potential counterpart of full realistic models can be used as a convenient reference for a successful perturbation expansion.
Collapse
Affiliation(s)
- Ariel A Chialvo
- Chemical Sciences Division, Aqueous Chemistry and Geochemistry Group, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831-6110, USA.
| | | | | |
Collapse
|
31
|
|
32
|
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
|
33
|
Vl[cbreve]ek L, Nezbeda I. From realistic to simple models of associating fluids. II. Primitive models of ammonia, ethanol and models of water revisited. Mol Phys 2004. [DOI: 10.1080/00268970410001668417] [Citation(s) in RCA: 20] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/26/2022]
|
34
|
Wierzchowski SJ, Kofke DA. Liquid-Phase Activity Coefficients for Saturated HF/H2O Mixtures with Vapor-Phase Nonidealities Described by Molecular Simulation. Ind Eng Chem Res 2003. [DOI: 10.1021/ie030437q] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Scott J. Wierzchowski
- Department of Chemical Engineering, University at Buffalo, The State University of New York, Buffalo, New York 14260-4200
| | - David A. Kofke
- Department of Chemical Engineering, University at Buffalo, The State University of New York, Buffalo, New York 14260-4200
| |
Collapse
|
35
|
|
36
|
Wierzchowski SJ, Kofke DA, Gao J. Hydrogen fluoride phase behavior and molecular structure: A QM/MM potential model approach. J Chem Phys 2003. [DOI: 10.1063/1.1607919] [Citation(s) in RCA: 26] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
37
|
Liem SY, Popelier PLA. High-rank quantum topological electrostatic potential: Molecular dynamics simulation of liquid hydrogen fluoride. J Chem Phys 2003. [DOI: 10.1063/1.1593012] [Citation(s) in RCA: 41] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
38
|
Vissers GWM, Groenenboom GC, van der Avoird A. Spectrum and vibrational predissociation of the HF dimer. I. Bound and quasibound states. J Chem Phys 2003. [DOI: 10.1063/1.1577111] [Citation(s) in RCA: 25] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
39
|
Toubin C, Picaud S, Hoang PNM, Girardet C, Lynden-Bell RM, Hynes JT. Adsorption of HF and HCl molecules on ice at 190 and 235 K from molecular dynamics simulations: Free energy profiles and residence times. J Chem Phys 2003. [DOI: 10.1063/1.1570408] [Citation(s) in RCA: 22] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
40
|
Muñoz-Losa A, Fdez.-Galván I, Martín ME, Aguilar MA. Theoretical Study of Liquid Hydrogen Fluoride. Application of the Averaged Solvent Electrostatic Potential/Molecular Dynamics Method. J Phys Chem B 2003. [DOI: 10.1021/jp022422w] [Citation(s) in RCA: 13] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- A. Muñoz-Losa
- Departamento de Química Física, Universidad de Extremadura, Avda de Elvas, s/n 06071 Badajoz, Spain
| | - I. Fdez.-Galván
- Departamento de Química Física, Universidad de Extremadura, Avda de Elvas, s/n 06071 Badajoz, Spain
| | - M. E. Martín
- Departamento de Química Física, Universidad de Extremadura, Avda de Elvas, s/n 06071 Badajoz, Spain
| | - M. A. Aguilar
- Departamento de Química Física, Universidad de Extremadura, Avda de Elvas, s/n 06071 Badajoz, Spain
| |
Collapse
|
41
|
Kreitmeir M, Bertagnolli H, Mortensen JJ, Parrinello M. Ab initio molecular dynamics simulation of hydrogen fluoride at several thermodynamic states. J Chem Phys 2003. [DOI: 10.1063/1.1539045] [Citation(s) in RCA: 44] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
42
|
Honda K. An effective potential function with enhanced charge-transfer-type interaction for hydrogen-bonding liquids. J Chem Phys 2002. [DOI: 10.1063/1.1495851] [Citation(s) in RCA: 14] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
43
|
Quack M, Stohner J, Suhm MA. Analytical three-body interaction potentials and hydrogen bond dynamics of hydrogen fluoride aggregates, (HF) n , n ≥3. J Mol Struct 2001. [DOI: 10.1016/s0022-2860(01)00825-0] [Citation(s) in RCA: 36] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/27/2022]
|
44
|
Jedlovszky P, Mezei M, Vallauri R. Comparison of polarizable and nonpolarizable models of hydrogen fluoride in liquid and supercritical states: A Monte Carlo simulation study. J Chem Phys 2001. [DOI: 10.1063/1.1413973] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
|
45
|
Chen B, Siepmann JI. Improving the Efficiency of the Aggregation−Volume−Bias Monte Carlo Algorithm. J Phys Chem B 2001. [DOI: 10.1021/jp012209k] [Citation(s) in RCA: 107] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Bin Chen
- Department of Chemistry and Department of Chemical Engineering and Materials Science, University of Minnesota, 207 Pleasant Street SE, Minneapolis, Minnesota 55455-0431
| | - J. Ilja Siepmann
- Department of Chemistry and Department of Chemical Engineering and Materials Science, University of Minnesota, 207 Pleasant Street SE, Minneapolis, Minnesota 55455-0431
| |
Collapse
|
46
|
Čuma M, Schmitt UW, Voth GA. A Multi-State Empirical Valence Bond Model for Weak Acid Dissociation in Aqueous Solution. J Phys Chem A 2001. [DOI: 10.1021/jp0038207] [Citation(s) in RCA: 73] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Martin Čuma
- Department of Chemistry and Henry Eyring Center for Theoretical Chemistry, 315 S. 1400 E. Rm 2020, University of Utah, Salt Lake City, Utah 84112-0850
| | - Udo W. Schmitt
- Department of Chemistry and Henry Eyring Center for Theoretical Chemistry, 315 S. 1400 E. Rm 2020, University of Utah, Salt Lake City, Utah 84112-0850
| | - Gregory A. Voth
- Department of Chemistry and Henry Eyring Center for Theoretical Chemistry, 315 S. 1400 E. Rm 2020, University of Utah, Salt Lake City, Utah 84112-0850
| |
Collapse
|
47
|
Martı́n C, Lombardero M, Anta JA, Lomba E. Integral equation study of liquid hydrogen fluoride. J Chem Phys 2001. [DOI: 10.1063/1.1328758] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
|
48
|
Pfleiderer T, Waldner I, Bertagnolli H, Tödheide K, Fischer HE. The structure of liquid and supercritical deuterium fluoride from neutron scattering using high-pressure techniques. J Chem Phys 2000. [DOI: 10.1063/1.1287427] [Citation(s) in RCA: 36] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
49
|
Chen B, Siepmann JI. A Novel Monte Carlo Algorithm for Simulating Strongly Associating Fluids: Applications to Water, Hydrogen Fluoride, and Acetic Acid. J Phys Chem B 2000. [DOI: 10.1021/jp001952u] [Citation(s) in RCA: 146] [Impact Index Per Article: 6.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Bin Chen
- Department of Chemistry and Department of Chemical Engineering and Materials Science, University of Minnesota, 207 Pleasant Street SE, Minneapolis, Minnesota 55455-0431
| | - J. Ilja Siepmann
- Department of Chemistry and Department of Chemical Engineering and Materials Science, University of Minnesota, 207 Pleasant Street SE, Minneapolis, Minnesota 55455-0431
| |
Collapse
|
50
|
|