We use the discrete element method, taking particle contact and hydrodynamic lubrication into account, to unveil the shear rheology of suspensions of frictionless non-Brownian rods in the dense packing fraction regime. We find that, analogously to the random close packing volume fraction, the shear-driven jamming point of this system varies in a nonmonotonic fashion as a function of the rod aspect ratio. The latter strongly influences how the addition of rodlike particles affects the rheological response of a suspension of frictionless non-Brownian spheres to an external shear flow. At fixed values of the total (rods plus spheres) packing fraction, the viscosity of the suspension is reduced by the addition of “short”($\le 2$) rods but is instead increased by the addition of “long”($\ge 2$) rods. A mechanistic interpretation is provided in terms of packing and excluded-volume arguments.

Despite extensive research, the fundamental physical mechanisms underlying the Mpemba effect, a phenomenon where a substance cools faster after initially being heated, remain elusive. Although historically linked with water, the Mpemba effect manifests across diverse systems, sparking heightened interest in Mpemba-like phenomena. Concurrently, the Kovacs effect, a memory phenomenon observed in materials such as polymers, involves rapid quenching and subsequent temperature changes, resulting in nonmonotonic relaxation behavior. This paper probes the intricacies of the Mpemba and Kovacs effects within the framework of the time-delayed Newton's law of cooling, recognized as a simplistic yet effective phenomenological model accommodating memory phenomena. This law allows for a nuanced comprehension of temperature variations, introducing a delay time ($\tau$) and incorporating specific protocols for the thermal bath temperature, contingent on a defined waiting time ($t_{\text{w}}$). Remarkably, the relevant parameter space is two-dimensional ($\tau$ and $t_{\text{w}}$), with bath temperatures exerting no influence on the presence or absence of the Mpemba effect or on the relative strength of the Kovacs effect. The findings enhance our understanding of these memory phenomena, providing valuable insights applicable to researchers across diverse fields, ranging from physics to materials science.

Swarmalators are oscillators that can swarm as well as sync via a dynamic balance between their spatial proximity and phase similarity. Swarmalator models employed so far in the literature comprise only one-dimensional phase variables to represent the intrinsic dynamics of the natural collectives. Nevertheless, the latter can indeed be represented more realistically by high-dimensional phase variables. For instance, the alignment of velocity vectors in a school of fish or a flock of birds can be more realistically set up in three-dimensional space, while the alignment of opinion formation in population dynamics could be multidimensional, in general. We present a generalized $D$-dimensional swarmalator model, which more accurately captures self-organizing behaviors of a plethora of real-world collectives by self-adaptation of high-dimensional spatial and phase variables. For a more sensible visualization and interpretation of the results, we restrict our simulations to three-dimensional spatial and phase variables. Our model provides a framework for modeling complicated processes such as flocking, schooling of fish, cell sorting during embryonic development, residential segregation, and opinion dynamics in social groups. We demonstrate its versatility by capturing the maneuvers of a school of fish, qualitatively and quantitatively, by a suitable extension of the original model to incorporate appropriate features besides a gallery of its intrinsic self-organizations for various interactions. We expect the proposed high-dimensional swarmalator model to be potentially useful in describing swarming systems and programmable and reconfigurable collectives in a wide range of disciplines, including the physics of active matter, developmental biology, sociology, and engineering.

Ensembles of coupled nonlinear oscillators are a popular paradigm and an ideal benchmark for analyzing complex collective behaviors. The onset of cluster synchronization is found to be at the core of various technological and biological processes. The current literature has investigated cluster synchronization by focusing mostly on the case of attractive coupling among the oscillators. However, the case of two coexisting competing interactions is of practical interest due to their relevance in diverse natural settings, including neuronal networks consisting of excitatory and inhibitory neurons, the coevolving social model with voters of opposite opinions, and ecological plant communities with both facilitation and competition, to name a few. In the present article, we investigate the impact of repulsive spanning trees on cluster formation within a connected network of attractively coupled limit-cycle oscillators. We successfully predict which nodes belong to each cluster and the emergent frustration of the connected networks independent of the particular local dynamics at the network nodes. We also determine local asymptotic stability of the cluster states using an approach based on the formulation of a master stability function. We additionally validate the emergence of solitary states and antisynchronization for some specific choices of spanning trees and networks.

We studied the single-particle Anderson localization problem for non-Hermitian systems on directed graphs. Random regular graph and various undirected standard random graph models were modified by controlling reciprocity and hopping asymmetry parameters. We found the emergence of left, biorthogonal, and right localized states depending on both parameters and graph structure properties such as node degree $d$. For directed random graphs, the occurrence of biorthogonal localization near exceptional points is described analytically and numerically. The clustering of localized states near the center of the spectrum and the corresponding mobility edge for left and right states are shown numerically. Structural features responsible for localization, such as topologically invariant nodes or drains and sources, were also described. Considering the diagonal disorder, we observed the disappearance of localization dependence on reciprocity around $W\sim 20$ for a random regular graph $d=4$. With a small diagonal disorder, the average biorthogonal fractal dimension drastically reduces. Around $W\sim 5$, localization scars occur within the spectrum, alternating as vertical bands of clustering of left and right localized states.

We investigate a two-dimensional system of interacting active Brownian particles. Using the Martin-Siggia-Rose-Janssen-de Dominicis formalism, we built up the generating functional for correlation functions. We study in detail the hydrodynamic regime with a constant density stationary state. Our findings reveal that, within a small density fluctuations regime, an emergent $U\left(1\right)$ gauge symmetry arises, originated from the conservation of fluid vorticity. Consequently, the interaction between the orientational order parameter and density fluctuations can be cast into a gauge theory, where the concept of “electric charge density” aligns with the local vorticity of the original fluid. We study in detail the case of a microscopic local two-body interaction. We show that, upon integrating out the gauge fields, the stationary states of the rotational degrees of freedom satisfy a nonlocal Frank free energy for a nematic fluid. We give explicit expressions for the splay and bend elastic constants as a function of the Péclet number ($\mathrm{Pe}$) and the diffusion interaction constant ($k_d$).

Rhythmic behaviors are generally observed in nonlinear chemical reactions such as the Belousov-Zhabotinsky reaction and enzymatic reactions. Similarly, a simple phase change can also lead to rhythmic behavior. It has been reported previously that camphor solid films alternate between generation and disappearance on ethanol (EtOH) solution, and a phenomenological mechanism has been suggested for this. The evaporation of EtOH decreases the temperature on the surface of the solution via vaporization heat and induces precipitation in the camphor solid film. At this time, the film prevents evaporation, and thus, the surface temperature increases due to thermal diffusion from the atmosphere, resulting in dissolution of the solid film. To verify the previously suggested phenomenological mechanism, we controlled the evaporation rate of EtOH using a porous plastic cover. As a result, the period of oscillation increased with decreasing pore diameter, and finally, the oscillation did not occur without pore in the cover, where the camphor solid film was not observed. Additionally, a new mathematical model was proposed, and the numerical calculations agreed well with experimental observations. Linear stability and bifurcation analyses revealed the detailed mechanism of this phenomenon, which agreed well with the phenomenological explanation mentioned above.

Understanding the mechanical instabilities of two-dimensional membranes has strong connection to the subjects of structure instabilities, morphology control, and materials failures. In this work, we investigate the plastic mechanism developed in the annular crystalline membrane system for adapting to the shrinking space, which is caused by the controllable gradual expansion of the inner boundary. In the process of plastic deformation, we find the continuous generation of dislocations at the inner boundary and their collective migration to the outer boundary; this neat dynamic scenario of dislocation current captures the complicated reorganization process of the particles. We also reveal the characteristic vortex structure arising from the interplay of topological defects and the displacement field. These results may find applications in the precise control of structural instabilities in packings of particulate matter and covalently bonded systems.

The hydraulic conveying of glass beads is studied in a horizontal tube. At low flow rates, plugs can be observed moving across the tube, whereas pseudoplugs can be seen at higher flow rates. A statistical analysis of the plugs' and pseudoplugs' velocities and the plugs' lengths observed is conducted. A transition of the propagation speed distribution is established when the crossing over from a plug to a pseudoplug regime is reached, where the peaked plug velocity distribution turns into a uniform pseudoplug velocity distribution. On the other hand, the statistical distribution of plug lengths exhibits a log-normal mathematical shape. The interpretation of the measured pressure drop evolution with the imposed flow rate by means of an effective viscosity shows an apparent shear-thinning effect coming from the dilution of the granular material. This approach provides a predictive tool for pressure drop calculation in the pseudoplug regime.

Precise modeling of shocks in inertial confinement fusion implosions is critical for obtaining the desired compression in experiments. Shock velocities and postshock conditions are determined by laser-energy deposition, heat conduction, and equations of state. This paper describes experiments at the National Ignition Facility (NIF) [E. M. Campbell and W. J. Hogan, Plasma Phys. Control. Fusion 41, B39 (1999)] where multiple shocks are launched into a cone-in-shell target made of polystyrene, using laser-pulse shapes with two or three pickets and varying on-target intensities. Shocks are diagnosed using the velocity interferometric system for any reflector (VISAR) diagnostic [P. M. Celliers et al., Rev. Sci. Instrum. 75, 4916 (2004)]. Simulated and inferred shock velocities agree well for the range of intensities studied in this work. These directly-driven shock-timing experiments on the NIF provide a good measure of early-time laser-energy coupling. The validated models add to the credibility of direct-drive-ignition designs at the megajoule scale.

For low-density plasmas, the ionization balance can be properly described by the normal Saha equation in the chemical picture. For dense plasmas, however, nonideal effects due to the interactions between the electrons and ions and among the electrons themselves affect the ionization potential depression and the ionization balance. With the increasing of plasma density, the pressure ionization starts to play a more obvious role and competes with the thermal ionization. Based on a local-density temperature-dependent ion-sphere model, we develop a unified and self-consistent theoretical formalism to simultaneously investigate the ionization potential depression, the ionization balance, and the charge states distributions of the dense plasmas. In this work, we choose Al and Au plasmas as examples as Al is a prototype light element and Au is an important heavy element in many research fields such as in the inertial confinement fusion. The nonideal effect of the free electrons in the plasmas is considered by the single-electron effective potential contributed by both the bound electrons of different charge states and the free electrons in the plasmas. For the Al plasmas, we can reconcile the results of two experiments on measuring the ionization potential depression, in which one experiment can be better explained by the Stewart-Pyatt model while the other fits better with the Ecker-Kröll model. For dense Au plasmas, the results show that the double peak structure of the charge state distribution appears to be a common phenomenon. In particular, the calculated ionization balance shows that the two- and three-peak structures can appear simultaneously for denser Au plasmas above $\sim 30\mathrm{g}/{\mathrm{cm}}^3$.

Aging in phase-ordering kinetics of the $d=3$ Ising model following a quench from infinite to zero temperature is studied by means of Monte Carlo simulations. In this model the two-time spin-spin autocorrelator $C_{\text{ag}}$ is expected to obey dynamical scaling and to follow asymptotically a power-law decay with the autocorrelation exponent $\lambda$. Previous work indicated that the lower Fisher-Huse bound of $\lambda \ge d/2=1.5$ is violated in this model. Using much larger systems than previously studied, the instantaneous exponent for $\lambda$ we obtain at late times does not disagree with this bound. By conducting systematic fits to the data of $C_{\text{ag}}$ using different Ansätze for the leading correction term, we find $\lambda =1.58\left(14\right)$, with most of the error attributed to the systematic uncertainty regarding the Ansätze. This result is in contrast to the recent report that below the roughening transition universality might be violated.

We present a multiscale stochastic analysis of foreign exchange rates using the H-theory formalism, which provides a hierarchical intermittency model for the information cascade in the currency market. We examine the distributions of returns and volatilities for the three most traded currency pairs: euro–U.S. dollar, U.S. dollar–Japanese yen, and British pound–U.S. dollar. We find that these markets have a hierarchy of timescales, with larger markets exhibiting more hierarchy levels. We provide a theoretical framework for understanding why the number of levels in the information cascade increases with market size, in analogy with similar behavior for the energy cascade in turbulence as a function of Reynolds number. We briefly argue that using turbulence-like models for financial markets can also provide valuable insights for developing efficient algorithmic trading strategies.

We present evidence for nonlocal contributions to the turbulent energy cascade in magnetohydrodynamic (MHD) plasmas. Therefore, we revisit a well-known result derived directly from the MHD equations, i.e., the Politano and Pouquet law for the transfer of kinetic and magnetic energy in scale. We propose adding a term that accounts for nonlocal transfer and represents the influence of fluctuations from large scales due to the Alfvén effect. Supported by direct numerical simulations of homogeneous and isotropic MHD turbulence, we verify that in some plasma configurations, neglecting the additional nonlocal term might consistently overestimate energy dissipation rates and, thus, the contributions of turbulent energy dissipation potentially affecting solar wind heating, i.e., a central puzzle in space plasma physics that motivates the present work.

The dynamics of open quantum systems connected with several reservoirs attract great attention due to their importance in quantum optics, biology, quantum thermodynamics, transport phenomena, etc. In many problems, the Born approximation is applicable, which implies that the influence of the open quantum system on the reservoirs can be neglected. However, in the case of long-time dynamics or mesoscopic reservoirs, the reverse influence can be crucial. In this paper, we investigate the transient dynamics of several bosonic reservoirs connected through an open quantum system. We use an adiabatic approach to study the temporal dynamics of temperatures of the reservoirs during relaxation to thermodynamic equilibrium. We show that there are various types of temperature dynamics that strongly depend on the values of dissipative rates and initial temperatures. We demonstrate that temperatures of the reservoirs, including the hottest and coldest ones, can exhibit nonmonotonic behavior. Moreover, there are moments of time during which the reservoir with an initially intermediate temperature becomes the hottest or coldest reservoir. The obtained results pave the way for managing energy flows in mesoscale and nanoscale systems.

An operating nuclear reactor is designed to maintain a sustained fission chain reaction in its core, which results from a delicate balance between neutron creations (i.e., fissions) and total absorptions. This balance is associated with random fluctuations that can have two, very different, origins. A distinction must thus be made between low-power noise, whose origin lies in the inherently stochastic nature of neutron interactions with matter, and high-power noise, whose origin lies in the particular thermomechanical constraints associated with the environment in which neutrons propagate. Modeling the behavior of this noisy neutron population with the help of stochastic differential equations, we first show how the Martin-Siggia-Rose-Janssen-De Dominicis (MSRJD) formalism, providing a field theoretical representation of the problem, reveals a convenient and adapted tool for the calculation of observable consequences of neutron noise. In particular, we show how the MSRJD approach is capable of encompassing both types of neutron noises in the same formalism. Emphasizing then on power noise, it is shown how the self-sustained chain reaction developing in a reactor core might be sensitive to noise-induced transitions. Establishing an unprecedented connection between the neutron population evolving in a reactor core and the celebrated Kardar-Parisi-Zhang (KPZ) equation, we indeed find evidence that a noisy reactor core power distribution might be subject to a process analogous to the roughening transition, well-known to occur in systems described by the KPZ equation.

We investigate the nature of the deconfinement transitions in three-dimensional lattice Abelian Higgs models, in which a complex scalar field of integer charge $Q\ge 2$ is minimally coupled with a compact $U\left(1\right)$ gauge field. Their phase diagram presents two phases separated by a transition line where static charges $q$, with $q<Q$, deconfine. We argue that these deconfinement transitions belong to the same universality class as transitions in generic three-dimensional ${\mathbb{Z}}_Q$ gauge models. In particular, they are Ising-like for $Q=2$, of first order for $Q=3$, and belong to the three-dimensional gauge $XY$ universality class for $Q\ge 4$. This general scenario is supported by numerical finite-size scaling analyses of the energy cumulants for $Q=2, Q=4$, and $Q=6$.

Stochastic resetting is a protocol of starting anew, which can be used to facilitate the escape kinetics. We demonstrate that restarting can accelerate the escape kinetics from a finite interval restricted by two absorbing boundaries also in the presence of heavy-tailed, Lévy-type, $\alpha$-stable noise. However, the width of the domain where resetting is beneficial depends on the value of the stability index $\alpha$ determining the power-law decay of the jump length distribution. For heavier (smaller $\alpha$) distributions, the domain becomes narrower in comparison to lighter tails. Additionally, we explore connections between Lévy flights (LFs) and Lévy walks (LWs) in the presence of stochastic resetting. First of all, we show that for Lévy walks, the stochastic resetting can also be beneficial in the domain where the coefficient of variation is smaller than 1. Moreover, we demonstrate that in the domain where LWs are characterized by a finite mean jump duration (length), with the increasing width of the interval, the LWs start to share similarities with LFs under stochastic resetting.

There are various research topics such as stochastic resonance, coherent resonance, and neuroavalanche in excitable systems under external noises. We perform numerical simulation of coupled noisy FitzHugh-Nagumo equations on the square lattice. Excitation waves are generated most efficiently at an intermediate noise strength. The cluster size distributions obey a power-law-like distribution at a certain parameter range. However, we consider that this is not a self-organized critical phenomenon, partly because the exponent of the power law is not constant. We have studied the propagation of excitation waves in the coupled noisy FitzHugh-Nagumo equations with a one-dimensional pacemaker region and found that there is a phase-transition-like phenomenon from the short-range propagation to the whole-system propagation by changing the noise strength $T$. The power-law distribution is observed most clearly near the phase transition of the propagation of excitation waves in the coupled noisy FitzHugh-Nagumo equations without the one-dimensional pacemaker.