Liquid-glass transition of a fluid confined in a disordered porous matrix: a mode-coupling theory

Quantitative Understanding of Liquid Dynamics at Interfaces from a Free-Energy Landscape Perspective

On the basis of molecular dynamics simulations of water and ethanol in nanopores, we devise a methodology to determine the free-energy landscape (FEL) imposed by an interface on an adjoining liquid directly from the particle trajectories. The methodology merely uses the statistical mechanical relation between occupancy and energy and, hence, is particularly suitable in complex situations, e.g., for disordered or rough atomistic interfaces and molecular liquids, as encountered in many biological, geological, and technological situations. Moreover, we show that the thus-obtained FEL enables a quantitative understanding of interface effects on liquid dynamics. Specifically, by determining the local minima and barriers of the FEL and using an Arrhenius-like relation, we reproduce the very strong spatial variation of the structural relaxation time of water and ethanol across nanopores over a broad temperature range. We anticipate that the proposed FEL approach is transferable to various other liquids and interfaces.

Deformable hard particles confined in a disordered porous matrix

With suitably designed Monte Carlo simulations, we have investigated the properties of mobile, impenetrable, yet deformable particles that are immersed into a porous matrix, the latter one realized by a frozen configuration of spherical particles. By virtue of a model put forward by Batista and Miller [Phys. Rev. Lett. 105, 088305 (2010)], the fluid particles can change in their surroundings, formed by other fluid particles or the matrix particles, their shape within the class of ellipsoids of revolution; such a change in shape is related to a change in energy, which is fed into suitably defined selection rules in the deformation "moves" of the Monte Carlo simulations. This concept represents a simple yet powerful model of realistic, deformable molecules with complex internal structures (such as dendrimers or polymers). For the evaluation of the properties of the system, we have used the well-known quenched-annealed protocol (with its characteristic double average prescription) and have analyzed the simulation data in terms of static properties (the radial distribution function and aspect ratio distribution of the ellipsoids) and dynamic features (notably the mean squared displacement). Our data provide evidence that the degree of deformability of the fluid particles has a distinct impact on the aforementioned properties of the system.

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.

Active topological glass

The glass transition in soft matter systems is generally triggered by an increase in packing fraction or a decrease in temperature. It has been conjectured that the internal topology of the constituent particles, such as polymers, can cause glassiness too. However, the conjecture relies on immobilizing a fraction of the particles and is therefore difficult to fulfill experimentally. Here we show that in dense solutions of circular polymers containing (active) segments of increased mobility, the interplay of the activity and the topology of the polymers generates an unprecedented glassy state of matter. The active isotropic driving enhances mutual ring threading to the extent that the rings can relax only in a cooperative way, which dramatically increases relaxation times. Moreover, the observed phenomena feature similarities with the conformation and dynamics of the DNA fibre in living nuclei of higher eukaryotes.

Glass-transition asymptotics in two theories of glassy dynamics: Self-consistent generalized Langevin equation and mode-coupling theory

We contrast the generic features of structural relaxation close to the idealized glass transition that are predicted by the self-consistent generalized Langevin equation theory (SCGLE) against those that are predicted by the mode-coupling theory of the glass transition (MCT). We present an asymptotic solution close to conditions of kinetic arrest that is valid for both theories, despite the different starting points that are adopted in deriving them. This in particular provides the same level of understanding of the asymptotic dynamics in the SCGLE as was previously done only for MCT. We discuss similarities and different predictions of the two theories for kinetic arrest in standard glass-forming models, as exemplified through the hard-sphere system. Qualitative differences are found for models where a decoupling of relaxation modes is predicted, such as the generalized Gaussian core model, or binary hard-sphere mixtures of particles with very disparate sizes. These differences, which arise in the distinct treatment of the memory kernels associated to self- and collective motion of particles, lead to distinct scenarios that are predicted by each theory for partially arrested states and in the vicinity of higher-order glass-transition singularities.

Dynamically asymmetric binary glass formers studied by dielectric and NMR spectroscopy

We investigate the component dynamics in asymmetric binary glass formers. Focusing on the dielectric spectra of the high-T_g components m-tricresyl phosphate and quinaldine mixed with toluene as low-T_g component, the broadend spectra cannot be described by Kohlrausch or Cole-Davidson (CD) functions. Instead, we apply a generalized CD function which allows to control the width of the susceptibility independently of its high-frequency flank. The spectra show a common broadening and failure of the frequency-temperature superposition with increasing toluene concentration. This is confirmed by stimulated echo experiments showing an increased stretching of the probed orientational correlation function. In analogy to the definition of T_g, we consider "isodynamic points". For each component, a different but linear concentration dependence of 1/T_iso is revealed, indicating different time scales. Qualitativly, we do not find significant differences for the present mixtures with T_g-contrasts of 63-89K compared to those with larger T_g-contrast ( [Formula: see text] K): Whereas the high-T_g component shows relaxation features similar to those of neat glass formers, yet, with "atypical" weak relaxation broadening, the faster low-T_g component displays pronounced dynamic heterogeneities. This is supported by scrutinizing NMR relaxation data of several mixtures investigated previously as a function of concentration. A universal evolution of the dynamics of the high-T_g as well as the low-T_g component is suggested for mixtures with high [Formula: see text]T_g .

Persistent Anti-Correlations in Brownian Dynamics Simulations of Dense Colloidal Suspensions Revealed by Noise Suppression

Transport properties of a hard-sphere colloidal fluid are investigated by Brownian dynamics simulations. We implement a novel algorithm for the time-dependent velocity-autocorrelation function (VACF) essentially eliminating the noise of the bare random motion. The measured VACF reveals persistent anti-correlations manifested by a negative algebraic power-law tail t^{-5/2} at all densities. At small packing fractions the simulations fully agree with the analytic low-density prediction, yet the amplitude of the tail becomes dramatically suppressed as the packing fraction is increased. The mode-coupling theory of the glass transition provides a qualitative explanation for the strong variation in terms of the static compressibility as well as the slowing down of the structural relaxation.

Connect the Thermodynamics of Bulk and Confined Fluids: Confinement-Adsorption Scaling

A fluid (a gas or a liquid) adsorbed in a porous material can behave very differently from its bulk counterpart. The advent of various synthesized materials with nanopores and their wide applications have provided strong impetus for studying fluids in confinement because our current understanding is still incomplete. From a large number of Monte Carlo simulations, we found a scaling relation that allows for connecting some thermodynamic properties (chemical potential, free energy per particle, and grand potential per particle) of a confined fluid to the bulk ones. Upon rescaling the adsorbed fluid density, the adsorption isotherms for many different confining environments collapse to the corresponding bulk curve. We also reveal the intimate connection of the reported scaling relation to Gibbs theory of inhomogeneous fluids and morphological thermodynamics. The advance in our understanding of confined fluids, gained from this study, also opens attractive perspectives for circumventing experimental difficulty for directly measuring some fluid thermodynamic properties in nanoporous materials.

Glassy relaxation slows down by increasing mobility

We find a striking trend reversal in the relaxation dynamics of mixtures with strong dynamical asymmetry. Simulations by both Brownian and Newtonian dynamics reveal that in mixtures of fast and slow hard spheres, above a critical density, the dynamics becomes slower upon increasing the mobility of the fast particles. Below that density, the same increase in mobility speeds up the dynamics. The critical density itself can be identified with the glass transition of the mode-coupling theory that does not depend on the dynamical asymmetry. The asymptotic dynamics close to the critical density is universal, but strong pre-asymptotic effects prevail in particular when the dynamical asymmetry also involves size asymmetry. Our observations reconcile earlier findings, where a strong dependence on kinetic parameters was found for the glassy dynamics, with the paradigm that the glass transition is determined by the properties of configuration space alone.

Different scenarios of dynamic coupling in glassy colloidal mixtures

The composition of mixtures determines the mechanism of glass formation and dynamic coupling of different species.

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.

Dynamics of fluids in quenched-random potential energy landscapes: a mode-coupling theory approach

Motivated by a number of recent experimental and computational studies of the dynamics of fluids plunged in quenched-disordered external fields, we report on a theoretical investigation of this topic within the framework of the mode-coupling theory, based on the simple model of the hard-sphere fluid in a Gaussian random field. The possible dynamical arrest scenarios driven by an increase of the disorder strength and/or of the fluid density are mapped, and the corresponding evolutions of time-dependent quantities typically used for the characterization of anomalous self-diffusion are illustrated with detailed computations. Overall, a fairly reasonable picture of the dynamics of the system at hand is outlined, which in particular involves a non-monotonicity of the self-diffusion coefficient with fluid density at fixed disorder strength, in agreement with experiments. The disorder correlation length is shown to have a strong influence on the latter feature.

Mode coupling theory for nonequilibrium glassy dynamics of thermal self-propelled particles

We present a mode coupling theory study for the relaxation and glassy dynamics of a system of strongly interacting self-propelled particles, wherein the self-propulsion force is described by Ornstein-Uhlenbeck colored noise and thermal noises are included. Our starting point is an effective Smoluchowski equation governing the distribution function of particle positions, from which we derive a memory function equation for the time dependence of density fluctuations in nonequilibrium steady states. With the basic assumption of the absence of macroscopic currents and standard mode coupling approximation, we can obtain expressions for the irreducible memory function and other relevant dynamic terms, wherein the nonequilibrium character of the active system is manifested through an averaged diffusion coefficient D[combining macron] and a nontrivial structural function S₂(q) with q being the magnitude of wave vector q. D[combining macron] and S₂(q) enter the frequency term and the vertex term for the memory function, and thus influence both the short time and the long time dynamics of the system. With these equations obtained, we study the glassy dynamics of this thermal self-propelled particle system by investigating the Debye-Waller factor f_q and relaxation time τ_α as functions of the persistence time τ_p of self-propulsion, the single particle effective temperature T_eff as well as the number density ρ. Consequently, we find the critical density ρ_c for given τ_p shifts to larger values with increasing magnitude of propulsion force or effective temperature, in good accordance with previously reported simulation work. In addition, the theory facilitates us to study the critical effective temperature T for fixed ρ as well as its dependence on τ_p. We find that T increases with τ_p and in the limit τ_p → 0, it approaches the value for a simple passive Brownian system as expected. Our theory also well recovers the results for passive systems and can be easily extended to more complex systems such as active-passive mixtures.

Fluids in porous media. IV. Quench effect on chemical potential

It appears to be a common sense to measure the crowdedness of a fluid system by the densities of the species constituting it. In the present work, we show that this ceases to be valid for confined fluids under some conditions. A quite thorough investigation is made for a hard sphere (HS) fluid adsorbed in a hard sphere matrix (a quench-annealed system) and its corresponding equilibrium binary mixture. When fluid particles are larger than matrix particles, the quench-annealed system can appear much more crowded than its corresponding equilibrium binary mixture, i.e., having a much higher fluid chemical potential, even when the density of each species is strictly the same in both systems, respectively. We believe that the insight gained from this study should be useful for the design of functionalized porous materials.

Dynamic heterogeneities and non-Gaussian behavior in two-dimensional randomly confined colloidal fluids

A binary mixture of superparamagnetic colloidal particles is confined between glass plates such that the large particles become fixed and provide a two-dimensional disordered matrix for the still mobile small particles, which form a fluid. By varying fluid and matrix area fractions and tuning the interactions between the superparamagnetic particles via an external magnetic field, different regions of the state diagram are explored. The mobile particles exhibit delocalized dynamics at small matrix area fractions and localized motion at high matrix area fractions, and the localization transition is rounded by the soft interactions [T. O. E. Skinner et al., Phys. Rev. Lett. 111, 128301 (2013)PRLTAO0031-900710.1103/PhysRevLett.111.128301]. Expanding on previous work, we find the dynamics of the tracers to be strongly heterogeneous and show that molecular dynamics simulations of an ideal gas confined in a fixed matrix exhibit similar behavior. The simulations show how these soft interactions make the dynamics more heterogeneous compared to the disordered Lorentz gas and lead to strong non-Gaussian fluctuations.

Diffusion of a Rouse chain in porous media: A mode-coupling-theory study

We use a kinetic mode-coupling theory (MCT) combining with generalized Langevin equation (GLE) to study the diffusion and conformational dynamics of a bead-spring Rouse chain (RC) dissolved in porous media. The media contains fluid particles and immobile matrix ones wherein the latter leads to the lack of translational invariance. The friction kernel ζ(t) used in the GLE can be obtained directly by adopting a simple density-functional approach in which the density correlators calculated by MCT equations of porous media serve as inputs. Due to cage effects generated by surrounding particles, ζ(t) shows a very long tail memory in the high volume fraction of fluid and matrix. It is found that the long-time center-of-mass diffusion constant D_{CM} of the RC decreases with the increment of volume fraction, influencing more strongly by the matrix particles than by the fluid ones. The auto-correlation function (ACF) of the end-to-end distance fluctuation can also be calculated theoretically based on GLE. Of particular interest is that the power-law region of ACF has a nearly fixed length in logarithmic scale when it shifts to longer time range, with increasing the volume fraction of media particles. Moreover, the effect of lack of translational invariance has been investigated by comparing the results between fluid-matrix and pure fluid cases under identical total volume fraction.

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.

Anomalous dynamics of intruders in a crowded environment of mobile obstacles

Many natural and industrial processes rely on constrained transport, such as proteins moving through cells, particles confined in nanocomposite materials or gels, individuals in highly dense collectives and vehicular traffic conditions. These are examples of motion through crowded environments, in which the host matrix may retain some glass-like dynamics. Here we investigate constrained transport in a colloidal model system, in which dilute small spheres move in a slowly rearranging, glassy matrix of large spheres. Using confocal differential dynamic microscopy and simulations, here we discover a critical size asymmetry, at which anomalous collective transport of the small particles appears, manifested as a logarithmic decay of the density autocorrelation functions. We demonstrate that the matrix mobility is central for the observed anomalous behaviour. These results, crucially depending on size-induced dynamic asymmetry, are of relevance for a wide range of phenomena ranging from glassy systems to cell biology.

The classical Lorentz gas model is widely used to describe constrained transport, but its assumption of an immobile environment is not applicable to many biological and industrial processes. Here, the authors show that the mobility of the matrix induces anomalous, logarithmic dynamics of the confined particles.

Characteristic length scales of the secondary relaxations in glass-forming glycerol

We investigate the secondary relaxations and their link to the main structural relaxation in glass-forming liquids using glycerol as a model system. We analyze the incoherent neutron scattering signal dependence on the scattering momentum transfer, Q , in order to obtain the characteristic length scale for different secondary relaxations. Such a capability of neutron scattering makes it somewhat unique and highly complementary to the traditional techniques of glass physics, such as light scattering and broadband dielectric spectroscopy, which provide information on the time scale, but not the length scales, of relaxation processes. The choice of suitable neutron scattering techniques depends on the time scale of the relaxation of interest. We use neutron backscattering to identify the characteristic length scale of 0.7 Å for the faster secondary relaxation described in the framework of the mode-coupling theory (MCT). Neutron spin-echo is employed to probe the slower secondary relaxation of the excess wing type at a low temperature ( ∼ 1.13T g . The characteristic length scale for this excess wing dynamics is approximately 4.7 Å. Besides the Q -dependence, the direct coupling of neutron scattering signal to density fluctuation makes this technique indispensable for measuring the length scale of the microscopic relaxation dynamics.

Splitting of the Universality Class of Anomalous Transport in Crowded Media

We investigate the emergence of subdiffusive transport by obstruction in continuum models for molecular crowding. While the underlying percolation transition for the accessible space displays universal behavior, the dynamic properties depend in a subtle nonuniversal way on the transport through narrow channels. At the same time, the different universality classes are robust with respect to introducing correlations in the obstacle matrix as we demonstrate for quenched hard-sphere liquids as underlying structures. Our results confirm that the microscopic dynamics can dominate the relaxational behavior even at long times, in striking contrast to glassy dynamics.

Light mediated emergence of surface patterns in azopolymers at low temperatures

Polymer thin films doped with azobenzene molecules do have the ability to organize themselves in spontaneous surface relief gratings (SRG) under irradiation using a single polarized beam. To shed some light on this still unexplained phenomenon, we use a new method that permits us to access experimentally the very first steps of the pattern formation process. By decreasing the temperature, we slow down the formation and organization of patterns, due to the large increase in the viscosity and relaxation time of the azopolymer. As a result, decreasing the temperature allows us to access and study much shorter time scales, in the physical mechanisms underlying the pattern formation, than those previously reported. We find that the patterns organize themselves in sub-structures which size increases with the temperature, following the diffusion coefficient evolution of the material. This result suggests that the pattern formation and organization are mainly governed by diffusive processes, in agreement with some theories of SRG formation. Upon decreasing the temperature further, we observe the emergence of small voids located at the junction of the sub-structures.

Classical Liquids in Fractal Dimension

We introduce fractal liquids by generalizing classical liquids of integer dimensions d=1,2,3 to a noninteger dimension dl. The particles composing the liquid are fractal objects and their configuration space is also fractal, with the same dimension. Realizations of our generic model system include microphase separated binary liquids in porous media, and highly branched liquid droplets confined to a fractal polymer backbone in a gel. Here, we study the thermodynamics and pair correlations of fractal liquids by computer simulation and semianalytical statistical mechanics. Our results are based on a model where fractal hard spheres move on a near-critical percolating lattice cluster. The predictions of the fractal Percus-Yevick liquid integral equation compare well with our simulation results.

Mode-coupling approach for the slow dynamics of a liquid on a spherical substrate

We study the dynamics of a one-component liquid constrained on a spherical substrate, a 2-sphere, and investigate how the mode-coupling theory (MCT) can describe the new features brought by the presence of curvature. To this end we have derived the MCT equations in a spherical geometry. We find that, as seen from the MCT, the slow dynamics of liquids in curved space at low temperature does not qualitatively differ from that of glass-forming liquids in Euclidean space. The MCT predicts the right trend for the evolution of the relaxation slowdown with curvature but is dramatically off at a quantitative level.

On polydispersity and the hard sphere glass transition

We investigate the dynamics of polydisperse hard spheres at high packing fractions ϕ. We use extensive numerical simulations based on an experimentally-realistic particle size distribution (PSD) and compare to commonly-used PSDs such as Gaussian or top hat distribution. We find that the mode of kinetic arrest depends on the PSD's shape and not only on its variance. For the experimentally-realistic PSD we find ageing dynamics even though the density correlators decay fully to zero for ϕ ≥ 0.59. We observe substantial decoupling of the dynamics of the smallest and largest particles. While the smallest particles remain diffusive in all our simulations, a power-law describes the largest-particle diffusion, suggesting an ideal arrest at ϕc ∼ 0.588. The latter is however averted just before ϕc, due to the presence of the mobile smallest particles. In addition, we identify that a partial aging mechanism is at work, whose effects are most pronounced for the largest particles. By comparing our results with recent experimental observations of ergodic behavior up to ϕ ∼ 0.6 in a hard-sphere system, we argue that this is an effect of polydispersity, which smears out the glass transition.

Critical dynamical heterogeneities close to continuous second-order glass transitions

We analyze, using inhomogeneous mode-coupling theory, the critical scaling behavior of the dynamical susceptibility at a distance ε from continuous second-order glass transitions. We find that the dynamical correlation length ξ behaves generically as ε(-1/3) and that the upper critical dimension is equal to six. More surprisingly, we find that ξ grows with time as ln²t exactly at criticality. All of these results suggest a deep analogy between the glassy behavior of attractive colloids or randomly pinned supercooled liquids and that of the random field Ising model.

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.

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.

Applicability of mode-coupling theory to polyisobutylene: a molecular dynamics simulation study

The applicability of Mode Coupling Theory (MCT) to the glass-forming polymer polyisobutylene (PIB) has been explored by using fully atomistic molecular dynamics simulations. MCT predictions for the so-called asymptotic regime have been successfully tested on the dynamic structure factor and the self-correlation function of PIB main-chain carbons calculated from the simulated cell. The factorization theorem and the time-temperature superposition principle are satisfied. A consistent fitting procedure of the simulation data to the MCT asymptotic power-laws predicted for the α-relaxation regime has delivered the dynamic exponents of the theory-in particular, the exponent parameter λ-the critical non-ergodicity parameters, and the critical temperature T(c). The obtained values of λ and T(c) agree, within the uncertainties involved in both studies, with those deduced from depolarized light scattering experiments [A. Kisliuk et al., J. Polym. Sci. Part B: Polym. Phys. 38, 2785 (2000)]. Both, λ and T(c)/T(g) values found for PIB are unusually large with respect to those commonly obtained in low molecular weight systems. Moreover, the high T(c)/T(g) value is compatible with a certain correlation of this parameter with the fragility in Angell's classification. Conversely, the value of λ is close to that reported for real polymers, simulated "realistic" polymers and simple polymer models with intramolecular barriers. In the framework of the MCT, such finding should be the signature of two different mechanisms for the glass-transition in real polymers: intermolecular packing and intramolecular barriers combined with chain connectivity.

Localization dynamics of fluids in random confinement

The dynamics of two-dimensional fluids confined within a random matrix of obstacles is investigated using both colloidal model experiments and molecular dynamics simulations. By varying fluid and matrix area fractions in the experiment, we find delocalized tracer particle dynamics at small matrix area fractions and localized motion of the tracers at high matrix area fractions. In the delocalized region, the dynamics is subdiffusive at intermediate times, and diffusive at long times, while in the localized regime, trapping in finite pockets of the matrix is observed. These observations are found to agree with the simulation of an ideal gas confined in a weakly correlated matrix. Our results show that Lorentz gas systems with soft interactions are exhibiting a smoothening of the critical dynamics and consequently a rounded delocalization-to-localization transition.

Anomalous transport in the crowded world of biological cells

A ubiquitous observation in cell biology is that the diffusive motion of macromolecules and organelles is anomalous, and a description simply based on the conventional diffusion equation with diffusion constants measured in dilute solution fails. This is commonly attributed to macromolecular crowding in the interior of cells and in cellular membranes, summarizing their densely packed and heterogeneous structures. The most familiar phenomenon is a sublinear, power-law increase of the mean-square displacement (MSD) as a function of the lag time, but there are other manifestations like strongly reduced and time-dependent diffusion coefficients, persistent correlations in time, non-Gaussian distributions of spatial displacements, heterogeneous diffusion and a fraction of immobile particles. After a general introduction to the statistical description of slow, anomalous transport, we summarize some widely used theoretical models: Gaussian models like fractional Brownian motion and Langevin equations for visco-elastic media, the continuous-time random walk model, and the Lorentz model describing obstructed transport in a heterogeneous environment. Particular emphasis is put on the spatio-temporal properties of the transport in terms of two-point correlation functions, dynamic scaling behaviour, and how the models are distinguished by their propagators even if the MSDs are identical. Then, we review the theory underlying commonly applied experimental techniques in the presence of anomalous transport like single-particle tracking, fluorescence correlation spectroscopy (FCS) and fluorescence recovery after photobleaching (FRAP). We report on the large body of recent experimental evidence for anomalous transport in crowded biological media: in cyto- and nucleoplasm as well as in cellular membranes, complemented by in vitro experiments where a variety of model systems mimic physiological crowding conditions. Finally, computer simulations are discussed which play an important role in testing the theoretical models and corroborating the experimental findings. The review is completed by a synthesis of the theoretical and experimental progress identifying open questions for future investigation.

Cooperative heterogeneous facilitation: multiple glassy states and glass-glass transition

The formal structure of glass singularities in the mode-coupling theory (MCT) of supercooled liquids dynamics is closely related to that appearing in the analysis of heterogeneous bootstrap percolation on Bethe lattices, random graphs, and complex networks. Starting from this observation one can build up microscopic on-lattice realizations of schematic MCT based on cooperative facilitated spin mixtures. I discuss a microscopic implementation of the F(13) schematic model including multiple glassy states and the glass-glass transition. Results suggest that our approach is flexible enough to bridge alternative theoretical descriptions of glassy matter based on the notions of quenched disorder and dynamic facilitation.

Mode-coupling theory of the glass transition for confined fluids

We present a detailed derivation of a microscopic theory for the glass transition of a liquid enclosed between two parallel walls relying on a mode-coupling approximation. This geometry lacks translational invariance perpendicular to the walls, which implies that the density profile and the density-density correlation function depends explicitly on the distances to the walls. We discuss the residual symmetry properties in slab geometry and introduce a symmetry adapted complete set of two-point correlation functions. Since the currents naturally split into components parallel and perpendicular to the walls the mathematical structure of the theory differs from the established mode-coupling equations in bulk. We prove that the equations for the nonergodicity parameters still display a covariance property similar to bulk liquids.

Signature of a type-A glass transition and intrinsic confinement effects in a binary glass-forming system

We study dynamically highly asymmetric binary mixtures comprised of small methyl tetrahydrofuran (MTHF) molecules and polystyrene. Combined use of dielectric spectroscopy, 2H nuclear magnetic resonance, incoherent quasielastic neutron scattering, and depolarized dynamic light scattering allows us to selectively probe the dynamics of the components in a broad dynamic range. It turns out that the mixtures exhibit two glass transitions in a wide concentration range although being fully miscible on a macroscopic scale. In between both glass transition temperatures, the dynamics of the small molecules show strong confinement effects, e.g., a crossover from Vogel-Fulcher to Arrhenius behavior of the time constants. Moreover, the dynamical behavior of small molecules close to the slow matrix is consistent with mode coupling theory predictions for a type-A glass transition, which was expected from recent theoretical and simulation studies in comparable systems.