We discuss how the synchronization in the Kuramoto model can be treated in terms of information geometry. We argue that the Fisher information is sensitive to synchronization transition; specifically, components of the Fisher metric diverge at the critical point. Our approach is based on the recently proposed relation between the Kuramoto model and geodesics in hyperbolic space.

Coupled nonlinear oscillators are ubiquitous in dynamical studies. A wealth of behaviors have been found mostly for globally coupled systems. From a complexity perspective, less studied have been systems with local coupling, which is the subject of this contribution. The phase approximation is used, as weak coupling is assumed. In particular, the so-called needle region, in parameter space, for Adler-type oscillators with nearest neighbors coupling is carefully characterized. The reason for this emphasis is that, in the border of this region to the surrounding chaotic one, computation enhancement at the edge of chaos has been reported. The present study shows that different behaviors within the needle region can be found and a smooth change of dynamics could be identified. Entropic measures further emphasize the region's heterogeneous nature with interesting features, as seen in the spatiotemporal diagrams. The occurrence of wave-like patterns in the spatiotemporal diagrams points to nontrivial correlations in both dimensions. The wave patterns change as the control parameters change without exiting the needle region. Spatial correlation is only achieved locally at the onset of chaos, with different clusters of oscillators behaving coherently while disordered boundaries appear between them.

Ultracold atoms confined to optical lattices provide a platform for simulation of phenomena not readily accessible in condensed matter and chemical systems. One area of growing interest is the mechanism by which isolated condensed matter systems can thermalize. The mechanism for thermalization of quantum systems has been directly linked to a transition to chaos in their classical counterpart. Here we show that the broken spatial symmetries of the honeycomb optical lattice leads to a transition to chaos in the single-particle dynamics which, in turn, causes mixing of the energy bands of the quantum honeycomb lattice. For systems with single-particle chaos, “soft” interactions between atoms can cause the system to thermalize (achieve a Fermi-Dirac distribution for fermions or a Bose-Einstein distribution for bosons).

Motivated by previous results showing that the addition of a linear dispersive term to the two-dimensional Kuramoto-Sivashinsky equation has a dramatic effect on the pattern formation, we study the Swift-Hohenberg equation with an added linear dispersive term, the dispersive Swift-Hohenberg equation (DSHE). The DSHE produces stripe patterns with spatially extended defects that we call seams. A seam is defined to be a dislocation that is smeared out along a line segment that is obliquely oriented relative to an axis of reflectional symmetry. In contrast to the dispersive Kuramoto-Sivashinsky equation, the DSHE has a narrow band of unstable wavelengths close to an instability threshold. This allows for analytical progress to be made. We show that the amplitude equation for the DSHE close to threshold is a special case of the anisotropic complex Ginzburg-Landau equation (ACGLE) and that seams in the DSHE correspond to spiral waves in the ACGLE. Seam defects and the corresponding spiral waves tend to organize themselves into chains, and we obtain formulas for the velocity of the spiral wave cores and for the spacing between them. In the limit of strong dispersion, a perturbative analysis yields a relationship between the amplitude and wavelength of a stripe pattern and its propagation velocity. Numerical integrations of the ACGLE and the DSHE confirm these analytical results.

Recurrently coupled oscillators that are sufficiently heterogeneous and/or randomly coupled can show an asynchronous activity in which there are no significant correlations among the units of the network. The asynchronous state can nevertheless exhibit a rich temporal correlation statistics that is generally difficult to capture theoretically. For randomly coupled rotator networks, it is possible to derive differential equations that determine the autocorrelation functions of the network noise and of the single elements in the network. So far, the theory has been restricted to statistically homogeneous networks, making it difficult to apply this framework to real-world networks, which are structured with respect to the properties of the single units and their connectivity. A particularly striking case are neural networks for which one has to distinguish between excitatory and inhibitory neurons, which drive their target neurons towards or away from the firing threshold. To take into account network structures like that, here we extend the theory for rotator networks to the case of multiple populations. Specifically, we derive a system of differential equations that govern the self-consistent autocorrelation functions of the network fluctuations in the respective populations. We then apply this general theory to the special but important case of recurrent networks of excitatory and inhibitory units in the balanced case and compare our theory to numerical simulations. We inspect the effect of the network structure on the noise statistics by comparing our results to the case of an equivalent homogeneous network devoid of internal structure. Our results show that structured connectivity and heterogeneity of the oscillator type can both enhance or reduce the overall strength of the generated network noise and shape its temporal correlations.

The contact stiffness of an aluminum bead confined between two slabs diminishes upon mechanical conditioning, and then recovers as $log\left(t\right)$ after the conditioning ceases. Here that structure is evaluated for its response to transient heating and cooling, with and without accompanying conditioning vibrations. It is found that, under heating or cooling alone, stiffness changes are mostly consistent with temperature-dependent material moduli; there is little or no slow dynamics. Hybrid tests in which vibration conditioning is followed by heating or cooling lead to recoveries that begin as $log\left(t\right)$ and then become more complex. On subtracting the known response to heating or cooling alone we discern the influence of higher or lower temperatures on slow dynamic recovery from vibrations. It is found that heating accelerates the initial $log\left(t\right)$ recovery, but by an amount more than predicted by an Arrhenius model of thermally activated barrier penetrations. Transient cooling has no discernible effect, in contrast to the Arrhenius prediction that it slows recovery.

A multiscale steady discrete unified gas kinetic scheme with macroscopic coarse mesh acceleration [accelerated steady discrete unified gas kinetic scheme (SDUGKS)] is proposed to improve the convergence of the original SDUGKS for an optically thick system in solving the multigroup neutron Boltzmann transport equation (NBTE) to analyze the distribution of fission energy in the reactor core. In the accelerated SDUGKS, by solving the coarse mesh macroscopic governing equations (MGEs) derived from the moment equations of the NBTE, the numerical solutions of the NBTE on fine meshes at the mesoscopic level can be rapidly obtained from the prolongation of the coarse mesh solutions of the MGE. Furthermore, the use of the coarse mesh can greatly reduce the computational variables and improve the computational efficiency of the MGE. The biconjugate gradient stabilized Krylov subspace method with the modified incomplete LU preconditioner and the lower-upper symmetric-Gauss-Seidel sweeping method are implemented to solve the discrete systems of the macroscopic coarse mesh acceleration model and mesoscopic SDUGKS to further improve the numerical efficiency. Numerical solutions validate good numerical accuracy and high acceleration efficiency of the proposed accelerated SDUGKS for the complicated multiscale neutron transport problems.

Single-file diffusion refers to the Brownian motion in narrow channels where particles cannot pass each other. In such processes, the diffusion of a tagged particle is typically normal at short times and becomes subdiffusive at long times. For hard-sphere interparticle interaction, the time-dependent mean squared displacement of a tracer is well understood. Here we develop a scaling theory for adhesive particles. It provides a full description of the time-dependent diffusive behavior with a scaling function that depends on an effective strength of adhesive interaction. Particle clustering induced by the adhesive interaction slows down the diffusion at short times, while it enhances subdiffusion at long times. The enhancement effect can be quantified in measurements irrespective of how tagged particles are injected into the system. Combined effects of pore structure and particle adhesiveness should speed up translocation of molecules through narrow pores.

In this work, we have studied the Ising model with one- and two-spin flip competing dynamics on a two-dimensional additive small-world network (A-SWN). The system model consists of an $L\times L$ square lattice where each site of the lattice is occupied by a spin variable that interacts with the nearest-neighbor spins and it has a certain probability $p$ of being additionally connected at random to one of its farther neighbors. The dynamics present in the system can be defined by the probability $q$ of being in contact with a heat bath at a given temperature $T$ and, at the same time, with a probability of $(1-q)$ the system is subjected to an external flux of energy into the system. The contact with the heat bath is simulated by one-spin flip according to the Metropolis prescription, while the input of energy is mimicked by the two-spin flip process, involving a simultaneous flipping of a pair of neighboring spins. We have employed Monte Carlo simulations to obtain the thermodynamic quantities of the system, such as the total $m_L^{\mathrm{F}}$ and staggered $m_L^{\mathrm{AF}}$ magnetizations per spin, the susceptibility ${\chi }_L$, and the reduced fourth-order Binder cumulant $U_L$. We have built the phase diagram for the stationary states of the model in the plane $T$ versus $q$, showing the existence of two continuous transition lines for each value of $p$: one line between the ferromagnetic F and paramagnetic P phases and the other line between the P and antiferromagnetic AF phases. Therefore, we have shown that the phase diagram topology changes when $p$ increases. Using the finite-size scaling analysis, we also obtained the critical exponents for the system, where, varying the parameter $p$, we have observed a different universality class from the Ising model in the regular square lattice to the A-SWN.

Via hydrodynamics-preserving molecular dynamics simulations we study growth phenomena in a phase-separating symmetric binary mixture model. We quench high-temperature homogeneous configurations to state points inside the miscibility gap, for various mixture compositions. For compositions at the symmetric or critical value we capture the rapid linear viscous hydrodynamic growth due to advective transport of material through tubelike interconnected domains. For state points very close to any of the branches of the coexistence curve, the growth in the system, following nucleation of disconnected droplets of the minority species, occurs via a coalescence mechanism. Using state-of-the-art techniques, we have identified that these droplets, between collisions, exhibit diffusive motion. The value of the exponent for the power-law growth, related to this diffusive coalescence mechanism, has been estimated. While the exponent nicely agrees with that for the growth via the well-known Lifshitz-Slyozov particle diffusion mechanism, the amplitude is stronger. For the intermediate compositions we observe initial rapid growth that matches the expectations for viscous or inertial hydrodynamic pictures. However, at later times these types of growth cross over to the exponent that is decided by the diffusive coalescence mechanism.

We study the mean first-passage time of a one-dimensional active fluctuating membrane that is stochastically returned to the same flat initial condition at a finite rate. We start with a Fokker-Planck equation to describe the evolution of the membrane coupled with an Ornstein-Uhlenbeck type of active noise. Using the method of characteristics, we solve the equation and obtain the joint distribution of the membrane height and active noise. In order to obtain the mean first-passage time (MFPT), we further obtain a relation between the MFPT and a propagator that includes stochastic resetting. The derived relation is then used to calculate it analytically. Our studies show that the MFPT increases with a larger resetting rate and decreases with a smaller rate, i.e., there is an optimal resetting rate. We compare the results in terms of MFPT of the membrane with active and thermal noises for different membrane properties. The optimal resetting rate is much smaller with active noise compared to thermal. When the resetting rate is much lower than the optimal rate, we demonstrate how the MFPT scales with resetting rates, distance to the target, and the properties of the membranes.

The Drazin inverse of the Liouvillian superoperator provides a solution to determine the dynamics of a time-dependent system governed by the Markovian master equation. Under the condition of slow driving, the perturbation expansion of the density operator of the system in powers of time can be derived. As an application, a finite-time cycle model of the quantum refrigerator driven by a time-dependent external field is established. The method of the Lagrange multiplier is adopted as a strategy to find the optimal cooling performance. The figure of merit given by the product of the coefficient of performance and the cooling rate is taken as a new objective function, and, consequently, the optimally operating state of the refrigerator is revealed. The effects of the frequency exponent determining dissipation characteristics on the optimal performance of the refrigerator are discussed systemically. The results obtained show that the adjacent areas of the state of the maximum figure of merit are the best operation region of low-dissipative quantum refrigerators.

The disorder-order transitions of layering and intralayer structural orders of three-dimensional Yukawa liquids, under the enhanced confinement effect with decreasing normal distance $z$ to the confinement boundary, is investigated numerically. The liquid between the two flat boundaries is segmented into many slabs parallel to the boundary, with the same slab width as the layer width. In each slab, particle sites are binarized into sites with layering order (LOSs)/ layering disorder (LDSs) and with intralayer structural order (SOSs)/disorder (SDSs). It is found that with decreasing $z$, a small fraction of LOSs starts to heterogeneously emerge in the form of small clusters in the slab, followed by the emergence of the large percolating LOS clusters spanning over the system. The smooth rapid rise of the fraction of LOSs from small values followed by their gradual saturations, and the scaling behavior of multiscale LOS clustering, are similar to those of the nonequilibrium systems governed by the percolation theory. The disorder-order transition of intraslab structural ordering also exhibits a similar generic behavior as that of layering with the same transition slab number. The spatial fluctuations of local layering order and local intralayer structural order are uncorrelated in the bulk liquid and the outmost layer next to the boundary. Approaching the percolating transition slab, their correlation gradually increases to the maximum.

We uncover the emergence of distinct sets of multistable chimera states in addition to chimera death and synchronized states in a smallest population of three globally coupled oscillators with mean-field diffusive coupling. Sequence of torus bifurcations result in the manifestation of distinct periodic orbits as a function of the coupling strength, which in turn result in the genesis of distinct chimera states constituted by two synchronized oscillators coexisting with an asynchronous oscillator. Two subsequent Hopf bifurcations result in homogeneous and inhomogeneous steady states resulting in desynchronized steady states and chimera death state among the coupled oscillators. The periodic orbits and the steady states lose their stability via a sequence of saddle-loop and saddle-node bifurcations finally resulting in a stable synchronized state. We have generalized these results to $N$ coupled oscillators and also deduced the variational equations corresponding to the perturbation transverse to the synchronization manifold and corroborated the synchronized state in the two-parameter phase diagrams using its largest eigenvalue. Chimera states in three coupled oscillators emerge as a solitary state in $N$ coupled oscillator ensemble.

We study numerically the vortex dynamics and vortex-lattice formation in a rotating density-dependent Bose-Einstein condensate (BEC), characterized by the presence of nonlinear rotation. By varying the strength of nonlinear rotation in density-dependent BECs, we calculate the critical frequency, ${\mathrm{\Omega }}_{\text{cr}}$, for vortex nucleation both in adiabatic and sudden external trap rotations. The nonlinear rotation modifies the extent of deformation experienced by the BEC due to the trap and shifts the ${\mathrm{\Omega }}_{\text{cr}}$ values for vortex nucleation. The critical frequencies, and thereby the transition to vortex-lattices in an adiabatic rotation ramp, depend on conventional $s$-wave scattering lengths through the strength of nonlinear rotation, $\mathit{C}$, such that ${\mathrm{\Omega }}_{\text{cr}}(\mathit{C}>0)<{\mathrm{\Omega }}_{\text{cr}}(\mathit{C}=0)<{\mathrm{\Omega }}_{\text{cr}}(\mathit{C}<0)$. In an analogous manner, the critical ellipticity (${\epsilon }_{\text{cr}}$) for vortex nucleation during an adiabatic introduction of trap ellipticity ($\epsilon$) depends on the nature of nonlinear rotation besides trap rotation frequency. The nonlinear rotation additionally affects the vortex-vortex interactions and the motion of the vortices through the condensate by altering the strength of Magnus force on them. The combined result of these nonlinear effects is the formation of the non-Abrikosov vortex-lattices and ring-vortex arrangements in the density-dependent BECs.

By interpreting a temporal network as a trajectory of a latent graph dynamical system, we introduce the concept of dynamical instability of a temporal network and construct a measure to estimate the network maximum Lyapunov exponent (nMLE) of a temporal network trajectory. Extending conventional algorithmic methods from nonlinear time-series analysis to networks, we show how to quantify sensitive dependence on initial conditions and estimate the nMLE directly from a single network trajectory. We validate our method for a range of synthetic generative network models displaying low- and high-dimensional chaos and finally discuss potential applications.

We explore the mechanics and damage of slide-ring gels by developing a discrete model for the mechanics of chain-ring polymer systems that accounts for both crosslink motion and internal chain sliding. The proposed framework utilizes an extendable Langevin chain model to describe the constitutive behavior of polymer chains undergoing large deformation and includes a rupture criterion to innately capture damage. Similarly, crosslinked rings are described as large molecules that also store enthalpic energy during deformation and thus have their own rupture criterion. Using this formalism, we show that the realized mode of damage in a slide-ring unit is a function of the loading rate, distribution of segments, and inclusion ratio (number of rings per chain). After analyzing an ensemble of representative units under different loading conditions, we find that failure is driven by damage to crosslinked rings at slow loading rates, but polymer chain scission at fast loading rates. Our results indicate that increasing the strength of the crosslinked rings may improve the toughness of the material.

In battery modeling, the electrode is discretized at the macroscopic scale with a single representative particle in each volume. This lacks the accurate physics to describe interparticle interactions in electrodes. To remedy this, we formulate a model that describes the evolution of degradation of a population of battery active material particles using ideas in population genetics of fitness evolution, where the state of a system depends on the health of each particle that contributes to the system. With the fitness formulation, the model incorporates effects of particle size and heterogeneous degradation effects which accumulate in the particles as the battery is cycled, accounting for different active material degradation mechanisms. At the particle scale, degradation progresses nonuniformly across the population of active particles, observed from the autocatalytic relationship between fitness and degradation. Electrode-level degradation is formed from various contributions of the particle-level degradation, especially from smaller particles. It is shown that specific mechanisms of particle-level degradation can be associated with characteristic signatures in the capacity-loss and voltage profiles. Conversely, certain features in the electrode-level phenomena can also provide insight into the relative importance of different particle-level degradation mechanisms.

Rod-shaped particles embedded in certain matrices have been reported to exhibit an increase in their center of mass diffusivity upon increasing the matrix density. This increase has been considered to be caused by a kinetic constraint in analogy with tube models. We investigate a mobile rodlike particle in a sea of immobile point obstacles using a kinetic Monte Carlo scheme equipped with a Markovian process, that generates gaslike collision statistics, so that such kinetic constraints do essentially not exist. Even in such a system, provided the particle's aspect ratio exceeds a threshold value of about 24, the unusual increase in the rod diffusivity emerges. This result implies that the kinetic constraint is not a necessary condition for the increase in the diffusivity.

We study size- and charge-asymmetric oppositely charged colloids driven by an external electric field. The large particles are connected by harmonic springs, forming a hexagonal-lattice network, while the small particles are free of bonds and exhibit fluidlike motion. We show that this model exhibits a cluster formation pattern when the external driving force exceeds a critical value. The clustering is accompanied with stable wave packets in vibrational motions of the large particles.