Structures and Dynamics of Glass-Forming Colloidal Liquids under Spherical Confinement

Orientational Order in Dense Colloidal Liquids and Glasses

We construct structural order parameters based on local angular and radial distribution functions in dense colloidal suspensions. All the order parameters show significant correlations to local dynamics in the supercooled and glass regime. In particular, the correlations between the orientational order and dynamical heterogeneity are consistently higher than those between the conventional two-body structural entropy and local dynamics. The structure-dynamics correlations can be explained by a excitation model with the energy barrier depending on local structural order. Our results suggest that in dense disordered packings, local orientational order is higher than translational order, and plays a more important role in determining the dynamics in glassy systems.

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.

Long-Wavelength Fluctuations in Quasi-2D Supercooled Liquids

We use molecular simulation to characterize the dynamics of supercooled liquids confined in quasi-2D slit geometries. Similar to bulk supercooled liquids, the confined systems exhibit subdiffusive dynamics on intermediate time scales arising from particle localization inside their neighbor cages, followed by an eventual crossover to diffusive behavior as cage rearrangement occurs. The quasi-2D confined liquids also exhibit signatures of long-wavelength fluctuations (LWFs) in the lateral directions parallel to the confining walls, reminiscent of the collective displacements observed in 2D but not 3D systems. The magnitude of the LWFs increases with the lateral dimensions of systems with the same particle volume fraction and confinement length scale, consistent with the logarithmic scaling predicted for 2D Mermin-Wagner fluctuations. The amplitude of the fluctuations is a nonmonotonic function of the confinement length scale because of a competition between caging and strengthening LWFs upon approaching the 2D limit. Our findings suggest that LWFs may play an important role in understanding the behavior of confined supercooled liquids due to their prevalence over a surprisingly broad range of particle densities and confinement length scales.

Nonmonotonic Dynamical Correlations beneath the Surface of Glass-Forming Liquids

Collective motion over increasing length scales is a signature of the vitrification process of liquids. We demonstrate how distinct static and dynamic length scales govern the dynamics of vitrifying films. In contrast to a monotonically growing static correlation length, the dynamical correlation length that measures the extent of surface-dynamics acceleration into the bulk displays a striking nonmonotonic temperature evolution that is robust also against changes in detailed interatomic interaction. This nonmonotonic change defines a crossover temperature T_{*} that is distinct from the critical temperature T_{c} of mode-coupling theory. We connect this nonmonotonic change to a morphological change of cooperative rearrangement regions of fast particles, and to the point where the decoupling of fast-particle motion from the bulk relaxation is most sensitive to fluctuations. We propose a rigorous definition of this new crossover temperature T_{*} within a recent extension of mode-coupling theory, the stochastic β-relaxation theory.

Evidence for growing structural correlation length in colloidal supercooled liquids

Using video microscopy, we measure the long-time diffusion coefficients of colloidal particles at different concentrations. The measured diffusion coefficients start to deviate from theoretical predictions based on random collision models upon entering the supercooled regime. The theoretical diffusion relation is recovered by assigning an effective mass proportional to the size of structurally correlated clusters to the diffusing particles, providing an indirect method to probe the growth of static correlation length scales approaching the glass transition. This method is tested and validated in the crystallization of mono-disperse colloids in quasi-two-dimensional experiments. The correlation length obtained for a binary colloidal liquid increases by a power law toward a critical packing fraction of ∼0.79. The system relaxation time exhibits a power-law dependence on the correlation length in agreement with dynamical facilitation theories.

Hierarchical glass transition of hard hemidisks with local assemblies

Using computer simulation, we investigate the glass transition of a two-dimensional hard-hemidisk system. Upon increasing the packing fraction of the system, we find that the system vitrifies into a glass with local assembled discal "dimers", which are free to rotate in a collective way. The rotational mean square displacement does not exhibit the typical plateau (slowdown) like what occurs in the translational mean square displacement. This effect induces a pronounced violation of the rotational Stokes-Einstein relationship compared with the translational degree of freedom at the supercooled region. However, the obtained glass transition points in these two freedom degrees are found to be the same within the numerical accuracy, which is due to the strong positive spatial and dynamic correlation between translational and rotational slow-moving particles. Moreover, we find that the locally assembled dimers can serve as fast rotating gears facilitating the orientational relaxation in the system, and this suggests that the locally favored finite structures play an important role in the hierarchical glass transition of anisotropic colloids.

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.

Quantification of spatio-temporal scales of dynamical heterogeneity of water near lipid membranes above supercooling

A hydrated 1,2-dimyristoyl-sn-glycero-3-phosphorylcholine (DMPC) lipid membrane is investigated using an all atom molecular dynamics simulation at 308 K to determine the physical sources of universal slow relaxations of hydration layers and length-scale of the spatially heterogeneous dynamics. Continuously residing interface water (IW) molecules hydrogen bonded to different moieties of lipid heads in the membrane are identified. The non-Gaussian parameters of all classes of IW molecules show a cross-over from cage vibration to translational diffusion. A significant non-Gaussianity is observed for the IW molecules exhibiting large length correlations in translational van Hove functions. Two time-scales for the ballistic motions and hopping transitions are obtained from the self intermediate scattering functions of the IW molecules with an additional long relaxation, which disappears for bulk water. The long relaxation time-scales for the IW molecules obtained from the self intermediate scattering functions are in good accordance with the hydrogen bond relaxation time-scales irrespective of the nature of the chemical confinement and the confinement lifetime. Employing a block analysis approach, the length-scale of dynamical heterogeneities is captured from a transition from non-Gaussianity to Gaussianity in van Hove correlation functions of the IW molecules. The heterogeneity length-scale is comparable to the wave-length of the small and weak undulations of the membrane calculated by Fourier transforms of lipid tilts. This opens up a new avenue towards a possible correlation between heterogeneity length-scale and membrane curvature more significant for rippled membranes. Thus, our analyses provide a measure towards the spatio-temporal scale of dynamical heterogeneity of confined water near membranes.

Six-fold director field configuration in amyloid nematic and cholesteric phases

Chiral liquid crystals, or cholesteric phases, have been widely studied in the last decades, leading to fundamental advances and a multitude of applications and technologies. In general, the rich phenomenology of these systems depends directly on the molecular traits and conditions of the system, imposing precise symmetry to the resulting nematic field. By selecting amyloid fibrils as model filamentous chiral colloids, we report an unprecedented breadth of liquid crystalline morphologies, where up to six distinct configurations of the nematic field are observed under identical conditions. Amyloid-rich droplets show homogeneous, bipolar, radial, uniaxial chiral and radial chiral nematic fields, with additional parabolic focal conics in bulk. Variational and scaling theories allow rationalizing the experimental evidence as a subtle interplay between surface and bulk energies. Our experimental and theoretical findings deepen the understanding of chiral liquid crystals under confinement, opening to a more comprehensive exploitation of these systems in related functional materials.

Long-wavelength fluctuations and static correlations in quasi-2D colloidal suspensions

Dimensionality strongly affects thermal fluctuations and critical dynamics of equilibrium systems. These influences persist in amorphous systems going through the nonequilibrium glass transition. Here, we experimentally study the glass transition of quasi-2D suspensions of spherical and ellipsoidal particles under different degrees of circular confinement. We show that the strength of the long-wavelength fluctuations increases logarithmically with system sizes and displays the signature of the Mermin-Wagner fluctuations. Moreover, using confinement as a tool, we also measure static structural correlations and extract a growing static correlation length in 2D supercooled liquids. Finally, we explore the influence of the Mermin-Wagner fluctuations on the translational and orientational relaxations of 2D ellipsoidal suspensions, which leads to a new interpretation of the two-step glass transition and the orientational glass phase of anisotropic particles. Our study reveals the importance of long-wavelength fluctuations in 2D supercooled liquids and provides new insights into the role of dimensionality in the glass transition.

Glass Transition in Supercooled Liquids with Medium-Range Crystalline Order

The origin of the rapid dynamical slowdown in glass forming liquids in the growth of static length scales, possibly associated with identifiable structural ordering, is a much debated issue. Growth of medium range crystalline order (MRCO) has been observed in various model systems to be associated with glassy behavior. Such observations raise the question of whether molecular mechanisms for the glass transition in liquids with and without MRCO are the same. In this study we perform extensive molecular dynamics simulations of a number of glass forming liquids and show that the static and dynamic properties of glasses with MRCO are different from those of other glass forming liquids with no predominant local order. We also resolve an important issue regarding the so-called point-to-set method for determining static length scales, and demonstrate it to be a robust method for determining static correlation lengths in glass formers.

Local structure in deeply supercooled liquids exhibits growing lengthscales and dynamical correlations

Glasses are among the most widely used of everyday materials, yet the process by which a liquid's viscosity increases by 14 decades to become a glass remains unclear, as often contradictory theories provide equally good descriptions of the available data. Knowledge of emergent lengthscales and higher-order structure could help resolve this, but this requires time-resolved measurements of dense particle coordinates-previously only obtained over a limited time interval. Here we present an experimental study of a model colloidal system over a dynamic window significantly larger than previous measurements, revealing structural ordering more strongly linked to dynamics than previously found. Furthermore we find that immobile regions and domains of local structure grow concurrently with density, and that these regions have low configurational entropy. We thus show that local structure plays an important role at deep supercooling, consistent with a thermodynamic interpretation of the glass transition rather than a principally dynamic description.

Statistical analysis of phase formation in 2D colloidal systems

Colloidal systems offer unique opportunities for the study of phase formation and structure since their characteristic length scales are accessible to visible light. As a model system the two-dimensional assembly of colloidal magnetic and non-magnetic particles dispersed in a ferrofluid (FF) matrix is studied by transmission optical microscopy. We present a method to statistically evaluate images with thousands of particles and map phases by extraction of local variables. Different lattice structures and long-range connected branching chains are observed, when tuning the effective magnetic interaction and varying particle ratios.

Block Analysis for the Calculation of Dynamic and Static Length Scales in Glass-Forming Liquids

We present block analysis, an efficient method of performing finite-size scaling for obtaining the length scale of dynamic heterogeneity and the point-to-set length scale for generic glass-forming liquids. This method involves considering blocks of varying sizes embedded in a system of a fixed (large) size. The length scale associated with dynamic heterogeneity is obtained from a finite-size scaling analysis of the dependence of the four-point dynamic susceptibility on the block size. The block size dependence of the variance of the α relaxation time yields the static point-to-set length scale. The values of the obtained length scales agree quantitatively with those obtained from other conventional methods. This method provides an efficient experimental tool for studying the growth of length scales in systems such as colloidal glasses for which performing finite-size scaling by carrying out experiments for varying system sizes may not be feasible.

Kinetic-contact-driven gigantic energy transfer in a two-dimensional Lennard-Jones fluid confined to a rotating pore

A two-dimensional Lennard-Jones system in a circular and rotating container has been studied by means of molecular dynamics technique. A nonequilibrium transition to the rotating stage has been detected in a delayed time since an instant switching of the frame rotation. This transition is attributed to the increase of the density at the wall because of the centrifugal force. At the same time the phase transition occurs, the inner system changes its configuration of the solid-state type into the liquid type. Impact of angular frequency and molecular roughness on the transport properties of the nonrotating and rotating systems is analyzed.

Pinning susceptibility: a novel method to study growth of amorphous order in glass-forming liquids

The existence and growth of amorphous order in supercooled liquids approaching glass transition is a subject of intense research. Even after decades of work, there is still no clear consensus on the molecular mechanisms that lead to a rapid slowing down of liquid dynamics approaching this putative transition. The existence of a correlation length associated with amorphous order has recently been postulated and has also been estimated using multi-point correlation functions which cannot be calculated easily in experiments. Thus the study of growing amorphous order remains mostly restricted to systems like colloidal glasses and simulations of model glass-forming liquids. In this work, we propose an experimentally realizable yet simple susceptibility to study the growth of amorphous order. We then demonstrate the validity of this approach for a few well-studied model supercooled liquids and obtain results which are consistent with other conventional methods.

Triple Junction at the Triple Point Resolved on the Individual Particle Level

At the triple point of a repulsive screened Coulomb system, a fcc crystal, a bcc crystal, and a fluid phase coexist. At their intersection, these three phases form a liquid groove, the triple junction. Using confocal microscopy, we resolve the triple junction on a single-particle level in a model system of charged PMMA colloids in a nonpolar solvent. The groove is found to be extremely deep and the incommensurate solid-solid interface to be very broad. Thermal fluctuations hence appear to dominate the solid-solid interface. This indicates a very low interfacial energy. The fcc-bcc interfacial energy is quantitatively determined based on Young's equation and, indeed, it is only about 1.3 times higher than the fcc-fluid interfacial energy close to the triple point.

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.

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.

Introduction to the Colloidal Glass Transition

Colloids are suspensions of small solid particles in a liquid and exhibit glassy behavior when the particle concentration is high. In these samples, the particles are roughly analogous to individual molecules in a traditional glass. This model system has been used to study the glass transition since the 1980s. In this Viewpoint I summarize some of the intriguing behaviors of the glass transition in colloids and discuss open questions.

Cooperative strings in glassy nanoparticles

Motivated by recent experimental results on glassy polymer nanoparticles, we develop a minimal theoretical framework for the glass transition in spherical confinement. This is accomplished using our cooperative-string model for supercooled dynamics, that was successful at recovering the bulk phenomenology and describing the thin-film anomalies. In particular, we obtain predictions for the mobile-layer thickness as a function of temperature, and for the effective glass-transition temperature as a function of the radius of the spherical nanoparticle - including the existence of a critical particle radius below which vitrification never occurs. Finally, we compare the theoretical results to experimental data on polystyrene from the recent literature, and we discuss the latter.

Colloidal cholesteric liquid crystal in spherical confinement

The organization of nanoparticles in constrained geometries is an area of fundamental and practical importance. Spherical confinement of nanocolloids leads to new modes of packing, self-assembly, phase separation and relaxation of colloidal liquids; however, it remains an unexplored area of research for colloidal liquid crystals. Here we report the organization of cholesteric liquid crystal formed by nanorods in spherical droplets. For cholesteric suspensions of cellulose nanocrystals, with progressive confinement, we observe phase separation into a micrometer-size isotropic droplet core and a cholesteric shell formed by concentric nanocrystal layers. Further confinement results in a transition to a bipolar planar cholesteric morphology. The distribution of polymer, metal, carbon or metal oxide nanoparticles in the droplets is governed by the nanoparticle size and yields cholesteric droplets exhibiting fluorescence, plasmonic properties and magnetic actuation. This work advances our understanding of how the interplay of order, confinement and topological defects affects the morphology of soft matter.