The additivity principle (AP) allows one to compute the current distribution in many one-dimensional nonequilibrium systems. Here we extend this conjecture to general $d$-dimensional driven diffusive systems, and validate its predictions against both numerical simulations of rare events and microscopic exact calculations of three paradigmatic models of diffusive transport in $d=2$. Crucially, the existence of a structured current vector field at the fluctuating level, coupled to the local mobility, turns out to be essential to understand current statistics in $d>1$. We prove that, when compared to the straightforward extension of the AP to high $d$, the so-called weak AP always yields a better minimizer of the macroscopic fluctuation theory action for current statistics.

We introduce a model of adaptive temporal networks whose evolution is regulated by an interplay between node activity and dynamic exchange of information through links. We study the model by using a master equation approach. Starting from a homogeneous initial configuration, we show that temporal and structural heterogeneities, characteristic of real-world networks, spontaneously emerge. This theoretically tractable model thus contributes to the understanding of the dynamics of human activity and interaction networks.

The excitation of precursor solitons ahead of a rapidly moving object in a fluid, a spectacular phenomenon in hydrodynamics that has often been observed ahead of moving ships, has surprisingly not been investigated in plasmas where the fluid model holds good for low frequency excitations such as ion acoustic waves. In this Rapid Communication we report an experimental observation of precursor solitons in a flowing dusty plasma. The nonlinear solitary dust acoustic waves (DAWs) are excited by a supersonic mass flow of the dust particles over an electrostatic potential hill. In a frame where the fluid is stationary and the hill is moving the solitons propagate in the upstream direction as precursors while wake structures consisting of linear DAWs are seen to propagate in the downstream region. A theoretical explanation of these excitations based on the forced Korteweg–deVries model equation is provided and their practical implications in situations involving a charged object moving in a plasma are discussed.

We consider the asymmetric simple exclusion process with Langmuir kinetics on a periodic lattice. We analytically obtain the exact time evolution of correlation functions with arbitrary length starting from the initial state with no particle in the system. The exact stationary state of this model has been known for the totally asymmetric case. We propose a basis transformation which simplifies the proof of the stationarity of this state and enables the generalization to the partially asymmetric case. Moreover, we construct low-energy excitations and obtain the exact relaxation time.

This paper extends the Lorentz-Abraham model of an electron (i.e., the equations of motion for a small spherical shell of charge, which is rigid in its proper frame) to treat a small, spherically symmetric charge distribution, allowing for small internal dipolar motion. This is done by dividing the distribution into thin spherical shells (in the continuum limit), and tracking the interactions between shells. Dipolar motion of each constituent spherical shell is allowed along the net dipole moment, but higher-order multipole moments are ignored. The amplitude of dipolar motion of each spherical shell is assumed to be linearly proportional to the net dipole moment. Under these assumptions, low-velocity equations of motion are determined for both the center-of-mass motion and net dipolar motion of the distribution. This is then generalized to arbitrary (relativistic) center-of-mass velocity and acceleration, assuming the motion of individual shells is completely in phase or out of phase with the net dipole moment.

We demonstrate that a three-terminal harmonic symmetric chain in the presence of a Coriolis force, produced by a rotating platform that is used to place the chain, can produce thermal rectification. The direction of heat flow is reconfigurable and controlled by the angular velocity $\mathrm{\Omega }$ of the rotating platform. A simple three-terminal triangular lattice is used to demonstrate the proposed principle.

We investigate a kinetic heat engine model composed of particles enclosed in a box where one side acts as a thermostat and the opposite side is a piston exerting a given pressure. Pressure and temperature are varied in a cyclical protocol of period $\tau$: their relative excursions, $\delta$ and $\epsilon$, respectively, constitute the thermodynamic forces dragging the system out of equilibrium. The analysis of the entropy production of the system allows us to define the conjugated fluxes, which are proportional to the extracted work and the consumed heat. In the limit of small $\delta$ and $\epsilon$ the fluxes are linear in the forces through a $\tau \text{-dependent}$ Onsager matrix whose off-diagonal elements satisfy a reciprocal relation. The dynamics of the piston can be approximated, through a coarse-graining procedure, by a Klein-Kramers equation which—in the linear regime—yields analytic expressions for the Onsager coefficients and the entropy production. A study of the efficiency at maximum power shows that the Curzon-Ahlborn formula is always an upper limit which is approached at increasing values of the thermodynamic forces, i.e., outside of the linear regime. In all our analysis the adiabatic limit $\tau \to \infty$ and the the small-force limit $\delta ,\epsilon \to 0$ are not directly related.

We derive and study two different formalisms used for nonequilibrium processes: the coherent-state path integral, and an effective, coarse-grained stochastic equation of motion. We first study the coherent-state path integral and the corresponding field theory, using the annihilation process $A+A\to A$ as an example. The field theory contains counterintuitive quartic vertices. We show how they can be interpreted in terms of a first-passage problem. Reformulating the coherent-state path integral as a stochastic equation of motion, the noise generically becomes imaginary. This renders it not only difficult to interpret, but leads to convergence problems at finite times. We then show how alternatively an effective coarse-grained stochastic equation of motion with real noise can be constructed. The procedure is similar in spirit to the derivation of the mean-field approximation for the Ising model, and the ensuing construction of its effective field theory. We finally apply our findings to stochastic Manna sandpiles. We show that the coherent-state path integral is inappropriate, or at least inconvenient. As an alternative, we derive and solve its mean-field approximation, which we then use to construct a coarse-grained stochastic equation of motion with real noise.

We study nonoscillating bifurcations of nonhomogeneous steady states of the Vlasov equation, a situation occurring in galactic models, or for Bernstein–Greene–Kruskal modes in plasma physics. Through an unstable manifold expansion, we show that in one spatial dimension the dynamics is very sensitive to the initial perturbation: the instability may saturate at small amplitude—generalizing the “trapping scaling” of plasma physics—or may grow to produce a large-scale modification of the system. Furthermore, resonances are strongly suppressed, leading to different phenomena with respect to the homogeneous case. These analytical findings are illustrated and extended by direct numerical simulations with a cosine interaction potential.

External dc bias-field effects on the nonlinear dielectric relaxation and dynamic Kerr effect of a system of permanent dipoles in a uniaxial mean-field potential are studied via the rotational Brownian motion model postulated in terms of the infinite hierarchy of differential-recurrence equations for the statistical moments $f_n\left(t\right)=\langle P_n\rangle \left(t\right)$ (the expectation value of the Legendre polynomials $P_n$). By solving these equations, the nonlinear dielectric and Kerr-effect ac stationary responses are evaluated for arbitrary dc field strength via perturbation theory in the ac field. Simple analytic equations based on the large separation of the time scales of the fast intrawell and slow overbarrier (interwell) relaxation processes are also derived.

Visibility algorithms transform time series into graphs and encode dynamical information in their topology, paving the way for graph-theoretical time series analysis as well as building a bridge between nonlinear dynamics and network science. In this work we introduce and study the concept of sequential visibility-graph motifs, smaller substructures of $n$ consecutive nodes that appear with characteristic frequencies. We develop a theory to compute in an exact way the motif profiles associated with general classes of deterministic and stochastic dynamics. We find that this simple property is indeed a highly informative and computationally efficient feature capable of distinguishing among different dynamics and robust against noise contamination. We finally confirm that it can be used in practice to perform unsupervised learning, by extracting motif profiles from experimental heart-rate series and being able, accordingly, to disentangle meditative from other relaxation states. Applications of this general theory include the automatic classification and description of physical, biological, and financial time series.

The equilibrium structure of a dinucleosome is studied using an elastic model that takes into account the force and torque balance conditions. Using the proper boundary conditions, it is found that the conformational energy of the problem does not depend on the length of the linker DNA. In addition it is shown that the two histone octamers are almost perpendicular to each other, and the linker DNA in short lengths is almost straight. These findings could shed some light on the role of DNA elasticity in the chromatin structure.

It is numerically demonstrated by means of a magnetohydrodynamic code that precession can trigger dynamo action in a cylindrical container. Fixing the angle between the spin and the precession axis to be $\frac{1}{2}\pi$, two limit configurations of the spinning axis are explored: either the symmetry axis of the cylinder is parallel to the spin axis (this configuration is henceforth referred to as the axial spin case), or it is perpendicular to the spin axis (this configuration is referred to as the equatorial spin case). In both cases, the centro-symmetry of the flow breaks when the kinetic Reynolds number increases. Equatorial spinning is found to be more efficient in breaking the centro-symmetry of the flow. In both cases, the average flow in the reference frame of the mantle converges to a counter-rotation with respect to the spin axis as the Reynolds number grows. We find a scaling law for the average kinetic energy in term of the Reynolds number in the axial spin case. In the equatorial spin case, the unsteady asymmetric flow is shown to be capable of sustaining dynamo action in the linear and nonlinear regimes. The magnetic field is mainly dipolar in the equatorial spin case, while it is is mainly quadrupolar in the axial spin case.

Three-dimensional simulations with fully resolved hydrodynamics are performed to study the collective motion of model swimmers in bulk and confinement. Calculating the dynamic structure factor, we clarified that the swarming in bulk systems can be understood as a pseudoacoustic mode. Under confinement between flat parallel walls, this pseudoacoustic mode leads to a traveling wavelike motion. This swarming behavior is due purely to the hydrodynamic interactions between the swimmers and depends strongly on the type and strength of swimming (i.e., pusher or puller).

This paper presents a numerical simulation study of the combustion structure of streamwise vortex breakdown at Mach number 2.48. Hydrogen fuel is injected into a combustor at sonic speed from the rear of a hypermixer strut that can generate streamwise vortices. The results show that the burning behavior is enhanced at the points of the shock waves that are incident on the vortex and therefore the vortex breakdown in the subsonic region occurs due to combustion. The breakdown domain in the mainstream is found to form a flame-holding region suited to combustion and to lead to a stable combustion field with detached flames. In this way, streamwise vortex breakdown has an essential role in combustion enhancement and the formation of flames that hold under supersonic inflow conditions. Finally, the combustion property defined here is shown to coincide with the produced-water mass flow. This property shows that the amount of combustion is saturated at equivalence ratios over $0.4$, although there is a slight increase beyond 1.

The Jarzynski equality (JE) method, which relates the work of a nonequilibrium process to the free-energy difference between its initial and final states, provides an efficient way to calculate free energies of thermodynamic systems in simulations or experiments. However, more extensive applications of the JE are hindered by the requirement that the initial state must be in equilibrium. In this work we extend the JE method to be the Jarzynski matrix equality (JME) method, which relates the work of trajectories connecting metastable conformational regions to their local free energies, and thus we can estimate the free energy from the nonequilibrium trajectories starting from an almost arbitrary initial distribution. We then apply the JME to toy models, Lennard-Jones fluids, and polymer chain models, demonstrating its efficiency in free-energy calculations with satisfactory accuracy. The JME extends the applicability of the nonequilibrium methods to complex systems whose initial equilibrium states are difficult to reach.

We consider the performance of periodically driven stochastic heat engines in the linear response regime. Reaching the theoretical bounds for efficiency and efficiency at maximum power typically requires full control over the design and the driving of the system. We develop a framework which allows us to quantify the role that limited control over the system has on the performance. Specifically, we show that optimizing the driving entering the work extraction for a given temperature protocol leads to a universal, one-parameter dependence for both maximum efficiency and maximum power as a function of efficiency. In particular, we show that reaching Carnot efficiency (and, hence, Curzon-Ahlborn efficiency at maximum power) requires to have control over the amplitude of the full Hamiltonian of the system. Since the kinetic energy cannot be controlled by an external parameter, heat engines based on underdamped dynamics can typically not reach Carnot efficiency. We illustrate our general theory with a paradigmatic case study of a heat engine consisting of an underdamped charged particle in a modulated two-dimensional harmonic trap in the presence of a magnetic field.

Transient dynamics are of large interest in many areas of science. Here, a generalization of basin stability (BS) is presented: constrained basin stability (CBS) that is sensitive to various different types of transients arising from finite size perturbations. CBS is applied to the paradigmatic Lorenz system for uncovering nonlinear precursory phenomena of a boundary crisis bifurcation. Further, CBS is used in a model of the Earth's carbon cycle as a return time-dependent stability measure of the system's global attractor. Both case studies illustrate how CBS's sensitivity to transients complements BS in its function as an early warning signal and as a stability measure. CBS is broadly applicable in systems where transients matter, from physics and engineering to sustainability science. Thus CBS complements stability analysis with BS as well as classical linear stability analysis and will be a useful tool for many applications.