Superadiabatic dynamical density functional study of Brownian hard-spheres in time-dependent external potentials

Superadiabatic dynamical density functional theory (superadiabatic-DDFT), a first-principles approach based on inhomogeneous two-body correlation functions, is employed to investigate the response of interacting Brownian particles to time-dependent external driving. Predictions for the superadiabatic dynamics of the one-body density are made directly from the underlying interparticle interactions without the need for either adjustable fit parameters or simulation input. The external potentials we investigate have been chosen to probe distinct aspects of structural relaxation in dense, strongly interacting liquid states. Nonequilibrium density profiles predicted by the superadiabatic theory are compared with those obtained from both adiabatic DDFT and event-driven Brownian dynamics simulation. Our findings show that superadiabatic-DDFT accurately predicts the time-evolution of the one-body density.

First-principles superadiabatic theory for the dynamics of inhomogeneous fluids

For classical many-body systems subject to Brownian dynamics, we develop a superadiabatic dynamical density functional theory (DDFT) for the description of inhomogeneous fluids out-of-equilibrium. By explicitly incorporating the dynamics of the inhomogeneous two-body correlation functions, we obtain superadiabatic forces directly from the microscopic interparticle interactions. We demonstrate the importance of these nonequilibrium forces for an accurate description of the one-body density by numerical implementation of our theory for three-dimensional hard-spheres in a time-dependent planar potential. The relaxation of the one-body density in superadiabatic-DDFT is found to be slower than that predicted by standard adiabatic DDFT and significantly improves the agreement with Brownian dynamics simulation data. We attribute this improved performance to the correct treatment of structural relaxation within the superadiabatic-DDFT. Our approach provides fundamental insight into the underlying structure of dynamical density functional theories and makes possible the study of situations for which standard approaches fail.

Fundamental measure theory of inhomogeneous two-body correlation functions

For the three-dimensional hard-sphere model we investigate the inhomogeneous two-body correlations predicted by Rosenfeld's fundamental measure theory. For the special cases in which the density has either planar or spherical symmetry we provide analytic formulas for the Hankel and Legendre transforms, respectively, of the inhomogeneous two-body direct correlation function as explicit functionals of the density. When combined with the Ornstein-Zernike relation our analytical results allow for rapid calculation of inhomogeneous hard-sphere density correlations in real space. These not only provide information about the packing structures of the hard-sphere system but also form an essential building-block for constructing perturbation theories of more realistic models.

Mean-field theory of inhomogeneous fluids

The Barker-Henderson perturbation theory is a bedrock of liquid-state physics, providing quantitative predictions for the bulk thermodynamic properties of realistic model systems. However, this successful method has not been exploited for the study of inhomogeneous systems. We develop and implement a first-principles "Barker-Henderson density functional," thus providing a robust and quantitatively accurate theory for classical fluids in external fields. Numerical results are presented for the hard-core Yukawa model in three dimensions. Our predictions for the density around a fixed test particle and between planar walls are in very good agreement with simulation data. The density profiles for the free liquid vapor interface show the expected oscillatory decay into the bulk liquid as the temperature is reduced toward the triple point, but with an amplitude much smaller than that predicted by the standard mean-field density functional.

Nondiffusive fluxes in a Brownian system with Lorentz force

The Fokker-Planck equation provides a complete statistical description of a particle undergoing random motion in a solvent. In the presence of Lorentz force due to an external magnetic field, the Fokker-Planck equation picks up a tensorial coefficient, which reflects the anisotropy of the particle's motion. This tensor, however, cannot be interpreted as a diffusion tensor; there are antisymmetric terms which give rise to fluxes perpendicular to the density gradients. Here, we show that for an inhomogeneous magnetic field these nondiffusive fluxes have finite divergence and therefore affect the density evolution of the system. Only in the special cases of a uniform magnetic field or carefully chosen initial condition with the same full rotational symmetry as the magnetic field can these fluxes be ignored in the density evolution.

Escape rate of transiently active Brownian particle in one dimension

Activity significantly enhances the escape rate of a Brownian particle over a potential barrier. Whereas constant activity has been extensively studied in the past, little is known about the effect of time-dependent activity on the escape rate of the particle. In this paper, we study the escape problem for a Brownian particle that is transiently active; the activity decreases rapidly during the escape process. Using the effective equilibrium approach, we analytically calculate the escape rate under the assumption that the particle is either completely passive or fully active when crossing the barrier. We perform numerical simulations of the escape process in one dimension and find good agreement with the theoretical predictions.

Particle-conserving dynamics on the single-particle level

We generalize the particle-conserving dynamics method of de las Heras et al. [J. Phys.: Condens. Matter 28, 244024 (2016)JCOMEL0953-898410.1088/0953-8984/28/24/244024] to binary mixtures and apply this to hard rods in one dimension. Considering the case of one species consisting of only one particle enables us to address the tagged-particle dynamics. The time-evolution of the species-labeled density profiles is compared to exact Brownian dynamics and (grand-canonical) dynamical density functional theory. The particle-conserving dynamics yields improved results over the dynamical density functional theory and well reproduces the simulation data at short and intermediate times. However, the neglect of a strict particle order (due to the fundamental statistical assumption of ergodicity) leads to errors at long times for our one-dimensional setup. The isolated study of that error makes clear the fundamental limitations of (adiabatic) density-based theoretical approaches when applied to systems of any dimension for which particle caging is a dominant physical mechanism.

Isotropic-nematic transition of self-propelled rods in three dimensions

Using overdamped Brownian dynamics simulations we investigate the isotropic-nematic (IN) transition of self-propelled rods in three spatial dimensions. For two well-known model systems (Gay-Berne potential and hard spherocylinders) we find that turning on activity moves to higher densities the phase boundary separating an isotropic phase from a (nonpolar) nematic phase. This active IN phase boundary is distinct from the boundary between isotropic and polar-cluster states previously reported in two-dimensional simulation studies and, unlike the latter, is not sensitive to the system size. We thus identify a generic feature of anisotropic active particles in three dimensions.

Effective equilibrium states in mixtures of active particles driven by colored noise

We consider the steady-state behavior of pairs of active particles having different persistence times and diffusivities. To this purpose we employ the active Ornstein-Uhlenbeck model, where the particles are driven by colored noises with exponential correlation functions whose intensities and correlation times vary from species to species. By extending Fox's theory to many components, we derive by functional calculus an approximate Fokker-Planck equation for the configurational distribution function of the system. After illustrating the predicted distribution in the solvable case of two particles interacting via a harmonic potential, we consider systems of particles repelling through inverse power-law potentials. We compare the analytic predictions to computer simulations for such soft-repulsive interactions in one dimension and show that at linear order in the persistence times the theory is satisfactory. This work provides the toolbox to qualitatively describe many-body phenomena, such as demixing and depletion, by means of effective pair potentials.

Dynamical density functional theory analysis of the laning instability in sheared soft matter

Using dynamical density functional theory (DDFT) methods we investigate the laning instability of a sheared colloidal suspension. The nonequilibrium ordering at the laning transition is driven by nonaffine particle motion arising from interparticle interactions. Starting from a DDFT which incorporates the nonaffine motion, we perform a linear stability analysis that enables identification of the regions of parameter space where lanes form. We illustrate our general approach by applying it to a simple one-component fluid of soft penetrable particles.

Brownian systems with spatially inhomogeneous activity

We generalize the Green-Kubo approach, previously applied to bulk systems of spherically symmetric active particles [J. Chem. Phys. 145, 161101 (2016)JCPSA60021-960610.1063/1.4966153], to include spatially inhomogeneous activity. The method is applied to predict the spatial dependence of the average orientation per particle and the density. The average orientation is given by an integral over the self part of the Van Hove function and a simple Gaussian approximation to this quantity yields an accurate analytical expression. Taking this analytical result as input to a dynamic density functional theory approximates the spatial dependence of the density in good agreement with simulation data. All theoretical predictions are validated using Brownian dynamics simulations.

Escape rate of active particles in the effective equilibrium approach

The escape rate of a Brownian particle over a potential barrier is accurately described by the Kramers theory. A quantitative theory explicitly taking the activity of Brownian particles into account has been lacking due to the inherently out-of-equilibrium nature of these particles. Using an effective equilibrium approach [Farage et al., Phys. Rev. E 91, 042310 (2015)PLEEE81539-375510.1103/PhysRevE.91.042310] we study the escape rate of active particles over a potential barrier and compare our analytical results with data from direct numerical simulation of the colored noise Langevin equation. The effective equilibrium approach generates an effective potential that, when used as input to Kramers rate theory, provides results in excellent agreement with the simulation data.

Phase separation on the sphere: Patchy particles and self-assembly

Motivated by observations of heterogeneous domain structure on the surface of cells, we consider a minimal model to describe the dynamics of phase separation on the surface of a spherical particle. Finite-size effects on the curved particle surface lead to the formation of long-lived, metastable states for which the density is distributed in patches over the particle surface. We study the time evolution and stability of these states as a function of both the particle size and the thermodynamic parameters. Finally, by connecting our findings with studies of patchy particles, we consider the implications for self-assembly in many-particle systems.

Effective interactions in active Brownian suspensions

Active colloids exhibit persistent motion, which can lead to motility-induced phase separation (MIPS). However, there currently exists no microscopic theory to account for this phenomenon. We report a first-principles theory, free of fit parameters, for active spherical colloids, which shows explicitly how an effective many-body interaction potential is generated by activity and how this can rationalize MIPS. For a passively repulsive system the theory predicts phase separation and pair correlations in quantitative agreement with simulation. For an attractive system the theory shows that phase separation becomes suppressed by moderate activity, consistent with recent experiments and simulations, and suggests a mechanism for reentrant cluster formation at high activity.

Normal-stress coefficients and rod climbing in colloidal dispersions

We calculate tractable microscopic expressions for the low-shear normal-stress coefficients of colloidal dispersions. Although restricted to the low rate regime, the presented formulas are valid for all volume fractions below the glass transition and for any interaction potential. Numerical results are presented for a system of colloids interacting via a hard-core attractive Yukawa potential, for which we explore the interplay between attraction strength and volume fraction. We show that the normal-stress coefficients exhibit nontrivial features close to the critical point and at high volume fractions in the vicinity of the reentrant glass transition. Finally, we exploit our formulas to make predictions about rod-climbing effects in attractive colloidal dispersions.

Residual stresses in glasses

The history dependence of glasses formed from flow-melted steady states by a sudden cessation of the shear rate γ[over ˙] is studied in colloidal suspensions, by molecular dynamics simulations and by mode-coupling theory. In an ideal glass, stresses relax only partially, leaving behind a finite persistent residual stress. For intermediate times, relaxation curves scale as a function of γ[over ˙]t, even though no flow is present. The macroscopic stress evolution is connected to a length scale of residual liquefaction displayed by microscopic mean-squared displacements. The theory describes this history dependence of glasses sharing the same thermodynamic state variables but differing static properties.

First-principles constitutive equation for suspension rheology

We provide a detailed derivation of a recently developed first-principles approach to calculating averages in systems of interacting, spherical Brownian particles under time-dependent flow. Although we restrict ourselves to flows which are both homogeneous and incompressible, the time dependence and geometry (e.g., shear and extension) are arbitrary. The approximations formulated within mode-coupling theory are particularly suited to dense colloidal suspensions and capture the slow relaxation arising from particle interactions and the resulting glass transition to an amorphous solid. The delicate interplay between slow structural relaxation and time-dependent external flow in colloidal suspensions thus may be studied within a fully tensorial theory.

Dynamics of localized particles from density functional theory

A fundamental assumption of the dynamical density functional theory (DDFT) of colloidal systems is that a grand-canonical free-energy functional may be employed to generate the thermodynamic driving forces. Using one-dimensional hard rods as a model system, we analyze the validity of this key assumption and show that unphysical self-interactions of the tagged particle density fields, arising from coupling to a particle reservoir, are responsible for the excessively fast relaxation predicted by the theory. Moreover, our findings suggest that even employing a canonical functional would not lead to an improvement for many-particle systems, if only the total density is considered. We present several possible schemes to suppress these effects by incorporating tagged densities. When applied to confined systems, we demonstrate, using a simple example, that DDFT necessarily leads to delocalized tagged particle density distributions, which do not respect the fundamental geometrical constraints apparent in Brownian dynamics simulation data. The implication of these results for possible applications of DDFT to treat the glass transition are discussed.

Nonlinear response of dense colloidal suspensions under oscillatory shear: mode-coupling theory and Fourier transform rheology experiments

Using a combination of theory, experiment, and simulation we investigate the nonlinear response of dense colloidal suspensions to large amplitude oscillatory shear flow. The time-dependent stress response is calculated using a recently developed schematic mode-coupling-type theory describing colloidal suspensions under externally applied flow. For finite strain amplitudes the theory generates a nonlinear response, characterized by significant higher harmonic contributions. An important feature of the theory is the prediction of an ideal glass transition at sufficiently strong coupling, which is accompanied by the discontinuous appearance of a dynamic yield stress. For the oscillatory shear flow under consideration we find that the yield stress plays an important role in determining the nonlinearity of the time-dependent stress response. Our theoretical findings are strongly supported by both large amplitude oscillatory experiments (with Fourier transform rheology analysis) on suspensions of thermosensitive core-shell particles dispersed in water and Brownian dynamics simulations performed on a two-dimensional binary hard-disk mixture. In particular, theory predicts nontrivial values of the exponents governing the final decay of the storage and loss moduli as a function of strain amplitude which are in good agreement with both simulation and experiment. A consistent set of parameters in the presented schematic model achieves to jointly describe linear moduli, nonlinear flow curves, and large amplitude oscillatory spectroscopy.

Nonlinear rheology of colloidal dispersions

Colloidal dispersions are commonly encountered in everyday life and represent an important class of complex fluid. Of particular significance for many commercial products and industrial processes is the ability to control and manipulate the macroscopic flow response of a dispersion by tuning the microscopic interactions between the constituents. An important step towards attaining this goal is the development of robust theoretical methods for predicting from first-principles the rheology and nonequilibrium microstructure of well defined model systems subject to external flow. In this review we give an overview of some promising theoretical approaches and the phenomena they seek to describe, focusing, for simplicity, on systems for which the colloidal particles interact via strongly repulsive, spherically symmetric interactions. In presenting the various theories, we will consider first low volume fraction systems, for which a number of exact results may be derived, before moving on to consider the intermediate and high volume fraction states which present both the most interesting physics and the most demanding technical challenges. In the high volume fraction regime particular emphasis will be given to the rheology of dynamically arrested states.

First-principles constitutive equation for suspension rheology

Using mode-coupling theory, we derive a constitutive equation for the nonlinear rheology of dense colloidal suspensions under arbitrary time-dependent homogeneous flow. Generalizing previous results for simple shear, this allows the full tensorial structure of the theory to be identified. Macroscopic deformation measures, such as the Cauchy-Green tensors, thereby emerge. So does a direct relation between the stress and the distorted microstructure, illuminating the interplay of slow structural relaxation and arbitrary imposed flow. We present flow curves for steady planar and uniaxial elongation and compare these to simple shear. The resulting nonlinear Trouton ratios point to a tensorially nontrivial dynamic yield condition for colloidal glasses.

Dynamic glass transition in two dimensions

The question of the existence of a structural glass transition in two dimensions is studied using mode coupling theory (MCT). We determine the explicit d dependence of the memory functional of mode coupling for one-component systems. Applied to two dimensions we solve the MCT equations numerically for monodisperse hard disks. A dynamic glass transition is found at a critical packing fraction phi(c)d=2 approximately equal 0.697 which is above phi(c)d=3 approximately equal 0.516 by about 35%. Phi(c)d scales approximately with phi(rcp)d, the value for random close packing, at least for d=2, 3. Quantities characterizing the local, cooperative "cage motion" do not differ much for d=2 and d=3, and we, e.g., find the Lindemann criterion for the localization length at the glass transition. The final relaxation obeys the superposition principle, collapsing remarkably well onto a Kohlrausch law. The d=2 MCT results are in qualitative agreement with existing results from Monte Carlo and molecular dynamics simulations. The mean-squared displacements measured experimentally for a quasi-two-dimensional binary system of dipolar hard spheres can be described satisfactorily by MCT for monodisperse hard disks over four decades in time provided the experimental control parameter Gamma (which measures the strength of dipolar interactions) and the packing fraction phi are properly related to each other.

Dense colloidal suspensions under time-dependent shear

We consider the nonlinear rheology of dense colloidal suspensions under a time-dependent simple shear flow. Starting from the Smoluchowski equation for interacting Brownian particles advected by shearing (ignoring fluctuations in fluid velocity), we develop a formalism which enables the calculation of time-dependent, far-from-equilibrium averages. Taking shear stress as an example, we derive exactly a generalized Green-Kubo relation and an equation of motion for the transient density correlator, involving a three-time memory function. Mode coupling approximations give a closed constitutive equation yielding the time-dependent stress for arbitrary shear rate history. We solve this equation numerically for the special case of a hard sphere glass subject to step strain.

Multiple views of the response of an ensemble of spectro-temporal features support concurrent classification of utterance, prosody, sex and speaker identity

Models of auditory processing, particularly of speech, face many difficulties. These difficulties include variability among speakers, variability in speech rate and robustness to moderate distortions such as time compression. In contrast to the 'invariance of percept' (across different speakers, of different sexes, using different intonation, and so on) is the observation that we are sensitive to the identity, sex and intonation of the speaker. In previous work we have reported that a model based on ensembles of spectro-temporal feature detectors, derived from onset sensitive pre-processing of a limited class of stimuli, preserves significant information about the stimulus class. We have also shown that this is robust with respect to the exact choice of feature set, moderate time compression in the stimulus and speaker variation. Here we extend these results to show a) that by using a classifier based on a network of spiking neurons with spike-driven plasticity, the output of the ensemble constitutes an effective rate coding representation of complex sounds; and b) that the same set of spectro-temporal features concurrently preserve information about a range of qualitatively different classes into which the stimulus might fall. We show that it is possible for multiple views of the same pattern of responses to generate different percepts. This is consistent with suggestions that multiple parallel processes exist within the auditory 'what' pathway with attentional modulation enhancing the task-relevant classification type. We also show that the responses of the ensemble are sparse in the sense that a small number of features respond for each stimulus type. This has implications for the ensembles' ability to generalise, and to respond differentially to a wide variety of stimulus classes.

Inhomogeneous model colloid-polymer mixtures: adsorption at a hard wall

We study the equilibrium properties of inhomogeneous model colloid-polymer mixtures. By integrating out the degrees of freedom of the ideal polymer coils, we derive a formal expression for the effective one-component Hamiltonian of the (hard sphere) colloids that is valid for arbitrary external potentials acting on both the colloids and the polymers. We show how one can recover information about the distribution of polymer in the mixture given knowledge of the colloid correlation functions calculated using the effective one-component Hamiltonian. This result is then used to furnish the connection between the free-volume and perturbation theory approaches to determining the bulk phase equilibria. For the special case of a planar hard wall the effective Hamiltonian takes an explicit form, consisting of zero-, one-, and two-body, but no higher-body, contributions provided the size ratio q=sigma(p)/sigma(c)<0.1547, where sigma(c) and sigma(p) denote the diameters of colloid and polymer respectively. We employ a simple density functional theory to calculate colloid density profiles from this effective Hamiltonian for q=0.1. The resulting profiles are found to agree well with those from Monte Carlo simulations for the same Hamiltonian. Adding very small amounts of polymer gives rise to strong depletion effects at the hard wall which lead to pronounced enhancement of the colloid density profile (close to the wall) over what is found for hard spheres at a hard wall.

Density functional for a model colloid-polymer mixture

We present a density functional theory for mixtures of (hard sphere) colloidal particles and ideal polymers. For this extreme nonadditive system we employ a fundamental measures approach to construct a functional which incorporates the correct dimensional crossover and the exact low density limit. In bulk fluid mixtures the functional yields the same free energy and, therefore, the same gas-liquid (demixing) transition as given by free-volume theory. It generates consistent pair correlation functions; the partial structure factors S(ij)(k) diverge, as k-->0, at the critical point obtained from the free energy. Our results for the structure agree well with those from simulation and Percus-Yevick theory.