Motility driven glassy dynamics in confluent epithelial monolayers.

As wounds heal, embryos develop, cancer spreads, or asthma progresses, the cellular monolayer undergoes a glass transition between solid-like jammed and fluid-like flowing states. During some of these processes, the cells undergo an epithelial-to-mesenchymal transition (EMT): they acquire in-plane polarity and become motile. Thus, how motility drives the glassy dynamics in epithelial systems is critical for the EMT process. However, no analytical framework that is indispensable for deeper insights exists. Here, we develop such a theory inspired by a well-known glass theory. One crucial result of this work is that the confluency affects the effective persistence time-scale of active force, described by its rotational diffusivity, Deffr. Deffr differs from the bare rotational diffusivity, Dr, of the motile force due to cell shape dynamics, which acts to rectify the force dynamics: Deffr is equal to Dr when Dr is small and saturates when Dr is large. We test the theoretical prediction of Deffr and how it affects the relaxation dynamics in our simulations of the active Vertex model. This novel effect of Deffr is crucial to understanding the new and previously published simulation data of active glassy dynamics in epithelial monolayers.

A comparative study of the correlation between the structure and the dynamics for systems interacting via attractive and repulsive potentials.

We present the study of the structure-dynamics correlation for systems interacting via attractive Lennard-Jones (LJ) and its repulsive counterpart, the Weeks-Chandler-Andersen (WCA) potentials. The structural order parameter (SOP) is related to the microscopic mean-field caging potential. At a particle level, the SOP shows a distribution. Although the two systems have similar pair structures, their average SOP differs. However, this difference alone is insufficient to explain the well known slowing down of the dynamics in the LJ system at low temperatures. The slowing down can be explained in terms of a stronger coupling between the SOP and the dynamics. To understand the origin of this system specific coupling, we study the difference in the microscopic structure between the hard and soft particles. We find that for the LJ system, the structural differences of the hard and soft particles are more significant and have a much stronger temperature dependence compared to the WCA system. Thus, the study suggests that attractive interaction creates more structurally different communities. This broader difference in the structural communities is probably responsible for stronger coupling between the structure and dynamics. Thus, the system specific structure-dynamics correlation, which also leads to a faster slowing down in the dynamics, appears to have a structural origin. A comparison of the predictive power of our SOP with the local energy and two body excess entropy in determining the dynamics shows that in the LJ system, the enthalpy plays a dominant role and in the WCA system, the entropy plays a dominant role, and our SOP can capture both these contributions.

Fluctuation-dissipation relations in the imbalanced Wilson-Cowan model.

The relation between spontaneous and stimulated brain activity is a fundamental question in neuroscience which has received wide attention in experimental studies. Recently, it has been suggested that the evoked response to external stimuli can be predicted from temporal correlations of spontaneous activity. Previous theoretical results, confirmed by the comparison with magnetoencephalography data for human brains, were obtained for the Wilson-Cowan model in the condition of balance of excitation and inhibition, a signature of a healthy brain. Here we extend previous studies to imbalanced conditions by examining a region of parameter space around the balanced fixed point. Analytical results are compared to numerical simulations of Wilson-Cowan networks. We evidence that in imbalanced conditions the functional form of the time correlation and response functions can show several behaviors, exhibiting also an oscillating regime caused by the emergence of complex eigenvalues. The analytical predictions are fully in agreement with numerical simulations, validating the role of cross-correlations in the response function. Furthermore, we identify the leading role of inhibitory neurons in controlling the overall activity of the system, tuning the level of excitability and imbalance.

Scaling of avalanche shape and activity power spectrum in neuronal networks.

Many systems in nature exhibit avalanche dynamics with scale-free features. A general scaling theory has been proposed for critical avalanche profiles in crackling noise, predicting the collapse onto a universal avalanche shape, as well as the scaling behavior of the activity power spectrum as Brown noise. Recently, much attention has been given to the profile of neuronal avalanches, measured in neuronal systems in vitro and in vivo. Although a universal profile was evidenced, confirming the validity of the general scaling theory, the parallel study of the power spectrum scaling under the same conditions was not performed. The puzzling observation is that in the majority of healthy neuronal systems the power spectrum exhibits a behavior close to 1/f, rather than Brown, noise. Here we perform a numerical study of the scaling behavior of the avalanche shape and the power spectrum for a model of integrate and fire neurons with a short-term plasticity parameter able to tune the system to criticality. We confirm that, at criticality, the average avalanche size and the avalanche profile fulfill the general avalanche scaling theory. However, the power spectrum consistently exhibits Brown noise behavior, for both fully excitatory networks and systems with 30% inhibitory networks. Conversely, a behavior closer to 1/f noise is observed in systems slightly off criticality. Results suggest that the power spectrum is a good indicator to determine how close neuronal activity is to criticality.

Identifying structural signature of dynamical heterogeneity via the local softness parameter.

In this work we study the relationship between the softness of a mean-field caging potential and dynamics at the local level. We first describe the local softness, which shows a distribution, thus identifying structural heterogeneity. We show that the lifetime of the softness parameter is connected to the lifetime of the well-known cage structure in supercooled liquids. Finally, our theory predicts that the local softness and the local dynamics is causal below the onset temperature where there is a decoupling between the short and long time dynamics, thus allowing a static description of the cage. With the decrease in temperature, the correlation between structure and dynamics increases. The study shows that at lower temperatures, the structural heterogeneity increases, and since the structure becomes a better predictor of the dynamics, it leads to an increase in the dynamical heterogeneity. We also find that the softness of a hard, immobile region evolves with time and becomes soft and eventually mobile due to the rearrangements in the neighborhood, confirming the well-known facilitation effect.

Thermodynamics and its correlation with dynamics in a mean-field model and pinned systems: A comparative study using two different methods of entropy calculation.

A recent study introduced a novel mean-field model system where each particle over and above the interaction with its regular neighbors interacts with k extra pseudo-neighbors. Here, we present an extensive study of thermodynamics and its correlation with the dynamics of this system. We surprisingly find that the well-known thermodynamic integration (TI) method of calculating the entropy provides unphysical results. It predicts vanishing of the configurational entropy at temperatures close to the onset temperature of the system and negative values of the configurational entropy at lower temperatures. Interestingly, well below the temperature at which the configurational entropy vanishes, both the collective and the single-particle dynamics of the system show complete relaxation. Negative values of the configurational entropy are unphysical, and complete relaxation when the configurational entropy is zero violates the prediction of the random first-order transition theory (RFOT). However, the entropy calculated using the two-phase thermodynamics (2PT) method remains positive at all temperatures for which we can equilibrate the system, and its values are consistent with RFOT predictions. We find that with an increase in k, the difference in the entropy computed using the two methods increases. A similar effect is also observed for a system where a randomly selected fraction of the particles are pinned in their positions in the equilibrated liquid. We show that the difference in entropy calculated via the 2PT and TI methods increases with pinning density.

Microscopic Theory of Softness in Supercooled Liquids.

We introduce a new measure of the structure of a liquid which is the softness of the mean-field potential developed by us earlier. We find that this softness is sensitive to small changes in the structure. Then, we study its correlation with the supercooled liquid dynamics. The study involves a wide range of liquids (fragile, strong, attractive, repulsive, and active) and predicts some universal behaviors such as the softness being linearly proportional to the temperature and inversely proportional to the activation barrier of the dynamics with system dependent proportionality constants. We establish a master equation between the dynamics and the softness parameter and show that, indeed, the dynamics, when scaled by the temperature and system dependent parameters, show a data collapse when plotted against softness. The dynamics of fragile liquids show a strong softness dependence, whereas that of strong liquids show a much weaker softness dependence. We also connect the present study with the earlier studies of softness involving machine learning (ML), thus, providing a theoretical framework for understanding the ML results.

Effective structure of a system with continuous polydispersity.

In a system of N particles, with continuous size polydispersity, there exists an N(N - 1) number of partial structure factors, making it analytically less tractable. A common practice is to treat the system as an effective one component system, which is known to exhibit an artificial softening of the structure. The aim of this study is to describe the system in terms of M pseudospecies such that we can avoid this artificial softening but, at the same time, have a value of M ≪ N. We use potential energy and pair excess entropy to estimate an optimum number of species, M0. We then define the maximum width of polydispersity, Δσ0, that can be treated as a monodisperse system. We show that M0 depends on the degree and type of polydispersity and also on the nature of the interaction potential, whereas Δσ0 weakly depends on the type of polydispersity but shows a stronger dependence on the type of interaction potential. Systems with a softer interaction potential have a higher tolerance with respect to polydispersity. Interestingly, M0 is independent of system size, making this study more relevant for bigger systems. Our study reveals that even 1% polydispersity cannot be treated as an effective monodisperse system. Thus, while studying the role of polydispersity by using the structure of an effective one component system, care must be taken in decoupling the role of polydispersity from that of the artificial softening of the structure.

Continuous time random walk concepts applied to extended mode coupling theory: a study of the Stokes-Einstein breakdown.

In an attempt to extend the mode coupling theory (MCT) to lower temperatures, some years back an Unified theory was proposed which within the MCT framework incorporated the activated dynamics via the random first order transition theory (RFOT). The theory successfully showed that there is hopping induced diffusive dynamics and the modified MCT coupled to the activated motion continues till low temperatures. Here we show that the theory although successful in describing other properties of supercooled liquids is unable to capture the Stokes-Einstein breakdown. We then show using continuous time random work (CTRW) formalism that the Unified theory is equivalent to a CTRW dynamics in presence of two waiting time distributions. It is known from earlier work on CTRW that in such cases the total dynamics is dominated by the fast motion. This explains the failure of the Unified theory in predicting the SE breakdown as both the structural relaxation and the diffusion process are described by the comparatively fast MCT like dynamics. The study also predicts that other forms of extended MCT with Markovian hopping kernel will face a similar issue. We next modify the Unified theory by applying the concept of renewal theory, usually used in CTRW models where the distribution has a long tail. According to this theory the first jump given by the persistent time is slower than the subsequent jumps given by the exchange time. We first show that for systems with two waiting time distributions even when both the distributions are exponential the persistent time is larger than the exchange time. We also identify the persistent time with the slower activated process. The extended Unified theory can now explain the SE breakdown. In this extended theory at low temperatures the structural relaxation is described by the activated dynamics whereas the diffusion is primarily determined by the MCT like dynamics leading to a decoupling between them. We also calculate a dynamic lengthscale from the wavenumber dependence of the relaxation time. We find that this dynamic length scale grows faster than the static length scale.

Analysis of the anomalous mean-field like properties of Gaussian core model in terms of entropy.

Studies of the Gaussian core model (GCM) have shown that it behaves like a mean-field model and the properties are quite different from standard glass former. In this work, we investigate the entropies, namely, the excess entropy (Sex) and the configurational entropy (Sc) and their different components to address these anomalies. Our study corroborates most of the earlier observations and also sheds new light on the high and low temperature dynamics. We find that unlike in standard glass former where high temperature dynamics is dominated by two-body correlation and low temperature by many-body correlations, in the GCM both high and low temperature dynamics are dominated by many-body correlations. We also find that the many-body entropy which is usually positive at low temperatures and is associated with activated dynamics is negative in the GCM suggesting suppression of activation. Interestingly despite the suppression of activation, the Adam-Gibbs (AG) relation that describes activated dynamics holds in the GCM, thus suggesting a non-activated contribution in AG relation. We also find an overlap between the AG relation and mode coupling power law regime leading to a power law behavior of Sc. From our analysis of this power law behavior, we predict that in the GCM the high temperature dynamics will disappear at dynamical transition temperature and below that there will be a transition to the activated regime. Our study further reveals that the activated regime in the GCM is quite narrow.

Determination of onset temperature from the entropy for fragile to strong liquids.

In this paper, we establish a connection between the onset temperature of glassy dynamics with the change in the entropy for a wide range of model systems. We identify the crossing temperature of pair and excess entropies as the onset temperature. Below the onset temperature, the residual multiparticle entropy, the difference between excess and pair entropies, becomes positive. The positive entropy can be viewed as equivalent to the larger phase space exploration of the system. The new method of onset temperature prediction from entropy is less ambiguous, as it does not depend on any fitting parameter like the existing methods.

Role of the Pair Correlation Function in the Dynamical Transition Predicted by Mode Coupling Theory.

In a recent study, we have found that for a large number of systems the configurational entropy at the pair level S_{c2}, which is primarily determined by the pair correlation function, vanishes at the dynamical transition temperature T_{c}. Thus, it appears that the information of the transition temperature is embedded in the structure of the liquid. In order to investigate this, we describe the dynamics of the system at the mean field level and, using the concepts of the dynamical density functional theory, show that the dynamical transition temperature depends only on the pair correlation function. Thus, this theory is similar in spirit to the microscopic mode coupling theory (MCT). However, unlike microscopic MCT, which predicts a very high transition temperature, the present theory predicts a transition temperature that is similar to T_{c}. This implies that the information of the dynamical transition temperature is embedded in the pair correlation function.

Effect of total and pair configurational entropy in determining dynamics of supercooled liquids over a range of densities.

In this paper, we present a study of supercooled liquids interacting with the Lennard Jones potential and the corresponding purely repulsive (Weeks-Chandler-Andersen) potential, over a range of densities and temperatures, in order to understand the origin of their different dynamics in spite of their structures being similar. Using the configurational entropy as the thermodynamic marker via the Adam Gibbs relation, we show that the difference in the dynamics of these two systems at low temperatures can be explained from thermodynamics. At higher densities both the thermodynamical and dynamical difference between these model systems decrease, which is quantitatively demonstrated in this paper by calculating different parameters. The study also reveals the origin of the difference in pair entropy despite the similarity in the structure. Although the maximum difference in structure is obtained in the partial radial distribution function of the B type of particles, the rdf of AA pairs and AB pairs gives rise to the differences in the entropy and dynamics. This work supports the observation made in an earlier study [A. Banerjee et al., Phys. Rev. Lett. 113, 225701 (2014)] and shows that they are generic in nature, independent of density.

Unraveling the success and failure of mode coupling theory from consideration of entropy.

We analyze the dynamics of model supercooled liquids in a temperature regime where predictions of mode coupling theory (MCT) are known to be valid qualitatively. In this regime, the Adam-Gibbs (AG) relation, based on an activation picture of dynamics, also describes the dynamics satisfactorily, and we explore the mutual consistency and interrelation of these descriptions. Although entropy and dynamics are related via phenomenological theories, the connection between MCT and entropy has not been argued for. In this work, we explore this connection and provide a microscopic derivation of the phenomenological Rosenfeld theory. At low temperatures, the overlap between the MCT power law regime and AG relation implies that the AG relation predicts an avoided divergence at Tc, the origin of which can be related to the vanishing of pair configurational entropy, which we find occurring at the same temperature. We also show that the residual multiparticle entropy plays an important role in describing the relaxation time.

Non-monotonic size dependence of diffusion and levitation effect: a mode-coupling theory analysis.

We present a study of diffusion of small tagged particles in a solvent, using mode coupling theory (MCT) analysis and computer simulations. The study is carried out for various interaction potentials. For the first time, using MCT, it is shown that only for strongly attractive interaction potential with allowing interpenetration between the solute-solvent pair the diffusion exhibits a non-monotonic solute size dependence which has earlier been reported in simulation studies [P. K. Ghorai and S. Yashonath, J. Phys. Chem. B 109, 5824-5835 (2005)]. For weak attractive and repulsive potential the solute size dependence of diffusion shows monotonic behaviour. It is also found that for systems where the interaction potential does not allow solute-solvent interpenetration, the solute cannot explore the neck of the solvent cage. Thus these systems even with strong attractive interaction will never show any non-monotonic size dependence of diffusion. This non-monotonic size dependence of diffusion has earlier been connected to levitation effect [S. Yashonath and P. Santikary, J. Phys. Chem. 98, 6368 (1994)]. We also show that although levitation is a dynamic phenomena, the effect of levitation can be obtained in the static radial distribution function.