A computer simulation study of fluids in model slit, tubular, and cubic micropores

Nanoconfined Fluids: Uniqueness of Water Compared to Other Liquids

Nanoconfinement can drastically change the behavior of liquids, puzzling us with counterintuitive properties. It is relevant in applications, including decontamination and crystallization control. However, it still lacks a systematic analysis for fluids with different bulk properties. Here we address this gap. We compare, by molecular dynamics simulations, three different liquids in a graphene slit pore: (1) A simple fluid, such as argon, described by a Lennard-Jones potential; (2) an anomalous fluid, such as a liquid metal, modeled with an isotropic core-softened potential; and (3) water, the prototypical anomalous liquid, with directional HBs. We study how the slit-pore width affects the structure, thermodynamics, and dynamics of the fluids. All the fluids show similar oscillating properties by changing the pore size. However, their free-energy minima are quite different in nature: (i) are energy-driven for the simple liquid; (ii) are entropy-driven for the isotropic core-softened potential; and (iii) have a changing nature for water. Indeed, for water, the monolayer minimum is entropy driven, at variance with the simple liquid, while the bilayer minimum is energy driven, at variance with the other anomalous liquid. Also, water has a large increase in diffusion for subnm slit pores, becoming faster than bulk. Instead, the other two fluids have diffusion oscillations much smaller than water, slowing down for decreasing slit-pore width. Our results, clarifying that water confined at the subnm scale behaves differently from other (simple or anomalous) fluids under similar confinement, are possibly relevant in nanopores applications, for example, in water purification from contaminants.

Osmosis in semi-permeable pores: an examination of the basic flow equations based on an experimental and molecular dynamics study

Classically ‘semi-permeable’ pores are generally considered to mediate osmotic flow at a rate dependent upon the hydraulic conductance of the pore and the difference in water potential. The shape or size of the solute molecules is not considered to exert a first-order effect on the flow rate nor is the hydraulic conductance thought to be solute dependent. By the experimental measurement of osmosis in the biological pore AQP (aquaporin) and hard-sphere molecular dynamics simulation of a model pore, we show here that the solute radius can have a profound effect on the osmotic flow rate, causing it to decline steeply with decreasing solute radius.

Using a simple non-equilibrium thermodynamic theory, we propose that an additional ‘osmotic flow coefficient’ is required to describe flows in semi-permeable structures such as AQPs, and that the fall in flow rate with radius represents a conversion from hydraulic to diffusive water flow due to increasing penetration of the pore by the solute. The interaction between the pore geometry and the solute size cannot, therefore, be overlooked, although for every solute the system obeys the criterion for semi-permeability required by basic thermodynamics. The osmotic pore theory therefore reveals a novel and potentially rich structure that remains to be explored in full.

Osmosis in small pores: a molecular dynamics study of the mechanism of solvent transport

Osmosis through semi–permeable pores is a complex process by which solvent is driven by its free energy gradient towards a solute–rich reservoir. We have studied osmotic flow across a semi–permeable cylindrical pore using hard–sphere molecular dynamics which simulates osmosis in the absence of attractive forces between solute and solvent. In addition, we recorded the rates of pressure–driven solvent flow and the diffusive flow of labelled solvent under concentration gradients. It is apparent that there are differences, which are radius dependent, between viscous and diffusive solvent permeabilities in small pores.

The osmotic flow rate is decreased by allowing solute entry into part of the pore, an effect which is not due to solute obstruction. The flow rate is dependent on the structure of the pore, which for asymmetric pores leads, surprisingly, to flow asymmetry or osmotic rectification. In the absence of any possible viscous rectification at these very low flow rates the effect correlates with changes between diffusive and pressure flows created by the presence of solute, an effect which has been predicted from thermodynamic arguments. The geometry of a semi–permeable pore in relation to the solute size is therefore required to predict the osmotic flow rate, a departure from the classical picture.

Finally, by extracting transport parameters from simulations with pure solvent, we examine the departure of observed flow rate from that predicted by continuum mechanics, obtaining drag coefficients which we compare with those derived from hydrodynamics alone.