Mode-coupling theory of the glass transition for confined fluids

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.

Structural properties of liquids in extreme confinement

We simulate a hard-sphere liquid in confined geometry where the separation of the two parallel, hard walls is smaller than two particle diameters. By systematically reducing the wall separation we analyze the behavior of structural and thermodynamic properties, such as inhomogeneous density profiles, structure factors, and compressibilities when approaching the two-dimensional limit. In agreement with asymptotic predictions, we find for quasi-two-dimensional fluids that the density profile becomes parabolic and the structure factor converges toward its two-dimensional counterpart. To extract the compressibility in polydisperse samples a perturbative expression is used which qualitatively influences the observed nonmonotonic dependence of the compressibility with wall separation. We also present theoretical calculations based on fundamental-measure theory and integral-equation theory, which are in very good agreement with the simulation results.

Layering and packing in confined colloidal suspensions

Confinement modifies the properties of a fluid. The particle density is no longer uniform but depends on the distance from the walls; parallel to the walls, layers with different particle densities form. This affects the particle packing in the layers. We investigated colloidal fluids with volume fractions between 0.19 and 0.32 confined between rough walls. The particle-particle interactions were dominated by hard-sphere interactions but also contained some electrostatic interactions. The particle locations were determined using confocal microscopy and served to calculate the density profile, radial distribution function, anisotropic and generalized structure factors but also to characterize the arrangement of the wall particles leading to the roughness of the walls. The experiments are complemented by molecular dynamics simulations and fundamental-measure theory. While the particle arrangements are mainly controlled by hard-core interactions, electrostatic interactions become more important as the volume fraction decreases. Furthermore, the structure of the rough walls was varied and found to have a significant effect on the fluid structure. An appropriate representation of the rough walls in the simulations is thus crucial to successfully mimic the experiments.

Isomorphs in nanoconfined liquids

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 good 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., 2013, 111, 235901] we find that these two nanoconfined liquids have isomorphs to a good approximation. More specifically, we show good invariance along the isomorphs of inhomogenous density profiles, mean-square displacements, and higher-order structures probed using the topological cluster classification algorithm. Our study thus provides an alternative framework for understanding nanoconfined liquids.

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.

Multiple character of non-monotonic size-dependence for relaxation dynamics in polymer-particle and binary mixtures

Adding plasticizers is a well-known procedure to reduce the glass transition temperature in polymers. It has been recently shown that this effect shows a non-monotonic dependence on the size of additive molecules (2019 J. Chem. Phys. 150 024903). In this work, we demonstrate that, as the size of the additive molecules is changed at fixed concentration, multiple extrema emerge in the dependence of the system's relaxation time on the size ratio. The effect occurs on all relevant length scales including single monomer dynamics, decay of Rouse modes and relaxation of the chain's end-to-end vector. A qualitatively similar trend is found within mode-coupling theoretical results for a binary hard-sphere mixture. An interpretation of the effect in terms of local packing efficiency and coupling between the dynamics of minority and majority species is provided.

Mode-coupling theory of the glass transition for colloidal liquids in slit geometry

We provide a detailed derivation of the mode-coupling equations for a colloidal liquid confined by two parallel smooth walls. We introduce irreducible memory kernels for the different relaxation channels thereby extending the projection operator technique to colloidal liquids in slit geometry. Investigating both the collective dynamics as well as the tagged-particle motion, we prove that the mode-coupling functional assumes the same form as in the Newtonian case corroborating the universality of the glass-transition singularity with respect to the microscopic dynamics.

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.

Study on the Glass Transition Process of Polymer System Using Differential Scanning Calorimetry and Fourier Transform Infrared Spectroscopy

The change in the infrared spectrum of polymer samples with temperature and their differential scanning calorimetry (DSC) experimental results are analyzed. According to the van't Hoff equation at constant pressure, the changes in the absorbance ratio corresponding to high and low vibrational states are calculated, and the apparent enthalpy differences of the vibration energy states transformation of the characteristic group can be obtained. From the experimental results, we can find that characteristic vibration modes of a chemical group in a polymer are under the influence of the glass transition process of the polymer with a different extent. The characteristic vibration modes of the same chemical group behave differently due to the influence of the polymer system at which the chemical moiety is situated.

Diverging Time Scale in the Dimensional Crossover for Liquids in Strong Confinement

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.

Colloidal diffusion in confined geometries

Colloidal diffusion in confined geometries is analysed at the level of anisotropic pair densities.

Impact of solvent granularity and layering on tracer hydrodynamics in confinement

Classic hydrodynamic arguments establish that when a spherical tracer particle is suspended between parallel walls, tracer-wall coupling mediated by the solvent will cause the tracer to exhibit position-dependent diffusivity. We investigate how the diffusivity profiles of confined tracers are impacted by the diameter size-ratio of the tracer to solvent: starting from the classic limit of infinite size-ratio (i.e., continuum solvent), we consider size-ratios of four or less to examine how hydrodynamic predictions are disrupted for systems where the tracer and solvent are of similar scale. We use computer simulations and techniques based on the Fokker-Planck formalism to calculate the diffusivity profiles of hard-sphere tracer particles in hard-sphere solvents, focusing on the dynamics perpendicular to the walls. Given wall separations of several tracer diameters, we first consider confinement between hard walls, where anisotropic structuring at the solvent lengthscale generates inhomogeneity in the tracer free-energy landscape and undermines hydrodynamic predictions locally. We then introduce confining planes that we term transparent walls, which restrict tracer and solvent center-accessibilities while completely eliminating static anisotropy, and reveal position-dependent signatures in tracer diffusivity solely attributable to confinement. With or without suppressing static heterogeneity, we find that tracer diffusivity increasingly deviates on a local basis from hydrodynamic predictions at smaller size-ratios. However, hydrodynamic theory still approximately captures spatially-averaged dynamics across the pores even for very small tracer-solvent size-ratios over a wide range of solvent densities and wall separations.

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.

Strongly confined fluids: Diverging time scales and slowing down of equilibration

The Newtonian dynamics of strongly confined fluids exhibits a rich behavior. Its confined and unconfined degrees of freedom decouple for confinement length L→0. In that case and for a slit geometry the intermediate scattering functions S_{μν}(q,t) simplify, resulting for (μ,ν)≠(0,0) in a Knudsen-gas-like behavior of the confined degrees of freedom, and otherwise in S_{∥}(q,t), describing the structural relaxation of the unconfined ones. Taking the coupling into account we prove that the energy fluctuations relax exponentially. For smooth potentials the relaxation times diverge as L^{-3} and L^{-4}, respectively, for the confined and unconfined degrees of freedom. The strength of the L^{-3} divergence can be calculated analytically. It depends on the pair potential and the two-dimensional pair distribution function. Experimental setups are suggested to test these predictions.

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.

Slow Water Dynamics near a Glass Transition or a Solid Interface: A Common Rationale

Performing molecular dynamics simulations, we investigate the enormous slowdowns of water dynamics when approaching a glass transition or a solid interface. We show that both effects can be described on common grounds within a theoretical framework, which was recently proposed by Schweizer et al. and considers coupled local hopping and elastic distortion. For confined water, we correctly describe the variation of the α-relaxation time, τα, as a function of both temperature and position with respect to the interface. Exploiting our knowledge of a cooperative length scale ξ(T) from the confinement studies, we quantitatively rationalize the glassy slowdown, τα(T), and the Stokes-Einstein breakdown of bulk water. For both confined and bulk liquid, variations of the α-relaxation time are intimately related to changes of the cage-rattling amplitude.

How Local and Average Particle Diffusivities of Inhomogeneous Fluids Depend on Microscopic Dynamics

Computer simulations and a stochastic Fokker-Planck equation based approach are used to compare the single-particle diffusion coefficients of equilibrium hard-sphere fluids exhibiting identical inhomogeneous static structure and governed by either Brownian (i.e., overdamped Langevin) or Newtonian microscopic dynamics. The physics of inhomogeneity is explored via the imposition of one-dimensional sinusoidal density profiles of different wavelengths and amplitudes. When imposed density variations are small in magnitude for distances on the scale of a particle diameter, bulk-like average correlations between local structure and mobility are observed. In contrast, when density variations are significant on that length scale, qualitatively different structure-mobility correlations emerge that are sensitive to the governing microscopic dynamics. Correspondingly, a previously proposed scaling between long-time diffusivities for bulk isotropic fluids of particles exhibiting Brownian versus Newtonian dynamics [Pond et al. Soft Matter 2011, 7, 9859-9862] cannot be generalized to describe the position-dependent behaviors of strongly inhomogeneous fluids. While average diffusivities in the inhomogeneous and homogeneous directions are coupled, their qualitative dependencies on inhomogeneity wavelength are sensitive to the details of the microscopic dynamics. Nonetheless, average diffusivities of the inhomogeneous fluids can be approximately predicted for either type of dynamics based on knowledge of bulk isotropic fluid behavior and how inhomogeneity modifies the distribution of available volume. Analogous predictions for average diffusivities of experimental, inhomogeneous colloidal dispersions (based on known bulk behavior) suggest that they will exhibit qualitatively different trends than those predicted by models governed by overdamped Langevin dynamics that do not account for hydrodynamic interactions.

Simulated glass-forming polymer melts: dynamic scattering functions, chain length effects, and mode-coupling theory analysis

We present molecular-dynamics simulations for a fully flexible model of polymer melts with different chain length N ranging from short oligomers (N = 4) to values near the entanglement length (N = 64). For these systems we explore the structural relaxation of the supercooled melt near the critical temperature T c of mode-coupling theory (MCT). Coherent and incoherent scattering functions are analyzed in terms of the idealized MCT. For temperatures T > T c we provide evidence for the space-time factorization property of the β relaxation and for the time-temperature superposition principle (TTSP) of the α relaxation, and we also discuss deviations from these predictions for T ≈ T c. For T larger than the smallest temperature where the TTSP holds we perform a quantitative analysis of the dynamics with the asymptotic MCT predictions for the late β regime. Within MCT a key quantity, in addition to T c, is the exponent parameter λ. For the fully flexible polymer models studied we find that λ is independent of N and has a value (λ = 0.735 ) typical of simple glass-forming liquids. On the other hand, the critical temperature increases with chain length toward an asymptotic value T c (∞) . This increase can be described by T c (∞) - T c(N) ∼ 1/N and may be interpreted in terms of the N dependence of the monomer density ρ, if we assume that the MCT glass transition is ruled by a soft-sphere-like constant coupling parameter Γ c = ρ c T c (-1/4), where ρ c is the monomer density at T c. In addition, we also estimate T c from a Hansen-Verlet-like criterion and MCT calculations based on structural input from the simulation. For our polymer model both the Hansen-Verlet criterion and the MCT calculations suggest T c to decrease with increasing chain length, in contrast to the direct analysis of the simulation data.

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.

Excess free energy of supercooled liquids at disordered walls

Using a thermodynamic integration scheme, we compute the excess free energy, γ, of a glass-forming, binary Lennard-Jones liquid in contact with a frozen amorphous wall, formed by particles frozen into a similar structure as the liquid. We find that γ is nonzero, becoming negative at low temperature. This indicates that the thermodynamics of the system is perturbed by the effect of the amorphous wall.

Percolation approach to glassy dynamics with continuously broken ergodicity

We show that the relaxation dynamics near a glass transition with continuous ergodicity breaking can be endowed with a geometric interpretation based on percolation theory. At the mean-field level this approach is consistent with the mode-coupling theory (MCT) of type-A liquid-glass transitions and allows one to disentangle the universal and nonuniversal contributions to MCT relaxation exponents. Scaling predictions for the time correlation function are successfully tested in the F(12) schematic model and facilitated spin systems on a Bethe lattice. Our approach immediately suggests the extension of MCT scaling laws to finite spatial dimensions and yields predictions for dynamic relaxation exponents below an upper critical dimension of 6.

The impact range for smooth wall-liquid interactions in nanoconfined liquids

Bulk and nanoconfined liquids have very different physics; for instance, nanoconfined liquids show stratification and position-dependent relaxation processes. A number of similarities between bulk and nanoconfined liquids have nevertheless been reported in computer simulations during the last decade. Inspired by these observations, we present results from molecular dynamics computer simulations of four nanoconfined liquids (the single-component Lennard-Jones liquid, the Kob-Andersen binary Lennard-Jones mixture, an asymmetric dumbbell model, and the Dzugutov liquid) demonstrating also a microscopic similarity between bulk and nanoconfined liquids. The results show that the interaction range for the wall-liquid and liquid-liquid interactions of the nanoconfined liquid is identical to that of the bulk liquid if the liquid is "Roskilde simple" in bulk as well as nanoconfinement, i.e., exhibits strong correlations between virial and potential-energy equilibrium fluctuations in the NVT ensemble. Under this condition, interactions beyond the first coordination shell may be ignored, in particular for the wall-liquid interactions. This is shown not to be the case for non-Roskilde-simple liquids.

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.

Predicting how nanoconfinement changes the relaxation time of a supercooled liquid

The properties of nanoconfined fluids can be strikingly different from those of bulk liquids. A basic unanswered question is whether the equilibrium and dynamic consequences of confinement are related to each other in a simple way. We study this question by simulation of a liquid comprising asymmetric dumbbell-shaped molecules, which can be deeply supercooled without crystallizing. We find that the dimensionless structural relaxation times-spanning six decades as a function of temperature, density, and degree of confinement-collapse when plotted versus excess entropy. The data also collapse when plotted versus excess isochoric heat capacity, a behavior consistent with the existence of isomorphs in the bulk and confined states.

Connection between thermodynamics and dynamics of simple fluids in highly attractive pores

Using molecular simulations, we investigate the structural and diffusive dynamics properties of a model fluid in highly absorptive cylindrical pores. At subcritical temperatures, self-diffusivity displays three distinct regimes as a function of average pore density ρ: (1) a decrease in self-diffusivity with increasing ρ at low ρ, (2) constant self-diffusivity with respect to varying ρ at moderate density, and (3) a decrease in self-diffusivity with increasing ρ at high density. These regimes are closely linked to the thermodynamic properties of the fluid in the pore, specifically, the adsorption isotherm, isosteric heat of adsorption, and the density profile. We show that these three diffusivity regimes qualitatively correspond to three distinct adsorption regimes: monolayer formation, multilayer adsorption, and pore filling, respectively. In addition, we find that the self-diffusivity is a universal function of the local film density in the monolayer formation regime at subcritical temperatures. The results of this work suggest a potential means to estimate the self-diffusivity over a broad pressure range using a limited number of experiments.