Glassy dynamics in a confined monatomic fluid

Computer simulations and mode-coupling theory of glass-forming confined hard-sphere fluids

We present mode-coupling theory (MCT) results for densely packed hard-sphere fluids confined between two parallel walls and compare them quantitatively to computer simulations. The numerical solution of MCT is calculated using the full system of matrix-valued integro-differential equations. We investigate several dynamical properties of supercooled liquids including scattering functions, frequency-dependent susceptibilities, and mean-square displacements. Close to the glass transition, we find quantitative agreement between the coherent scattering function predicted from theory and that evaluated from simulations, which enables us to make quantitative statements on caging and relaxation dynamics of the confined hard-sphere fluid.

Viscoelastic response of fluid trapped between two dissimilar van der Waals surfaces

Employing grand canonical Monte-Carlo and molecular dynamics simulations, the viscoelastic response of trapped fluid under molecularly thin confinement by walls having different wall-fluid interaction strengths, is investigated. With increase in slit asymmetry, given by the ratio of interaction strengths of the wall having strong wall-fluid interaction to that of the wall with weak wall-fluid interaction, a crossover in effective density of the fluid film, from rarer (R) to denser (D) than the bulk density is observed. Upon increasing asymmetry further, the dense fluid (F) layers undergo bond-orientational (S) ordering. The variation of viscoelastic relaxation time with scaled asymmetry shows a universal behavior, independent of slit width, with two distinct regimes. Below a critical value of asymmetry, the viscoelastic relaxation time is a slowly varying function of asymmetry, comparable with the structural relaxation time. Beyond the critical asymmetry, on the other hand, viscoelastic response time shows a sharp increase upon increasing asymmetry, deviating markedly from the structural relaxation time. Interestingly the critical asymmetry value is found to correlate with R to D crossover. The microscopic origin of the two-regime universal behavior of viscoelastic response time is found to stem from the fact that below critical asymmetry, the overall viscoelastic behaviour of the slit is dominated by that of the fast relaxing layer close to the weakly attracting surface, while above the critical asymmetry, the relaxation behaviour is guided by the dense fluid layer adjacent to the strongly attracting wall. In vicinity of fluid to ordering transition, the loss and storage moduli merge for low frequencies as in gel-like mechanical behaviour. The storage modulus takes over the loss modulus in the phase co-existence region even before the long ranged order sets in. Our findings bear important implications for fluid transport in hetero-structured geometry in nanotechnology.

Tagged-particle dynamics in confined colloidal liquids

We present numerical results for the tagged-particle dynamics by solving the mode-coupling theory in confined geometry for colloidal liquids (cMCT). We show that neither the microscopic dynamics nor the type of intermediate scattering function qualitatively changes the asymptotic dynamics in vicinity of the glass transition. In particular, we find similar characteristics of confinement in the low-frequency susceptibility spectrum which we interpret as footprints of parallel relaxation. We derive predictions for the localization length and the scaling of the diffusion coefficient in the supercooled regime and discover a pronounced nonmonotonic dependence on the confinement length. For dilute liquids in the hydrodynamic limit we calculate an analytical expression for the intermediate scattering functions, which is in perfect agreement with event-driven Brownian dynamics simulations. From this, we derive an expression for persistent anticorrelations in the velocity autocorrelation function (VACF) for confined motion. Using numerical results of the cMCT equations for the VACF we also identify a crossover between different scalings corresponding to a transition from unconfined to confined behavior.

Dynamical properties of densely packed confined hard-sphere fluids

Numerical solutions of the mode-coupling theory (MCT) equations for a hard-sphere fluid confined between two parallel hard walls are elaborated. The governing equations feature multiple parallel relaxation channels which significantly complicate their numerical integration. We investigate the intermediate scattering functions and the susceptibility spectra close to structural arrest and compare to an asymptotic analysis of the MCT equations. We corroborate that the data converge in the β-scaling regime to two asymptotic power laws, viz. the critical decay and the von Schweidler law. The numerical results reveal a nonmonotonic dependence of the power-law exponents on the slab width and a nontrivial kink in the low-frequency susceptibility spectra. We also find qualitative agreement of these theoretical results to event-driven molecular dynamics simulations of polydisperse hard-sphere systems. In particular, the nontrivial dependence of the dynamical properties on the slab width is well reproduced.

Nonergodicity parameters of confined hard-sphere glasses

Within a recently developed mode-coupling theory for fluids confined to a slit we elaborate numerical results for the long-time limits of suitably generalized intermediate scattering functions. The theory requires as input the density profile perpendicular to the plates, which we obtain from density functional theory within the fundamental-measure framework, as well as symmetry-adapted static structure factors, which can be calculated relying on the inhomogeneous Percus-Yevick closure. Our calculations for the nonergodicity parameters for both the collective as well as for the self motion are in qualitative agreement with our extensive event-driven molecular dynamics simulations for the intermediate scattering functions for slightly polydisperse hard-sphere systems at high packing fraction. We show that the variation of the nonergodicity parameters as a function of the wavenumber correlates with the in-plane static structure factors, while subtle effects become apparent in the structure factors and relaxation times of higher mode indices. A criterion to predict the multiple reentrant from the variation of the in-plane static structure is presented.

Role of density modulation in the spatially resolved dynamics of strongly confined liquids

Confinement by walls usually produces a strong modulation in the density of dense liquids near the walls. Using molecular dynamics simulations, we examine the effects of the density modulation on the spatially resolved dynamics of a liquid confined between two parallel walls, using a resolution of a fraction of the interparticle distance in the liquid. The local dynamics is quantified by the relaxation time associated with the temporal autocorrelation function of the local density. We find that this local relaxation time varies in phase with the density modulation. The amplitude of the spatial modulation of the relaxation time can be quite large, depending on the characteristics of the wall and thermodynamic parameters of the liquid. To disentangle the effects of confinement and density modulation on the spatially resolved dynamics, we compare the dynamics of a confined liquid with that of an unconfined one in which a similar density modulation is induced by an external potential. We find several differences indicating that density modulation alone cannot account for all the features seen in the spatially resolved dynamics of confined liquids. We also examine how the dynamics near a wall depends on the separation between the two walls and show that the features seen in our simulations persist in the limit of large wall separation.

Non-monotonic effect of confinement on the glass transition

The relaxation dynamics of glass forming liquids and their structure are influenced in the vicinity of confining walls. This effect has mostly been observed to be a monotonic function of the slit width. Recently, a qualitatively new behaviour has been uncovered by Mittal and coworkers, who reported that the single particle dynamics in a hard-sphere fluid confined in a planar slit varies in a non-monotonic way as the slit width is decreased from five to roughly two particle diametres (Mittal et al 2008 Phys. Rev. Lett. 100 145901). In view of the great potential of this effect for applications in those fields of science and industry, where liquids occur under strong confinement (e.g. nano-technology), the number of researchers studying various aspects and consequences of this non-monotonic behaviour has been rapidly growing. This review aims at providing an overview of the research activity in this newly emerging field. We first briefly discuss how competing mechanisms such as packing effects and short-range attraction may lead to a non-monotonic glass transition scenario in the bulk. We then analyse confinement effects on the dynamics of fluids using a thermodynamic route which relates the single particle dynamics to the excess entropy. Moreover, relating the diffusive dynamics to the Widom's insertion probability, the oscillations of the local dynamics with density at moderate densities are fairly well described. At high densities belonging to the supercooled regime, however, this approach breaks down signaling the onset of strongly collective effects. Indeed, confinement introduces a new length scale which in the limit of high densities and small pore sizes competes with the short-range local order of the fluid. This gives rise to a non-monotonic dependence of the packing structure on confinement, with a corresponding effect on the dynamics of structural relaxation. This non-monotonic effect occurs also in the case of a cone-plate type channel, where the degree of confinement varies with distance from the apex. This is a very promising issue for future research with the possibility of uncovering the existence of alternating glassy and liquid-like domains.

Structure and Dynamics of Octamethylcyclotetrasiloxane Confined between Mica Surfaces

Using a molecular model for octamethylcyclotetrasiloxane (OMCTS), molecular dynamics simulations are carried out to probe the phase state of OMCTS confined between two mica surfaces in equilibrium with a reservoir. Molecular dynamics simulations are carried out for elevations ranging from 5 to 35 K above the melting point for the OMCTS model used in this study. The Helmholtz free energy is computed for a specific confinement using the two-phase thermodynamic (2PT) method. Analysis of the in-plane pair correlation functions did not reveal signatures of freezing even under an extreme confinement of two layers. OMCTS is found to orient with a wide distribution of orientations with respect to the mica surface, with a distinct preference for the surface parallel configuration in the contact layers. The self-intermediate scattering function is found to decay with increasing relaxation times as the surface separation is decreased, and the two-step relaxation in the scattering function, a signature of glassy dynamics, distinctly evolves as the temperature is lowered. However, even at 5 K above the melting point, we did not observe a freezing transition and the self-intermediate scattering functions relax within 200 ps for the seven-layered confined system. The self-diffusivity and relaxation times obtained from the Kohlrausch-Williams-Watts stretched exponential fits to the late α-relaxation exhibit power law scalings with the packing fraction as predicted by mode coupling theory. A distinct discontinuity in the Helmholtz free energy, potential energy, and a sharp change in the local bond order parameter, Q4, was observed at 230 K for a five-layered system upon cooling, indicative of a first-order transition. A freezing point depression of about 30 K was observed for this five-layered confined system, and at the lower temperatures, contact layers were found to be disordered with long-range order present only in the inner layers. These dynamical signatures indicate that confined OMCTS undergoes a slowdown akin to a fluid approaching a glass transition upon increasing confinement, and freezing under confinement would require substantial subcooling below the bulk melting point of OMCTS.

Glassy dynamics in confinement: planar and bulk limits of the mode-coupling theory

We demonstrate how the matrix-valued mode-coupling theory of the glass transition and glassy dynamics in planar confinement converges to the corresponding theory for two-dimensional (2D) planar and the three-dimensional bulk liquid, provided the wall potential satisfies certain conditions. Since the mode-coupling theory relies on the static properties as input, the emergence of a homogeneous limit for the matrix-valued intermediate scattering functions is directly connected to the convergence of the corresponding static quantities to their conventional counterparts. We show that the 2D limit is more subtle than the bulk limit, in particular, the in-planar dynamics decouples from the motion perpendicular to the walls. We investigate the frozen-in parts of the intermediate scattering function in the glass state and find that the limits time t→∞ and effective wall separation L→0 do not commute due to the mutual coupling of the residual transversal and lateral force kernels.

Tagged-particle motion in a dense confined liquid

We investigate the dynamics of a tagged particle embedded in a strongly interacting confined liquid enclosed between two opposing flat walls. Using the Zwanzig-Mori projection operator formalism we obtain an equation of motion for the incoherent scattering function suitably generalized to account for the lack of translational symmetry. We close the equations of motion by a self-consistent mode-coupling ansatz. The interaction of the tracer with the surrounding liquid is encoded in generalized direct correlation functions. We extract the in-plane dynamics and provide a microscopic expression for the diffusion coefficient parallel to the walls. The solute particle may differ in size or interaction from the surrounding host-liquid constituents offering the possibility of a systematic analysis of dynamic effects on the tagged-particle motion in confinement.

Distributions of single-molecule properties as tools for the study of dynamical heterogeneities in nanoconfined water

The explicit trend of the distribution functions of single-molecule rotational relaxation constants and atomic mean-square displacement are used to study the dynamical heterogeneities in nanoconfined water. The trend of the single-molecule properties distributions is related to the dynamic heterogeneities, and to the dynamic crossovers found in water clusters of different shapes and sizes and confined in a variety of zeolites. This was true in all the cases that were considered, in spite of the various shapes and sizes of the clusters. It is confirmed that the high temperature dynamical crossover occurring in the temperature range 200-230 K can be interpreted at a molecular level as the formation of almost translationally rigid clusters, characterized by some rotational freedom, hydrogen bond exchange and translational jumps as cage-to-cage processes. We also suggest a mechanism for the low temperature dynamical crossover (LTDC), falling in the temperature range 150-185 K, through which the adsorbed water clusters are made of nearly rigid sub-clusters, slightly mismatched, and thus permitting a relatively free librational motion at their borders. It appears that the condition required for LTDC to occur is the presence of highly heterogeneous environments for the adsorbed molecules, with some dangling hydrogen bonds or weaker than water-water hydrogen bonds. Under these conditions some dynamics are permitted at very low temperature, although most rotational motion is frozen. Therefore, it is unlikely, though not entirely excluded, that LTDC will be found in supercooled bulk water where no heterogeneous interface is present.