Dynamic equivalences in the hard-sphere dynamic universality class

Induced phase transformation in ionizable colloidal nanoparticles

Acid-base equilibria directly influence the functionality and behavior of particles in a system. Due to the ionizing effects of acid-base functional groups, particles will undergo charge exchange. The degree of ionization and their intermolecular and electrostatic interactions are controlled by varying the pH and salt concentration of the solution in a system. Although the pH can be tuned in experiments, it is hard to model this effect using simulations or theoretical approaches. This is due to the difficulty in treating charge regulation and capturing the cooperative effects in a colloidal suspension with Coulombic interaction. In this work, we analyze a suspension of ionizable colloidal particles via molecular dynamics (MD) simulations, along with Monte Carlo simulations for charge regulation (MC-CR) and derive a phase diagram of the system as a function of pH. It is observed that as pH increases, particles functionalized with acid groups change their arrangement from face-centered cubic (FCC) packing to a disordered state. We attribute these transitions to an increase in the degree of charge polydispersity arising from an increase in pH. Our work shows that charge regulation leads to amorphous solids in colloids when the mean nanoparticle charge is sufficiently high.

The origin of the density scaling exponent for polyatomic molecules and the estimation of its value from the liquid structure

In this article, we unravel the problem of interpreting the density scaling exponent for the polyatomic molecules representing the real van der Waals liquids. Our studies show that the density scaling exponent is a weighted average of the exponents of the repulsive terms of all interatomic interactions that occur between molecules, where the potential energy of a given interaction represents its weight. It implies that potential energy is a key quantity required to calculate the density scaling exponent value for real molecules. Finally, we use the well-known method for potential energy estimation and show that the density scaling exponent could be successfully predicted from the liquid structure for fair representatives of the real systems.

Non-equilibrium view of the amorphous solidification of liquids with competing interactions

The interplay between short-range attractions and long-range repulsions (SALR) characterizes the so-called liquids with competing interactions, which are known to exhibit a variety of equilibrium and non-equilibrium phases. The theoretical description of the phenomenology associated with glassy or gel states in these systems has to take into account both the presence of thermodynamic instabilities (such as those defining the spinodal line and the so called λ line) and the limited capability to describe genuine non-equilibrium processes from first principles. Here, we report the first application of the non-equilibrium self-consistent generalized Langevin equation theory to the description of the dynamical arrest processes that occur in SALR systems after being instantaneously quenched into a state point in the regions of thermodynamic instability. The physical scenario predicted by this theory reveals an amazing interplay between the thermodynamically driven instabilities, favoring equilibrium macro- and micro-phase separation, and the kinetic arrest mechanisms, favoring non-equilibrium amorphous solidification of the liquid into an unexpected variety of glass and gel states.

Structural relaxation, dynamical arrest, and aging in soft-sphere liquids

We investigate the structural relaxation of a soft-sphere liquid quenched isochorically (ϕ = 0.7) and instantaneously to different temperatures T_f above and below the glass transition. For this, we combine extensive Brownian dynamics simulations and theoretical calculations based on the non-equilibrium self-consistent generalized Langevin equation (NE-SCGLE) theory. The response of the liquid to a quench generally consists of a sub-linear increase of the α-relaxation time with system's age. Approaching the ideal glass-transition temperature from above (T_f > T^a), sub-aging appears as a transient process describing a broad equilibration crossover for quenches to nearly arrested states. This allows us to empirically determine an equilibration timescale t^eq(T_f) that becomes increasingly longer as T_f approaches T^a. For quenches inside the glass (T_f ≤ T^a), the growth rate of the structural relaxation time becomes progressively larger as T_f decreases and, unlike the equilibration scenario, τ_α remains evolving within the whole observation time-window. These features are consistently found in theory and simulations with remarkable semi-quantitative agreement and coincide with those revealed in a previous and complementary study [P. Mendoza-Méndez et al., Phys. Rev. 96, 022608 (2017)] that considered a sequence of quenches with fixed final temperature T_f = 0 but increasing ϕ toward the hard-sphere dynamical arrest volume fraction ϕ_HS ^a=0.582. The NE-SCGLE analysis, however, unveils various fundamental aspects of the glass transition, involving the abrupt passage from the ordinary equilibration scenario to the persistent aging effects that are characteristic of glass-forming liquids. The theory also explains that, within the time window of any experimental observation, this can only be observed as a continuous crossover.

Structural relaxation dynamics of colloidal nanotrimers

By Molecular Dynamics simulation, we investigate the dynamics of isotropic fluids of colloidal nanotrimers whose interactions are described by varying the strength of attractive and repulsive terms of the Mie potential. To provide a consistent comparison between the systems described by different force fields, we determine the phase diagram and critical points of each system, characterize the morphology of high-density liquid phases at the same reduced temperature and density, and finally investigate their long-time relaxation dynamics. In particular, we detect an especially complex dynamics that reveals the existence of slow and fast nanotrimers and the resulting occurrence of non-Gaussianity, which develops at intermediate timescales. Deviations from Gaussianity are temporary and vanish within the timescales of the system's density fluctuations decay, when a Fickian-like diffusion regime is eventually observed.

Softness mapping of the concentration dependence of the dynamics in model soft colloidal systems

The dynamics of a series of soft colloids comprised of polystyrene cores with poly(N-isopropylacrylamide) (PNIPAM) coronas was investigated by diffusing wave spectroscopy (DWS). The modulus of the coronas was varied by changing the cross-link density and we were able to interpret the results within a hard-soft mapping framework. The soft, swellable particle properties were modeled using an extended Flory-Rehner theory and a Hertzian pair potential. Following volume fraction jumps, softness effects on the concentration dependence of dynamics were determined, with a 'soft colloids make strong glass-forming liquid'-type of behavior observed close to the nominal glass transition volume fraction, φ_g. Such behavior from the current systems cannot be fully explained by the osmotic deswelling model alone. However, inspired by the soft-hard mapping from Schmiedeberg et al, [Europhys. Lett. 2011, 96(3), 36010] we estimated effective hard-sphere diameters and achieved a successful mapping of the α-relaxation times to a master curve below φ_g. Above φ_g, the curves no longer collapse but show strong deviations from a Vogel-Fulcher type of divergence onto soft jamming plateaux. Our results provide evidence that osmotic deswelling itself cannot fully explain the observed dynamics. Softness also plays an important role in the dynamics of soft, concentrated colloids.

Effect of attractions on the hard-sphere dynamic universality class

The fundamental understanding of the dynamic and transport properties of liquids is crucial for the better processing of most materials. The usefulness of this understanding increases when it involves general scaling rules, such as the concept of the hard-sphere dynamic universality class, which provides a unifying scaling of the dynamics of soft-sphere repulsive systems. A relevant question is how far this concept extends to systems that also involve attractive interactions. To answer this question, in this work we performed systematic molecular and Brownian dynamics simulations with the Lennard-Jones system in a wide range of temperatures and densities and verify the extent to which its static and dynamic properties map onto those of the hard-sphere system. We determine that in most of the fluid regime, the Lennard-Jones liquid exhibits the same dynamic equivalence with the hard-sphere system as most purely repulsive fluids, thus establishing the degree of its inclusion in the hard-sphere dynamic universality class.

Understanding simple liquids through statistical and deep learning approaches

Statistical and deep learning-based methods are employed to obtain insights into the quasi-universal properties of simple liquids. In the first part, a statistical model is employed to provide a probabilistic explanation for the similarity in the structure of simple liquids interacting with different pair potential forms, collectively known as simple liquids. The methodology works by sampling the radial distribution function and the number of interacting particles within the cutoff distance, and it produces the probability density function of the net force. We show that matching the probability distribution of the net force can be a direct route to parameterize simple liquid pair potentials with a similar structure, as the net force is the main component of the Newtonian equations of motion. The statistical model is assessed and validated against various cases. In the second part, we exploit DeepILST [A. Moradzadeh and N. R. Aluru, J. Phys. Chem. Lett. 10, 1242-1250 (2019)], a data-driven and deep-learning assisted framework to parameterize the standard 12-6 Lennard-Jones (LJ) pair potential, to find structurally equivalent/isomorphic LJ liquids that identify constant order parameter [τ=∫₀ ^ξ_cf gξ-1ξ²dξ, where gξ and ξ(=rρ¹³) are the reduced radial distribution function and radial distance, respectively] systems in the space of non-dimensional temperature and density of the LJ liquids. We also investigate the consistency of DeepILST in reproducibility of radial distribution functions of various quasi-universal potentials, e.g., exponential, inverse-power-law, and Yukawa pair potentials, quantified based on the radial distribution functions and Kullback-Leibler errors. Our results provide insights into the quasi-universality of simple liquids using the statistical and deep learning methods.

Transport properties of Lennard-Jones fluids: Freezing density scaling along isotherms

It is demonstrated that properly reduced transport coefficients (self-diffusion, shear viscosity, and thermal conductivity) of Lennard-Jones fluids along isotherms exhibit quasi-universal scaling on the density divided by its value at the freezing point. Moreover, this scaling is closely related to the density scaling of transport coefficients of hard-sphere fluids. The Stokes-Einstein relation without the hydrodynamic diameter is valid in the dense fluid regime. The lower density boundary of its validity can serve as a practical demarcation line between gaslike and liquidlike regimes.

Universality in Driven and Equilibrium Hard Sphere Liquid Dynamics

We demonstrate that the time evolution of the van Hove dynamical pair correlation function is governed by adiabatic forces that arise from the free energy and by superadiabatic forces that are induced by the flow of the van Hove function. The superadiabatic forces consist of drag, viscous, and structural contributions, as occur in active Brownian particles, in liquids under shear and in lane forming mixtures. For hard sphere liquids, we present a power functional theory that predicts these universal force fields in quantitative agreement with our Brownian dynamics simulation results.

Brownian motion of ellipsoidal particles on a granular magnetic bath

We study the Brownian motion of ellipsoidal particles lying on an agitated granular bath composed of magnetic particles. We quantify the mobility of different floating ellipsoidal particles using the mean-square displacement and the mean-square angular displacement, and relate the diffusion coefficients to the bath particle motion. In terms of the particle major radius R, we find the translational diffusion coefficient scales roughly as 1/R^{2} and the rotational diffusion coefficient scales as roughly 1/R^{4}; this is consistent with the assumption that diffusion arises from random kicks of the bath particles underneath the floating particle. By varying the magnetic forcing, the bath particles' diffusivity changes by a factor of ten; over this range, the translational and rotational diffusion of the floating particles change by a factor of 50. However, the ratio of the two diffusion constants for the floating particles is forcing-independent. Unusual aspects of the floating particle motion include non-Gaussian statistics for their displacements.

Exploring the connection between the density-scaling exponent and the intermolecular potential for liquids on the basis of computer simulations of quasireal model systems

In this paper, based on the molecular dynamics simulations of quasireal model systems, we propose a method for determination of the effective intermolecular potential for real materials. We show that in contrast to the simple liquids, the effective intermolecular potential for the studied systems depends on the thermodynamic conditions. Nevertheless, the previously established relationship for simple liquids between the exponent of the inverse power law approximation of intermolecular potential and the density-scaling exponent is still preserved when small enough intermolecular distances are considered. However, our studies show that molecules approach each other at these very short distances relatively rarely. Consequently, only sparse interactions between extremely close molecules determine the value of the scaling exponent and then strongly influence the connection between dynamics and thermodynamics of the whole system.

Temperature dependence of the transition packing fraction of thermal jamming in a harmonic soft sphere system

The glassy dynamics of soft harmonic spheres are often mapped onto the dynamics of hard spheres by considering an effective diameter for the soft particles and therefore an effective packing fraction. While in this approach the thermal fluctuations within valleys of the energy landscape are covered, the crossing of energy barriers from one valley into another usually is neglected. Here we argue-motivated by studies of the glass transition based on explorations of the energy landscape-that the crossing of energy barriers can be attributed by an effective decrease of the glass transition packing fraction with increasing temperature T according to T ^0.2. Furthermore, we reanalyzing data of soft sphere simulations. Since fitting scaling laws to simulation data always allows for some arbitrariness, we cannot prove based on the simulation data that our idea of a shift of the glass transition packing fraction due to barrier crossings is the only possible way to explain the discrepancies that have been observed previously. However, we show that a possible explanation of the simulation data is given by our approach to characterize the dynamics of soft spheres by both, the previously-considered temperature-dependent effective packing fraction due to the increase of the mean overlap between neighboring particles with stronger thermal fluctuations and the newly introduced increase of the glass transition packing with an increasing probability of barrier crossings.

Glassy dynamics in asymmetric binary mixtures of hard spheres

We perform a systematic and detailed study of the glass transition in highly asymmetric binary mixtures of colloidal hard spheres, combining differential dynamic microscopy experiments, event-driven molecular dynamics simulations, and theoretical calculations, exploring the whole state diagram and determining the self-dynamics and collective dynamics of both species. Two distinct glassy states involving different dynamical arrest transitions are consistently described, namely, a double glass with the simultaneous arrest of the self-dynamics and collective dynamics of both species, and a single glass of large particles in which the self-dynamics of the small species remains ergodic. In the single-glass scenario, spatial modulations in the collective dynamics of both species occur due to the structure of the large spheres, a feature not observed in the double-glass domain. The theoretical results, obtained within the self-consistent generalized Langevin equation formalism, are in agreement with both simulations and experimental data, thus providing a stringent validation of this theoretical framework in the description of dynamical arrest in highly asymmetric mixtures. Our findings are summarized in a state diagram that classifies the various amorphous states of highly asymmetric mixtures by their dynamical arrest mechanisms.

The effect of molecular architecture on the physical properties of supercooled liquids studied by MD simulations: Density scaling and its relation to the equation of state

Theoretical concepts in condensed matter physics are typically verified and also developed by exploiting computer simulations mostly in simple models. Predictions based on these usually isotropic models are often at odds with measurement results obtained for real materials. One of the examples is an intriguing problem within the density scaling idea that has attracted attention in recent decades due to its hallmarks of universality, i.e., the fact that the difference between the density scaling exponent and the exponent of the equation of state is observed for real materials, whereas it has not been reported for the model system. In this paper, we use new model molecules of simple but anisotropic architecture to study the effect of molecular anisotropy on the dynamic and thermodynamic properties of the system. We identify the applicable range of intermolecular interactions for a given physical process, and then we explain the reason for observed differences between the behavior of the model and real systems. It demonstrates that the new model systems open broad perspectives for simulation and theoretical research, for example, into unifying concepts in the glass transition physics.

Crystallization of hard spheres revisited. I. Extracting kinetics and free energy landscape from forward flux sampling

We investigate the kinetics and the free energy landscape of the crystallization of hard spheres from a supersaturated metastable liquid though direct simulations and forward flux sampling. In this first paper, we describe and test two different ways to reconstruct the free energy barriers from the sampled steady state probability distribution of cluster sizes without sampling the equilibrium distribution. The first method is based on mean first passage times, and the second method is based on splitting probabilities. We verify both methods for a single particle moving in a double-well potential. For the nucleation of hard spheres, these methods allow us to probe a wide range of supersaturations and to reconstruct the kinetics and the free energy landscape from the same simulation. Results are consistent with the scaling predicted by classical nucleation theory although a quantitative fit requires a rather large effective interfacial tension.

Crossover from equilibration to aging: Nonequilibrium theory versus simulations

Understanding glasses and the glass transition requires comprehending the nature of the crossover from the ergodic (or equilibrium) regime, in which the stationary properties of the system have no history dependence, to the mysterious glass transition region, where the measured properties are nonstationary and depend on the protocol of preparation. In this work we use nonequilibrium molecular dynamics simulations to test the main features of the crossover predicted by the molecular version of the recently developed multicomponent nonequilibrium self-consistent generalized Langevin equation theory. According to this theory, the glass transition involves the abrupt passage from the ordinary pattern of full equilibration to the aging scenario characteristic of glass-forming liquids. The same theory explains that this abrupt transition will always be observed as a blurred crossover due to the unavoidable finiteness of the time window of any experimental observation. We find that within their finite waiting-time window, the simulations confirm the general trends predicted by the theory.

Cage Size and Jump Precursors in Glass-Forming Liquids: Experiment and Simulations

Glassy dynamics is intermittent, as particles suddenly jump out of the cage formed by their neighbors, and heterogeneous, as these jumps are not uniformly distributed across the system. Relating these features of the dynamics to the diverse local environments explored by the particles is essential to rationalize the relaxation process. Here we investigate this issue characterizing the local environment of a particle with the amplitude of its short time vibrational motion, as determined by segmenting in cages and jumps the particle trajectories. Both simulations of supercooled liquids and experiments on colloidal suspensions show that particles in large cages are likely to jump after a small time-lag, and that, on average, the cage enlarges shortly before the particle jumps. At large time-lags, the cage has essentially a constant size, which is smaller for longer-lasting cages. Finally, we clarify how this coupling between cage size and duration controls the average behavior and opens the way to a better understanding of the relaxation process in glass-forming liquids.

Anisotropy and memory during cage breaking events close to a wall

The slow dynamics in a glassy hard-sphere system is dominated by cage breaking events, i.e. rearrangements where a particle escapes from the cage formed by its neighboring particles. We study such events for an overdamped colloidal system by the means of Brownian dynamics simulations. While it is difficult to relate cage breaking events to structural mean field results in bulk, we show that the microscopic dynamics of particles close to a wall can be related to the anisotropic two-particle density. In particular, we study cage-breaking trajectories, mean forces on a tracked particle, and the impact of the history of trajectories. Based on our simulation results, we further construct two different one-particle random-walk models-one without and one with memory incorporated-and find the local anisotropy and the history-dependence of particles as crucial ingredients to describe the escape from a cage. Finally, our detailed study of a rearrangement event close to a wall not only reveals the memory effect of cages, but leads to a deeper insight into the fundamental mechanisms of glassy dynamics.

Nonvibrating granular model for a glass-forming liquid: Equilibration and aging

We studied experimentally a model of a glass-forming liquid on the basis of a nonvibrating magnetic granular system under an unsteady magnetic field. A sudden quenching was produced that drove the system from a liquid state to a different final state with lower temperature; the latter could be a liquid state or a solid state. We determined the mean-squared displacement in temporal windows to obtain the dynamic evolution of the system, and we determined the radial distribution function to obtain its structural characteristics. The results were analyzed using the intermediate scattering function and the effective potential between two particles. We observed that when quenching drives the system to a final state in the liquid phase far from the glass-transition temperature, equilibration occurs very quickly. When the final state has a temperature far below the glass-transition temperature, the system reaches its equilibrium state very quickly. In contrast, when the final state has an intermediate temperature but is below that corresponding to the glass transition, the system falls into a state that evolves slowly, presenting aging. The system evolves by an aging process toward more ordered states. However, after a waiting time, the dynamic behavior changes. It was observed that some particles get close enough to overpass the repulsive interactions and form small stable aggregates. In the effective potential curves, it was observed that the emergence of a second effective well due to the attraction quickly evolves and results in a deeper well than the initial effective well due to the repulsion. With the increase in time, more particles fall in the attractive well forming inhomogeneities, which produce a frustration in the aging process.

Simple liquids' quasiuniversality and the hard-sphere paradigm

This topical review discusses the quasiuniversality of simple liquids' structure and dynamics and two possible justifications of it. The traditional one is based on the van der Waals picture of liquids in which the hard-sphere system reflects the basic physics. An alternative explanation argues that all quasiuniversal liquids to a good approximation conform to the same equation of motion, referring to the exponentially repulsive pair-potential system as the basic reference system. The paper, which is aimed at non-experts, ends by listing a number of open problems in the field.

Short-time dynamics of monomers and dimers in quasi-two-dimensional colloidal mixtures

We report on the short-time dynamics in colloidal mixtures made up of monomers and dimers highly confined between two glass plates. At low concentrations, the experimental measurements of colloidal motion agree well with the solution of the Navier-Stokes equation at low Reynolds numbers; the latter takes into account the increase in the drag force on a colloidal particle due to wall-particle hydrodynamic forces. More importantly, we find that the ratio of the short-time diffusion coefficient of the monomer and that of the center of mass of the dimmer is almost independent of both the dimer molar fraction, x_{d}, and the total packing fraction, ϕ, up to ϕ≈0.5. At higher concentrations, this ratio displays a small but systematic increase. A similar physical scenario is observed for the ratio between the parallel and the perpendicular components of the short-time diffusion coefficients of the dimer. This dynamical behavior is corroborated by means of molecular dynamics computer simulations that include explicitly the particle-particle hydrodynamic forces induced by the solvent. Our results suggest that the effects of colloid-colloid hydrodynamic interactions on the short-time diffusion coefficients are almost identical and factorable in both species.

Glassy dynamics of model colloidal polymers: The effect of "monomer" size

In recent years, attempts have been made to assemble colloidal particles into chains, which are termed "colloidal polymers." An apparent difference between molecular and colloidal polymers is the "monomer" size. Here, we propose a model to represent the variation from molecular polymer to colloidal polymer and study the quantitative differences in their glassy dynamics. For chains, two incompatible local length scales, i.e., monomer size and bond length, are manifested in the radial distribution function and intramolecular correlation function. The mean square displacement of monomers exhibits Rouse-like sub-diffusion at intermediate time/length scale and the corresponding exponent depends on the volume fraction and the monomer size. We find that the threshold volume fraction at which the caging regime emerges can be used as a rescaling unit so that the data of localization length versus volume fraction for different monomer sizes can gather close to an exponential curve. The increase of monomer size effectively increases the hardness of monomers and thus makes the colloidal polymers vitrify at lower volume fraction. Static and dynamic equivalences between colloidal polymers of different monomer sizes have been discussed. In the case of having the same peak time of the non-Gaussian parameter, the motion of monomers of larger size is much less non-Gaussian. The mode-coupling critical exponents for colloidal polymers are in agreement with that of flexible bead-spring chains.

Anisotropic pair correlations in binary and multicomponent hard-sphere mixtures in the vicinity of a hard wall: A combined density functional theory and simulation study

The fundamental measure approach to classical density functional theory has been shown to be a powerful tool to predict various thermodynamic properties of hard-sphere systems. We employ this approach to determine not only one-particle densities but also two-particle correlations in binary and six-component mixtures of hard spheres in the vicinity of a hard wall. The broken isotropy enables us to carefully test a large variety of theoretically predicted two-particle features by quantitatively comparing them to the results of Brownian dynamics simulations. Specifically, we determine and compare the one-particle density, the total correlation functions, their contact values, and the force distributions acting on a particle. For this purpose, we follow the compressibility route and theoretically calculate the direct correlation functions by taking functional derivatives. We usually observe an excellent agreement between theory and simulations, except for small deviations in cases where local crystal-like order sets in. Our results set the course for further investigations on the consistency of functionals as well as for structural analysis on, e.g., the primitive model. In addition, we demonstrate that due to the suppression of local crystallization, the predictions of six-component mixtures are better than those in bidisperse or monodisperse systems. Finally, we are confident that our results of the structural modulations induced by the wall lead to a deeper understanding of ordering in anisotropic systems in general, the onset of heterogeneous crystallization, caging effects, and glassy dynamics close to a wall, as well as structural properties in systems with confinement.

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.

Self-consistent generalized Langevin-equation theory for liquids of nonspherically interacting particles

A self-consistent generalized Langevin-equation theory is proposed to describe the self- and collective dynamics of a liquid of linear Brownian particles. The equations of motion for the spherical harmonics projections of the collective and self-intermediate-scattering functions, F_{lm,lm}(k,t) and F_{lm,lm}^{S}(k,t), are derived as a contraction of the description involving the stochastic equations of the corresponding tensorial one-particle density n_{lm}(k,t) and the translational (α=T) and rotational (α=R) current densities j_{lm}^{α}(k,t). Similar to the spherical case, these dynamic equations require as an external input the equilibrium structural properties of the system contained in the projections of the static structure factor, denoted by S_{lm,lm}(k). Complementing these exact equations with simple (Vineyard-like) approximate relations for the collective and the self-memory functions we propose a closed self-consistent set of equations for the dynamic properties involved. In the long-time asymptotic limit, these equations become the so-called bifurcation equations, whose solutions (the nonergodicity parameters) can be written, extending the spherical case, in terms of one translational and one orientational scalar dynamic order parameter, γ_{T} and γ_{R}, which characterize the possible dynamical arrest transitions of the system. As a concrete illustrative application of this theory we determine the dynamic arrest diagram of the dipolar hard-sphere fluid. In qualitative agreement with mode coupling theory, the present self-consistent equations also predict three different regions in the state space spanned by the macroscopic control parameters η (volume fraction) and T* (scaled temperature): a region of fully ergodic states, a region of mixed states, in which the translational degrees of freedom become arrested while the orientational degrees of freedom remain ergodic, and a region of fully nonergodic states.

Hidden scale invariance in condensed matter

Recent developments show that many liquids and solids have an approximate "hidden" scale invariance that implies the existence of lines in the thermodynamic phase diagram, so-called isomorphs, along which structure and dynamics in properly reduced units are invariant to a good approximation. This means that the phase diagram becomes effectively one-dimensional with regard to several physical properties. Liquids and solids with isomorphs include most or all van der Waals bonded systems and metals, as well as weakly ionic or dipolar systems. On the other hand, systems with directional bonding (hydrogen bonds or covalent bonds) or strong Coulomb forces generally do not exhibit hidden scale invariance. The article reviews the theory behind this picture of condensed matter and the evidence for it coming from computer simulations and experiments.