We study the breakdown of Anderson localization in the one-dimensional nonlinear Klein-Gordon chain, a prototypical example of a disordered classical many-body system. A series of numerical works indicate that an initially localized wave packet spreads polynomially in time, while analytical studies rather suggest a much slower spreading. Here, we focus on the decorrelation time in equilibrium. On the one hand, we provide a mathematical theorem establishing that this time is larger than any inverse power law in the effective anharmonicity parameter $\lambda$, and on the other hand our numerics show that it follows a power law for a broad range of values of $\lambda$. This numerical behavior is fully consistent with the power law observed numerically in spreading experiments, and we conclude that the state-of-the-art numerics may well be unable to capture the long-time behavior of such classical disordered systems.

We investigate the convergence of chemical reaction networks (CRNs), aiming to establish an upper bound on their reaction rates. The nonlinear characteristics and discrete composition of CRNs pose significant challenges in this endeavor. To circumvent these complexities, we adopt an information geometric perspective, utilizing the natural gradient to formulate a nonlinear system. This system effectively determines an upper bound for the dynamics of CRNs. We corroborate our methodology through numerical simulations, which reveal that our constructed system converges more rapidly than CRNs within a particular class of reactions. This class is defined by the count of chemicals, the highest stoichiometric coefficients in the reactions, and the total number of reactions involved. Further, we juxtapose our approach with traditional methods, illustrating that the latter falls short in providing an upper bound for CRN reaction rates. Although our investigation centers on CRNs, the widespread presence of hypergraphs across various disciplines, ranging from natural sciences to engineering, indicates potential wider applications of our method, including in the realm of information science.

We study effects of the mutant's degree on the fixation probability, extinction, and fixation times in Moran processes on Erdös-Rényi and Barabási-Albert graphs. We performed stochastic simulations and used mean-field-type approximations to obtain analytical formulas. We showed that the initial placement of a mutant has a significant impact on the fixation probability and extinction time, while it has no effect on the fixation time. In both types of graphs, an increase in the degree of an initial mutant results in a decreased fixation probability and a shorter time to extinction. Our results extend previous ones to arbitrary fitness values.

This paper is concerned with the processes of spatial propagation and penetration of turbulence from the regions where it is locally excited into initially laminar regions. The phenomenon has come to be known as “turbulence spreading” and witnessed a renewed attention in the literature recently. Here, we propose a comprehensive theory of turbulence spreading based on fractional kinetics. We argue that the use of fractional-derivative equations permits a general approach focusing on fundamentals of the spreading process regardless of a specific turbulence model and/or specific instability type. The starting point is the Hamiltonian of resonant wave-wave interactions, from which a family of scaling laws for the asymptotic spreading is derived. Both three- and four-wave interactions are considered. The results span from a subdiffusive spreading in the parameter range of weak chaos to avalanche propagation in regimes with population inversion. Attention is paid to how nonergodicity introduces weak mixing, memory and intermittency into spreading dynamics, and how the properties of non-Markovianity and nonlocality emerge from the presence of islands of regular dynamics in phase space. Also we resolve an existing question concerning turbulence spillover into gap regions, where the instability growth is locally suppressed, and show that the spillover occurs through exponential (Anderson-like) localization in case of four-wave interactions and through an algebraic (weak) localization in case of triad interactions. In the latter case an inverse-cubic behavior of the spillover function is found. Wherever relevant, we contrast our findings against the available observational and numerical evidence, and we also commit ourselves to establish connections with the models of turbulence spreading proposed previously.

The Higgs RNA model with an added term for a coupling to an external force is studied in regard to finite-time force-driving protocols with a minimal-work requirement. In this paper, RNA sequences which at low temperature exhibit hairpins are considered, which are often cited as typical template systems in stochastic thermodynamics. The optimized work protocols for this glassy many-particle system are determined numerically using the parallel tempering method. The protocols show distinct jumps at the beginning and end, which have been observed for single-particle systems and are proven to be optimal in the fast protocol limit generally. Optimality seems to be achieved by staying close to the equilibrium unfolding transition point, in agreement with experimental and theoretical observations. The change of work distributions, compared to those resulting from a naive linear driving protocol, are discussed generally and in terms of free energy estimation as well as the effect of optimized protocols on rare work process starting conditions.

We study quantum Otto thermal machines with a two-spin working system coupled by anisotropic interaction. Depending on the choice of different parameters, the quantum Otto cycle can function as different thermal machines, including a heat engine, refrigerator, accelerator, and heater. We aim to investigate how the anisotropy plays a fundamental role in the performance of the quantum Otto engine (QOE) operating in different timescales. We find that while the engine's efficiency increases with the increase in anisotropy for the quasistatic operation, quantum internal friction and incomplete thermalization degrade the performance in a finite-time cycle. Further, we study the quantum heat engine (QHE) with one of the spins (local spin) as the working system. We show that the efficiency of such an engine can surpass the standard quantum Otto limit, along with maximum power, thanks to the anisotropy. This can be attributed to quantum interference effects. We demonstrate that the enhanced performance of a local-spin QHE originates from the same interference effects, as in a measurement-based QOE for their finite-time operation.

The COVID-19 pandemic has underscored the importance of understanding, forecasting, and avoiding infectious processes, as well as the necessity for understanding the diffusion and acceptance of preventative measures. Simple contagions, like virus transmission, can spread with a single encounter, while complex contagions, such as preventive social measures (e.g., wearing masks, social distancing), may require multiple interactions to propagate. This disparity in transmission mechanisms results in differing contagion rates and contagion patterns between viruses and preventive measures. Furthermore, the dynamics of complex contagions are significantly less understood than those of simple contagions. Stochastic models, integrating inherent variability and randomness, offer a way to elucidate complex contagion dynamics. This paper introduces a stochastic model for both simple and complex contagions and assesses its efficacy against ensemble simulations for homogeneous and heterogeneous threshold configurations. The model provides a unified framework for analyzing both types of contagions, demonstrating promising outcomes across various threshold setups on Erds-Rényi graphs.

We report on experiments where solid-density Mg plasmas are created by heating with the focused output of the Linac Coherent Light Source x-ray free-electron laser. We study the K-shell emission from the helium- and lithium-like ions using Bragg crystal spectroscopy. Observation of the dielectronic satellites in lithium-like ions confirms that the M-shell electrons appear bound for these high charge states. An analysis of the intensity of these satellites indicates that when modeled with an atomic-kinetics code, the ionization potential depression model employed needs to produce depressions for these ions which lie between those predicted by the well known Stewart-Pyatt and Ecker-Kroll models. These results are largely consistent with recent density functional theory calculations.

The Lagrange-mesh method is an approximate variational method which provides accurate solutions of the Schrödinger equation for bound-state and scattering few-body problems. The stationary Klein-Gordon equation depends quadratically on the energy. For a central potential, it is solved on a Lagrange-Laguerre mesh by iteration. Results are tested with the Coulomb potential for which exact solutions are available. A high accuracy is obtained with a rather small number of mesh points. For various potentials and levels, few iterations provide accurate energies and mean values in short computer times. Analytical expressions of the wave functions are available.

We study heat transfer in one-dimensional Fermi-Pasta-Ulam-Tsingou-type systems with long-range (LR) interactions. The strength of the LR interaction between two lattice sites decays as a power $\sigma$ of the inverse of their distance. We focus on the strong LR regime ($0\le \sigma \le 1$) and show that the thermal transport behaviors are remarkably nuanced. Specifically, we observe that the antipersistent (negative) energy current correlation in this regime is intricately dependent on $\sigma$, displaying a nonmonotonic variation. Notably, a significant qualitative change occurs at ${\sigma }_c=0.5$, where with respect to other $\sigma$ values the correlation shows a minimum negative value. Furthermore, our findings also demonstrate that within the long-time range considered, these antipersistent correlations will eventually vanish for certain $\sigma >0.5$. The underlying mechanisms behind these intriguing phenomena are related to the crossover of two diverse space-time scaling properties of equilibrium heat correlations and the various scattering processes of phonons and discrete breathers.

The probability distribution (PD) of spin configurations in kinetic Ising models has been cast in the form of the canonical Boltzmann PD with a time-dependent effective Hamiltonian (EH). It has been argued that in systems with extensive energy EH depends linearly on the number of spins $N$ leading to the exponential dependence of PD on the system size. In macroscopic systems the argument of the exponential function may reach values of the order of the Avogadro number which is impossible to deal with computationally, thus making unusable the linear master equation (ME) governing the PD evolution. To overcome the difficulty, it has been suggested to use instead the nonlinear ME (NLME) for the EH density per spin. It has been shown that in spatially homogeneous systems NLME contains only terms of order unity even in the thermodynamic limit. The approach has been illustrated with the kinetic Husimi-Temperley model (HTM) evolving under the Glauber dynamics. At finite $N$ the known numerical results has been reproduced and extended to broader parameter ranges. In the thermodynamic limit an exact nonlinear partial differential equation of the Hamilton-Jacobi type for EH has been derived. It has been shown that the average magnetization in HTM evolves according to the conventional kinetic mean field equation.

Thermal conduction force plays a crucial role in manipulating the local thermal conductivity of crystals. However, due to the diffusive nature of thermal conduction, investigating the force effect is challenging. Recently, researchers have explored the force effect based on the wavelike behavior of thermal conduction, specifically second sound. However, their focus has been primarily on the progressive case, neglecting the more complex standing temperature field case. Additionally, establishing a connection between the results obtained from the progressive case and the standing case poses a challenging problem. In this study, we investigate the force effect of standing and quasistanding temperature fields, revealing distinct characteristics of thermal conduction force. Moreover, we establish a link between the progressive and standing cases through the quasistanding case. Our findings pave the way for research in more intricate scenarios and provide an additional degree of freedom for manipulating the local thermal conductivity of dielectric crystals.

Non-Markovian dynamics of a charged particle in a two-dimensional harmonic oscillator linearly coupled to a neutral bosonic heat bath is investigated in an external uniform magnetic field and two perpendicular time-dependent electric fields. The analytical expressions for the time-dependent and asymptotic angular momentum are derived for the Markovian and non-Markovian dynamics. The dependence of the angular momentum on the frequency of the electric field, cyclotron frequency, collective frequency, and anisotropy of the heat bath is studied. The angular momentum (or magnetization) of a charged particle can be ruled by varying the frequency of the electric field.

Bifurcations in dynamical systems are often studied experimentally and numerically using a slow parameter sweep. Focusing on the cases of period-doubling and pitchfork bifurcations in maps, we show that the adiabatic approximation always breaks down sufficiently close to the bifurcation, so the upsweep and downsweep dynamics diverge from one another, disobeying standard bifurcation theory. Nevertheless, we demonstrate universal upsweep and downsweep trajectories for sufficiently slow sweep rates, revealing that the slow trajectories depend essentially on a structural asymmetry parameter, whose effect is negligible for the stationary dynamics. We obtain explicit asymptotic expressions for the universal trajectories and use them to calculate the area of the hysteresis loop enclosed between the upsweep and downsweep trajectories as a function of the asymmetry parameter and the sweep rate.

Heider's structural balance theory has proven invaluable in comprehending the dynamics of social groups characterized by both friendly and hostile relationships. Since people's relations are rarely single faceted, we investigate Heider balance dynamics on a multiplex network, consisting of several copies of the same agent displaying correlated relations at different layers building the multiplex. Intralayer interactions in our model adhere to Heider dynamics, while interlayer correlations stem from Ising interactions, with the heat-bath dynamics of link signs. Our investigation reveals a multifaceted system with a diverse equilibrium landscape contingent on the coexistence of distinct phases across layers. We observe that, starting from a paradise state with positive links in all layers, an increase in temperature triggers a discontinuous transition to a disordered state akin to single-layer scenarios. The critical temperature surpasses that of the single-layer case, a fact verified through extended mean-field analysis and agent-based simulations. Furthermore, the scenario shifts when one layer exhibits a two-clique configuration instead of a paradise state. This change introduces additional transitions: synchronization of interlayer relations and a transition to the disorder, appearing at a different, lower temperature compared to matching paradise states. This exploration shows the intricate interplay of Heider balance and multiplex interactions.

Active propulsion, as performed by bacteria and Janus particles, in combination with hydrodynamic interaction results in the accumulation of bacteria at a flat wall. However, in microfluidic devices with cylindrical pillars of sufficiently small radius, self-propelled particles can slide along and scatter off the surface of a pillar, without becoming trapped over long times. This nonequilibrium scattering process has been predicted to result in large diffusivities, even at high obstacle density, unlike particles that undergo classical specular reflection. Here, we test this prediction by experimentally studying the nonequilibrium scattering of pusherlike swimmers in microfluidic obstacle lattices. To explore the role of tumbles in the scattering process, we microscopically tracked wild-type (run and tumble) and smooth-swimming (run only) mutants of the bacterium Escherichia coli scattering off microfluidic pillars. We quantified key scattering parameters and related them to previously proposed models that included a prediction for the diffusivity, discussing their relevance. Finally, we discuss potential interpretations of the role of tumbles in the scattering process and connect our work to the broader study of swimmers in porous media.

Chromatin polymer dynamics are commonly described using the classical Rouse model. The subsequent discovery, however, of intermediate-scale chromatin organization known as topologically associating domains (TADs) in experimental Hi-C contact maps for chromosomes across the tree of life, together with the success of loop extrusion factor (LEF) model in explaining TAD formation, motivates efforts to understand the effect of loops and loop extrusion on chromatin dynamics. This paper seeks to fulfill this need by combining LEF-model simulations with extended Rouse-model polymer simulations to investigate the dynamics of chromatin with loops and dynamic loop extrusion. We show that loops significantly suppress the averaged mean-square displacement (MSD) of a gene locus, consistent with recent experiments that track fluorescently labeled chromatin loci. We also find that loops reduce the MSD's stretching exponent from the classical Rouse-model value of 1/2 to a loop-density-dependent value in the 0.45–0.40 range. Remarkably, stretching exponent values in this range have also been observed in recent experiments [Weber et al., Phys. Rev. Lett. 104, 238102 (2010); Bailey et al., Mol. Biol. Cell 34, ar78 (2023)]. We also show that the dynamics of loop extrusion itself negligibly affects chromatin mobility. By studying static “rosette” loop configurations, we also demonstrate that chromatin MSDs and stretching exponents depend on the location of the locus in question relative to the position of the loops and on the local friction environment.

We investigate the formation of wrinkling instabilities at the interface between layers of hydrogel and water, which arise to relieve horizontal compressive stresses caused by either differential swelling or confinement. Modelling the gel using a linear-elastic-nonlinear-swelling approach, we determine both a criterion for marginal stability and the growth rates of normal modes. Furthermore, our formalism allows us to understand the influence of differential swelling on the stability of hydrogels brought into contact with water, and we find three distinct phases of the instability. Initially, when only a thin skin layer of gel has swollen, buckles grow rapidly and the gel deforms as an incompressible material. A balance between normal elastic stress and pore pressure selects a wavelength for these buckles that increases with the square root of time. At late times, when the gel approaches a uniformly swollen state, buckles can only grow by differential swelling on much slower timescales determined by solvent transport. At intermediate times, growth is driven by the same fluid transport process as at late times but gradients in fluid pressure in the gel as it swells destabilize the interface, driving faster growth of wrinkles. We also explain why some instabilities can be transient, “healing” as time progresses, while others must remain for all time.

We present a simple model of driven matter in a 1D medium with pinning impurities, applicable to magnetic domains walls, confined colloids, and other systems. We find rich dynamics, including hysteresis, reentrance, quasiperiodicity, and two distinct routes to chaos. In contrast to other minimal models of driven matter, the model is solvable: we derive the full phase diagram for small $N$, and for large $N$, we derive expressions for order parameters and several bifurcation curves. The model is also realistic. Its collective states match those seen in the experiments of magnetic domain walls.