1
|
Prosmiti R, González-Lezana T. Computational Modeling: Up-to-Date Approaches and Cutting-Edge Applications from Clusters, Nanostructures to Bulk Systems. Chemphyschem 2024:e202400207. [PMID: 38837591 DOI: 10.1002/cphc.202400207] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/26/2024] [Indexed: 06/07/2024]
Abstract
The contributions in this special theme collection, in honor to Prof. P. Villarreal, cover a broad variety of computational methodologies and experimental techniques, containing studies on gas phase, clusters and condensed phase systems.
Collapse
Affiliation(s)
- Rita Prosmiti
- Instituto de Fìsica Fundamental, IFF-CSIC, Serrano 123, 28006, Madrid, Spain
| | | |
Collapse
|
2
|
Yuan H, Zhang Y, Huang X, Zhang X, Li J, Huang Y, Li K, Weng H, Xu Y, Zhang Y. Exploration of the Existence Forms and Patterns of Dissolved Oxygen Molecules in Water. NANO-MICRO LETTERS 2024; 16:208. [PMID: 38833205 DOI: 10.1007/s40820-024-01427-z] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/27/2024] [Accepted: 04/25/2024] [Indexed: 06/06/2024]
Abstract
The structure of liquid water is primarily composed of three-dimensional networks of water clusters formed by hydrogen bonds, and dissolved oxygen is one of the most important indicators for assessing water quality. In this work, distilled water with different concentration of dissolved oxygen were prepared, and a clear negative correlation between the size of water clusters and dissolved oxygen concentration was observed. Besides, a phenomenon of rapid absorption and release of oxygen at the water interfaces was unveiled, suggesting that oxygen molecules predominantly exist at the interfaces of water clusters. Oxygen molecules can move rapidly through the interfaces among water clusters, allowing dissolved oxygen to quickly reach a saturation level at certain partial pressure of oxygen and temperature. Further exploration into the mechanism by molecular dynamics simulations of oxygen and water clusters found that oxygen molecules can only exist stably at the interfaces among water clusters. A semi-empirical formula relating the average number of water molecules in a cluster (n) to 17O NMR half-peak width (W) was summarized: n = 0.1 W + 0.85. These findings provide a foundation for exploring the structure and properties of water.
Collapse
Affiliation(s)
- Hewei Yuan
- Key Laboratory of Thin Film and Microfabrication (Ministry of Education), Department of Micro/Nano Electronics, School of Electronic Information and Electrical Engineering, Shanghai Jiao Tong University, Shanghai, 200240, People's Republic of China
| | - Yaozhong Zhang
- School of Materials Science and Engineering, Shanghai Jiao Tong University, Shanghai, 200240, People's Republic of China.
| | - Xiaolu Huang
- Key Laboratory of Thin Film and Microfabrication (Ministry of Education), Department of Micro/Nano Electronics, School of Electronic Information and Electrical Engineering, Shanghai Jiao Tong University, Shanghai, 200240, People's Republic of China
| | - Xiwu Zhang
- Jinduo Yuchen Water Environment Engineering Co., Ltd, Shanghai, 201702, People's Republic of China
| | - Jinjin Li
- Key Laboratory of Thin Film and Microfabrication (Ministry of Education), Department of Micro/Nano Electronics, School of Electronic Information and Electrical Engineering, Shanghai Jiao Tong University, Shanghai, 200240, People's Republic of China
| | - Yufeng Huang
- Jinduo Yuchen Water Environment Engineering Co., Ltd, Shanghai, 201702, People's Republic of China
| | - Kun Li
- Key Laboratory of Thin Film and Microfabrication (Ministry of Education), Department of Micro/Nano Electronics, School of Electronic Information and Electrical Engineering, Shanghai Jiao Tong University, Shanghai, 200240, People's Republic of China
| | - Haotian Weng
- Key Laboratory of Thin Film and Microfabrication (Ministry of Education), Department of Micro/Nano Electronics, School of Electronic Information and Electrical Engineering, Shanghai Jiao Tong University, Shanghai, 200240, People's Republic of China
| | - Yang Xu
- Key Laboratory of Thin Film and Microfabrication (Ministry of Education), Department of Micro/Nano Electronics, School of Electronic Information and Electrical Engineering, Shanghai Jiao Tong University, Shanghai, 200240, People's Republic of China
| | - Yafei Zhang
- Key Laboratory of Thin Film and Microfabrication (Ministry of Education), Department of Micro/Nano Electronics, School of Electronic Information and Electrical Engineering, Shanghai Jiao Tong University, Shanghai, 200240, People's Republic of China.
| |
Collapse
|
3
|
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
|
4
|
Siani P, Frigerio G, Donadoni E, Di Valentin C. Modeling Zeta Potential for Nanoparticles in Solution: Water Flexibility Matters. THE JOURNAL OF PHYSICAL CHEMISTRY. C, NANOMATERIALS AND INTERFACES 2023; 127:9236-9247. [PMID: 37223652 PMCID: PMC10201526 DOI: 10.1021/acs.jpcc.2c08988] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 12/23/2022] [Revised: 04/17/2023] [Indexed: 05/25/2023]
Abstract
Nonequilibrium molecular dynamics simulations were performed to study the electrokinetic properties of five mainstream TIPxP water models (namely, TIP3P-FB, TIP3Pm, TIP4P-FB, TIP4P-Ew, and TIP4P/2005) in NaCl aqueous solutions in the presence of a negatively charged TiO2 surface. The impact of solvent flexibility and system geometry on the electro-osmotic (EO) mobility and flow direction was systematically assessed and compared. We found that lack of water flexibility decelerates the forward EO flow of aqueous solutions at moderate (0.15 M) or high (0.30 M) NaCl concentrations, in some special cases to such an extent that EO flow reversal occurs. Zeta potential (ZP) values were then determined from the bulk EO mobilities using the Helmholtz-Smoluchowski formula. The straight comparison against available experimental data strongly suggests that water flexibility improves the ZP determination of NaCl solutions adjacent to a realistic TiO2 surface under neutral pH conditions.
Collapse
Affiliation(s)
- Paulo Siani
- Dipartimento
di Scienza dei Materiali, Università
di Milano Bicocca, via
R. Cozzi 55, 20125 Milano, Italy
- BioNanoMedicine
Center NANOMIB, University of Milano-Bicocca, 20126 Milano, Italy
| | - Giulia Frigerio
- Dipartimento
di Scienza dei Materiali, Università
di Milano Bicocca, via
R. Cozzi 55, 20125 Milano, Italy
- BioNanoMedicine
Center NANOMIB, University of Milano-Bicocca, 20126 Milano, Italy
| | - Edoardo Donadoni
- Dipartimento
di Scienza dei Materiali, Università
di Milano Bicocca, via
R. Cozzi 55, 20125 Milano, Italy
- BioNanoMedicine
Center NANOMIB, University of Milano-Bicocca, 20126 Milano, Italy
| | - Cristiana Di Valentin
- Dipartimento
di Scienza dei Materiali, Università
di Milano Bicocca, via
R. Cozzi 55, 20125 Milano, Italy
- BioNanoMedicine
Center NANOMIB, University of Milano-Bicocca, 20126 Milano, Italy
| |
Collapse
|
5
|
Weldon R, Wang F. Simulating a flexible water model as rigid: Best practices and lessons learned. J Chem Phys 2023; 158:134506. [PMID: 37031157 PMCID: PMC10076064 DOI: 10.1063/5.0143836] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/26/2023] [Accepted: 03/09/2023] [Indexed: 03/14/2023] Open
Abstract
Two ways to create rigid versions of flexible models are explored. The rigid model can assume the Model's Geometry (MG) as if the molecule is not interacting with any other molecules or the ensemble averaged geometry (EG) under a particular thermodynamic condition. Although the MG model is more straightforward to create, it leads to relatively poor performance. The EG model behaves similarly to the corresponding flexible model (the FL model) and, in some cases, agrees even better with experiments. While the difference between the EG and the FL models is mostly a result of flexibility, the MG and EG models have different dipole moments as a result of an effective induction in the condensed phase. For the three water models studied, the property that shows the most difference is the temperature dependence of density. The MG version of the water model by adaptive force matching for ice and liquid does not possess a temperature of maximum density, which is attributed to a downshift of the putative liquid-liquid phase transition line, leading to the hypothesized second critical point of liquid water to manifest at negative pressure. A new three-phase coexistence method for determining the melting temperature of ice is also presented.
Collapse
Affiliation(s)
- Raymond Weldon
- Department of Chemistry and Biochemistry, University of Arkansas, Fayetteville, Arkansas 72701, USA
| | - Feng Wang
- Department of Chemistry and Biochemistry, University of Arkansas, Fayetteville, Arkansas 72701, USA
| |
Collapse
|
6
|
Hong SN, Ri JH, Mun SY, Yu CJ. Revealing the influence of porosity and temperature on transport properties of nanobubble solution with molecular dynamics simulations. J Mol Liq 2022. [DOI: 10.1016/j.molliq.2022.120518] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022]
|
7
|
Grigoriev A, Kuzovkov Y, Markov I. Bulk viscosity of hydrocarbon solutions at extreme state parameters. I. Linear alkane solutions (C6H14-C16H34). J Mol Liq 2022. [DOI: 10.1016/j.molliq.2021.118328] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
|
8
|
A Facile Strategy to Prepare Small Water Clusters via Interacting with Functional Molecules. Int J Mol Sci 2021; 22:ijms22158250. [PMID: 34361016 PMCID: PMC8347634 DOI: 10.3390/ijms22158250] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/28/2021] [Revised: 07/17/2021] [Accepted: 07/28/2021] [Indexed: 12/03/2022] Open
Abstract
Although small water clusters (SWCs) are important in many research fields, efficient methods of preparing SWCs are still rarely reported, which is mainly due to the lack of related materials and understanding of the molecular interaction mechanisms. In this study, a series of functional molecules were added in water to obtain small water cluster systems. The decreasing rate of the half-peak width in a sodium dodecyl sulfate (SDS)–water system reaches ≈20% at 0.05 mM from 17O nuclear magnetic resonance (NMR) results. Based on density functional theory (DFT) and molecular dynamics (MD) simulation calculation, it can be concluded that functional molecules with stronger negative electrostatic potential (ESP) and higher hydrophilicity have a stronger ability to destroy big water clusters. Notably, the concentrations of our selected molecule systems are one to two magnitudes lower than that of previous reports. This study provides a promising way to optimize aqueous systems in various fields such as oilfield development, protein stability, and metal anti-corrosion.
Collapse
|
9
|
Does uranyl-TBP complex formation happen at the aqueous-organic interface? Revelation by molecular dynamics simulations. J Mol Liq 2021. [DOI: 10.1016/j.molliq.2021.115621] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
10
|
Hashimoto K, Amano KI, Nishi N, Onishi H, Sakka T. Comparison of atomic force microscopy force curve and solvation structure studied by integral equation theory. J Chem Phys 2021; 154:164702. [PMID: 33940841 DOI: 10.1063/5.0046600] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Atomic force microscopy can observe structures of liquids (solvents) on solid surfaces as oscillating force curves. The oscillation originates from the solvation force, which is affected by the interaction between the probe, substrate, and solvents. To investigate the effects of the interactions on the force curve, we calculated the force curves by integral equation theory with various probe and substrate conditions. The probe solvophilicity affected the force curves more than the substrate solvophilicity in our calculation, and its reason is qualitatively explained by the amount of the desolvated solvents. We evaluated the probes and parameters in terms of the qualitative estimation of the number density distribution of the solvent on the wall. The negative of the force curve's derivative with respect to the surface separation reflected the number density distribution better than the force curve. This parameter is based on the method that is proposed previously by Amano et al. [Phys. Chem. Chem. Phys. 18, 15534 (2016)]. The normalized frequency shift can also be used for the qualitative estimation of the number density distribution if the cantilever amplitude is small. Solvophobic probes reflected the number density distribution better than the solvophilic probes. Solvophilic probes resulted in larger oscillation amplitudes than solvophobic probes and are suitable for measurements with a high S/N ratio.
Collapse
Affiliation(s)
- Kota Hashimoto
- Department of Energy and Hydrocarbon Chemistry, Graduate School of Engineering, Kyoto University, Kyoto 615-8510, Japan
| | - Ken-Ichi Amano
- Faculty of Agriculture, Meijo University, 1-501 Shiogamaguchi, Tenpaku, Nagoya 468-8502, Japan
| | - Naoya Nishi
- Department of Energy and Hydrocarbon Chemistry, Graduate School of Engineering, Kyoto University, Kyoto 615-8510, Japan
| | - Hiroshi Onishi
- Department of Chemistry, Graduate School of Science, Kobe University, Nada, Kobe, Hyogo 657-8501, Japan
| | - Tetsuo Sakka
- Department of Energy and Hydrocarbon Chemistry, Graduate School of Engineering, Kyoto University, Kyoto 615-8510, Japan
| |
Collapse
|
11
|
Raucci U, Perrella F, Donati G, Zoppi M, Petrone A, Rega N. Ab-initio molecular dynamics and hybrid explicit-implicit solvation model for aqueous and nonaqueous solvents: GFP chromophore in water and methanol solution as case study. J Comput Chem 2020; 41:2228-2239. [PMID: 32770577 DOI: 10.1002/jcc.26384] [Citation(s) in RCA: 16] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/16/2020] [Revised: 06/21/2020] [Accepted: 06/27/2020] [Indexed: 12/15/2022]
Abstract
Solute-solvent interactions are proxies for understanding how the electronic density of a chromophore interacts with the environment in a more exhaustive way. The subtle balance between polarization, electrostatic, and non-bonded interactions need to be accurately described to obtain good agreement between simulations and experiments. First principles approaches providing accurate configurational sampling through molecular dynamics may be a suitable choice to describe solvent effects on solute chemical-physical properties and spectroscopic features, such as optical absorption of dyes. In this context, accurate energy potentials, obtained by hybrid implicit/explicit solvation methods along with employing nonperiodic boundary conditions, are required to represent bulk solvent around a large solute-solvent cluster. In this work, a novel strategy to simulate methanol solutions is proposed combining ab initio molecular dynamics, a hybrid implicit/explicit flexible solvent model, nonperiodic boundary conditions, and time dependent density functional theory. As case study, the robustness of the proposed protocol has been gauged by investigating the microsolvation and electronic absorption of the anionic green fluorescent protein chromophore in methanol and aqueous solution. Satisfactory results are obtained, reproducing the microsolvation layout of the chromophore and, as a consequence, the experimental trends shown by the optical absorption in different solvents.
Collapse
Affiliation(s)
- Umberto Raucci
- Dipartimento di Scienze Chimiche, Università di Napoli Federico II, Complesso Universitario di M.S. Angelo, Naples, Italy
| | - Fulvio Perrella
- Dipartimento di Scienze Chimiche, Università di Napoli Federico II, Complesso Universitario di M.S. Angelo, Naples, Italy
| | - Greta Donati
- Dipartimento di Scienze Chimiche, Università di Napoli Federico II, Complesso Universitario di M.S. Angelo, Naples, Italy.,Dipartimento di Chimica e Biologia "Adolfo Zambelli", Università di Salerno, Fisciano, Italy
| | - Maria Zoppi
- Dipartimento di Scienze Chimiche, Università di Napoli Federico II, Complesso Universitario di M.S. Angelo, Naples, Italy
| | - Alessio Petrone
- Dipartimento di Scienze Chimiche, Università di Napoli Federico II, Complesso Universitario di M.S. Angelo, Naples, Italy
| | - Nadia Rega
- Dipartimento di Scienze Chimiche, Università di Napoli Federico II, Complesso Universitario di M.S. Angelo, Naples, Italy.,Center for Advanced Biomaterials for Healthcare@CRIB, Naples, Italy
| |
Collapse
|
12
|
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
|
13
|
Abstract
An equation of state for the bulk viscosity of liquid noble gases is proposed. On the basis of dedicated equilibrium molecular dynamics simulations, a multi-mode relaxation ansatz is used to obtain precise bulk viscosity data over a wide range of liquid states. From this dataset, the equation of state emerges as a two-parametric power function with both parameters showing a conspicuous saturation behavior over temperature. After passing a temperature threshold, the bulk viscosity is found to vary significantly over density, a behavior that resembles the frequency response of a one pole low-pass filter. The proposed equation of state is in good agreement with available experimental sound attenuation data.
Collapse
Affiliation(s)
- René Spencer Chatwell
- Thermodynamics and Process Engineering, Technische Universität Berlin, 10587 Berlin, Germany
| | - Jadran Vrabec
- Thermodynamics and Process Engineering, Technische Universität Berlin, 10587 Berlin, Germany
| |
Collapse
|
14
|
Baz J, Hansen N, Gross J. Transferable Anisotropic Mie-Potential Force Field for n-Alcohols: Static and Dynamic Fluid Properties of Pure Substances and Binary Mixtures. Ind Eng Chem Res 2019. [DOI: 10.1021/acs.iecr.9b05323] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Jörg Baz
- Institute of Thermodynamics and Thermal Process Engineering, University of Stuttgart, 70569 Stuttgart, Germany
| | - Niels Hansen
- Institute of Thermodynamics and Thermal Process Engineering, University of Stuttgart, 70569 Stuttgart, Germany
| | - Joachim Gross
- Institute of Thermodynamics and Thermal Process Engineering, University of Stuttgart, 70569 Stuttgart, Germany
| |
Collapse
|
15
|
Celebi AT, Vlugt TJH, Moultos OA. Structural, Thermodynamic, and Transport Properties of Aqueous Reline and Ethaline Solutions from Molecular Dynamics Simulations. J Phys Chem B 2019; 123:11014-11025. [PMID: 31794220 PMCID: PMC6935864 DOI: 10.1021/acs.jpcb.9b09729] [Citation(s) in RCA: 48] [Impact Index Per Article: 9.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/25/2023]
Abstract
Deep eutectic solvents (DESs) are a new generation of green solvents, which are considered an environmentally friendly alternative to ionic liquids and volatile organic compounds. The addition of controlled amounts of water to DESs has a significant effect on their microscopic structure and thus on their thermodynamic and transport properties. In this way, DESs can be modified, leading to solvents with improved characteristics. In this work, molecular dynamics (MD) simulations are performed to obtain a better understanding of the relation between the microscopic structure, molecular interactions, and thermophysical properties of aqueous reline and ethaline solutions at temperatures ranging from 303.15 to 363.15 K. For both reline and ethaline solutions, the hydrogen bond (HB) networks disappear with increasing mass fraction of water, and the intensity of radial distribution function (RDF) peaks decreases. For a mass fraction of water of 40%, most of the HBs between the compounds of reline and ethaline are broken, and DESs are fully dissolved in water. Consequently, a monotonic decrease in viscosities and an increase in self-diffusion coefficients are observed. Ionic conductivities show a nonmonotonic behavior with increasing water content. Up to 60% water mass fraction, the ionic conductivities increase with increasing water content. A further increase in the mass fraction of water decreases conductivities. For all studied systems, the HB network and the peaks of RDFs show relatively small changes for water mass fractions below 5% and beyond 40%. The MD results show that viscosities decrease with temperature, while diffusivities and ionic conductivities increase. The effect of the temperature on the structure of DES-water mixtures is negligible.
Collapse
Affiliation(s)
- Alper T Celebi
- Engineering Thermodynamics, Process & Energy Department, Faculty of Mechanical, Maritime and Materials Engineering , Delft University of Technology , Leeghwaterstraat 39 , 2628CB Delft , The Netherlands
| | - Thijs J H Vlugt
- Engineering Thermodynamics, Process & Energy Department, Faculty of Mechanical, Maritime and Materials Engineering , Delft University of Technology , Leeghwaterstraat 39 , 2628CB Delft , The Netherlands
| | - Othonas A Moultos
- Engineering Thermodynamics, Process & Energy Department, Faculty of Mechanical, Maritime and Materials Engineering , Delft University of Technology , Leeghwaterstraat 39 , 2628CB Delft , The Netherlands
| |
Collapse
|
16
|
Celebi AT, Nguyen CT, Hartkamp R, Beskok A. The role of water models on the prediction of slip length of water in graphene nanochannels. J Chem Phys 2019; 151:174705. [DOI: 10.1063/1.5123713] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Alper Tunga Celebi
- Process and Energy Department, Delft University of Technology, Leeghwaterstraat 39, 2628 CB Delft, The Netherlands
| | - Chinh Thanh Nguyen
- Lyle School of Engineering, Southern Methodist University, 3101 Dyer Street, Dallas, Texas 75205, USA
| | - Remco Hartkamp
- Process and Energy Department, Delft University of Technology, Leeghwaterstraat 39, 2628 CB Delft, The Netherlands
| | - Ali Beskok
- Lyle School of Engineering, Southern Methodist University, 3101 Dyer Street, Dallas, Texas 75205, USA
| |
Collapse
|
17
|
Mendis CH, Piskulich ZA, Thompson WH. Tests of the Stokes–Einstein Relation through the Shear Viscosity Activation Energy of Water. J Phys Chem B 2019; 123:5857-5865. [DOI: 10.1021/acs.jpcb.9b04647] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/22/2023]
Affiliation(s)
- Camina H. Mendis
- Department of Chemistry, University of Kansas, Lawrence, Kansas 66045, United States
| | - Zeke A. Piskulich
- 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
|
18
|
Baer A, Miličević Z, Smith DM, Smith AS. Water in an electric field does not dance alone: The relation between equilibrium structure, time dependent viscosity and molecular motions. J Mol Liq 2019. [DOI: 10.1016/j.molliq.2019.02.055] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/31/2023]
|
19
|
Jamali SH, Westen TV, Moultos OA, Vlugt TJH. Optimizing Nonbonded Interactions of the OPLS Force Field for Aqueous Solutions of Carbohydrates: How to Capture Both Thermodynamics and Dynamics. J Chem Theory Comput 2018; 14:6690-6700. [PMID: 30407814 PMCID: PMC6293444 DOI: 10.1021/acs.jctc.8b00909] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Abstract
![]()
Knowledge on thermodynamic and transport
properties of aqueous
solutions of carbohydrates is of great interest for process and product
design in the food, pharmaceutical, and biotechnological industries.
Molecular simulation is a powerful tool to calculate these properties,
but current classical force fields cannot provide accurate estimates
for all properties of interest. The poor performance of the force
fields is mainly observed for concentrated solutions, where solute–solute
interactions are overestimated. In this study, we propose a method
to refine force fields, such that solute–solute interactions
are more accurately described. The OPLS force field combined with
the SPC/Fw water model is used as a basis. We scale the nonbonded
interaction parameters of sucrose, a disaccharide. The scaling factors
are chosen in such a way that experimental thermodynamic and transport
properties of aqueous solutions of sucrose are accurately reproduced.
Using a scaling factor of 0.8 for Lennard-Jones energy parameters
(ϵ) and a scaling factor of 0.95 for partial atomic charges
(q), we find excellent agreement between experiments
and computed liquid densities, thermodynamic factors, shear viscosities,
self-diffusion coefficients, and Fick (mutual) diffusion coefficients.
The transferability of these optimum scaling factors to other carbohydrates
is verified by computing thermodynamic and transport properties of
aqueous solutions of d-glucose, a monosaccharide. The good
agreement between computed properties and experiments suggests that
the scaled interaction parameters are transferable to other carbohydrates,
especially for concentrated solutions.
Collapse
Affiliation(s)
- Seyed Hossein Jamali
- Engineering Thermodynamics, Process & Energy Department, Faculty of Mechanical, Maritime and Materials Engineering , Delft University of Technology , Leeghwaterstraat 39 , 2628CB Delft , The Netherlands
| | - Thijs van Westen
- Institute AMOLF , Science Park 104 , 1098XG , Amsterdam , The Netherlands.,Institute of Thermodynamics and Thermal Process Engineering , University of Stuttgart , Pfaffenwaldring 9 , D-70569 Stuttgart , Germany
| | - Othonas A Moultos
- Engineering Thermodynamics, Process & Energy Department, Faculty of Mechanical, Maritime and Materials Engineering , Delft University of Technology , Leeghwaterstraat 39 , 2628CB Delft , The Netherlands
| | - Thijs J H Vlugt
- Engineering Thermodynamics, Process & Energy Department, Faculty of Mechanical, Maritime and Materials Engineering , Delft University of Technology , Leeghwaterstraat 39 , 2628CB Delft , The Netherlands
| |
Collapse
|
20
|
In-situ discrimination of the water cluster size distribution in aqueous solution by ToF-SIMS. Sci China Chem 2017. [DOI: 10.1007/s11426-017-9180-1] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
|
21
|
Medina JS, Arismendi-Arrieta DJ, Alemán JV, Prosmiti R. Developing time to frequency-domain descriptors for relaxation processes: Local trends. J Mol Liq 2017. [DOI: 10.1016/j.molliq.2017.08.034] [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]
|
22
|
Köhler MH, Bordin JR, da Silva LB, Barbosa MC. Breakdown of the Stokes–Einstein water transport through narrow hydrophobic nanotubes. Phys Chem Chem Phys 2017; 19:12921-12927. [DOI: 10.1039/c7cp02058a] [Citation(s) in RCA: 34] [Impact Index Per Article: 4.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
As water density is increased inside narrow hydrophobic nanotubes, the viscosity shows a huge increase associated with a small increase in the diffusion, which violates the Stokes–Einstein relation.
Collapse
Affiliation(s)
- Mateus Henrique Köhler
- Instituto de Física
- Universidade Federal do Rio Grande do Sul
- Caixa Postal 15051
- Porto Alegre
- Brazil
| | - José Rafael Bordin
- Campus Caapava do Sul
- Universidade Federal do Pampa
- v. Pedro Anunciacao 111
- CEP 96570-000
- Ca apava do Sul
| | | | - Marcia C. Barbosa
- Instituto de Física
- Universidade Federal do Rio Grande do Sul
- Caixa Postal 15051
- Porto Alegre
- Brazil
| |
Collapse
|
23
|
Parameterization and optimization of the menthol force field for molecular dynamics simulations. J Mol Model 2016; 22:234. [PMID: 27604277 PMCID: PMC5014899 DOI: 10.1007/s00894-016-3082-1] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/14/2015] [Accepted: 08/05/2016] [Indexed: 11/30/2022]
Abstract
Menthol’s various biological properties render it a useful component for medical and cosmetological applications, while its three centers of asymmetry mean that it can be used in a range of organic reactions. Menthol-substituted ionic liquids (ILs) have been found to exhibit promising antimicrobial and antielectrostatic properties, as well as being useful in organic catalysis and biochemical studies. However, so far, a force field designed and validated specifically for the menthol molecule has not been constructed. In the present work, the validation and optimization of force field parameters with regard to the ability to reproduce the macroscopic properties of menthol is presented. The set of optimized potentials for liquid simulations all atom (OPLS-AA) compatible parameters was tested and carefully tuned. The refinement of parameters included fitting of partial atomic charges, optimization of Lennard-Jones parameters, and recalculation of the dihedral angle parameters needed to reproduce quantum energy profiles. To validate the force field, a variety of physicochemical properties were calculated for liquid menthol. Both thermodynamic and kinetic properties were taken into account, including density, surface tension, enthalpy of vaporization, and shear viscosity. The obtained force field was proven to accurately reproduce the properties of the investigated compound while being fully compatible with the OPLS-AA force field.
Collapse
|
24
|
Zong D, Hu H, Duan Y, Sun Y. Viscosity of Water under Electric Field: Anisotropy Induced by Redistribution of Hydrogen Bonds. J Phys Chem B 2016; 120:4818-27. [DOI: 10.1021/acs.jpcb.6b01686] [Citation(s) in RCA: 23] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Diyuan Zong
- Department
of Thermal Engineering, Tsinghua University, Beijing 100084, P. R. China
| | - Han Hu
- Department
of Mechanical Engineering and Mechanics, Drexel University, Philadelphia, Pennsylvania 19104, United States
| | - Yuanyuan Duan
- Department
of Thermal Engineering, Tsinghua University, Beijing 100084, P. R. China
| | - Ying Sun
- Department
of Thermal Engineering, Tsinghua University, Beijing 100084, P. R. China
- Department
of Mechanical Engineering and Mechanics, Drexel University, Philadelphia, Pennsylvania 19104, United States
| |
Collapse
|
25
|
Sahu P, Ali SM, Shenoy KT. Passage of TBP–uranyl complexes from aqueous–organic interface to the organic phase: insights from molecular dynamics simulation. Phys Chem Chem Phys 2016; 18:23769-84. [DOI: 10.1039/c6cp02194h] [Citation(s) in RCA: 30] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Water/organic interface representing TBP orientation for neutral versus acidic interface and occurrence of UO22+–TBP–NO3− species in various stoichiometry.
Collapse
Affiliation(s)
- Pooja Sahu
- Chemical Engineering Division
- Bhabha Atomic Research Centre
- Mumbai-400085
- India
| | - Sk. Musharaf Ali
- Chemical Engineering Division
- Bhabha Atomic Research Centre
- Mumbai-400085
- India
| | | |
Collapse
|
26
|
Nieto-Draghi C, Fayet G, Creton B, Rozanska X, Rotureau P, de Hemptinne JC, Ungerer P, Rousseau B, Adamo C. A General Guidebook for the Theoretical Prediction of Physicochemical Properties of Chemicals for Regulatory Purposes. Chem Rev 2015; 115:13093-164. [PMID: 26624238 DOI: 10.1021/acs.chemrev.5b00215] [Citation(s) in RCA: 70] [Impact Index Per Article: 7.8] [Reference Citation Analysis] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Affiliation(s)
- Carlos Nieto-Draghi
- IFP Energies nouvelles , 1 et 4 avenue de Bois-Préau, 92852 Rueil-Malmaison, France
| | - Guillaume Fayet
- INERIS, Parc Technologique Alata, BP2 , 60550 Verneuil-en-Halatte, France
| | - Benoit Creton
- IFP Energies nouvelles , 1 et 4 avenue de Bois-Préau, 92852 Rueil-Malmaison, France
| | - Xavier Rozanska
- Materials Design S.A.R.L. , 18, rue de Saisset, 92120 Montrouge, France
| | - Patricia Rotureau
- INERIS, Parc Technologique Alata, BP2 , 60550 Verneuil-en-Halatte, France
| | | | - Philippe Ungerer
- Materials Design S.A.R.L. , 18, rue de Saisset, 92120 Montrouge, France
| | - Bernard Rousseau
- Laboratoire de Chimie-Physique, Université Paris Sud , UMR 8000 CNRS, Bât. 349, 91405 Orsay Cedex, France
| | - Carlo Adamo
- Institut de Recherche Chimie Paris, PSL Research University, CNRS, Chimie Paristech , 11 rue P. et M. Curie, F-75005 Paris, France.,Institut Universitaire de France , 103 Boulevard Saint Michel, F-75005 Paris, France
| |
Collapse
|
27
|
Pham TT, Lemaire T, Capiez-Lernout E, Lewerenz M, To QD, Christie JK, Di Tommaso D, de Leeuw NH, Naili S. Properties of water confined in hydroxyapatite nanopores as derived from molecular dynamics simulations. Theor Chem Acc 2015. [DOI: 10.1007/s00214-015-1653-3] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/23/2022]
|
28
|
Arismendi-Arrieta D, Medina JS, Fanourgakis GS, Prosmiti R, Delgado-Barrio G. Simulating liquid water for determining its structural and transport properties. Appl Radiat Isot 2014; 83 Pt B:115-21. [DOI: 10.1016/j.apradiso.2013.01.020] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/20/2012] [Revised: 01/10/2013] [Accepted: 01/10/2013] [Indexed: 11/29/2022]
|
29
|
Kuang S, Gezelter JD. Velocity shearing and scaling RNEMD: a minimally perturbing method for simulating temperature and momentum gradients. Mol Phys 2012. [DOI: 10.1080/00268976.2012.680512] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/28/2022]
|
30
|
Villarreal P, Miret-Artés S, Roncero O, Campos-Martínez J. New trends in atomic and molecular clusters (in honour of Gerardo Delgado-Barrio). Chem Phys 2012. [DOI: 10.1016/j.chemphys.2011.10.012] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/15/2022]
|
31
|
Fanourgakis GS, Medina JS, Prosmiti R. Determining the Bulk Viscosity of Rigid Water Models. J Phys Chem A 2012; 116:2564-70. [DOI: 10.1021/jp211952y] [Citation(s) in RCA: 42] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
| | - J. S. Medina
- Instituto de Física Fundamental (CSIC), Serrano 123, 28006 Madrid, Spain
| | - R. Prosmiti
- Instituto de Física Fundamental (CSIC), Serrano 123, 28006 Madrid, Spain
| |
Collapse
|