Here we study standard and higher-order birth-death processes on fully connected networks, within the perspective of large-deviation theory [also referred to as the Wentzel-Kramers-Brillouin (WKB) method in some contexts]. We obtain a general expression for the leading and next-to-leading terms of the stationary probability distribution of the fraction of “active” sites as a function of parameters and network size $N$. We reproduce several results from the literature and, in particular, we derive all the moments of the stationary distribution for the $q$-susceptible-infected-susceptible ($q\text{-SIS}$) model, i.e., a high-order epidemic model requiring $q$ active (“infected”) sites to activate an additional one. We uncover a very rich scenario for the fluctuations of the fraction of active sites, with nontrivial finite-size-scaling properties. In particular, we show that the variance-to-mean ratio diverges at criticality for $[1\le q\le 3]$, with a maximal variability at $q=2$, confirming that complex-contagion processes can exhibit peculiar scaling features including wild variability. Moreover, the leading order in a large-deviation approach does not suffice to describe them: next-to-leading terms are essential to capture the intrinsic singularity at the origin of systems with absorbing states. Some possible extensions of this work are also discussed.

Traditional stochastic modeling of reactive systems limits the domain of applicability of the associated path thermodynamics to systems involving a single elementary reaction at the origin of each observed change in composition. An alternative stochastic modeling has recently been proposed to overcome this limitation. These two ways of modeling reactive systems are in principle incompatible. The question thus arises about choosing the appropriate type of modeling to be used in practical situations. In the absence of sufficiently accurate experimental results, one way to address this issue is through the microscopic simulation of reactive fluids, usually based on hard-sphere dynamics in the Boltzmann limit. In this paper, we show that results obtained through such simulations unambiguously confirm the predictions of traditional stochastic modeling, invalidating a recently proposed alternative.

A type of chaos called laminar chaos was found in singularly perturbed dynamical systems with periodic time-varying delay [Phys. Rev. Lett. 120, 084102 (2018)]. It is characterized by nearly constant laminar phases, which are periodically interrupted by irregular bursts, where the intensity levels of the laminar phases vary chaotically from phase to phase. In this paper, we demonstrate that laminar chaos can also be observed in systems with quasiperiodic delay, where we generalize the concept of conservative and dissipative delays to such systems. It turns out that the durations of the laminar phases vary quasiperiodically and follow the dynamics of a torus map in contrast to the periodic variation observed for periodic delay. Theoretical and numerical results indicate that introducing a quasiperiodic delay modulation into a time-delay system can lead to a giant reduction of the dimension of the chaotic attractors. By varying the mean delay and keeping other parameters fixed, we found that the Kaplan-Yorke dimension is modulated quasiperiodically over several orders of magnitudes, where the dynamics switches quasiperiodically between different types of high- and low-dimensional types of chaos.

Cells are dynamic systems characterized by temporal variations in biophysical properties such as stiffness and contractility. Recent studies show that the recruitment and release of actin filaments into and out of the cell cortex—a network of proteins underneath the cell membrane—leads to cell stiffening prior to division and softening immediately afterward. In three-dimensional (3D) cell collectives, it is unclear whether the stiffness change during division at the single-cell scale controls the spatial structure and dynamics at the multicellular scale. This is an important question to understand because cell stiffness variations impact cell spatial organization and cancer progression. Using a minimal 3D model incorporating cell birth, death, and cell-to-cell elastic and adhesive interactions, we investigate the effect of mechanical heterogeneity—variations in individual cell stiffnesses that make up the cell collective—on tumor spatial organization and cell dynamics. We discover that spatial mechanical heterogeneity characterized by a spheroid core composed of stiffer cells and softer cells in the periphery emerges within dense 3D cell collectives, which may be a general feature of multicellular tumor growth. We show that heightened spatial mechanical heterogeneity enhances single-cell dynamics and volumetric tumor growth driven by fluctuations in cell elasticity. Our results could have important implications in understanding how spatiotemporal variations in single-cell stiffness determine tumor growth and spread.

When a system is in equilibrium, external perturbations yield a time series of nonequilibrium distributions, and recent experimental techniques give access to the nonequilibrium data that may contain critical information. Jinwoo and Tanaka [Sci. Rep. 5, 7832 (2015)] have provided mathematical proof that such a process's nonequilibrium free energy profile over a system's substates has Jarzynski's work as content, which spontaneously dissipates while molecules perform their tasks. Here we numerically verify this fact and give a practical example where we analyze a computer simulation of RNA translocation by a ring-shaped ATPase motor. By interpreting the cyclic process of substrate translocation as a series of quenching, relaxation, and second quenching, the theory gives how the ATPase motor allocates the hydrolysis energy between individual substates until the end of the process. It turns out that the efficiency of RNA translocation is $48\%\text{–}60\%$ for most molecules, but $12\%$ of molecules achieve $80\%\text{–}100\%$ efficiency, which is consistent with the literature. This theory would be a valuable tool for extracting quantitative information about molecular nonequilibrium behavior from experimental observations.

We study a classical model of thermally fluctuating polymers confined to two dimensions, experiencing a grooved periodic potential, and subject to pulling forces both along and transverse to the grooves. The equilibrium polymer conformations are described by a mapping to a quantum system with a non-Hermitian Hamiltonian and with fermionic statistics generated by noncrossing interactions among polymers. Using molecular dynamics simulations and analytical calculations, we identify a localized and a delocalized phase of the polymer conformations, separated by a delocalization transition which corresponds (in the quantum description) to the breakdown of a band insulator when driven by an imaginary vector potential. We calculate the average tilt of the many-body system, at arbitrary shear values and filling density of polymer chains, in terms of the complex-valued non-Hermitian band structure. We find the critical shear value, the localization length, and the critical exponent by which the shear modulus diverges in terms of the branch points (exceptional points) in the band structure at which the bandgap closes. We also investigate the combined effects of non-Hermitian delocalization and localization due to both periodicity and disorder, uncovering preliminary evidence that while disorder favors localization at high values, it encourages delocalization at lower values.

We experimentally investigate the microfluidic flow of oil-in-water nanoemulsions in aqueous sodium dodecyl sulfate (SDS) solutions having different concentrations and injection flow rates. A coaxial microfluidic device is employed to explore the behavior of nanoemulsion threads in these sheathing SDS solutions. Using two high-speed cameras, which simultaneously capture both top and side views, we reveal a variety of flow phenomena, ranging from simple core-annular flow to complex flows, such as gravitational, inertial, and buckling thread flows. By analyzing these complex flows, we develop a methodology that elucidates the relationship of core-annular and gravitational flows at low flow rates. Further, we examine the off-axis displacements and bending of core threads at large flow rates, and we study the buckling dynamics of nanoemulsion threads subjected to osmotic stresses caused by large SDS concentrations in the sheathing fluid.

The single-component multiphase fluids can be described by a single equation of state (EOS), and various EOSs have been employed in the multiphase lattice Boltzmann (LB) method. In this work, we revisit five commonly used EOSs, including the van der Waals EOS, the Redlich-Kwong EOS, the Redlich-Kwong-Soave EOS, the Peng-Robinson EOS, and the Carnahan-Starling EOS. The recent multiphase LB model with self-tuning EOS is employed because of its thermodynamic consistency in a strict sense and clear physical picture at the microscopic level. First, the way to incorporate these multiphase EOSs is proposed. Two scaling factors are introduced to independently adjust the surface tension and interface thickness, and the lattice sound speed is EOS-dependent to ensure the numerical stability. Then, numerical tests are conducted to validate the incorporations of these EOSs and compare their numerical performances. The surface tension and interface thickness are set to the same values for different EOSs in the comparisons. The liquid and gas densities, surface tension, and interface thickness by the LB simulation agree well with the thermodynamic results. The maximum density ratios achieved with different EOSs are at the same level and could be very close to each other when the interface thickness is relatively small. The effects of multiphase EOS, density ratio, and dimensionless relaxation time on the spurious current are discussed in detail. It is interesting to find the van der Waals EOS shows the best numerical performance in reducing the spurious current.

The macroscopic hydrodynamic equations are derived for many-body systems in the local-equilibrium approach, using the Schrödinger picture of quantum mechanics. In this approach, statistical operators are defined in terms of microscopic densities associated with the fundamentally conserved quantities and other slow modes possibly emerging from continuous symmetry breaking, as well as macrofields conjugated to these densities. Functional identities can be deduced, allowing us to identify the reversible and dissipative parts of the mean current densities, to obtain general equations for the time evolution of the conjugate macrofields, and to establish the relationship to projection-operator methods. The entropy production is shown to be nonnegative by applying the Peierls-Bogoliubov inequality to a quantum integral fluctuation theorem. Using the expansion in the gradients of the conjugate macrofields, the transport coefficients are given by Green-Kubo formulas and the entropy production rate can be expressed in terms of quantum Einstein-Helfand formulas, implying its nonnegativity in agreement with the second law of thermodynamics. The results apply to multicomponent fluids and can be extended to condensed matter phases with broken continuous symmetries.

Accelerating controlled thermodynamic processes requires an auxiliary Hamiltonian to steer the system into instantaneous equilibrium states. An extra energy cost is inevitably needed in such finite-time operation. We recently developed a geodesic approach to minimize such energy cost for the shortcut to isothermal process. The auxiliary control typically contains momentum-dependent terms, which are hard to be experimentally implemented due to the requirement of constantly monitoring the speed. In this work, we employ a variational auxiliary control without the momentum-dependent force to approximate the exact control. Following the geometric approach, we obtain the optimal control protocol with variational minimum energy cost. We demonstrate the construction of such protocol via an example of Brownian motion with a controllable harmonic potential.

The nonequilibrium Fokker-Planck dynamics in an arbitrary force field $\vec{f}\left(\vec{x}\right)$ in dimension $N$ is revisited via the correspondence with the non-Hermitian quantum mechanics in a real scalar potential $V\left(\vec{x}\right)$ and in a purely imaginary vector potential $[-i\vec{A}\left(\vec{x}\right)]$ of real amplitude $\vec{A}\left(\vec{x}\right)$. The relevant parameters of irreversibility are then the $\frac{N(N-1)}{2}$ magnetic matrix elements $B_{nm}\left(\vec{x}\right)=-B_{mn}\left(\vec{x}\right)={\partial }_n{A}_m\left(\vec{x}\right)-{\partial }_m{A}_n\left(\vec{x}\right)$, while it is enlightening to explore the corresponding gauge transformations of the vector potential $\vec{A}\left(\vec{x}\right)$. This quantum interpretation is even more fruitful to study the statistics of all the time-additive observables of the stochastic trajectories, since their generating functions correspond to the same quantum problem with additional scalar and/or vector potentials. Our main conclusion is that the analysis of their large deviations properties and the construction of the corresponding Doob conditioned processes can be drastically simplified via the choice of an appropriate gauge for each purpose. This general framework is then applied to the special time-additive observables of Ornstein-Uhlenbeck trajectories in dimension $N$, whose generating functions correspond to quantum propagators involving quadratic scalar potentials and linear vector potentials, i.e., to quantum harmonic oscillators in constant magnetic matrices. As simple illustrative example, we finally focus on the Brownian gyrator in dimension $N=2$ to compute the large deviations properties of the entropy production of its stochastic trajectories and to construct the corresponding conditioned processes having a given value of the entropy production per unit time.

The well known nonlinear fluctuating hydrodynamics theory has grouped diffusions in anharmonic chains into two universality classes: one is the Kardar-Parisi-Zhang (KPZ) class for chains with either asymmetric potential or nonzero static pressure and the other is the Gaussian class for chains with symmetric potential at zero static pressure, such as Fermi-Pasta-Ulam-Tsingou (FPUT)-$\beta$ chains. However, little is known of the nonequilibrium transient diffusion in anharmonic chains. Here, we reveal that the KPZ class is the only universality class for nonequilibrium transient diffusion, manifested as the KPZ scaling of the side peaks of momentum correlation (corresponding to the sound modes correlation), which was completely unexpected in equilibrium FPUT-$\beta$ chains. The underlying mechanism is that the nonequilibrium soliton dynamics cause nonzero transient pressure so that the sound modes satisfy approximately the noisy Burgers equation, in which the collisions of solitons was proved to yield the KPZ dynamic exponent of the soliton dispersion. Therefore, the unexpected KPZ universality class is obtained in the nonequilibrium transient diffusion in FPUT-$\beta$ chains and the corresponding carriers of nonequilibrium transient diffusion are attributed to solitons.

Experiment, theory, and simulation are employed to understand the dispersion of colloidal particles in a periodic array of oscillating harmonic traps generated by optical tweezers. In the presence of trap oscillation, a nonmonotonic and anisotropic dispersion is observed. Surprisingly, the stiffest traps produce the largest dispersion at a critical frequency, and the particles diffuse significantly faster in the direction of oscillation than those undergoing passive Stokes-Einstein-Sutherland diffusion. Theoretical predictions for the effective diffusivity of the particles as a function of trap stiffness and oscillation frequency are developed using generalized Taylor dispersion theory and Brownian dynamics simulations. Both theory and simulation demonstrate excellent agreement with the experiments, and reveal a “slingshot” mechanism that predicts a significant enhancement of colloidal diffusion in dynamic external fields.

We report here the results of comprehensive plasma characterization and diagnostics by analyzing time-resolved absorption spectra of short ultrabroadband (0.1–50 THz) pulses propagated through the test plasma. Spectral analysis of plasma-induced absorption of such THz pulses provides very direct, in situ, high dynamical range, potentially single-shot access to the plasma density, plasma decay time, electron temperature, and ballistic dynamics of the plasma expansion. We have demonstrated a proof-of-principle measurement of plasma created by an intense laser beam. In particular, we showed a reliable measurement of plasma densities from around ${10}^{16}$ to ${10}^{20}{\mathrm{cm}}^{-3}$. Apart from the plasma parameters, this method allowed us to reconstruct peak intensity inside the plasma spot and to observe a very early stage of plasma evolution after its excitation.

A long-standing expectation is that two repulsively coupled oscillators tend to oscillate in opposite directions. It has been difficult to achieve complete synchrony in coupled identical oscillators with purely repulsive coupling. Here, we introduce a general coupling condition based on the linear matrix of dynamical systems for the emergence of the complete synchronization in pure repulsively coupled oscillators. The proposed coupling profiles (coupling matrices) define a bidirectional cross-coupling link that plays the role of indicator for the onset of complete synchrony between identical oscillators. We illustrate the proposed coupling scheme on several paradigmatic two-coupled chaotic oscillators and validate its effectiveness through the linear stability analysis of the synchronous solution based on the master stability function approach. We further demonstrate that the proposed general condition for the selection of coupling profiles to achieve synchronization even works perfectly for a large ensemble of oscillators.

We identify a class of potentials for which the scattering of flat-top solitons and thin-top solitons of the nonlinear Schrödinger equation with dual nonlinearity can be reflectionless. The scattering is characterized by sharp resonances between regimes of full transmission and full quantum reflection. Perturbative expansion in terms of the magnitude of radiation losses leads to the general form of reflectionless potentials. Simulating the scattering of flat-top solitons and thin-top solitons confirms the reflectionless feature of these potentials.

We develop a rigorous theory of external influences on finite discrete dynamical systems, going beyond the perturbation paradigm, in that the external influence need not be a small contribution. Indeed, the covariance condition can be stated as follows: If we evolve the dynamical system for $n$ time steps and then disturb it, then it is the same as first disturbing the system with the same influence and then letting the system evolve for $n$ time steps. Applying the powerful machinery of resource theories, we develop a theory of covariant influences both when there is a purely deterministic evolution and when randomness is involved. Subsequently, we provide necessary and sufficient conditions for the transition between states under deterministic covariant influences and necessary conditions in the presence of stochastic covariant influences, predicting which transitions between states are forbidden. Our approach, for the first time, employs the framework of resource theories, borrowed from quantum information theory, to the study of finite discrete dynamical systems. The laws we articulate unify the behavior of different types of finite discrete dynamical systems, and their mathematical flavor makes them rigorous and checkable.

The nonbacktracking matrix and the related nonbacktracking centrality (NBC) play a crucial role in models of percolation-type processes on networks, such as nonrecurrent epidemics. Here we study the localization of NBC in infinite sparse networks that contain an arbitrary finite subgraph. Assuming the local tree likeness of the enclosing network, and that branches emanating from the finite subgraph do not intersect at finite distances, we show that the largest eigenvalue of the nonbacktracking matrix of the composite network is equal to the highest of the two largest eigenvalues: that of the finite subgraph and of the enclosing network. In the localized state, when the largest eigenvalue of the subgraph is the highest of the two, we derive explicit expressions for the NBCs of nodes in the subgraph and other nodes in the network. In this state, nonbacktracking centrality is concentrated on the subgraph and its immediate neighborhood in the enclosing network. We obtain simple, exact formulas in the case where the enclosing network is uncorrelated. We find that the mean NBC decays exponentially around the finite subgraph, at a rate which is independent of the structure of the enclosing network, contrary to what was found for the localization of the principal eigenvector of the adjacency matrix. Numerical simulations confirm that our results provide good approximations even in moderately sized, loopy, real-world networks.

A prominent spatiotemporal failure mode of frictional systems is self-healing slip pulses, which are propagating solitonic structures that feature a characteristic length. Here, we numerically derive a family of steady state slip pulse solutions along generic and realistic rate-and-state dependent frictional interfaces, separating large deformable bodies in contact. Such nonlinear interfaces feature a nonmonotonic frictional strength as a function of the slip velocity, with a local minimum. The solutions exhibit a diverging length and strongly inertial propagation velocities, when the driving stress approaches the frictional strength characterizing the local minimum from above, and change their character when it is away from it. An approximate scaling theory quantitatively explains these observations. The derived pulse solutions also exhibit significant spatially-extended dissipation in excess of the edge-localized dissipation (the effective fracture energy) and an unconventional edge singularity. The relevance of our findings for available observations is discussed.

The radiation model is a parameter-free model of human mobility that has been applied primarily for short-distance moves, such as commuting. When applied to migration, it underestimates the number of long-range moves, such as between different US states. Here we show that it additionally suffers from a conceptual inconsistency that can have substantial numerical effects on long-distance moves. We propose a modification of the radiation model that introduces a dependence on the angle between any two alternative potential destinations, accounting for the possibility that migrants may have preferences about the approximate direction of their move. We demonstrate that this modification mitigates the conceptual inconsistency and improves the model fit to observational migration data, without introducing any fitting parameters.