Capillary condensation in cylindrical pores: Monte Carlo study of the interplay of surface and finite size effects

Polymer solution structure and dynamics within pores of hexagonally close-packed nanoparticles

Using coarse-grained molecular dynamics simulations, we examine structure and dynamics of polymer solutions under confinement within the pores of a hexagonally close-packed (HCP) nanoparticle system with nanoparticle diameter fifty times that of the polymer Kuhn segment size. We model a condition where the polymer chain is in a good solvent (i.e., polymer-polymer interaction is purely repulsive and polymer-solvent and solvent-solvent interactions are attractive) and the polymer-nanoparticle and solvent-nanoparticle interactions are purely repulsive. We probe three polymer lengths (N = 10, 114, and 228 Kuhn segments) and three solution concentrations (1, 10, and 25%v) to understand how the polymer chain conformations and chain center-of-mass diffusion change under confinement within the pores of the HCP nanoparticle structure from those seen in bulk. The known trend of bulk polymer R_g² decreasing with increasing concentration no longer holds when confined in the pores of HCP nanoparticle structure; for example, for the 114-mer, the HCP 〈R_g²〉 at 1%v concentration is lower than HCP 〈R_g²〉 at 10%v concentration. The 〈R_g²〉 of the 114-mer and 228-mer exhibit the largest percent decline going from bulk to HCP at the 1%v concentration and the smallest percent decline at the 25%v concentration. We also provide insight into how the confinement ratio (CR) of polymer chain size to pore size within tetrahedral and octahedral pores in the HCP arrangement of nanoparticles affects the chain conformation and diffusion at various concentrations. At the same concentration, the N = 114 has significantly more movement between pores than the N = 228 chains. For the N = 114 polymer, the diffusion between pores (i.e., inter-pore diffusion) accelerates the overall diffusion rate for the confined HCP system while for the N = 228 polymer, the polymer diffusion in the entire HCP is dominated by the diffusion within the tetrahedral or octahedral pores with minor contributions from inter-pore diffusion. These findings augment the fundamental understanding of macromolecular diffusion through large, densely packed nanoparticle assemblies and are relevant to research focused on fabrication of polymer composite materials for chemical separations, storage, optics, and photonics. We perform coarse-grained molecular dynamics simulations to understand structure and dynamics of polymer solutions under confinement within hexagonal close packed nanoparticles with radii much larger than the polymer chain's bulk radius of gyration.

Machine-Learned Free Energy Surfaces for Capillary Condensation and Evaporation in Mesopores

Using molecular simulations, we study the processes of capillary condensation and capillary evaporation in model mesopores. To determine the phase transition pathway, as well as the corresponding free energy profile, we carry out enhanced sampling molecular simulations using entropy as a reaction coordinate to map the onset of order during the condensation process and of disorder during the evaporation process. The structural analysis shows the role played by intermediate states, characterized by the onset of capillary liquid bridges and bubbles. We also analyze the dependence of the free energy barrier on the pore width. Furthermore, we propose a method to build a machine learning model for the prediction of the free energy surfaces underlying capillary phase transition processes in mesopores.

Phase behavior of colloid-polymer mixtures in planar, spherical, and cylindrical confinement: A density functional theory study

The Asakura-Oosawa (AO) model of colloid-polymer mixtures has been extensively studied over the past several decades both via computer simulations and Density Functional Theory (DFT). At this point, its structural and thermodynamic properties both in the bulk and in contact with flat structureless walls are well understood. At the same time, the phase behavior of AO mixtures in spherical cavities and cylindrical pores, while thoroughly investigated by simulations, has not received a comparably detailed DFT treatment. In this paper, we use the DFT results for the AO model in the bulk and under planar confinement as a point of reference for studying its thermodynamic and structural properties in cavities and pores. The accuracy of the DFT approach is assessed by comparing its predictions with the available extensive simulation data; good overall agreement is generally found with some notable exceptions in the vicinity of wetting and drying transitions. The deviations of the phase behavior in confinement from the bulk phase diagram are analyzed using the Kelvin equation, which is seen to work reasonably well under moderate confinement, i.e., for sufficiently large radii of confining cavities and pores.

Cylindrical confinement of solutions containing semiflexible macromolecules: surface-induced nematic order versus phase separation

Solutions of semiflexible polymers confined in cylindrical pores with repulsive walls are studied by Molecular Dynamics simulations for a wide range of polymer concentrations. Both the case where both lengths are of the same order and the case when the persistence length by far exceeds the contour length are considered, and the enhancement of nematic order along the cylinder axis is characterized. With increasing density the character of the surface effect changes from depletion to the formation of a layered structure. For binary 50 : 50 mixtures of the two types of polymers an interplay between surface enrichment of the stiffer component and the isotropic-nematic transition is found, and a phase separated structure with cylindrical symmetry occurs, with the isotropic phase located around the cylinder axis. For melt densities the mixed nematic phase forms at the wall a layer with a screw-like structure of a tilted smectic phase. The observed behavior is tentatively interpreted in terms of the competition of the chain orientational entropy with entropy of mixing and excluded volume due to the wall.

Nucleation of Capillary Bridges and Bubbles in Nanoconfined CO2

Using molecular simulation, we examine the capillary condensation and the capillary evaporation of CO₂ in cylindrical nanopores. More specifically, we employ the recently developed μV T-S method to determine the microscopic mechanism associated with these processes and the corresponding free energy profiles. We calculate the free energy barrier for capillary condensation and identify that the key step consists in the nucleation of a liquid bridge of a critical size. Similarly, the free energy maximum found for the capillary evaporation process is found to correspond to the nucleation of a vapor bubble of a critical size. In addition, we assess the impact of the strength of the wall-fluid on the height of the free energy barrier and on the critical size of liquid bridges (condensation process) and vapor bubbles (evaporation process). We observe that the height of the free energy barrier increases with the strength of the wall-fluid interactions. Finally, we build a theoretical model, based on capillary theory, to rationalize our findings. In particular, the simulation results reveal a linear scaling of the free energy barrier with the critical size, in excellent agreement with the theoretical predictions.

Hysteresis of liquid adsorption in porous media by coarse-grained Monte Carlo with direct experimental validation

The effects of path-dependent wetting and drying manifest themselves in many types of physical systems, including nanomaterials, biological systems, and porous media such as soil. It is desirable to better understand how these hysteretic macroscopic properties result from a complex interplay between gasses, liquids, and solids at the pore scale. Coarse-Grained Monte Carlo (CGMC) is an appealing approach to model these phenomena in complex pore spaces, including ones determined experimentally. We present two-dimensional CGMC simulations of wetting and drying in two systems with pore spaces determined by sections from micro X-ray computed tomography: a system of randomly distributed spheres and a system of Ottawa sand. Results for the phase distribution, water uptake, and matric suction when corrected for extending to three dimensions show excellent agreement with experimental measurements on the same systems. This supports the hypothesis that CGMC can generate metastable configurations representative of experimental hysteresis and can also be used to predict hysteretic constitutive properties of particular experimental systems, given pore space images.

Soft matter in hard confinement: phase transition thermodynamics, structure, texture, diffusion and flow in nanoporous media

Spatial confinement in nanoporous media affects the structure, thermodynamics and mobility of molecular soft matter often markedly. This article reviews thermodynamic equilibrium phenomena, such as physisorption, capillary condensation, crystallisation, self-diffusion, and structural phase transitions as well as selected aspects of the emerging field of spatially confined, non-equilibrium physics, i.e. the rheology of liquids, capillarity-driven flow phenomena, and imbibition front broadening in nanoporous materials. The observations in the nanoscale systems are related to the corresponding bulk phenomenologies. The complexity of the confined molecular species is varied from simple building blocks, like noble gas atoms, normal alkanes and alcohols to liquid crystals, polymers, ionic liquids, proteins and water. Mostly, experiments with mesoporous solids of alumina, gold, carbon, silica, and silicon with pore diameters ranging from a few up to 50 nm are presented. The observed peculiarities of nanopore-confined condensed matter are also discussed with regard to applications. A particular emphasis is put on texture formation upon crystallisation in nanoporous media, a topic both of high fundamental interest and of increasing nanotechnological importance, e.g. for the synthesis of organic/inorganic hybrid materials by melt infiltration, the usage of nanoporous solids in crystal nucleation or in template-assisted electrochemical deposition of nano structures.

Molecular dynamics simulation of heterogeneous nucleation on nanotubes

Condensation rate inside and outside the tube depends on the pore diameter and length.

Logarithmic finite-size effects on interfacial free energies: phenomenological theory and Monte Carlo studies

The computation of interfacial free energies between coexisting phases (e.g., saturated vapor and liquid) by computer simulation methods is still a challenging problem due to the difficulty of an atomistic identification of an interface and interfacial fluctuations on all length scales. The approach to estimate the interfacial tension from the free-energy excess of a system with interfaces relative to corresponding single-phase systems does not suffer from the first problem but still suffers from the latter. Considering d-dimensional systems with interfacial area L(d-1) and linear dimension L(z) in the direction perpendicular to the interface, it is argued that the interfacial fluctuations cause logarithmic finite-size effects of order ln(L)/L(d-1) and order ln(L(z))/L(d-1), in addition to regular corrections (with leading-order const/L(d-1)). A phenomenological theory predicts that the prefactors of the logarithmic terms are universal (but depend on the applied boundary conditions and the considered statistical ensemble). The physical origin of these corrections are the translational entropy of the interface as a whole, "domain breathing" (coupling of interfacial fluctuations to the bulk order parameter fluctuations of the coexisting domains), and capillary waves. Using a new variant of the ensemble switch method, interfacial tensions are found from Monte Carlo simulations of d = 2 and d = 3 Ising models and a Lennard-Jones fluid. The simulation results are fully consistent with the theoretical predictions.

Critical behaviour of the Ising ferromagnet confined in quasi-cylindrical pores: a Monte Carlo study

The critical behaviour of the Ising ferromagnet confined in pores of radius R and length L is studied by means of Monte Carlo computer simulations. Quasi-cylindrical pores are obtained by replicating n-times a triangular lattice disc of radius R, where L = na and a is the spacing between consecutive replications. So, spins placed at the surface of the pores have less nearest-neighbours (NN) as compared to 8 NN for spins in the bulk. These "missing neighbour" effects undergone by surface spins cause a strong suppression of surface ordering, leading to an ordinary surface transition. Also, the effect propagates into the bulk for small tubes (R ≤ 12) and the effective critical temperature of the pores is shifted towards lower values than in the bulk case. By applying the standard finite-size scaling theory, subsequently supported by numerical data, we concluded that data collapse of relevant observables, e.g., magnetization (m), susceptibility, specific heat, etc., can only be observed by comparing simulation results obtained by keeping the aspect ratio C ≡ R∕L constant. Also, by extrapolating "effective" R-dependent critical temperatures to the thermodynamic limit (R → ∞, C fixed), we obtained T(C)(∞) = 6.208(4). As suggested by finite-size scaling arguments, the magnetization is measured at the critical point scales according to [|m|]Tc R(β/ν) is proportional to [R/L](1/2), where β and ν are the standard exponents for the order parameter and the correlation length, respectively. Furthermore, it is shown that close to criticality the axial correlation length decreases exponentially with the distance. That result is the signature of the formation of (randomly distributed) alternating domains of different magnetization, which can be directly observed by means of snapshot configurations, whose typical length (ξ) is given by the characteristic length of the exponential decay of correlations. Moreover, we show that at criticality ξ = 0.43(2)R.

Mixing-demixing phase diagram for simple liquids in nonuniform electric fields

We deduce the mixing-demixing phase diagram for binary liquid mixtures in an electric field for various electrode geometries and arbitrary constitutive relation for the dielectric constant. By focusing on the behavior of the liquid-liquid interface, we produce simple analytic expressions for the dependence of the interface location on experimental parameters. We also show that the phase diagram contains regions where liquid separation cannot occur under any applied field. The analytic expression for the boundary "electrostatic binodal" line reveals that the regions' size and shape depend strongly on the dielectric relation between the liquids. Moreover, we predict the existence of an "electrostatic spinodal" line that identifies conditions where the liquids are in a metastable state. We finally construct the phase diagram for closed systems by mapping solutions onto those of an open system via an effective liquid composition. For closed systems at a fixed temperature and mixture composition, liquid separation occurs in a finite "window" of surface potential (or charge density). Higher potentials or charge densities counterintuitively destroy the interface, leading to liquid mixing. These results give valuable guides for experiments by providing easily testable predictions for how liquids behave in nonuniform electric fields.

Effect of capsid confinement on the chromatin organization of the SV40 minichromosome

Using small-angle X-ray scattering, we determined the three-dimensional packing architecture of the minichromosome confined within the SV40 virus. In solution, the minichromosome, composed of closed circular dsDNA complexed in nucleosomes, was shown to be structurally similar to cellular chromatin. In contrast, we find a unique organization of the nanometrically encapsidated chromatin, whereby minichromosomal density is somewhat higher at the center of the capsid and decreases towards the walls. This organization is in excellent agreement with a coarse-grained computer model, accounting for tethered nucleosomal interactions under viral capsid confinement. With analogy to confined liquid crystals, but contrary to the solenoid structure of cellular chromatin, our simulations indicate that the nucleosomes within the capsid lack orientational order. Nucleosomes in the layer adjacent to the capsid wall, however, align with the boundary, thereby inducing a 'molten droplet' state of the chromatin. These findings indicate that nucleosomal interactions suffice to predict the genome organization in polyomavirus capsids and underscore the adaptable nature of the eukaryotic chromatin architecture to nanoscale confinement.

Controlling the wetting properties of the Asakura-Oosawa model and applications to spherical confinement

We demonstrate for the Asakura-Oosawa model and an extension of this model that uses continuous rather than hard potentials, how wetting properties at walls can be easily controlled. By increasing the interaction range of the repulsive wall potential acting on the colloids (while keeping the polymer-wall interactions constant) polymers begin to substitute colloids at walls and the system can be driven from complete wetting of colloids via partial wetting to complete wetting of polymers. As an application, we discuss the morphology and wetting behavior of colloid-polymer mixtures in spherical confinement. We apply the recently developed 'ensemble switch method' where the Hamiltonian is extended to a combination of a system with walls and of a system without walls to calculate the surface excess free energies of colloid-rich and polymer-rich phases. The contact angle then is inferred from Young's equation.

Wetting transition in the two-dimensional Blume-Capel model: a Monte Carlo study

The wetting transition of the Blume-Capel model is studied by a finite-size scaling analysis of L×M lattices where competing boundary fields ±H_{1} act on the first or last row of the L rows in the strip, respectively. We show that using the appropriate anisotropic version of finite-size scaling, critical wetting in d=2 is equivalent to a "bulk" critical phenomenon with exponents α=-1, β=0, and γ=3. These concepts are also verified for the Ising model. For the Blume-Capel model, it is found that the field strength H_{1c}(T) where critical wetting occurs goes to zero when the bulk second-order transition is approached, while H_{1c}(T) stays nonzero in the region where in the bulk a first-order transition from the ordered phase, with nonzero spontaneous magnetization, to the disordered phase occurs. Interfaces between coexisting phases then show interfacial enrichment of a layer of the disordered phase which exhibits in the second-order case a finite thickness only. A tentative discussion of the scaling behavior of the wetting phase diagram near the tricritical point is also given.

Formation of metastable structures by phase separation triggered by initial composition gradients in thin films

Phase separation kinetics of a binary (A,B) mixture contained in a thin film of thickness D induced by a quench from the one-phase region into the miscibility gap is studied by simulations using a Cahn-Hilliard-Cook model. The initial randomly mixed state (50% A, 50% B) contains a concentration gradient perpendicular to the film, while the surfaces of the film are "neutral" (no preference for either A or B). In thermal equilibrium, a pattern of large A-rich and B-rich domains must result, separated by domain walls oriented perpendicularly to the external surfaces of the thin film. However, it is shown that for many choices of D and the strength of the initial gradient Ψ(g), instead a very long-lived metastable layered structure forms, with two domains separated by a single interface parallel to the external walls. The transient time evolution that leads to this structure is interpreted in terms of a competition between domain growth in the bulk and surface-directed spinodal decomposition caused by the gradient during the initial stages. A surprising and potentially useful finding is that a moderate concentration gradient perpendicular to the film does not favor the layered structure but facilitates the approach toward the true equilibrium with just two domain walls perpendicular to the film. This mechanism may have useful applications in producing layered materials.

Phase separation in fluids exposed to spatially periodic external fields

When a fluid is confined within a spatially periodic external field, the liquid-vapor transition is replaced by a different transition called laser-induced condensation (LIC) [Götze et al., Mol. Phys. 101, 1651 (2003)]. In d=3 dimensions, the periodic field induces an additional phase, characterized by large density modulations along the field direction. At the triple point, all three phases (modulated, vapor, and liquid) coexist. At temperatures slightly above the triple point and for low (high) values of the chemical potential, two-phase coexistence between the modulated phase and the vapor (liquid) is observed; by increasing the temperature further, both coexistence regions terminate in critical points. In this paper, we reconsider LIC using the Ising model to resolve a number of open issues. To be specific, we (1) determine the universality class of the LIC critical points and elucidate the nature of the correlations along the field direction, (2) present a mean-field analysis to show how the LIC phase diagram changes as a function of the field wavelength and amplitude, (3) develop a simulation method by which the extremely low tension of the interface between modulated and vapor or liquid phase can be measured, (4) present a finite-size scaling analysis to accurately extract the LIC triple point from finite-size simulation data, and (5) consider the fate of LIC in d=2 dimensions.

Single-file water in nanopores

Water molecules confined to pores with sub-nanometre diameters form single-file hydrogen-bonded chains. In such nanoscale confinement, water has unusual physical properties that are exploited in biology and hold promise for a wide range of biomimetic and nanotechnological applications. The latter can be realized by carbon and boron nitride nanotubes which confine water in a relatively non-specific way and lend themselves to the study of intrinsic properties of single-file water. As a consequence of strong water-water hydrogen bonds, many characteristics of single-file water are conserved in biological and synthetic pores despite differences in their atomistic structures. Charge transport and orientational order in water chains depend sensitively on and are mainly determined by electrostatic effects. Thus, mimicking functions of biological pores with apolar pores and corresponding external fields gives insight into the structure-function relation of biological pores and allows the development of technical applications beyond the molecular devices found in living systems. In this Perspective, we revisit results for single-file water in apolar pores, and examine the similarities and the differences between these simple systems and water in more complex pores.