Restricted Boltzmann machines are described by the Gibbs measure of a bipartite spin glass, which in turn can be seen as a generalized Hopfield network. This equivalence allows us to characterize the state of these systems in terms of their retrieval capabilities, both at low and high load, of pure states. We study the paramagnetic-spin glass and the spin glass-retrieval phase transitions, as the pattern (i.e., weight) distribution and spin (i.e., unit) priors vary smoothly from Gaussian real variables to Boolean discrete variables. Our analysis shows that the presence of a retrieval phase is robust and not peculiar to the standard Hopfield model with Boolean patterns. The retrieval region becomes larger when the pattern entries and retrieval units get more peaked and, conversely, when the hidden units acquire a broader prior and therefore have a stronger response to high fields. Moreover, at low load retrieval always exists below some critical temperature, for every pattern distribution ranging from the Boolean to the Gaussian case.

We use Langevin dynamics simulations to study dynamical behavior of a dense planar layer of active semiflexible filaments. Using the strength of active force and the thermal persistence length as parameters, we map a detailed phase diagram and identify several nonequilibrium phases in this system. In addition to a slowly flowing melt phase, we observe that, for sufficiently high activity, collective flow accompanied by signatures of local polar and nematic order appears in the system. This state is also characterized by strong density fluctuations. Furthermore, we identify an activity-driven crossover from this state of coherently flowing bundles of filaments to a phase with no global flow, formed by individual filaments coiled into rotating spirals. This suggests a mechanism where the system responds to activity by changing the shape of active agents, an effect with no analog in systems of active particles without internal degrees of freedom.

We use computer simulations to study the behavior of a mixture of large passive charged colloids in a suspension of smaller active dipolar Janus particles. We find that when a single charged colloid is present in solution, it acquires a rotational or translational motion depending on how the active dipoles self-assemble on its surface to form active complexes. The collective behavior of these complexes is quite remarkable, and includes swarming behavior and coherent macroscopic motion. We detail how the variety of different phenomenologies emerging in this system can ultimately be controlled by the strength of the active forces and the relative concentration of the two species.

We study the quasistatic penetration of a flexible beam into a two-dimensional dense granular medium lying on a horizontal plate. Rather than a buckling-like behavior we observe a transition between a regime of crack-like penetration in which the fiber only shows small fluctuations around a stable straight geometry and a bending regime in which the fiber fully bends and advances through series of loading and unloading steps. We show that the shape reconfiguration of the fiber is controlled by a single nondimensional parameter $L/L_c$, which is the ratio of the length of the flexible beam $L$ to $L_c$, a bending elastogranular length scale that depends on the rigidity of the fiber and on the departure from the jamming packing fraction of the granular medium. We show, moreover, that the dynamics of the bending transition in the course of the penetration experiment is gradual and is accompanied by a symmetry breaking of the granular packing fraction in the vicinity of the fiber. Together with the progressive bending of the fiber, a cavity grows downstream of the fiber and the accumulation of grains upstream of the fiber leads to the development of a jammed cluster of grains. We discuss our experimental results in the framework of a simple model of bending-induced compaction and we show that the rate of the bending transition only depends on the control parameter $L/L_c$.

We consider a collection of self-driven apolar particles on a substrate that organize into an active nematic phase at sufficiently high density or low noise. Using the dynamical renormalization group, we systematically study the two-dimensional fluctuating ordered phase in a coarse-grained hydrodynamic description involving both the nematic director and the conserved density field. In the presence of noise, we show that the system always displays only quasi-long-ranged orientational order beyond a crossover scale. A careful analysis of the nonlinearities permitted by symmetry reveals that activity is dangerously irrelevant over the linearized description, allowing giant number fluctuations to persist although now with strong finite-size effects and a nonuniversal scaling exponent. Nonlinear effects from the active currents lead to power-law correlations in the density field, thereby preventing macroscopic phase separation in the thermodynamic limit.

When a flat elastic strip is compressed along its axis, it is bent in one of two possible directions via spontaneous symmetry breaking, forming a cylindrical arc. This is a phenomenon well known as Euler buckling. When this cylindrical section is pushed in the other direction, the bending direction can suddenly reverse. This instability is called “snap-through buckling” and is one of the elementary shape transitions in a prestressed thin structure. Combining experiments and theory, we study snap-buckling of an elastic strip with one end hinged and the other end clamped. These asymmetric boundary constraints break the intrinsic symmetry of the strip, generating mechanical behaviors, including largely hysteretic but reproducible force responses and switchlike discontinuous shape changes. We establish the set of exact analytical solutions to fully explain all our major experimental and numerical findings. Asymmetric boundary conditions arise naturally in diverse situations when a thin object is in contact with a solid surface at one end. The introduction of asymmetry through boundary conditions yields new insight into complex and programmable functionalities in material and industrial design.

The complex zeros of partition functions were originally investigated by Lee and Yang to explain the behavior of condensing gases. Since then, Lee-Yang zeros have become a powerful tool to describe phase transitions in interacting systems. Today, Lee-Yang zeros are no longer just a theoretical concept; they have been determined in recent experiments. In one approach, the Lee-Yang zeros are extracted from the high cumulants of thermodynamic observables at finite size. Here we employ this method to investigate a phase transition in a molecular zipper. From the energy fluctuations in small zippers, we can predict the temperature at which a phase transition occurs in the thermodynamic limit. Even when the system does not undergo a sharp transition, the Lee-Yang zeros carry important information about the large-deviation statistics and its symmetry properties. Our work suggests an interesting duality between fluctuations in small systems and their phase behavior in the thermodynamic limit. These predictions may be tested in future experiments.

Confining a liquid crystal imposes topological constraints on the orientational order, allowing global control of equilibrium systems by manipulation of anchoring boundary conditions. In this article, we investigate whether a similar strategy allows control of active liquid crystals. We study a hydrodynamic model of an extensile active nematic confined in containers, with different anchoring conditions that impose different net topological charges on the nematic director. We show that the dynamics are controlled by a complex interplay between topological defects in the director and their induced vortical flows. We find three distinct states by varying confinement and the strength of the active stress: A topologically minimal state, a circulating defect state, and a turbulent state. In contrast to equilibrium systems, we find that anchoring conditions are screened by the active flow, preserving system behavior across different topological constraints. This observation identifies a fundamental difference between active and equilibrium materials.

The Totally Asymmetric Exclusion Process (TASEP) is a classical stochastic model for describing the transport of interacting particles, such as ribosomes moving along the messenger ribonucleic acid (mRNA) during translation. Although this model has been widely studied in the past, the extent of collision between particles and the average distance between a particle to its nearest neighbor have not been quantified explicitly. We provide here a theoretical analysis of such quantities via the distribution of isolated particles. In the classical form of the model in which each particle occupies only a single site, we obtain an exact analytic solution using the matrix ansatz. We then employ a refined mean-field approach to extend the analysis to a generalized TASEP with particles of an arbitrary size. Our theoretical study has direct applications in mRNA translation and the interpretation of experimental ribosome profiling data. In particular, our analysis of data from Saccharomyces cerevisiae suggests a potential bias against the detection of nearby ribosomes with a gap distance of less than approximately three codons, which leads to some ambiguity in estimating the initiation rate and protein production flux for a substantial fraction of genes. Despite such ambiguity, however, we demonstrate theoretically that the interference rate associated with collisions can be robustly estimated and show that approximately 1% of the translating ribosomes get obstructed.

The fight against biological aging (bio-aging) is long-standing, with the focus of intense research aimed at maintaining high rates of tissue regeneration to promote health and longevity. Nevertheless, there are overwhelming complexities associated with the quantitative analysis of aging. In this study, we sought to quantify bio-aging based on physical aging, by mapping instances of multicellular regeneration to the relaxation of physical systems. An experiment of delayed wound healing assays was devised to obtain delay-dependent healing data. The experiment confirmed the slowdown of healing events, which fitted dynamical scaling just as relaxation events do in physical aging. The scaling exponent, which describes the aging rate in physics, is here similarly proposed as an indicator of the deterioration rate of tissue-regenerative power. Parallel equation-based and cell-based simulations also revealed that asymmetric cell cycle-regulatory mechanisms under strong growth-inhibitory conditions predominantly control the critical slowdown of healing analogous to physical criticality. By establishing a direct link between physical aging and biological aging, we are able to estimate the aging rate of tissues and to achieve an integrated understanding of bio-aging mechanism which may improve the modulation of regeneration for clinical use.

We study a simple stochastic model of neuronal excitatory and inhibitory interactions. The model is defined on a directed lattice and internodes couplings are modulated by a nonlinear function that mimics the process of synaptic activation. We prove that such a system behaves as a fully tunable amplifier: the endogenous component of noise, stemming from finite size effects, seeds a coherent (exponential) amplification across the chain generating giant oscillations with tunable frequencies, a process that the brain could exploit to enhance, and eventually encode, different signals. On a wider perspective, the characterized amplification process could provide a reliable pacemaking mechanism for biological systems. The device extracts energy from the finite size bath and operates as an out of equilibrium thermal machine, under stationary conditions.

We expand the theory of Hawkes processes to the nonstationary case, in which the mutually exciting point processes receive time-dependent inputs. We derive an analytical expression for the time-dependent correlations, which can be applied to networks with arbitrary connectivity, and inputs with arbitrary statistics. The expression shows how the network correlations are determined by the interplay between the network topology, the transfer functions relating units within the network, and the pattern and statistics of the external inputs. We illustrate the correlation structure using several examples in which neural network dynamics are modeled as a Hawkes process. In particular, we focus on the interplay between internally and externally generated oscillations and their signatures in the spike and rate correlation functions.

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

A colloidal particle placed in an inhomogeneous solution of smaller nonadsorbing polymers will move towards regions of lower polymer concentration in order to reduce the free energy of the interface between the surface of the particle and the solution. This phenomenon is known as diffusiophoresis. Treating the polymer as penetrable hard spheres, as in the Asakura-Oosawa model, a simple analytic expression for the diffusiophoretic drift velocity can be obtained. In the context of drying films we show that diffusiophoresis by this mechanism can lead to stratification under easily accessible experimental conditions. By stratification we mean spontaneous formation of a layer of polymer on top of a layer of the colloid. Transposed to the case of binary colloidal mixtures, this offers an explanation for the stratification observed recently in these systems [A. Fortini et al., Phys. Rev. Lett. 116, 118301 (2016)]. Our results emphasize the importance of treating solvent dynamics explicitly in these problems and caution against the neglect of hydrodynamic interactions or the use of implicit solvent models in which the absence of solvent backflow results in an unbalanced osmotic force that gives rise to large but unphysical effects.

The Bayesian data analysis framework has been proven to be a systematic and effective method of parameter inference and model selection for stochastic processes. In this work, we introduce an information content model check that may serve as a goodness-of-fit, like the ${\chi }^2$ procedure, to complement conventional Bayesian analysis. We demonstrate this extended Bayesian framework on a system of Langevin equations, where coordinate-dependent mobilities and measurement noise hinder the normal mean-squared displacement approach.

Chemically active Brownian particles with surface catalytic reactions may repel each other due to diffusiophoretic interactions in the reaction and product concentration fields. The system behavior can be described by a “chemical” coupling parameter ${\mathrm{\Gamma }}_c$ that compares the strength of diffusiophoretic repulsion to Brownian motion, and by a mapping to the classical electrostatic one component plasma (OCP) system. When confined to a constant-volume domain, body-centered cubic (bcc) crystals spontaneously form from random initial configurations when the repulsion is strong enough to overcome Brownian motion. Face-centered cubic (fcc) crystals may also be stable. The “melting point” of the “liquid-to-crystal transition” occurs at ${\mathrm{\Gamma }}_c\approx 140$ for both bcc and fcc lattices.

The generation of the plasma current resulting from Bremsstrahlung absorption is considered. It is shown that the electric current is higher than the naive estimates assuming that electrons absorb only the photon momentum and using the Spitzer conductivity would suggest. The current enhancement is in part because electrons get the recoil momentum from the Coulomb field of ions during the absorption and in part because the electromagnetic power is absorbed asymmetrically within the electron velocity distribution space.

The dynamic and static critical behaviors of driven and equilibrium lattice gas models are studied in two spatial dimensions. We show that in the short-time regime immediately following a critical quench, the dynamics of the transverse anisotropic order parameter, its autocorrelation, and Binder cumulant are consistent with the prediction of a Gaussian, i.e., noninteracting, effective theory, both for the nonequilibrium lattice gases and, to some extent, their equilibrium counterpart. Such a superuniversal behavior is observed only at short times after a critical quench, while the various models display their distinct behaviors in the stationary states, described by the corresponding, known universality classes.

The prediction of residue contacts in a protein solely from sequence information is a promising approach to computational structure prediction. Recent developments use statistical or information theoretic methods to extract contact information from a multiple sequence alignment. Despite good results, accuracy is limited due to usage of two-body interactions within a Potts model. In this paper we generalize this approach and propose a Hamiltonian with an additional three-body interaction term. We derive a mean-field approximation for inference of three-body couplings within a Potts model which is fast enough on modern computers. Finally, we show that our model has a higher accuracy in predicting residue contacts in comparison with the plain two-body-interaction model.

The phase transition of finite Yukawa balls (ordered systems of microspheres in a gaseous plasma environment) with less than 100 particles is studied experimentally by means of configurational entropies. We have developed cylindrical two- and three-particle-correlation functions to measure these entropies for multiple cluster sizes over a wide temperature range. The cluster temperature is finely tuned using a stochastic laser heating setup. It is shown that the correlation functions give a detailed insight into the structural properties of the cluster. The derived configurational entropies give a clear indication of the transition temperature from a solid-like to a fluid-like state. Comparing the transition temperatures of different sized clusters it is found that the transition temperature increases with cluster size in general agreement with theoretical predictions.