1
|
Chaboksavar M, Soltanabadi A. Experimental and computational study of binary mixture ethanolamine and 2-amino-2-methyl-1-propanol. J Mol Graph Model 2024; 131:108816. [PMID: 38909381 DOI: 10.1016/j.jmgm.2024.108816] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/26/2024] [Revised: 06/19/2024] [Accepted: 06/19/2024] [Indexed: 06/25/2024]
Abstract
The present work involves experimental and computational investigations into the density of pure and mixed states of ethanolamine (ET) and 2-amino-2-methyl-1-propanol (AMP) under a pressure of 1 atm and temperatures ranging from 293.15 K to 333.15 K The density data were used to derive the excess molar volume, thermal expansion coefficient, and isothermal coefficient of pressure excess molar enthalpy. The Redlich-Kister equation was employed to calculate the excess molar and its accompanying coefficients. In the gas phase, density functional theory (DFT) was utilized to explore the most stable structures of ET … ET, AMP … AMP, and the ET … AMP mixture. Molecular dynamics simulation (MD) was used to calculate the structural properties of these mixtures in the liquid phase. Radial distribution function (RDFs) combined distribution function (CDF) and spatial distribution function (SDF) in different mole fractions calculated in the liquid phase. The intramolecular and intermolecular interactions of ethanolamine and AMP were obtained using the radial distribution function in different molar fractions. It was found that the ethanolamine molecule has a greater tendency to form intramolecular hydrogen bonds, while the AMP molecule has a greater tendency to form intermolecular hydrogen bonds.
Collapse
Affiliation(s)
- Mastaneh Chaboksavar
- Department of Physical Chemistry, Faculty of Chemistry, Razi University, Kermanshah, 6714414971, Iran.
| | - Azim Soltanabadi
- Department of Physical Chemistry, Faculty of Chemistry, Razi University, Kermanshah, 6714414971, Iran.
| |
Collapse
|
2
|
Hulikal Chakrapani T, Hajibeygi H, Moultos OA, Vlugt TJH. Mutual Diffusivities of Mixtures of Carbon Dioxide and Hydrogen and Their Solubilities in Brine: Insight from Molecular Simulations. Ind Eng Chem Res 2024; 63:10456-10481. [PMID: 38882502 PMCID: PMC11177264 DOI: 10.1021/acs.iecr.4c01078] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/21/2024] [Revised: 05/14/2024] [Accepted: 05/20/2024] [Indexed: 06/18/2024]
Abstract
H2-CO2 mixtures find wide-ranging applications, including their growing significance as synthetic fuels in the transportation industry, relevance in capture technologies for carbon capture and storage, occurrence in subsurface storage of hydrogen, and hydrogenation of carbon dioxide to form hydrocarbons and alcohols. Here, we focus on the thermodynamic properties of H2-CO2 mixtures pertinent to underground hydrogen storage in depleted gas reservoirs. Molecular dynamics simulations are used to compute mutual (Fick) diffusivities for a wide range of pressures (5 to 50 MPa), temperatures (323.15 to 423.15 K), and mixture compositions (hydrogen mole fraction from 0 to 1). At 5 MPa, the computed mutual diffusivities agree within 5% with the kinetic theory of Chapman and Enskog at 423.15 K, albeit exhibiting deviations of up to 25% between 323.15 and 373.15 K. Even at 50 MPa, kinetic theory predictions match computed diffusivities within 15% for mixtures comprising over 80% H2 due to the ideal-gas-like behavior. In mixtures with higher concentrations of CO2, the Moggridge correlation emerges as a dependable substitute for the kinetic theory. Specifically, when the CO2 content reaches 50%, the Moggridge correlation achieves predictions within 10% of the computed Fick diffusivities. Phase equilibria of ternary mixtures involving CO2-H2-NaCl were explored using Gibbs Ensemble (GE) simulations with the Continuous Fractional Component Monte Carlo (CFCMC) technique. The computed solubilities of CO2 and H2 in NaCl brine increased with the fugacity of the respective component but decreased with NaCl concentration (salting out effect). While the solubility of CO2 in NaCl brine decreased in the ternary system compared to the binary CO2-NaCl brine system, the solubility of H2 in NaCl brine increased less in the ternary system compared to the binary H2-NaCl brine system. The cooperative effect of H2-CO2 enhances the H2 solubility while suppressing the CO2 solubility. The water content in the gas phase was found to be intermediate between H2-NaCl brine and CO2-NaCl brine systems. Our findings have implications for hydrogen storage and chemical technologies dealing with CO2-H2 mixtures, particularly where experimental data are lacking, emphasizing the need for reliable thermodynamic data on H2-CO2 mixtures.
Collapse
Affiliation(s)
- Thejas Hulikal Chakrapani
- Reservoir Engineering, Geoscience and Engineering Department, Faculty of Civil Engineering and Geosciences, Delft University of Technology, Delft 2628 CN, The Netherlands
| | - Hadi Hajibeygi
- Reservoir Engineering, Geoscience and Engineering Department, Faculty of Civil Engineering and Geosciences, Delft University of Technology, Delft 2628 CN, The Netherlands
| | - Othonas A Moultos
- Engineering Thermodynamics, Process and Energy Department, Faculty of Mechanical Engineering, Delft University of Technology, Delft 2628 CB, The Netherlands
| | - Thijs J H Vlugt
- Engineering Thermodynamics, Process and Energy Department, Faculty of Mechanical Engineering, Delft University of Technology, Delft 2628 CB, The Netherlands
| |
Collapse
|
3
|
Chaimovich M, Chaimovich A. Relative Resolution: An Analysis with the Kullback-Leibler Entropy. J Chem Theory Comput 2024; 20:2074-2087. [PMID: 38416535 DOI: 10.1021/acs.jctc.3c01052] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/29/2024]
Abstract
A novel type of a multiscale approach, called Relative Resolution (RelRes), can correctly retrieve the behavior of various nonpolar liquids while speeding up molecular simulations by almost an order of magnitude. In this approach in a single system, molecules switch their resolution in terms of their relative separation, with near neighbors interacting via fine-grained potentials, yet far neighbors interacting via coarse-grained potentials; notably, these two potentials are analytically parametrized by a multipole approximation. Our current work focuses on analyzing RelRes by relating it with the Kullback-Leibler (KL) entropy, which is a useful metric for multiscale errors. In particular, we thoroughly examine the exact and approximate versions of this informatic measure for several alkane systems. By analyzing its dependency on the system size, we devise a formula for predicting the exact KL entropy of an "infinite" system via the computation of the approximate KL entropy of an "infinitesimal" system. Demonstrating that the KL entropy can holistically capture many multiscale errors, we settle bounds for the KL entropy that ensure a sufficient representation of the structural and thermal behavior by the RelRes algorithm. This, in turn, allows the scientific community to readily determine the ideal switching distance for an arbitrary RelRes system.
Collapse
Affiliation(s)
- Mark Chaimovich
- Russian School of Mathematics, North Bethesda, Maryland 20852, United States
| | - Aviel Chaimovich
- Department of Chemical and Biological Engineering, Drexel University, Philadelphia, Pennsylvania 19104, United States
| |
Collapse
|
4
|
Fields B, Schelling PK. Connection between partial pressure, volatility, and the Soret effect elucidated using simulations of nonideal supercritical fluid mixtures. J Chem Phys 2024; 160:084501. [PMID: 38385515 DOI: 10.1063/5.0185603] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/31/2023] [Accepted: 01/31/2024] [Indexed: 02/23/2024] Open
Abstract
Building on recent simulation work, it is demonstrated using molecular dynamics simulations of two-component fluid mixtures that the chemical contribution to the Soret effect in two-component nonideal fluid mixtures arises due to differences in how the partial pressures of the components respond to temperature and density gradients. Further insight is obtained by reviewing the connection between activity and deviations from Raoult's law in the measurement of the vapor pressure of a liquid mixture. A new parameter γsS, defined in a manner similar to the activity coefficient, is used to characterize differences deviations from "ideal" behavior. It is then shown that the difference γ2S-γ1S is predictive of the sign of the Soret coefficient and is correlated to its magnitude. We hence connect the Soret effect to the relative volatility of the components of a fluid mixture, with the more volatile component enriched in the low-density, high-temperature region, and the less volatile component enriched in the high-density, low-temperature region. Because γsS is closely connected to the activity coefficient, this suggests the possibility that measurement of partial vapor pressures might be used to indirectly determine the Soret coefficient. It is proposed that the insight obtained here is quite general and should be applicable to a wide range of materials systems. An attempt is made to understand how these results might apply to other materials systems including interstitials in solids and multicomponent solids with interdiffusion occurring via a vacancy mechanism.
Collapse
Affiliation(s)
- Brandon Fields
- Department of Physics, University of Central Florida, Orlando, Florida 32816-2385, USA
| | - Patrick K Schelling
- Department of Physics, University of Central Florida, Orlando, Florida 32816-2385, USA
- Advanced Materials Processing and Analysis Center, University of Central Florida, Orlando, Florida 32816-2385, USA
- Renewable Energy and Chemical Transformations (REACT) Cluster, University of Central Florida, Orlando, Florida 32816-2385, USA
| |
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
|
Shoemaker BA, Haji-Akbari A. Ideal conductor/dielectric model (ICDM): A generalized technique to correct for finite-size effects in molecular simulations of hindered ion transport. J Chem Phys 2024; 160:024116. [PMID: 38197447 DOI: 10.1063/5.0180029] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/06/2023] [Accepted: 12/18/2023] [Indexed: 01/11/2024] Open
Abstract
Molecular simulations serve as indispensable tools for investigating the kinetics and elucidating the mechanism of hindered ion transport across nanoporous membranes. In particular, recent advancements in advanced sampling techniques have made it possible to access translocation timescales spanning several orders of magnitude. In our prior study [Shoemaker et al., J. Chem. Theory Comput. 18, 7142 (2022)], we identified significant finite size artifacts in simulations of pressure-driven hindered ion transport through nanoporous graphitic membranes. We introduced the ideal conductor model, which effectively corrects for such artifacts by assuming the feed to be an ideal conductor. In the present work, we introduce the ideal conductor dielectric model (Icdm), a generalization of our earlier model, which accounts for the dielectric properties of both the membrane and the filtrate. Using the Icdm model substantially enhances the agreement among corrected free energy profiles obtained from systems of varying sizes, with notable improvements observed in regions proximate to the pore exit. Moreover, the model has the capability to consider secondary ion passage events, including the transport of a co-ion subsequent to the traversal of a counter-ion, a feature that is absent in our original model. We also investigate the sensitivity of the new model to various implementation details. The Icdm model offers a universally applicable framework for addressing finite size artifacts in molecular simulations of ion transport. It stands as a significant advancement in our quest to use molecular simulations to comprehensively understand and manipulate ion transport processes through nanoporous membranes.
Collapse
Affiliation(s)
- Brian A Shoemaker
- Department of Chemical and Environmental Engineering, Yale University, New Haven, Connecticut 06520, USA
| | - Amir Haji-Akbari
- Department of Chemical and Environmental Engineering, Yale University, New Haven, Connecticut 06520, USA
| |
Collapse
|
7
|
Olowookere FV, Turner CH. Predicting Gaseous Solute Diffusion in Viscous Multivalent Ionic Liquid Solvents. J Phys Chem B 2023; 127:9144-9154. [PMID: 37831616 DOI: 10.1021/acs.jpcb.3c03858] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/15/2023]
Abstract
Calculating solute diffusion in dense, viscous solvents can be particularly challenging in molecular dynamics simulations due to the long time scales involved. Here, a new scaling approach is developed for predicting solute diffusion based on analyses of CO2 and SO2 diffusion in two different multivalent ionic liquid solvents. Various scaling approaches are initially evaluated, including single and separate thermostats for the solute and solvent, as well as the application of the Arrhenius relationship and the Speedy-Angell power law. A very strong logarithmic correlation is established between the solvent-accessible surface area and solute diffusion. This relationship, reflecting Danckwerts' surface renewal theory and the Vrentas-Duda free volume model, presents a valuable method for estimating diffusion behavior from short simulation trajectories at elevated temperatures. The approach may be beneficial for enhancing predictive modeling in similar challenging systems and should be more broadly evaluated.
Collapse
Affiliation(s)
- Feranmi V Olowookere
- Department of Chemical and Biological Engineering, The University of Alabama, Tuscaloosa, Alabama 35487-0203, United States
| | - C Heath Turner
- Department of Chemical and Biological Engineering, The University of Alabama, Tuscaloosa, Alabama 35487-0203, United States
| |
Collapse
|
8
|
Wei L, Bai Q, Li X, Liu Z, Li C, Cui Y, Shen L, Zhu C, Fang W. Puckered Zigzag Monolayer Ice: Does a Confined Flat Four-Coordinated Monolayer Ice Always Have a Corresponding Puckered Phase? J Phys Chem Lett 2023; 14:8890-8895. [PMID: 37767947 DOI: 10.1021/acs.jpclett.3c02065] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/29/2023]
Abstract
We note that a flat, four-coordinated monolayer ice under confinement always has a corresponding puckered phase. Recently, a monolayer ice consisting of an array of zigzag water chains (ZZMI) predicted by first-principles calculations of water under confinement is a flat four-coordinated monolayer ice. Herein, to investigate whether puckered ZZMI exists stably, we perform molecular dynamics simulations of two-dimensional (2D) ice formation for water constrained in graphene nanocapillaries. We find a novel monolayer ice structure that can be viewed as the ZZMI puckered along the direction perpendicular to the zigzag chain (pZZMI). Unlike ZZMI that does not satisfy the ice rule, each water molecule in pZZMI can form four hydrogen bonds (HBs) via forming two stable intersublayer HBs and two intrasublayer HBs. This work provides a fresh perspective on 2D confined ice, highlighting the intrinsic connections between 2D confined ices.
Collapse
Affiliation(s)
- Laiyang Wei
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, People's Republic of China
| | - Qi Bai
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, People's Republic of China
| | - Xiaojiao Li
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, People's Republic of China
| | - Ziyuan Liu
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, People's Republic of China
| | - Chenruyuan Li
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, People's Republic of China
| | - Yanhong Cui
- College of Chemical Engineering, Zhejiang University of Technology, Hangzhou, Zhejiang 310023, People's Republic of China
| | - Lin Shen
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, People's Republic of China
| | - Chongqin Zhu
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, People's Republic of China
| | - Weihai Fang
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, People's Republic of China
| |
Collapse
|
9
|
Saric D, Guevara-Carrion G, Gaponenko Y, Shevtsova V, Vrabec J. Diffusion of hydrocarbons diluted in supercritical carbon dioxide. Sci Rep 2023; 13:16107. [PMID: 37752219 PMCID: PMC10522683 DOI: 10.1038/s41598-023-42892-7] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/24/2023] [Accepted: 09/15/2023] [Indexed: 09/28/2023] Open
Abstract
Mutual diffusion of six hydrocarbons (methane, ethane, isobutane, benzene, toluene or naphthalene) diluted in supercritical carbon dioxide ([Formula: see text]) is studied by molecular dynamics simulation near the Widom line, i.e., in the temperature range from 290 to 345 K along the isobar 9 MPa. The [Formula: see text] + aromatics mixtures are additionally sampled at 10 and 12 MPa and an experimental database with Fick diffusion coefficient data for those systems is provided. Taylor dispersion experiments of [Formula: see text] with benzene, toluene, n-dodecane and 1,2,3,4-tetrahydronaphthalene are conducted along the [Formula: see text] 10 MPa isobar. Maxwell-Stefan and Fick diffusion coefficients are analyzed, together with the thermodynamic factor that relates them. It is found that the peculiar behavior of the Fick diffusion coefficient of some [Formula: see text] mixtures in the extended critical region is a consequence of the thermodynamic factor minimum due to pronounced clustering on the molecular scale. Further, the strong dependence of the Fick diffusion coefficient on the molecular mass of the solute as well as the breakdown of the Stokes-Einstein relation near the Widom line are confirmed. Eleven correlations for the prediction of the Fick diffusion coefficient of [Formula: see text] mixtures are assessed. An alternative two-step approach for the prediction of the infinite dilution Fick diffusion coefficient of supercritical [Formula: see text] mixtures is proposed. It requires only the state point in terms of temperature and pressure (or density) as well as the molecular solute mass as input parameters. First, entropy scaling is applied to estimate the self-diffusion coefficient of [Formula: see text]. Subsequently, this coefficient is used to determine the infinite dilution Fick diffusion coefficient of the mixture, based on the finding that these two diffusion coefficients exhibit a linear relationship, where the slope depends only on the molecular solute mass.
Collapse
Affiliation(s)
- Denis Saric
- Thermodynamics, Technical University of Berlin, Ernst-Reuter-Platz 1, 10587, Berlin, Germany
| | | | - Yury Gaponenko
- MRC, CP-165/62, Université libre de Bruxelles (ULB), Ave. F.D. Roosevelt 50, B-1050, Brussels, Belgium
| | - Valentina Shevtsova
- Fluid Mechanics Group, Faculty of Engineering, Mondragon University, 20500, Mondragon, Spain
- IKERBASQUE, Basque Foundation for Science, Plaza Euskadi 5, 48009, Bilbao, Spain
| | - Jadran Vrabec
- Thermodynamics, Technical University of Berlin, Ernst-Reuter-Platz 1, 10587, Berlin, Germany.
| |
Collapse
|
10
|
Gittus OR, Bresme F. Mass dipole contribution to the isotopic Soret effect in molecular mixtures. J Chem Phys 2023; 159:114503. [PMID: 37724736 DOI: 10.1063/5.0164253] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/21/2023] [Accepted: 08/03/2023] [Indexed: 09/21/2023] Open
Abstract
Temperature gradients induce mass separation in mixtures in a process called thermal diffusion and are quantified by the Soret coefficient ST. Thermal diffusion in fluid mixtures has been interpreted recently in terms of the so-called (pseudo-)isotopic Soret effect but only considering the mass and moment of inertia differences of the molecules. We demonstrate that the first moment of the molecular mass distribution, the mass dipole, contributes significantly to the isotopic Soret effect. To probe this physical effect, we investigate fluid mixtures consisting of rigid linear molecules that differ only by the first moment of their mass distributions. We demonstrate that such mixtures have non-zero Soret coefficients in contrast with ST = 0 predicted by current formulations. For the isotopic mixtures investigated in this work, the dependence of ST on the mass dipole arises mainly through the thermal diffusion coefficient DT. In turn, DT is correlated with the dependence of the molecular librational modes on the mass dipole. We examine the interplay of the mass dipole and the moment of inertia in defining the isotopic Soret effect and propose empirical equations that include the mass dipole contribution.
Collapse
Affiliation(s)
- Oliver R Gittus
- Department of Chemistry, Molecular Sciences Research Hub, Imperial College London, London W12 0BZ, United Kingdom
| | - Fernando Bresme
- Department of Chemistry, Molecular Sciences Research Hub, Imperial College London, London W12 0BZ, United Kingdom
| |
Collapse
|
11
|
Blazquez S, Abascal JLF, Lagerweij J, Habibi P, Dey P, Vlugt TJH, Moultos OA, Vega C. Computation of Electrical Conductivities of Aqueous Electrolyte Solutions: Two Surfaces, One Property. J Chem Theory Comput 2023; 19:5380-5393. [PMID: 37506381 PMCID: PMC10448725 DOI: 10.1021/acs.jctc.3c00562] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/26/2023] [Indexed: 07/30/2023]
Abstract
In this work, we computed electrical conductivities under ambient conditions of aqueous NaCl and KCl solutions by using the Einstein-Helfand equation. Common force fields (charge q = ±1 e) do not reproduce the experimental values of electrical conductivities, viscosities, and diffusion coefficients. Recently, we proposed the idea of using different charges to describe the potential energy surface (PES) and the dipole moment surface (DMS). In this work, we implement this concept. The equilibrium trajectories required to evaluate electrical conductivities (within linear response theory) were obtained by using scaled charges (with the value q = ±0.75 e) to describe the PES. The potential parameters were those of the Madrid-Transport force field, which accurately describe viscosities and diffusion coefficients of these ionic solutions. However, integer charges were used to compute the conductivities (thus describing the DMS). The basic idea is that although the scaled charge describes the ion-water interaction better, the integer charge reflects the value of the charge that is transported due to the electric field. The agreement obtained with experiments is excellent, as for the first time electrical conductivities (and the other transport properties) of NaCl and KCl electrolyte solutions are described with high accuracy for the whole concentration range up to their solubility limit. Finally, we propose an easy way to obtain a rough estimate of the actual electrical conductivity of the potential model under consideration using the approximate Nernst-Einstein equation, which neglects correlations between different ions.
Collapse
Affiliation(s)
- Samuel Blazquez
- Dpto.
Química Física I, Fac. Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - Jose L. F. Abascal
- Dpto.
Química Física I, Fac. Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - Jelle Lagerweij
- Engineering
Thermodynamics, Process and Energy Department, Faculty of Mechanical,
Maritime and Materials Engineering, Delft
University of Technology, Leeghwaterstraat 39, 2628CB Delft, The Netherlands
| | - Parsa 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
| | - Poulumi Dey
- Department
of Materials Science and Engineering, Faculty of Mechanical, Maritime
and Materials Engineering, Delft University
of Technology, Mekelweg
2, 2628CD Delft, The Netherlands
| | - Thijs 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
| | - Othonas 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
| | - Carlos Vega
- Dpto.
Química Física I, Fac. Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| |
Collapse
|
12
|
Antolović I, Staubach J, Stephan S, Vrabec J. Phase equilibria of symmetric Lennard-Jones mixtures and a look at the transport properties near the upper critical solution temperature. Phys Chem Chem Phys 2023. [PMID: 37365979 DOI: 10.1039/d3cp01434g] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/28/2023]
Abstract
This study investigates phase equilibria and transport properties of five symmetric binary Lennard-Jones mixtures using molecular simulation and equation of state models. The mixtures are selected for their representation of different types of phase behavior and the research contributes to the development of simulation techniques, mixture theories and understanding of thermophysical mixture properties. A novel method is introduced for determining the critical end point (CEP) and critical azeotropic end point (CAEP) by molecular simulation. The van der Waals one-fluid theory is assessed for its performance in conjunction with Lennard-Jones equation of state models, while addressing different phase equilibrium types simultaneously. An empirical correlation is introduced to account for deviations between the equation of state and simulation that arise when using the same binary interaction parameter. This study also investigates the influence of the liquid-liquid critical point on thermophysical properties, which are found to exhibit no significant anomalies or singularities. System-size effects of diffusion coefficients are addressed by extrapolating simulation data to the thermodynamic limit and applying analytical finite-size corrections.
Collapse
Affiliation(s)
- Ivan Antolović
- Thermodynamics, Technical University Berlin, 10587 Berlin, Germany.
| | - Jens Staubach
- Laboratory of Engineering Thermodynamics (LTD), RPTU Kaiserslautern, Kaiserslautern, Germany
| | - Simon Stephan
- Laboratory of Engineering Thermodynamics (LTD), RPTU Kaiserslautern, Kaiserslautern, Germany
| | - Jadran Vrabec
- Thermodynamics, Technical University Berlin, 10587 Berlin, Germany.
| |
Collapse
|
13
|
Halder R, Nissley DA, Sitarik I, Jiang Y, Rao Y, Vu QV, Li MS, Pritchard J, O'Brien EP. How soluble misfolded proteins bypass chaperones at the molecular level. Nat Commun 2023; 14:3689. [PMID: 37344452 DOI: 10.1038/s41467-023-38962-z] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/24/2022] [Accepted: 05/24/2023] [Indexed: 06/23/2023] Open
Abstract
Subpopulations of soluble, misfolded proteins can bypass chaperones within cells. The extent of this phenomenon and how it happens at the molecular level are unknown. Through a meta-analysis of the experimental literature we find that in all quantitative protein refolding studies there is always a subpopulation of soluble but misfolded protein that does not fold in the presence of one or more chaperones, and can take days or longer to do so. Thus, some misfolded subpopulations commonly bypass chaperones. Using multi-scale simulation models we observe that the misfolded structures that bypass various chaperones can do so because their structures are highly native like, leading to a situation where chaperones do not distinguish between the folded and near-native-misfolded states. More broadly, these results provide a mechanism by which long-time scale changes in protein structure and function can persist in cells because some misfolded states can bypass components of the proteostasis machinery.
Collapse
Affiliation(s)
- Ritaban Halder
- Department of Chemistry, Pennsylvania State University, University Park, PA, 16802, USA
| | - Daniel A Nissley
- Department of Chemistry, Pennsylvania State University, University Park, PA, 16802, USA
- Department of Statistics, University of Oxford, Oxford, OX1 3LB, UK
| | - Ian Sitarik
- Department of Chemistry, Pennsylvania State University, University Park, PA, 16802, USA
| | - Yang Jiang
- Department of Chemistry, Pennsylvania State University, University Park, PA, 16802, USA
| | - Yiyun Rao
- Molecular, Cellular and Integrative Biosciences Program, The Huck Institutes of the Life Sciences, Pennsylvania State University, University Park, PA, 16802, USA
| | - Quyen V Vu
- Institute of Physics, Polish Academy of Sciences; Al. Lotnikow 32/46, 02-668, Warsaw, Poland
| | - Mai Suan Li
- Institute of Physics, Polish Academy of Sciences; Al. Lotnikow 32/46, 02-668, Warsaw, Poland
- Institute for Computational Sciences and Technology; Quang Trung Software City, Tan Chanh Hiep Ward, District 12, Ho Chi Minh City, Vietnam
| | - Justin Pritchard
- Department of Biomedical Engineering, Pennsylvania State University, State College, PA, 16802, USA
- Huck Institute for the Life Sciences, Pennsylvania State University, State College, PA, 16802, USA
| | - Edward P O'Brien
- Department of Chemistry, Pennsylvania State University, University Park, PA, 16802, USA.
- Bioinformatics and Genomics Graduate Program, The Huck Institutes of the Life Sciences, Pennsylvania State University, University Park, PA, 16802, USA.
- Institute for Computational and Data Sciences, Pennsylvania State University, University Park, PA, 16802, USA.
| |
Collapse
|
14
|
Farshad M, DelloStritto MJ, Suma A, Carnevale V. Detecting Liquid-Liquid Phase Separations Using Molecular Dynamics Simulations and Spectral Clustering. J Phys Chem B 2023; 127:3682-3689. [PMID: 37053472 DOI: 10.1021/acs.jpcb.3c00805] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/15/2023]
Abstract
A stringent test of the accuracy of empirical force fields is reproducing the phase diagram of bulk phases and mixtures. Exploring the phase diagram of mixtures requires the detection of phase boundaries and critical points. In contrast to most solid-liquid transitions, in which a global order parameter (average density) can be used to discriminate between two phases, some demixing transitions entail relatively subtle changes in the local environment of each molecule. In such cases, finite sampling errors and finite-size effects make the identification of trends in local order parameters extremely challenging. Here we analyze one such example, namely a methanol/hexane mixture, and compute several local and global structural properties. We simulate the system at various temperatures and study the structural changes associated with demixing. We show that despite a seemingly continuous transformation between mixed and demixed states, the topological properties of the H-bond network change abruptly as the system crosses the demixing line. In particular, by using spectral clustering, we show that the distribution of cluster sizes develops a fat tail (as expected from percolation theory) in the vicinity of the critical point. We illustrate a simple criterion to identify this behavior, which results from the emergence of large system-spanning clusters from a collection of aggregates. We further tested the spectral clustering analysis on a Lennard-Jones system as a standard example of a system with no H-bonds, and also, in this case, we were able to detect the demixing transition.
Collapse
Affiliation(s)
- Mohsen Farshad
- Institute for Computational Molecular Science, Temple University, Philadelphia, Pennsylvania 19122, United States
- Department of Chemistry, Temple University, Philadelphia, Pennsylvania 19122, United States
| | - Mark J DelloStritto
- Institute for Computational Molecular Science, Temple University, Philadelphia, Pennsylvania 19122, United States
- Department of Chemistry, Temple University, Philadelphia, Pennsylvania 19122, United States
| | - Antonio Suma
- Institute for Computational Molecular Science, Temple University, Philadelphia, Pennsylvania 19122, United States
- Dipartimento di Fisica, Università di Bari, 70121 Bari, Italy
| | - Vincenzo Carnevale
- Institute for Computational Molecular Science, Temple University, Philadelphia, Pennsylvania 19122, United States
- Institute for Genomics and Evolutionary Medicine, Temple University, Philadelphia, Pennsylvania 19122, United States
- Department of Biology, Temple University, Philadelphia, Pennsylvania 19122, United States
| |
Collapse
|
15
|
Kondratyuk N, Ryltsev R, Ankudinov V, Chtchelkatchev N. First-principles calculations of the viscosity in multicomponent metallic melts: Al-Cu-Ni as a test case. J Mol Liq 2023. [DOI: 10.1016/j.molliq.2023.121751] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/03/2023]
|
16
|
Bhatta T, Khanal P, Khanal SP, Adhikari NP. Thermodynamics and Transport Properties of Valine and Cysteine Peptides in Water. J Mol Liq 2023. [DOI: 10.1016/j.molliq.2023.121472] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/18/2023]
|
17
|
Wang L, Zhong W, Yu A. Simulation of O2/N2 behaviors on multi-component polymeric membranes in oxy-fuel combustion system. POWDER TECHNOL 2023. [DOI: 10.1016/j.powtec.2023.118288] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/26/2023]
|
18
|
Shoemaker BA, Domingues TS, Haji-Akbari A. Ideal Conductor Model: An Analytical Finite-Size Correction for Nonequilibrium Molecular Dynamics Simulations of Ion Transport through Nanoporous Membranes. J Chem Theory Comput 2022; 18:7142-7154. [PMID: 36327152 DOI: 10.1021/acs.jctc.2c00375] [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/06/2022]
Abstract
Modulating ion transport through nanoporous membranes is critical to many important chemical and biological separation processes. The corresponding transport timescales, however, are often too long to capture accurately using conventional molecular dynamics (MD). Recently, path sampling techniques, such as forward-flux sampling (FFS), have emerged as attractive alternatives for efficiently and accurately estimating arbitrarily long ionic passage times. Here, we use non-equilibrium MD and FFS to explore how the kinetics and mechanisms of pressure-driven chloride transport through a nanoporous graphitic membrane are affected by its lateral dimensions. We not only find ionic passage times and free energy barriers to decrease dramatically upon increasing the membrane surface area but also observe an abrupt and discontinuous change in the locus of the transition state. These strong finite size effects arise due to the cumulative effect of the periodic images of the leading ion entering the pore on the distribution of the induced excess charge at the membrane surface in the feed. By assuming that the feed is an ideal conductor, we analytically derive a finite size correction term that can be computed from the information obtained from a single simulation and successfully use it to obtain corrected free energy profiles with no dependence on the system size. We then estimate ionic passage times in the thermodynamic limit by assuming an Eyring-type dependence of rates on barriers with a size-independent prefactor. This approach constitutes a universal framework for removing finite size artifacts in molecular simulations of ion transport through nanoporous membranes and biological channel proteins.
Collapse
Affiliation(s)
- Brian A Shoemaker
- Department of Chemical and Environmental Engineering, Yale University, New Haven, Connecticut06520, United States
| | - Tiago S Domingues
- Department of Chemical and Environmental Engineering, Yale University, New Haven, Connecticut06520, United States
| | - Amir Haji-Akbari
- Department of Chemical and Environmental Engineering, Yale University, New Haven, Connecticut06520, United States
| |
Collapse
|
19
|
Gupta O, Roy S, Rao L, Mitra S. Graphene Oxide-Carbon Nanotube (GO-CNT) Hybrid Mixed Matrix Membrane for Pervaporative Dehydration of Ethanol. MEMBRANES 2022; 12:membranes12121227. [PMID: 36557134 PMCID: PMC9783890 DOI: 10.3390/membranes12121227] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/20/2022] [Revised: 11/24/2022] [Accepted: 11/29/2022] [Indexed: 05/26/2023]
Abstract
The pervaporation process is an energy-conservative and environmentally sustainable way for dehydration studies. It efficiently separates close boiling point and azeotrope mixtures unlike the distillation process. The separation of ethanol and water is challenging as ethanol and water form an azeotrope at 95.6 wt.% of ethanol. In the last few decades, various polymers have been used as candidates in membrane preparation for pervaporation (PV) application, which are currently used in the preparation of mixed matrix membranes (MMMs) for ethanol recovery and ethanol dehydration but have not been able to achieve an enhanced performance both in terms of flux and selectivity. Composite membranes comprising of poly (vinyl alcohol) (PVA) incorporated with carboxylated carbon nanotubes (CNT-COOH), graphene oxide (GO) and GO-CNT-COOH mixtures were fabricated for the dehydration of ethanol by pervaporation (PV). The membranes were characterized with Scanning Electron Microscopy (SEM), Fourier Transform Infrared Spectroscopy (FTIR), Thermogravimetric Analysis (TGA), Differential Scanning Calorimetry (DSC), Raman spectroscopy, Raman imaging, contact angle measurement, and water sorption to determine the effects of various nanocarbons on the intermolecular interactions, surface hydrophilicity, and degrees of swelling. The effects of feed water concentration and temperature on the dehydration performance were investigated. The incorporation of nanocarbons led to an increase in the permeation flux and separation factor. At a feed water concentration of 10 wt.%, a permeation flux of 0.87 kg/m2.h and a separation factor of 523 were achieved at 23 °C using a PVA-GO-CNT-COOH hybrid membrane.
Collapse
Affiliation(s)
| | | | | | - Somenath Mitra
- Correspondence: ; Tel.: +1-973-596-5611; Fax: +1-973-596-3586
| |
Collapse
|
20
|
Yang W, Xu D, Diao Y, Zhao J, Jing Z, Guo Y. Molecular dynamics simulations on K2SO4 nucleation in supercritical water. J Mol Liq 2022. [DOI: 10.1016/j.molliq.2022.120565] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/08/2022]
|
21
|
Habibi P, Rahbari A, Blazquez S, Vega C, Dey P, Vlugt TJH, Moultos OA. A New Force Field for OH – for Computing Thermodynamic and Transport Properties of H 2 and O 2 in Aqueous NaOH and KOH Solutions. J Phys Chem B 2022; 126:9376-9387. [DOI: 10.1021/acs.jpcb.2c06381] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/05/2022]
Affiliation(s)
- Parsa Habibi
- Engineering Thermodynamics, Process & Energy Department, Faculty of Mechanical, Maritime and Materials Engineering, Delft University of Technology, Leeghwaterstraat 39, 2628 CBDelft, The Netherlands
- Department of Materials Science and Engineering, Faculty of Mechanical, Maritime and Materials Engineering, Delft University of Technology, Mekelweg 2, 2628 CDDelft, The Netherlands
| | - Ahmadreza Rahbari
- Engineering Thermodynamics, Process & Energy Department, Faculty of Mechanical, Maritime and Materials Engineering, Delft University of Technology, Leeghwaterstraat 39, 2628 CBDelft, The Netherlands
| | - Samuel Blazquez
- Depto. Química Física, Fac. Ciencias Químicas, Universidad Complutense de Madrid, 28040Madrid, Spain
| | - Carlos Vega
- Depto. Química Física, Fac. Ciencias Químicas, Universidad Complutense de Madrid, 28040Madrid, Spain
| | - Poulumi Dey
- Department of Materials Science and Engineering, Faculty of Mechanical, Maritime and Materials Engineering, Delft University of Technology, Mekelweg 2, 2628 CDDelft, The Netherlands
| | - Thijs J. H. Vlugt
- Engineering Thermodynamics, Process & Energy Department, Faculty of Mechanical, Maritime and Materials Engineering, Delft University of Technology, Leeghwaterstraat 39, 2628 CBDelft, The Netherlands
| | - Othonas A. Moultos
- Engineering Thermodynamics, Process & Energy Department, Faculty of Mechanical, Maritime and Materials Engineering, Delft University of Technology, Leeghwaterstraat 39, 2628 CBDelft, The Netherlands
| |
Collapse
|
22
|
Simon JM, Krüger P, Schnell SK, Vlugt TJH, Kjelstrup S, Bedeaux D. Kirkwood-Buff integrals: From fluctuations in finite volumes to the thermodynamic limit. J Chem Phys 2022; 157:130901. [PMID: 36209013 DOI: 10.1063/5.0106162] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
The Kirkwood-Buff theory is a cornerstone of the statistical mechanics of liquids and solutions. It relates volume integrals over the radial distribution function, so-called Kirkwood-Buff integrals (KBIs), to particle number fluctuations and thereby to various macroscopic thermodynamic quantities such as the isothermal compressibility and partial molar volumes. Recently, the field has seen a strong revival with breakthroughs in the numerical computation of KBIs and applications to complex systems such as bio-molecules. One of the main emergent results is the possibility to use the finite volume KBIs as a tool to access finite volume thermodynamic quantities. The purpose of this Perspective is to shed new light on the latest developments and discuss future avenues.
Collapse
Affiliation(s)
- J-M Simon
- Laboratoire Interdisciplinaire Carnot de Bourgogne, UMR-6303 CNRS - Université de Bourgogne Franche-Comté, F-21078 Dijon, France
| | - P Krüger
- Graduate School of Science and Engineering, Molecular Chirality Research Center, Chiba University, Chiba 263-8522, Japan
| | - S K Schnell
- Department of Materials Science and Engineering, NTNU - Norwegian University of Science and Technology, N-7491 Trondheim, Norway
| | - T J H Vlugt
- Process & Energy Department, Delft University of Technology, Leeghwaterstraat 39, 2628CB Delft, The Netherlands
| | - S Kjelstrup
- Center of Excellence PoreLab, Department of Chemistry, Faculty of Natural Sciences, NTNU - Norwegian University of Science and Technology, Trondheim, Norway
| | - D Bedeaux
- Center of Excellence PoreLab, Department of Chemistry, Faculty of Natural Sciences, NTNU - Norwegian University of Science and Technology, Trondheim, Norway
| |
Collapse
|
23
|
Iwashita T, Nagao M, Yoshimori A, Terazima M, Akiyama R. Usefulness of higher-order system-size correction for macromolecule diffusion coefficients: A molecular dynamics study. Chem Phys Lett 2022. [DOI: 10.1016/j.cplett.2022.140096] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
|
24
|
Deshchenya V, Kondratyuk N, Lankin A, Norman G. Molecular dynamics study of sucrose aqueous solutions: from solution structure to transport coefficients. J Mol Liq 2022. [DOI: 10.1016/j.molliq.2022.120456] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2022]
|
25
|
Ricci E, Minelli M, De Angelis MG. Modelling Sorption and Transport of Gases in Polymeric Membranes across Different Scales: A Review. MEMBRANES 2022; 12:membranes12090857. [PMID: 36135877 PMCID: PMC9502097 DOI: 10.3390/membranes12090857] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/29/2022] [Revised: 08/24/2022] [Accepted: 08/27/2022] [Indexed: 06/02/2023]
Abstract
Professor Giulio C. Sarti has provided outstanding contributions to the modelling of fluid sorption and transport in polymeric materials, with a special eye on industrial applications such as membrane separation, due to his Chemical Engineering background. He was the co-creator of innovative theories such as the Non-Equilibrium Theory for Glassy Polymers (NET-GP), a flexible tool to estimate the solubility of pure and mixed fluids in a wide range of polymers, and of the Standard Transport Model (STM) for estimating membrane permeability and selectivity. In this review, inspired by his rigorous and original approach to representing membrane fundamentals, we provide an overview of the most significant and up-to-date modeling tools available to estimate the main properties governing polymeric membranes in fluid separation, namely solubility and diffusivity. The paper is not meant to be comprehensive, but it focuses on those contributions that are most relevant or that show the potential to be relevant in the future. We do not restrict our view to the field of macroscopic modelling, which was the main playground of professor Sarti, but also devote our attention to Molecular and Multiscale Hierarchical Modeling. This work proposes a critical evaluation of the different approaches considered, along with their limitations and potentiality.
Collapse
Affiliation(s)
- Eleonora Ricci
- Department of Civil, Chemical, Environmental and Materials Engineering (DICAM), Alma Mater Studiorum—University of Bologna, 40126 Bologna, Italy
| | - Matteo Minelli
- Department of Civil, Chemical, Environmental and Materials Engineering (DICAM), Alma Mater Studiorum—University of Bologna, 40126 Bologna, Italy
| | - Maria Grazia De Angelis
- Institute for Materials and Processes, School of Engineering, University of Edinburgh, Edinburgh EH9 3FB, UK
| |
Collapse
|
26
|
Mistry A, Yu Z, Peters BL, Fang C, Wang R, Curtiss LA, Balsara NP, Cheng L, Srinivasan V. Toward Bottom-Up Understanding of Transport in Concentrated Battery Electrolytes. ACS CENTRAL SCIENCE 2022; 8:880-890. [PMID: 35912355 PMCID: PMC9335914 DOI: 10.1021/acscentsci.2c00348] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/05/2023]
Abstract
Bottom-up understanding of transport describes how molecular changes alter species concentrations and electrolyte voltage drops in operating batteries. Such an understanding is essential to predictively design electrolytes for desired transport behavior. We herein advocate building a structure-property-performance relationship as a systematic approach to accurate bottom-up understanding. To ensure generalization across salt concentrations as well as different electrolyte types and cell configurations, the property-performance relation must be described using Newman's concentrated solution theory. It uses Stefan-Maxwell diffusivity, ij , to describe the role of molecular motions at the continuum scale. The key challenge is to connect ij to the structure. We discuss existing methods for making such a connection, their peculiarities, and future directions to advance our understanding of electrolyte transport.
Collapse
Affiliation(s)
- Aashutosh Mistry
- Chemical
Sciences and Engineering, Argonne National
Laboratory, Lemont, Illinois 60439, United States
- Joint
Center for Energy Storage Research, Argonne
National Laboratory, Lemont, Illinois 60439, United States
| | - Zhou Yu
- Joint
Center for Energy Storage Research, Argonne
National Laboratory, Lemont, Illinois 60439, United States
- Materials
Science Division, Argonne National Laboratory, Lemont, Illinois 60439, United States
| | - Brandon L. Peters
- Joint
Center for Energy Storage Research, Argonne
National Laboratory, Lemont, Illinois 60439, United States
- Materials
Science Division, Argonne National Laboratory, Lemont, Illinois 60439, United States
| | - Chao Fang
- Department
of Chemical and Biomolecular Engineering, University of California Berkeley, Berkeley, California 94720, United States
- Materials
Sciences Division, Lawrence Berkeley National
Laboratory, Berkeley, California 94720, United States
- Joint Center
for Energy Storage Research, Lawrence Berkeley
National Laboratory, Berkeley, California 94720, United States
| | - Rui Wang
- Department
of Chemical and Biomolecular Engineering, University of California Berkeley, Berkeley, California 94720, United States
- Materials
Sciences Division, Lawrence Berkeley National
Laboratory, Berkeley, California 94720, United States
- Joint Center
for Energy Storage Research, Lawrence Berkeley
National Laboratory, Berkeley, California 94720, United States
| | - Larry A. Curtiss
- Joint
Center for Energy Storage Research, Argonne
National Laboratory, Lemont, Illinois 60439, United States
- Materials
Science Division, Argonne National Laboratory, Lemont, Illinois 60439, United States
| | - Nitash P. Balsara
- Department
of Chemical and Biomolecular Engineering, University of California Berkeley, Berkeley, California 94720, United States
- Materials
Sciences Division, Lawrence Berkeley National
Laboratory, Berkeley, California 94720, United States
- Joint Center
for Energy Storage Research, Lawrence Berkeley
National Laboratory, Berkeley, California 94720, United States
| | - Lei Cheng
- Joint
Center for Energy Storage Research, Argonne
National Laboratory, Lemont, Illinois 60439, United States
- Materials
Science Division, Argonne National Laboratory, Lemont, Illinois 60439, United States
| | - Venkat Srinivasan
- Chemical
Sciences and Engineering, Argonne National
Laboratory, Lemont, Illinois 60439, United States
- Joint
Center for Energy Storage Research, Argonne
National Laboratory, Lemont, Illinois 60439, United States
| |
Collapse
|
27
|
Alam TM, Allers JP, Leverant CJ, Harvey JA. Symbolic regression development of empirical equations for diffusion in Lennard-Jones fluids. J Chem Phys 2022; 157:014503. [DOI: 10.1063/5.0093658] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022] Open
Abstract
Symbolic regression (SR) with a multi-gene genetic program has been used to elucidate new empirical equations describing diffusion in Lennard-Jones (LJ) fluids. Examples include equations to predict self-diffusion in pure LJ fluids and equations describing the finite-size correction for self-diffusion in binary LJ fluids. The performance of the SR-obtained equations was compared to that of both the existing empirical equations in the literature and to the results from artificial neural net (ANN) models recently reported. It is found that the SR equations have improved predictive performance in comparison to the existing empirical equations, even though employing a smaller number of adjustable parameters, but show an overall reduced performance in comparison to more extensive ANNs.
Collapse
Affiliation(s)
- Todd M. Alam
- ACC Consulting New Mexico, Cedar Crest, New Mexico 87008, USA
- Department of Organic Materials Science, Sandia National Laboratories, Albuquerque, New Mexico 87185, USA
| | - Joshua P. Allers
- Department of Organic Materials Science, Sandia National Laboratories, Albuquerque, New Mexico 87185, USA
| | - Calen J. Leverant
- Department of WMD Threats and Aerosol Science, Sandia National Laboratories, Albuquerque, New Mexico 87185, USA
- Department of Chemical Engineering, University of Florida, Gainesville, Florida 32611, USA
| | - Jacob A. Harvey
- Geochemistry Department, Sandia National Laboratories, Albuquerque, New Mexico 87185, USA
| |
Collapse
|
28
|
Allers JP, Keth J, Alam TM. Prediction of Self-Diffusion in Binary Fluid Mixtures Using Artificial Neural Networks. J Phys Chem B 2022; 126:4555-4564. [PMID: 35675158 DOI: 10.1021/acs.jpcb.2c01723] [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
Artificial neural networks (ANNs) were developed to accurately predict the self-diffusion constants for individual components in binary fluid mixtures. The ANNs were tested on an experimental database of 4328 self-diffusion constants from 131 mixtures containing 75 unique compounds. The presence of strong hydrogen bonding molecules may lead to clustering or dimerization resulting in non-linear diffusive behavior. To address this, self- and binary association energies were calculated for each molecule and mixture to provide information on intermolecular interaction strength and were used as input features to the ANN. An accurate, generalized ANN model was developed with an overall average absolute deviation of 4.1%. Forward input feature selection reveals the importance of critical properties and self-association energies along with other fluid properties. Additional ANNs were developed with subsets of the full input feature set to further investigate the impact of various properties on model performance. The results from two specific mixtures are discussed in additional detail: one providing an example of strong hydrogen bonding and the other an example of extreme pressure changes, with the ANN models predicting self-diffusion well in both cases.
Collapse
Affiliation(s)
- Joshua P Allers
- Department of Organic Materials Science, Sandia National Laboratories, Albuquerque, New Mexico 87185, United States.,Virtual Technologies and Engineering, Sandia National Laboratories, Albuquerque, New Mexico 87185, United States
| | - Jane Keth
- Department of Organic Materials Science, Sandia National Laboratories, Albuquerque, New Mexico 87185, United States
| | - Todd M Alam
- Department of Organic Materials Science, Sandia National Laboratories, Albuquerque, New Mexico 87185, United States.,ACC Consulting New Mexico, Cedar Crest, New Mexico 87008, United States
| |
Collapse
|
29
|
Dragulet F, Goyal A, Ioannidou K, Pellenq RJM, Del Gado E. Ion Specificity of Confined Ion-Water Structuring and Nanoscale Surface Forces in Clays. J Phys Chem B 2022; 126:4977-4989. [PMID: 35731697 DOI: 10.1021/acs.jpcb.2c01738] [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
Ion specificity and related Hofmeister effects, which are ubiquitous in aqueous systems, can have spectacular consequences in hydrated clays, where ion-specific nanoscale surface forces can determine large-scale cohesive swelling and shrinkage behaviors of soil and sediments. We have used a semiatomistic computational approach and examined sodium, calcium, and aluminum counterions confined with water between charged surfaces representative of clay materials to show that ion-water structuring in nanoscale confinement is at the origin of surface forces between clay particles which are intrinsically ion-specific. When charged surfaces strongly confine ions and water, the amplitude and oscillations of the net pressure naturally emerge from the interplay of electrostatics and steric effects, which cannot be captured by existing theories. Increasing confinement and surface charge densities promote ion-water structures that increasingly deviate from the ions' bulk hydration shells, being strongly anisotropic, persistent, and self-organizing into optimized, nearly solid-like assemblies where hardly any free water is left. Under these conditions, strongly attractive interactions can prevail between charged surfaces because of the dramatically reduced dielectric screening of water and the highly organized water-ion structures. By unravelling the ion-specific nature of these nanoscale interactions, we provide evidence that ion-specific solvation structures determined by confinement are at the origin of ion specificity in clays and potentially a broader range of confined aqueous systems.
Collapse
Affiliation(s)
- Francis Dragulet
- Department of Physics, Institute for Soft Matter Synthesis and Metrology, Georgetown University, 37th and O Streets NW, Washington, D.C. 20057, United States
| | - Abhay Goyal
- Department of Physics, Institute for Soft Matter Synthesis and Metrology, Georgetown University, 37th and O Streets NW, Washington, D.C. 20057, United States.,Infrastructure Materials Group, Engineering Laboratory, National Institute of Standards and Technology, 100 Bureau Drive, Gaithersburg, Maryland 20899, United States
| | - Katerina Ioannidou
- Laboratoire de Mécanique et Génie Civil, CNRS Université de Montpellier, Montpellier 34090, France
| | - Roland J-M Pellenq
- EPiDaPo, The Joint CNRS and George Washington University Laboratory, Children's National Medical Center, Children's Research Institute, 111 Michigan Avenue NW, Washington, D.C. 20010, United States
| | - Emanuela Del Gado
- Department of Physics, Institute for Soft Matter Synthesis and Metrology, Georgetown University, 37th and O Streets NW, Washington, D.C. 20057, United States
| |
Collapse
|
30
|
Jia Q, Venton BJ, DuBay KH. Structure and Dynamics of Adsorbed Dopamine on Solvated Carbon Nanotubes and in a CNT Groove. MOLECULES (BASEL, SWITZERLAND) 2022; 27:molecules27123768. [PMID: 35744896 PMCID: PMC9228466 DOI: 10.3390/molecules27123768] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 05/05/2022] [Revised: 06/01/2022] [Accepted: 06/06/2022] [Indexed: 11/29/2022]
Abstract
Advanced carbon microelectrodes, including many carbon-nanotube (CNT)-based electrodes, are being developed for the in vivo detection of neurotransmitters such as dopamine (DA). Our prior simulations of DA and dopamine-o-quinone (DOQ) on pristine, flat graphene showed rapid surface diffusion for all adsorbed species, but it is not known how CNT surfaces affect dopamine adsorption and surface diffusivity. In this work, we use molecular dynamics simulations to investigate the adsorbed structures and surface diffusion dynamics of DA and DOQ on CNTs of varying curvature and helicity. In addition, we study DA dynamics in a groove between two aligned CNTs to model the spatial constraints at the junctions within CNT assemblies. We find that the adsorbate diffusion on a solvated CNT surface depends upon curvature. However, this effect cannot be attributed to changes in the surface energy roughness because the lateral distributions of the molecular adsorbates are similar across curvatures, diffusivities on zigzag and armchair CNTs are indistinguishable, and the curvature dependence disappears in the absence of solvent. Instead, adsorbate diffusivities correlate with the vertical placement of the adsorbate’s moieties, its tilt angle, its orientation along the CNT axis, and the number of waters in its first hydration shell, all of which will influence its effective hydrodynamic radius. Finally, DA diffuses into and remains in the groove between a pair of aligned and solvated CNTs, enhancing diffusivity along the CNT axis. These first studies of surface diffusion on a CNT electrode surface are important for understanding the changes in diffusion dynamics of dopamine on nanostructured carbon electrode surfaces.
Collapse
|
31
|
Khanal SP, Adhikari NP. Thermodynamic and transport properties of amoxicillin. J Mol Liq 2022. [DOI: 10.1016/j.molliq.2022.118865] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
|
32
|
Grasselli F. Investigating finite-size effects in molecular dynamics simulations of ion diffusion, heat transport, and thermal motion in superionic materials. J Chem Phys 2022; 156:134705. [PMID: 35395883 DOI: 10.1063/5.0087382] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The effects of the finite size of the simulation box in equilibrium molecular dynamics simulations are investigated for prototypical superionic conductors of different types, namely, the fluorite-structure materials PbF2, CaF2, and UO2 (type II), and the α phase of AgI (type I). Largely validated empirical force-fields are employed to run ns-long simulations and extract general trends for several properties, at increasing size and in a wide temperature range. This work shows that, for the considered type-II superionic conductors, the diffusivity dramatically depends on the system size and that the superionic regime is shifted to larger temperatures in smaller cells. Furthermore, only simulations of several hundred atoms are able to capture the experimentally observed, characteristic change in the activation energy of the diffusion process, occurring at the order-disorder transition to the superionic regime. Finite-size effects on ion diffusion are instead much weaker in α-AgI. The thermal conductivity is found generally smaller for smaller cells, where the temperature-independent (Allen-Feldman) regime is also reached at significantly lower temperatures. The finite-size effects on the thermal motion of the non-mobile ions composing the solid matrix follow the simple law that holds for solids.
Collapse
Affiliation(s)
- Federico Grasselli
- COSMO-Laboratory of Computational Science and Modelling, IMX, École Polytechnique Fédérale de Lausanne, 1015 Lausanne, Switzerland
| |
Collapse
|
33
|
Deshchenya VI, Kondratyuk ND, Lankin AV, Norman GE. Modeling of Transport Properties of Aqueous Sucrose Solutions by the Molecular Dynamics Method. RUSSIAN JOURNAL OF PHYSICAL CHEMISTRY A 2022. [DOI: 10.1134/s0036024422030086] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/23/2022]
|
34
|
Hussain S, Haji-Akbari A. How to quantify and avoid finite size effects in computational studies of crystal nucleation: The case of homogeneous crystal nucleation. J Chem Phys 2022; 156:054503. [DOI: 10.1063/5.0079702] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023] Open
Affiliation(s)
- Sarwar Hussain
- Department of Chemical and Environmental Engineering, Yale University, New Haven, Connecticut 06520, USA
| | - Amir Haji-Akbari
- Department of Chemical and Environmental Engineering, Yale University, New Haven, Connecticut 06520, USA
| |
Collapse
|
35
|
Deng X, Zhang Q, Zhang Z, Li Q, Liu X. Adsorption and diffusion behavior of CO2/H2 mixture in calcite slit pores: A molecular simulation study. J Mol Liq 2022. [DOI: 10.1016/j.molliq.2021.118306] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/26/2022]
|
36
|
Omrani S, Ghasemi M, Mahmoodpour S, Shafiei A, Rostami B. Insights from molecular dynamics on CO2 diffusion coefficient in saline water over a wide range of temperatures, pressures, and salinity: CO2 geological storage implications. J Mol Liq 2022. [DOI: 10.1016/j.molliq.2021.117868] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
|
37
|
Desbiens N, Arnault P, Weens W, Dubois V, Perrin G. Bootstrapping time correlation functions of molecular dynamics. Phys Rev E 2021; 104:055310. [PMID: 34942746 DOI: 10.1103/physreve.104.055310] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/20/2021] [Accepted: 11/04/2021] [Indexed: 11/07/2022]
Abstract
Molecular dynamics is often considered as a numerical experiment. The error bars on the results are therefore mandatory, but sometimes difficult to determine and computationally demanding. As a low-cost approach, we describe the application of the bootstrap (BS) method to the quantification of uncertainties pertaining to the time correlation functions. We chose the autocorrelation functions of velocity and interdiffusion current for a binary ionic mixture as a test bed, and we assessed the merit of the Darken approximation relating both of them. The intrinsic errors related to phase space sampling is investigated comparing the BS method with the reference method of replica. We also study how the BS method can assist in addressing the finite-size effects.
Collapse
Affiliation(s)
| | | | - William Weens
- CEA, DAM, DIF, 91297 Arpajon, France.,Laboratoire en Informatique Haute Performance pour le Calcul et la Simulation, 91680 Bruyères-le-Châtel, France
| | | | - Guillaume Perrin
- COSYS, Université Gustave Eiffel, 77420 Champs-sur-Marne, France
| |
Collapse
|
38
|
Shi W, Tiwari SP, Thompson RL, Culp JT, Hong L, Hopkinson DP, Smith K, Resnik K, Steckel JA, Siefert NS. Computational Screening of Physical Solvents for CO 2 Pre-combustion Capture. J Phys Chem B 2021; 125:13467-13481. [PMID: 34734716 DOI: 10.1021/acs.jpcb.1c07268] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
A computational scheme was used to screen physical solvents for CO2 pre-combustion capture by integrating the commercial NIST database, an in-house computational database, chem-informatics, and molecular modeling. A commercially available screened hydrophobic solvent, diethyl sebacate, was identified from the screening with favorable physical properties and promising absorption performance. The promising performance to use diethyl sebacate in CO2 pre-combustion capture has also been confirmed from experiments. Water loading in diethyl sebacate is very low, and therefore, water is kept with H2 in the gas stream. The favorable CO2 interaction with diethyl sebacate and the intermediate solvent free volume fraction leads to both high CO2 solubility and high CO2/H2 solubility selectivity in diethyl sebacate. An in-house NETL computational database was built to characterize CO2, H2, N2, and H2O interactions with 202 different chemical functional groups. It was found that 13% of the functional groups belong to the strong interaction category with the CO2 interaction energy between -15 and -21 kJ/mol; 62% of the functional groups interact intermediately with CO2 (-8 to -15 kJ/mol). The remaining 25% of functional groups interact weakly with CO2 (below -8 kJ/mol). In addition, calculations show that CO2 interactions with the functional groups are stronger than N2 and H2 interactions but are weaker than H2O interactions. The CO2 and H2O interactions with the same functional groups exhibit a very strong linear positive correlation coefficient of 0.92. The relationship between CO2 and H2 gas solubilities and solvent fractional free volume (FFV) has been systematically studied for seven solvents at 298.2 K. A skewed bell-shaped relation was obtained between CO2 solubility and solvent FFV. When an organic compound has a density approximately 10% lower than its density at 298.2 K and 1 bar, it exhibits the highest CO2 loading at that specific solvent density and FFV. Note that the solvent densities were varied using simulations, which are difficult to be obtained from the experiment. In contrast, H2 solubility results exhibit an almost perfect positive linear correlation with the solvent FFV. The theoretical maximum and minimum physical CO2 solubilities in any organic compound at 298.2 K were estimated to be 11 and 0.4 mol/MPa L, respectively. An examination of 182 experimental CO2 physical solubility data and 29 simulated CO2 physical solubilities shows that all the CO2 physical solubility data are within the maximum and minimum with only a few exceptions. Finally, simulations suggest that in order to develop physical solvents with both high CO2 solubility and high CO2/H2 solubility selectivity, the solvents should contain functional groups which are available to interact strongly with CO2 while minimizing FFV.
Collapse
Affiliation(s)
- Wei Shi
- National Energy Technology Laboratory, 626 Cochrans Mill Road, P.O. Box 10940, Pittsburgh, Pennsylvania 15236-0940, United States.,NETL Support Contractor, 626 Cochrans Mill Road, P.O. Box 10940, Pittsburgh, Pennsylvania 15236-0940, United States
| | - Surya P Tiwari
- National Energy Technology Laboratory, 626 Cochrans Mill Road, P.O. Box 10940, Pittsburgh, Pennsylvania 15236-0940, United States.,NETL Support Contractor, 626 Cochrans Mill Road, P.O. Box 10940, Pittsburgh, Pennsylvania 15236-0940, United States
| | - Robert L Thompson
- National Energy Technology Laboratory, 626 Cochrans Mill Road, P.O. Box 10940, Pittsburgh, Pennsylvania 15236-0940, United States.,NETL Support Contractor, 626 Cochrans Mill Road, P.O. Box 10940, Pittsburgh, Pennsylvania 15236-0940, United States
| | - Jeffrey T Culp
- National Energy Technology Laboratory, 626 Cochrans Mill Road, P.O. Box 10940, Pittsburgh, Pennsylvania 15236-0940, United States.,NETL Support Contractor, 626 Cochrans Mill Road, P.O. Box 10940, Pittsburgh, Pennsylvania 15236-0940, United States
| | - Lei Hong
- National Energy Technology Laboratory, 626 Cochrans Mill Road, P.O. Box 10940, Pittsburgh, Pennsylvania 15236-0940, United States.,NETL Support Contractor, 626 Cochrans Mill Road, P.O. Box 10940, Pittsburgh, Pennsylvania 15236-0940, United States
| | - David P Hopkinson
- National Energy Technology Laboratory, 626 Cochrans Mill Road, P.O. Box 10940, Pittsburgh, Pennsylvania 15236-0940, United States
| | - Kathryn Smith
- National Energy Technology Laboratory, 626 Cochrans Mill Road, P.O. Box 10940, Pittsburgh, Pennsylvania 15236-0940, United States.,Carbon Capture Scientific, 4000 Brownsville Road, South Park, Pennsylvania 15129, United States
| | - Kevin Resnik
- National Energy Technology Laboratory, 626 Cochrans Mill Road, P.O. Box 10940, Pittsburgh, Pennsylvania 15236-0940, United States.,NETL Support Contractor, 626 Cochrans Mill Road, P.O. Box 10940, Pittsburgh, Pennsylvania 15236-0940, United States
| | - Janice A Steckel
- National Energy Technology Laboratory, 626 Cochrans Mill Road, P.O. Box 10940, Pittsburgh, Pennsylvania 15236-0940, United States
| | - Nicholas S Siefert
- National Energy Technology Laboratory, 626 Cochrans Mill Road, P.O. Box 10940, Pittsburgh, Pennsylvania 15236-0940, United States
| |
Collapse
|
39
|
Chen X, Wang Y, Wu L, Zhang W, Hu Y. Testing and validation of a self-diffusion coefficient model based on molecular dynamics simulations. Chin J Chem Eng 2021. [DOI: 10.1016/j.cjche.2021.04.036] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
|
40
|
Srinivasan H, Sharma VK, Mitra S. Water accelerates the hydrogen-bond dynamics and abates heterogeneity in deep eutectic solvent based on acetamide and lithium perchlorate. J Chem Phys 2021; 155:024505. [PMID: 34266283 DOI: 10.1063/5.0054942] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022] Open
Abstract
Deep eutectic solvents (DESs) have become a prevalent and promising medium in various industrial applications. The addition of water to DESs has attracted a lot of attention as a scheme to modulate their functionalities and improve their physicochemical properties. In this work, we study the effects of water on an acetamide based DES by probing its microscopic structure and dynamics using classical molecular dynamics simulation. It is observed that, at low water content, acetamide still remains the dominant solvate in the first solvation shell of lithium ions, however, beyond 10 wt. %, it is replaced by water. The increase in the water content in the solvent accelerates the H-bond dynamics by drastically decreasing the lifetimes of acetamide-lithium H-bond complexes. Additionally, water-lithium H-bond complexes are also found to form, with systematically longer lifetimes in comparison to acetamide-lithium complexes. Consequently, the diffusivity and ionic conductivity of all the species in the DES are found to increase substantially. Non-Gaussianity parameters for translational motions of acetamide and water in the DES show a conspicuous decrease with addition of water in the system. The signature of jump-like reorientation of acetamide is observed in the DES by quantifying the deviation from rotational Brownian motion. However, a notable decrease in the deviation is observed with an increase in the water content in the DES. This study demonstrates the intricate connection between H-bond dynamics and various microscopic dynamical parameters in the DES, by investigating the modulation of the former with addition of water.
Collapse
Affiliation(s)
- H Srinivasan
- Solid State Physics Division, Bhabha Atomic Research Centre, Mumbai 400085, India
| | - V K Sharma
- Solid State Physics Division, Bhabha Atomic Research Centre, Mumbai 400085, India
| | - S Mitra
- Solid State Physics Division, Bhabha Atomic Research Centre, Mumbai 400085, India
| |
Collapse
|
41
|
Celebi AT, Dawass N, Moultos OA, Vlugt TJH. How sensitive are physical properties of choline chloride-urea mixtures to composition changes: Molecular dynamics simulations and Kirkwood-Buff theory. J Chem Phys 2021; 154:184502. [PMID: 34241035 DOI: 10.1063/5.0049064] [Citation(s) in RCA: 16] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
Deep eutectic solvents (DESs) have emerged as a cheaper and greener alternative to conventional organic solvents. Choline chloride (ChCl) mixed with urea at a molar ratio of 1:2 is one of the most common DESs for a wide range of applications such as electrochemistry, material science, and biochemistry. In this study, molecular dynamics simulations are performed to study the effect of urea content on the thermodynamic and transport properties of ChCl and urea mixtures. With increased mole fraction of urea, the number of hydrogen bonds (HBs) between cation-anion and ion-urea decreases, while the number of HBs between urea-urea increases. Radial distribution functions (RDFs) for ChCl-urea and ChCl-ChCl pairs shows a significant decrease as the mole fraction of urea increases. Using the computed RDFs, Kirkwood-Buff Integrals (KBIs) are computed. KBIs show that interactions of urea-urea become stronger, while interactions of urea-ChCl and ChCl-ChCl pairs become slightly weaker with increasing mole fraction of urea. All thermodynamic factors are found larger than one, indicating a non-ideal mixture. Our results also show that self- and collective diffusivities increase, while viscosities decrease with increasing urea content. This is mainly due to the weaker interactions between ions and urea, resulting in enhanced mobilities. Ionic conductivities exhibit a non-monotonic behavior. Up to a mole fraction of 0.5, the ionic conductivities increase with increasing urea content and then reach a plateau.
Collapse
Affiliation(s)
- Alper T Celebi
- Engineering Thermodynamics, Process and Energy Department, Faculty of Mechanical, Maritime, and Materials Engineering, Delft University of Technology, Leeghwaterstraat 39, 2628CB Delft, The Netherlands
| | - Noura Dawass
- Engineering Thermodynamics, Process and 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 and 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 and Energy Department, Faculty of Mechanical, Maritime, and Materials Engineering, Delft University of Technology, Leeghwaterstraat 39, 2628CB Delft, The Netherlands
| |
Collapse
|
42
|
Gu HY, Gao W, Gong XG. Hyperdynamics simulations with ab initio forces. J Chem Phys 2021; 154:214112. [PMID: 34240996 DOI: 10.1063/5.0047669] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
By applying the locally optimal rotation method to deal with the lowest eigenvalue of a Hessian matrix, we have efficiently incorporated the hyperdynamics method into the ab initio scheme. In the present method, we only need to calculate the first derivative of the potential and several more force calls in each molecular dynamics (MD) step, which makes hyperdynamics simulation applicable in ab initio MD simulations. With this implementation, we are able to simulate defect diffusion in silicon with boost factors up to 105. We utilized both direct MD and the hyperdynamics method to investigate diffusion of lithium atoms and silicon vacancies in silicon. We identified the complex diffusion process. The obtained diffusion coefficients of Li atoms and Si vacancies are in good agreement with the direct MD results.
Collapse
Affiliation(s)
- Hong-Yang Gu
- Key Laboratory for Computational Physical Sciences (MOE), State Key Laboratory of Surface Physics, Department of Physics, Fudan University, Shanghai 200433, China
| | - Weiguo Gao
- School of Mathematical Sciences, Fudan University, Shanghai 200433, China
| | - Xin-Gao Gong
- Key Laboratory for Computational Physical Sciences (MOE), State Key Laboratory of Surface Physics, Department of Physics, Fudan University, Shanghai 200433, China
| |
Collapse
|
43
|
Klein T, Piszko M, Kankanamge CJ, Kasapis G, Fröba AP. Fick Diffusion Coefficient in Binary Mixtures of [HMIM][NTf 2] and Carbon Dioxide by Dynamic Light Scattering and Molecular Dynamics Simulations. J Phys Chem B 2021; 125:5100-5113. [PMID: 33960781 DOI: 10.1021/acs.jpcb.1c01616] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Dynamic light scattering (DLS) experiments and equilibrium molecular dynamics (EMD) simulations were performed in the saturated liquid phase of the binary mixture of 1-hexyl-3-methylimidazolium bis(trifluormethylsulfonyl)imide ([HMIM][NTf2]) and carbon dioxide (CO2) to access the Fick diffusion coefficient (D11). The investigations were performed within or close to saturation conditions at temperatures between (298.15 and 348.15) K and CO2 mole fractions (xCO2) up to 0.81. The DLS experiments were combined with polarization-difference Raman spectroscopy (PDRS) to simultaneously access the composition of the liquid phase. For the first time in an electrolyte-based system, D11 was directly calculated from EMD simulations by accessing the Maxwell-Stefan (MS) diffusion coefficient and the thermodynamic factor. Agreement within combined uncertainties was found between D11 from DLS and EMD simulations for CO2 mole fractions up to 0.5. In general, an increasing D11 with increasing xCO2 could be observed, with a local maximum present at a CO2 mole fraction of about 0.75. The local maximum could be explained by an increasing MS diffusion coefficient with increasing xCO2 over the entire studied composition range and a decreasing thermodynamic factor at xCO2 above 0.7. Finally, PDRS and EMD simulations were combined to investigate the influence of the fluid structure on the diffusive process.
Collapse
Affiliation(s)
- Tobias Klein
- Institute of Advanced Optical Technologies - Thermophysical Properties (AOT-TP), Department of Chemical and Biological Engineering (CBI) and Erlangen Graduate School in Advanced Optical Technologies (SAOT), Friedrich-Alexander-University Erlangen-Nürnberg (FAU), Paul-Gordan-Straße 8, 91052 Erlangen, Germany
| | - Maximilian Piszko
- Institute of Advanced Optical Technologies - Thermophysical Properties (AOT-TP), Department of Chemical and Biological Engineering (CBI) and Erlangen Graduate School in Advanced Optical Technologies (SAOT), Friedrich-Alexander-University Erlangen-Nürnberg (FAU), Paul-Gordan-Straße 8, 91052 Erlangen, Germany
| | - Chathura J Kankanamge
- Institute of Advanced Optical Technologies - Thermophysical Properties (AOT-TP), Department of Chemical and Biological Engineering (CBI) and Erlangen Graduate School in Advanced Optical Technologies (SAOT), Friedrich-Alexander-University Erlangen-Nürnberg (FAU), Paul-Gordan-Straße 8, 91052 Erlangen, Germany
| | - Georgios Kasapis
- Institute of Advanced Optical Technologies - Thermophysical Properties (AOT-TP), Department of Chemical and Biological Engineering (CBI) and Erlangen Graduate School in Advanced Optical Technologies (SAOT), Friedrich-Alexander-University Erlangen-Nürnberg (FAU), Paul-Gordan-Straße 8, 91052 Erlangen, Germany.,Institute for Multiscale Thermofluids, School of Engineering, University of Edinburgh, Robert Stevenson Road, Edinburgh EH9 3FB, U.K
| | - Andreas P Fröba
- Institute of Advanced Optical Technologies - Thermophysical Properties (AOT-TP), Department of Chemical and Biological Engineering (CBI) and Erlangen Graduate School in Advanced Optical Technologies (SAOT), Friedrich-Alexander-University Erlangen-Nürnberg (FAU), Paul-Gordan-Straße 8, 91052 Erlangen, Germany
| |
Collapse
|
44
|
Fick diffusion coefficients via molecular dynamics: An alternative approach in the Fourier domain. J Mol Liq 2021. [DOI: 10.1016/j.molliq.2021.115460] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/17/2022]
|
45
|
Jahanbakhsh Bonab P, Rastkar Ebrahimzadeh A, Jahanbin Sardroodi J. Insights into the interactions and dynamics of a DES formed by phenyl propionic acid and choline chloride. Sci Rep 2021; 11:6384. [PMID: 33737540 PMCID: PMC7973556 DOI: 10.1038/s41598-021-85260-z] [Citation(s) in RCA: 22] [Impact Index Per Article: 7.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/07/2020] [Accepted: 02/28/2021] [Indexed: 01/31/2023] Open
Abstract
Deep eutectic solvents (DESs) have received much attention in modern green chemistry as inexpensive and easy to handle analogous ionic liquids. This work employed molecular dynamics techniques to investigate the structure and dynamics of a DES system composed of choline chloride and phenyl propionic acid as a hydrogen bond donor and acceptor, respectively. Dynamical parameters such as mean square displacement, liquid phase self-diffusion coefficient and viscosity are calculated at the pressure of 0.1 MPa and temperatures 293, 321 and 400 K. The system size effect on the self-diffusion coefficient of DES species was also examined. Structural parameters such as liquid phase densities, hydrogen bonds, molecular dipole moment of species, and radial and spatial distribution functions (RDF and SDF) were investigated. The viscosity of the studied system was compared with the experimental values recently reported in the literature. A good agreement was observed between simulated and experimental values. The electrostatic and van der Waals nonbonding interaction energies between species were also evaluated and interpreted in terms of temperature. These investigations could play a vital role in the future development of these designer solvents.
Collapse
Affiliation(s)
- Parisa Jahanbakhsh Bonab
- grid.411468.e0000 0004 0417 5692Department of Chemistry, Faculty of Basic Sciences, Azarbaijan Shahid Madani University, Tabriz, Iran ,grid.411468.e0000 0004 0417 5692Molecular Simulation Laboratory (MSL), Azarbaijan Shahid Madani University, Tabriz, Iran
| | - Alireza Rastkar Ebrahimzadeh
- grid.411468.e0000 0004 0417 5692Molecular Simulation Laboratory (MSL), Azarbaijan Shahid Madani University, Tabriz, Iran ,grid.411468.e0000 0004 0417 5692Computational Nanomaterials Research Group (CNRG), Department of Physics, Faculty of Basic Sciences, Azarbaijan Shahid Madani University, Tabriz, Iran
| | - Jaber Jahanbin Sardroodi
- grid.411468.e0000 0004 0417 5692Department of Chemistry, Faculty of Basic Sciences, Azarbaijan Shahid Madani University, Tabriz, Iran ,grid.411468.e0000 0004 0417 5692Molecular Simulation Laboratory (MSL), Azarbaijan Shahid Madani University, Tabriz, Iran
| |
Collapse
|
46
|
Fong KD, Self J, McCloskey BD, Persson KA. Ion Correlations and Their Impact on Transport in Polymer-Based Electrolytes. Macromolecules 2021. [DOI: 10.1021/acs.macromol.0c02545] [Citation(s) in RCA: 20] [Impact Index Per Article: 6.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/08/2023]
Affiliation(s)
- Kara D. Fong
- Department of Chemical and Biomolecular Engineering, University of California, Berkeley, Berkeley, California 94720, United States
- Energy Technologies Area, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| | - Julian Self
- Energy Technologies Area, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
- Department of Materials Science and Engineering, University of California, Berkeley, Berkeley, California 94720, United States
| | - Bryan D. McCloskey
- Department of Chemical and Biomolecular Engineering, University of California, Berkeley, Berkeley, California 94720, United States
- Energy Technologies Area, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| | - Kristin A. Persson
- Department of Materials Science and Engineering, University of California, Berkeley, Berkeley, California 94720, United States
- The Molecular Foundry, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| |
Collapse
|
47
|
Chatwell RS, Guevara-Carrion G, Gaponenko Y, Shevtsova V, Vrabec J. Diffusion of the carbon dioxide-ethanol mixture in the extended critical region. Phys Chem Chem Phys 2021; 23:3106-3115. [PMID: 33491706 DOI: 10.1039/d0cp04985a] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The effect of traces of ethanol in supercritical carbon dioxide on the mixture's thermodynamic properties is studied by molecular simulations and Taylor dispersion measurements. This mixture is investigated along the isobar p = 10 MPa in the temperature range between T = 304 and 343 K. Along this path, the mixture undergoes two transitions: First, the Widom line is crossed, marking the transition from liquid-like to gas-like conditions. A second transition occurs from the supercritical gas-like domain to a subcritical gas. The Widom line crossover entails inflection points for most of the studied properties, i.e. density, enthalpy, shear viscosity, Maxwell-Stefan and intradiffusion coefficients. On the other hand, the transition between the super- and subcritical regions is found to be generally smooth, an observation that is qualitatively confirmed by experimental Taylor dispersion measurements. Dedicated atomistic simulations show the presence of microheterogeneities due to ethanol self-association along the investigated path, which lead to the mixture's anomalous behavior in its extended critical region.
Collapse
Affiliation(s)
- René Spencer Chatwell
- Thermodynamics and Process Engineering, Technische Universität Berlin, 10587 Berlin, Germany.
| | | | - Yuri Gaponenko
- Microgravity Research Center, Université Libre de Bruxelles, 1050 Bruxelles, Belgium
| | - Valentina Shevtsova
- Microgravity Research Center, Université Libre de Bruxelles, 1050 Bruxelles, Belgium
| | - Jadran Vrabec
- Thermodynamics and Process Engineering, Technische Universität Berlin, 10587 Berlin, Germany.
| |
Collapse
|
48
|
Zhao X, Luo T, Jin H. A predictive model for self-, Maxwell-Stefan, and Fick diffusion coefficients of binary supercritical water mixtures. J Mol Liq 2021. [DOI: 10.1016/j.molliq.2020.114735] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2022]
|
49
|
Chaimovich M, Chaimovich A. Relative Resolution: A Computationally Efficient Implementation in LAMMPS. J Chem Theory Comput 2021; 17:1045-1059. [PMID: 33512166 DOI: 10.1021/acs.jctc.0c01003] [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/13/2022]
Abstract
Recently, a novel type of multiscale simulation, called Relative Resolution (RelRes), was introduced. In a single system, molecules switch their resolution in terms of their relative separation, with near neighbors interacting via fine-grained potentials yet far neighbors interacting via coarse-grained potentials; notably, these two potentials are analytically parametrized by a multipole approximation. This multiscale approach is consequently able to correctly retrieve across state space the structural and thermal, as well as static and dynamic, behavior of various nonpolar mixtures. Our current work focuses on the practical implementation of RelRes in LAMMPS, specifically for the commonly used Lennard-Jones potential. By examining various correlations and properties of several alkane liquids, including complex solutions of alternate cooligomers and block copolymers, we confirm the validity of this automated LAMMPS algorithm. Most importantly, we demonstrate that this RelRes implementation gains almost an order of magnitude in computational efficiency, as compared with conventional simulations. We thus recommend this novel LAMMPS algorithm for anyone studying systems governed by Lennard-Jones interactions.
Collapse
Affiliation(s)
- Mark Chaimovich
- Russian School of Mathematics, North Bethesda, Maryland 20852, United States
| | - Aviel Chaimovich
- Department of Chemical and Biological Engineering, Drexel University, Philadelphia, Pennsylvania 19104, United States
| |
Collapse
|
50
|
Celebi AT, Vlugt TJH, Moultos OA. Thermal conductivity of aqueous solutions of reline, ethaline, and glyceline deep eutectic solvents; a molecular dynamics simulation study. Mol Phys 2021. [DOI: 10.1080/00268976.2021.1876263] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022]
Affiliation(s)
- Alper T. Celebi
- Engineering Thermodynamics, Process & Energy Department, Faculty of Mechanical, Maritime and Materials Engineering, Delft University of Technology, Delft, The Netherlands
| | - Thijs J. H. Vlugt
- Engineering Thermodynamics, Process & Energy Department, Faculty of Mechanical, Maritime and Materials Engineering, Delft University of Technology, Delft, The Netherlands
| | - Othonas A. Moultos
- Engineering Thermodynamics, Process & Energy Department, Faculty of Mechanical, Maritime and Materials Engineering, Delft University of Technology, Delft, The Netherlands
| |
Collapse
|