We describe the steady state of the annihilation process of a one-dimensional system of two initially separated reactants A and B. The parameters that define the dynamical behavior of the system are the diffusion constant, the reaction rate, and the deposition rate. Depending on the ratio between those parameters, the system exhibits a crossover between a diffusion-limited (DL) regime and a reaction-limited (RL) regime. We found that a key quantity to describe the reaction process in the system is the probability $p(x_A,x_B)$ to find the rightmost A (RMA) particle and the leftmost B (LMB) particle at the positions $x_A$ and $x_B$, respectively. The statistical behavior of the system in both regimes is described using the density of particles, the gap length distribution $x_B-x_A$, the marginal probabilities $p_A\left(x_A\right)$ and $p_B\left(x_B\right)$, and the reaction kernel. For both regimes, this kernel can be approximated by using $p(x_A,x_B)$. We found an excellent agreement between the numerical and analytical results for all calculated quantities despite the reaction process being quite different in both regimes. In the DL regime, the reaction kernel can be approximated by the probability to find the RMA and LMB particles in adjacent sites. In the RL regime, the kernel depends on the marginal probabilities $p_A\left(x_A\right)$ and $p_B\left(x_B\right)$.

We introduce and demonstrate the usage of the origin-fate map (OFM) as a tool for the detailed investigation of phase space transport in reactant-product-type systems. For these systems, which exhibit clearly defined start and end states, it is possible to build a comprehensive picture of the lobe dynamics by considering backward and forward integration of sets of initial conditions to index their origin and fate. We illustrate the method and its utility in the study of a two degrees of freedom caldera potential with four exits, demonstrating that the OFM not only recapitulates results from classical manifold theory but even provides more detailed information about complex lobe structures. The OFM allows the detection of dynamically significant transitions caused by the creation of new lobes and is also able to guide the prediction of the position of unstable periodic orbits (UPOs). Further, we compute the OFM on the periodic orbit dividing surface (PODS) associated with the transition state of a caldera entrance, which allows for a powerful analysis of reactive trajectories. The intersection of the manifolds corresponding to this UPO with other manifolds in the phase space results in the appearance of lobes on the PODS, which are directly classified by the OFM. This allows computations of branching ratios and the exploration of a fractal cascade of lobes as the caldera is stretched, which results in fluctuations in the branching ratio and chaotic selectivity. The OFM is found to be a simple and very useful tool with a vast range of descriptive and quantitative applications.

The purpose of this work is to use active particles to investigate the effect of facilitation on supercooled liquids. To this end we examine the behavior of a model supercooled liquid that is doped with a mixture of active particles and slowed particles. To simulate the facilitation mechanism, the activated particles are subjected to a force that follows the mobility of their most mobile neighboring molecule, while the slowed particles experience a friction force. Upon activation, we observe a fluidization of the entire medium along with a significant increase in dynamic heterogeneity. This effect is reminiscent of the fluidization observed experimentally when introducing molecular motors into soft materials. Interestingly, when the characteristic time ${\tau }_{\mu }$, used to define the mobility in the facilitation mechanism, matches the physical time $t^{*}$ that characterizes the spontaneous cooperativity of the material, we observe a phase transition accompanied by structural aggregation of the active molecules. This transition is characterized by a sharp increase in fluidization and dynamic heterogeneity.

The study by Semenova et al. [Phys. Rev. Lett. 117, 014102 (2016)] discovered a type of chimera state known as coherence-resonance chimera (CRC), which combines the effects of coherence resonance (CR) and the spatial property of classical chimeras. In this Letter, we present yet another form of chimera, which we refer to as self-induced-stochastic-resonance breathing chimera (SISR-BC), which differs fundamentally from the CRC in that it combines the mechanism and effects of self-induced stochastic resonance (SISR, previously shown by DeVille et al. [Phys. Rev. E 72, 031105 (2005)] to be intrinsically different from CR), the symmetry breaking in the rotational coupling between the slow and fast subsystems of the coupled oscillators, and the property of breathing chimera—a form of chimera state characterized by nonstationary periodic dynamics of coherent-incoherent patterns with a periodically oscillating global order parameter. Unlike other types of chimeras, including CRC, SISR-BC demonstrates remarkable resilience to a relatively wide range of stochastic perturbations and persists even when the purely excitable system is significantly distant from the Hopf bifurcation threshold—thanks to the mechanism of SISR—and globally attracts random distributions of initial conditions. Considering its potential impact on information processing in neuronal networks, SISR-BC could have special significance and applications.

Designing a protocol to efficiently drive a stochastic system is an active field of research. Here we extend such control theory to an active Ornstein-Uhlenbeck particle (AOUP) in a bistable potential, driven by a harmonic trap. We find that protocols designed to minimize the excess work (up to linear response) perform better than naive protocols with constant velocity for a wide range of protocol durations.

We present the results of linear stability of a damped coplanar double pendulum and its nonlinear motion, when the point of suspension is vibrated sinusoidally in the vertical direction with amplitude $a$ and frequency $\omega$. A double pendulum has two pairs of Floquet multipliers, which have been calculated for various driving parameters. We have considered the stability of a double pendulum when it is in any of its possible stationary states: (i) both pendulums are either vertically downward or upward and (ii) one pendulum is downward and other is upward. The damping is considered to be velocity-dependent, and the driving frequency is taken in a wide range. A double pendulum excited from its stable state shows both periodic and chaotic motion. The periodic motion about its pivot may be either oscillatory or rotational. The periodic swings of a driven double pendulum may be either harmonic or subharmonic for lower values of $a$. The limit cycles corresponding to the normal mode oscillations of a double pendulum of two equal masses are squeezed into a line in its configuration space. For unequal masses, the pendulum shows multi-period swings for smaller values of $a$ and damping, while chaotic swings or rotational motion at relatively higher values of $a$. The parametric driving may lead to stabilization of a partially or fully inverted double pendulum.

The elastic response of dense suspensions under an impact is studied using coupled lattice Boltzmann method and discrete element method (LBM-DEM) and its reduced model. We succeed to extract the elastic force acting on the impactor in dense suspensions, which can exist even in the absence of percolating clusters of suspended particles. We then propose a reduced model to describe the motion of the impactor and demonstrate its relevancy through the comparison of the solution of the reduced model and that of LBM-DEM. Furthermore, we illustrate that the perturbation analysis of the reduced model captures the short-time behavior of the impactor motion quantitatively. We apply this reduced model to the impact of a foot-spring-body system on a dense suspension, which is the minimal model to realize walking on the suspension. Due to the spring force of the system and the stiffness of the suspension, the foot undergoes multiple bounces. We also study the parameter dependencies of the hopping motion and find that multiple bounces are suppressed as the spring stiffness increases.

Thermodynamic properties of fluids confined in nanopores differ from those observed in the bulk. To investigate the effect of nanoconfinement on water compressibility, we perform water sorption experiments on two nanoporous glass samples while concomitantly measuring the speed of longitudinal and shear ultrasonic waves in these samples. These measurements yield the longitudinal and shear moduli of the water-laden nanoporous glass as a function of relative humidity that we utilize in the Gassmann theory to infer the bulk modulus of the confined water. This analysis shows that the bulk modulus (inverse of compressibility) of confined water is noticeably higher than that of the bulk water at the same temperature. Moreover, the modulus exhibits a linear dependence on the Laplace pressure. The results for water, which is a polar fluid, agree with previous experimental and numerical data reported for nonpolar fluids. This similarity suggests that irrespective of intermolecular forces, confined fluids are stiffer than bulk fluids. Accounting for fluid stiffening in nanopores may be important for accurate interpretation of wave propagation measurements in fluid-filled nanoporous media, including in petrophysics, catalysis, and other applications, such as in porous materials characterization.

We report an instability of a slider slowly dragged at the surface of a granular bed in a quasistatic regime. The boat-shaped slider sits on the granular medium under its own weight and is free to translate vertically and to rotate around the pitch axis while a constant horizontal speed is imposed. For a wide range of parameters (mass, length, shape, velocity) a regular pattern of peaks and troughs spontaneously emerges as the slider travels forward. This instability is studied through experiments using a conveyor belt and by means of two-dimensional discrete elements method simulations. We show that the wavelength and amplitude of the pattern scale as the length of the slider. We also observe that the ripples disappear for low and high masses, indicating an optimal confining pressure. The effect of the shape, more specifically the inclination of the front spatula, is studied and found to drastically influence both the wavelength and the amplitude. Finally, we show that the mechanical details (friction, cohesion) of the contact point between the slider and the pulling device is critical and remains to be fully understood.

Understanding the organization of matter under the long-range electrostatic force is a fundamental problem in multiple fields. In this work, based on the electrically charged tethered membrane model, we reveal regular structures underlying the lowest-energy states of inhomogeneously stretched planar lattices by a combination of numerical simulation and analytical geometric analysis. Specifically we show the conformal order characterized by the preserved bond angle in the lattice deformation and reveal the Poincaré-Klein mapping underlying the electrostatics-driven inhomogeneity. The discovery of the Poincaré-Klein mapping, which connects the Poincaré disk and the Klein disk for the hyperbolic plane, implies the connection of long-range electrostatic force and hyperbolic geometry. We also discuss lattices with patterned charges of opposite signs for modulating in-plane inhomogeneity and even creating 3D shapes, which may have a connection to metamaterials design. This work suggests the geometric analysis as a promising approach for elucidating the organization of matter under the long-range force.

Preionization is believed to play an important role on the implosion of gas-puff Z pinches. Some experiments have used an external preionization source, e.g., UV light or electron beam. In contrast, other experiments rely completely on over voltage breakdown by the own generator's voltage pulse. However, this approach lacks shot-to-shot reproducibility since self-breakdown is mainly a stochastic process. In this work, we performed a systematic study on self-breakdown using two different cathode geometries: (i) a smooth, round cathode to provide a homogeneous electric field, (ii) a sharp, knife-edge-like geometry to enhance the electric field locally and eventually electron emission. The experiments were carried out on the Llampudken current generator, which provides a current pulse of $\sim 400\text{kA}$ amplitude and 200 ns rise time ($10\%\text{−}90\%$). We implemented gated XUV imaging, filtered diodes and time-integrated x-ray imaging to obtain information about the implosion as well as the stagnation phase for the two cathode geometries. We found that erosion of the knife-edge cathode might be a serious problem, and we had to replace it every 15 shots. On the other hand, the round cathode lasted for the whole series of experiments. We also measured a more reproducible and larger peak current for the knife cathode. From the photo-conductive detectors we observed that even if the round cathode might present shots with higher x-ray yield compared to the knife cathode, dispersion is almost twice as large. Moreover, after a statistic analysis, it is demonstrated that the dispersion in the yield is due solely to differences imposed by the cathodes and not to variations in the driver, as no correlation was found between them. We found that in order to fit the experimental data with the snowplow model, only $\sim 60\%$ of the total mass is compressed in the knife cathode while $\sim 20\%$ for the round one, highlighting the importance of the cathode and preionization. Therefore, we conclude that the use of the knife cathode increases the reproducibility of the experiment in comparison with the round cathode.

Degradation in the energy spread of accelerated protons due to the transverse instability induced transparency is one of the critical issues in the laser-driven radiation pressure acceleration (RPA) scheme. This issue is more severe for linearly polarized lasers due to enhanced heating of electrons. Therefore, in spite of being experimentally challenging, most of the numerical studies are performed with circularly polarized lasers. In this work, through particle-in-cell simulations, we demonstrate a significant improvement in the energy spread of the accelerated protons when a multilayered target is irradiated by a linearly polarized laser. This multilayered target consists of a near-critical-density (NCD) layer, sandwiched between a thick metallic foil and a thin RPA target. The role of the NCD target is to suppress the laser transparency to increase the coupling of laser momentum to the RPA protons. On the other hand, the metallic foil utilizes the kinetic energy of the escaping fast electrons to form an electrostatic sheath to filter the low-energy RPA protons. This results in significant improvement in the accelerated proton spectrum, even with a linearly polarized laser.

Understanding the flow of yield stress fluids in porous media is a major challenge. In particular, experiments and extensive numerical simulations report a nonlinear Darcy law as a function of the pressure gradient. In this letter we consider a treelike porous structure for which the problem of the flow can be resolved exactly due to a mapping with the directed polymer (DP) with disordered bond energies on the Cayley tree. Our results confirm the nonlinear behavior of the flow and expresses its full pressure dependence via the density of low-energy paths of DP restricted to vanishing overlap. These universal predictions are confirmed by extensive numerical simulations.

We determine the full distribution and moments of the first passage time for a wide class of stochastic search processes in the limit of frequent stochastic resetting. Our results apply to any system whose short-time behavior of the search process without resetting can be estimated. In addition to the typical case of exponentially distributed resetting times, we prove our results for a variety of resetting time distributions. We illustrate our results in several examples and show that the errors of our approximations vanish exponentially fast in typical scenarios of diffusive search.

We perform a numerical study of transport properties of a one-dimensional chain with couplings decaying as an inverse power $r^{-(1+\sigma )}$ of the intersite distance $r$ and open boundary conditions, interacting with two heat reservoirs. Despite its simplicity, the model displays highly nontrivial features in the strong long-range regime $-1<\sigma <0$. At weak coupling with the reservoirs, the energy flux departs from the predictions of perturbative theory and displays anomalous superdiffusive scaling of the heat current with the chain size. We trace this behavior back to the transmission spectrum of the chain, which displays a self-similar structure with a characteristic $\sigma$-dependent fractal dimension.

A spin system is studied with simultaneous permutation-symmetric Potts and spin-rotation-symmetric clock interactions in spatial dimensions $d=2$ and 3. The global phase diagram is calculated from the renormalization-group solution with the recently improved (spontaneous first-order detecting) Migdal-Kadanoff approximation or, equivalently, with hierarchical lattices with the inclusion of effective vacancies. Five different ordered phases are found: Conventionally ordered ferromagnetic, quadrupolar, antiferromagnetic phases and algebraically ordered antiferromagnetic, antiquadrupolar phases. These five different ordered phases and the disordered phase are mutually bounded by first- and second-order phase transitions, themselves delimited by multicritical points: Inverted bicritical, zero-temperature bicritical, tricritical, second-order bifurcation, and zero-temperature highly degenerate multicritical points. One rich phase diagram topology exhibits all of these phenomena.

The out-of-time order correlator (OTOC) has been widely studied in closed quantum systems. However, there are very few studies for open systems and they are mainly focused on isolating the effects of scrambling from those of decoherence. Adopting a different point of view, we study the interplay between these two processes. This proves crucial in order to explain the OTOC behavior when a phase space contracting dissipation is present, ubiquitous not only in real life quantum devices but in the dynamical systems area. The OTOC decay rate is closely related to the classical Lyapunov exponent—with some differences—and more sensitive in order to distinguish the chaotic from the regular behavior than other measures. On the other hand, it is revealed as a generally simple function of the longest lived eigenvalues of the quantum evolution operator. We find no simple connection with the Ruelle-Pollicott resonances, but by adding Gaussian noise of ${\hslash }_{\text{eff}}$ size to the classical system we recover the OTOC decay rate, which is a consequence of the correspondence principle put forward in Phys. Rev. Lett. 108, 210605 (2012) and Phys. Rev. E 99, 042214 (2019).

This study presents an improved solution procedure for determining the transient probability density function (PDF) solutions of the nonlinear oscillators with both odd and even nonlinearity under modulated random stimulation, which is an extension of the exponential-polynomial-closure approach. An evolutionary exponential-polynomial function with time-varying undetermined variables is considered as the transient probabilistic solution. By selecting a set of independent evolutionary base functions spanning a $R^n$ space as weight functions, a set of ordinary differential equations can be formulated by integrating the weighted residual error. The undetermined variables can be determined numerically by solving those ordinary differential equations. Three numerical examples illustrate that the improved solution procedure can acquire the transient probabilistic responses of the stochastic dynamic systems effectively and efficiently even in the PDF solution tails when compared with Monte Carlo simulation. Moreover, the results indicate that the PDF solutions of the oscillators are asymmetrical at their nonzero means due to the influence of even nonlinearity. The nonstationary behaviors of the system responses are also investigated along with the behavior of modulated Gaussian white noise.

The contact process is an emblematic model of a nonequilibrium system, containing a phase transition between inactive and active dynamical regimes. In the epidemiological context, the model is known as the susceptible-infected-susceptible model, and it is widely used to describe contagious spreading. In this work, we demonstrate how accurate and efficient representations of the full probability distribution over all configurations of the contact process on a one-dimensional chain can be obtained by means of matrix product states (MPSs). We modify and adapt MPS methods from many-body quantum systems to study the classical distributions of the driven contact process at late times. We give accurate and efficient results for the distribution of large gaps, and illustrate the advantage of our methods over Monte Carlo simulations. Furthermore, we study the large deviation statistics of the dynamical activity, defined as the total number of configuration changes along a trajectory, and investigate quantum-inspired entropic measures, based on the second Rényi entropy.

Achieving perfect synchronization in a complex network, specially in the presence of higher-order interactions (HOIs) at a targeted point in the parameter space, is an interesting, yet challenging task. Here we present a theoretical framework to achieve the same under the paradigm of the Sakaguchi-Kuramoto (SK) model. We analytically derive a frequency set to achieve perfect synchrony at some desired point in a complex network of SK oscillators with higher-order interactions. Considering the SK model with HOIs on top of the scale-free, random, and small world networks, we perform extensive numerical simulations to verify the proposed theory. Numerical simulations show that the analytically derived frequency set not only provides stable perfect synchronization in the network at a desired point but also proves to be very effective in achieving a high level of synchronization around it compared to the other choices of frequency sets. The stability and the robustness of the perfect synchronization state of the system are determined using the low-dimensional reduction of the network and by introducing a Gaussian noise around the derived frequency set, respectively.