1
|
Sevilla M, Baptista LA, Kremer K, Cortes-Huerto R. Density fluctuations, solvation thermodynamics, and coexistence curves in grand canonical molecular dynamics simulations. J Chem Phys 2025; 162:080901. [PMID: 39991994 DOI: 10.1063/5.0243895] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/14/2024] [Accepted: 01/31/2025] [Indexed: 02/25/2025] Open
Abstract
Fluid transport across nanometric channels induced by electric, pressure, and concentration gradients is ubiquitous in biological systems and fosters various applications. In this context, computer simulation setups with well-defined open-boundary equilibrium starting states are essential in understanding and assisting experimental studies. However, open-boundary computational methods are scarce and do not typically satisfy all the equilibrium conditions imposed by reality. Namely, in the absence of external gradients, (1) the system of interest (SoI) must be at thermodynamic and chemical equilibrium with an infinite reservoir of particles; (2) the fluctuations of the SoI in equilibrium should sample the grand canonical ensemble; (3) the local solvation thermodynamics, which is extremely sensitive to finite-size effects due to solvent depletion, should be correctly described. This point is particularly relevant for out-of-equilibrium systems; and (4) finally, the method should be robust enough to deal with phase transitions and coexistence conditions in the SoI. In this study, we demonstrate with prototypical liquid systems embedded into a reservoir of ideal gas particles that the adaptive resolution simulation (AdResS) method, coupled with particle insertion/deletion steps (AdResS+PI), satisfies all these requirements. Therefore, the AdResS+PI setup is suitable for performing grand canonical and stationary non-equilibrium simulations of open systems.
Collapse
Affiliation(s)
- Mauricio Sevilla
- Max Planck Institute for Polymer Research, Ackermannweg 10, 55128 Mainz, Germany
| | - Luis A Baptista
- Max Planck Institute for Polymer Research, Ackermannweg 10, 55128 Mainz, Germany
| | - Kurt Kremer
- Max Planck Institute for Polymer Research, Ackermannweg 10, 55128 Mainz, Germany
| | | |
Collapse
|
2
|
Yang Y, Anwari Che Ruslan MF, Zhu W, Zhao G, Sun S. Interfacial Behaviors of the H2O+CO2+CH4+C10H22 System in Three Phase Equilibrium: A Combined Molecular Dynamics Simulation and Density Gradient Theory Investigation. J Mol Liq 2022. [DOI: 10.1016/j.molliq.2022.121031] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
|
3
|
Xue P, Qu M, Shi J, Jiang Y, He N, Zhao T, Luo S, Zhou S, Zhang JJ, Luo Y, Chu G, Li H, Chen JF, Sun SG, Liao HG. In Situ TEM Observation of Stagnant Liquid Layer Activation in Nanochannel. NANO LETTERS 2022; 22:6958-6963. [PMID: 36037446 DOI: 10.1021/acs.nanolett.2c01762] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/15/2023]
Abstract
The kinetics of mass transfer in a stagnant fluid layer next to an interface govern numerous dynamic reactions in diffusional micro/nanopores, such as catalysis, fuel cells, and chemical separation. However, the effect of the interplay between stagnant liquid and flowing fluid on the micro/nanoscopic mass transfer dynamics remains poorly understood. Here, by using liquid cell transmission electron microscopy (TEM), we directly tracked microfluid unit migration at the nanoscale. By tracking the trajectories, an unexpected mass transfer phenomenon in which fluid units in the stagnant liquid layer migrated two orders faster during gas-liquid interface updating was identified. Molecular dynamics (MD) simulations indicated that the chemical potential difference between nanoscale liquid layers led to convective flow, which greatly enhanced mass transfer on the surface. Our study opens up a pathway toward research on mass transfer in the surface liquid layers at high spatial and temporal resolutions.
Collapse
Affiliation(s)
- Peng Xue
- State Key Laboratory of Physical Chemistry of Solid Surfaces, Department of Chemistry, College of Chemistry and Chemical Engineering, Xiamen University, Xiamen 361005, People's Republic of China
| | - Mi Qu
- State Key Laboratory of Physical Chemistry of Solid Surfaces, Department of Chemistry, College of Chemistry and Chemical Engineering, Xiamen University, Xiamen 361005, People's Republic of China
| | - Jie Shi
- Beijing Advanced Innovation Center for Soft Matter Science and Engineering, Beijing University of Chemical Technology, Beijing 100029, People's Republic of China
- State Key Laboratory of Organic Inorganic Composites, Beijing University of Chemical Technology, Beijing 100029, People's Republic of China
| | - Youhong Jiang
- State Key Laboratory of Physical Chemistry of Solid Surfaces, Department of Chemistry, College of Chemistry and Chemical Engineering, Xiamen University, Xiamen 361005, People's Republic of China
| | - Nana He
- State Key Laboratory of Physical Chemistry of Solid Surfaces, Department of Chemistry, College of Chemistry and Chemical Engineering, Xiamen University, Xiamen 361005, People's Republic of China
| | - Tiqing Zhao
- State Key Laboratory of Physical Chemistry of Solid Surfaces, Department of Chemistry, College of Chemistry and Chemical Engineering, Xiamen University, Xiamen 361005, People's Republic of China
| | - Shiwen Luo
- State Key Laboratory of Physical Chemistry of Solid Surfaces, Department of Chemistry, College of Chemistry and Chemical Engineering, Xiamen University, Xiamen 361005, People's Republic of China
| | - Shiyuan Zhou
- State Key Laboratory of Physical Chemistry of Solid Surfaces, Department of Chemistry, College of Chemistry and Chemical Engineering, Xiamen University, Xiamen 361005, People's Republic of China
| | - Jia-Jun Zhang
- Xiamen Chip-Nova Technology Co., Ltd., Xiamen 361005, People's Republic of China
| | - Yong Luo
- State Key Laboratory of Organic Inorganic Composites, Beijing University of Chemical Technology, Beijing 100029, People's Republic of China
| | - Guangwen Chu
- State Key Laboratory of Organic Inorganic Composites, Beijing University of Chemical Technology, Beijing 100029, People's Republic of China
| | - Hui Li
- Beijing Advanced Innovation Center for Soft Matter Science and Engineering, Beijing University of Chemical Technology, Beijing 100029, People's Republic of China
| | - Jian-Feng Chen
- Beijing Advanced Innovation Center for Soft Matter Science and Engineering, Beijing University of Chemical Technology, Beijing 100029, People's Republic of China
- State Key Laboratory of Organic Inorganic Composites, Beijing University of Chemical Technology, Beijing 100029, People's Republic of China
| | - Shi-Gang Sun
- State Key Laboratory of Physical Chemistry of Solid Surfaces, Department of Chemistry, College of Chemistry and Chemical Engineering, Xiamen University, Xiamen 361005, People's Republic of China
| | - Hong-Gang Liao
- State Key Laboratory of Physical Chemistry of Solid Surfaces, Department of Chemistry, College of Chemistry and Chemical Engineering, Xiamen University, Xiamen 361005, People's Republic of China
| |
Collapse
|
4
|
Study of interfacial properties of water + methane + oil three-phase systems by a simple molecular simulation protocol. J Mol Liq 2022. [DOI: 10.1016/j.molliq.2022.118951] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
5
|
Smith ER. The importance of reference frame for pressure at the liquid–vapour interface. MOLECULAR SIMULATION 2021. [DOI: 10.1080/08927022.2021.1953697] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/20/2022]
Affiliation(s)
- Edward R. Smith
- Department of Mechanical and Aerospace Engineering, Brunel University London, Uxbridge, London, UK
| |
Collapse
|
6
|
Squarcini A, Malgaretti P. Inhomogeneous surface tension of chemically active fluid interfaces. J Chem Phys 2020; 153:234903. [PMID: 33353340 DOI: 10.1063/5.0025989] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We study the dependence of the surface tension of a fluid interface on the density profile of a third suspended phase. By means of an approximated model for the binary mixture and of a perturbative approach, we derive closed-form expressions for the free energy of the system and for the surface tension of the interface. Our results show a remarkable non-monotonous dependence of the surface tension on the spatial separation between the peaks of the density of the suspended phase. Our results also predict the local value of the surface tension in the case in which the density of the suspended phase is not homogeneous along the interface.
Collapse
Affiliation(s)
- Alessio Squarcini
- Max-Planck-Institut für Intelligente Systeme, Heisenbergstr. 3, D-70569 Stuttgart, Germany
| | - Paolo Malgaretti
- Max-Planck-Institut für Intelligente Systeme, Heisenbergstr. 3, D-70569 Stuttgart, Germany
| |
Collapse
|
7
|
Wei J. Simulation of diffusio-phoretic motion of colloidal particle suppressed by bound solutes within adsorption shell. Mol Phys 2020. [DOI: 10.1080/00268976.2020.1718226] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/25/2022]
Affiliation(s)
- Jiachen Wei
- Institute of Mechanics, Chinese Academy of Sciences, Beijing, People's Republic of China
- Department of Chemistry, University of Cambridge, Cambridge, UK
| |
Collapse
|
8
|
Heidari M, Kremer K, Golestanian R, Potestio R, Cortes-Huerto R. Open-boundary Hamiltonian adaptive resolution. From grand canonical to non-equilibrium molecular dynamics simulations. J Chem Phys 2020; 152:194104. [DOI: 10.1063/1.5143268] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022] Open
Affiliation(s)
- Maziar Heidari
- Max Planck Institute for Polymer Research, Ackermannweg 10, 55128 Mainz, Germany
| | - Kurt Kremer
- Max Planck Institute for Polymer Research, Ackermannweg 10, 55128 Mainz, Germany
| | - Ramin Golestanian
- Max Planck Institute for Dynamics and Self-Organization (MPIDS), 37077 Göttingen, Germany
- Rudolf Peierls Centre for Theoretical Physics, University of Oxford, Oxford OX1 3PU, United Kingdom
| | - Raffaello Potestio
- Physics Department, University of Trento, Via Sommarive, 14, I-38123 Trento, Italy
- INFN-TIFPA, Trento Institute for Fundamental Physics and Applications, I-38123 Trento, Italy
| | | |
Collapse
|
9
|
Wei J, Ramírez-Hinestrosa S, Dobnikar J, Frenkel D. Effect of the interaction strength and anisotropy on the diffusio-phoresis of spherical colloids. SOFT MATTER 2020; 16:3621-3627. [PMID: 32101215 DOI: 10.1039/c9sm02053e] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/10/2023]
Abstract
Gradients in temperature, concentration or electrostatic potential cannot exert forces on a bulk fluid; they can, however, exert forces on a fluid in a microscopic boundary layer surrounding a (nano)colloidal solute, resulting in so-called phoretic flow. Here we present a simulation study of phoretic flow around a spherical colloid held fixed in a concentration gradient. We show that the resulting flow velocity depends non-monotonically on the strength of the colloid-fluid interaction. The reason for this non-monotonic dependence is that solute particles are effectively trapped in a shell around the colloid and cannot contribute to diffusio-phoresis. We also observe that the flow depends sensitively on the anisotropy of solute-colloid interaction.
Collapse
Affiliation(s)
- Jiachen Wei
- Institute of Mechanics, Chinese Academy of Sciences, Beijing 100190, China and Department of Chemistry, University of Cambridge, CB21EW Cambridge, UK.
| | | | - Jure Dobnikar
- Department of Chemistry, University of Cambridge, CB21EW Cambridge, UK. and Institute of Physics, Chinese Academy of Sciences, Beijing 100190, China and Songshan Lake Materials Laboratory, Dongguan 523808, China
| | - Daan Frenkel
- Department of Chemistry, University of Cambridge, CB21EW Cambridge, UK.
| |
Collapse
|
10
|
Marbach S, Bocquet L. Osmosis, from molecular insights to large-scale applications. Chem Soc Rev 2019; 48:3102-3144. [PMID: 31114820 DOI: 10.1039/c8cs00420j] [Citation(s) in RCA: 139] [Impact Index Per Article: 23.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Abstract
Osmosis is a universal phenomenon occurring in a broad variety of processes and fields. It is the archetype of entropic forces, both trivial in its fundamental expression - the van 't Hoff perfect gas law - and highly subtle in its physical roots. While osmosis is intimately linked with transport across membranes, it also manifests itself as an interfacial transport phenomenon: the so-called diffusio-osmosis and -phoresis, whose consequences are presently actively explored for example for the manipulation of colloidal suspensions or the development of active colloidal swimmers. Here we give a global and unifying view of the phenomenon of osmosis and its consequences with a multi-disciplinary perspective. Pushing the fundamental understanding of osmosis allows one to propose new perspectives for different fields and we highlight a number of examples along these lines, for example introducing the concepts of osmotic diodes, active separation and far from equilibrium osmosis, raising in turn fundamental questions in the thermodynamics of separation. The applications of osmosis are also obviously considerable and span very diverse fields. Here we discuss a selection of phenomena and applications where osmosis shows great promises: osmotic phenomena in membrane science (with recent developments in separation, desalination, reverse osmosis for water purification thanks in particular to the emergence of new nanomaterials); applications in biology and health (in particular discussing the kidney filtration process); osmosis and energy harvesting (in particular, osmotic power and blue energy as well as capacitive mixing); applications in detergency and cleaning, as well as for oil recovery in porous media.
Collapse
Affiliation(s)
- Sophie Marbach
- Laboratoire de Physique de l'Ecole Normale Supérieure, ENS, Université PSL, CNRS, Sorbonne Université, Université Paris-Diderot, Sorbonne Paris Cité, Paris, France.
| | | |
Collapse
|
11
|
Ganti R, Liu Y, Frenkel D. Hamiltonian Transformation to Compute Thermo-osmotic Forces. PHYSICAL REVIEW LETTERS 2018; 121:068002. [PMID: 30141645 DOI: 10.1103/physrevlett.121.068002] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/04/2017] [Revised: 04/16/2018] [Indexed: 06/08/2023]
Abstract
If a thermal gradient is applied along a fluid-solid interface, the fluid experiences a thermo-osmotic force. In the steady state, this force is balanced by the gradient of the shear stress. Surprisingly, there appears to be no unique microscopic expression that can be used for computing the magnitude of the thermo-osmotic force. Here we report how, by treating the mass M of the fluid particles as a tensor in the Hamiltonian, we can eliminate the balancing shear force in a nonequilibrium simulation and therefore compute the thermo-osmotic force at simple solid-fluid interfaces. We compare the nonequilibrium force measurement with estimates of the thermo-osmotic force based on computing gradients of the stress tensor. We find that the thermo-osmotic force as measured in our simulations cannot be derived from the most common microscopic definitions of the stress tensor.
Collapse
Affiliation(s)
- Raman Ganti
- Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| | - Yawei Liu
- Beijing University of Chemical Technology, Beijing, People's Republic of China
| | - Daan Frenkel
- Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| |
Collapse
|
12
|
Liu Y, Ganti R, Frenkel D. Pressure gradients fail to predict diffusio-osmosis. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2018; 30:205002. [PMID: 29637907 DOI: 10.1088/1361-648x/aabd58] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/08/2023]
Abstract
We present numerical simulations of diffusio-osmotic flow, i.e. the fluid flow generated by a concentration gradient along a solid-fluid interface. In our study, we compare a number of distinct approaches that have been proposed for computing such flows and compare them with a reference calculation based on direct, non-equilibrium molecular dynamics simulations. As alternatives, we consider schemes that compute diffusio-osmotic flow from the gradient of the chemical potentials of the constituent species and from the gradient of the component of the pressure tensor parallel to the interface. We find that the approach based on treating chemical potential gradients as external forces acting on various species agrees with the direct simulations, thereby supporting the approach of Marbach et al (2017 J. Chem. Phys. 146 194701). In contrast, an approach based on computing the gradients of the microscopic pressure tensor does not reproduce the direct non-equilibrium results.
Collapse
Affiliation(s)
- Yawei Liu
- Beijing Advanced Innovation Center for Soft Matter Science and Engineering, Beijing University of Chemical Technology, Beijing 100029, People's Republic of China
| | | | | |
Collapse
|