We consider the effect of global repulsive coupling on an ensemble of identical excitable elements. An increase of the coupling strength destabilizes the synchronous equilibrium and replaces it with many attracting oscillatory states, created in the transcritical heteroclinic bifurcation. The period of oscillations is inversely proportional to the distance from the critical parameter value. If the elements interact with the global field via the first Fourier harmonics of their phases, the stable equilibrium is in one step replaced by the attracting continuum of periodic motions.

For a general sensory system following an external stochastic signal, we introduce the sensory capacity. This quantity characterizes the performance of a sensor: sensory capacity is maximal if the instantaneous state of the sensor has as much information about a signal as the whole time series of the sensor. We show that adding a memory to the sensor increases the sensory capacity. This increase quantifies the improvement of the sensor with the addition of the memory. Our results are obtained with the framework of stochastic thermodynamics of bipartite systems, which allows for the definition of an efficiency that relates the rate with which the sensor learns about the signal with the energy dissipated by the sensor, which is given by the thermodynamic entropy production. We demonstrate a general trade-off between sensory capacity and efficiency: if the sensory capacity is equal to its maximum 1, then the efficiency must be less than 1/2. As a physical realization of a sensor we consider a two-component cellular network estimating a fluctuating external ligand concentration as signal. This model leads to coupled linear Langevin equations that allow us to obtain explicit analytical results.

We investigate the transition radiation on a periodically deformed interface between two dielectric media. Under the assumption that the dielectric permittivities of the media are close, a formula is derived for the spectral-angular distribution of the radiated energy in the general case of a nonstatic profile function for the separating boundary. In particular, the latter includes the case of surface waves propagating along the boundary. The numerical examples are given for triangular grating and for sinusoidal profile. We show that instead of a single peak in the backward transition radiation on a flat interface, for periodic interface one has a set of peaks. The number and the locations of the peaks depend on the incidence angle of the charge and on the period of the interface. The conditions are specified for their appearance.

The binding of clusters of metal nanoparticles is partly electrostatic. We address difficulties in calculating the electrostatic energy when high charging energies limit the total charge to a single quantum, entailing unequal potentials on the particles. We show that the energy at small separation $h$ has a singular logarithmic dependence on $h$. We derive a general form for this energy in terms of the singular capacitance of two spheres in near contact $c\left(h\right)$, together with nonsingular geometric features of the cluster. Using this form, we determine the energies of various clusters, finding that more compact clusters are more stable. These energies are proposed to be significant for metal-semiconductor binary nanoparticle lattices found experimentally. We sketch how these effects should dictate the relative abundances of metal nanoparticle clusters in nonpolar solvents.

We derive from first principles a one-way radiative transfer equation for the wave intensity resolved over directions (Wigner transform of the wave field) in random media. It is an initial value problem with excitation from a source which emits waves in a preferred, forward direction. The equation is derived in a regime with small random fluctuations of the wave speed but long distances of propagation with respect to the wavelength, so that cumulative scattering is significant. The correlation length of the medium and the scale of the support of the source are slightly larger than the wavelength, and the waves propagate in a wide cone with an opening angle less than ${180}^{\circ }$, so that the backward and evanescent waves are negligible. The scattering regime is a bridge between that of radiative transfer, where the waves propagate in all directions, and the paraxial regime, where the waves propagate in a narrow angular cone. We connect the one-way radiative transfer equation with the equations satisfied by the Wigner transform of the wave field in these regimes.

In this paper we establish a connection between particle trajectories subject to a nonholonomic constraint and light ray trajectories in a variable index of refraction. In particular, we extend the analysis of systems with linear nonholonomic constraints to the dynamics of particles in a potential subject to nonlinear velocity constraints. We contrast the long time behavior of particles subject to a constant kinetic energy constraint (a thermostat) to particles with the constraint of parallel velocities. We show that, while in the former case the velocities of each particle equalize in the limit, in the latter case all the kinetic energies of each particle remain the same.

We investigate nonequilibrium fluctuations of a solid surface governed by the stochastic Mullins-Herring equation with conserved noise. This equation describes surface diffusion of adatoms accompanied by their exchange between the surface and the bulk of the solid, when desorption of adatoms is negligible. Previous works dealt with dynamic scaling behavior of the fluctuating interface. Here we determine the probability that the interface first reaches a large given height at a specified time. We also find the optimal time history of the interface conditional on this nonequilibrium fluctuation. We obtain these results by developing a macroscopic fluctuation theory of surface diffusion.

Granular self-assembly of confined non-Brownian spheres under gravity is studied by molecular dynamics simulations. Starting from a disordered phase, dry or cohesive spheres organize, by vibrational annealing, into body-centered-tetragonal or face-centered-cubic structures, respectively. During the self-assembling process, isothermal and isodense points are observed. The existence of such points indicates that both granular temperature and packing fraction undergo an inversion process that may be in the core of crystal nucleation. Around the isothermal point, a sudden growth of granular clusters having the maximum coordination number takes place, indicating the outcome of a first-order phase transition. We propose a heuristic equation that successfully describes the dynamic evolution of the local packing fraction in terms of the local granular temperature, along the entire crystallization process.

Mechanical instability takes different forms in various ordered and disordered systems and little is known about how thermal fluctuations affect different classes of mechanical instabilities. We develop an analytic theory involving renormalization of rigidity and coherent potential approximation that can be used to understand finite-temperature mechanical stabilities in various disordered systems. We use this theory to study two disordered lattices: a randomly diluted triangular lattice and a randomly braced square lattice. These two lattices belong to two different universality classes as they approach mechanical instability at $T=0$. We show that thermal fluctuations stabilize both lattices. In particular, the triangular lattice displays a critical regime in which the shear modulus scales as $G\sim T^{1/2}$, whereas the square lattice shows $G\sim T^{2/3}$. We discuss generic scaling laws for finite-$T$ mechanical instabilities and relate them to experimental systems.

We revisit the universal behavior of crystalline membranes at and below the crumpling transition, which pertains to the mechanical properties of important soft and hard matter materials, such as the cytoskeleton of red blood cells or graphene. Specifically, we perform large-scale Monte Carlo simulations of a triangulated two-dimensional phantom network which is freely fluctuating in three-dimensional space. We obtain a continuous crumpling transition characterized by critical exponents which we estimate accurately through the use of finite-size techniques. By controlling the scaling corrections, we additionally compute with high accuracy the asymptotic value of the Poisson ratio in the flat phase, thus characterizing the auxetic properties of this class of systems. We obtain agreement with the value which is universally expected for polymerized membranes with a fixed connectivity.

Bento et al. [Phys. Rev. E 91, 022105 (2015)] recently stated that the Tsallis entropy violates the third law of thermodynamics for $0<q<1$ in the subadditive regime. We first show that the division between the regimes $q<1$ and $q>1$ is already inherent in the fundamental incomplete structure of the deformed logarithms and exponentials underlying the Tsallis entropy. Then, we provide the complete deformed functions and show that the Tsallis entropy conforms to the third law of thermodynamics for both superadditive $q<1$ and subadditive $q>1$ regimes. Finally, we remark that the Tsallis entropy does not require the use of an escort-averaging scheme once it is expressed in terms of the complete deformed functions.

Critical dynamics in two-dimension Ising lattices up to $2048\times 2048$ is simulated on field-programmable-gate-array- based computing devices. Linear relaxation times are measured from extremely long Monte Carlo simulations. The longest simulation has $7.1\times {10}^{16}$ spin updates, which would take over 37 years to simulate on a general purpose computer. The linear relaxation time of the Ising lattices is found to follow the dynamic scaling law for correlation lengths as long as 2048. The dynamic exponent $z$ of the system is found to be 2.179(12), which is consistent with previous studies of Ising lattices with shorter correlation lengths. It is also found that Monte Carlo simulations of critical dynamics in Ising lattices larger than $512\times 512$ are very sensitive to the statistical correlations between pseudorandom numbers, making it even more difficult to study such large systems.

From the Horowitz-Esposito stochastic thermodynamical description of information flows in dynamical systems [J. M. Horowitz and M. Esposito, Phys. Rev. X 4, 031015 (2014)], it is known that while the second law of thermodynamics is satisfied by a joint system, the entropic balance for the subsystems is adjusted by a term related to the mutual information exchange rate between the two subsystems. In this article, we present a quantitative discussion of the conceptual link between the Horowitz-Esposito analysis and the Liang-Kleeman work on information transfer between dynamical system components [X. S. Liang and R. Kleeman, Phys. Rev. Lett. 95, 244101 (2005)]. In particular, the entropic balance arguments employed in the two approaches are compared. Notwithstanding all differences between the two formalisms, our work strengthens the Liang-Kleeman heuristic balance reasoning by showing its formal analogy with the recent Horowitz-Esposito thermodynamic balance arguments.

Chaos synchronization may arise in networks of nonlinear units with delayed couplings. We study complete and sublattice synchronization generated by resonance of two large time delays with a specific ratio. As it is known for single-delay networks, the number of synchronized sublattices is determined by the greatest common divisor (GCD) of the network loop lengths. We demonstrate analytically the GCD condition in networks of iterated Bernoulli maps with multiple delay times and complement our analytic results by numerical phase diagrams, providing parameter regions showing complete and sublattice synchronization by resonance for Tent and Bernoulli maps. We compare networks with the same GCD with single and multiple delays, and we investigate the sensitivity of the correlation to a detuning between the delays in a network of coupled Stuart-Landau oscillators. Moreover, the GCD condition also allows detection of time-delay resonances, leading to high correlations in nonsynchronizable networks. Specifically, GCD-induced resonances are observed both in a chaotic asymmetric network and in doubly connected rings of delay-coupled noisy linear oscillators.

We introduce new quantities for exploratory causal inference between bivariate time series. The quantities, called penchants and leanings, are computationally straightforward to apply, follow directly from assumptions of probabilistic causality, do not depend on any assumed models for the time series generating process, and do not rely on any embedding procedures; these features may provide a clearer interpretation of the results than those from existing time series causality tools. The penchant and leaning are computed based on a structured method for computing probabilities.

In many adult tissues, stem cells and differentiated cells are not homogeneously distributed: stem cells are arranged in periodic “niches,” and differentiated cells are constantly produced and migrate out of these niches. In this article, we provide a general theoretical framework to study mixtures of dividing and actively migrating particles, which we apply to biological tissues. We show in particular that the interplay between the stresses arising from active cell migration and stem cell division give rise to robust stem cell patterns. The instability of the tissue leads to spatial patterns which are either steady or oscillating in time. The wavelength of the instability has an order of magnitude consistent with the biological observations. We also discuss the implications of these results for future in vitro and in vivo experiments.

We discuss the plastic behavior of an amorphous matrix reinforced by hard particles. A mesoscopic depinning-like model accounting for Eshelby elastic interactions is implemented. Only the effect of a plastic disorder is considered. Numerical results show a complex size dependence of the effective flow stress of the amorphous composite. In particular, the departure from the mixing law shows opposite trends associated to the competing effects of the matrix and the reinforcing particles, respectively. The reinforcing mechanisms and their effects on localization are discussed. Plastic strain is shown to gradually concentrate on the weakest band of the system. This correlation of the plastic behavior with the material structure is used to design a simple analytical model. The latter nicely captures reinforcement size effects in ${(logN/N)}^{1/2}$, where $N$ is the linear size of the system, observed numerically. Predictions of the effective flow stress accounting for further logarithmic corrections show a very good agreement with numerical results.

We consider viscous fluids in spherical geometry, a lighter fluid supporting a heavier one. Chandrasekhar [Q. J. Mech. Appl. Math. 8, 1 (1955)] analyzed this unstable configuration providing the equations needed to find, numerically, the exact growth rates for the ensuing Rayleigh-Taylor instability. He also derived an analytic but approximate solution. We point out a weakness in his approximate dispersion relation (DR) and offer a somewhat improved one. A third DR, based on transforming a planar DR into a spherical one, suffers no unphysical predictions and compares reasonably well with the exact work of Chandrasekhar and a more recent numerical analysis of the problem [Terrones and Carrara, Phys. Fluids 27, 054105 (2015)].

Equal size air bubbles that are entrapped by a Taylor vortex of the secondary flow in a Couette device, thereby defying buoyancy, slowly form a stable ordered ring with equal separation distances between all neighbors. We present two models of the process dynamics based on force balance on a bubble in the presence of other bubbles positioned on the same streamline in a simple shear flow. The forces taken into account are the viscous resistance, the added mass force, and the inertia-induced repulsing force between two bubbles in a low-Reynolds-number shear flow obtained in Prakash et al. [J. Prakash et al., Phys. Rev. E 87, 043002 (2013)]. The first model of the process assumes that each bubble interacts solely with its nearest neighbors. The second model takes into account pairwise interactions among all the bubbles in the ring. The performed dynamic simulations were compared to the experimental results reported in Prakash et al. [J. Prakash et al., Phys. Rev. E 87, 043002 (2013)] and to the results of quasistationary models (ignoring the added mass effect) suggested in that paper. It is demonstrated that taking into account the effect of added mass, the models describe the major effect of the bubbles’ ordering, provide good estimation of the relaxation time, and also predict nonmonotonic behavior of the separation distance between the bubbles, which exhibit over- and undershooting of equilibrium separations. The latter effects were observed in experiments, but are not predicted by the quasistationary models.