1
|
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
|
2
|
Akbarishandiz S, Khani S, Maia J. Adhesion dynamics of Janus nanocarriers to endothelial cells: A dissipative particle dynamics study. Phys Rev E 2024; 109:064408. [PMID: 39020963 DOI: 10.1103/physreve.109.064408] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/02/2023] [Accepted: 03/28/2024] [Indexed: 07/20/2024]
Abstract
Janus nanocarriers (NCs) provide promising features in interfacial applications such as targeted drug delivery. Herein, we use dissipative particle dynamics simulations to study the adhesion dynamics of NCs with Janus ligand compositions to the endothelial cell as a function of a series of effects, such as the initial orientation, ligand density, shape, and size of Janus NCs. The Janus NCs, with its long axis parallel to the endothelial glycocalyx (EG) layer, has the best penetration depth due to its lower potential energy and the lowest shell entropy loss. Among different shapes of Janus NCs, both the potential energy and the EG entropy loss control the penetration. In fact, at the parallel orientations, Janus shapes with a robust mechanical strength and larger surface area at the EG/water interface can rotate and penetrate more efficiently. An increase in the ligand density of Janus NCs increases entropy losses of both the hydrophilic and the hydrophobic ligands and decreases the potential energy. Thus, for a specific Janus NCs, functionalizing with an appropriate ligand density would help driving forces prevail over barriers of penetration into the EG layer. For a particular ligand density, once the radius of the Janus NCs exceeds the appropriate size, barriers such as hydrophobic ligands and shell entropy losses are also reinforced significantly and surpass driving forces. Our observations reveal that entropy losses for hydrophobic ligands of Janus NCs and for the shell of NCs are decisive for the adhesion and penetration of Janus NCs to endothelial cells.
Collapse
|
3
|
Saric D, Bell IH, Guevara-Carrion G, Vrabec J. Influence of repulsion on entropy scaling and density scaling of monatomic fluids. J Chem Phys 2024; 160:104503. [PMID: 38456532 DOI: 10.1063/5.0196592] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/08/2024] [Accepted: 02/18/2024] [Indexed: 03/09/2024] Open
Abstract
Entropy scaling is applied to the shear viscosity, self-diffusion coefficient, and thermal conductivity of simple monatomic fluids. An extensive molecular dynamics simulation series is performed to obtain these transport properties and the residual entropy of three potential model classes with variable repulsive exponents: n, 6 Mie (n = 9, 12, 15, and 18), Buckingham's exponential-six (α = 12, 14, 18, and 30), and Tang-Toennies (αT = 4.051, 4.275, and 4.600). A wide range of liquid and supercritical gas- and liquid-like states is covered with a total of 1120 state points. Comparisons to equations of state, literature data, and transport property correlations are made. Although the absolute transport property values within a given potential model class may strongly depend on the repulsive exponent, it is found that the repulsive steepness plays a negligible role when entropy scaling is applied. Hence, the plus-scaled transport properties of n, 6 Mie, exponential-six, and Tang-Toennies fluids lie basically on one master curve, which closely corresponds with entropy scaling correlations for the Lennard-Jones fluid. This trend is confirmed by literature data of n, 6 Mie, and exponential-six fluids. Furthermore, entropy scaling holds for state points where the Pearson correlation coefficient R is well below 0.9. The condition R > 0.9 for strongly correlating liquids is thus not necessary for the successful application of entropy scaling, pointing out that isomorph theory may be a part of a more general framework that is behind the success of entropy scaling. Density scaling reveals a strong influence of the repulsive exponent on this particular approach.
Collapse
Affiliation(s)
- Denis Saric
- Thermodynamics, Technical University of Berlin, Ernst-Reuter-Platz 1, 10587 Berlin, Germany
| | - Ian H Bell
- Applied Chemicals and Materials Division, National Institute of Standards and Technology, Boulder, Colorado 80305, USA
| | | | - Jadran Vrabec
- Thermodynamics, Technical University of Berlin, Ernst-Reuter-Platz 1, 10587 Berlin, Germany
| |
Collapse
|
4
|
Akbarishandiz S, Khani S, Maia J. Adhesion dynamics of functionalized nanocarriers to endothelial cells: a dissipative particle dynamics study. SOFT MATTER 2023; 19:9254-9268. [PMID: 38009071 DOI: 10.1039/d3sm00865g] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/28/2023]
Abstract
Targeted drug delivery to endothelial cells utilizing functionalized nanocarriers (NCs) is an essential procedure in therapeutic and diagnosis therapies. Using dissipative particle dynamics simulation, NCs have been designed and combined with an endothelial environment, such as the endothelial glycocalyx (EG) layer, receptors, water, and cell wall. Furthermore, the energy landscapes of the functionalized NC with the endothelial cell have been analyzed as a function of properties such as the shape, size, initial orientation, and ligand density of NCs. Our results show that an appropriate higher ligand density for each particular NC provides more driving forces than barriers for the penetration of the NCs. Herein we report the importance of shell entropy loss for the NC shape effect on the adhesion and penetration into the EG layer. Moreover, the rotation of the disc shape NC as a wheel during the penetration is an extra driving force for its further inclusion. By increasing the NCs' size larger than the appropriate size for each particular ligand density, due to an increase in the NCs' shell entropy loss, the barriers surpass the driving forces for NC penetration. Furthermore, the parallel orientation provides the NCs with the best penetration capabilities. However, the rotation of the disc shape NCs enhances their diffusion in the perpendicular orientation too. Overall, our findings highlight the crucial role of the shell entropy loss in governing the penetration of NCs. Besides, studying NCs with a homogeneous ligand composition enabled us to cross barriers and probe energetics after the complete inclusion of the NCs.
Collapse
Affiliation(s)
- Saeed Akbarishandiz
- Department of Macromolecular Science and Engineering, Case Western Reserve University, USA.
| | - Shaghayegh Khani
- Department of Macromolecular Science and Engineering, Case Western Reserve University, USA.
| | - Joao Maia
- Department of Macromolecular Science and Engineering, Case Western Reserve University, USA.
| |
Collapse
|
5
|
Sorkin B, Diamant H, Ariel G. Universal Relation between Entropy and Kinetics. PHYSICAL REVIEW LETTERS 2023; 131:147101. [PMID: 37862659 DOI: 10.1103/physrevlett.131.147101] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/15/2023] [Accepted: 08/25/2023] [Indexed: 10/22/2023]
Abstract
Relating thermodynamic and kinetic properties is a conceptual challenge with many practical benefits. Here, based on first principles, we derive a rigorous inequality relating the entropy and the dynamic propagator of particle configurations. It is universal and applicable to steady states arbitrarily far from thermodynamic equilibrium. Applying the general relation to diffusive dynamics yields a relation between the entropy and the (normal or anomalous) diffusion coefficient. The relation can be used to obtain useful bounds for the late-time diffusion coefficient from the calculated steady-state entropy or, conversely, to estimate the entropy based on measured diffusion coefficients. We demonstrate the validity and usefulness of the relation through several examples and discuss its broad range of applications, in particular, for systems far from equilibrium.
Collapse
Affiliation(s)
- Benjamin Sorkin
- School of Chemistry and Center for Physics and Chemistry of Living Systems, Tel Aviv University, 69978 Tel Aviv, Israel
| | - Haim Diamant
- School of Chemistry and Center for Physics and Chemistry of Living Systems, Tel Aviv University, 69978 Tel Aviv, Israel
| | - Gil Ariel
- Department of Mathematics, Bar-Ilan University, 52000 Ramat Gan, Israel
| |
Collapse
|
6
|
Jin J, Schweizer KS, Voth GA. Understanding dynamics in coarse-grained models. II. Coarse-grained diffusion modeled using hard sphere theory. J Chem Phys 2023; 158:034104. [PMID: 36681632 DOI: 10.1063/5.0116300] [Citation(s) in RCA: 10] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022] Open
Abstract
The first paper of this series [J. Chem. Phys. 158, 034103 (2023)] demonstrated that excess entropy scaling holds for both fine-grained and corresponding coarse-grained (CG) systems. Despite its universality, a more exact determination of the scaling relationship was not possible due to the semi-empirical nature. In this second paper, an analytical excess entropy scaling relation is derived for bottom-up CG systems. At the single-site CG resolution, effective hard sphere systems are constructed that yield near-identical dynamical properties as the target CG systems by taking advantage of how hard sphere dynamics and excess entropy can be analytically expressed in terms of the liquid packing fraction. Inspired by classical equilibrium perturbation theories and recent advances in constructing hard sphere models for predicting activated dynamics of supercooled liquids, we propose a new approach for understanding the diffusion of molecular liquids in the normal regime using hard sphere reference fluids. The proposed "fluctuation matching" is designed to have the same amplitude of long wavelength density fluctuations (dimensionless compressibility) as the CG system. Utilizing the Enskog theory to derive an expression for hard sphere diffusion coefficients, a bridge between the CG dynamics and excess entropy is then established. The CG diffusion coefficient can be roughly estimated using various equations of the state, and an accurate prediction of accelerated CG dynamics at different temperatures is also possible in advance of running any CG simulation. By introducing another layer of coarsening, these findings provide a more rigorous method to assess excess entropy scaling and understand the accelerated CG dynamics of molecular fluids.
Collapse
Affiliation(s)
- Jaehyeok Jin
- Department of Chemistry, Chicago Center for Theoretical Chemistry, Institute for Biophysical Dynamics, and James Franck Institute, The University of Chicago, Chicago, Illinois 60637, USA
| | - Kenneth S Schweizer
- Department of Material Science, Department of Chemistry, Department of Chemical and Biomolecular Engineering, and Materials Research Laboratory, University of Illinois, Urbana, Illinois 61801, USA
| | - Gregory A Voth
- Department of Chemistry, Chicago Center for Theoretical Chemistry, Institute for Biophysical Dynamics, and James Franck Institute, The University of Chicago, Chicago, Illinois 60637, USA
| |
Collapse
|