We explore the case of a group of random walkers looking for a target randomly located in space, such that the number of walkers is not constant but new ones can join the search, or those that are active can abandon it, with constant rates $r_b$ and $r_d$, respectively. Exact analytical solutions are provided both for the fastest-first-passage time and for the collective time cost required to reach the target, for the exemplifying case of Brownian walkers with $r_d=0$. We prove that even for such a simple situation there exists an optimal rate $r_b$ at which walkers should join the search to minimize the collective search costs. We discuss how these results open a new line to understand the optimal regulation in searches conducted through multiparticle random walks, e.g., in chemical or biological processes.

It is known that two-dimensional two-component fundamental solitons of the semivortex (SV) type, with vorticities $(s_{+},s_{-})=(0,1)$ in their components, are stable ground states (GSs) in the spin-orbit-coupled (SOC) binary Bose-Einstein condensate with the contact self-attraction acting in both components, in spite of the possibility of the critical collapse in the system. However, excited states (ESs) of the SV solitons, with the vorticity set $(s_{+},s_{-})=(S_{+},S_{+}+1)$ and $S_{+}=1,2,3,...$, are unstable in the same system. We construct ESs of SV solitons in the SOC system with opposite signs of the self-interaction in the two components. The main finding is stability of the ES-SV solitons, with the extra vorticity (at least) up to $S_{+}=6$. The threshold value of the norm for the onset of the critical collapse, $N_{\mathrm{thr}}$, in these excited states is higher than the commonly known critical value, $N_c\approx 5.85$, associated with the single-component Townes solitons, $N_{\mathrm{thr}}$ increasing with the growth of $S_{+}$. A velocity interval for stable motion of the GS-SV solitons is found too. The results suggest a solution for the challenging problem of the creation of stable vortex solitons with high topological charges.

We study the $\beta$ model ($\beta$-NG) and the Bayesian Naming Game (BNG) as dynamical systems. By applying linear stability analysis to the dynamical system associated with the $\beta$ model, we demonstrate the existence of a nongeneric bifurcation with a bifurcation point ${\beta }_c=1/3$. As $\beta$ passes through ${\beta }_c$, the stability of isolated fixed points changes, giving rise to a one-dimensional manifold of fixed points. Notably, this attracting invariant manifold forms an arc of an ellipse. In the context of the BNG, we propose modeling the Bayesian learning probabilities $p_A$ and $p_B$ as logistic functions. This modeling approach allows us to establish the existence of fixed points without relying on the overly strong assumption that $p_A=p_B=p$, where $p$ is a constant.

Random matrix theory allows for the deduction of stability criteria for complex systems using only a summary knowledge of the statistics of the interactions between components. As such, results like the well-known elliptical law are applicable in a myriad of different contexts. However, it is often assumed that all components of the complex system in question are statistically equivalent, which is unrealistic in many applications. Here we introduce the concept of a finely structured random matrix. These are random matrices with element-specific statistics, which can be used to model systems in which the individual components are statistically distinct. By supposing that the degree of “fine structure” in the matrix is small, we arrive at a succinct “modified” elliptical law. We demonstrate the direct applicability of our results to the niche and cascade models in theoretical ecology, as well as a model of a neural network, and a directed network with arbitrary degree distribution. The simple closed form of our central results allow us to draw broad qualitative conclusions about the effect of fine structure on stability.

We present a theoretical analysis of phase-separated compartments to facilitate enzymatic chemical reactions. While phase separation can facilitate reactions by increasing local concentration, it can also hinder the mobility of reactants. In particular, we find that the attractive interactions that concentrate reactants within the dense phase can inhibit reactions by lowering the mobility of the reactants. This mobility loss severely limits the potential to enhance reaction rates. Phase separation provides greater benefit in situations where multiple sequential reactions occur and/or high order reactions, provided the enzymes are unsaturated, transport to the condensate is not limiting, and the reactants are mobile. We show that mobility can be maintained if recruitment to the condensed phase is driven by multiple attractive moieties that can bind and release independently. However, the spacers necessary to ensure independence between stickers are prone to entangle with the dense phase scaffold. We find an optimal sticker affinity that balances the need for rapid binding/unbinding kinetics and minimal entanglement. Reaction rates can be accelerated by shrinking the size of the dense phase with a corresponding increase in the number of stickers. Our results showcase the potential capabilities of phase-separated compartments to act as biochemical reaction crucibles within living cells.

The effect of shear flows on the thermal conductivity of polymer melts is investigated using a reversed nonequilibrium molecular-dynamics (RNEMD) method. We extended the original RNEMD method to simultaneously produce spatial gradients of temperature and flow velocity in a single direction. This method enables accurate measurement of the thermal conductivity in the direction transverse to shear flow. The Weissenberg number defined with the shear rate and the relaxation time of the polymer conformation can uniformly differentiate the occurrence of shear rate dependence of the thermal conductivity across different chain lengths. The stress-thermal rule (STR) (i.e., the linear relationship between anisotropic parts of the stress tensor and the thermal conductivity tensor) holds for entangled polymer melts even under shear flows but not for unentangled polymer melts. Furthermore, once entanglements form in polymer chains, the stress-thermal coefficient in the STR remains independent of the polymer chain length. These observations align with the theoretical foundation of the STR, which focuses on energy transmission along the network structure of entangled polymer chains [B. van den Brule, Rheol. Acta 28, 257 (1989)]. However, under driven shear flows, the stress-thermal coefficient is notably smaller than that measured in the literature for a quasiquiescent state without external forces. Although the mechanism of the STR in shear flows has yet to be fully elucidated, our study confirmed the validity of the STR in shear flows. This allows us to use the STR as a constitutive equation for computational thermofluid dynamics of polymer melts, thus having broad engineering applications.

Brownian dynamics simulations are utilized to study segregation phenomena far from thermodynamic equilibrium. In the present study, we expand upon the analysis of binary colloid mixtures and introduce a third particle species to further our understanding of colloidal systems. Gravitationally driven, spherical colloids immersed in an implicit solvent are confined in two-dimensional linear microchannels. The interaction between the colloids is modeled by the Weeks-Chandler-Andersen potential, and the confinement of the colloids is realized by hard walls based on the solution of the Smoluchowski equation in half space. In binary and ternary colloidal systems, a difference in the driving force is achieved by differing colloid sizes but fixed mass density. We observe for both the binary and ternary systems that a driving force difference induces a nonequilibrium phase transition to lanes. For ternary systems, we study the tendency of lane formation to depend on the diameter of the medium-sized colloids. Here we find a sweet spot for lane formation in ternary systems. Furthermore, we study the interaction of two differently sized colloids at the channel walls. Recently we observed that driven large colloids push smaller colloids to the walls. This results in small particle lanes at the walls at early simulation times. In this work we additionally find that thin lanes are unstable and dissolve over very long time frames. Furthermore, we observe a connection between lane formation and the nonuniform distribution of particles along the channel length. This nonuniform distribution occurs either alongside lane formation or in shared lanes (i.e., lanes consisting of two colloid types).

We numerically investigate the effect of an asymmetric periodic obstacle array in a two-dimensional active nematic. We find that activity in conjunction with the asymmetry leads to a ratchet effect or unidirectional flow of the fluid along the asymmetry direction. The directional flow is still present even in the active turbulent phase when the gap between obstacles is sufficiently small. We demonstrate that the dynamics of the topological defects transition from flow mirroring to smectic-like as the gap between obstacles is made smaller, and explain this transition in terms of the pinning of negative winding number defects between obstacles. This also leads to a nonmonotonic ratchet effect magnitude as a function of obstacle size, so that there is an optimal obstacle size for ratcheting at fixed activity.

Non-Newtonian transport properties of a dilute gas of inelastic hard spheres immersed in a molecular gas are determined. We assume that the granular gas is sufficiently rarefied, and hence the state of the molecular gas is not disturbed by the presence of the solid particles. In this situation, one can treat the molecular gas as a bath (or thermostat) of elastic hard spheres at a given temperature. Moreover, in spite of the fact that the number density of grains is quite small, we take into account their inelastic collisions among themselves in its kinetic equation. The system (granular gas plus a bath of elastic hard spheres) is subjected to a simple (or uniform) shear flow. In the low-density regime, the rheological properties of the granular gas are determined by solving the Boltzmann kinetic equation by means of Grad's moment method. These properties turn out to be highly nonlinear functions of the shear rate and the remaining parameters of the system. Our results show that the kinetic granular temperature and the non-Newtonian viscosity present a discontinuous shear thickening effect for sufficiently high values of the mass ratio $m/m_g$ ($m$ and $m_g$ being the mass of grains and gas particles, respectively). This effect becomes more pronounced as the mass ratio $m/m_g$ increases. In particular, in the Brownian limit ($m/m_g\to \infty$) the expressions of the non-Newtonian transport properties derived here are consistent with those previously obtained by considering a coarse-grained approach where the effect of gas phase on grains is through an effective force. Theoretical results are compared against computer simulations, showing an excellent agreement.

The dense active matter exhibits characteristics reminiscent of traditional glassy phenomena, yet the role of rotational inertia in glass dynamics remains elusive. In this study, we investigate the glass dynamics of chiral active particles influenced by rotational inertia. Rotational inertia endows exponential memory to particle orientation, restricting its alteration and amplifying the effective persistence time. At lower spinning frequencies, the diffusion coefficient exhibits a peak function relative to rotational inertia for shorter persistence times, while it steadily increases with rotational inertia for longer persistence times. In the realm of high-frequency spinning, the impact of rotational inertia on diffusion behavior becomes more pronounced, resulting in a nonmonotonic and intricate relationship between the diffusion coefficient and rotational inertia. Consequently, the introduction of rotational inertia significantly alters the glassy dynamics of chiral active particles, allowing for the control over transitions between fluid and glassy states by modulating rotational inertia. Moreover, our findings indicate that at a specific spinning temperature, there exists an optimal spinning frequency at which the diffusion coefficient attains its maximum value.

Dislocation motion under cyclic loading is of great interest from theoretical and practical viewpoints. In this paper, we develop a random walk model for the purpose of evaluating the diffusion coefficient of dislocation under cyclic loading condition. The dislocation behavior was modeled as a series of binomial stochastic processes (one-dimensional random walk), where dislocations are randomly driven by the external load. The probability distribution of dislocation motion and the diffusion coefficient per cycle were analytically derived from the random-walk description as a function of the loading condition and the microscopic material properties. The derived equation was validated by comparing the predicted diffusion coefficient with the molecular dynamics simulation result copper under cyclic deformation. As a result, we confirmed fairly good agreement between the random walk model and the molecular dynamics simulation results.

Shock-driven implosions with 100% deuterium (${\mathrm{D}}_2$) gas fill compared to implosions with 50:50 nitrogen-deuterium (${\mathrm{N}}_2{\mathrm{D}}_2$) gas fill have been performed at the OMEGA laser facility to test the impact of the added mid-$Z$ fill gas on implosion performance. Ion temperature ($T_{\mathrm{ion}}$) as inferred from the width of measured DD-neutron spectra is seen to be $34\%\pm 6\%$ higher for the ${\mathrm{N}}_2{\mathrm{D}}_2$ implosions than for the ${\mathrm{D}}_2$-only case, while the DD-neutron yield from the ${\mathrm{D}}_2$-only implosion is $7.2\pm 0.5$ times higher than from the ${\mathrm{N}}_2{\mathrm{D}}_2$ gas fill. The $T_{\mathrm{ion}}$ enhancement for ${\mathrm{N}}_2{\mathrm{D}}_2$ is observed in spite of the higher $Z$, which might be expected to lead to higher radiative loss, and higher shock strength for the ${\mathrm{D}}_2$-only versus ${\mathrm{N}}_2{\mathrm{D}}_2$ implosions due to lower mass, and is understood in terms of increased shock heating of N compared to D, heat transfer from N to D prior to burn, and limited amount of ion-electron-equilibration-mediated additional radiative loss due to the added higher-$Z$ material. This picture is supported by interspecies equilibration timescales for these implosions, constrained by experimental observables. The one-dimensional (1D) kinetic Vlasov-Fokker-Planck code ifp and the radiation hydrodynamic simulation codes hyades (1D) and xrage [1D, two-dimensional (2D)] are brought to bear to understand the observed yield ratio. Comparing measurements and simulations, the yield loss in the ${\mathrm{N}}_2{\mathrm{D}}_2$ implosions relative to the pure ${\mathrm{D}}_2$-fill implosion is determined to result from the reduced amount of ${\mathrm{D}}_2$ in the fill (fourfold effect on yield) combined with a lower fraction of the ${\mathrm{D}}_2$ fuel being hot enough to burn in the ${\mathrm{N}}_2{\mathrm{D}}_2$ case. The experimental yield and $T_{\mathrm{ion}}$ ratio observations are relatively well matched by the kinetic simulations, which suggest interspecies diffusion is responsible for the lower fraction of hot ${\mathrm{D}}_2$ in the ${\mathrm{N}}_2{\mathrm{D}}_2$ relative to the ${\mathrm{D}}_2$-only case. The simulated absolute yields are higher than measured; a comparison of 1D versus 2D xrage simulations suggest that this can be explained by dimensional effects. The hydrodynamic simulations suggest that radiative losses primarily impact the implosion edges, with ion-electron equilibration times being too long in the implosion cores. The observations of increased $T_{\mathrm{ion}}$ and limited additional yield loss (on top of the fourfold expected from the difference in D content) for the ${\mathrm{N}}_2{\mathrm{D}}_2$ versus ${\mathrm{D}}_2$-only fill suggest it is feasible to develop the platform for studying CNO-cycle-relevant nuclear reactions in a plasma environment.

We study melting in a two-dimensional system of classical particles with Gaussian-core interactions in disordered environments. The pure system validates the conventional two-step melting with a hexatic phase intervening between the solid and the liquid. This picture is modified in the presence of pinning impurities. A random distribution of pinning centers forces a hexaticlike low-temperature phase that transits into a liquid at a single melting temperature $T_{\mathrm{m}}^{\mathrm{RP}}$. In contrast, pinning centers located at randomly chosen sites of a perfect crystal anchor a solid at low temperatures which undergoes a direct transition to the liquid at $T_{\mathrm{m}}^{\mathrm{CP}}$. Thus, the two-step melting is lost in either case of disorder. We discuss the characteristics of melting depending on the nature of the impurities.

It is demonstrated that quantum systems classically exhibiting strong and homogeneous chaos in a bounded region of the phase space can induce a global quantum diffusion. As an ideal model system, a small quantum chaos with finite Hilbert space dimension $N$ weakly coupled with $M$ additional degrees of freedom which is approximated by linear systems is proposed. By twinning the system the diffusion process in the additional modes can be numerically investigated without taking the unbounded diffusion space into account explicitly. Even though $N$ is not very large, diffusion occurs in the additional modes as the coupling strength increases if $M\ge 3$. If $N$ is large enough, a definite quantum transition to diffusion takes place through a critical subdiffusion characterized by an anomalous diffusion exponent.

We report on the unexpected precessional motion that occurs when a small, rigid ring is rotated on a vertical smooth rod. Using high-speed imaging, two distinct regimes of motion are observed experimentally. (i) Oscillatory motion when the ring has a single contact with the rod (transient motion). (ii) Steady state (terminal velocity) when the ring has two instantaneously stationary contacts with the rod. These two regimes and the transition between them are analyzed through qualitative, analytical, and numerical means. One key feature of the steady-state motion is energy dissipation through rolling resistance and air drag. Our model predicts that steady-state motion occurs only when certain geometric conditions (the ratio of the radii of a ring and a rod, and the tilt angle) are fulfilled. As the steady-state regime is often found to be the preferred final state, the dynamic stability of the single contact configuration is investigated through bifurcation, phase-space, and linear stability analyses. Our numerical simulation accounts for the stick-slip transition at the contact point, and is in good agreement with experimental data across all experiments. Our findings have potential applications in various engineering fields, such as studying the vibrational detachment of threaded fasteners and developing gravity-driven centrifuges.

We consider random hyperbolic graphs in hyperbolic spaces of any dimension $d+1\ge 2$. We present a rescaling of model parameters that casts the random hyperbolic graph model of any dimension to a unified mathematical framework, leaving the degree distribution invariant with respect to the dimension. Unlike the degree distribution, clustering does depend on the dimension, decreasing to 0 at $d\to \infty$. We analyze all of the other limiting regimes of the model, and we release a software package that generates random hyperbolic graphs and their limits in hyperbolic spaces of any dimension.

One of the main challenges in developing high-performance quantum batteries is the self-discharging process, where energy is dissipated from a quantum battery into the environment. In this work, we investigate the influence of non-Markovian noises on the performance of a quantum battery. Our results demonstrate that adding auxiliary qubits to a quantum battery system can effectively suppress the self-discharging process, leading to an improvement in both the steady-state energy and extractable work. We reveal that the physical mechanism inhibiting the self-discharging process is the formation of system-environment bound states, rather than an increase in non-Markovianity. Our results could be of both theoretical and experimental interest in exploring the ability of quantum batteries to maintain long stored energy in the environment.

Kagome spin ice is one of the canonical examples of highly frustrated magnets. The effective magnetic degrees of freedom in kagome spin ice are Ising spins residing on a two-dimensional network of corner-sharing triangles. Due to strong geometrical frustration, nearest-neighbor antiferromagnetic interactions on the kagome lattice give rise to a macroscopic number of degenerate classical ground states characterized by ice rules. Elementary excitations at low temperatures are defect-triangles that violate the ice rules and carry an additional net magnetic charge relative to the background. We perform large-scale Glauber dynamics simulations to study the nonequilibrium dynamics of kagome ice under slow cooling. We show that the density of residual charge defects exhibits a power-law dependence on the quench rate for the class of algebraic cooling protocols. The numerical results are well captured by the rate equation for the charge defects based on the reaction kinetics theory. As the relaxation time of the kagome ice phase remains finite, there is no dynamical freezing as in the Kibble-Zurek scenario. Instead, we show that the power-law behavior originates from a thermal excitation that decays algebraically with time at the late stage of the cooling schedule. Similarities and differences in quench dynamics of other spin ice systems are also discussed.

The conformational and dynamical properties of active ring polymers are studied by numerical simulations. The two-dimensionally confined polymer is modeled as a closed bead-spring chain, driven by tangential forces, put in contact with a heat bath described by the Brownian multiparticle collision dynamics. Both phantom polymers and chains comprising excluded volume interactions are considered for different bending rigidities. The size and shape are found to be dependent on persistence length, driving force, and bead mutual exclusion. The lack of excluded volume interactions is responsible for a shrinkage of active rings when increasing driving force in the flexible limit, while the presence induces a moderate swelling of chains. The internal dynamics of flexible phantom active rings shows activity-enhanced diffusive behavior at large activity values while, in the case of self-avoiding active chains, it is characterized by active ballistic motion not depending on stiffness. The long-time dynamics of active rings is marked by rotational motion whose period scales as the inverse of the applied tangential force, irrespective of persistence length and beads' self-exclusion.

We consider discrete models of kinetic rough interfaces that exhibit space-time scale invariance in height-height correlation. We use the generic scaling theory of Ramasco et al. [Phys. Rev. Lett. 84, 2199 (2000)] to confirm that the dynamical structure factor of the height profile can uniquely characterize the underlying dynamics. We apply both finite-size and finite-time scaling methods that systematically allow an estimation of the critical exponents and the scaling functions, eventually establishing the universality class accurately. The finite-size scaling analysis offers an alternative way to characterize the anomalous rough interfaces. As an illustration, we investigate a class of self-organized interface models in random media with extremal dynamics. The isotropic version shows a faceted pattern and belongs to the same universality class (as shown numerically) as the Sneppen model (version A). We also examine an anisotropic version of the Sneppen model and suggest that the model belongs to the universality class of the tensionless Kardar-Parisi-Zhang (tKPZ) equation in one dimension.