Properties of a soft-core model of methanol: an integral equation theory and computer simulation study

Size-Pore-Dependent Methanol Sequestration from Water-Methanol Mixtures by an Embedded Graphene Slit

The separation of liquid mixture components is relevant to many applications-ranging from water purification to biofuel production-and is a growing concern related to the UN Sustainable Development Goals (SDGs), such as "Clean water and Sanitation" and "Affordable and clean energy". One promising technique is using graphene slit-pores as filters, or sponges, because the confinement potentially affects the properties of the mixture components in different ways, favoring their separation. However, no systematic study has shown how the size of a pore changes the thermodynamics of the surrounding mixture. Here, we focus on water-methanol mixtures and explore, using Molecular Dynamics simulations, the effects of a graphene pore, with size ranging from 6.5 to 13 Å, for three compositions: pure water, 90%-10%, and 75%-25% water-methanol. We show that tuning the pore size can change the mixture pressure, density and composition in bulk due to the size-dependent methanol sequestration within the pore. Our results can help in optimizing the graphene pore size for filtering applications.

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.

Evidence of Structural Inhomogeneities in Hard-Soft Dimeric Particles without Attractive Interactions

We perform Monte Carlo simulations of a simple hard-soft dimeric model constituted by two tangent spheres experiencing different interactions. Specifically, two hard spheres belonging to different dimers interact via a bare hard-core repulsion, whereas two soft spheres experience a softly repulsive Hertzian interaction. The cross correlations are soft as well. By exploring a wide range of temperatures and densities we investigate the capability of this model to document the existence of structural inhomogeneities indicating the possible onset of aggregates, even if no attraction is set. The fluid phase behavior is studied by analyzing structural and thermodynamical properties of the observed structures, in particular by computing radial distribution functions, structure factors and cluster size distributions. The numerical results are supported by integral equation theories of molecular liquids which allow for a finer and faster spanning of the temperature-density diagram. Our results may serve as a framework for a more systematic investigation of self-assembled structures of functionalized hard-soft dimers able to aggregate in a variety of structures widely oberved in colloidal dispersion.

Monte Carlo simulation and integral equation study of Hertzian spheres in the low-temperature regime

We investigate the behavior of Hertzian spheres in the fluid phase and in proximity of the freezing threshold by using Monte Carlo (MC) simulations and integral equation theories, based on the Ornstein-Zernike (OZ) approach. The study is motivated by the importance of the Hertzian model in representing a large class of systems interacting via soft interactions, such as star polymers or microgels. Radial distribution functions, structure factors, and excess entropy clearly show the reentrant behavior typical of the Hertzian fluid, well caught by both MC simulations and OZ theory. Then, we make use of some phenomenological one-phase criteria for testing their reliability in detecting the freezing threshold. All criteria provide evidence of the fluid-solid transition with different degrees of accuracy. This suggests the possibility to adopt these empirical rules to provide a quick and reasonable estimate of the freezing transition in model potentials of wide interest for soft matter systems.

Structural properties of the Jagla fluid

The structural properties of the Jagla fluid are studied by Monte Carlo (MC) simulations, numerical solutions of integral equation theories, and the (semi-analytical) rational-function approximation (RFA) method. In the latter case, the results are obtained from the assumption (supported by our MC simulations) that the Jagla potential and a potential with a hard core plus an appropriate piecewise constant function lead to practically the same cavity function. The predictions obtained for the radial distribution function, g(r), from this approach are compared against MC simulations and integral equations for the Jagla model, and also for the limiting cases of the triangle-well potential and the ramp potential, with a general good agreement. The analytical form of the RFA in Laplace space allows us to describe the asymptotic behavior of g(r) in a clean way and compare it with MC simulations for representative states with oscillatory or monotonic decay. The RFA predictions for the Fisher-Widom and Widom lines of the Jagla fluid are obtained.

Temperature of maximum density and excess properties of short-chain alcohol aqueous solutions: A simplified model simulation study

We perform an extensive computational study of binary mixtures of water and short-chain alcohols resorting to two-scale potential models to account for the singularities of hydrogen bonded liquids. Water molecules are represented by a well studied core softened potential which is known to qualitatively account for a large number of water's characteristic anomalies. Along the same lines, alcohol molecules are idealized by dimers in which the hydroxyl groups interact with each other and with water with a core softened potential as well. Interactions involving non-polar groups are all deemed purely repulsive. We find that the qualitative behavior of excess properties (excess volume, enthalpy, and constant pressure heat capacity) agrees with that found experimentally for alcohols such as t-butanol in water. Moreover, we observe that our simple solute under certain conditions acts as a "structure-maker," in the sense that the temperature of maximum density of the bulk water model increases as the solute is added, i.e., the anomalous behavior of the solvent is enhanced by the solute.

Integral equation study of soft-repulsive dimeric fluids

We study fluid structure and water-like anomalies of a system constituted by dimeric particles interacting via a purely repulsive core-softened potential by means of integral equation theories. In our model, dimers interact through a repulsive pair potential of inverse-power form with a softened repulsion strength. By employing the Ornstein-Zernike approach and the reference interaction site model (RISM) theory, we study the behavior of water-like anomalies upon progressively increasing the elongation λ of the dimers from the monomeric case ([Formula: see text]) to the tangent configuration ([Formula: see text]). For each value of the elongation we consider two different values of the interaction potential, corresponding to one and two length scales, with the aim to provide a comprehensive description of the possible fluid scenarios of this model. Our theoretical results are systematically compared with already existing or newly generated Monte Carlo data: we find that theories and simulations agree in providing the picture of a fluid exhibiting density and structural anomalies for low values of λ and for both the two values of the interaction potential. Integral equation theories give accurate predictions for pressure and radial distribution functions, whereas the temperatures where anomalies occur are underestimated. Upon increasing the elongation, the RISM theory still predicts the existence of anomalies; the latter are no longer observed in simulations, since their development is likely precluded by the onset of crystallization. We discuss our results in terms of the reliability of integral equation theories in predicting the existence of water-like anomalies in core-softened fluids.

Lattice model for water-solute mixtures

A lattice model for the study of mixtures of associating liquids is proposed. Solvent and solute are modeled by adapting the associating lattice gas (ALG) model. The nature of interaction of solute/solvent is controlled by tuning the energy interactions between the patches of ALG model. We have studied three set of parameters, resulting in, hydrophilic, inert, and hydrophobic interactions. Extensive Monte Carlo simulations were carried out, and the behavior of pure components and the excess properties of the mixtures have been studied. The pure components, water (solvent) and solute, have quite similar phase diagrams, presenting gas, low density liquid, and high density liquid phases. In the case of solute, the regions of coexistence are substantially reduced when compared with both the water and the standard ALG models. A numerical procedure has been developed in order to attain series of results at constant pressure from simulations of the lattice gas model in the grand canonical ensemble. The excess properties of the mixtures, volume and enthalpy as the function of the solute fraction, have been studied for different interaction parameters of the model. Our model is able to reproduce qualitatively well the excess volume and enthalpy for different aqueous solutions. For the hydrophilic case, we show that the model is able to reproduce the excess volume and enthalpy of mixtures of small alcohols and amines. The inert case reproduces the behavior of large alcohols such as propanol, butanol, and pentanol. For the last case (hydrophobic), the excess properties reproduce the behavior of ionic liquids in aqueous solution.

Effects of confinement between attractive and repulsive walls on the thermodynamics of an anomalous fluid

We study via molecular-dynamics simulations the thermodynamics of an anomalous fluid confined in a slit pore with one wall structured and attractive and another unstructured and repulsive. We find that the phase diagram of the homogeneous part of the confined fluid is shifted to higher temperatures, densities, and pressures with respect to the bulk, but it can be rescaled on the bulk case. We calculate a moderate increase of mobility of the homogeneous confined fluid that we interpret as a consequence of the layering due to confinement and the collective modes due to long-range correlations. We show that, as in bulk, the confined fluid has structural, diffusion, and density anomalies that order in the waterlike hierarchy, and a liquid-liquid critical point (LLCP). The overall anomalous region moves to higher temperatures, densities, and pressure, and the LLCP displaces to higher temperature compared to bulk. Motivated by experiments, we calculate also the phase diagram not just for the homogeneous part of the confined fluid but for the entire fluid in the pore, and we show that it is shifted toward higher pressures but preserves the thermodynamics, including the LLCP. Because our model has waterlike properties, we argue that in experiments with supercooled water confined in slit pores with a width of >3 nm if hydrophilic and of >1.5 nm if hydrophobic, the existence of the LLCP could be easier to test than in bulk, where it is not directly accessible.

Density and structural anomalies in soft-repulsive dimeric fluids

We performed a simulation study of the fluid structure of dimeric particles interacting via a core-softened potential and shed light on their anomalous behaviours upon varying both geometrical and interaction parameters.

Existence of a liquid-liquid phase transition in methanol

A simple model is constructed to study the phase diagram and thermodynamic properties of methanol, which is described as a dimer of an apolar sphere mimicking the methyl group and a sphere with core-softened potential as the hydroxyl group. Performing classical Monte Carlo simulations, we obtained the phase diagram, showing a second critical point between two different liquid phases. Evaluating systems with a different number of particles, we extrapolate to infinite size in accordance with Ising universality class to obtain bulk values for critical temperature, pressure, and density. Strong evidence that the structure of the liquid changes upon transition from high- to low-density phase was provided. From the experimentally determined hydrogen bond strength and length in methanol and water, we propose where the second critical point of methanol should be.