Recent large deviation results have provided general lower bounds for the fluctuations of time-integrated currents in the steady state of stochastic systems. A corollary are so-called thermodynamic uncertainty relations connecting precision of estimation to average dissipation. Here we consider this problem but for counting observables, i.e., trajectory observables which, in contrast to currents, are non-negative and nondecreasing in time (and possibly symmetric under time reversal). In the steady state, their fluctuations to all orders are bound from below by a Conway-Maxwell-Poisson distribution dependent only on the averages of the observable and of the dynamical activity. We show how to obtain the corresponding bounds for first-passage times (times when a certain value of the counting variable is first reached) and their uncertainty relations. Just like entropy production does for currents, dynamical activity controls the bounds on fluctuations of counting observables.

We consider diffusion of a particle in rearranging environment, so that the diffusivity of the particle is a stochastic function of time. In our previous model of “diffusing diffusivity” [Jain and Sebastian, J. Phys. Chem. B 120, 3988 (2016)], it was shown that the mean square displacement of particle remains Fickian, i.e., $〈x^2\left(T\right)〉\propto T$ at all times, but the probability distribution of particle displacement is not Gaussian at all times. It is exponential at short times and crosses over to become Gaussian only in a large time limit in the case where the distribution of $D$ in that model has a steady state limit which is exponential, i.e., ${\pi }_e\left(D\right)\sim e^{-D/D_0}$. In the present study, we model the diffusivity of a particle as a Lévy flight process so that $D$ has a power-law tailed distribution, viz., ${\pi }_e\left(D\right)\sim D^{-1-\alpha }$ with $0<\alpha <1$. We find that in the short time limit, the width of displacement distribution is proportional to $\sqrt{T}$, implying that the diffusion is Fickian. But for long times, the width is proportional to $T^{1/2\alpha }$ which is a characteristic of anomalous diffusion. The distribution function for the displacement of the particle is found to be a symmetric stable distribution with a stability index $2\alpha$ which preserves its shape at all times.

A two-dimensional concentrated emulsion exhibits spontaneous rapid destabilization through an avalanche of coalescence events which propagate through the assembly stochastically. We propose a deterministic model to explain the average dynamics of the avalanching process. The dynamics of the avalanche phenomenon is studied as a function of a composite parameter, the decay time ratio, which characterizes the ratio of the propensity of coalescence to cease propagation to that of propagation. When this ratio is small, the avalanche grows autocatalytically to destabilize the emulsion. Using a scaling analysis, we unravel the relation between a local characteristic of the system and a global system wide effect. The anisotropic nature of local coalescence results in a system size dependent transition from nonautocatalytic to autocatalytic behavior. By incorporating uncertainty into the parameters in the model, several possible realizations of the coalescence avalanche are generated. The results are compared with the Monte Carlo simulations to derive insights into how the uncertainty propagates in the system.

The dynamics of an initially nonspherical liquid droplet falling in air under the action of gravity is investigated via three-dimensional numerical simulations of the Navier-Stokes and continuity equations in the inertial regime. The surface tension is considered to be high enough so that a droplet does not undergo breakup. Vertically symmetric oscillations which decay with time are observed for low inertia. The amplitude of these oscillations increases for high Gallilei numbers and the shape asymmetry in the vertical direction becomes prominent. The reason for this asymmetry has been attributed to the higher aerodynamic inertia. Moreover, even for large inertia, no path deviations or oscillations are observed.

A general scheme is derived to connect transitions in configuration space with features in the heat capacity. A formulation in terms of occupation probabilities for local minima that define the potential energy landscape provides a quantitative description of how contributions arise from competition between different states. The theory does not rely on a structural interpretation for the local minima, so it is equally applicable to molecular energy landscapes and the landscapes defined by abstract functions. Applications are presented for low-temperature solid-solid transitions in atomic clusters, which involve just a few local minima with different morphologies, and for cluster melting, which is driven by the landscape entropy associated with the more numerous high-energy minima. Analyzing these features in terms of the balance between states with increasing and decreasing occupation probabilities provides a direct interpretation of the underlying transitions. This approach enables us to identify a qualitatively different transition that is caused by a single local minimum associated with an exceptionally large catchment volume in configuration space for a machine learning landscape.

We consider a spring-block model with both dry and viscous frictions, subjected to a periodic driving allowing mechanically stable configurations to be sampled. We show that under strong driving, the scaling of the correlation length with the energy density is incompatible with the prediction of the Edwards statistical approach, which assumes a uniform sampling of mechanically stable configurations. A crossover between the Edwards scaling and nonstandard high-energy scaling is observed at energy scales that depend on the viscous friction coefficient. Generalizing Edwards thermodynamics, we propose a statistical framework, based on a sampling of marginally stable states, that is able to describe the scaling of the correlation length in the highly viscous regime.

Conventional digital computation is rapidly approaching physical limits for speed and energy dissipation. Here we fabricate and test a simple neuromorphic circuit that models neuronal somas, axons, and synapses with superconducting Josephson junctions. The circuit models two mutually coupled excitatory neurons. In some regions of parameter space the neurons are desynchronized. In others, the Josephson neurons synchronize in one of two states, in-phase or antiphase. An experimental alteration of the delay and strength of the connecting synapses can toggle the system back and forth in a phase-flip bifurcation. Firing synchronization states are calculated >70 000 times faster than conventional digital approaches. With their speed and low energy dissipation (${10}^{-17}\mathrm{J}/\mathrm{spike}$), this set of proof-of-concept experiments establishes Josephson junction neurons as a viable approach for improvements in neuronal computation as well as applications in neuromorphic computing.

An ellipsoid, the simplest nonspherical shape, has been extensively used as a model for elongated building blocks for a wide spectrum of molecular, colloidal, and granular systems. Yet the densest packing of congruent hard ellipsoids, which is intimately related to the high-density phase of many condensed matter systems, is still an open problem. We discover an unusual family of dense crystalline packings of self-dual ellipsoids (ratios of the semiaxes $\alpha :\sqrt{\alpha }:1$), containing 24 particles with a quasi-square-triangular (SQ-TR) tiling arrangement in the fundamental cell. The associated packing density $\varphi$ exceeds that of the densest known SM2 crystal [ A. Donev et al., Phys. Rev. Lett. 92, 255506 (2004)] for aspect ratios α in (1.365, 1.5625), attaining a maximal $\varphi \approx 0.75806...$ at $\alpha =93/64$. We show that the SQ-TR phase derived from these dense packings is thermodynamically stable at high densities over the aforementioned α range and report a phase diagram for self-dual ellipsoids. The discovery of the SQ-TR crystal suggests organizing principles for nonspherical particles and self-assembly of colloidal systems.

The migration of a neutrally buoyant droplet in a tube containing another immiscible liquid is investigated numerically in the creeping flow regime. A fully developed velocity profile is imposed at the inlet of the tube. The interface between the two immiscible fluids is captured using a coupled level-set and volume-of-fluid approach. The deformation and breakup dynamics of the droplet are investigated in terms of three dimensionless parameters, namely, the ratio between the radius of the undeformed droplet and the radius of the capillary tube, the viscosity ratio between the dispersed and the continuous phases, and the capillary number that measures the relative importance of the viscous force over the surface tension force. It has been observed that the droplet, while traversing through the tube, either approaches a steady bulletlike shape or develops a prominent reentrant cavity at its rear. Depending on the initial droplet size, there exists a critical capillary number for every flow configuration beyond which the drop fails to maintain a steady shape and breaks into fragments. The deformation and breakup phenomena depend primarily on the droplet size, the viscosity ratio, and the capillary number. Special attention has been given to the case where the drop diameter is comparable with the tube diameter. A thorough computational study has been conducted to find the critical capillary number for a range of droplets of varied sizes suspended in systems having different viscosity ratios.

Observation of complex systems produces time series with specific dynamics at different time scales. The majority of the existing numerical methods for multiscale analysis first decompose the time series into several simpler components and the multiscale structure is given by the properties of their components. We present a numerical method which describes the multiscale structure of arbitrary time series without decomposing them. It is based on the monotony spectrum defined as the variation of the mean amplitude of the monotonic segments with respect to the mean local time scale during successive averagings of the time series, the local time scales being the durations of the monotonic segments. The maxima of the monotony spectrum indicate the time scales which dominate the variations of the time series. We show that the monotony spectrum can correctly analyze a diversity of artificial time series and can discriminate the existence of deterministic variations at large time scales from the random fluctuations. As an application we analyze the multifractal structure of some hydrological time series.

Recently, it was shown that the probability distribution function (PDF) of the free energy of a single continuum directed polymer (DP) in a random potential, equivalently to the height of a growing interface described by the Kardar-Parisi-Zhang (KPZ) equation, converges at large scale to the Tracy-Widom distribution. The latter describes the fluctuations of the largest eigenvalue of a random matrix, drawn from the Gaussian unitary ensemble (GUE), and the result holds for a DP with fixed end points, i.e., for the KPZ equation with droplet initial conditions. A more general conjecture can be put forward, relating the free energies of $N>1$ noncrossing continuum DP in a random potential, to the sum of the $N$th largest eigenvalues of the GUE. Here, using replica methods, we provide an important test of this conjecture by calculating exactly the right tails of both PDFs and showing that they coincide for arbitrary $N$.

Motivated by the demand for efficient quantum devices to engineer energy transport, we analyze some inhomogeneous quantum spin systems, including $XXZ$ chains, with magnetization baths at the ends. With a goal of finding general properties, we study the effects of suitable transformations on the boundary-driven Lindblad master equation associated with the dynamics of the systems. For asymmetric models with target polarization at the edges or twisted $XY$ boundary gradients, we show the properties of the steady state, which establish the features of the energy current irrespective of the system size and the regime of transport. We show the ubiquitous occurrence of energy rectification and, more interestingly, of an unusual phenomenon: in the absence of an external magnetic field, there is a one-way street for the energy current, i.e., the direction of the energy current does not change as we invert the magnetization baths at the boundaries. Given the extensiveness of the procedures, which essentially involve the properties of the Lindblad master equation, our results certainly follow for other interactions and other boundary conditions. Moreover, our results indicate graded spin chains as genuine quantum rectifiers.

The random walk process in a nonhomogeneous medium, characterized by a Lévy stable distribution of jump length, is discussed. The width depends on a position: either before the jump or after that. In the latter case, the density slope is affected by the variable width and the variance may be finite; then all kinds of the anomalous diffusion are predicted. In the former case, only the time characteristics are sensitive to the variable width. The corresponding Langevin equation with different interpretations of the multiplicative noise is discussed. The dependence of the distribution width on position after jump is interpreted in terms of cognitive abilities and related to such problems as migration in a human population and foraging habits of animals.

The classical setting of stochastic resonance is extended to account for parameter variations leading to transitions between a unique stable state, bistability, and multistability regimes, across singularities of various kinds. Analytic expressions for the amplitude and the phase of the response in terms of key parameters are obtained. The conditions for optimal responses are derived in terms of the bifurcation parameter, the driving frequency, and the noise strength.

We investigate the formation of stable ecological networks where many species share the same resource. We show that such a stable ecosystem naturally occurs as a result of extinctions. We obtain an analytical relation for the number of coexisting species, and we find a relation describing how many species that may become extinct as a result of a sharp environmental change. We introduce a special parameter that is a combination of species traits and resource characteristics used in the model formulation. This parameter describes the pressure on the system to converge, by extinctions. When that stress parameter is large, we obtain that the species traits are concentrated at certain values. This stress parameter is thereby a parameter that determines the level of final biodiversity of the system. Moreover, we show that the dynamics of this limit system can be described by simple differential equations.

We have developed a numerical simulation method for predicting the time dependence (wave form) of pressure at any location in the systemic arterial system in humans. The method uses the matlab-Simulink environment. The input data include explicitly the geometry of the arterial tree, treated up to an arbitrary bifurcation level, and the elastic properties of arteries as well as rheological parameters of blood. Thus, the impact of anatomic details of an individual subject can be studied. The method is applied here to reveal the earliest stages of mechanical reaction of the pressure profiles to sudden local blockages (thromboses or embolisms) of selected arteries. The results obtained with a purely passive model provide reference data indispensable for studies of longer-term effects due to neural and humoral mechanisms. The reliability of the results has been checked by comparison of two available sets of anatomic, elastic, and rheological data involving (i) 55 and (ii) 138 arterial segments. The remaining arteries have been replaced with the appropriate resistive elements. Both models are efficient in predicting an overall shift of pressure, whereas the accuracy of the 55-segment model in reproducing the detailed wave forms and stabilization times turns out dependent on the location of the blockage and the observation point.

We theoretically study the dispersion relation for axially propagating electromagnetic waves throughout a one-dimensional helical structure whose pitch and dielectric and magnetic properties are spatial random functions with specific statistical characteristics. In the system of coordinates rotating with the helix, by using a matrix formalism, we write the set of differential equations that governs the expected value of the electromagnetic field amplitudes and we obtain the corresponding dispersion relation. We show that the dispersion relation depends strongly on the noise intensity introduced in the system and the autocorrelation length. When the autocorrelation length increases at fixed fluctuation and when the fluctuation augments at fixed autocorrelation length, the band gap widens and the attenuation coefficient of electromagnetic waves propagating in the random medium gets larger. By virtue of the degeneracy in the imaginary part of the eigenvalues associated with the propagating modes, the random medium acts as a filter for circularly polarized electromagnetic waves, in which only the propagating backward circularly polarized wave can propagate with no attenuation. Our results are valid for any kind of dielectric and magnetic structures which possess a helical-like symmetry such as cholesteric and chiral smectic-C liquid crystals, structurally chiral materials, and stressed cholesteric elastomers.

The electrical responses of a nematic liquid crystal cell are investigated by means of the elastic continuum theory. The nematic medium is considered as a parallel circuit of a resistance and a capacitance and the electric current profile across the sample is determined as a function of the elastic constants. In the reorientation process of the nematic director, the resistance and capacitance of the sample are determined by taking into account the elastic anisotropy. A nonmonotonic profile for the current is observed in which a minimum value of the current may be used to estimate the elastic constants values. This scenario suggests a theoretical method to determine the values of the bulk elastic constants in a single planar aligned cell just by changing the direction of applied electrical field and measuring the resulting electrical current.

The collective behavior of vibrated hexagonal disks confined in a monolayer is investigated experimentally. Due to the broken circular symmetry, hexagons prefer to rotate upon sufficiently strong driving. Due to the formation of liquid bridges, short-ranged cohesive interactions are introduced upon wetting. Consequently, a nonequilibrium stationary state with the rotating disks self-organized in a hexagonal structure arises. The bond length of the hexagonal structure is slightly smaller than the circumdiameter of a hexagon, indicating geometric frustration. This investigation provides an example where the collective behavior of granular matter is tuned by the shape of individual particles.