1
|
Valle JVL, Mendonça BHS, Barbosa MC, Chacham H, de Moraes EE. Accuracy of TIP4P/2005 and SPC/Fw Water Models. J Phys Chem B 2024; 128:1091-1097. [PMID: 38253517 DOI: 10.1021/acs.jpcb.3c07044] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/24/2024]
Abstract
Water is used as the main solvent in model systems containing bioorganic molecules. Choosing the right water model is an important step in the study of the biophysical and biochemical processes that occur in cells. In the present work, we perform molecular dynamics simulations using two distinct force fields for water: the rigid model TIP4P/2005, where only intermolecular interactions are considered, and the flexible model SPC/Fw, where intramolecular interactions are also taken into account. The simulations aim to determine the effect of the inclusion of intramolecular interactions on the accuracy of calculated properties of bulk water (density and thermal expansion coefficient, self-diffusion coefficients, shear viscosity, radial distribution functions, and dielectric constant), as compared to experimental results, over a temperature range between 250 and 370 K. We find that the results of the rigid model present the smallest deviations relative to experiments for most of the calculated quantities, except for the shear viscosity of supercooled water and the water dielectric constant, where the flexible model presents better agreement with experiments.
Collapse
Affiliation(s)
- João V L Valle
- Instituto de Física, Universidade Federal da Bahia, Campus Universitário de Ondina, Salvador 40210-340, BA, Brazil
| | - Bruno H S Mendonça
- Departamento de Física, ICEX, Universidade Federal de Minas Gerais, CP 702, Belo Horizonte 30123-970, MG, Brazil
| | - Marcia C Barbosa
- Instituto de Física, Universidade Federal do Rio Grande do Sul, Porto Alegre 91501-970, RS, Brazil
| | - Helio Chacham
- Departamento de Física, ICEX, Universidade Federal de Minas Gerais, CP 702, 30123-970 Belo Horizonte, MG, Brazil
| | - Elizane E de Moraes
- Instituto de Física, Universidade Federal da Bahia, Campus Universitário de Ondina, Salvador 40210-340, BA, Brazil
| |
Collapse
|
2
|
Mendonça BHS, de Moraes EE, Kirch A, Batista RJC, de Oliveira AB, Barbosa MC, Chacham H. Flow through Deformed Carbon Nanotubes Predicted by Rigid and Flexible Water Models. J Phys Chem B 2023; 127:8634-8643. [PMID: 37754781 DOI: 10.1021/acs.jpcb.3c02889] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/28/2023]
Abstract
In this study, using nonequilibrium molecular dynamics simulation, the flow of water in deformed carbon nanotubes is studied for two water models TIP4P/2005 and simple point charge/FH (SPC/FH). The results demonstrated a nonuniform dependence of the flow on the tube deformation and the flexibility imposed on the water molecules, leading to an unexpected increase in the flow in some cases. The effects of the tube diameter and pressure gradient are investigated to explain the abnormal flow behavior with different degrees of structural deformation.
Collapse
Affiliation(s)
- Bruno H S Mendonça
- Departamento de Física, ICEX, Universidade Federal de Minas Gerais, CP 702, Belo Horizonte 30123-970, MG, Brazil
| | - Elizane E de Moraes
- Instituto de Física, Universidade Federal da Bahia, Campus Universitário de Ondina, Salvador 40210-340, BA, Brazil
| | - Alexsandro Kirch
- Instituto de Física, Universidade de São Paulo, CP 66318, São Paulo 05315-970, SP, Brazil
| | - Ronaldo J C Batista
- Departamento de Física, Universidade Federal de Ouro Preto, Campus Morro do Cruzeiro, Ouro Preto 35400-000, MG, Brazil
| | - Alan B de Oliveira
- Departamento de Física, Universidade Federal de Ouro Preto, Campus Morro do Cruzeiro, Ouro Preto 35400-000, MG, Brazil
| | - Marcia C Barbosa
- Instituto de Física, Universidade Federal do Rio Grande do Sul, Porto Alegre 91501-970, RS, Brazil
| | - Hélio Chacham
- Departamento de Física, ICEX, Universidade Federal de Minas Gerais, CP 702, Belo Horizonte 30123-970, MG, Brazil
| |
Collapse
|
3
|
Wilson BA, Nielsen SO, Randrianalisoa J, Qin Z. Curvature and temperature-dependent thermal interface conductance between nanoscale-gold and water. J Chem Phys 2022; 157:054703. [PMID: 35933210 PMCID: PMC9355664 DOI: 10.1063/5.0090683] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
<p>Plasmonic gold nanoparticles (AuNPs) can convert laser irradiation into thermal energy for a variety of applications. Although heat transfer through the AuNP-water interface is considered an essential part of the plasmonic heating process, there is a lack of mechanistic understanding of how interface curvature and the heating itself impact interfacial heat transfer. Here, we report atomistic molecular dynamics simulations that investigate heat transfer through nanoscale gold-water interfaces. We simulated four nanoscale gold structures under various applied heat flux to evaluate how gold-water interface curvature and temperature affect the interfacial heat transfer. We also considered a case in which we artificially reduced wetting at the gold surfaces by tuning the gold-water interactions to determine if such a perturbation alters the curvature and temperature dependence of the gold-water interfacial heat transfer. We first confirmed that interfacial heat transfer is particularly important for small particles (diameter {less than or equal to} 10 nm). We found that the thermal interface conductance increases linearly with interface curvature regardless of the gold wettability, while it increases non-linearly with the applied heat flux under normal wetting and remains constant under reduced wetting. Our analysis suggests the curvature dependence of the interface conductance coincides with changes in interfacial water adsorption, while the temperature dependence may arise from temperature-induced shifts in the distribution of water vibrational states. Our study advances the current understanding of interface thermal conductance for a broad range of applications.
Collapse
Affiliation(s)
- Blake A Wilson
- Chemistry, The University of Texas at Dallas, United States of America
| | - Steven O. Nielsen
- Department of Chemistry, University of Texas at Dallas, United States of America
| | | | - Zhenpeng Qin
- Mechanical Engineering, The University of Texas at Dallas, United States of America
| |
Collapse
|
4
|
Chen Y, Schultz AJ, Errington JR. Coupled Monte Carlo and Molecular Dynamics Simulations on Interfacial Properties of Antifouling Polymer Membranes. J Phys Chem B 2021; 125:8193-8204. [PMID: 34259529 DOI: 10.1021/acs.jpcb.1c01966] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
We use molecular simulation to study the wetting behavior of antifouling polymer-tethered membranes. We obtain the interfacial properties (e.g., contact angle) of water at various temperatures for five polymer membranes, including a base polysulfone (PSF) membrane and four other PSF membranes grafted with antifouling polymers (two poly(ethylene glycol) (PEG) tethers and two zwitterionic tethers). We implement a coupled Monte Carlo (MC)/molecular dynamics (MD) approach to determine the interface potentials of water on the membrane surfaces in an efficient manner. Within this method, short MC and MD simulations are performed in cycles to collect the surface excess free energy of a thin water film on polymer membrane surfaces. Simulation results show that the grafting of zwitterionic tethers provides a more significant enhancement in the hydrophilicity of the PSF membrane than that of the PEG tethers. Water completely wets the surface of zwitterionic polymer membranes.
Collapse
Affiliation(s)
- Yiqi Chen
- Department of Chemical and Biological Engineering, University at Buffalo, Buffalo, New York 14260-4200, United States
| | - Andrew J Schultz
- Department of Chemical and Biological Engineering, University at Buffalo, Buffalo, New York 14260-4200, United States
| | - Jeffrey R Errington
- Department of Chemical and Biological Engineering, University at Buffalo, Buffalo, New York 14260-4200, United States
| |
Collapse
|
5
|
Ta TD, Tieu AK, Tran BH. Hydroxyl Influence on Adsorption and Lubrication of an Ultrathin Aqueous Triblock Copolymer Lubricant. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2021; 37:1465-1479. [PMID: 33476165 DOI: 10.1021/acs.langmuir.0c02987] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/12/2023]
Abstract
This research aims to provide insights into the adsorption behaviors of two monomers of triblock copolymers (1,2-dimethoxyethane (1,2-DME) and 1,2-dimethoxypropane (1,2-DMP)) on a TiO2 surface in aqueous solution. A multiscale theoretical framework by means of the density functional theory (DFT), ab initio molecular dynamics (AIMD), and classical molecular dynamics (MD) simulations is established. The DFT calculation confirms that these molecules adsorb more energetically on a hydroxylated surface than pure oxide. There is a difference in adsorption behaviors between 1,2-DMP and 1,2-DME molecules due to the covalent bonding between carbons and oxygen of the hydroxylated TiO2 surface. The AIMD simulation reveals that the adsorption of both copolymers to the TiO2 surface is hindered by the presence of water with 1,2-DME exhibiting a weaker adsorption than 1,2-DMP. The presence of 1,2-DME on the TiO2 surface with water produced a smaller number of hydroxyl groups on the surface than 1,2-DMP. Moreover, the dissociative adsorption of water onto the rutile surface is the main cause for a chemical formation of terminating hydroxyl groups. The number of associated bonds is insignificant compared to the dissociated one since the dissociative adsorption is more favored than the associative one. MD simulation indicates that triblock copolymers adsorb stronger on the hydroxylated surface with a thinner adsorbed film thickness than that on the pure rutile. The presence of terminal hydroxyl groups on the rutile surface helps reducing the friction for aqueous 17R2 triblock copolymers, while it results in an increase of friction for normal copolymer L62.
Collapse
Affiliation(s)
- Thi D Ta
- School of Mechanical, Materials Mechatronic and Biomedical Engineering, Faculty of Engineering and Information Sciences (EIS), University of Wollongong, Northfield Avenue, Wollongong, New South Wales 2522, Australia
| | - A Kiet Tieu
- School of Mechanical, Materials Mechatronic and Biomedical Engineering, Faculty of Engineering and Information Sciences (EIS), University of Wollongong, Northfield Avenue, Wollongong, New South Wales 2522, Australia
| | - Bach H Tran
- School of Mechanical, Materials Mechatronic and Biomedical Engineering, Faculty of Engineering and Information Sciences (EIS), University of Wollongong, Northfield Avenue, Wollongong, New South Wales 2522, Australia
| |
Collapse
|
6
|
Mosallanejad S, Oluwoye I, Altarawneh M, Gore J, Dlugogorski BZ. Interfacial and bulk properties of concentrated solutions of ammonium nitrate. Phys Chem Chem Phys 2020; 22:27698-27712. [PMID: 33242055 DOI: 10.1039/d0cp04874g] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
We conducted molecular dynamics (MD) simulations to calculate the density and surface tension of concentrated ammonium nitrate (AN) solutions up to the solubility limit of ammonium nitrate in water, by combining the SPC/E, SPCE/F and TIP4P/2005 water models with OPLS model for ammonium and nitrate ions. This is the first time that the properties of concentrated solutions of nitrates, especially AN, have been studied by molecular dynamics. We effectively account for the polarisation effects by the electronic continuum correction (ECC), practically realised via rescaling of the ionic charges. We found that, the full-charge force field MD simulations overestimate the experimental results, as the ions experience repulsion from the interface and prefer to remain in the subsurface layer and the bulk solution. In contrast, reducing the ionic charges results in the behaviour that fits well with the experimental data. The nitrate anions display a greater propensity for the interface than the ammonium cations. We accurately predict both the density and the rise in the surface tension of concentrated solutions of AN, recommending TIP4P/2005 for water and the scaled-charge OPLS model (OPLS/ECC) for the ions in the solutions. We observe that, the adsorption of anions to the interface accompanies their depletion in the subsurface layer, which is preferentially occupied by cations, resulting in the formation of the electric double layer. We demonstrate the ion deficiency for up to 3 Å below the surface and establish the requirement to include the polarisability effects in the OPLS model for AN. While these results confirmed the findings of the previous studies for dilute solutions, they are new in the solubility limit. Concentrated solutions exhibit a strong effect of the abundance of solute on the coordination numbers of ions and on the degree of ion pairing. Surprisingly, ion pairing decreases significantly at the interface compared with the bulk. The present study identifies OPLS/ECC, along with TIP4P/2005, to yield accurate predictions of physical properties of concentrated AN, with precision required for industrial applications, such as a formulation of emulsion and fuel-oil explosives that now predominate the civilian use of AN. An application of this model will allow one to predict the surface properties of supersaturated solutions of AN which fall outside the capability of the present laboratory experiments but are important industrially.
Collapse
Affiliation(s)
- Sara Mosallanejad
- Discipline of Chemistry and Physics, College of Science, Health, Engineering and Education, Murdoch University, 90 South Street, Murdoch, WA 6150, Australia
| | | | | | | | | |
Collapse
|
7
|
Lynch C, Rao S, Sansom MSP. Water in Nanopores and Biological Channels: A Molecular Simulation Perspective. Chem Rev 2020; 120:10298-10335. [PMID: 32841020 PMCID: PMC7517714 DOI: 10.1021/acs.chemrev.9b00830] [Citation(s) in RCA: 88] [Impact Index Per Article: 22.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/23/2019] [Indexed: 12/18/2022]
Abstract
This Review explores the dynamic behavior of water within nanopores and biological channels in lipid bilayer membranes. We focus on molecular simulation studies, alongside selected structural and other experimental investigations. Structures of biological nanopores and channels are reviewed, emphasizing those high-resolution crystal structures, which reveal water molecules within the transmembrane pores, which can be used to aid the interpretation of simulation studies. Different levels of molecular simulations of water within nanopores are described, with a focus on molecular dynamics (MD). In particular, models of water for MD simulations are discussed in detail to provide an evaluation of their use in simulations of water in nanopores. Simulation studies of the behavior of water in idealized models of nanopores have revealed aspects of the organization and dynamics of nanoconfined water, including wetting/dewetting in narrow hydrophobic nanopores. A survey of simulation studies in a range of nonbiological nanopores is presented, including carbon nanotubes, synthetic nanopores, model peptide nanopores, track-etched nanopores in polymer membranes, and hydroxylated and functionalized nanoporous silica. These reveal a complex relationship between pore size/geometry, the nature of the pore lining, and rates of water transport. Wider nanopores with hydrophobic linings favor water flow whereas narrower hydrophobic pores may show dewetting. Simulation studies over the past decade of the behavior of water in a range of biological nanopores are described, including porins and β-barrel protein nanopores, aquaporins and related polar solute pores, and a number of different classes of ion channels. Water is shown to play a key role in proton transport in biological channels and in hydrophobic gating of ion channels. An overall picture emerges, whereby the behavior of water in a nanopore may be predicted as a function of its hydrophobicity and radius. This informs our understanding of the functions of diverse channel structures and will aid the design of novel nanopores. Thus, our current level of understanding allows for the design of a nanopore which promotes wetting over dewetting or vice versa. However, to design a novel nanopore, which enables fast, selective, and gated flow of water de novo would remain challenging, suggesting a need for further detailed simulations alongside experimental evaluation of more complex nanopore systems.
Collapse
Affiliation(s)
- Charlotte
I. Lynch
- Department of Biochemistry, University of Oxford, South Parks Road, Oxford, OX1 3QU, U.K.
| | - Shanlin Rao
- Department of Biochemistry, University of Oxford, South Parks Road, Oxford, OX1 3QU, U.K.
| | - Mark S. P. Sansom
- Department of Biochemistry, University of Oxford, South Parks Road, Oxford, OX1 3QU, U.K.
| |
Collapse
|
8
|
Ranathunga DTS, Shamir A, Dai X, Nielsen SO. Molecular Dynamics Simulations of Water Condensation on Surfaces with Tunable Wettability. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2020; 36:7383-7391. [PMID: 32498521 DOI: 10.1021/acs.langmuir.0c00915] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/11/2023]
Abstract
Water condensation plays a major role in a wide range of industrial applications. Over the past few years, many studies have shown interest in designing surfaces with enhanced water condensation and removal properties. It is well known that heterogeneous nucleation outperforms homogeneous nucleation in the condensation process. Because heterogeneous nucleation initiates on a surface at a small scale, it is highly desirable to characterize water-surface interactions at the molecular level. Molecular dynamics (MD) simulations can provide direct insight into heterogeneous nucleation and advance surface designs. Existing MD simulations of water condensation on surfaces were conducted by tuning the solid-water van der Waals interaction energy as a substitute for modeling surfaces with different wettabilities. However, this approach cannot reflect the real intermolecular interactions between the surface and water molecules. Here, we report MD simulations of water condensation on realistic surfaces of alkanethiol self-assembled monolayers with different head group chemistries. We show that decreasing surface hydrophobicity significantly increases the electrostatic forces between water molecules and the surface, thus increasing the water condensation rate. We observe a strong correlation between our rate of condensation results and the results from other surface characterization metrics, such as the interfacial thermal conductance, contact angle, and the molecular-scale wettability metric of Garde and co-workers. This work provides insight into the water condensation process at the molecular scale on surfaces with tunable wettability.
Collapse
Affiliation(s)
- Dineli T S Ranathunga
- Department of Chemistry and Biochemistry, The University of Texas at Dallas, 800 West Campbell Road, Richardson, Texas 75080, United States
| | - Alexandra Shamir
- Department of Chemistry and Biochemistry, The University of Texas at Dallas, 800 West Campbell Road, Richardson, Texas 75080, United States
| | - Xianming Dai
- Department of Mechanical Engineering, The University of Texas at Dallas, 800 West Campbell Road, Richardson, Texas 75080, United States
| | - Steven O Nielsen
- Department of Chemistry and Biochemistry, The University of Texas at Dallas, 800 West Campbell Road, Richardson, Texas 75080, United States
| |
Collapse
|
9
|
Zhang C, Giberti F, Sevgen E, de Pablo JJ, Gygi F, Galli G. Dissociation of salts in water under pressure. Nat Commun 2020; 11:3037. [PMID: 32546791 PMCID: PMC7298052 DOI: 10.1038/s41467-020-16704-9] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/28/2019] [Accepted: 05/15/2020] [Indexed: 11/09/2022] Open
Abstract
The investigation of salts in water at extreme conditions is crucial to understanding the properties of aqueous fluids in the Earth. We report first principles (FP) and classical molecular dynamics simulations of NaCl in the dilute limit, at temperatures and pressures relevant to the Earth’s upper mantle. Similar to ambient conditions, we observe two metastable states of the salt: the contact (CIP) and the solvent-shared ion-pair (SIP), which are entropically and enthalpically favored, respectively. We find that the free energy barrier between the CIP and SIP minima increases at extreme conditions, and that the stability of the CIP is enhanced in FP simulations, consistent with the decrease of the dielectric constant of water. The minimum free energy path between the CIP and SIP becomes smoother at high pressure, and the relative stability of the two configurations is affected by water self-dissociation, which can only be described properly by FP simulations. Salts in water at extreme conditions play a fundamental role in determining the properties of the Earthʼs mantle constituents. Here the authors shed light on ion-water and ion-ion interactions for NaCl dissolved in water at conditions relevant to the Earthʼs upper mantle by molecular dynamics simulations.
Collapse
Affiliation(s)
- Cunzhi Zhang
- Department of Materials Science and Engineering, COE, Peking University, 100871, Beijing, China
| | - Federico Giberti
- University of Chicago, 5640 S. Ellis Ave., Chicago, IL, 60637, USA
| | - Emre Sevgen
- University of Chicago, 5640 S. Ellis Ave., Chicago, IL, 60637, USA
| | - Juan J de Pablo
- University of Chicago, 5640 S. Ellis Ave., Chicago, IL, 60637, USA.,Materials Science Division, Argonne National Laboratory, Argonne, IL, 60439, USA
| | - Francois Gygi
- University of California Davis, Davis, CA, 95616, USA
| | - Giulia Galli
- University of Chicago, 5640 S. Ellis Ave., Chicago, IL, 60637, USA. .,Materials Science Division, Argonne National Laboratory, Argonne, IL, 60439, USA.
| |
Collapse
|
10
|
Li F, Korotkin IA, Karabasov SA. Rheology of Water Flows Confined between Multilayer Graphene Walls. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2020; 36:5633-5646. [PMID: 32370511 DOI: 10.1021/acs.langmuir.0c01049] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/11/2023]
Abstract
Water confined by hydrophilic materials shows unique transport properties compared to bulk water, thereby offering new opportunities for the development of nanofluidic devices. Recent experimental and numerical studies showed that nanoconfined water undergoes liquid- to solid-phase-like transitions depending on the degree of confinement. In the case of water confined by graphene layers, the van der Waals forces are known to deform the graphene layers, whose bending leads to further nonuniform confinement effects. Despite the extensive studies of nanoconfined water under equilibrium conditions, the interplay between the confinement and rheological water properties, such as viscosity, slip length, and normal stress differences under shear flow conditions, is poorly understood. The current investigation uses a validated all-atom nonequilibrium molecular dynamics model to simultaneously analyze the continuum transport and atomistic structural properties of water in a slit between two moving graphene walls under Couette flow conditions. A range of different slit widths and velocity strain rates are considered. It is shown that under subnanometer confinement, water loses the rotational symmetry of a Newtonian fluid. Under such conditions, water transforms into ice, where the atomistic structure is completely insensitive to the applied shear force and behaves like a frozen slab sliding between the graphene walls. This leads to the shear viscosity increase, although it is not as dramatic as the normal force increase that contributes to the increased friction force reported in previous experimental studies. On the other end of the spectrum, for flows at large velocity strain rates in moderate to large slits between the graphene walls, water is in the liquid state and reveals shear thinning behavior. In this case, water exhibits a constant slip length on the wall, which is typical of liquids in the vicinity of hydrophobic surfaces.
Collapse
Affiliation(s)
- F Li
- The School of Engineering and Materials Science, Queen Mary University of London, Mile End Road, E1 4NS London, United Kingdom
| | - I A Korotkin
- The School of Engineering and Materials Science, Queen Mary University of London, Mile End Road, E1 4NS London, United Kingdom
- Mathematical Sciences, University of Southampton, University Road, SO17 1BJ Southampton, United Kingdom
| | - S A Karabasov
- The School of Engineering and Materials Science, Queen Mary University of London, Mile End Road, E1 4NS London, United Kingdom
| |
Collapse
|
11
|
Silvestri A, Ataman E, Budi A, Stipp SLS, Gale JD, Raiteri P. Wetting Properties of the CO 2-Water-Calcite System via Molecular Simulations: Shape and Size Effects. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2019; 35:16669-16678. [PMID: 31714788 DOI: 10.1021/acs.langmuir.9b02881] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/10/2023]
Abstract
Assessment of the risks and environmental impacts of carbon geosequestration requires knowledge about the wetting behavior of mineral surfaces in the presence of CO2 and the pore fluids. In this context, the interfacial tension (IFT) between CO2 and the aqueous fluid and the contact angle, θ, with the pore mineral surfaces are the two key parameters that control the capillary pressure in the pores of the candidate host rock. Knowledge of these two parameters and their dependence on the local conditions of pressure, temperature, and salinity is essential for the correct prediction of structural and residual trapping. We have performed classical molecular dynamics simulations to predict the CO2-water IFT and the CO2-water-calcite contact angle. The IFT results are consistent with previous simulations, where simple point charge water models have been shown to underestimate the water surface tension, thus affecting the simulated IFT values. When combined with the EPM2 CO2 model, the SPC/Fw water model indeed underestimates the IFT in the low-pressure region at all temperatures studied. On the other hand, at high pressure and low temperature, the IFT is overestimated by ∼5 mN/m. Literature data regarding the CO2/water/calcite contact angle on calcite are contradictory. Using our new set of force field parameters, we performed NVT simulations at 323 K and 20 MPa to calculate the contact angle of a water droplet on the calcite {10.4} surface in a CO2 atmosphere. We performed simulations for both spherical and cylindrical droplet configurations for different initial radii to study the size dependence of the water contact angle on calcite in the presence of CO2. Our results suggest that the contact angle of a cylindrical droplet, is independent of droplet size, for droplets with a radius of 50 Å or more. On the contrary, spherical droplets make a contact angle that is strongly influenced by their size. At the largest size explored in this study, both spherical and cylindrical droplets converge to the same contact angle, 38°, indicating that calcite is strongly wetted by water.
Collapse
Affiliation(s)
- A Silvestri
- Curtin Institute for Computation, The Institute for Geoscience Research (TIGeR), School of Molecular and Life Sciences , Curtin University , PO Box U1987, Perth , WA 6845 , Australia
| | - E Ataman
- Nano-Science Center, Department of Chemistry , University of Copenhagen , Universitetsparken 5 , København Ø DK-2100 , Denmark
| | - A Budi
- Institute for Frontier Materials , Deakin University , Geelong , VIC 3216 , Australia
| | - S L S Stipp
- Department of Physics , Technical University of Denmark , Fysikvej , DK-2800 Kongens Lyngby , Denmark
| | - J D Gale
- Curtin Institute for Computation, The Institute for Geoscience Research (TIGeR), School of Molecular and Life Sciences , Curtin University , PO Box U1987, Perth , WA 6845 , Australia
| | - P Raiteri
- Curtin Institute for Computation, The Institute for Geoscience Research (TIGeR), School of Molecular and Life Sciences , Curtin University , PO Box U1987, Perth , WA 6845 , Australia
| |
Collapse
|
12
|
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
|
13
|
Nikzad M, Azimian AR, Rezaei M, Nikzad S. Water liquid-vapor interface subjected to various electric fields: A molecular dynamics study. J Chem Phys 2017; 147:204701. [DOI: 10.1063/1.4985875] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Mohammadreza Nikzad
- Department of Mechanical Engineering, Khomeinishahr Branch, Islamic Azad University, 84181-48499 Khomeinishahr/Isfahan, Iran
| | - Ahmad Reza Azimian
- Department of Mechanical Engineering, Khomeinishahr Branch, Islamic Azad University, 84181-48499 Khomeinishahr/Isfahan, Iran
| | - Majid Rezaei
- Mechanical Engineering Department, Isfahan University of Technology, 84156-8311 Isfahan, Iran
| | - Safoora Nikzad
- Department of Medical Physics, Hamadan University of Medical Sciences, 65176-19654 Hamadan, Iran
| |
Collapse
|
14
|
Giberti F, Hassanali AA. The excess proton at the air-water interface: The role of instantaneous liquid interfaces. J Chem Phys 2017; 146:244703. [DOI: 10.1063/1.4986082] [Citation(s) in RCA: 32] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Federico Giberti
- Institute for Molecular Engineering, The University of Chicago, 5640 South Ellis Avenue, Chicago, Illinois 60637, USA
| | - Ali A. Hassanali
- Condensed Matter and Statistical Physics, The Abdus Salam International Centre for Theoretical Physics, Strada Costiera 11, 34151 Trieste, Italy
| |
Collapse
|
15
|
Na-Montmorillonite Edge Structure and Surface Complexes: An Atomistic Perspective. MINERALS 2017. [DOI: 10.3390/min7050078] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/17/2022]
|
16
|
Zhao L, Ji J, Tao L, Lin S. Ionic Effects on Supercritical CO2-Brine Interfacial Tensions: Molecular Dynamics Simulations and a Universal Correlation with Ionic Strength, Temperature, and Pressure. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2016; 32:9188-9196. [PMID: 27564433 DOI: 10.1021/acs.langmuir.6b02485] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/06/2023]
Abstract
For geological CO2 storage in deep saline aquifers, the interfacial tension (IFT) between supercritical CO2 and brine is critical for the storage security and design of the storage capacitance. However, currently, no predictive model exists to determine the IFT of supercritical CO2 against complex electrolyte solutions involving various mixed salt species at different concentrations and compositions. In this paper, we use molecular dynamics (MD) simulations to investigate the effect of salt ions on the incremental IFT at the supercritical CO2-brine interface with respect to that at the reference supercritical CO2-water interface. Supercritical CO2-NaCl solution, CO2-CaCl2 solution and CO2-(NaCl+CaCl2) mixed solution systems are simulated at 343 K and 20 MPa under different salinities and salt compositions. We find that the valence of the cations is the primary contributor to the variation in IFT, while the Lennard-Jones potentials for the cations pose a smaller impact on the IFT. Interestingly, the incremental IFT exhibits a general linear correlation with the ionic strength in the above three electrolyte systems, and the slopes are almost identical and independent of the solution types. Based on this finding, a universal predictive formula for IFTs of CO2-complex electrolyte solution systems is established, as a function of ionic strength, temperature, and pressure. The predicted IFTs using the established formula agree perfectly (with a high statistical confidence level of ∼96%) with a wide range of experimental data for CO2 interfacing with different electrolyte solutions, such as those involving MgCl2 and Na2SO4. This work provides an efficient and accurate route to directly predict IFTs in supercritical CO2-complex electrolyte solution systems for practical engineering applications, such as geological CO2 sequestration in deep saline aquifers and other interfacial systems involving complex electrolyte solutions.
Collapse
Affiliation(s)
- Lingling Zhao
- Key Laboratory of Energy Thermal Conversion and Control of Ministry of Education, School of Energy & Environment, Southeast University , Nanjing, Jiangsu 210096, China
| | - Jiayuan Ji
- Key Laboratory of Energy Thermal Conversion and Control of Ministry of Education, School of Energy & Environment, Southeast University , Nanjing, Jiangsu 210096, China
| | - Lu Tao
- Key Laboratory of Energy Thermal Conversion and Control of Ministry of Education, School of Energy & Environment, Southeast University , Nanjing, Jiangsu 210096, China
| | - Shangchao Lin
- Department of Mechanical Engineering, Materials Science & Engineering Program, FAMU-FSU College of Engineering, Florida State University , Tallahassee, Florida 32310, United States
| |
Collapse
|
17
|
|
18
|
Ta TD, Tieu AK, Zhu H, Zhu Q, Kosasih PB, Zhang J, Deng G. Tribological Behavior of Aqueous Copolymer Lubricant in Mixed Lubrication Regime. ACS APPLIED MATERIALS & INTERFACES 2016; 8:5641-5652. [PMID: 26828119 DOI: 10.1021/acsami.5b10905] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/05/2023]
Abstract
Although a number of experiments have been attempted to investigate the lubrication of aqueous copolymer lubricant, which is applied widely in metalworking operations, a comprehensive theoretical investigation at atomistic level is still lacking. This study addresses the influence of loading pressure and copolymer concentration on the structural properties and tribological performance of aqueous copolymer solution of poly(propylene oxide)-poly(ethylene oxide)-poly(propylene oxide) (PPO-PEO-PPO) at mixed lubrication using a molecular dynamic (MD) simulation. An effective interfacial potential, which has been derived from density functional theory (DFT) calculations, was employed for the interactions between the fluid's molecules and iron surface. The simulation results have indicated that the triblock copolymer is physisorption on iron surface. Under confinement by iron surfaces, the copolymer molecules form lamellar structure in aqueous solution and behave differently from its bulk state. The lubrication performance of aqueous copolymer lubricant increases with concentration, but the friction reduction is insignificant at high loading pressure. Additionally, the plastic deformation of asperity is dependent on both copolymer concentration and loading pressure, and the wear behavior shows a linear dependence of friction force on the number of transferred atoms between contacting asperities.
Collapse
Affiliation(s)
- Thi D Ta
- School of Mechanical, Materials and Mechatronic Engineering, Faculty of Engineering and Information Sciences (EIS), University of Wollongong , Northfield Avenue, Wollongong, New South Wales 2522, Australia
| | - A Kiet Tieu
- School of Mechanical, Materials and Mechatronic Engineering, Faculty of Engineering and Information Sciences (EIS), University of Wollongong , Northfield Avenue, Wollongong, New South Wales 2522, Australia
| | - Hongtao Zhu
- School of Mechanical, Materials and Mechatronic Engineering, Faculty of Engineering and Information Sciences (EIS), University of Wollongong , Northfield Avenue, Wollongong, New South Wales 2522, Australia
| | - Qiang Zhu
- School of Mechanical, Materials and Mechatronic Engineering, Faculty of Engineering and Information Sciences (EIS), University of Wollongong , Northfield Avenue, Wollongong, New South Wales 2522, Australia
| | - Prabouno B Kosasih
- School of Mechanical, Materials and Mechatronic Engineering, Faculty of Engineering and Information Sciences (EIS), University of Wollongong , Northfield Avenue, Wollongong, New South Wales 2522, Australia
| | - Jie Zhang
- School of Mechanical, Materials and Mechatronic Engineering, Faculty of Engineering and Information Sciences (EIS), University of Wollongong , Northfield Avenue, Wollongong, New South Wales 2522, Australia
| | - Guanyu Deng
- School of Mechanical, Materials and Mechatronic Engineering, Faculty of Engineering and Information Sciences (EIS), University of Wollongong , Northfield Avenue, Wollongong, New South Wales 2522, Australia
| |
Collapse
|
19
|
Stinson JL, Kathmann SM, Ford IJ. A classical reactive potential for molecular clusters of sulphuric acid and water. Mol Phys 2015. [DOI: 10.1080/00268976.2015.1090027] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/23/2022]
Affiliation(s)
- Jake L. Stinson
- Department of Physics and Astronomy and London Centre for Nanotechnology, University College London, London, United Kingdom
| | - Shawn M. Kathmann
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, WA, United States
| | - Ian J. Ford
- Department of Physics and Astronomy and London Centre for Nanotechnology, University College London, London, United Kingdom
| |
Collapse
|
20
|
Jiang H, Mester Z, Moultos OA, Economou IG, Panagiotopoulos AZ. Thermodynamic and Transport Properties of H2O + NaCl from Polarizable Force Fields. J Chem Theory Comput 2015; 11:3802-10. [DOI: 10.1021/acs.jctc.5b00421] [Citation(s) in RCA: 53] [Impact Index Per Article: 5.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Hao Jiang
- Department
of Chemical and Biological Engineering, Princeton University, Princeton, New Jersey 08544, United States
| | - Zoltan Mester
- Department
of Chemical and Biological Engineering, Princeton University, Princeton, New Jersey 08544, United States
| | - Othonas A. Moultos
- Chemical Engineering Program, Texas A&M University at Qatar, P.O. Box 23874, Doha, Qatar
| | - Ioannis G. Economou
- Chemical Engineering Program, Texas A&M University at Qatar, P.O. Box 23874, Doha, Qatar
| | | |
Collapse
|
21
|
Shvab I, Sadus RJ. Thermophysical properties of supercritical water and bond flexibility. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2015; 92:012124. [PMID: 26274141 DOI: 10.1103/physreve.92.012124] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/15/2015] [Indexed: 06/04/2023]
Abstract
Molecular dynamics results are reported for the thermodynamic properties of supercritical water using examples of both rigid (TIP4P/2005) and flexible (TIP4P/2005f) transferable interaction potentials. Data are reported for pressure, isochoric and isobaric heat capacities, the thermal expansion coefficient, isothermal and adiabatic compressibilities, Joule-Thomson coefficient, speed of sound, self-diffusion coefficient, viscosities, and thermal conductivity. Many of these properties have unusual behavior in the supercritical phase such as maximum and minimum values. The effectiveness of bond flexibility on predicting these properties is determined by comparing the results to experimental data. The influence of the intermolecular potential on these properties is both variable and state point dependent. In the vicinity of the critical density, the rigid and flexible potentials yield very different values for the compressibilities, heat capacities, and thermal expansion coefficient, whereas the self-diffusion coefficient, viscosities, and thermal conductivities are much less potential dependent. Although the introduction of bond flexibility is a computationally expedient way to improve the accuracy of an intermolecular potential, it can be counterproductive in some cases and it is not an adequate replacement for incorporating the effects of polarization.
Collapse
Affiliation(s)
- I Shvab
- School of Civil Engineering and Geosciences, Newcastle University, Newcastle upon Tyne, NE1 7RU, United Kingdom
| | - Richard J Sadus
- Centre for Molecular Simulation, Swinburne University of Technology, PO Box 218, Hawthorn, Victoria 3122, Australia
| |
Collapse
|
22
|
Orozco GA, Moultos OA, Jiang H, Economou IG, Panagiotopoulos AZ. Molecular simulation of thermodynamic and transport properties for the H2O+NaCl system. J Chem Phys 2014; 141:234507. [DOI: 10.1063/1.4903928] [Citation(s) in RCA: 51] [Impact Index Per Article: 5.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022] Open
Affiliation(s)
- Gustavo A. Orozco
- Department of Chemical and Biological Engineering, Princeton University, Princeton, New Jersey 08544, USA
| | - Othonas A. Moultos
- Chemical Engineering Program, Texas A&M University at Qatar, P.O. Box 23874, Doha, Qatar
| | - Hao Jiang
- Department of Chemical and Biological Engineering, Princeton University, Princeton, New Jersey 08544, USA
| | - Ioannis G. Economou
- Chemical Engineering Program, Texas A&M University at Qatar, P.O. Box 23874, Doha, Qatar
| | | |
Collapse
|
23
|
Orozco GA, Economou IG, Panagiotopoulos AZ. Optimization of Intermolecular Potential Parameters for the CO2/H2O Mixture. J Phys Chem B 2014; 118:11504-11. [DOI: 10.1021/jp5067023] [Citation(s) in RCA: 27] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/19/2023]
Affiliation(s)
- Gustavo A. Orozco
- Department
of Chemical and Biological Engineering, Princeton University, Princeton, New Jersey 08544, United States
| | | | | |
Collapse
|
24
|
Kamath G, Deshmukh SA, Sankaranarayanan SKRS. Comparison of select polarizable and non-polarizable water models in predicting solvation dynamics of water confined between MgO slabs. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2013; 25:305003. [PMID: 23819970 DOI: 10.1088/0953-8984/25/30/305003] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/02/2023]
Abstract
We present a molecular dynamics simulation study in which we compare and contrast the performance of a polarizable shell water potential model and non-polarizable water force field-extended simple point charge (SPC/EF) model in predicting the solvation dynamics of confined water molecules sandwiched between MgO(100) slabs. Structural features based on radial distribution functions, atomic density profiles, adsorption patterns, orientational ordering and dynamical correlations such as diffusional characteristics, hydrogen bonding lifetimes and residence probabilities are used as metrics for comparison. The simulations yield significant ordering of water molecules in the two layers adjacent to the oxide interface and the extent of ordering decreases with increasing distance from the oxide-water interface. These results elucidate that the dependence of local ordering and solvation dynamics on the molecular geometry and charge distribution, observed for typical three- and four-site water models, is generally lost for confined water if polarization is explicitly included. While the interfacial water structure predicted by the polarizable and non-polarizable models are similar, the confinement and interface proximity effects on the solvation dynamics are seen to be more pronounced for polarizable water models in comparison to non-polarizable ones. The study also shows that the polarizable water model over predicts the orientational order and under predicts the transport properties of confined water. In addition, analysis of the orientational preferences and hydrogen bonding characteristics of water near oxide interfaces suggests a higher degree of tetrahedral disorder in the polarizable shell compared to the non-polarizable SPC/E flexible model. The origin of the differences in solvation behavior of confined water between oxide slabs is analyzed based on the energetic contributions of the dispersive and electrostatic terms in the two force fields. Our findings suggest some new considerations regarding the role of polarization terms in predicting confinement and interface proximity effects that may guide future development of reliable polarizable water models for confined liquids.
Collapse
Affiliation(s)
- Ganesh Kamath
- Department of Chemistry, University of Missouri-Columbia, Columbia 65211, USA
| | | | | |
Collapse
|
25
|
Raabe G, Sadus RJ. Molecular dynamics simulation of the effect of bond flexibility on the transport properties of water. J Chem Phys 2012; 137:104512. [DOI: 10.1063/1.4749382] [Citation(s) in RCA: 50] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022] Open
|
26
|
Abstract
A new flexible water model, TIP4P/2005f, is developed. The idea was to add intramolecular degrees of freedom to the successful rigid model TIP4P/2005 in order to try to improve the predictions for some properties, and to enable the calculation of new ones. The new model incorporates flexibility by means of a Morse potential for the bond stretching and a harmonic term for the angle bending. The parameters have been fitted to account for the peaks of the infrared spectrum of liquid water and to produce an averaged geometry close to that of TIP4P/2005. As for the intermolecular interactions, only a small change in the σ parameter of the Lennard-Jones potential has been introduced. The overall predictions are very close to those of TIP4P/2005. This ensures that the new model may be used with the same confidence as its predecessor in studies where a flexible model is advisable.
Collapse
Affiliation(s)
- Miguel A González
- Departamento de Química Física, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | | |
Collapse
|
27
|
Deshmukh SA, Sankaranarayanan SKRS. Atomic scale characterization of interfacial water near an oxide surface using molecular dynamics simulations. Phys Chem Chem Phys 2012; 14:15593-605. [DOI: 10.1039/c2cp42308a] [Citation(s) in RCA: 27] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
28
|
WANG JUN, ZENG XIAOCHENG. COMPUTER SIMULATION OF LIQUID–VAPOR INTERFACIAL TENSION: LENNARD-JONES FLUID AND WATER REVISITED. JOURNAL OF THEORETICAL & COMPUTATIONAL CHEMISTRY 2011. [DOI: 10.1142/s0219633609005027] [Citation(s) in RCA: 12] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
Abstract
We review several commonly used simulation methods for computing liquid–vapor surface tension and associated theoretical treatments of the long-range correction for inhomogeneous systems. Prototype model systems considered in this review are the Lennard-Jones (LJ) fluid and the SPC/E model water. In addition, we examine a variety of factors that can affect calculation of the surface tension γ via the mechanical approach (i.e. using either KB or IK method). It is found that for the LJ fluid, the size of simulation box and the number of particles in the system can have notable effects on the computed surface tension. For SPC/E water, the Ewald parameters can influence computed surface tensions (γ) as well, e.g., very small Ewald parameters tend to overestimate γ. It is also found that the IK method consistently gives γ that are 0.6 - 0.9 mN/m greater than γ computed based on the KB method. When computing the first reciprocal–space contribution to the surface tension, the Ghoufi's strategy gives rise to more sensible profile of pressure difference PN(z)-PT(z) than the Alejandre's strategy although both strategies result in nearly the same average surface tension through the integration of PN(z)-PT(z).
Collapse
Affiliation(s)
- JUN WANG
- Department of Chemistry, University of Nebraska-Lincoln, Lincoln, Nebraska 68588, USA
| | - XIAO CHENG ZENG
- Department of Chemistry, University of Nebraska-Lincoln, Lincoln, Nebraska 68588, USA
| |
Collapse
|
29
|
Raabe G, Sadus RJ. Molecular dynamics simulation of the dielectric constant of water: The effect of bond flexibility. J Chem Phys 2011; 134:234501. [DOI: 10.1063/1.3600337] [Citation(s) in RCA: 71] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022] Open
|
30
|
Zhao L, Lin S, Mendenhall JD, Yuet PK, Blankschtein D. Molecular Dynamics Investigation of the Various Atomic Force Contributions to the Interfacial Tension at the Supercritical CO2–Water Interface. J Phys Chem B 2011; 115:6076-87. [DOI: 10.1021/jp201190g] [Citation(s) in RCA: 43] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Lingling Zhao
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
- School of Energy & Environment, Southeast University, Nanjing, China 210096
| | - Shangchao Lin
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
- Department of Mechanical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Jonathan D. Mendenhall
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Pak K. Yuet
- Department of Process Engineering and Applied Science, Dalhousie University, Canada B3J 2X4
| | - Daniel Blankschtein
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| |
Collapse
|
31
|
Sakamaki R, Sum AK, Narumi T, Yasuoka K. Molecular dynamics simulations of vapor/liquid coexistence using the nonpolarizable water models. J Chem Phys 2011; 134:124708. [DOI: 10.1063/1.3574038] [Citation(s) in RCA: 51] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
32
|
da Silva JAB, Moreira FGB, dos Santos VML, Longo RL. Hydrogen bond networks in water and methanol with varying interaction strengths. Phys Chem Chem Phys 2011; 13:593-603. [DOI: 10.1039/c0cp01204a] [Citation(s) in RCA: 30] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
33
|
Vega C, Abascal JLF. Simulating water with rigid non-polarizable models: a general perspective. Phys Chem Chem Phys 2011; 13:19663-88. [DOI: 10.1039/c1cp22168j] [Citation(s) in RCA: 658] [Impact Index Per Article: 50.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
34
|
Yuet PK, Blankschtein D. Molecular Dynamics Simulation Study of Water Surfaces: Comparison of Flexible Water Models. J Phys Chem B 2010; 114:13786-95. [DOI: 10.1021/jp1067022] [Citation(s) in RCA: 86] [Impact Index Per Article: 6.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Pak K. Yuet
- Department of Process Engineering and Applied Science, Dalhousie University, P.O. Box 1000, Halifax, Nova Scotia, Canada B3J 2X4, and Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Daniel Blankschtein
- Department of Process Engineering and Applied Science, Dalhousie University, P.O. Box 1000, Halifax, Nova Scotia, Canada B3J 2X4, and Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| |
Collapse
|
35
|
Alejandre J, Chapela GA. The surface tension of TIP4P/2005 water model using the Ewald sums for the dispersion interactions. J Chem Phys 2010; 132:014701. [PMID: 20078174 DOI: 10.1063/1.3279128] [Citation(s) in RCA: 72] [Impact Index Per Article: 5.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The liquid-vapor phase equilibria and surface tension of the TIP4P/2005 water model is obtained by using the Ewald summation method to determine the long range Lennard-Jones and electrostatic interactions. The method is implemented in a straightforward manner into standard simulation programs. The computational cost of using Ewald sums in dispersion interactions of water is estimated in direct simulation of interfaces. The results of this work at 300 K show a dramatic change in surface tension with an oscillatory behavior for surface areas smaller than 5x5sigma(2), where sigma is the Lennard-Jones oxygen diameter. The amplitude of such oscillations substantially decreases with temperature. Finite size effects are less important on coexisting densities. Phase equilibria and interfacial properties can be determined using a small number of water molecules; their fluctuations are around the same size of simulation error at all temperatures, even in systems where the interfaces are separated a few molecular diameters only. The difference in surface tension of this work compared to the results of other authors is not significant (on the contrary, there is a good agreement). What should be stressed is the different and more consistent approach to obtain the surface tension using the Ewald sums for dispersion interactions. There are two relevant aspects at the interface: An adsorption of water molecules is observed at small surface areas and its thickness systematically increases with system size.
Collapse
Affiliation(s)
- José Alejandre
- Departamento de Química, Universidad Autónoma Metropolitana-Iztapalapa, Av. San Rafael Atlixco 186, Col. Vicentina, México D.F. 09340, Mexico.
| | | |
Collapse
|
36
|
Alejandre J, Chapela GA. The dipole moment distribution on water is improved by using large flexibility and large bending angle. Mol Phys 2010. [DOI: 10.1080/00268970903563451] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/19/2022]
Affiliation(s)
| | - Gustavo A. Chapela
- b Departamento de Física , Universidad Autónoma Metropolitana-Iztapalapa , Av. San Rafael Atlixco 186, Col. Vicentina, 09340 México D.F., México
| |
Collapse
|
37
|
González BS, Noya EG, Vega C, Sesé LM. Nuclear Quantum Effects in Water Clusters: The Role of the Molecular Flexibility. J Phys Chem B 2010; 114:2484-92. [DOI: 10.1021/jp910770y] [Citation(s) in RCA: 28] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Briesta S. González
- Departamento Química Física I, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - Eva G. Noya
- Departamento Química Física I, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - Carlos Vega
- Departamento Química Física I, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - Luis M. Sesé
- Departamento Ciencias y Técnicas Fisicoquímicas, Facultad de Ciencias, UNED, Paseo Senda del Rey 9, 28040 Madrid, Spain
| |
Collapse
|
38
|
Luis D, López-Lemus J, Mayorga M, Romero-Salazar L. Performance of rigid water models in the phase transition of clathrates†. MOLECULAR SIMULATION 2010. [DOI: 10.1080/08927020903096072] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/20/2022]
|
39
|
Lopes PEM, Roux B, MacKerell AD. Molecular modeling and dynamics studies with explicit inclusion of electronic polarizability. Theory and applications. Theor Chem Acc 2009; 124:11-28. [PMID: 20577578 PMCID: PMC2888514 DOI: 10.1007/s00214-009-0617-x] [Citation(s) in RCA: 265] [Impact Index Per Article: 17.7] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/24/2022]
Abstract
A current emphasis in empirical force fields is on the development of potential functions that explicitly treat electronic polarizability. In the present article, the commonly used methodologies for modelling electronic polarization are presented along with an overview of selected application studies. Models presented include induced point-dipoles, classical Drude oscillators, and fluctuating charge methods. The theoretical background of each method is followed by an introduction to extended Langrangian integrators required for computationally tractable molecular dynamics simulations using polarizable force fields. The remainder of the review focuses on application studies using these methods. Emphasis is placed on water models, for which numerous examples exist, with a more thorough discussion presented on the recently published models associated with the Drude-based CHARMM and the AMOEBA force fields. The utility of polarizable models for the study of ion solvation is then presented followed by an overview of studies of small molecules (e.g. CCl(4), alkanes, etc) and macromolecule (proteins, nucleic acids and lipid bilayers) application studies. The review is written with the goal of providing a general overview of the current status of the field and to facilitate future application and developments.
Collapse
Affiliation(s)
- Pedro E. M. Lopes
- Department of Pharmaceutical Sciences, School of Pharmacy, University of Maryland, 20 Penn Street, Baltimore, MD 21230, USA
| | - Benoit Roux
- Institute of Molecular Pediatric Sciences, Gordon Center for Integrative Science, University of Chicago 929 E. 57th St. Chicago, IL 60637
| | - Alexander D. MacKerell
- Department of Pharmaceutical Sciences, School of Pharmacy, University of Maryland, 20 Penn Street, Baltimore, MD 21230, USA
| |
Collapse
|
40
|
Habershon S, Markland TE, Manolopoulos DE. Competing quantum effects in the dynamics of a flexible water model. J Chem Phys 2009; 131:024501. [DOI: 10.1063/1.3167790] [Citation(s) in RCA: 373] [Impact Index Per Article: 24.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
41
|
Alejandre J, Chapela GA, Bresme F, Hansen JP. The short range anion-H interaction is the driving force for crystal formation of ions in water. J Chem Phys 2009; 130:174505. [DOI: 10.1063/1.3124184] [Citation(s) in RCA: 60] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
42
|
Vega C, Abascal JLF, Conde MM, Aragones JL. What ice can teach us about water interactions: a critical comparison of the performance of different water models. Faraday Discuss 2009; 141:251-76; discussion 309-46. [DOI: 10.1039/b805531a] [Citation(s) in RCA: 328] [Impact Index Per Article: 21.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/21/2023]
|
43
|
Mendoza FN, López-Lemus J, Chapela GA, Alejandre J. The Wolf method applied to the liquid-vapor interface of water. J Chem Phys 2008; 129:024706. [DOI: 10.1063/1.2948951] [Citation(s) in RCA: 27] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|