No Arabic abstract
We study in this paper the possible existence of Roskilde-simple liquids and their isomorphs in a rough-wall nanoconfinement. Isomorphs are curves in the thermodynamic phase diagram along which structure and dynamics are invariant in suitable nondimensionalized units. Two model liquids using molecular dynamics computer simulations are considered: the single-component Lennard-Jones (LJ) liquid and the Kob-Andersen binary LJ mixture, both of which in the bulk phases are known to have isomorphs. Nanoconfinement is implemented by adopting a slit-pore geometry with fcc crystalline walls; this implies inhomogenous density profiles both parallel and perpendicular to the confining walls. Despite this fact and consistent with an earlier study [Ingebrigtsen et. al, Phys. Rev. Lett. 111, 235901 (2013)] we find that these nanoconfined liquids have isomorphs to a good approximation. More specifically, we show good scaling of inhomogenous density profiles, mean-square displacements, and higher-order structures probed using the topological cluster classification algorithm along the isomorphs. From this study, we conjecture that in experiments, Roskilde-simple liquids may exhibit isomorphs if confined in a suitable manner, for example with carbon nanotubes. Our study thus provides an alternative framework for understanding nanoconfined liquids.
This paper continues the investigation of the exponentially repulsive EXP pair-potential system of Paper I with a focus on isomorphs in the low-temperature gas and liquid phases. As expected from the EXP systems strong virial potential-energy correlations, the systems reduced-unit structure and dynamics are isomorph invariant to a good approximation. Three methods for generating isomorphs are compared: the small-step method that is exact in the limit of small density changes and t
Nanoconfinement can drastically change the behavior of liquids, puzzling us with counterintuitive properties. Moreover, it is relevant in applications, including decontamination and crystallization control. It still lacks a systematic analysis for fluids with different bulk properties. Here we fill this gap. We compare, by molecular dynamics simulations, three different liquids in a graphene slit pore: (A) A simple fluid, such as argon, described by a Lennard-Jones potential; (B) An anomalous fluid, such as a liquid metal, modeled with an isotropic core-softened potential; (C) Water, the prototypical anomalous liquid, with directional hydrogen bonds. We study how the slit-pore width affects the structure, thermodynamics, and dynamics of the fluids. We check that all the fluids, as expected, show similar oscillating properties by changing the pore size. However, the nature of the free-energy minima for the three fluids is quite different: i) only for the simple liquid all the minima are energy-driven, while their structural order increases with decreasing slit-pore width; ii) only for the isotropic core-softened potential all the minima are entropy-driven, while the energy in the minima increases with decreasing slit-pore width; iii) only the water has a changing nature of the minima: 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 diffusion has a large increase for sub-nm 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 clarify that nanoconfined water is unique compared to other (simple or anomalous) fluids under similar confinement, and are possibly relevant in nanopores applications, e.g., in water purification from contaminants.
We review the works devoted to third and fifth harmonic susceptibilities in glasses, namely $chi$ (3) 3 and $chi$ (5) 5. We explain why these nonlinear responses are especially well adapted to test whether or not some amorphous correlations develop upon cooling. We show that the experimental frequency and temperature dependences of $chi$ (3) 3 and of $chi$ (5) 5 have anomalous features, since their behavior is qualitatively different to that of an ideal gas, which is the high-temperature limit of a fluid. Most of the works have interpreted this anomalous behavior as reflecting the growth, upon cooling, of amorphously ordered domains, as predicted by the general framework of Bouchaud and Biroli (BB). We explain why most-if not all-of the challenging interpretations can be recast in a way which is consistent with that of Bouchaud and Biroli. Finally, the comparison of the anomalous features of $chi$ (5) 5 and of $chi$ (3) 3 shows that the amorphously ordered domains are compact, i.e., the fractal dimension d f is close to the dimension d of space. This suggests that the glass transition of molecular liquids corresponds to a new universality class of critical phenomena.
We study a strongly interacting dense hard-sphere system confined between two parallel plates by event-driven molecular dynamics simulations to address the fundamental question of the nature of the 3D to 2D crossover. As the fluid becomes more and more confined the dynamics of the transverse and lateral degrees of freedom decouple, which is accompanied by a diverging time scale separating 2D from 3D behavior. Relying on the time-correlation function of the transversal kinetic energy the scaling behavior and its density-dependence is explored. Surprisingly, our simulations reveal that its time-dependence becomes purely exponential such that memory effects can be ignored. We rationalize our findings quantitatively in terms of an analytic theory which becomes exact in the limit of strong confinement.
A first-principle multiscale modeling approach is presented, which is derived from the solution of the Ornstein-Zernike equation for the coarse-grained representation of polymer liquids. The approach is analytical, and for this reason is transferable. It is here applied to determine the structure of several polymeric systems, which have different parameter values, such as molecular length, monomeric structure, local flexibility, and thermodynamic conditions. When the pair distribution function obtained from this procedure is compared with the results from a full atomistic simulation, it shows quantitative agreement. Moreover, the multiscale procedure accurately captures both large and local scale properties while remaining computationally advantageous.