BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Peter Sollich

Peter Sollich appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
cond-mat.stat-mechcond-mat.softcond-mat.dis-nnMachine LearningMolecular Networkscond-mat.mtrl-sciphysics.data-anQuantitative MethodsBiological PhysicsComputation and LanguageGenomicsSoundComputer Visioneess.ASmath-phmath.MPnlin.CDphysics.app-phphysics.chem-phTissues and Organs

Catalog footprint

What is connected

60works
20topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

60 published item(s)

preprint2022arXiv

Delayed elastic contributions to the viscoelastic response of foams

We show that the slow viscoelastic response of a foam is that of a power-law fluid with a terminal relaxation. Investigations of the foam mechanics in creep and recovery tests reveal that the power-law contribution is fully reversible, indicative of a delayed elastic response. We demonstrate how this contribution fully accounts for the non-Maxwellian features observed in all tests, probing the linear mechanical response function. The associated power-law spectrum is consistent with soft glassy rheology of systems with mechanical noise temperatures just above the glass transition [Fielding et al., J.Rheol. 44, 323(2000)] and originates from a combination of superdiffusive bubble dynamics and stress diffusion, as recently evidenced in simulations of coarsening foam [Hwang et al., Nat.Mater. 15, 1031 (2016)].

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Exponential increase of transition rates in metastable systems driven by non-Gaussian noise

Non-Gaussian noise influences many complex out-of-equilibrium systems on a wide range of scales such as quantum devices, active and living matter, and financial markets. Despite the ubiquitous nature of non-Gaussian noise, its effect on activated transitions between metastable states has so far not been understood in generality, notwithstanding prior work focusing on specific noise types and scaling regimes. Here, we present a unified framework for a general class of non-Gaussian noise, which we take as any finite-intensity noise with independent and stationary increments. Our framework identifies optimal escape paths as minima of a stochastic action, which enables us to derive analytical results for the dominant scaling of the escape rates in the weak-noise regime generalizing the conventional Arrhenius law. We show that non-Gaussian noise always induces a more efficient escape, by reducing the effective potential barrier compared to the Gaussian case with the same noise intensity. Surprisingly, for a broad class of amplitude distributions even noise of infinitesimally small intensity can induce an exponentially larger escape rate. As the underlying reason we identify the appearance of discontinuous minimal action paths, for which escape from the metastable state involves a finite jump. We confirm the existence of such paths by calculating the prefactor of the escape rate, as well as by numerical simulations. Our results highlight fundamental differences in the escape behaviour of systems subject to thermal and non-thermal fluctuations, which can be tuned to optimize switching behaviour in metastable systems.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Localization properties of the sparse Barrat-Mézard trap model

Inspired by works on the Anderson model on sparse graphs, we devise a method to analyze the localization properties of sparse systems that may be solved using cavity theory. We apply this method to study the properties of the eigenvectors of the master operator of the sparse Barrat-Mézard trap model, with an emphasis on the extended phase. As probes for localization, we consider the inverse participation ratio and the correlation volume, both dependent on the distribution of the diagonal elements of the resolvent. Our results reveal a rich and non-trivial behavior of the estimators across the spectrum of relaxation rates and an interplay between entropic and activation mechanisms of relaxation that give rise to localized modes embedded in the bulk of extended states. We characterize this route to localization and find it to be distinct from the paradigmatic Anderson model or standard random matrix systems.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Mean field description of aging linear response in athermal amorphous solids

We study the linear response to strain in a mean field elastoplastic model for athermal amorphous solids, incorporating the power-law mechanical noise spectrum arising from plastic events. In the "jammed" regime of the model, where the plastic activity exhibits a non-trivial slow relaxation referred to as aging, we find that the stress relaxes incompletely to an age-dependent plateau, on a timescale which grows with material age. We determine the scaling behaviour of this aging linear response analytically, finding that key scaling exponents are universal and independent of the noise exponent $μ$. For $μ>1$, we find simple aging, where the stress relaxation timescale scales linearly with the age $t_{\rm w}$ of the material. At $μ=1$, which corresponds to interactions mediated by the physical elastic propagator, we find instead a $t_{\rm w}^{1/2}$ scaling arising from the stretched exponential decay of the plastic activity. We compare these predictions with measurements of the linear response in computer simulations of a model jammed system of repulsive soft athermal particles, during its slow dissipative relaxation towards mechanical equilibrium, and find good agreement with the theory.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Mean field theory of yielding under oscillatory shear

We study a mean field elastoplastic model, embedded within a disordered landscape of local yield barriers, to shed light on the behaviour of athermal amorphous solids subject to oscillatory shear. We show that the model presents a genuine dynamical transition between an elastic and a yielded state, and qualitatively reproduces the dependence on the initial degree of annealing found in particle simulations. For initial conditions prepared below the analytically derived threshold energy, we observe a non-trivial, non-monotonic approach to the yielded state. The timescale diverges as one approaches the yielding point from above, which we identify with the fatigue limit. We finally discuss the connections to brittle yielding under uniform shear.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Towards Robust Waveform-Based Acoustic Models

We study the problem of learning robust acoustic models in adverse environments, characterized by a significant mismatch between training and test conditions. This problem is of paramount importance for the deployment of speech recognition systems that need to perform well in unseen environments. First, we characterize data augmentation theoretically as an instance of vicinal risk minimization, which aims at improving risk estimates during training by replacing the delta functions that define the empirical density over the input space with an approximation of the marginal population density in the vicinity of the training samples. More specifically, we assume that local neighborhoods centered at training samples can be approximated using a mixture of Gaussians, and demonstrate theoretically that this can incorporate robust inductive bias into the learning process. We then specify the individual mixture components implicitly via data augmentation schemes, designed to address common sources of spurious correlations in acoustic models. To avoid potential confounding effects on robustness due to information loss, which has been associated with standard feature extraction techniques (e.g., FBANK and MFCC features), we focus on the waveform-based setting. Our empirical results show that the approach can generalize to unseen noise conditions, with 150% relative improvement in out-of-distribution generalization compared to training using the standard risk minimization principle. Moreover, the results demonstrate competitive performance relative to models learned using a training sample designed to match the acoustic conditions characteristic of test utterances.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Coarse graining of biochemical systems described by discrete stochastic dynamics

Many biological systems can be described by finite Markov models. A general method for simplifying master equations is presented that is based on merging adjacent states. The approach preserves the steady-state probability distribution and all steady-state fluxes except the one between the merged states. Different levels of coarse graining of the underlying microscopic dynamics can be obtained by iteration, with the result being independent of the order in which states are merged. A criterion for the optimal level of coarse graining or resolution of the process is proposed via a tradeoff between the simplicity of the coarse-grained model and the information loss relative to the original model. As a case study, the method is applied to the cycle kinetics of the molecular motor kinesin.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Robust Prediction of Force Chains in Jammed Solids using Graph Neural Networks

Force chains, which are quasi-linear self-organised structures carrying large stresses, are ubiquitous in jammed amorphous materials, such as granular materials, foams, emulsions or even assemblies of cells. Predicting where they will form upon mechanical deformation is crucial in order to describe the physical properties of such materials, but remains an open question. Here we demonstrate that graph neural networks (GNN) can accurately infer the location of these force chains in frictionless materials from the local structure prior to deformation, without receiving any information about the inter-particle forces. Once trained on a prototypical system, the GNN prediction accuracy proves to be robust to changes in packing fraction, mixture composition, amount of deformation, and the form of the interaction potential. The GNN is also scalable, as it can make predictions for systems much larger than those it was trained on. Our results and methodology will be of interest for experimental realizations of granular matter and jammed disordered systems, e.g. in cases where direct visualisation of force chains is not possible or contact forces cannot be measured.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Shear Induced Orientational Ordering in Active Glass

Dense assemblies of self propelled particles, also known as active or living glasses are abundantaround us, covering different length and time scales: from the cytoplasm to tissues, from bacterialbio-films to vehicular traffic jams, from Janus colloids to animal herds. Being structurally disorderedas well as strongly out of equilibrium, these systems show fascinating dynamical and mechanicalproperties. Using extensive molecular dynamics simulation and a number of different dynamicaland mechanical order parameters we differentiate three dynamical steady states in a sheared modelactive glassy system: (a) a disordered phase, (b) a propulsion-induced ordered phase, and (c) ashear-induced ordered phase. We supplement these observations with an analytical theory based onan effective single particle Fokker-Planck description to rationalise the existence of the novel shear-induced orientational ordering behaviour in our model active glassy system that has no explicitaligning interactions,e.g.of Vicsek-type. This ordering phenomenon occurs in the large persistencetime limit and is made possible only by the applied steady shear. Using a Fokker-Planck descriptionwe make testable predictions without any fit parameters for the joint distribution of single particleposition and orientation. These predictions match well with the joint distribution measured fromdirect numerical simulation. Our results are of relevance for experiments exploring the rheologicalresponse of dense active colloids and jammed active granular matter systems.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

The random first-order transition theory of active glass in the high-activity regime

Dense active matter, in the fluid or amorphous-solid form, has generated intense interest as a model for the dynamics inside living cells and multicellular systems. An extension of the random first-order transition theory (RFOT) to include activity was developed, whereby the activity of the individual particles was added to the free energy of the system in the form of the potential energy of an active particle, trapped by a harmonic potential that describes the effective confinement by the surrounding medium. This active-RFOT model was shown to successfully account for the dependence of the structural relaxation time in the active glass, extracted from simulations, as a function of the activity parameters: the magnitude of the active force ($f_0$) and its persistence time ($τ_p$). However, significant deviations were found in the limit of large activity (large $f_0$ and/or $τ_p$). Here we extend the active-RFOT model to high activity using an activity-dependent harmonic confining potential, which we solve self-consistently. The extended model predicts qualitative changes in the high activity regime, which agree with the results of simulations in both three-dimensional and two-dimensional models of active glass.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Systematic model reduction captures the dynamics of extrinsic noise in biochemical subnetworks

We consider the general problem of describing the dynamics of subnetworks of larger biochemical reaction networks, e.g. protein interaction networks involving complex formation and dissociation reactions. We propose the use of model reduction strategies to understand the 'extrinsic' sources of stochasticity arising from the rest of the network. Our approaches are based on subnetwork dynamical equations derived by projection methods and by path integrals. The results provide a principled derivation of the different components of the extrinsic noise that is observed experimentally in cellular biochemical reactions, over and above the intrinsic noise from the stochasticity of biochemical events in the subnetwork. We explore several intermediate approximations to assess systematically the relative importance of different extrinsic noise components, including initial transients, long-time plateaus, temporal correlations, multiplicative noise terms and nonlinear noise propagation. The best approximations achieve excellent accuracy in quantitative tests on a simple protein network and on the epidermal growth factor receptor signalling network.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Wrinkles, folds and ripplocations: unusual deformation structures of confined elastic sheets at non-zero temperatures

We study the deformation of a fluctuating crystalline sheet confined between two flat rigid walls as a simple model for layered solids where bonds among atoms {\it within} the same layer are much stronger than those {\it between} layers. When subjected to sufficiently high loads in an appropriate geometry, these solids deform and fail in unconventional ways. Recent experiments suggest that configurations named {\it ripplocations}, where a layer folds backwards over itself, are involved. These structures are distinct and separated by large free energy barriers from smooth {\it ripples} of the atomic layers that are always present at any non-zero temperature. We use Monte Carlo simulation in combination with an umbrella sampling technique to obtain conditions under which such structures form and study their specific experimental signatures.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

A nucleation theory for yielding of nearly defect-free crystals: understanding rate dependent yield points

Experiments and simulations show that when an initially defect free rigid crystal is subjected to deformation at a constant rate, irreversible plastic flow commences at the so-called {\em yield point}. The yield point is a weak function of the deformation rate, which is usually expressed as a power law with an extremely small non-universal exponent. We re-analyze a representative set of published data on nanometer sized, mostly defect free, Cu, Ni and Au crystals in the light of a recently proposed theory of yielding based on nucleation of stable stress-free regions inside the metastable rigid solid. The single relation derived here, which is {\em not} a power law, explains data covering {\em fifteen} orders of magnitude in time scales.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Glassy dynamics on networks: local spectra and return probabilities

The slow relaxation and aging of glassy systems can be modelled as a Markov process on a simplified rough energy landscape: energy minima where the system tends to get trapped are taken as nodes of a random network, and the dynamics are governed by the transition rates among these. In this work we consider the case of purely activated dynamics, where the transition rates only depend on the depth of the departing trap. The random connectivity and the disorder in the trap depths make it impossible to solve the model analytically, so we base our analysis on the spectrum of eigenvalues $λ$ of the master operator. We compute the local density of states $ρ(λ|τ)$ for traps with a fixed lifetime $τ$ by means of the cavity method. This exhibits a power law behaviour $ρ(λ|τ)\simτ|λ|^T$ in the regime of small relaxation rates $|λ|$, which we rationalize using a simple analytical approximation. In the time domain, we find that the probabilities of return to a starting node have a power law-tail that is determined by the distribution of excursion times $F(t)\sim t^{-(T+1)}$. We show that these results arise only by the combination of finite configuration space connectivity and glassy disorder, and interpret them in a simple physical picture dominated by jumps to deep neighbouring traps.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Amorphization of a crystalline solid by plastic deformation

We demonstrate that plastic deformation in solids is associated with a dynamic transition that is reminiscent to the transition from a superconducting to a mixed phase in Type II superconductors. We report analytic calculations, extensive molecular dynamics and sequential umbrella sampling Monte Carlo simulations of a two dimensional triangular crystalline solid undergoing plastic deformation under strain. The solid consists of particles connected by harmonic springs. Non-affine displacement fluctuations of the solid are enhanced using an external field, causing a rich deformation behaviour. The external field leads to a mixed phase, where defect and stress-free crystallites are surrounded by a network of amorphous regions with large local internal stress --- a "stress Meissner" effect. The transition shows slow ageing dynamics caused by the presence of many competing, non-crystalline free-energy minima. Under uniform uniaxial strain, this amorphization transition gives rise to irreversible plastic deformation with the amorphous regions appearing as band-like structures. Our results may be checked by careful experiments on colloidal crystals using holographic optical tweezers.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Michaelis-Menten dynamics in protein subnetworks

To understand the behaviour of complex systems it is often necessary to use models that describe the dynamics of subnetworks. It has previously been established using projection methods that such subnetwork dynamics generically involves memory of the past, and that the memory functions can be calculated explicitly for biochemical reaction networks made up of unary and binary reactions. However, many established network models involve also Michaelis-Menten kinetics, to describe e.g. enzymatic reactions. We show that the projection approach to subnetwork dynamics can be extended to such networks, thus significantly broadening its range of applicability. To derive the extension we construct a larger network that represents enzymes and enzyme complexes explicitly, obtain the projected equations, and finally take the limit of fast enzyme reactions that gives back Michaelis-Menten kinetics. The crucial point is that this limit can be taken in closed form. The outcome is a simple procedure that allows one to obtain a description of subnetwork dynamics, including memory functions, starting directly from any given network of unary, binary and Michaelis-Menten reactions. Numerical tests show that this closed form enzyme elimination gives a much more accurate description of the subnetwork dynamics than the simpler method that represents enzymes explicitly, and is also more efficient computationally.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Dynamics and instantaneous normal modes in a liquid with density anomalies

We investigate the relation between the dynamical features of a supercooled liquid and those of its potential energy landscape, focusing on a model liquid with density anomalies. We consider, at fixed temperature, pairs of state points with different density but the same diffusion constant, and find that surprisingly they have identical dynamical features at all length and time scales. This is shown by the collapse of their mean square displacements and of their self--intermediate scattering functions at different wavevectors. We then investigate how the features of the energy landscape change with density, and establish that state points with equal diffusion constant have different landscapes. In particular, we find a correlation between the fraction of instantaneous normal modes connecting different energy minima and the diffusion constant, but unlike in other systems these two quantities are not in one--to--one correspondence with each other, showing that additional landscape features must be relevant in determining the diffusion constant.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Effective interactions and large deviations in stochastic processes

We discuss the relationships between large deviations in stochastic systems, and "effective interactions" that induce particular rare events. We focus on the nature of these effective interactions in physical systems with many interacting degrees of freedom, which we illustrate by reviewing several recent studies. We describe the connections between effective interactions, large deviations at "level 2.5", and the theory of optimal control. Finally, we discuss possible physical applications of variational results associated with those theories.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Elastic models of the glass transition applied to a liquid with density anomalies

Elastic models of the glass transition relate the relaxation dynamics and the elastic properties of structural glasses. They are based on the assumption that the relaxation dynamics occurs through activated events in the energy landscape whose energy scale is set by the elasticity of the material. Here we investigate whether such elastic models describe the relaxation dynamics of systems of particles interacting via a purely repulsive harmonic potential, focusing on a volume fraction and temperature range that is characterized by entropy--driven water--like density anomalies. We do find clear correlations between relaxation time and diffusivity on the one hand, and plateau shear modulus and Debye--Waller factor on the other, thus supporting the validity of elastic models of the glass transition. However, we also show that the plateau shear modulus is not related to the features of the underlying energy landscape of the system, at variance with recent results for power--law potentials. This challenges the common potential energy landscape interpretation of elastic models.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Increased signaling entropy in cancer requires the scale-free property of protein interaction networks

One of the key characteristics of cancer cells is an increased phenotypic plasticity, driven by underlying genetic and epigenetic perturbations. However, at a systems-level it is unclear how these perturbations give rise to the observed increased plasticity. Elucidating such systems-level principles is key for an improved understanding of cancer. Recently, it has been shown that signaling entropy, an overall measure of signaling pathway promiscuity, and computable from integrating a sample's gene expression profile with a protein interaction network, correlates with phenotypic plasticity and is increased in cancer compared to normal tissue. Here we develop a computational framework for studying the effects of network perturbations on signaling entropy. We demonstrate that the increased signaling entropy of cancer is driven by two factors: (i) the scale-free (or near scale-free) topology of the interaction network, and (ii) a subtle positive correlation between differential gene expression and node connectivity. Indeed, we show that if protein interaction networks were random graphs, described by Poisson degree distributions, that cancer would generally not exhibit an increased signaling entropy. In summary, this work exposes a deep connection between cancer, signaling entropy and interaction network topology.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Large-scale fluctuations of the largest Lyapunov exponent in diffusive systems

We present a general formalism for computing the largest Lyapunov exponent and its fluctuations in spatially extended systems described by diffusive fluctuating hydrodynamics, thus extending the concepts of dynamical system theory to a broad range of non-equilibrium systems. Our analytical results compare favourably with simulations of a lattice model of heat conduction. We further show how the computation of the Lyapunov exponent for the Symmetric Simple Exclusion Process relates to damage spreading and to a two-species pair annihilation process, for which our formalism yields new finite size results.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Modeling the emergence of polarity patterns for the intercellular transport of auxin in plants

The hormone auxin is actively transported throughout plants via protein machineries including the dedicated transporter known as PIN. The associated transport is ordered with nearby cells driving auxin flux in similar directions. Here we provide a model of both the auxin transport and of the dynamics of cellular polarisation based on flux sensing. Our main findings are: (i) spontaneous intracellular PIN polarisation arises if PIN recycling dynamics are sufficiently non-linear, (ii) there is no need for an auxin concentration gradient, and (iii) ordered multi-cellular patterns of PIN polarisation are favored by molecular noise.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Non-affine fluctuations and the Statistics of Defect Precursors in the Planar Honeycomb Lattice

Certain localised displacement fluctuations in the planar honeycomb lattice may be identified as precursors to topological defects. We show that these fluctuations are among the most pronounced {\em non-affine} distortions of an elemental coarse graining volume of the honeycomb structure at non zero temperatures. We obtain the statistics of these precursor modes in the canonical ensemble, evaluating exactly their single point and two-point spatio-temporal distributions, for a lattice with harmonic nearest neighbour and next near neighbour bonds. As the solid is destabilised by tuning interactions, the precursor fluctuations diverge and correlations become long-lived and long-ranged.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Speech Recognition Front End Without Information Loss

Speech representation and modelling in high-dimensional spaces of acoustic waveforms, or a linear transformation thereof, is investigated with the aim of improving the robustness of automatic speech recognition to additive noise. The motivation behind this approach is twofold: (i) the information in acoustic waveforms that is usually removed in the process of extracting low-dimensional features might aid robust recognition by virtue of structured redundancy analogous to channel coding, (ii) linear feature domains allow for exact noise adaptation, as opposed to representations that involve non-linear processing which makes noise adaptation challenging. Thus, we develop a generative framework for phoneme modelling in high-dimensional linear feature domains, and use it in phoneme classification and recognition tasks. Results show that classification and recognition in this framework perform better than analogous PLP and MFCC classifiers below 18 dB SNR. A combination of the high-dimensional and MFCC features at the likelihood level performs uniformly better than either of the individual representations across all noise levels.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The average number of distinct sites visited by a random walker on random graphs

We study the linear large $n$ behavior of the average number of distinct sites $S(n)$ visited by a random walker after $n$ steps on a large random graph. An expression for the graph topology dependent prefactor $B$ in $S(n) = Bn$ is proposed. We use generating function techniques to relate this prefactor to the graph adjacency matrix and then devise message-passing equations to calculate its value. Numerical simulations are performed to evaluate the agreement between the message passing predictions and random walk simulations on random graphs. Scaling with system size and average graph connectivity are also analysed.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Demixing cascades in cluster crystals

In a cluster crystal, each lattice site is occupied by multiple soft-core particles. As the number density is increased at zero temperature, a `cascade' of isostructural phase transitions can occur between states whose site occupancy differs by unity. For low but finite temperature, each of these transitions terminates in a critical point. Using tailored Monte Carlo simulation techniques we have studied such demixing cascades in systems of soft particles interacting via potentials of the generalized exponential form $u(r)=ε\exp[-(r/σ)^n]$. We have estimated the critical parameters of the first few transitions in the cascade as a function of the softness parameter $n$. The critical temperature and pressure exhibit non-monotonic behaviour as $n$ is varied, although the critical chemical potential remains monotonic. The trends for the pressure and chemical potential are confirmed by cell model calculations at zero temperature. As $n\to 2^+$, all the transitions that we have observed are preempted by melting although we cannot rule out that clustering transitions survive at high density.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Extensive load in multitasking associative networks

We use belief-propagation techniques to study the equilibrium behavior of a bipartite spin-glass, with interactions between two sets of $N$ and $P = αN$ spins. Each spin has a finite degree, i.e.\ number of interaction partners in the opposite set; an equivalent view is then of a system of $N$ neurons storing $P$ diluted patterns. We show that in a large part of the parameter space of noise, dilution and storage load, delimited by a critical surface, the network behaves as an extensive parallel processor, retrieving all $P$ patterns {\it in parallel} without falling into spurious states due to pattern cross-talk and typical of the structural glassiness built into the network. Our approach allows us to consider effects beyond those studied in replica theory so far, including pattern asymmetry and heterogeneous dilution. Parallel extensive retrieval is more robust for homogeneous degree distributions, and is not disrupted by biases in the distributions of the spin-glass links.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Hyperuniformity and phase separation in biased ensembles of trajectories for diffusive systems

We analyse biased ensembles of trajectories for diffusive systems. In trajectories biased either by the total activity or the total current, we use fluctuating hydrodynamics to show that these systems exhibit phase transtions into `hyperuniform' states, where large-wavelength density fluctuations are strongly suppressed. We illustrate this behaviour numerically for a system of hard particles in one dimension and we discuss how it appears in simple exclusion processes. We argue that these diffusive systems generically respond very strongly to any non-zero bias, so that homogeneous states with "normal" fluctuations (finite compressibility) exist only when the bias is very weak.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Memory effects in biochemical networks as the natural counterpart of extrinsic noise

We show that in the generic situation where a biological network, e.g. a protein interaction network, is in fact a subnetwork embedded in a larger "bulk" network, the presence of the bulk causes not just extrinsic noise but also memory effects. This means that the dynamics of the subnetwork will depend not only on its present state, but also its past. We use projection techniques to get explicit expressions for the memory functions that encode such memory effects, for generic protein interaction networks involving binary and unary reactions such as complex formation and phosphorylation, respectively. Remarkably, in the limit of low intrinsic copy-number noise such expressions can be obtained even for nonlinear dependences on the past. We illustrate the method with examples from a protein interaction network around epidermal growth factor receptor (EGFR), which is relevant to cancer signalling. These examples demonstrate that inclusion of memory terms is not only important conceptually but also leads to substantially higher quantitative accuracy in the predicted subnetwork dynamics.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Signalling Entropy: a novel network-theoretical framework for systems analysis and interpretation of functional omic data

A key challenge in systems biology is the elucidation of the underlying principles, or fundamental laws, which determine the cellular phenotype. Understanding how these fundamental principles are altered in diseases like cancer is important for translating basic scientific knowledge into clinical advances. While significant progress is being made, with the identification of novel drug targets and treatments by means of systems biological methods, our fundamental systems level understanding of why certain treatments succeed and others fail is still lacking. We here advocate a novel methodological framework for systems analysis and interpretation of molecular omic data, which is based on statistical mechanical principles. Specifically, we propose the notion of cellular signalling entropy (or uncertainty), as a novel means of analysing and interpreting omic data, and more fundamentally, as a means of elucidating systems-level principles underlying basic biology and disease. We describe the power of signalling entropy to discriminate cells according to differentiation potential and cancer status. We further argue the case for an empirical cellular entropy-robustness correlation theorem and demonstrate its existence in cancer cell line drug sensitivity data. Specifically, we find that high signalling entropy correlates with drug resistance and further describe how entropy could be used to identify the achilles heels of cancer cells. In summary, signalling entropy is a deep and powerful concept, based on rigorous statistical mechanical principles, which, with improved data quality and coverage, will allow a much deeper understanding of the systems biological principles underlying normal and disease physiology.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A Subband-Based SVM Front-End for Robust ASR

This work proposes a novel support vector machine (SVM) based robust automatic speech recognition (ASR) front-end that operates on an ensemble of the subband components of high-dimensional acoustic waveforms. The key issues of selecting the appropriate SVM kernels for classification in frequency subbands and the combination of individual subband classifiers using ensemble methods are addressed. The proposed front-end is compared with state-of-the-art ASR front-ends in terms of robustness to additive noise and linear filtering. Experiments performed on the TIMIT phoneme classification task demonstrate the benefits of the proposed subband based SVM front-end: it outperforms the standard cepstral front-end in the presence of noise and linear filtering for signal-to-noise ratio (SNR) below 12-dB. A combination of the proposed front-end with a conventional front-end such as MFCC yields further improvements over the individual front ends across the full range of noise levels.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Density anomalies and high-order jamming crossovers

Jamming crossovers occur at zero temperature in assemblies of particles interacting via finite range repulsive potentials, when on increasing the density particles make contacts with those of subsequent coordination shells. Density anomalies, including an increased diffusivity upon isothermal compression and a negative thermal expansion coefficient, are the finite temperature signatures of the jamming crossovers. In this manuscript we show that the jamming crossovers are correlated to an increase of the non affine response of the system to density changes, and clarify that jammed systems evolve upon compression through subsequent Eshlby-like plastic instabilities.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Disentangling glass and jamming physics in the rheology of soft materials

The shear rheology of soft particles systems becomes complex at large density because crowding effects may induce a glass transition for Brownian particles, or a jamming transition for non-Brownian systems. Here we successfully explore the hypothesis that the shear stress contributions from glass and jamming physics are `additive'. We show that the experimental flow curves measured in a large variety of soft materials (colloidal hard spheres, microgel suspensions, emulsions, aqueous foams) as well as numerical flow curves obtained for soft repulsive particles in both thermal and athermal limits are well described by a simple model assuming that glass and jamming rheologies contribute linearly to the shear stress, provided that the relevant scales for time and stress are correctly identified in both sectors. Our analysis confirms that the dynamics of colloidal hard spheres is uniquely controlled by glass physics while aqueous foams are only sensitive to jamming effects. We show that for micron-sized emulsions both contributions are needed to successfully account for the flow curves, which reveal distinct signatures of both phenomena. Finally, for two systems of soft microgel particles we show that the flow curves are representative of the glass transition of colloidal systems, and deduce that microgel particles are not well suited to studying the jamming transition experimentally.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Large deviations of the dynamical activity in the East model: analysing structure in biased trajectories

We consider large deviations of the dynamical activity in the East model. We bias this system to larger than average activity and investigate the structure that emerges. To best characterise this structure, we exploit the fact that there are effective interactions that would reproduce the same behaviour in an equilibrium system. We combine numerical results with linear response theory and variational estimates of these effective interactions, giving the first insights into such interactions in a many-body system, across a wide range of biases. The system exhibits a hierarchy of responses to the bias, remaining quasi-equilibrated on short length scales, but deviating far from equilibrium on large length scales. We discuss the connection between this hierarchy and the hierarchical aging behaviour of the system.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Random walk kernels and learning curves for Gaussian process regression on random graphs

We consider learning on graphs, guided by kernels that encode similarity between vertices. Our focus is on random walk kernels, the analogues of squared exponential kernels in Euclidean spaces. We show that on large, locally treelike, graphs these have some counter-intuitive properties, specifically in the limit of large kernel lengthscales. We consider using these kernels as covariance matrices of e.g.\ Gaussian processes (GPs). In this situation one typically scales the prior globally to normalise the average of the prior variance across vertices. We demonstrate that, in contrast to the Euclidean case, this generically leads to significant variation in the prior variance across vertices, which is undesirable from the probabilistic modelling point of view. We suggest the random walk kernel should be normalised locally, so that each vertex has the same prior variance, and analyse the consequences of this by studying learning curves for Gaussian process regression. Numerical calculations as well as novel theoretical predictions for the learning curves using belief propagation make it clear that one obtains distinctly different probabilistic models depending on the choice of normalisation. Our method for predicting the learning curves using belief propagation is significantly more accurate than previous approximations and should become exact in the limit of large random graphs.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A Monte Carlo method for chemical potential determination in single and multiple occupancy crystals

We describe a Monte Carlo scheme which, in a single simulation, yields a measurement of the chemical potential of a crystalline solid. Within the isobaric ensemble, this immediately provides an estimate of the system free energy, with statistical uncertainties that are determined precisely and transparently. An extension to multiple occupancy ("cluster") solids permits the direct determination of the cluster chemical potential and hence the equilibrium conditions. We apply the method to a model exhibiting cluster crystalline phases, where we find evidence for an infinite cascade of critical points terminating coexistence between crystals of differing site occupancies.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

High--order jamming crossovers and density anomalies

We demonstrate the existence of high--order jamming crossovers in systems of particles with repulsive contact interactions, which originate from the collapse of successive coordination shells. At zero temperature, these crossovers induce an anomalous behavior of the bulk modulus, which varies non--monotonically with the density, while at finite temperature they induce density anomalies consisting in an increased diffusivity upon isothermal compression and in a negative thermal expansion coefficient. We rationalize the dependence of these crossovers on the softness of the interaction potential, and relate the jamming crossovers and the anomalous diffusivity through the investigation of the vibrational spectrum.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Learning curves for multi-task Gaussian process regression

We study the average case performance of multi-task Gaussian process (GP) regression as captured in the learning curve, i.e. the average Bayes error for a chosen task versus the total number of examples $n$ for all tasks. For GP covariances that are the product of an input-dependent covariance function and a free-form inter-task covariance matrix, we show that accurate approximations for the learning curve can be obtained for an arbitrary number of tasks $T$. We use these to study the asymptotic learning behaviour for large $n$. Surprisingly, multi-task learning can be asymptotically essentially useless, in the sense that examples from other tasks help only when the degree of inter-task correlation, $ρ$, is near its maximal value $ρ=1$. This effect is most extreme for learning of smooth target functions as described by e.g. squared exponential kernels. We also demonstrate that when learning many tasks, the learning curves separate into an initial phase, where the Bayes error on each task is reduced down to a plateau value by "collective learning" even though most tasks have not seen examples, and a final decay that occurs once the number of examples is proportional to the number of tasks.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Non-affine displacements in crystalline solids in the harmonic limit

A systematic coarse graining of microscopic atomic displacements generates a local elastic deformation tensor ${\mathsf D}$ as well as a positive definite scalar $χ$ measuring non-affinity, i.e. the extent to which the displacements are not representable as affine deformations of a reference crystal. We perform an exact calculation of the statistics of $χ$ and ${\mathsf D}$ and their spatial correlations for solids at low temperatures, within a harmonic approximation and in one and two dimensions. We obtain the joint distribution $P(χ, {\mathsf D})$ and the two point spatial correlation functions for $χ$ and ${\mathsf D}$. We show that non-affine and affine deformations are coupled even in a harmonic solid, with a strength that depends on the size of the coarse graining volume $Ω$ and dimensionality. As a corollary to our work, we identify the field, $h_χ$, conjugate to $χ$ and show that this field may be tuned to produce a transition to a state where the ensemble average, $<χ>$, and the correlation length of $χ$ diverge. Our work should be useful as a template for understanding non-affine displacements in realistic systems with or without disorder and as a means for developing computational tools for studying the effects of non-affine displacements in melting, plastic flow and the glass transition.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Notes on the polynomial identities in random overlap structures

In these notes we review first in some detail the concept of random overlap structure (ROSt) applied to fully connected and diluted spin glasses. We then sketch how to write down the general term of the expansion of the energy part from the Boltzmann ROSt (for the Sherrington-Kirkpatrick model) and the corresponding term from the RaMOSt, which is the diluted extension suitable for the Viana-Bray model. From the ROSt energy term, a set of polynomial identities (often known as Aizenman-Contucci or AC relations) is shown to hold rigorously at every order because of a recursive structure of these polynomials that we prove. We show also, however, that this set is smaller than the full set of AC identities that is already known. Furthermore, when investigating the RaMOSt energy for the diluted counterpart, at higher orders, combinations of such AC identities appear, ultimately suggesting a crucial role for the entropy in generating these constraints in spin glasses.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Rectification of asymmetric surface vibrations with dry friction: an exactly solvable model

We consider a stochastic model for the directed motion of a solid object due to the rectification of asymmetric surface vibrations with Poissonian shot-noise statistics. The friction between the object and the surface is given by a piecewise-linear friction force. This models the combined effect of dynamic friction and singular dry friction. We derive an exact solution of the stationary Kolmogorov-Feller (KF) equation in the case of two-sided exponentially distributed amplitudes. The stationary density of the velocity exhibits singular features such as a discontinuity and a delta-peak singularity at zero velocity, and also contains contributions from non-integrable solutions of the KF equation. The mean velocity in our model generally varies non-monotonically as the strength of the dry friction is increased, indicating that transport improves for increased dissipation.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Replica theory for learning curves for Gaussian processes on random graphs

Statistical physics approaches can be used to derive accurate predictions for the performance of inference methods learning from potentially noisy data, as quantified by the learning curve defined as the average error versus number of training examples. We analyse a challenging problem in the area of non-parametric inference where an effectively infinite number of parameters has to be learned, specifically Gaussian process regression. When the inputs are vertices on a random graph and the outputs noisy function values, we show that replica techniques can be used to obtain exact performance predictions in the limit of large graphs. The covariance of the Gaussian process prior is defined by a random walk kernel, the discrete analogue of squared exponential kernels on continuous spaces. Conventionally this kernel is normalised only globally, so that the prior variance can differ between vertices; as a more principled alternative we consider local normalisation, where the prior variance is uniform.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Singular features in noise-induced transport with dry friction

We present an exactly solvable nonlinear model for the directed motion of an object due to zero-mean fluctuations on a uniform featureless surface. Directed motion results from the effect of dry (Coulombic) friction coupled to asymmetric surface vibrations with Poissonian shot noise statistics. We find that the transport of the object exhibits striking non-monotonic and singular features: transport actually improves for increasing dry friction up to a critical dry friction strength $Δ^*$ and undergoes a transition to a unidirectional mode of motion at $Δ^*$. This transition is indicated by a cusp singularity in the mean velocity of the object. Moreover, the stationary velocity distribution also contains singular features, such as a discontinuity and a delta peak at zero velocity. Our results highlight that dissipation can in fact enhance transport, which might be exploited in artificial small scale systems.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Spectra of Empirical Auto-Covariance Matrices

We compute spectra of sample auto-covariance matrices of second order stationary stochastic processes. We look at a limit in which both the matrix dimension $N$ and the sample size $M$ used to define empirical averages diverge, with their ratio $α=N/M$ kept fixed. We find a remarkable scaling relation which expresses the spectral density $ρ(λ)$ of sample auto-covariance matrices for processes with dynamical correlations as a continuous superposition of appropriately rescaled copies of the spectral density $ρ^{(0)}_α(λ)$ for a sequence of uncorrelated random variables. The rescaling factors are given by the Fourier transform $\hat C(q)$ of the auto-covariance function of the stochastic process. We also obtain a closed-form approximation for the scaling function $ρ^{(0)}_α(λ)$. This depends on the shape parameter $α$, but is otherwise universal: it is independent of the details of the underlying random variables, provided only they have finite variance. Our results are corroborated by numerical simulations using auto-regressive processes.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Spin glass polynomial identities from entropic constraints

The core idea of stochastic stability is that thermodynamic observables must be robust under small (random) perturbations of the quenched Gibbs measure. Combining this idea with the cavity field technique, which aims to measure the free energy increment under addition of a spin to the system, we sketch how to write a stochastic stability approach to diluted mean field spin glasses which explicitly gives overlap constraints as the outcome. We then show that, under minimal mathematical assumptions and for gauge invariant systems (namely those with even Ising interactions), it is possible to "reverse" the idea of stochastic stability and use it to derive a broad class of constraints on the unperturbed quenched Gibbs measure. This paper extends a previous study where we showed how to derive (linear) polynomial identities from the "energy" contribution to the free energy, while here we focus on the consequences of "entropic" constraints. Interestingly, in diluted spin glasses, the entropic approach generates more identities than those found by the energy route or other techniques. The two sets of identities become identical on a fully connected topology, where they reduce to the ones derived by Aizenman and Contucci.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Thermodynamic Interpretation of Soft Glassy Rheology Models

Mesoscopic models play an important role in our understanding of the deformation and flow of amorphous materials. One such description, based on the Shear Transformation Zone (STZ) theory, has recently been re-formulated within a non-equilibrium thermodynamics framework, and found to be consistent with it. We show here that a similar interpretation can be made for the Soft Glassy Rheology (SGR) model. Conceptually this means that the "noise temperature" x, proposed phenomenologically in the SGR model to control the dynamics of a set of slow mesoscopic degrees of freedom, can consistently be interpreted as their actual thermodynamic temperature. (Because such modes are slow to equilibrate, this generally does not coincide with the temperature of the fast degrees of freedom and/or heat bath.) If one chooses to make this interpretation, the thermodynamic framework significantly constrains extensions of the SGR approach to models in which x is a dynamical variable. We assess in this light some such extensions recently proposed in the context of shear banding.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Unified study of glass and jamming rheology in soft particle systems

We explore numerically the shear rheology of soft repulsive particles at large volume fraction. The interplay between viscous dissipation and thermal motion results in multiple rheological regimes encompassing Newtonian, shear-thinning and yield stress regimes near the `colloidal' glass transition when thermal fluctuations are important, crossing over to qualitatively similar regimes near the `jamming' transition when dissipation dominates. In the crossover regime, glass and jamming sectors coexist and give complex flow curves. Although glass and jamming limits are characterized by similar macroscopic flow curves, we show that they occur over distinct time and stress scales and correspond to distinct microscopic dynamics. We propose a simple rheological model describing the glass to jamming crossover in the flow curves, and discuss the experimental implications of our results.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Polydispersity induced solid-solid transitions in model colloids

Specialized Monte Carlo simulation techniques and moment free energy method calculations, capable of treating fractionation exactly, are deployed to study the crystalline phase behaviour of an assembly of spherical particles described by a top-hat "parent" distribution of particle sizes. An increase in either the overall density or the degree of polydispersity is shown to generate a succession of phase transitions in which the system demixes into an ever greater number of face-centred cubic "daughter" phases. Each of these phases is strongly fractionated: it contains a much narrower distribution of particle sizes than is present in the system overall. Certain of the demixing transitions are found to be nearly continuous, accompanied by fluctuations in local particle size correlated over many lattice spacings. We explore possible factors controlling the stability of the phases and the character of the demixing transitions.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Suppressed compressibility at large scale in jammed packings of size disperse spheres

We analyze the large scale structure and fluctuations of jammed packings of size disperse spheres, produced in a granular experiment as well as numerically. While the structure factor of the packings reveals no unusual behavior for small wavevectors, the compressibility displays an anomalous linear dependence at low wavectors and vanishes when q -> 0. We show that such behavior occurs because jammed packings of size disperse spheres have no bulk fluctuations of the volume fraction and are thus hyperuniform, a property not observed experimentally before. Our results apply to arbitrary particle size distributions. For continuous distributions, we derive a perturbative expression for the compressibility that is accurate for polydispersity up to about 30%.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Crystalline phases of polydisperse spheres

We use specialized Monte Carlo simulation methods and moment free energy calculations to provide conclusive evidence that dense polydisperse spheres at equilibrium demix into coexisting fcc phases, with more phases appearing as the spread of diameters increases. We manage to track up to four coexisting phases. Each of these is fractionated: it contains a narrower distribution of particle sizes than is present in the system overall. We also demonstrate that, surprisingly, demixing transitions can be nearly continuous, accompanied by fluctuations in local particle size correlated over many lattice spacings.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Duality symmetries in driven one-dimensional hopping models

We consider some duality relations for models of non-interacting particles hopping on disordered one-dimensional chains. In particular, we discuss symmetries of bulk-driven barrier and trap models, and relations between boundary-driven and equilibrium models with related energy landscapes. We discuss the relationships between these duality relations and similar results for interacting many-body systems.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Kinetically Constrained Models

In this chapter we summarize recent developments in the study of kinetically constrained models (KCMs) as models for glass formers. After recalling the definition of the KCMs which we cover we study the possible occurrence of ergodicity breaking transitions and discuss in some detail how, before any such transition occurs, relaxation timescales depend on the relevant control parameter (density or temperature). Then we turn to the main issue: the prediction of KCMs for dynamical heterogeneities. We focus in particular on multipoint correlation functions and susceptibilities, and decoupling in the transport coefficients. Finally we discuss the recent view of KCMs as being at first order coexistence between an active and an inactive space-time phase.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Large deviations and ensembles of trajectories in stochastic models

We consider ensembles of trajectories associated with large deviations of time-integrated quantities in stochastic models. Motivated by proposals that these ensembles are relevant for physical processes such as shearing and glassy relaxation, we show how they can be generated directly using auxiliary stochastic processes. We illustrate our results using the Glauber-Ising chain, for which biased ensembles of trajectories can exhibit ferromagnetic ordering. We discuss the relation between such biased ensembles and quantum phase transitions.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Phase behaviour of polydisperse spheres: simulation strategies and an application to the freezing transition

The statistical mechanics of phase transitions in dense systems of polydisperse particles presents distinctive challenges to computer simulation and analytical theory alike. The core difficulty, namely dealing correctly with particle size fractionation between coexisting phases, is set out in the context of a critique of previous simulation work on such systems. Specialized Monte Carlo simulation techniques and moment free energy method calculations, capable of treating fractionation exactly, are then described and deployed to study the fluid-solid transition of an assembly of repulsive spherical particles described by a top-hat "parent" distribution of particle sizes. The cloud curve delineating the solid-fluid coexistence region is mapped as a function of the degree of polydispersity $δ$, and the properties of the incipient "shadow" phases are presented. The coexistence region is found to shift to higher densities as $δ$ increases, but does not exhibit the sharp narrowing predicted by many theories and some simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Fluid phase coexistence and critical behaviour from simulations in the restricted Gibbs ensemble

The symmetrical restricted Gibbs ensemble (RGE) is a version of the Gibbs ensemble in which particles are exchanged between two boxes of fixed equal volumes. It has recently come to prominence because -- when combined with specialized algorithms -- it provides for the study of near-coexistence density fluctuations in highly size-asymmetric binary mixtures. Hitherto, however, a detailed framework for extracting accurate estimates of critical point and coexistence curve parameters from RGE density fluctuations has been lacking. Here we address this problem by exploiting an exact link between the RGE density fluctuations and those of the grand canonical ensemble. In the sub-critical region we propose and test a simple method for obtaining accurate estimates of coexistence densities. In the critical region we identify an observable that serves as a finite system size estimator for the critical point parameters, and present a finite-size scaling theory that allows extrapolation to the thermodynamic limit.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Phase Transition in a Random Minima Model: Mean Field Theory and Exact Solution on the Bethe Lattice

We consider the number and distribution of minima in random landscapes defined on non-Euclidean lattices. Using an ensemble where random landscapes are reweighted by a fugacity factor $z$ for each minimum they contain, we construct first a `two-box' mean field theory. This exhibits an ordering phase transition at $z\c=2$ above which one box contains an extensive number of minima. The onset of order is governed by an unusual order parameter exponent $β=1$, motivating us to study the same model on the Bethe lattice. Here we find from an exact solution that for any connectivity $μ+1>2$ there is an ordering transition with a conventional mean field order parameter exponent $β=1/2$, but with the region where this behaviour is observable shrinking in size as $1/μ$ in the mean field limit of large $μ$. We show that the behaviour in the transition region can also be understood directly within a mean field approach, by making the assignment of minima `soft'. Finally we demonstrate, in the simplest mean field case, how the analysis can be generalized to include both maxima and minima. In this case an additional first order phase transition appears, to a landscape in which essentially all sites are either minima or maxima.

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

Nematic-nematic demixing in polydisperse thermotropic liquid crystals

We consider the effects of polydispersity on isotropic-nematic phase equilibria in thermotropic liquid crystals, using a Maier-Saupe theory with factorized interactions. A sufficient spread (approx. 50%) in the interaction strengths of the particles leads to phase separation into two or more nematic phases, which can in addition coexist with an isotropic phase. The isotropic-nematic coexistence region widens dramatically as polydispersity is increased, leading to re-entrant isotropic-nematic phase separation in some regions of the phase diagram. We show that similar phenomena will occur also for non-factorized interactions as long as the interaction strength between any two particle species is lower than the mean of the intra-species interactions.

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

Trap models with slowly decorrelating observables

We study the correlation and response dynamics of trap models of glassy dynamics, considering observables that only partially decorrelate with every jump. This is inspired by recent work on a microscopic realization of such models, which found strikingly simple linear out-of-equilibrium fluctuation-dissipation relations in the limit of slow decorrelation. For the Barrat-Mezard model with its entropic barriers we obtain exact results at zero temperature $T$ for arbitrary decorrelation factor $κ$. These are then extended to nonzero $T$, where the qualitative scaling behaviour and all scaling exponents can still be found analytically. Unexpectedly, the choice of transition rates (Glauber versus Metropolis) affects not just prefactors but also some exponents. In the limit of slow decorrelation even complete scaling functions are accessible in closed form. The results show that slowly decorrelating observables detect persistently slow out-of-equilibrium dynamics, as opposed to intermittent behaviour punctuated by excursions into fast, effectively equilibrated states.

0 citations0 reviewsNo code linkedOpen access
preprint2002arXiv

Isotropic-nematic phase equilibria in the Onsager theory of hard rods with length polydispersity

We analyse the effect of a continuous spread of particle lengths on the phase behavior of rodlike particles, using the Onsager theory of hard rods. Our aim is to establish whether ``unusual'' effects such as isotropic-nematic-nematic (I-N-N) phase separation can occur even for length distributions with a single peak. We focus on the onset of I-N coexistence. For a log-normal distribution we find that a finite upper cutoff on rod lengths is required to make this problem well-posed. The cloud curve, which tracks the density at the onset of I-N coexistence as a function of the width of the length distribution, exhibits a kink; this demonstrates that the phase diagram must contain a three-phase I-N-N region. Theoretical analysis shows that in the limit of large cutoff the cloud point density actually converges to zero, so that phase separation results at any nonzero density; this conclusion applies to all length distributions with fatter-than-exponentail tails. Finally we consider the case of a Schulz distribution, with its exponential tail. Surprisingly, even here the long rods (and hence the cutoff) can dominate the phase behaviour, and a kink in the cloud curve and I-N-N coexistence again result. Theory establishes that there is a nonzero threshold for the width of the length distribution above which these long rod effects occur, and shows that the cloud and shadow curves approach nonzero limits for large cutoff, both in good agreement with the numerical results.

0 citations0 reviewsNo code linkedOpen access
preprint2002arXiv

Simplified Onsager theory for isotropic-nematic phase equilibria of length polydisperse hard rods

Polydispersity is believed to have important effects on the formation of liquid crystal phases in suspensions of rod-like particles. To understand such effects, we analyse the phase behaviour of thin hard rods with length polydispersity. Our treatment is based on a simplified Onsager theory, obtained by truncating the series expansion of the angular dependence of the excluded volume. We describe the model and give the full phase equilibrium equations; these are then solved numerically using the moment free energy method which reduces the problem from one with an infinite number of conserved densities to one with a finite number of effective densities that are moments of the full density distribution. The method yields exactly the onset of nematic ordering. Beyond this, results are approximate but we show that they can be made essentially arbitrarily precise by adding adaptively chosen extra moments, while still avoiding the numerical complications of a direct solution of the full phase equilibrium conditions. We investigate in detail the phase behaviour of systems with three different length distributions: a (unimodal) Schulz distribution, a bidisperse distribution and a bimodal mixture of two Schulz distributions which interpolates between these two cases. A three-phase isotropic-nematic-nematic coexistence region is shown to exist for the bimodal and bidisperse length distributions if the ratio of long and short rod lengths is sufficiently large, but not for the unimodal one. We systematically explore the topology of the phase diagram as a function of the width of the length distribution and of the rod length ratio in the bidisperse and bimodal cases.

0 citations0 reviewsNo code linkedOpen access