1
|
Hrahsheh F, Jum'h I, Wilemski G. Second inflection point of supercooled water surface tension induced by hydrogen bonds: A molecular-dynamics study. J Chem Phys 2024; 160:114504. [PMID: 38506292 DOI: 10.1063/5.0185832] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/02/2023] [Accepted: 02/26/2024] [Indexed: 03/21/2024] Open
Abstract
Surface tension of supercooled water is a fundamental property in various scientific processes. In this study, we perform molecular dynamics simulations with the TIP4P-2005 model to investigate the surface tension of supercooled water down to 220 K. Our results show a second inflection point (SIP) in the surface tension at temperature TSIP ≈ 267.5 ± 2.3 K. Using an extended IAPWS-E functional fit for the water surface tension, we calculate the surface excess internal-energy and entropy terms of the excess Helmholtz free energy. Similar to prior studies [Wang et al., Phys. Chem. Chem. Phys. 21, 3360 (2019); Gorfer et al., J. Chem. Phys. 158, 054503 (2023)], our results show that the surface tension is governed by two driving forces: a surface excess entropy change above the SIP and a surface excess internal-energy change below it. We study hydrogen-bonding near the SIP because it is the main cause of water's anomalous properties. With decreasing temperature, our results show that the entropy contribution to the surface tension reaches a maximum slightly below the SIP and then decreases. This is because the number of hydrogen bonds increases more slowly below the SIP. Moreover, the strengths and lifetimes of the hydrogen bonds also rise dramatically below the SIP, causing the internal-energy term to dominate the excess surface free energy. Thus, the SIP in the surface tension of supercooled TIP4P-2005 water is associated with an increase in the strengths and lifetimes of hydrogen bonds, along with a decrease in the formation rate (#/K) of new hydrogen bonds.
Collapse
Affiliation(s)
- Fawaz Hrahsheh
- Higher Colleges of Technology, ETS, MZWC, Abu Dhabi 25026, United Arab Emirates
| | - Inshad Jum'h
- School of Basic Sciences and Humanities, German Jordanian University, Amman 11180, Jordan
| | - Gerald Wilemski
- Department of Physics, Missouri University of Science and Technology, Rolla, Missouri 65409, USA
| |
Collapse
|
2
|
Burrows SA, Shon JW, Peychev B, Slavchov RI, Smoukov SK. Phase transitions of fluorotelomer alcohols at the water|alkane interface studied via molecular dynamics simulation. SOFT MATTER 2024; 20:2243-2257. [PMID: 38351894 DOI: 10.1039/d3sm01444d] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 03/07/2024]
Abstract
Fluorosurfactants are long-lasting environmental pollutants that accumulate at interfaces ranging from aerosol droplet surfaces to cell membranes. Modeling of adsorption-based removal technologies for fluorosurfactants requires accurate simulation methods which can predict their adsorption isotherm and monolayer structure. Fluorotelomer alcohols with one or two methylene groups adjacent to the alcohol (7 : 1 FTOH and 6 : 2 FTOH, respectively) are investigated using the OPLS-AA force field at the water|hexane interface, varying the interfacial area per surfactant. The acquired interfacial pressure isotherms and monolayer phase behavior are compared with previous experimental results. The results are consistent with the experimental data inasmuch as, at realistic adsorption densities, only 7 : 1 FTOH shows a phase transition between liquid-expanded (LE) and 2D crystalline phases. Structures of the LE and crystalline phases are in good agreement with the sticky disc and Langmuir defective crystal models, respectively, used previously to interpret experimental data. Interfacial pressure of the LE phase agrees well with experiment, and sticky disc interaction parameters indicate no 2D LE-gas transition is present for either molecule. Conformation analysis reveals 7 : 1 FTOH favors conformers where the OH dipole is perpendicular to the molecular backbone, such that the crystalline phase is stabilized when these dipoles align.
Collapse
Affiliation(s)
- Stephen A Burrows
- Centre for Sustainable Engineering, School of Engineering and Materials Science, Queen Mary University of London, Mile End Road, London, E1 4NS, UK.
| | - Jang Won Shon
- Centre for Sustainable Engineering, School of Engineering and Materials Science, Queen Mary University of London, Mile End Road, London, E1 4NS, UK.
| | - Boyan Peychev
- Centre for Sustainable Engineering, School of Engineering and Materials Science, Queen Mary University of London, Mile End Road, London, E1 4NS, UK.
| | - Radomir I Slavchov
- Centre for Sustainable Engineering, School of Engineering and Materials Science, Queen Mary University of London, Mile End Road, London, E1 4NS, UK.
| | - Stoyan K Smoukov
- Centre for Sustainable Engineering, School of Engineering and Materials Science, Queen Mary University of London, Mile End Road, London, E1 4NS, UK.
| |
Collapse
|
3
|
Wei Z, Elliott JD, Papaderakis AA, Dryfe RA, Carbone P. Relation between Double Layer Structure, Capacitance, and Surface Tension in Electrowetting of Graphene and Aqueous Electrolytes. J Am Chem Soc 2024; 146:760-772. [PMID: 38153698 PMCID: PMC10785801 DOI: 10.1021/jacs.3c10814] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/30/2023] [Revised: 12/03/2023] [Accepted: 12/05/2023] [Indexed: 12/29/2023]
Abstract
Deciphering the mechanisms of charge storage on carbon-based materials is pivotal for the development of next-generation electrochemical energy storage systems. Graphene, the building block of graphitic electrodes, is an ideal model for probing such processes on a fundamental level. Herein, we investigate the thermodynamics of the graphene/aqueous electrolyte interface by utilizing a multiscale quantum mechanics-classical molecular dynamics (QM/MD) approach to provide insights into the effect of alkali metal ion (Li+) concentration on the interfacial tension (γSL) of the charged graphene/electrolyte interface. We demonstrate that the dependence of γSL on the applied surface charge exhibits an asymmetric behavior relative to the neutral surface. At the positively charged graphene sheet, the electrowetting response is amplified by electrolyte concentration, resulting in a strongly hydrophilic surface. On the contrary, at negative potential bias, γSL shows a weaker response to the charging of the electrode. Changes in γSL greatly affect the total areal capacitance predicted by the Young-Lippmann equation but have a negligible impact on the simulated total areal capacitance, indicating that the EDL structure is not directly correlated with the wettability of the surface and different interfacial mechanisms drive the two phenomena. The proposed model is validated experimentally by studying the electrowetting response of highly oriented pyrolytic graphite over a wide range of electrolyte concentrations. Our work presents the first combined theoretical and experimental study on electrowetting using carbon surfaces, introducing new conceptual routes for the investigation of wetting phenomena under potential bias.
Collapse
Affiliation(s)
- Zixuan Wei
- Department
of Chemical Engineering, The University
of Manchester, Oxford Road, Manchester M13 9PL, United Kingdom
| | - Joshua D. Elliott
- Diamond
Light Source, Diamond House, Harwell Science
and Innovation Park, Oxfordshire, Didcot OX11 ODE, United Kingdom
| | - Athanasios A. Papaderakis
- Department
of Chemistry and Henry Royce Institute, The University of Manchester, Oxford Road, Manchester M13 9PL, United
Kingdom
| | - Robert A.W. Dryfe
- Department
of Chemistry and Henry Royce Institute, The University of Manchester, Oxford Road, Manchester M13 9PL, United
Kingdom
| | - Paola Carbone
- Department
of Chemical Engineering, The University
of Manchester, Oxford Road, Manchester M13 9PL, United Kingdom
| |
Collapse
|
4
|
Le Donne A, Littlefair JD, Tortora M, Merchiori S, Bartolomé L, Grosu Y, Meloni S. Hydrophobicity of molecular-scale textured surfaces: The case of zeolitic imidazolate frameworks, an atomistic perspective. J Chem Phys 2023; 159:184709. [PMID: 37955326 DOI: 10.1063/5.0173110] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/20/2023] [Accepted: 10/19/2023] [Indexed: 11/14/2023] Open
Abstract
Hydrophobicity has proven fundamental in an inexhaustible amount of everyday applications. Material hydrophobicity is determined by chemical composition and geometrical characteristics of its macroscopic surface. Surface roughness or texturing enhances intrinsic hydrophilic or hydrophobic characteristics of a material. Here we consider crystalline surfaces presenting molecular-scale texturing typical of crystalline porous materials, e.g., metal-organic frameworks. In particular, we investigate one such material with remarkable hydrophobic qualities, ZIF-8. We show that ZIF-8 hydrophobicity is driven not only by its chemical composition but also its sub-nanoscale surface corrugations, a physical enhancement rare amongst hydrophobes. Studying ZIF-8's hydrophobic properties is challenging as experimentally it is difficult to distinguish between the materials' and the macroscopic corrugations' contributions to the hydrophobicity. The computational contact angle determination is also difficult as the standard "geometric" technique of liquid nanodroplet deposition is prone to many artifacts. Here, we characterise ZIF-8 hydrophobicity via: (i) the "geometric" approach and (ii) the "energetic" method, utilising the Young-Dupré formula and computationally determining the liquid-solid adhesion energy. Both approaches reveal nanoscale Wenzel-like bathing of the corrugated surface. Moreover, we illustrate the importance of surface linker termination in ZIF-8 hydrophobicity, which reduces when varied from sp3 N to sp2 N termination. We also consider halogenated analogues of the methyl-imidazole linker, which promote the transition from nanoWenzel-like to nanoCassie-Baxter-like states, further enhancing surface hydrophobicity. Present results reveal the complex interface physics and chemistry between water and complex porous, molecular crystalline surfaces, providing a hint to tune their hydrophobicity.
Collapse
Affiliation(s)
- Andrea Le Donne
- Dipartimento di Scienze Chimiche, Farmaceutiche ed Agrarie (DOCPAS), Università degli Studi di Ferrara (Unife), Via Luigi Borsari 46, I-44121 Ferrara, Italy
| | - Josh D Littlefair
- Dipartimento di Scienze Chimiche, Farmaceutiche ed Agrarie (DOCPAS), Università degli Studi di Ferrara (Unife), Via Luigi Borsari 46, I-44121 Ferrara, Italy
| | - Marco Tortora
- Dipartimento di Ingegneria Meccanica e Aerospaziale, Universitá di Roma "Sapienza," Via Eudossiana 18, 00184 Rome, Italy
| | - Sebastiano Merchiori
- Dipartimento di Scienze Chimiche, Farmaceutiche ed Agrarie (DOCPAS), Università degli Studi di Ferrara (Unife), Via Luigi Borsari 46, I-44121 Ferrara, Italy
| | - Luis Bartolomé
- Centre for Cooperative Research on Alternative Energies (CIC energiGUNE), Basque Research and Technology Alliance (BRTA), Alava Technology Park, Albert Einstein 48, 01510 Vitoria-Gasteiz, Spain
| | - Yaroslav Grosu
- Centre for Cooperative Research on Alternative Energies (CIC energiGUNE), Basque Research and Technology Alliance (BRTA), Alava Technology Park, Albert Einstein 48, 01510 Vitoria-Gasteiz, Spain
- Institute of Chemistry, University of Silesia in Katowice, Szkolna 9, 40-006 Katowice, Poland
| | - Simone Meloni
- Dipartimento di Scienze Chimiche, Farmaceutiche ed Agrarie (DOCPAS), Università degli Studi di Ferrara (Unife), Via Luigi Borsari 46, I-44121 Ferrara, Italy
| |
Collapse
|
5
|
R Leivas F, Barbosa MC. Functionalized carbon nanocones performance in water harvesting. J Chem Phys 2023; 158:2890471. [PMID: 37184010 DOI: 10.1063/5.0142718] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/16/2023] [Accepted: 04/30/2023] [Indexed: 05/16/2023] Open
Abstract
In this work, we investigate the water capture process for functionalized carbon nanocones (CNCs) through molecular dynamic simulations in the following three scenarios: a single CNC in contact with a reservoir containing liquid water, a single CNC in contact with a water vapor reservoir, and a combination of more than one CNC in contact with vapor. We found that water flows through the nanocones when in contact with the liquid reservoir if the nanocone tip presents hydrophilic functionalization. In contact with steam, we observed the formation of droplets at the base of the nanocone only when hydrophilic functionalization is present. Then, water flows through in a linear manner, a process that is more efficient than that in the liquid reservoir regime. The scalability of the process is tested by analyzing the water flow through more than one nanocone. The results suggest that the distance between the nanocones is a fundamental ingredient for the efficiency of water harvesting.
Collapse
Affiliation(s)
- Fernanda R Leivas
- Instituto de Física, Universidade Federal do Rio Grande do Sul, CP 15051, 91501-970 Porto Alegre, RS, Brazil
| | - Marcia C Barbosa
- Instituto de Física, Universidade Federal do Rio Grande do Sul, CP 15051, 91501-970 Porto Alegre, RS, Brazil
| |
Collapse
|
6
|
Sanchez-Burgos I, Muniz MC, Espinosa JR, Panagiotopoulos AZ. A Deep Potential model for liquid-vapor equilibrium and cavitation rates of water. J Chem Phys 2023; 158:2889532. [PMID: 37158636 DOI: 10.1063/5.0144500] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/30/2023] [Accepted: 03/30/2023] [Indexed: 05/10/2023] Open
Abstract
Computational studies of liquid water and its phase transition into vapor have traditionally been performed using classical water models. Here, we utilize the Deep Potential methodology-a machine learning approach-to study this ubiquitous phase transition, starting from the phase diagram in the liquid-vapor coexistence regime. The machine learning model is trained on ab initio energies and forces based on the SCAN density functional, which has been previously shown to reproduce solid phases and other properties of water. Here, we compute the surface tension, saturation pressure, and enthalpy of vaporization for a range of temperatures spanning from 300 to 600 K and evaluate the Deep Potential model performance against experimental results and the semiempirical TIP4P/2005 classical model. Moreover, by employing the seeding technique, we evaluate the free energy barrier and nucleation rate at negative pressures for the isotherm of 296.4 K. We find that the nucleation rates obtained from the Deep Potential model deviate from those computed for the TIP4P/2005 water model due to an underestimation in the surface tension from the Deep Potential model. From analysis of the seeding simulations, we also evaluate the Tolman length for the Deep Potential water model, which is (0.091 ± 0.008) nm at 296.4 K. Finally, we identify that water molecules display a preferential orientation in the liquid-vapor interface, in which H atoms tend to point toward the vapor phase to maximize the enthalpic gain of interfacial molecules. We find that this behavior is more pronounced for planar interfaces than for the curved interfaces in bubbles. This work represents the first application of Deep Potential models to the study of liquid-vapor coexistence and water cavitation.
Collapse
Affiliation(s)
- Ignacio Sanchez-Burgos
- Maxwell Centre, Cavendish Laboratory, Department of Physics, University of Cambridge, J J Thomson Avenue,Cambridge CB3 0HE, United Kingdom
- Department of Chemical and Biological Engineering, Princeton University, Princeton, New Jersey 08544, USA
| | - Maria Carolina Muniz
- Department of Chemical and Biological Engineering, Princeton University, Princeton, New Jersey 08544, USA
| | - Jorge R Espinosa
- Maxwell Centre, Cavendish Laboratory, Department of Physics, University of Cambridge, J J Thomson Avenue,Cambridge CB3 0HE, United Kingdom
- Departamento de Química Fisica, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | | |
Collapse
|
7
|
Johnson LJW, Paulo G, Bartolomé L, Amayuelas E, Gubbiotti A, Mirani D, Le Donne A, López GA, Grancini G, Zajdel P, Meloni S, Giacomello A, Grosu Y. Optimization of the wetting-drying characteristics of hydrophobic metal organic frameworks via crystallite size: The role of hydrogen bonding between intruded and bulk liquid. J Colloid Interface Sci 2023; 645:775-783. [PMID: 37172487 DOI: 10.1016/j.jcis.2023.04.059] [Citation(s) in RCA: 7] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/16/2023] [Revised: 03/31/2023] [Accepted: 04/13/2023] [Indexed: 05/15/2023]
Abstract
HYPOTHESIS The behavior of Heterogeneous Lyophobic Systems (HLSs) comprised of a lyophobic porous material and a corresponding non-wetting liquid is affected by a variety of different structural parameters of the porous material. Dependence on exogenic properties such as crystallite size is desirable for system tuning as they are much more facilely modified. We explore the dependence of intrusion pressure and intruded volume on crystallite size, testing the hypothesis that the connection between internal cavities and bulk water facilitates intrusion via hydrogen bonding, a phenomenon that is magnified in smaller crystallites with a larger surface/volume ratio. EXPERIMENTS Water intrusion/extrusion pressures and intrusion volume were experimentally measured for ZIF-8 samples of various crystallite sizes and compared to previously reported values. Alongside the practical research, molecular dynamics simulations and stochastic modeling were performed to illustrate the effect of crystallite size on the properties of the HLSs and uncover the important role of hydrogen bonding within this phenomenon. FINDINGS A reduction in crystallite size led to a significant decrease of intrusion and extrusion pressures below 100 nm. Simulations indicate that this behavior is due to a greater number of cages being in proximity to bulk water for smaller crystallites, allowing cross-cage hydrogen bonds to stabilize the intruded state and lower the threshold pressure of intrusion and extrusion. This is accompanied by a reduction in the overall intruded volume. Simulations demonstrate that this phenomenon is linked to ZIF-8 surface half-cages exposed to water being occupied by water due to non-trivial termination of the crystallites, even at atmospheric pressure.
Collapse
Affiliation(s)
- Liam J W Johnson
- Centre for Cooperative Research on Alternative Energies (CIC energiGUNE), Basque Research and Technology Alliance (BRTA), Calle Albert Einstein, 48, Vitoria-Gasteiz, 01510, Araba/Alava, Spain; Department of Physics, Faculty of Science and Technology, University of the Basque Country (UPV/EHU), Barrio Sarriena s/n, Bilbao, 48490, Leioa, Spain
| | - Gonçalo Paulo
- Dipartimento di Ingegneria Meccanica e Aerospaziale, Sapienza Università di Roma, Roma, Italy
| | - Luis Bartolomé
- Centre for Cooperative Research on Alternative Energies (CIC energiGUNE), Basque Research and Technology Alliance (BRTA), Calle Albert Einstein, 48, Vitoria-Gasteiz, 01510, Araba/Alava, Spain
| | - Eder Amayuelas
- Centre for Cooperative Research on Alternative Energies (CIC energiGUNE), Basque Research and Technology Alliance (BRTA), Calle Albert Einstein, 48, Vitoria-Gasteiz, 01510, Araba/Alava, Spain
| | - Alberto Gubbiotti
- Dipartimento di Ingegneria Meccanica e Aerospaziale, Sapienza Università di Roma, Roma, Italy
| | - Diego Mirani
- Department of Chemistry & INSTM University of Pavia, Via Taramelli 14, Pavia, I-27100, Italy
| | - Andrea Le Donne
- Dipartimento di Scienze Chimiche e Farmaceutiche (DipSCF), Università degli Studi di Ferrara (Unife) Via Luigi Borsari 46, Ferrara, I-44121, Italy
| | - Gabriel A López
- Department of Physics, Faculty of Science and Technology, University of the Basque Country (UPV/EHU), Barrio Sarriena s/n, Bilbao, 48490, Leioa, Spain
| | - Giulia Grancini
- Department of Chemistry & INSTM University of Pavia, Via Taramelli 14, Pavia, I-27100, Italy
| | - Paweł Zajdel
- Institute of Physics, University of Silesia, 75 Pulku Piechoty 1, Chorzow, 41-500, Poland
| | - Simone Meloni
- Dipartimento di Scienze Chimiche e Farmaceutiche (DipSCF), Università degli Studi di Ferrara (Unife) Via Luigi Borsari 46, Ferrara, I-44121, Italy.
| | - Alberto Giacomello
- Dipartimento di Ingegneria Meccanica e Aerospaziale, Sapienza Università di Roma, Roma, Italy.
| | - Yaroslav Grosu
- Centre for Cooperative Research on Alternative Energies (CIC energiGUNE), Basque Research and Technology Alliance (BRTA), Calle Albert Einstein, 48, Vitoria-Gasteiz, 01510, Araba/Alava, Spain; Institute of Chemistry, University of Silesia, Szkolna 9, Katowice, 40-006, Poland.
| |
Collapse
|
8
|
Takayama T, Otosu T, Yamaguchi S. Transferability of vibrational spectroscopic map from TIP4P to TIP4P-like water models. J Chem Phys 2023; 158:136101. [PMID: 37031108 DOI: 10.1063/5.0146084] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/05/2023] Open
Abstract
We computed the IR, Raman, and sum frequency generation spectra of water in the OH-stretch region by employing the quantum/classical mixed approach that consists of a vibrational spectroscopic map and molecular dynamics (MD) simulation. We carried out the MD simulation with the TIP4P, TIP4P/2005, and TIP4P/Ice models and applied the map designed for TIP4P by Skinner et al. to each MD trajectory. Although the map is not tuned for TIP4P-like models, TIP4P/2005 and TIP4P/Ice provide the best reproduction of the experimental vibrational spectra of liquid water and crystalline ice, respectively. This result demonstrates the transferability of the map from TIP4P to TIP4P/2005 and TIP4P/Ice, meaning that one can choose an appropriate TIP4P-like model to calculate the vibrational spectra of an aqueous system without rebuilding the map.
Collapse
Affiliation(s)
- Tetsuyuki Takayama
- Department of Applied Chemistry, Graduate School of Science and Engineering, Saitama University, 255 Shimo-Okubo, Sakura, Saitama 338-8570, Japan
| | - Takuhiro Otosu
- Department of Applied Chemistry, Graduate School of Science and Engineering, Saitama University, 255 Shimo-Okubo, Sakura, Saitama 338-8570, Japan
| | - Shoichi Yamaguchi
- Department of Applied Chemistry, Graduate School of Science and Engineering, Saitama University, 255 Shimo-Okubo, Sakura, Saitama 338-8570, Japan
| |
Collapse
|
9
|
P Lamas C, Vega C, G Noya E, Sanz E. The water cavitation line as predicted by the TIP4P/2005 model. J Chem Phys 2023; 158:124504. [PMID: 37003766 DOI: 10.1063/5.0139470] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/03/2023] Open
Abstract
The formation of vapor bubbles in a metastable liquid, cavitation, is an activated process due to the free energy cost of having both phases at contact. Such an energetic penalty enables the existence of the liquid beyond its thermodynamic borders. Establishing the stability limits of a liquid as ubiquitous as water has important practical implications and has thereby attracted a lot of attention. Different experimental strategies and theoretical analyses have been employed to measure and predict the cavitation line, or the pressure-temperature kinetic stability border of liquid water. Understanding the location of the cavitation line requires knowing the cavitation rate dependence on pressure and temperature. Such dependency is difficult to obtain in experiments, and we use molecular simulations with the TIP4P/2005 model to fill this gap. By deeply overstretching liquid water below the saturation pressure, we are able to observe and quantify spontaneous cavitation. To deal with a lower overstretching regime, we resort to the Seeding technique, which consists of analyzing simulations of a liquid containing a vapor bubble under the theoretical framework of Classical Nucleation Theory. Combining spontaneous cavitation with Seeding, we get a wide overview of the cavitation rate. We study two different temperatures (450 and 550 K) and complement our perspective with the results previously obtained at 296.4 K [Menzl et al., Proc. Natl. Acad. Sci. 113, 13582 (2016)] to establish a broad simulation-experiment comparison. We find a good agreement between simulations and both isobaric heating and isochoric cooling experiments using quartz inclusions. We are, however, unable to reconcile simulations with other experimental techniques. Our results predict a decrease in the solid-liquid interfacial free energy as the liquid becomes increasingly overstretched with a temperature independent Tolman length of 0.1 nm. Therefore, the capillarity approximation underestimates the nucleation rate. Nonetheless, it provides a fair indication of the location of the cavitation line given the steep rate vs pressure dependence. Overall, our work provides a comprehensive view of the water cavitation phenomenon and sets an efficient strategy to investigate it with molecular simulations.
Collapse
Affiliation(s)
- Cintia P Lamas
- Departamento de Química-Física I (Unidad de I+D+i Asociada al CSIC), Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - Carlos Vega
- Departamento de Química-Física I (Unidad de I+D+i Asociada al CSIC), Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - Eva G Noya
- Instituto de Química Física Rocasolano, CSIC, C/ Serrano 119, 28006 Madrid, Spain
| | - Eduardo Sanz
- Departamento de Química-Física I (Unidad de I+D+i Asociada al CSIC), Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| |
Collapse
|
10
|
Hrahsheh F, Wilemski G. Effects of molecular size and orientation on the interfacial properties and wetting behavior of water/ n-alkane systems: a molecular-dynamics study. Phys Chem Chem Phys 2023; 25:5808-5816. [PMID: 36744733 DOI: 10.1039/d2cp05735b] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/28/2023]
Abstract
Molecular dynamics simulations (MD) are performed to study the interfacial structure/tension and wetting behavior of water/n-alkane systems (water/nC5 to water/nC16 where nCx = CxH(2x + 2)). In particular, we study complete-to-partial wetting transitions by changing the n-alkane chain length (NC) at a constant temperature, T = 295 K. Simulations are carried out with a united-atom TraPPE model for n-alkanes and the TIP4P-2005 model of water. Simulation results are in excellent agreement with the initial spreading coefficients and contact angles calculated using experimental values of the surface and interfacial tensions. In addition, it has been determined that water/(nC5-nC7) and water/(nC8-nC16), respectively, exhibit complete and partial initial wetting modes. Simulations show that the interfacial structures of water/(nC5-nC7) are different from water/(nC8-nC16) systems. In the latter, water preferentially orients near the interface to increase the number of hydrogen bonds and the charge and mass densities. Moreover, the orientation of n-alkane molecules at water/(nC8-nC16) interfaces has a long-range persistence, resulting in layered structures that increase with NC. In addition, simulation results of the orientational order parameter Sz show alignment behavior of the n-alkane molecules with respect to the interfaces. Simulations predict that the central segments of n-alkane are strongly packed in the interfaces while the end segments (methyl groups) form smaller peaks in the outer edge of the layer. This observation confirms the "horseshoe" or "C-shaped" structure of n-alkane molecules in the water/n-alkane interfaces. At constant temperature, the interface widths of both water and the n-alkanes decrease with increasing n-alkane molecular length. These results suggest that increasing the n-alkane chain length affects the water/n-alkane interfacial properties in a manner similar to that of cooling.
Collapse
Affiliation(s)
- Fawaz Hrahsheh
- Higher Colleges of Technology, ETS, MZWC, Abu Dhabi, 58855, United Arab Emirates.
| | - Gerald Wilemski
- Department of Physics, Missouri University of Science and Technology, Rolla, MO 65409, USA.
| |
Collapse
|
11
|
Blazquez S, Conde MM, Vega C. Scaled charges for ions: An improvement but not the final word for modeling electrolytes in water. J Chem Phys 2023; 158:054505. [PMID: 36754806 DOI: 10.1063/5.0136498] [Citation(s) in RCA: 15] [Impact Index Per Article: 15.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/09/2023] Open
Abstract
In this work, we discuss the use of scaled charges when developing force fields for NaCl in water. We shall develop force fields for Na+ and Cl- using the following values for the scaled charge (in electron units): ±0.75, ±0.80, ±0.85, and ±0.92 along with the TIP4P/2005 model of water (for which previous force fields were proposed for q = ±0.85 and q = ±1). The properties considered in this work are densities, structural properties, transport properties, surface tension, freezing point depression, and maximum in density. All the developed models were able to describe quite well the experimental values of the densities. Structural properties were well described by models with charges equal to or larger than ±0.85, surface tension by the charge ±0.92, maximum in density by the charge ±0.85, and transport properties by the charge ±0.75. The use of a scaled charge of ±0.75 is able to reproduce with high accuracy the viscosities and diffusion coefficients of NaCl solutions for the first time. We have also considered the case of KCl in water, and the results obtained were fully consistent with those of NaCl. There is no value of the scaled charge able to reproduce all the properties considered in this work. Although certainly scaled charges are not the final word in the development of force fields for electrolytes in water, its use may have some practical advantages. Certain values of the scaled charge could be the best option when the interest is to describe certain experimental properties.
Collapse
Affiliation(s)
- S Blazquez
- Dpto. Química Física I, Fac. Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - M M Conde
- Departamento de Ingeniería Química Industrial y Medio Ambiente, Escuela Técnica Superior de Ingenieros Industriales, Universidad Politécnica de Madrid, 28006 Madrid, Spain
| | - C Vega
- Dpto. Química Física I, Fac. Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| |
Collapse
|
12
|
Leivas FR, Barbosa MC. Atmospheric water harvesting using functionalized carbon nanocones. BEILSTEIN JOURNAL OF NANOTECHNOLOGY 2023; 14:1-10. [PMID: 36703909 PMCID: PMC9830493 DOI: 10.3762/bjnano.14.1] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/27/2022] [Accepted: 12/14/2022] [Indexed: 05/28/2023]
Abstract
In this work, we propose a method to harvest liquid water from water vapor using carbon nanocones. The condensation occurs due to the presence of hydrophilic sites at the nanocone entrance. The functionalization, together with the high mobility of water inside nanostructures, leads to a fast water flow through the nanostructure. We show using molecular dynamics simulations that this device is able to collect water if the surface functionalization is properly selected.
Collapse
Affiliation(s)
- Fernanda R Leivas
- Instituto de Física, Universidade Federal do Rio Grande do Sul, CP 15051, 91501-970, Porto Alegre, RS, Brazil
| | - Marcia C Barbosa
- Instituto de Física, Universidade Federal do Rio Grande do Sul, CP 15051, 91501-970, Porto Alegre, RS, Brazil
| |
Collapse
|
13
|
The lung surfactant activity probed with molecular dynamics simulations. Adv Colloid Interface Sci 2022; 304:102659. [PMID: 35421637 DOI: 10.1016/j.cis.2022.102659] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/20/2021] [Revised: 03/18/2022] [Accepted: 03/31/2022] [Indexed: 01/17/2023]
Abstract
The surface of pulmonary alveolar subphase is covered with a mixture of lipids and proteins. This lung surfactant plays a crucial role in lung functioning. It shows a complex phase behavior which can be altered by the interaction with third molecules such as drugs or pollutants. For studying multicomponent biological systems, it is of interest to couple experimental approach with computational modelling yielding atomic-scale information. Simple two, three, or four-component model systems showed to be useful for getting more insight in the interaction between lipids, lipids and proteins or lipids and proteins with drugs and impurities. These systems were studied theoretically using molecular dynamic simulations and experimentally by means of the Langmuir technique. A better understanding of the structure and behavior of lung surfactants obtained from this research is relevant for developing new synthetic surfactants for efficient therapies, and may contribute to public health protection.
Collapse
|
14
|
Atherton D, Michaelides A, Cox SJ. Can molecular simulations reliably compare homogeneous and heterogeneous ice nucleation? J Chem Phys 2022; 156:164501. [PMID: 35490004 DOI: 10.1063/5.0085750] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
In principle, the answer to the posed titular question is undoubtedly "yes." But in practice, requisite reference data for homogeneous systems have been obtained with a treatment of intermolecular interactions that is different from that typically employed for heterogeneous systems. In this article, we assess the impact of the choice of truncation scheme when comparing water in homogeneous and inhomogeneous environments. Specifically, we use explicit free energy calculations and a simple mean field analysis to demonstrate that using the "cut-and-shift" version of the Lennard-Jones potential (common to most simple point charge models of water) results in a systematic increase in the melting temperature of ice Ih. In addition, by drawing an analogy between a change in cutoff and a change in pressure, we use existing literature data for homogeneous ice nucleation at negative pressures to suggest that enhancements due to heterogeneous nucleation may have been overestimated by several orders of magnitude.
Collapse
Affiliation(s)
- Dominic Atherton
- Yusuf Hamied Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| | - Angelos Michaelides
- Yusuf Hamied Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| | - Stephen J Cox
- Yusuf Hamied Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| |
Collapse
|
15
|
Bhakat S, Söderhjelm P. Flap Dynamics in Pepsin-Like Aspartic Proteases: A Computational Perspective Using Plasmepsin-II and BACE-1 as Model Systems. J Chem Inf Model 2022; 62:914-926. [PMID: 35138093 PMCID: PMC8889585 DOI: 10.1021/acs.jcim.1c00840] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/13/2021] [Indexed: 11/28/2022]
Abstract
The flexibility of β hairpin structure known as the flap plays a key role in catalytic activity and substrate intake in pepsin-like aspartic proteases. Most of these enzymes share structural and sequential similarity. In this study, we have used apo Plm-II and BACE-1 as model systems. In the apo form of the proteases, a conserved tyrosine residue in the flap region remains in a dynamic equilibrium between the normal and flipped states through rotation of the χ1 and χ2 angles. Independent MD simulations of Plm-II and BACE-1 remained stuck either in the normal or flipped state. Metadynamics simulations using side-chain torsion angles (χ1 and χ2 of tyrosine) as collective variables sampled the transition between the normal and flipped states. Qualitatively, the two states were predicted to be equally populated. The normal and flipped states were stabilized by H-bond interactions to a tryptophan residue and to the catalytic aspartate, respectively. Further, mutation of tyrosine to an amino-acid with smaller side-chain, such as alanine, reduced the flexibility of the flap and resulted in a flap collapse (flap loses flexibility and remains stuck in a particular state). This is in accordance with previous experimental studies, which showed that mutation to alanine resulted in loss of activity in pepsin-like aspartic proteases. Our results suggest that the ring flipping associated with the tyrosine side-chain is the key order parameter that governs flap dynamics and opening of the binding pocket in most pepsin-like aspartic proteases.
Collapse
Affiliation(s)
- Soumendranath Bhakat
- Division
of Biophysical Chemistry, Center for Molecular Protein Science, Department
of Chemistry, Lund University, P.O. Box 124, SE-22100 Lund, Sweden
- Department
of Biochemistry and Molecular Biophysics, Washington University, School of Medicine, St. Louis, Missouri 63110, United States
| | - Pär Söderhjelm
- Division
of Biophysical Chemistry, Center for Molecular Protein Science, Department
of Chemistry, Lund University, P.O. Box 124, SE-22100 Lund, Sweden
| |
Collapse
|
16
|
Papadopoulou E, Megaridis CM, Walther JH, Koumoutsakos P. Nanopumps without Pressure Gradients: Ultrafast Transport of Water in Patterned Nanotubes. J Phys Chem B 2022; 126:660-669. [DOI: 10.1021/acs.jpcb.1c07562] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Ermioni Papadopoulou
- Computational Science and Engineering Laboratory, ETH Zürich, Zürich CH-8092, Switzerland
| | - Constantine M. Megaridis
- Department of Mechanical and Industrial Engineering, University of Illinois at Chicago, Chicago, Illinois 60607, United States
| | - Jens H. Walther
- Computational Science and Engineering Laboratory, ETH Zürich, Zürich CH-8092, Switzerland
- Department of Mechanical Engineering, Technical University of Denmark, 2800 Kongens Lyngby, Denmark
| | - Petros Koumoutsakos
- Computational Science and Engineering Laboratory, ETH Zürich, Zürich CH-8092, Switzerland
- John A. Paulson School of Engineering and Applied Sciences, Harvard University, Cambridge, Massachusetts 02138, United States
| |
Collapse
|
17
|
Gao Z, Wu W, Sun W, Wang B. Understanding the Stabilization of a Bulk Nanobubble: A Molecular Dynamics Analysis. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2021; 37:11281-11291. [PMID: 34520212 DOI: 10.1021/acs.langmuir.1c01796] [Citation(s) in RCA: 13] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/13/2023]
Abstract
Bulk nanobubbles (NBs) have received considerable attention because of their extensive potential applications, such as in ultrasound imaging and water management. Although multiple types of experimental evidence have supported the existence and stabilization of bulk NBs, the underlying mechanism remains unclear. This study numerically investigates the bulk NB stabilization with molecular dynamics (MD) methods: the all-atom (AA) MD simulation is used for NBs of several nanometers diameter; the coarse-grained (CG) MD simulation is for the NBs of about 100 nm. The NB properties are statistically obtained and analyzed, including the inner density, inner pressure, surface charge, interfacial hydrogen bond (HB), and gaseous diffusion. The results show that the gas inside an NB has ultrahigh density (tens of kilograms per cubic meter). A double-layer surface charge exists on the NB. The inner/outer layer is positively/negatively charged, and the electrostatic stress can counteract part of the surface tension. In addition, the interfacial HB is weakened by the interaction between gas and water molecules, causing less surface tension. The above features are beneficial to NB stabilization. The NB equilibrium radii solved by the interfacial mechanical equilibrium equation agree with the MD results, indicating that this equation can describe the force balance of an NB as small as several nanometers. Besides, supersaturation appears to be necessary for the NB thermodynamic equilibrium. Based on Henry's law and the ideal gas law, the theoretical analysis suggests that the stability of the NB thermodynamic equilibrium is conditional: the number of gas molecules in NBs should be more than half that dissolved in liquid. This study unravels a stabilized bulk NB's properties and discusses the NB equilibrium and stabilization mechanism, which will advance the understanding and application of bulk NBs.
Collapse
Affiliation(s)
- Zhan Gao
- School of Aerospace Engineering, Tsinghua University, Beijing 100084, People's Republic of China
| | - Wangxia Wu
- School of Aerospace Engineering, Beijing Institute of Technology, Beijing 100081, People's Republic of China
| | - Weitao Sun
- School of Aerospace Engineering, Tsinghua University, Beijing 100084, People's Republic of China
| | - Bing Wang
- School of Aerospace Engineering, Tsinghua University, Beijing 100084, People's Republic of China
| |
Collapse
|
18
|
Gittus OR, Bresme F. Thermophysical properties of water using reactive force fields. J Chem Phys 2021; 155:114501. [PMID: 34551553 DOI: 10.1063/5.0057868] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
The widescale importance and rich phenomenology of water continue to motivate the development of computational models. ReaxFF force fields incorporate many characteristics desirable for modeling aqueous systems: molecular flexibility, polarization, and chemical reactivity (bond formation and breaking). However, their ability to model the general properties of water has not been evaluated in detail. We present comprehensive benchmarks of the thermophysical properties of water for two ReaxFF models, the water-2017 and CHON-2017_weak force fields. These include structural, electrostatic, vibrational, thermodynamic, coexistence, and transport properties at ambient conditions (300 K and 0.997 g cm-3) and along the standard pressure (1 bar) isobar. Overall, CHON-2017_weak predicts more accurate thermophysical properties than the water-2017 force field. Based on our results, we recommend potential avenues for improvement: the dipole moment to quadrupole moment ratio, the self-diffusion coefficient, especially for water-2017, and the gas phase vibrational frequencies with the aim to improve the vibrational properties of liquid water.
Collapse
Affiliation(s)
- Oliver R Gittus
- Department of Chemistry, Molecular Sciences Research Hub Imperial College, London W12 0BZ, United Kingdom
| | - Fernando Bresme
- Department of Chemistry, Molecular Sciences Research Hub Imperial College, London W12 0BZ, United Kingdom
| |
Collapse
|
19
|
Naeiji P, Woo TK, Alavi S, Ohmura R. Molecular dynamics simulations of interfacial properties of the CO 2-water and CO 2-CH 4-water systems. J Chem Phys 2020; 153:044701. [PMID: 32752701 DOI: 10.1063/5.0008114] [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/29/2022] Open
Abstract
Molecular dynamics simulations were performed to study the interfacial behavior of the pure carbon dioxide-water system and a binary 40:60 mol. % gas mixture of (carbon dioxide + methane)-water at the temperatures of 275.15 K and 298.15 K and pressures near 4 MPa for CO2 and up to 10 MPa for methane. The simulations are used to study the dynamic equilibrium of the gases at the water-gas interface, to determine the z-density profiles for the gases and water, and calculate the interfacial tension γ under the different temperature/pressure conditions close to those of the formation of clathrate hydrates of these gases. At the same hydrostatic gas phase pressure, the CO2-water interface has a lower interfacial tension than the CH4-water interface. A greater number of CO2 molecules, as much as three times more than methane at the same pressure, were adsorbed at the interfacial layer, which reflects the stronger electrostatic quadrupolar and van der Waals interactions between CO2 and water molecules at the interface. The water surfaces are covered by less than a monolayer of gas even when the pressure of the system goes near the saturation pressure of CO2. The surface adsorbed molecules are in dynamic equilibrium with the bulk gas and with exchange between the gas and interface regions occurring repeatedly within the timescale of the simulations. The effects of the changes in the CO2-water interfacial tension with external temperature and pressure conditions on the formation of the clathrate hydrates and other CO2 capture and sequestration processes are discussed.
Collapse
Affiliation(s)
- Parisa Naeiji
- Department of Chemistry and Biomolecular Sciences, University of Ottawa, Ottawa, Ontario K1A 0R6, Canada
| | - Tom K Woo
- Department of Chemistry and Biomolecular Sciences, University of Ottawa, Ottawa, Ontario K1A 0R6, Canada
| | - Saman Alavi
- Department of Chemistry and Biomolecular Sciences, University of Ottawa, Ottawa, Ontario K1A 0R6, Canada
| | - Ryo Ohmura
- Department of Mechanical Engineering, Keio University, 3-14-1 Hiyoshi, Kohoku-Ku, Yokohama 223-8522, Japan
| |
Collapse
|
20
|
Le Breton G, Joly L. Molecular modeling of aqueous electrolytes at interfaces: Effects of long-range dispersion forces and of ionic charge rescaling. J Chem Phys 2020; 152:241102. [PMID: 32610967 DOI: 10.1063/5.0011058] [Citation(s) in RCA: 17] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2023] Open
Abstract
Molecular dynamics simulations of aqueous electrolytes generally rely on empirical force fields, combining dispersion interactions-described by a truncated Lennard-Jones (LJ) potential-and electrostatic interactions-described by a Coulomb potential computed with a long-range solver. Recently, force fields using rescaled ionic charges [electronic continuum correction (ECC)], possibly complemented with rescaling of LJ parameters [ECC rescaled (ECCR)], have shown promising results in bulk, but their performance at interfaces has been less explored. Here, we started by exploring the impact of the LJ potential truncation on the surface tension of a sodium chloride aqueous solution. We show a discrepancy between the numerical predictions for truncated LJ interactions with a large cutoff and for untruncated LJ interactions computed with a long-range solver, which can bias comparison of force field predictions with experiments. Using a long-range solver for LJ interactions, we then show that an ionic charge rescaling factor chosen to correct long-range electrostatic interactions in bulk accurately describes image charge repulsion at the liquid-vapor interface, and the rescaling of LJ parameters in ECCR models-aimed at capturing local ion-ion and ion-water interactions in bulk- describes well the formation of an ionic double layer at the liquid-vapor interface. Overall, these results suggest that the molecular modeling of aqueous electrolytes at interfaces would benefit from using long-range solvers for dispersion forces and from using ECCR models, where the charge rescaling factor should be chosen to correct long-range electrostatic interactions.
Collapse
Affiliation(s)
- Guillaume Le Breton
- Département de Physique, École Normale Supérieure de Lyon, 46 Allée d'Italie, Lyon Cedex 07, France
| | - Laurent Joly
- Univ. Lyon, Univ. Claude Bernard Lyon 1, CNRS, Institut Lumière Matière, F-69622 Villeurbanne, France
| |
Collapse
|
21
|
Gravelle S, Kamal C, Botto L. Liquid exfoliation of multilayer graphene in sheared solvents: A molecular dynamics investigation. J Chem Phys 2020; 152:104701. [PMID: 32171224 DOI: 10.1063/1.5141515] [Citation(s) in RCA: 16] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Liquid-phase exfoliation, the use of a sheared liquid to delaminate graphite into few-layer graphene, is a promising technique for the large-scale production of graphene. However, the microscale and nanoscale fluid-structure processes controlling the exfoliation are not fully understood. Here, we perform non-equilibrium molecular dynamics simulations of a defect-free graphite nanoplatelet suspended in a shear flow and measure the critical shear rate γ̇c needed for the exfoliation to occur. We compare γ̇c for different solvents, including water and N-methyl-pyrrolidone, and nanoplatelets of different lengths. Using a theoretical model based on a balance between the work done by viscous shearing forces and the change in interfacial energies upon layer sliding, we are able to predict the critical shear rates γ̇c measured in simulations. We find that an accurate prediction of the exfoliation of short graphite nanoplatelets is possible only if both hydrodynamic slip and the fluid forces on the graphene edges are considered and if an accurate value of the solid-liquid surface energy is used. The commonly used "geometric-mean" approximation for the solid-liquid energy leads to grossly incorrect predictions.
Collapse
Affiliation(s)
- Simon Gravelle
- School of Engineering and Material Science, Queen Mary University of London, London, United Kingdom
| | - Catherine Kamal
- School of Engineering and Material Science, Queen Mary University of London, London, United Kingdom
| | - Lorenzo Botto
- School of Engineering and Material Science, Queen Mary University of London, London, United Kingdom
| |
Collapse
|
22
|
Charge transfer as a ubiquitous mechanism in determining the negative charge at hydrophobic interfaces. Nat Commun 2020; 11:901. [PMID: 32060273 PMCID: PMC7021814 DOI: 10.1038/s41467-020-14659-5] [Citation(s) in RCA: 47] [Impact Index Per Article: 11.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/28/2019] [Accepted: 01/22/2020] [Indexed: 01/17/2023] Open
Abstract
The origin of the apparent negative charge at hydrophobic–water interfaces has fueled debates in the physical chemistry community for decades. The most common interpretation given to explain this observation is that negatively charged hydroxide ions (OH–) bind strongly to the interfaces. Using first principles calculations of extended air–water and oil–water interfaces, we unravel a mechanism that does not require the presence of OH–. Small amounts of charge transfer along hydrogen bonds and asymmetries in the hydrogen bond network due to topological defects can lead to the accumulation of negative surface charge at both interfaces. For water near oil, some spillage of electron density into the oil phase is also observed. The computed surface charge densities at both interfaces is approximately \documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$-0.015\ {\rm{e}}/{{\rm{nm}}}^{2}$$\end{document}−0.015e∕nm2 in agreement with electrophoretic experiments. We also show, using an energy decomposition analysis, that the electronic origin of this phenomena is rooted in a collective polarization/charge transfer effect. The accumulation of negative charge at hydrophobic–water interfaces has been a source of debate for a long time. Here the authors use ab initio calculations to show that the charge accumulation at air–water and oil–water interfaces is caused by subtle charge transfer processes.
Collapse
|
23
|
Yuan Y, Ma Z, Wang F. Leveraging local MP2 to reduce basis set superposition errors: An efficient first-principles based force-field for carbon dioxide. J Chem Phys 2019; 151:184501. [PMID: 31731863 DOI: 10.1063/1.5124811] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Pairwise additive model potentials for CO2 were developed by fitting to gradients computed with the local second order Møller Plesset Perturbation theory (LMP2) method, with and without consideration of 3-body dispersion using adaptive force matching. Without fitting to experiments, all models gave good predictions of properties of CO2, such as the density-temperature diagram, diffusion constants, and radial distribution functions. For the prediction of vibrational spectra, the inclusion of a bond-bond coupling term has been shown to be important. The CO2 models developed only have pairwise additive terms, thus allowing microsecond time scale simulations to be performed with practical computational cost. LMP2 performed significantly better than second order Møller Plesset Perturbation theory (MP2) for the development of the CO2 model. This is attributed to the appreciable reduction in the basis set superposition error when the localized method was used. It is argued that LMP2 is a more appropriate method than MP2 for force matching for systems where the basis set superposition error is large.
Collapse
Affiliation(s)
- Ying Yuan
- Department of Chemistry and Biochemistry, University of Arkansas, Fayetteville, Arkansas 72701, USA
| | - Zhonghua Ma
- Department of Chemistry and Biochemistry, University of Arkansas, Fayetteville, Arkansas 72701, USA
| | - Feng Wang
- Department of Chemistry and Biochemistry, University of Arkansas, Fayetteville, Arkansas 72701, USA
| |
Collapse
|
24
|
Qiu Y, Nerenberg PS, Head-Gordon T, Wang LP. Systematic Optimization of Water Models Using Liquid/Vapor Surface Tension Data. J Phys Chem B 2019; 123:7061-7073. [DOI: 10.1021/acs.jpcb.9b05455] [Citation(s) in RCA: 21] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Yudong Qiu
- Chemistry Department, University of California, Davis, Davis, California 95616, United States
| | - Paul S. Nerenberg
- Departments of Physics & Astronomy and Biological Sciences, California State University, Los Angeles, California 90032, United States
| | - Teresa Head-Gordon
- Pitzer Theory Center and Departments of Chemistry, Bioengineering and Chemical and Biomolecular Engineering, University of California, Berkeley, Berkeley, California 94720, United States
| | - Lee-Ping Wang
- Chemistry Department, University of California, Davis, Davis, California 95616, United States
| |
Collapse
|
25
|
Papadopoulou E, Megaridis CM, Walther JH, Koumoutsakos P. Ultrafast Propulsion of Water Nanodroplets on Patterned Graphene. ACS NANO 2019; 13:5465-5472. [PMID: 31025854 DOI: 10.1021/acsnano.9b00252] [Citation(s) in RCA: 24] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/12/2023]
Abstract
The directed transport of liquids at the nanoscale is of great importance for nanotechnology applications ranging from water filtration to the cooling of electronics and precision medicine. Here we demonstrate such unidirectional, pumpless transport of water nanodroplets on graphene sheets patterned with hydrophilic/phobic areas inspired by natural systems. We find that spatially varying patterning of the graphene surfaces can lead to water transport at ultrafast velocities, far exceeding macroscale estimates. We perform extensive molecular dynamics simulations to show that such high transport velocities ( O(102 m/s)) are due to differences of the advancing and receding contact angles of the moving droplet. This contact angle hysteresis and the ensuing transport depend on the surface pattern and the droplet size. We present a scaling law for the driving capillary and resisting friction forces on the water droplet and use it to predict nanodroplet trajectories on a wedge-patterned graphene sheet. The present results demonstrate that graphene with spatially variable wettability is a potent material for fast and precise transport of nanodroplets with significant potential for directed nanoscale liquid transport and precision drug delivery.
Collapse
Affiliation(s)
- Ermioni Papadopoulou
- Computational Science and Engineering Laboratory , ETH Zürich , Zürich CH-8092 , Switzerland
| | - Constantine M Megaridis
- Department of Mechanical and Industrial Engineering , University of Illinois at Chicago , Chicago , Illinois 60607 , United States
| | - Jens H Walther
- Computational Science and Engineering Laboratory , ETH Zürich , Zürich CH-8092 , Switzerland
- Department of Mechanical Engineering , Technical University of Denmark , 2800 Kongens Lyngby , Denmark
| | - Petros Koumoutsakos
- Computational Science and Engineering Laboratory , ETH Zürich , Zürich CH-8092 , Switzerland
| |
Collapse
|
26
|
Naeiji P, Woo TK, Alavi S, Varaminian F, Ohmura R. Interfacial properties of hydrocarbon/water systems predicted by molecular dynamic simulations. J Chem Phys 2019; 150:114703. [DOI: 10.1063/1.5078739] [Citation(s) in RCA: 20] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/21/2023] Open
Affiliation(s)
- Parisa Naeiji
- Department of Chemistry and Biomolecular Sciences, University of Ottawa, Ottawa, Ontario K1A 0R6, Canada
- School of Chemical, Petroleum, and Gas Engineering, Semnan University, Semnan, Iran
| | - Tom K. Woo
- Department of Chemistry and Biomolecular Sciences, University of Ottawa, Ottawa, Ontario K1A 0R6, Canada
| | - Saman Alavi
- Department of Chemistry and Biomolecular Sciences, University of Ottawa, Ottawa, Ontario K1A 0R6, Canada
- Department of Mechanical Engineering, Keio University, 3-14-1 Hiyoshi, Kohoku-Ku, Yokohama 223-8522, Japan
| | - Farshad Varaminian
- School of Chemical, Petroleum, and Gas Engineering, Semnan University, Semnan, Iran
| | - Ryo Ohmura
- Department of Mechanical Engineering, Keio University, 3-14-1 Hiyoshi, Kohoku-Ku, Yokohama 223-8522, Japan
| |
Collapse
|
27
|
Li B, Gu Y, Chen M. Cavitation inception of water with solid nanoparticles: A molecular dynamics study. ULTRASONICS SONOCHEMISTRY 2019; 51:120-128. [PMID: 30420302 DOI: 10.1016/j.ultsonch.2018.10.036] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/28/2018] [Revised: 10/28/2018] [Accepted: 10/29/2018] [Indexed: 05/24/2023]
Abstract
Cavitation in liquid with impurities is important in heterogeneous nucleation applications. One of the most widely existing kinds of impurities is solid particles, which can be found in natural water from rivers and specially prepared water such as nanofluids. Understanding the effects caused by the existence of nanoparticles on cavitation in water is vital to the rapidly developed nanotechnologies and medical researches. In this study, cavitation in water with nanoparticles is investigated through molecular dynamics simulations. The effects by nanoparticle materials and sizes on cavitation are discussed by using SiO2 and polyethylene spherical nanoparticles with different diameters. The nucleation rate and the formation of critical bubbles in cavitation are studied via the Voronoi tessellation and the mean first passage time methods. The hydrogen bond network in water is also analyzed. Results reveal that SiO2 and polyethylene nanoparticles may destabilize the hydrogen bond network in water. With the same particle size, cavitation in water with polyethylene nanoparticles is promoted to a greater extent than that with SiO2 nanoparticles. With the same nanoparticle material, cavitation is promoted with the increase in particle size in a range spanning half to ten times the critical bubble radius. Beyond this range, particle size has little influence on cavitation. Reasons for those effects on cavitation due to the presence of solid nanoparticles are discussed by analysing the changes of hydrogen bonds network in water.
Collapse
Affiliation(s)
- Buxuan Li
- Department of Engineering Mechanics, Tsinghua University, Beijing 100084, China
| | - Youwei Gu
- Department of Engineering Mechanics, Tsinghua University, Beijing 100084, China
| | - Min Chen
- Department of Engineering Mechanics, Tsinghua University, Beijing 100084, China.
| |
Collapse
|
28
|
Zhou Y, Li B, Gu Y, Chen M. A molecular dynamics simulation study on the cavitation inception of water with dissolved gases. Mol Phys 2018. [DOI: 10.1080/00268976.2018.1559371] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/27/2022]
Affiliation(s)
- Yusi Zhou
- Department of Engineering Mechanics, Center for Nano and Micro Mechanics, Tsinghua University, Beijing, People’s Republic of China
| | - Buxuan Li
- Department of Engineering Mechanics, Center for Nano and Micro Mechanics, Tsinghua University, Beijing, People’s Republic of China
| | - Youwei Gu
- Department of Engineering Mechanics, Center for Nano and Micro Mechanics, Tsinghua University, Beijing, People’s Republic of China
| | - Min Chen
- Department of Engineering Mechanics, Center for Nano and Micro Mechanics, Tsinghua University, Beijing, People’s Republic of China
| |
Collapse
|
29
|
Ghoufi A, Malfreyt P. Calculation of the surface tension of water: 40 years of molecular simulations. MOLECULAR SIMULATION 2018. [DOI: 10.1080/08927022.2018.1513648] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/05/2023]
Affiliation(s)
- Aziz Ghoufi
- Institut de Physique de Rennes, Université Rennes 1, Rennes, France
| | - Patrice Malfreyt
- Institut de Chimie de Clermont-Ferrand (ICCF), Université Clermont Auvergne, CNRS, SIGMA Clermont, Clermont-Ferrand, France
| |
Collapse
|
30
|
Essafri I, Le breton JC, Saint-Jalmes A, Soldera A, Szymczyk A, Malfreyt P, Ghoufi A. Contact angle and surface tension of water on a hexagonal boron nitride monolayer: a methodological investigation. MOLECULAR SIMULATION 2018. [DOI: 10.1080/08927022.2018.1502427] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/28/2022]
Affiliation(s)
- Ilham Essafri
- Univ Rennes, CNRS, IPR (Institut de Physique de Rennes), Rennes, France
| | | | | | - Armand Soldera
- Laboratory of Physical Chemistry of Matter (LPCM), Department of Chemistry, Université de Sherbrooke, Sherbrooke, Québec, Canada
| | - Anthony Szymczyk
- Univ Rennes, CNRS, ISCR (Institut des Sciences Chimiques de Rennes), Rennes, France
| | - Patrice Malfreyt
- Universite Clermont Auvergne, CNRS, SIGMA Clermont, Institut de Chimie de Clermont-Ferrand (ICCF), Clermont-Ferrand, France
| | - Aziz Ghoufi
- Univ Rennes, CNRS, IPR (Institut de Physique de Rennes), Rennes, France
| |
Collapse
|
31
|
Longford FGJ, Essex JW, Skylaris CK, Frey JG. Unexpected finite size effects in interfacial systems: Why bigger is not always better-Increase in uncertainty of surface tension with bulk phase width. J Chem Phys 2018; 148:214704. [PMID: 29884027 DOI: 10.1063/1.5025887] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022] Open
Abstract
We present an unexpected finite size effect affecting interfacial molecular simulations that is proportional to the width-to-surface-area ratio of the bulk phase Ll/A. This finite size effect has a significant impact on the variance of surface tension values calculated using the virial summation method. A theoretical derivation of the origin of the effect is proposed, giving a new insight into the importance of optimising system dimensions in interfacial simulations. We demonstrate the consequences of this finite size effect via a new way to estimate the surface energetic and entropic properties of simulated air-liquid interfaces. Our method is based on macroscopic thermodynamic theory and involves comparing the internal energies of systems with varying dimensions. We present the testing of these methods using simulations of the TIP4P/2005 water forcefield and a Lennard-Jones fluid model of argon. Finally, we provide suggestions of additional situations, in which this finite size effect is expected to be significant, as well as possible ways to avoid its impact.
Collapse
Affiliation(s)
| | | | | | - Jeremy G Frey
- University of Southampton, Southampton SO17 1BJ, United Kingdom
| |
Collapse
|
32
|
Javanainen M, Lamberg A, Cwiklik L, Vattulainen I, Ollila OHS. Atomistic Model for Nearly Quantitative Simulations of Langmuir Monolayers. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2018; 34:2565-2572. [PMID: 28945973 DOI: 10.1021/acs.langmuir.7b02855] [Citation(s) in RCA: 46] [Impact Index Per Article: 7.7] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/03/2023]
Abstract
Lung surfactant and a tear film lipid layer are examples of biologically relevant macromolecular structures found at the air-water interface. Because of their complexity, they are often studied in terms of simplified lipid layers, the simplest example being a Langmuir monolayer. Given the profound biological significance of these lipid assemblies, there is a need to understand their structure and dynamics on the nanoscale, yet there are not many techniques able to provide this information. Atomistic molecular dynamics simulations would be a tool fit for this purpose; however, the simulation models suggested until now have been qualitative instead of quantitative. This limitation has mainly stemmed from the challenge to correctly describe the surface tension of water with simulation parameters compatible with other biomolecules. In this work, we show that this limitation can be overcome by using the recently introduced four-point OPC water model, whose surface tension for water is demonstrated to be quantitatively consistent with experimental data and which is also shown to be compatible with the commonly employed lipid models. We further establish that the approach of combining the OPC four-point water model with the CHARMM36 lipid force field provides nearly quantitative agreement with experiments for the surface pressure-area isotherm for POPC and DPPC monolayers, also including the experimentally observed phase coexistence in a DPPC monolayer. The simulation models reported in this work pave the way for nearly quantitative atomistic studies of lipid-rich biological structures at air-water interfaces.
Collapse
Affiliation(s)
- Matti Javanainen
- Laboratory of Physics, Tampere University of Technology , 33101 Tampere, Finland
| | | | - Lukasz Cwiklik
- J. Heyrovský Institute of Physical Chemistry of the Czech Academy of Sciences , 182 23 Prague 8, Czech Republic
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences , 166 10 Prague 6, Czech Republic
| | - Ilpo Vattulainen
- Laboratory of Physics, Tampere University of Technology , 33101 Tampere, Finland
- MEMPHYS - Center for Biomembrane Physics, www.memphys.dk
| | - O H Samuli Ollila
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences , 166 10 Prague 6, Czech Republic
| |
Collapse
|
33
|
Nandi PK, Burnham CJ, English NJ. Electro-suppression of water nano-droplets’ solidification in no man’s land: Electromagnetic fields’ entropic trapping of supercooled water. J Chem Phys 2018; 148:044503. [DOI: 10.1063/1.5004509] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Prithwish K. Nandi
- Irish Centre for High-End Computing, Grand Canal Quay, Dublin 2, Ireland
- School of Chemical and Bioprocess Engineering, University College Dublin, Belfield, Dublin 4, Ireland
| | - Christian J. Burnham
- School of Chemical and Bioprocess Engineering, University College Dublin, Belfield, Dublin 4, Ireland
| | - Niall J. English
- School of Chemical and Bioprocess Engineering, University College Dublin, Belfield, Dublin 4, Ireland
| |
Collapse
|
34
|
Sinha N, Singh JK. Effect of nanoparticles on vapour-liquid surface tension of water: A molecular dynamics study. J Mol Liq 2017. [DOI: 10.1016/j.molliq.2017.09.059] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
|
35
|
Sega M, Dellago C. Long-Range Dispersion Effects on the Water/Vapor Interface Simulated Using the Most Common Models. J Phys Chem B 2017; 121:3798-3803. [DOI: 10.1021/acs.jpcb.6b12437] [Citation(s) in RCA: 28] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Marcello Sega
- Faculty of Physics, University of Vienna, Boltzmanngasse 5, A-1090, Vienna, Austria
| | - Christoph Dellago
- Faculty of Physics, University of Vienna, Boltzmanngasse 5, A-1090, Vienna, Austria
| |
Collapse
|
36
|
Fugel M, Weiss VC. A corresponding-states analysis of the liquid-vapor equilibrium properties of common water models. J Chem Phys 2017; 146:064505. [DOI: 10.1063/1.4975778] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022] Open
|
37
|
Reddy SK, Straight SC, Bajaj P, Huy Pham C, Riera M, Moberg DR, Morales MA, Knight C, Götz AW, Paesani F. On the accuracy of the MB-pol many-body potential for water: Interaction energies, vibrational frequencies, and classical thermodynamic and dynamical properties from clusters to liquid water and ice. J Chem Phys 2016; 145:194504. [DOI: 10.1063/1.4967719] [Citation(s) in RCA: 161] [Impact Index Per Article: 20.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022] Open
Affiliation(s)
- Sandeep K. Reddy
- Department of Chemistry and Biochemistry, University of California, San Diego, La Jolla, California 92093, USA
| | - Shelby C. Straight
- Department of Chemistry and Biochemistry, University of California, San Diego, La Jolla, California 92093, USA
| | - Pushp Bajaj
- Department of Chemistry and Biochemistry, University of California, San Diego, La Jolla, California 92093, USA
| | - C. Huy Pham
- Department of Chemistry and Biochemistry, University of California, San Diego, La Jolla, California 92093, USA
| | - Marc Riera
- Department of Chemistry and Biochemistry, University of California, San Diego, La Jolla, California 92093, USA
| | - Daniel R. Moberg
- Department of Chemistry and Biochemistry, University of California, San Diego, La Jolla, California 92093, USA
| | - Miguel A. Morales
- Lawrence Livermore National Laboratory, 7000 East Avenue, Livermore, California 94550, USA
| | - Chris Knight
- Leadership Computing Facility, Argonne National Laboratory, 9700 South Cass Avenue, Argonne, Illinois 60439, USA
| | - Andreas W. Götz
- San Diego Supercomputer Center, University of California, San Diego, La Jolla, California 92093, USA
| | - Francesco Paesani
- Department of Chemistry and Biochemistry, University of California, San Diego, La Jolla, California 92093, USA
| |
Collapse
|
38
|
Wang H, Gao X, Fang J. Multiple Staggered Mesh Ewald: Boosting the Accuracy of the Smooth Particle Mesh Ewald Method. J Chem Theory Comput 2016; 12:5596-5608. [PMID: 27760290 DOI: 10.1021/acs.jctc.6b00701] [Citation(s) in RCA: 28] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The smooth particle mesh Ewald (SPME) method is the standard method for computing the electrostatic interactions in the molecular simulations. In this work, we develop the multiple staggered mesh Ewald (MSME) method, which averages the SPME forces computed on, for example, M, staggered meshes. We prove, from a theoretical perspective, that the MSME is as accurate as the SPME, but uses M2 times less mesh points in a certain parameter range. In the complementary parameter range, the MSME is as accurate as the SPME with twice the interpolation order. The theoretical conclusions are numerically validated both by a uniform and uncorrelated charge system, and by a three-point-charge water system that is widely used as solvent for the biomacromolecules.
Collapse
Affiliation(s)
- Han Wang
- Institute of Applied Physics and Computational Mathematics , Fenghao East Road 2, Beijing 100094, P.R. China.,CAEP Software Center for High Performance Numerical Simulation , Huayuan Road 6, Beijing 100088, P.R. China
| | - Xingyu Gao
- Institute of Applied Physics and Computational Mathematics , Fenghao East Road 2, Beijing 100094, P.R. China.,CAEP Software Center for High Performance Numerical Simulation , Huayuan Road 6, Beijing 100088, P.R. China.,Laboratory of Computational Physics, Huayuan Road 6, Beijing 100088, P.R. China
| | - Jun Fang
- Institute of Applied Physics and Computational Mathematics , Fenghao East Road 2, Beijing 100094, P.R. China.,CAEP Software Center for High Performance Numerical Simulation , Huayuan Road 6, Beijing 100088, P.R. China
| |
Collapse
|
39
|
Sega M, Fábián B, Jedlovszky P. Pressure Profile Calculation with Mesh Ewald Methods. J Chem Theory Comput 2016; 12:4509-15. [DOI: 10.1021/acs.jctc.6b00576] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Marcello Sega
- Computational
Physics Group, University of Vienna, Sensengasse 8/9, 1090 Vienna, Austria
| | - Balázs Fábián
- Institut
UTINAM (CNRS UMR 6213), Université Bourgogne Franche Comté, 16 route de Gray, F-25030 Besançon, France
- Department
of Inorganic and Analytical Chemistry, Budapest University of Technology and Economics, Szt. Gellért tér 4, H-1111 Budapest, Hungary
| | - Pál Jedlovszky
- Department
of Chemistry, Eszterházy Károly University, Leányka
u. 6, H-3300 Eger, Hungary
- Laboratory
of Interfaces and Nanosize Systems, Institute of Chemistry, Eötvös Loránd University, Pázmány P. Stny 1/A, H-1117 Budapest, Hungary
- MTA-BME Research
Group of Technical Analytical Chemistry, Szt. Gellért tér 4, H-1111 Budapest, Hungary
| |
Collapse
|
40
|
Rick SW. A polarizable, charge transfer model of water using the drude oscillator. J Comput Chem 2016; 37:2060-6. [DOI: 10.1002/jcc.24426] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/27/2015] [Revised: 04/14/2016] [Accepted: 05/17/2016] [Indexed: 12/24/2022]
Affiliation(s)
- Steven W. Rick
- Department of ChemistryUniversity of New OrleansNew Orleans70148 Los Angeles
| |
Collapse
|
41
|
Hu H, Wang F. The liquid-vapor equilibria of TIP4P/2005 and BLYPSP-4F water models determined through direct simulations of the liquid-vapor interface. J Chem Phys 2016; 142:214507. [PMID: 26049508 DOI: 10.1063/1.4922166] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022] Open
Abstract
In this paper, the surface tension and critical properties for the TIP4P/2005 and BLYPSP-4F models are reported. A clear dependence of surface tension on the van der Waals cutoff radius (rvdw) is shown when van der Waals interactions are modeled with a simple cutoff scheme. A linear extrapolation formula is proposed that can be used to determine the infinite rvdw surface tension through a few simulations with finite rvdw. A procedure for determining liquid and vapor densities is proposed that does not require fitting to a profile function. Although the critical temperature of water is also found to depend on the choice of rvdw, the dependence is weaker. We argue that a rvdw of 1.75 nm is a good compromise for water simulations when long-range van der Waals correction is not applied. Since the majority of computational programs do not support rigorous treatment of long-range dispersion, the establishment of a minimal acceptable rvdw is important for the simulation of a variety of inhomogeneous systems, such as water bubbles, and water in confined environments. The BLYPSP-4F model predicts room temperature surface tension marginally better than TIP4P/2005 but overestimates the critical temperature. This is expected since only liquid configurations were fit during the development of the BLYPSP-4F potential. The potential is expected to underestimate the stability of vapor and thus overestimate the region of stability for the liquid.
Collapse
Affiliation(s)
- Hongyi Hu
- Department of Chemistry and Biochemistry, University of Arkansas, Fayetteville, Arkansas 72701, USA
| | - Feng Wang
- Department of Chemistry and Biochemistry, University of Arkansas, Fayetteville, Arkansas 72701, USA
| |
Collapse
|
42
|
Wilhelmsen Ø, Trinh TT, Lervik A, Badam VK, Kjelstrup S, Bedeaux D. Coherent description of transport across the water interface: From nanodroplets to climate models. Phys Rev E 2016; 93:032801. [PMID: 27078427 DOI: 10.1103/physreve.93.032801] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/28/2015] [Indexed: 11/07/2022]
Abstract
Transport of mass and energy across the vapor-liquid interface of water is of central importance in a variety of contexts such as climate models, weather forecasts, and power plants. We provide a complete description of the transport properties of the vapor-liquid interface of water with the framework of nonequilibrium thermodynamics. Transport across the planar interface is then described by 3 interface transfer coefficients where 9 more coefficients extend the description to curved interfaces. We obtain all coefficients in the range 260-560 K by taking advantage of water evaporation experiments at low temperatures, nonequilibrium molecular dynamics with the TIP4P/2005 rigid-water-molecule model at high temperatures, and square gradient theory to represent the whole range. Square gradient theory is used to link the region where experiments are possible (low vapor pressures) to the region where nonequilibrium molecular dynamics can be done (high vapor pressures). This enables a description of transport across the planar water interface, interfaces of bubbles, and droplets, as well as interfaces of water structures with complex geometries. The results are likely to improve the description of evaporation and condensation of water at widely different scales; they open a route to improve the understanding of nanodroplets on a small scale and the precision of climate models on a large scale.
Collapse
Affiliation(s)
- Øivind Wilhelmsen
- Department of Chemistry, Norwegian University of Science and Technology, Trondheim, Norway
| | - Thuat T Trinh
- Department of Chemistry, Norwegian University of Science and Technology, Trondheim, Norway
| | - Anders Lervik
- Department of Chemistry, Norwegian University of Science and Technology, Trondheim, Norway
| | - Vijay Kumar Badam
- Institute of Fluid Mechanics (LSTM), Friedrich-Alexander-Universität Erlangen-Nürnberg, Cauerstrasse 4, D-91058 Erlangen, Germany
| | - Signe Kjelstrup
- Department of Chemistry, Norwegian University of Science and Technology, Trondheim, Norway
| | - Dick Bedeaux
- Department of Chemistry, Norwegian University of Science and Technology, Trondheim, Norway
| |
Collapse
|
43
|
Ghoufi A, Malfreyt P, Tildesley DJ. Computer modelling of the surface tension of the gas–liquid and liquid–liquid interface. Chem Soc Rev 2016; 45:1387-409. [DOI: 10.1039/c5cs00736d] [Citation(s) in RCA: 134] [Impact Index Per Article: 16.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/03/2023]
Abstract
This review presents the state of the art in molecular simulations of interfacial systems and of the calculation of the surface tension from the underlying intermolecular potential.
Collapse
Affiliation(s)
- Aziz Ghoufi
- Institut de Physique de Rennes
- UMR CNRS 6251
- 35042 Rennes
- France
| | - Patrice Malfreyt
- Institut de Chimie de Clermont-Ferrand
- ICCF
- CNRS
- UMR 6296
- F-63000 Clermont-Ferrand
| | | |
Collapse
|
44
|
Zhang J, Borg MK, Sefiane K, Reese JM. Wetting and evaporation of salt-water nanodroplets: A molecular dynamics investigation. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2015; 92:052403. [PMID: 26651708 DOI: 10.1103/physreve.92.052403] [Citation(s) in RCA: 23] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/30/2015] [Indexed: 06/05/2023]
Abstract
We employ molecular dynamics simulations to study the wetting and evaporation of salt-water nanodroplets on platinum surfaces. Our results show that the contact angle of the droplets increases with the salt concentration. To verify this, a second simulation system of a thin salt-water film on a platinum surface is used to calculate the various surface tensions. We find that both the solid-liquid and liquid-vapor surface tensions increase with salt concentration and as a result these cause an increase in the contact angle. However, the evaporation rate of salt-water droplets decreases as the salt concentration increases, due to the hydration of salt ions. When the water molecules have all evaporated from the droplet, two forms of salt crystals are deposited, clump and ringlike, depending on the solid-liquid interaction strength and the evaporation rate. To form salt crystals in a ring, it is crucial that there is a pinned stage in the evaporation process, during which salt ions can move from the center to the rim of the droplets. With a stronger solid-liquid interaction strength, a slower evaporation rate, and a higher salt concentration, a complete salt crystal ring can be deposited on the surface.
Collapse
Affiliation(s)
- Jun Zhang
- James Weir Fluids Laboratory, Department of Mechanical and Aerospace Engineering, University of Strathclyde, Glasgow G1 1XJ, United Kingdom
| | - Matthew K Borg
- School of Engineering, University of Edinburgh, Edinburgh EH9 3FB, United Kingdom
| | - Khellil Sefiane
- School of Engineering, University of Edinburgh, Edinburgh EH9 3FB, United Kingdom
- Tianjin Key Laboratory of Refrigeration Technology, Tianjin University of Commerce, Tianjin, China
| | - Jason M Reese
- School of Engineering, University of Edinburgh, Edinburgh EH9 3FB, United Kingdom
| |
Collapse
|
45
|
Goujon F, Ghoufi A, Malfreyt P, Tildesley DJ. Controlling the Long-Range Corrections in Atomistic Monte Carlo Simulations of Two-Phase Systems. J Chem Theory Comput 2015; 11:4573-85. [DOI: 10.1021/acs.jctc.5b00377] [Citation(s) in RCA: 23] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/27/2022]
Affiliation(s)
- Florent Goujon
- Institut
de Chimie de Clermont-Ferrand, Université Clermont Auvergne, Université Blaise Pascal, BP 10448, F-63000 Clermont-Ferrand, France
| | - Aziz Ghoufi
- Institut
de Physique de Rennes, Université Rennes 1, 35042 Rennes, France
| | - Patrice Malfreyt
- Institut
de Chimie de Clermont-Ferrand, Université Clermont Auvergne, Université Blaise Pascal, BP 10448, F-63000 Clermont-Ferrand, France
| | | |
Collapse
|
46
|
Martínez-Ruiz FJ, Blas FJ, Mendiboure B, Moreno-Ventas Bravo AI. Effect of dispersive long-range corrections to the pressure tensor: The vapour-liquid interfacial properties of the Lennard-Jones system revisited. J Chem Phys 2014; 141:184701. [DOI: 10.1063/1.4900773] [Citation(s) in RCA: 28] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Affiliation(s)
- F. J. Martínez-Ruiz
- Departamento de Física Aplicada, Universidad de Huelva, 21071 Huelva, Spain
- Centro de Investigación de Física Teórica y Matemática, Universidad de Huelva, 21071 Huelva, Spain
| | - F. J. Blas
- Departamento de Física Aplicada, Universidad de Huelva, 21071 Huelva, Spain
- Centro de Investigación de Física Teórica y Matemática, Universidad de Huelva, 21071 Huelva, Spain
| | - B. Mendiboure
- Laboratoire des Fluides Complexes et leurs Réservoirs, UMR5150, Université de Pau et des Pays de l’Adour, B. P. 1155, Pau Cedex 64014, France
| | - A. I. Moreno-Ventas Bravo
- Centro de Investigación de Física Teórica y Matemática, Universidad de Huelva, 21071 Huelva, Spain
- Departamento de Geología, Facultad de Ciencias Experimentales, Universidad de Huelva, 21071 Huelva, Spain
| |
Collapse
|
47
|
Pascal TA, Goddard WA. Interfacial Thermodynamics of Water and Six Other Liquid Solvents. J Phys Chem B 2014; 118:5943-56. [DOI: 10.1021/jp410861h] [Citation(s) in RCA: 29] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Affiliation(s)
- Tod A. Pascal
- Materials and Process Simulation
Center, California Institute of Technology, Pasadena, California 91125, United States
| | - William A. Goddard
- Materials and Process Simulation
Center, California Institute of Technology, Pasadena, California 91125, United States
| |
Collapse
|
48
|
Fuentes-Azcatl R, Alejandre J. Non-Polarizable Force Field of Water Based on the Dielectric Constant: TIP4P/ε. J Phys Chem B 2014; 118:1263-72. [DOI: 10.1021/jp410865y] [Citation(s) in RCA: 85] [Impact Index Per Article: 8.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Raúl Fuentes-Azcatl
- Departamento de Química, Universidad Autónoma Metropolitana-Iztapalapa, Av. San Rafael Atlixco 186, Col.
Vicentina, 09340 México D.F., México
| | - José Alejandre
- Departamento de Química, Universidad Autónoma Metropolitana-Iztapalapa, Av. San Rafael Atlixco 186, Col.
Vicentina, 09340 México D.F., México
| |
Collapse
|
49
|
Affiliation(s)
- Mario Orsi
- School of Engineering and Materials Science, Queen Mary University of London, London, UK
| |
Collapse
|
50
|
Bárcenas M, Odriozola G, Orea P. Structure and coexistence properties of shoulder–square well fluids. J Mol Liq 2013. [DOI: 10.1016/j.molliq.2012.10.027] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/27/2022]
|