1
|
Johnson E, Haussener S. Contrasting Views of the Electric Double Layer in Electrochemical CO 2 Reduction: Continuum Models vs Molecular Dynamics. THE JOURNAL OF PHYSICAL CHEMISTRY. C, NANOMATERIALS AND INTERFACES 2024; 128:10450-10464. [PMID: 38957368 PMCID: PMC11215773 DOI: 10.1021/acs.jpcc.4c03469] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 05/26/2024] [Revised: 06/05/2024] [Accepted: 06/05/2024] [Indexed: 07/04/2024]
Abstract
In the field of electrochemical CO2 reduction, both continuum models and molecular dynamics (MD) models have been used to understand the electric double layer (EDL). MD often focuses on the region within a few nm of the electrode, while continuum models can span up to the device level (cm). Still, both methods model the EDL, and for a cohesive picture of the CO2 electrolysis system, the two methods should agree in the regions where they overlap length scales. To this end, we make a direct comparison between state-of-the-art continuum models and classical MD simulations under the conditions of CO2 reduction on a Ag electrode. For continuum modeling, this includes the Poisson-Nernst-Planck formulation with steric (finite ion size) effects, and in MD the electrode is modeled with the constant potential method. The comparison yields numerous differences between the two modeling methods. MD shows cations forming two adsorbed layers, including a fully hydrated outer layer and a partial hydration layer closer to the electrode surface. The strength of the inner adsorbed layer increases with cation size (Li+ < Na+ < K+ < Cs+) and with more negative applied potentials. Continuum models that include steric effects predict CO2 to be mostly excluded within 1 nm of the cathode due to tightly packed cations, yet we find little evidence to support these predictions from the MD results. In fact, MD shows that the concentration of CO2 increases within a few Å of the cathode surface due to interactions with the Ag electrode, a factor not included in continuum models. The EDL capacitance is computed from the MD results, showing values in the range of 7-9 μF cm-2, irrespective of the electrolyte concentration, cation identity, or applied potential. The direct comparison between the two modeling methods is meant to show the areas of agreement and disagreement between the two views of the EDL, so as to improve and better align these models.
Collapse
Affiliation(s)
- Evan Johnson
- Laboratory of Renewable Energy
Science and Engineering, École Polytechnique
Fédérale de Lausanne, Station 9, 1015 Lausanne, Switzerland
| | - Sophia Haussener
- Laboratory of Renewable Energy
Science and Engineering, École Polytechnique
Fédérale de Lausanne, Station 9, 1015 Lausanne, Switzerland
| |
Collapse
|
2
|
Khammari A, Jeon J, Cho M. Solvation Structures and Transport Mechanisms of Cations in Water-in-Bisalt Electrolytes for High-Concentration Lithium-Ion Batteries Revealed by Molecular Dynamics Simulations. J Phys Chem B 2024; 128:5735-5745. [PMID: 38820050 DOI: 10.1021/acs.jpcb.4c02001] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/02/2024]
Abstract
The development of safe and cost-effective electrolytes for rechargeable batteries is currently underway. While water-based electrolytes hold promise, their restricted electrochemical stability window poses a challenge. Combining multiple ionic species emerges as a promising strategy to broaden this stability window and optimize Li-ion battery performance. This study focuses on dual-cation electrolytes, which blend lithium and potassium acetates to enhance the electrochemical characteristics of the solution at high concentrations. We investigated the solvation structure of each ion and its interactions on a molecular level. Our analysis reveals that ion clusters and aggregates are formed through shared acetate and water molecules at high salt concentrations. Furthermore, the residence time analyses of atom pairs indicate that cations diffuse in vehicular mode at low concentrations. In contrast, they switch to a structural mode at high concentrations due to diminishing water content. This study offers a comprehensive model for exploring diverse solvation structures of cations and gaining insights into their diffusion mechanisms within water-in-bisalt electrolytes for aqueous Li-ion batteries.
Collapse
Affiliation(s)
- Anahita Khammari
- Center for Molecular Spectroscopy and Dynamics, Institute for Basic Science (IBS), Seoul 02841, Republic of Korea
- Korea University, Seoul 02841, Republic of Korea
| | - Jonggu Jeon
- Center for Molecular Spectroscopy and Dynamics, Institute for Basic Science (IBS), Seoul 02841, Republic of Korea
- Korea University, Seoul 02841, Republic of Korea
| | - Minhaeng Cho
- Center for Molecular Spectroscopy and Dynamics, Institute for Basic Science (IBS), Seoul 02841, Republic of Korea
- Department of Chemistry, Korea University, Seoul 02841, Republic of Korea
| |
Collapse
|
3
|
Habibi P, Polat HM, Blazquez S, Vega C, Dey P, Vlugt TJH, Moultos OA. Accurate Free Energies of Aqueous Electrolyte Solutions from Molecular Simulations with Non-polarizable Force Fields. J Phys Chem Lett 2024; 15:4477-4485. [PMID: 38634502 PMCID: PMC11057036 DOI: 10.1021/acs.jpclett.4c00428] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/09/2024] [Revised: 04/11/2024] [Accepted: 04/12/2024] [Indexed: 04/19/2024]
Abstract
Non-polarizable force fields fail to accurately predict free energies of aqueous electrolytes without compromising the predictive ability for densities and transport properties. A new approach is presented in which (1) TIP4P/2005 water and scaled charge force fields are used to describe the interactions in the liquid phase and (2) an additional Effective Charge Surface (ECS) is used to compute free energies at zero additional computational expense. The ECS is obtained using a single temperature-independent charge scaling parameter per species. Thereby, the chemical potential of water and the free energies of hydration of various aqueous salts (e.g., NaCl and LiCl) are accurately described (deviations less than 5% from experiments), in sharp contrast to calculations where the ECS is omitted (deviations larger than 20%). This approach enables accurate predictions of free energies of aqueous electrolyte solutions using non-polarizable force fields, without compromising liquid-phase properties.
Collapse
Affiliation(s)
- Parsa Habibi
- Engineering
Thermodynamics, Process & Energy Department, Faculty of Mechanical
Engineering, Delft University of Technology, Leeghwaterstraat 39, 2628 CB Delft, Netherlands
- Department
of Materials Science and Engineering, Faculty of Mechanical Engineering, Delft University of Technology, Mekelweg 2, 2628 CD Delft, Netherlands
| | - H. Mert Polat
- Engineering
Thermodynamics, Process & Energy Department, Faculty of Mechanical
Engineering, Delft University of Technology, Leeghwaterstraat 39, 2628 CB Delft, Netherlands
| | - Samuel Blazquez
- Departamento
de Química Física, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - Carlos Vega
- Departamento
de Química Física, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - Poulumi Dey
- Department
of Materials Science and Engineering, Faculty of Mechanical Engineering, Delft University of Technology, Mekelweg 2, 2628 CD Delft, Netherlands
| | - Thijs J. H. Vlugt
- Engineering
Thermodynamics, Process & Energy Department, Faculty of Mechanical
Engineering, Delft University of Technology, Leeghwaterstraat 39, 2628 CB Delft, Netherlands
| | - Othonas A. Moultos
- Engineering
Thermodynamics, Process & Energy Department, Faculty of Mechanical
Engineering, Delft University of Technology, Leeghwaterstraat 39, 2628 CB Delft, Netherlands
| |
Collapse
|
4
|
Wagemann E, Kim NY, Mitra SK. Cs + conductance in graphene membranes with Ångström-scale pores: the role of pore entrance geometry. Phys Chem Chem Phys 2024; 26:11311-11319. [PMID: 38444318 DOI: 10.1039/d4cp00400k] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/07/2024]
Abstract
Recently, 2D materials, such as graphene, have been successfully implemented as artificial conduits of molecular sizes. The extreme precision with which these structures can be fabricated provides an unprecedented framework for the development of highly specific and efficient devices. In this work, we study the electrophoretic transport of Cs+ ions in a graphene membrane with effective pore heights of 3.4 Å by conducting molecular dynamics simulations. The entrance of the pore is systematically modified to investigate the effect of pore geometry on ionic conductance. Simulation results suggest a significant correlation between ionic conductance and entrance geometry, with a variation of the conductance up to 100% across the studied cases. To explain the observed correlation, two mechanisms involving an intimate relationship between ion dehydration and edge functional groups are proposed. The present study provides theoretical insights that can aid the design of graphene-based membranes with tunable ionic transport properties.
Collapse
Affiliation(s)
- Enrique Wagemann
- Departamento de Ingeniería Mecanica, Facultad de Ingeniería, Universidad de Concepción, Concepción 4070409, Chile
| | - Na Young Kim
- Institute for Quantum Computing, University of Waterloo, 200 University Ave. West, Waterloo, Ontario N2L 3G1, Canada
- Department of Electrical and Computer Engineering, University of Waterloo, 200 University Ave. West, Waterloo, Ontario N2L 3G1, Canada
- Waterloo Institute for Nanotechnology, University of Waterloo, 200 University Ave. West, Waterloo, Ontario N2L 3G1, Canada
| | - Sushanta K Mitra
- Micro & Nano-Scale Transport Laboratory, Waterloo Institute for Nanotechnology, Department of Mechanical and Mechatronics Engineering, University of Waterloo, Waterloo, Ontario N2L 3G1, Canada.
| |
Collapse
|
5
|
van Rooijen WA, Habibi P, Xu K, Dey P, Vlugt TJH, Hajibeygi H, Moultos OA. Interfacial Tensions, Solubilities, and Transport Properties of the H 2/H 2O/NaCl System: A Molecular Simulation Study. JOURNAL OF CHEMICAL AND ENGINEERING DATA 2024; 69:307-319. [PMID: 38352074 PMCID: PMC10859954 DOI: 10.1021/acs.jced.2c00707] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 11/12/2022] [Accepted: 12/23/2022] [Indexed: 02/16/2024]
Abstract
Data for several key thermodynamic and transport properties needed for technologies using hydrogen (H2), such as underground H2 storage and H2O electrolysis are scarce or completely missing. Force field-based Molecular Dynamics (MD) and Continuous Fractional Component Monte Carlo (CFCMC) simulations are carried out in this work to cover this gap. Extensive new data sets are provided for (a) interfacial tensions of H2 gas in contact with aqueous NaCl solutions for temperatures of (298 to 523) K, pressures of (1 to 600) bar, and molalities of (0 to 6) mol NaCl/kg H2O, (b) self-diffusivities of infinitely diluted H2 in aqueous NaCl solutions for temperatures of (298 to 723) K, pressures of (1 to 1000) bar, and molalities of (0 to 6) mol NaCl/kg H2O, and (c) solubilities of H2 in aqueous NaCl solutions for temperatures of (298 to 363) K, pressures of (1 to 1000) bar, and molalities of (0 to 6) mol NaCl/kg H2O. The force fields used are the TIP4P/2005 for H2O, the Madrid-2019 and the Madrid-Transport for NaCl, and the Vrabec and Marx for H2. Excellent agreement between the simulation results and available experimental data is found with average deviations lower than 10%.
Collapse
Affiliation(s)
- W. A. van Rooijen
- Reservoir
Engineering, Geoscience and Engineering Department, Faculty of Civil
Engineering and Geosciences, Delft University
of Technology, Stevinweg 1, 2628CN, Delft, The Netherlands
| | - P. Habibi
- Engineering
Thermodynamics, Process and Energy Department, Faculty of Mechanical,
Maritime and Materials Engineering, Delft
University of Technology, Leeghwaterstraat 39, 2628CB, Delft, The Netherlands
- Department
of Materials Science and Engineering, Faculty of Mechanical, Maritime
and Materials Engineering, Delft University
of Technology, Mekelweg
2, 2628CD, Delft, The Netherlands
| | - K. Xu
- Department
of Materials Science and Engineering, Faculty of Mechanical, Maritime
and Materials Engineering, Delft University
of Technology, Mekelweg
2, 2628CD, Delft, The Netherlands
| | - P. Dey
- Department
of Materials Science and Engineering, Faculty of Mechanical, Maritime
and Materials Engineering, Delft University
of Technology, Mekelweg
2, 2628CD, Delft, The Netherlands
| | - T. J. H. Vlugt
- Engineering
Thermodynamics, Process and Energy Department, Faculty of Mechanical,
Maritime and Materials Engineering, Delft
University of Technology, Leeghwaterstraat 39, 2628CB, Delft, The Netherlands
| | - H. Hajibeygi
- Reservoir
Engineering, Geoscience and Engineering Department, Faculty of Civil
Engineering and Geosciences, Delft University
of Technology, Stevinweg 1, 2628CN, Delft, The Netherlands
| | - O. A. Moultos
- Engineering
Thermodynamics, Process and Energy Department, Faculty of Mechanical,
Maritime and Materials Engineering, Delft
University of Technology, Leeghwaterstraat 39, 2628CB, Delft, The Netherlands
| |
Collapse
|
6
|
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
|
7
|
Hosseni A, Ashbaugh HS. Osmotic Force Balance Evaluation of Aqueous Electrolyte Osmotic Pressures and Chemical Potentials. J Chem Theory Comput 2023; 19:8826-8838. [PMID: 37978934 PMCID: PMC10720338 DOI: 10.1021/acs.jctc.3c00982] [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/07/2023] [Revised: 10/26/2023] [Accepted: 10/27/2023] [Indexed: 11/19/2023]
Abstract
Concentrated aqueous salt solutions are ubiquitous in problems of biological and environmental relevance. The development of accurate force fields that capture the interactions between dissolved species in solution is crucial to simulating these systems to gain molecular insights into the underlying processes under saline conditions. The osmotic pressure is a relatively simple thermodynamic property connecting the experimental and simulation measurements of the associative properties of the ions in solution. Milner [C. Gillespie and S. T. Milner, Soft Matter, 16, 9816 (2020)] proposed a simulation approach to evaluate the osmotic pressures of salts in solution by applying a restraint potential to the ions alone in solution and determining the resulting pressure required to balance that potential, referred to here as the osmotic force balance. Here, we expand Milner's approach, demonstrating that the chemical potentials of the salts in solution as a function of concentration can be fitted to the concentration profiles determined from simulation, additionally providing an analytical expression for the osmotic pressure. This approach is used to determine the osmotic pressures of 15 alkali halide salts in water from simulations. The cross interactions between cations and anions in solution are subsequently optimized to capture their experimental osmotic pressures.
Collapse
Affiliation(s)
- Alireza Hosseni
- Department of Chemical and
Biomolecular Engineering, Tulane University, New Orleans, Louisiana 70118, United States
| | - Henry S. Ashbaugh
- Department of Chemical and
Biomolecular Engineering, Tulane University, New Orleans, Louisiana 70118, United States
| |
Collapse
|
8
|
Almeida NMS, Bali SK, James D, Wang C, Wilson AK. Binding of Per- and Polyfluoroalkyl Substances (PFAS) to the PPARγ/RXRα-DNA Complex. J Chem Inf Model 2023; 63:7423-7443. [PMID: 37990410 DOI: 10.1021/acs.jcim.3c01384] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/23/2023]
Abstract
Nuclear receptors are the fundamental building blocks of gene expression regulation and the focus of many drug targets. While binding to DNA, nuclear receptors act as transcription factors, governing a multitude of functions in the human body. Peroxisome proliferator-activator receptor γ (PPARγ) and the retinoid X receptor α (RXRα) form heterodimers with unique properties and have a primordial role in insulin sensitization. This PPARγ/RXRα heterodimer has been shown to be impacted by per- and polyfluoroalkyl substances (PFAS) and linked to a variety of significant health conditions in humans. Herein, a selection of the most common PFAS (legacy and emerging) was studied utilizing molecular dynamics simulations for PPARγ/RXRα. The local and global structural effects of PFAS binding on the known ligand binding pockets of PPARγ and RXRα as well as the DNA binding domain (DBD) of RXRα were inspected. The binding free energies were predicted computationally and were compared between the different binding pockets. In addition, two electronic structure approaches were utilized to model the interaction of PFAS within the DNA binding domain, density functional theory (DFT) and domain-based pair natural orbital coupled cluster with perturbative triples (DLPNO-CCSD(T)) approaches, with implicit solvation. Residue decomposition and hydrogen-bonding analysis were also performed, detailing the role of prominent residues in molecular recognition. The role of l-carnitine is explored as a potential in vivo remediation strategy for PFAS interaction with the PPARγ/RXRα heterodimer. In this work, it was found that PFAS can bind and act as agonists for all of the investigated pockets. For the first time in the literature, PFAS are postulated to bind to the DNA binding domain in a nonspecific manner. In addition, for the PPARγ ligand binding domain, l-carnitine shows promise in replacing smaller PFAS from the pocket.
Collapse
Affiliation(s)
- Nuno M S Almeida
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48864, United States
| | - Semiha Kevser Bali
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48864, United States
| | - Deepak James
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48864, United States
| | - Cong Wang
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48864, United States
| | - Angela K Wilson
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48864, United States
| |
Collapse
|
9
|
Omrani S, Ghasemi M, Singh M, Mahmoodpour S, Zhou T, Babaei M, Niasar V. Interfacial Tension-Temperature-Pressure-Salinity Relationship for the Hydrogen-Brine System under Reservoir Conditions: Integration of Molecular Dynamics and Machine Learning. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2023; 39:12680-12691. [PMID: 37650690 PMCID: PMC10501201 DOI: 10.1021/acs.langmuir.3c01424] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/26/2023] [Revised: 08/14/2023] [Indexed: 09/01/2023]
Abstract
Hydrogen (H2) underground storage has attracted considerable attention as a potentially efficient strategy for the large-scale storage of H2. Nevertheless, successful execution and long-term storage and withdrawal of H2 necessitate a thorough understanding of the physical and chemical properties of H2 in contact with the resident fluids. As capillary forces control H2 migration and trapping in a subsurface environment, quantifying the interfacial tension (IFT) between H2 and the resident fluids in the subsurface is important. In this study, molecular dynamics (MD) simulation was employed to develop a data set for the IFT of H2-brine systems under a wide range of thermodynamic conditions (298-373 K temperatures and 1-30 MPa pressures) and NaCl salinities (0-5.02 mol·kg-1). For the first time to our knowledge, a comprehensive assessment was carried out to introduce the most accurate force field combination for H2-brine systems in predicting interfacial properties with an absolute relative deviation (ARD) of less than 3% compared with the experimental data. In addition, the effect of the cation type was investigated for brines containing NaCl, KCl, CaCl2, and MgCl2. Our results show that H2-brine IFT decreases with increasing temperature under any pressure condition, while higher NaCl salinity increases the IFT. A slight decrease in IFT occurs when the pressure increases. Under the impact of cation type, Ca2+ can increase IFT values more than others, i.e., up to 12% with respect to KCl. In the last step, the predicted IFT data set was used to provide a reliable correlation using machine learning (ML). Three white-box ML approaches of the group method of data handling (GMDH), gene expression programming (GEP), and genetic programming (GP) were applied. GP demonstrates the most accurate correlation with a coefficient of determination (R2) and absolute average relative deviation (AARD) of 0.9783 and 0.9767%, respectively.
Collapse
Affiliation(s)
- Sina Omrani
- Department
of Chemical Engineering, The University
of Manchester, Manchester M13 9PL, United
Kingdom
| | - Mehdi Ghasemi
- Department
of Chemical Engineering, The University
of Manchester, Manchester M13 9PL, United
Kingdom
| | - Mrityunjay Singh
- Institute
of Applied Geosciences, Geothermal Science and Technology, Technische Universität Darmstadt, 64289 Darmstadt, Germany
| | - Saeed Mahmoodpour
- Group
of Geothermal Technologies, Technische Universität
Munchen, 80333 Munich, Germany
| | - Tianhang Zhou
- College
of Carbon Neutrality Future Technology, China University of Petroleum (Beijing), 102249 Beijing, China
| | - Masoud Babaei
- Department
of Chemical Engineering, The University
of Manchester, Manchester M13 9PL, United
Kingdom
| | - Vahid Niasar
- Department
of Chemical Engineering, The University
of Manchester, Manchester M13 9PL, United
Kingdom
| |
Collapse
|
10
|
Korede V, Nagalingam N, Penha FM, van der Linden N, Padding JT, Hartkamp R, Eral HB. A Review of Laser-Induced Crystallization from Solution. CRYSTAL GROWTH & DESIGN 2023; 23:3873-3916. [PMID: 37159656 PMCID: PMC10161235 DOI: 10.1021/acs.cgd.2c01526] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 01/04/2023] [Indexed: 05/11/2023]
Abstract
Crystallization abounds in nature and industrial practice. A plethora of indispensable products ranging from agrochemicals and pharmaceuticals to battery materials are produced in crystalline form in industrial practice. Yet, our control over the crystallization process across scales, from molecular to macroscopic, is far from complete. This bottleneck not only hinders our ability to engineer the properties of crystalline products essential for maintaining our quality of life but also hampers progress toward a sustainable circular economy in resource recovery. In recent years, approaches leveraging light fields have emerged as promising alternatives to manipulate crystallization. In this review article, we classify laser-induced crystallization approaches where light-material interactions are utilized to influence crystallization phenomena according to proposed underlying mechanisms and experimental setups. We discuss nonphotochemical laser-induced nucleation, high-intensity laser-induced nucleation, laser trapping-induced crystallization, and indirect methods in detail. Throughout the review, we highlight connections among these separately evolving subfields to encourage the interdisciplinary exchange of ideas.
Collapse
Affiliation(s)
- Vikram Korede
- Process
& Energy Department, Delft University
of Technology, Leeghwaterstraat 39, 2628 CB Delft, The Netherlands
| | - Nagaraj Nagalingam
- Process
& Energy Department, Delft University
of Technology, Leeghwaterstraat 39, 2628 CB Delft, The Netherlands
| | - Frederico Marques Penha
- Department
of Chemical Engineering, KTH Royal Institute
of Technology, Teknikringen
42, 114-28 Stockholm, Sweden
| | - Noah van der Linden
- Process
& Energy Department, Delft University
of Technology, Leeghwaterstraat 39, 2628 CB Delft, The Netherlands
| | - Johan T. Padding
- Process
& Energy Department, Delft University
of Technology, Leeghwaterstraat 39, 2628 CB Delft, The Netherlands
| | - Remco Hartkamp
- Process
& Energy Department, Delft University
of Technology, Leeghwaterstraat 39, 2628 CB Delft, The Netherlands
| | - Huseyin Burak Eral
- Process
& Energy Department, Delft University
of Technology, Leeghwaterstraat 39, 2628 CB Delft, The Netherlands
| |
Collapse
|
11
|
Di Pasquale N, Finney AR, Elliott JD, Carbone P, Salvalaglio M. Constant chemical potential-quantum mechanical-molecular dynamics simulations of the graphene-electrolyte double layer. J Chem Phys 2023; 158:134714. [PMID: 37031135 DOI: 10.1063/5.0138267] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/10/2023] Open
Abstract
We present the coupling of two frameworks-the pseudo-open boundary simulation method known as constant potential molecular dynamics simulations (CμMD), combined with quantum mechanics/molecular dynamics (QMMD) calculations-to describe the properties of graphene electrodes in contact with electrolytes. The resulting CμQMMD model was then applied to three ionic solutions (LiCl, NaCl, and KCl in water) at bulk solution concentrations ranging from 0.5 M to 6 M in contact with a charged graphene electrode. The new approach we are describing here provides a simulation protocol to control the concentration of electrolyte solutions while including the effects of a fully polarizable electrode surface. Thanks to this coupling, we are able to accurately model both the electrode and solution side of the double layer and provide a thorough analysis of the properties of electrolytes at charged interfaces, such as the screening ability of the electrolyte and the electrostatic potential profile. We also report the calculation of the integral electrochemical double layer capacitance in the whole range of concentrations analyzed for each ionic species, while the quantum mechanical simulations provide access to the differential and integral quantum capacitance. We highlight how subtle features, such as the adsorption of potassium graphene or the tendency of the ions to form clusters contribute to the ability of graphene to store charge, and suggest implications for desalination.
Collapse
Affiliation(s)
- Nicodemo Di Pasquale
- Department of Chemical Engineering, Brunel University London, Uxbridge UB8 3PH, United Kingdom
| | - Aaron R Finney
- Department of Chemical Engineering, University College London, London WC1E 7JE, United Kingdom
| | - Joshua D Elliott
- Department of Chemical Engineering, University of Manchester, Manchester M13 9PL, United Kingdom
| | - Paola Carbone
- Department of Chemical Engineering, University of Manchester, Manchester M13 9PL, United Kingdom
| | - Matteo Salvalaglio
- Department of Chemical Engineering, University College London, London WC1E 7JE, United Kingdom
| |
Collapse
|
12
|
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: 11] [Impact Index Per Article: 11.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
|
13
|
Upadhyaya A, Panthi B, Verma S, Kumar S, Rajouria SK, Srivastava HK, Chandra P. Analogue and structure based approaches for modelling HIV-1 integrase inhibitors. J Biomol Struct Dyn 2023; 41:11946-11956. [PMID: 36734646 DOI: 10.1080/07391102.2023.2171129] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/27/2021] [Accepted: 12/25/2022] [Indexed: 02/04/2023]
Abstract
A set of 220 inhibitors belonging to different structure classes and having HIV-1 integrase activity were collected along with their experimental pIC50 values. Geometries of all the inhibitors were fully optimized using B3LYP/6-31 + G(d) level of theory. These ligands were docked against 4 different HIV-1 integrase receptors (PDB IDs: 4LH5, 5KRS, 3ZSQ and 3ZSV). 30 docked poses were generated for all 220 inhibitors and ligand interaction of the first docked pose and the docked pose with the highest score were analysed. Residue GLU170 of 4LH5 receptor shows the highest number of interactions followed by ALA169, GLN168, HIS171 and ASP167 residues. Hydrogen bonding and stacking are mainly responsible for the interactions of these inhibitors with the receptor. We performed Molecular Dynamics (MD) simulation to observe the root-mean-square deviation (RMSD), for measure the average change of displacement between the atoms for a particular frame with respect to a reference and The Root Mean Square Fluctuation (RMSF) for characterization of local changes along the protein chain of the docked complexes. Analogue based models were generated to predict the pIC50 values for integrase inhibitors using various types of descriptors such as constitutional, geometrical, topological, quantum chemical and docking based descriptors. The best models were selected on the basis of statistical parameters and were validated by training and test set division. A few new inhibitors were designed on the basis of structure activity relationship and their pIC50 values were predicted using the generated models. All the designed new inhibitors a very high potential and may be used as potent inhibitors of HIV integrase. These models may be useful for further design and development of new and potent HIV integrase inhibitors.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Anurag Upadhyaya
- Department of Physics, Institute of Science, Banaras Hindu University, Varanasi, Uttar Pradesh, India
| | - Bhavana Panthi
- Department of Chemistry, Indian Institute of Technology Kanpur, Kalyanpur Kanpur, Uttar Pradesh, India
| | - Shubham Verma
- Department of Medicinal Chemistry, National Institute of Pharmaceutical Education and Research, Changsari, Guwahati, Assam, India
- Department of Biosciences and Bioengineering, Indian Institute of Technology Guwahati, Guwahati, Assam, India
| | - Suresh Kumar
- Department of Physics, Institute of Science, Banaras Hindu University, Varanasi, Uttar Pradesh, India
- Department of Physics, Dyal Singh College, University of Delhi, Delhi, India
| | - Satish Kumar Rajouria
- Department of Physics, Zakir Husain Delhi College, University of Delhi, Delhi, India
| | - Hemant Kumar Srivastava
- Department of Medicinal Chemistry, National Institute of Pharmaceutical Education and Research, Changsari, Guwahati, Assam, India
| | - Pranjal Chandra
- Department of Biosciences and Bioengineering, Indian Institute of Technology Guwahati, Guwahati, Assam, India
| |
Collapse
|
14
|
Borkowski AK, Thompson WH. Shining (Infrared) Light on the Hofmeister Series: Driving Forces for Changes in the Water Vibrational Spectra in Alkali-Halide Salt Solutions. J Phys Chem B 2022; 126:6700-6712. [PMID: 36004804 DOI: 10.1021/acs.jpcb.2c03957] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The Hofmeister series is frequently used to rank ions based on their behavior from chaotropes ("structure breakers"), which weaken the surrounding hydrogen-bond network, to kosmotropes ("structure makers"), which enhance it. Here, we use fluctuation theory to investigate the energetic and entropic driving forces underlying the Hofmeister series for aqueous alkali-halide solutions. Specifically, we exploit the OH stretch infrared (IR) spectrum in isotopically dilute HOD/D2O solutions as a probe of the effect of the salt on the water properties for different concentrations and choice of halide anion. Fluctuation theory is used to calculate the temperature derivative of these IR spectra, including decomposition of the derivative into different energetic contributions. These contributions are used to determine the thermodynamic driving forces in terms of effective internal energy and entropic contributions. This analysis implicates entropic contributions as the key factor in the Hofmeister series behavior of the OH stretch IR spectra, while the effective internal energy is nearly ion-independent.
Collapse
Affiliation(s)
- Ashley K Borkowski
- Department of Chemistry, University of Kansas, Lawrence, Kansas 66045, United States
| | - Ward H Thompson
- Department of Chemistry, University of Kansas, Lawrence, Kansas 66045, United States
| |
Collapse
|
15
|
Biriukov D, Wang HW, Rampal N, Tempra C, Kula P, Neuefeind JC, Stack AG, Předota M. The "good," the "bad," and the "hidden" in neutron scattering and molecular dynamics of ionic aqueous solutions. J Chem Phys 2022; 156:194505. [PMID: 35597655 DOI: 10.1063/5.0093643] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022] Open
Abstract
We characterize a concentrated 7.3 m CaCl2 solution, combining neutron diffraction with chloride isotopic substitution (Cl-NDIS) in null water and molecular dynamics (MD) simulations. We elucidate the solution structure, thermodynamic properties, and extent of ion pairing previously suggested as concentration-dependent and often not observed at lower concentrations. Our Cl-NDIS measurements designate the solvent-shared ion pairing as dominant and the contact ion pairing (CIP) as insignificant even under conditions close to the solubility limit. The MD models parameterized against neutron diffraction with calcium isotopic substitution (Ca-NDIS) overestimate CIP despite successfully reproducing most of the Cl-NDIS signal. This drawback originates from the fact that Ca2+-Cl- interactions were primarily "hidden" in the Ca-NDIS signal due to overlapping with Ca2+-Ow and Ca2+-Hw contributions to the total scattering. Contrary, MD models with moderate CIP and possessing generally good performance at high concentrations fail to reproduce the NDIS measurements accurately. Therefore, the electronic polarization, introduced in most of the recent MD models via scaling ionic charges, resolves some but not all parameterization drawbacks. We conclude that despite improving the quality of MD models "on average," the question "which model is the best" has not been answered but replaced by the question "which model is better for a given research." An overall "good" model can still be inappropriate or, in some instances, "bad" and, unfortunately, produce erroneous results. The accurate interpretation of several NDIS datasets, complemented by MD simulations, can prevent such mistakes and help identify the strengths, weaknesses, and convenient applications for corresponding computational models.
Collapse
Affiliation(s)
- Denys Biriukov
- Institute of Organic Chemistry and Biochemistry, Czech Academy of Sciences, Flemingovo nám. 2, 16610 Prague 6, Czech Republic
| | - Hsiu-Wen Wang
- Chemical Sciences Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, USA
| | - Nikhil Rampal
- Chemical Sciences Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, USA
| | - Carmelo Tempra
- Institute of Organic Chemistry and Biochemistry, Czech Academy of Sciences, Flemingovo nám. 2, 16610 Prague 6, Czech Republic
| | - Patrik Kula
- Institute of Organic Chemistry and Biochemistry, Czech Academy of Sciences, Flemingovo nám. 2, 16610 Prague 6, Czech Republic
| | - Joerg C Neuefeind
- Neutron Scattering Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, USA
| | - Andrew G Stack
- Chemical Sciences Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, USA
| | - Milan Předota
- Department of Physics, Faculty of Science, University of South Bohemia, Branišovská 1760, 370 05 České Budějovice, Czech Republic
| |
Collapse
|
16
|
Zofchak ES, Zhang Z, Marioni N, Duncan TJ, Sachar HS, Chamseddine A, Freeman BD, Ganesan V. Cation–Ligand Interactions Dictate Salt Partitioning and Diffusivity in Ligand-Functionalized Polymer Membranes. Macromolecules 2022. [DOI: 10.1021/acs.macromol.2c00035] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Affiliation(s)
- Everett S. Zofchak
- McKetta Department of Chemical Engineering, University of Texas at Austin, Austin, Texas 78712, United States
| | - Zidan Zhang
- McKetta Department of Chemical Engineering, University of Texas at Austin, Austin, Texas 78712, United States
| | - Nico Marioni
- McKetta Department of Chemical Engineering, University of Texas at Austin, Austin, Texas 78712, United States
| | - Tyler J. Duncan
- McKetta Department of Chemical Engineering, University of Texas at Austin, Austin, Texas 78712, United States
| | - Harnoor S. Sachar
- McKetta Department of Chemical Engineering, University of Texas at Austin, Austin, Texas 78712, United States
| | - Alyssa Chamseddine
- McKetta Department of Chemical Engineering, University of Texas at Austin, Austin, Texas 78712, United States
| | - Benny D. Freeman
- McKetta Department of Chemical Engineering, University of Texas at Austin, Austin, Texas 78712, United States
| | - Venkat Ganesan
- McKetta Department of Chemical Engineering, University of Texas at Austin, Austin, Texas 78712, United States
| |
Collapse
|
17
|
Grotz KK, Schwierz N. Magnesium Force Fields for OPC Water with Accurate Solvation, Ion-Binding, and Water-Exchange Properties: Successful Transfer from SPC/E. J Chem Phys 2022; 156:114501. [DOI: 10.1063/5.0087292] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Magnesium plays a vital role in a large variety of biological processes. To model such processes by molecular dynamics simulations, researchers rely on accurate force field parameters for Mg2+ and water. OPC is one of the most promising water models yielding an improved description of biomolecules in water. The aim of this work is to provide force field parameters for Mg2+ that lead to accurate simulation results in combination with OPC water. Using twelve different Mg2+ parameter sets, that were previously optimized with different water models, we systematically assess the transferability to OPC based on a large variety of experimental properties. The results show that the Mg2+ parameters for SPC/E are transferable to OPC and closely reproduce the experimental solvation free energy, radius of the first hydration shell, coordination number, activity derivative, and binding affinity toward the phosphate oxygens on RNA. Two optimal parameter sets are presented: MicroMg yields water exchange in OPC on the microsecond timescale in agreement with experiments. NanoMg yields accelerated exchange on the nanosecond timescale and facilitates the direct observation of ion binding events for enhanced sampling purposes.
Collapse
Affiliation(s)
- Kara K. Grotz
- Theoretical Biophysics, Max Planck Institute of Biophysics, Germany
| | - Nadine Schwierz
- Theoretical Biophysics, Max Planck Institute of Biophysics, Germany
| |
Collapse
|
18
|
Döpke MF, Westerbaan van der Meij F, Coasne B, Hartkamp R. Surface Protolysis and Its Kinetics Impact the Electrical Double Layer. PHYSICAL REVIEW LETTERS 2022; 128:056001. [PMID: 35179914 DOI: 10.1103/physrevlett.128.056001] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/14/2021] [Revised: 12/08/2021] [Accepted: 01/05/2022] [Indexed: 06/14/2023]
Abstract
Surface conductivity in the electrical double layer (EDL) is known to be affected by proton hopping and diffusion at solid-liquid interfaces. Yet, the role of surface protolysis and its kinetics on the thermodynamic and transport properties of the EDL are usually ignored as physical models consider static surfaces. Here, using a novel molecular dynamics method mimicking surface protolysis, we unveil the impact of such chemical events on the system's response. Protolysis is found to strongly affect the EDL and electrokinetic aspects with major changes in ζ potential and electro-osmotic flow.
Collapse
Affiliation(s)
- Max F Döpke
- Process & Energy Department, Delft University of Technology, Leeghwaterstraat 39, 2628 CB Delft, Netherlands
| | | | - Benoit Coasne
- Université Grenoble Alpes, CNRS, LIPhy, 38000 Grenoble, France
| | - Remco Hartkamp
- Process & Energy Department, Delft University of Technology, Leeghwaterstraat 39, 2628 CB Delft, Netherlands
| |
Collapse
|
19
|
Application of Imidazolium-based polyionic liquids to separate the 1,3,5-Trioxane-Water/Ethanol-Water system based on experimental verification and molecular mechanism analysis. J Mol Liq 2022. [DOI: 10.1016/j.molliq.2021.118079] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/24/2023]
|
20
|
Blazquez S, Conde MM, Abascal JLF, Vega C. The Madrid-2019 force field for electrolytes in water using TIP4P/2005 and scaled charges: Extension to the ions F−, Br−, I−, Rb+, and Cs+. J Chem Phys 2022; 156:044505. [DOI: 10.1063/5.0077716] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022] Open
Affiliation(s)
- S. Blazquez
- Departamento Química Física I, Facultad de 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
| | - J. L. F. Abascal
- Departamento Química Física I, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - C. Vega
- Departamento Química Física I, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| |
Collapse
|
21
|
Grotz KK, Schwierz N. Optimized Magnesium Force Field Parameters for Biomolecular Simulations with Accurate Solvation, Ion-Binding, and Water-Exchange Properties in SPC/E, TIP3P-fb, TIP4P/2005, TIP4P-Ew, and TIP4P-D. J Chem Theory Comput 2022; 18:526-537. [PMID: 34881568 PMCID: PMC8757469 DOI: 10.1021/acs.jctc.1c00791] [Citation(s) in RCA: 13] [Impact Index Per Article: 6.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/06/2021] [Indexed: 12/14/2022]
Abstract
Magnesium is essential in many vital processes. To correctly describe Mg2+ in physiological processes by molecular dynamics simulations, accurate force fields are fundamental. Despite the importance, force fields based on the commonly used 12-6 Lennard-Jones potential showed significant shortcomings. Recently progress was made by an optimization procedure that implicitly accounts for polarizability. The resulting microMg and nanoMg force fields (J. Chem. Theory Comput. 2021, 17, 2530-2540) accurately reproduce a broad range of experimental solution properties and the binding affinity to nucleic acids in TIP3P water. Since countless simulation studies rely on available water models and ion force fields, we here extend the optimization and provide Mg2+ parameters in combination with the SPC/E, TIP3P-fb, TIP4P/2005, TIP4P-Ew, and TIP4P-D water models. For each water model, the Mg2+ force fields reproduce the solvation free energy, the distance to oxygens in the first hydration shell, the hydration number, the activity coefficient derivative in MgCl2 solutions, and the binding affinity and distance to the phosphate oxygens on nucleic acids. We present two parameter sets: MicroMg yields water exchange on the microsecond time scale and matches the experimental exchange rate. Depending on the water model, nanoMg yields accelerated water exchange in the range of 106 to 108 exchanges per second. The nanoMg parameters can be used to enhance the sampling of binding events, to obtain converged distributions of Mg2+, or to predict ion binding sites in biomolecular simulations. The parameter files are freely available at https://github.com/bio-phys/optimizedMgFFs.
Collapse
Affiliation(s)
- Kara K. Grotz
- Department of Theoretical
Biophysics, Max-Planck-Institute of Biophysics, Frankfurt am Main 60438, Germany
| | - Nadine Schwierz
- Department of Theoretical
Biophysics, Max-Planck-Institute of Biophysics, Frankfurt am Main 60438, Germany
| |
Collapse
|
22
|
Zhu R, Li G, Lei Z, Gui C. Mechanistic insight into absorption performance assessment for SO2 by mixed ionic liquids. J Mol Liq 2021. [DOI: 10.1016/j.molliq.2021.117927] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/20/2022]
|
23
|
Chen Z, Dai Y, Chi S, Su Z, Xing J, Wang Y, Lu Y. Analysis and intensification of energy saving process for separation of azeotrope by ionic liquid extractive distillation based on molecular dynamics simulation. Sep Purif Technol 2021. [DOI: 10.1016/j.seppur.2021.119254] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2023]
|
24
|
Warnock SJ, Sujanani R, Zofchak ES, Zhao S, Dilenschneider TJ, Hanson KG, Mukherjee S, Ganesan V, Freeman BD, Abu-Omar MM, Bates CM. Engineering Li/Na selectivity in 12-Crown-4-functionalized polymer membranes. Proc Natl Acad Sci U S A 2021; 118:e2022197118. [PMID: 34493651 PMCID: PMC8449368 DOI: 10.1073/pnas.2022197118] [Citation(s) in RCA: 39] [Impact Index Per Article: 13.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
Lithium is widely used in contemporary energy applications, but its isolation from natural reserves is plagued by time-consuming and costly processes. While polymer membranes could, in principle, circumvent these challenges by efficiently extracting lithium from aqueous solutions, they usually exhibit poor ion-specific selectivity. Toward this end, we have incorporated host-guest interactions into a tunable polynorbornene network by copolymerizing 1) 12-crown-4 ligands to impart ion selectivity, 2) poly(ethylene oxide) side chains to control water content, and 3) a crosslinker to form robust solids at room temperature. Single salt transport measurements indicate these materials exhibit unprecedented reverse permeability selectivity (∼2.3) for LiCl over NaCl-the highest documented to date for a dense, water-swollen polymer. As demonstrated by molecular dynamics simulations, this behavior originates from the ability of 12-crown-4 to bind Na+ ions more strongly than Li+ in an aqueous environment, which reduces Na+ mobility (relative to Li+) and offsets the increase in Na+ solubility due to binding with crown ethers. Under mixed salt conditions, 12-crown-4 functionalized membranes showed identical solubility selectivity relative to single salt conditions; however, the permeability and diffusivity selectivity of LiCl over NaCl decreased, presumably due to flux coupling. These results reveal insights for designing advanced membranes with solute-specific selectivity by utilizing host-guest interactions.
Collapse
Affiliation(s)
- Samuel J Warnock
- Materials Department, University of California, Santa Barbara, CA 93106
| | - Rahul Sujanani
- McKetta Department of Chemical Engineering, The University of Texas at Austin, Austin, TX 78712
| | - Everett S Zofchak
- McKetta Department of Chemical Engineering, The University of Texas at Austin, Austin, TX 78712
| | - Shou Zhao
- Department of Chemistry & Biochemistry, University of California, Santa Barbara, CA 93106
| | | | - Kalin G Hanson
- Department of Chemistry & Biochemistry, University of California, Santa Barbara, CA 93106
| | - Sanjoy Mukherjee
- Materials Research Laboratory, University of California, Santa Barbara, CA 93106
| | - Venkat Ganesan
- McKetta Department of Chemical Engineering, The University of Texas at Austin, Austin, TX 78712;
| | - Benny D Freeman
- McKetta Department of Chemical Engineering, The University of Texas at Austin, Austin, TX 78712;
| | - Mahdi M Abu-Omar
- Department of Chemistry & Biochemistry, University of California, Santa Barbara, CA 93106;
- Department of Chemical Engineering, University of California, Santa Barbara, CA 93106
| | - Christopher M Bates
- Materials Department, University of California, Santa Barbara, CA 93106;
- Department of Chemistry & Biochemistry, University of California, Santa Barbara, CA 93106
- Materials Research Laboratory, University of California, Santa Barbara, CA 93106
- Department of Chemical Engineering, University of California, Santa Barbara, CA 93106
| |
Collapse
|
25
|
Zhang Y, Jiang Y, Peng J, Zhang H. Rational Design of Nonbonded Point Charge Models for Divalent Metal Cations with Lennard-Jones 12-6 Potential. J Chem Inf Model 2021; 61:4031-4044. [PMID: 34313132 DOI: 10.1021/acs.jcim.1c00580] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Abstract
Exploring a metal-involved biochemical process at a molecular level often requires a reliable description of metal properties in aqueous solution by classical nonbonded models. An additional C4 term for considering ion-induced dipole interactions was previously proposed to supplement the widely used Lennard-Jones 12-6 potential (known as the 12-6-4 LJ-type model) with good accuracy. Here, we demonstrate an alternative to modeling divalent metal cations (M2+) with the traditional 12-6 LJ potential by developing nonbonded point charge models for use with 11 water models: TIP3P, SPC/E, SPC/Eb, TIP4P-Ew, TIP4P-D, and TIP4P/2005 and the more recent OPC3, TIP3P-FB, OPC, TIP4P-FB, and a99SB-disp. Our designed models simultaneously reproduce the experimental hydration free energy, ion-oxygen distance, and coordination number in the first hydration shell accurately for most of the metal cations, an accuracy equivalent to that of the complex 12-6-4 LJ-type and double exponential potential models. A systematic comparison with the existing M2+ models is presented as well in terms of effective ion radii, diffusion constants, water exchange rates, and ion-water interactions. Molecular dynamics simulations of metal substitution in Escherichia coli glyoxalase I variants show the great potential of our new models for metalloproteins.
Collapse
Affiliation(s)
- Yongguang Zhang
- Department of Biological Science and Engineering, School of Chemistry and Biological Engineering, University of Science and Technology Beijing, Beijing 100083, China
| | - Yang Jiang
- Department of Chemistry, Pennsylvania State University, University Park 16802, Pennsylvania, United States
| | - Jiarong Peng
- Department of Biological Science and Engineering, School of Chemistry and Biological Engineering, University of Science and Technology Beijing, Beijing 100083, China
| | - Haiyang Zhang
- Department of Biological Science and Engineering, School of Chemistry and Biological Engineering, University of Science and Technology Beijing, Beijing 100083, China
| |
Collapse
|
26
|
Abstract
Many important processes affecting the earth's climate are determined by the physical properties of seawater. In addition, desalination of seawater is a significant source of drinking water for the human population living in coastal areas. Since the physical properties of seawater governing these processes depend on the molecular interactions among its components, a deeper knowledge of seawater at the molecular level would contribute to a better understanding of these phenomena. However, in strong contrast with the situation in other areas such as biomolecules or materials science, molecular simulation studies reporting the physical properties of seawater are currently lacking. This is probably due to the usual perception of the seawater composition being too complex to approach. This point of view ignores the fact that physical properties of seawater are dependent on a single parameter representing the composition, namely the salinity. This is because the relative proportions of any two major constituents of seasalt are always the same. Another obstacle to performing molecular simulations of seawater could have been the unavailability of a satisfactory force field representing the interactions between water molecules and dissolved substances. However, this drawback has recently been overcome with the proposal of the Madrid-2019 force field. In this work we show for the first time that molecular simulation of seawater is feasible. We have performed molecular dynamics simulations of a system, the composition of which is close to the average composition of standard seawater and with the molecular interactions given by the Madrid-2019 force field. In this way we are able to provide quantitative or semiquantitative predictions for a number of relevant physical properties of seawater for temperatures and salinities from the oceanographic range to those relevant to desalination processes. The computed magnitudes include static (density), dynamical (viscosity and diffusion coefficients), structural (ionic hydration, ion-ion distribution functions), and interfacial (surface tension) properties.
Collapse
Affiliation(s)
- Ivan M Zeron
- Departamento de Química Física, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - Miguel A Gonzalez
- Departamento de Química Física, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - Edoardo Errani
- Departamento de Química Física, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - Carlos Vega
- Departamento de Química Física, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - Jose L F Abascal
- Departamento de Química Física, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| |
Collapse
|
27
|
Döpke MF, Hartkamp R. The importance of specifically adsorbed ions for electrokinetic phenomena: Bridging the gap between experiments and MD simulations. J Chem Phys 2021; 154:094701. [DOI: 10.1063/5.0038161] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023] Open
Affiliation(s)
- Max F. Döpke
- Process & Energy Department, Delft University of Technology, Leeghwaterstraat 39, 2628 CB Delft, The Netherlands
| | - Remco Hartkamp
- Process & Energy Department, Delft University of Technology, Leeghwaterstraat 39, 2628 CB Delft, The Netherlands
| |
Collapse
|
28
|
Patel LA, Yoon TJ, Currier RP, Maerzke KA. NaCl aggregation in water at elevated temperatures and pressures: Comparison of classical force fields. J Chem Phys 2021; 154:064503. [PMID: 33588550 DOI: 10.1063/5.0030962] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022] Open
Abstract
The properties of water vary dramatically with temperature and density. This can be exploited to control its effectiveness as a solvent. Thus, supercritical water is of keen interest as solvent in many extraction processes. The low solubility of salts in lower density supercritical water has even been suggested as a means of desalination. The high temperatures and pressures required to reach supercritical conditions can present experimental challenges during collection of required physical property and phase equilibria data, especially in salt-containing systems. Molecular simulations have the potential to be a valuable tool for examining the behavior of solvated ions at these high temperatures and pressures. However, the accuracy of classical force fields under these conditions is unclear. We have, therefore, undertaken a parametric study of NaCl in water, comparing several salt and water models at 200 bar-600 bar and 450 K-750 K for a range of salt concentrations. We report a comparison of structural properties including ion aggregation, hydrogen bonding, density, and static dielectric constants. All of the force fields qualitatively reproduce the trends in the liquid phase density. An increase in ion aggregation with decreasing density holds true for all of the force fields. The propensity to aggregate is primarily determined by the salt force field rather than the water force field. This coincides with a decrease in the water static dielectric constant and reduced charge screening. While a decrease in the static dielectric constant with increasing NaCl concentration is consistent across all model combinations, the salt force fields that exhibit more ionic aggregation yield a slightly smaller dielectric decrement.
Collapse
Affiliation(s)
- Lara A Patel
- Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA
| | - Tae Jun Yoon
- Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA
| | - Robert P Currier
- Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA
| | - Katie A Maerzke
- Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA
| |
Collapse
|
29
|
Erdős M, Frangou M, Vlugt TJH, Moultos OA. Diffusivity of α-, β-, γ-cyclodextrin and the inclusion complex of β-cyclodextrin: Ibuprofen in aqueous solutions; A molecular dynamics simulation study. FLUID PHASE EQUILIBRIA 2021; 528:112842. [PMID: 33024350 PMCID: PMC7529625 DOI: 10.1016/j.fluid.2020.112842] [Citation(s) in RCA: 22] [Impact Index Per Article: 7.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/25/2020] [Revised: 09/03/2020] [Accepted: 09/18/2020] [Indexed: 05/28/2023]
Abstract
Cyclodextrins (CDs) are widely used in drug delivery, catalysis, food and separation processes. In this work, a comprehensive simulation study on the diffusion of the native α-, β- and γ-CDs in aqueous solutions is carried out using Molecular Dynamics simulations. The effect of the system size on the computed self-diffusivity is investigated and it is found that the required correction can be as much as 75% of the final value. The effect of the water force field is examined and it is shown that the q4md-CD/TIP4P/2005 force field combination predicts the experimentally measured self-diffusion coefficients of CDs very accurately. The self-diffusion coefficients of the three native CDs were also computed in aqueous-NaCl solutions using the Joung and Cheatham (JC) and the Madrid-2019 force fields. It is found that Na+ ions have higher affinity towards the CDs when the JC force field is used and for this reason the predicted diffusivity of CDs is lower compared to simulations using the Madrid-2019 force field. As a model system for drug delivery and waste-water treatment applications, the diffusion of the β-CD:Ibuprofen inclusion complex in water is studied. In agreement with experiments for similar components, it is shown that the inclusion complex and the free β-CD have almost equal self-diffusion coefficients. Our analysis revealed that this is most likely caused by the almost full inclusion of the ibuprofen in the cavity of the β-CD. Our findings show that Molecular Dynamics simulation can be used to provide reasonable diffusivity predictions, and to obtain molecular-level understanding useful for industrial applications of CDs.
Collapse
Affiliation(s)
- Máté Erdős
- Engineering Thermodynamics, Process & Energy Department, Faculty of Mechanical, Maritime and Materials Engineering, Delft University of Technology, Leeghwaterstraat 39, Delft 2628CB, The Netherlands
| | - Michalis Frangou
- Engineering Thermodynamics, Process & Energy Department, Faculty of Mechanical, Maritime and Materials Engineering, Delft University of Technology, Leeghwaterstraat 39, Delft 2628CB, The Netherlands
| | - Thijs J H Vlugt
- Engineering Thermodynamics, Process & Energy Department, Faculty of Mechanical, Maritime and Materials Engineering, Delft University of Technology, Leeghwaterstraat 39, Delft 2628CB, The Netherlands
| | - Othonas A Moultos
- Engineering Thermodynamics, Process & Energy Department, Faculty of Mechanical, Maritime and Materials Engineering, Delft University of Technology, Leeghwaterstraat 39, Delft 2628CB, The Netherlands
| |
Collapse
|
30
|
Zhang W, Ma Z, Liu X, Liu Y, Hou Z, Qi J, Ma Y, Wang L, Wang Y. Molecular Mechanism and Absorption Performance Evaluation of CO 2 Capture from the PCC Process by Monoethanolamine-Based Deep Eutectic Solvents. Ind Eng Chem Res 2021. [DOI: 10.1021/acs.iecr.0c04970] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022]
Affiliation(s)
- Wanxiang Zhang
- College of Chemical Engineering, Qingdao University of Science and Technology, 53 Zhengzhou Road, Qingdao 266042, P. R. China
| | - Zhaoyuan Ma
- College of Chemical Engineering, Qingdao University of Science and Technology, 53 Zhengzhou Road, Qingdao 266042, P. R. China
| | - Xingyi Liu
- College of Chemical Engineering, Qingdao University of Science and Technology, 53 Zhengzhou Road, Qingdao 266042, P. R. China
| | - Yigang Liu
- College of Chemical Engineering, Qingdao University of Science and Technology, 53 Zhengzhou Road, Qingdao 266042, P. R. China
| | - Zhengkun Hou
- College of Chemical Engineering, Qingdao University of Science and Technology, 53 Zhengzhou Road, Qingdao 266042, P. R. China
| | - Jianguang Qi
- College of Chemical Engineering, Qingdao University of Science and Technology, 53 Zhengzhou Road, Qingdao 266042, P. R. China
| | - Yixin Ma
- College of Chemical and Environmental Engineering, Shandong University of Science and Technology, 579 Qianwangang Road, Qingdao 266590, P. R. China
| | - Lei Wang
- Key Laboratory of Eco-chemical Engineering, Taishan Scholar Advantage and Characteristic Discipline Team of Eco Chemical Process and Technology, College of Chemistry and Molecular Engineering, Qingdao University of Science and Technology, Qingdao 266042, P. R. China
| | - Yinglong Wang
- College of Chemical Engineering, Qingdao University of Science and Technology, 53 Zhengzhou Road, Qingdao 266042, P. R. China
| |
Collapse
|
31
|
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
|
32
|
Eken Y, Almeida NMS, Wang C, Wilson AK. SAMPL7: Host-guest binding prediction by molecular dynamics and quantum mechanics. J Comput Aided Mol Des 2020; 35:63-77. [PMID: 33150463 DOI: 10.1007/s10822-020-00357-3] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/03/2020] [Accepted: 10/28/2020] [Indexed: 01/11/2023]
Abstract
Statistical Assessment of Modeling of Proteins and Ligands (SAMPL) challenges provide routes to compare chemical quantities determined using computational chemistry approaches to experimental measurements that are shared after the competition. For this effort, several computational methods have been used to calculate the binding energies of Octa Acid (OA) and exo-Octa Acid (exoOA) host-guest systems for SAMPL7. The initial poses for molecular dynamics (MD) were generated by molecular docking. Binding free energy calculations were performed using molecular mechanics combined with Poisson-Boltzmann or generalized Born surface area solvation (MMPBSA/MMGBSA) approaches. The factors that affect the utility of the MMPBSA/MMGBSA approaches including solvation, partial charge, and solute entropy models were also analyzed. In addition to MD calculations, quantum mechanics (QM) calculations were performed using several different density functional theory (DFT) approaches. From SAMPL6 results, B3PW91-D3 was found to overestimate binding energies though it was effective for geometry optimizations, so it was considered for the DFT geometry optimizations in the current study, with single-point energy calculations carried out with B2PLYP-D3 with double-, triple-, and quadruple-ζ level basis sets. Accounting for dispersion effects, and solvation models was deemed essential for the predictions. MMGBSA and MMPBSA correlated better to experiment when used in conjunction with an empirical/linear correction.
Collapse
Affiliation(s)
- Yiğitcan Eken
- Department of Chemistry, Michigan State University, East Lansing, MI, 48864, USA
| | - Nuno M S Almeida
- Department of Chemistry, Michigan State University, East Lansing, MI, 48864, USA
| | - Cong Wang
- Department of Chemistry, Michigan State University, East Lansing, MI, 48864, USA
| | - Angela K Wilson
- Department of Chemistry, Michigan State University, East Lansing, MI, 48864, USA.
- Department of Chemistry, University of North Texas, Denton, TX, 76201, USA.
| |
Collapse
|
33
|
Kumar N, Naik PK, Banerjee T. Molecular modeling insights in the extraction of benzene from hydrocarbon stream using deep eutectic solvent. J Mol Liq 2020. [DOI: 10.1016/j.molliq.2020.113909] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
|
34
|
Panagiotopoulos AZ. Simulations of activities, solubilities, transport properties, and nucleation rates for aqueous electrolyte solutions. J Chem Phys 2020; 153:010903. [DOI: 10.1063/5.0012102] [Citation(s) in RCA: 17] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/23/2023] Open
|
35
|
Urano R, Shinoda W, Yoshii N, Okazaki S. Exact long-range Coulombic energy calculation for net charged systems neutralized by uniformly distributed background charge using fast multipole method and its application to efficient free energy calculation. J Chem Phys 2020; 152:244115. [DOI: 10.1063/5.0007957] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Ryo Urano
- Department of Materials Chemistry, Nagoya University, Nagoya 464-8603, Japan
| | - Wataru Shinoda
- Department of Materials Chemistry, Nagoya University, Nagoya 464-8603, Japan
| | - Noriyuki Yoshii
- Department of Materials Chemistry, Nagoya University, Nagoya 464-8603, Japan
- Center for Computational Science, Graduate School of Engineering, Nagoya University, Nagoya 464-8603, Japan
| | - Susumu Okazaki
- Department of Materials Chemistry, Nagoya University, Nagoya 464-8603, Japan
- Center for Computational Science, Graduate School of Engineering, Nagoya University, Nagoya 464-8603, Japan
| |
Collapse
|
36
|
Abal JPK, Bordin JR, Barbosa MC. Salt parameterization can drastically affect the results from classical atomistic simulations of water desalination by MoS 2 nanopores. Phys Chem Chem Phys 2020; 22:11053-11061. [PMID: 32373906 DOI: 10.1039/d0cp00484g] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/30/2023]
Abstract
Water scarcity is a reality in our world, and scenarios predicted by leading scientists in this area indicate that it will worsen in the next decades. However, new technologies based on low-cost seawater desalination can prevent the worst scenarios, providing fresh water for humanity. With this goal, membranes based on nanoporous materials have been suggested in recent years. One of the materials suggested is MoS2, and classical Molecular Dynamics (MD) simulation is one of the most powerful tools to explore these nanomaterials. However, distinct force fields employed in MD simulations are parameterized based on distinct experimental quantities. In this paper, we compare two models of salt that were built based on distinct properties of water-salt mixtures. One model fits the hydration free energy and lattice properties, and the second fits the crystal density and the density and the dielectric constant of water and salt mixtures. To compare the models, MD simulations for salty water flow through nanopores of two sizes were used - one pore big enough to accommodate hydrated ions, and one smaller in which the ion has to dehydrate to enter - and two rigid water models from the TIP4P family - TIP4P/2005 and TIP4P/ε. Our results indicate that the water permeability and salt rejection by the membrane are more influenced by the salt model than by the water model, especially for the narrow pore. In fact, completely distinct mechanisms were observed, and they are related to the characteristics employed in the ion model parameterization. The results show that not only can the water model influence the outcomes, but the ion model plays a crucial role when the pore is small enough.
Collapse
Affiliation(s)
- João P K Abal
- Institute of Physics, Federal University of Rio Grande do Sul, 91501-970, Porto Alegre, Brazil.
| | | | | |
Collapse
|