We study a three-wave truncation of the high-order nonlinear Schrödinger equation for deep-water waves (also named Dysthe equation). We validate the model by comparing it to numerical simulation; we distinguish the impact of the different fourth-order terms and classify the solutions according to their topology. This allows us to properly define the temporary spectral upshift occurring in the nonlinear stage of Benjamin-Feir instability and provides a tool for studying further generalizations of this model.

We consider optimization of the linear stability of synchronized states between a pair of weakly coupled limit-cycle oscillators with cross coupling, where different components of state variables of the oscillators are allowed to interact. On the basis of the phase reduction theory, we derive the coupling matrix between different components of the oscillator states that maximizes the linear stability of the synchronized state under given constraints on the overall coupling intensity and the stationary phase difference. The improvement in the linear stability is illustrated by using several types of limit-cycle oscillators as examples.

Optimization of the stability of synchronized states between a pair of symmetrically coupled reaction-diffusion systems exhibiting rhythmic spatiotemporal patterns is studied in the framework of the phase reduction theory. The optimal linear filter that maximizes the linear stability of the in-phase synchronized state is derived for the case in which the two systems are nonlocally coupled. The optimal nonlinear interaction function that theoretically gives the largest linear stability of the in-phase synchronized state is also derived. The theory is illustrated by using typical rhythmic patterns in FitzHugh-Nagumo systems as examples.

By introducing evolving disorder in the binary kagome ribbons, we study the establishment of diffusive spreading of flat band states characterized by diffractionless propagation in regular periodic ribbons. Our numerical analysis relies on controlling strength and rate of change of disorder during light propagation while tailoring binarism of the kagome ribbons in order to isolate the flat band with the gap from the rest of the ribbon's eigenvalue spectrum and study systematically its influence on diffusion. We show that the flat band plays a dominant role in the establishment of the diffusion for a given strength and rate of change of disorder, whereas the rest of the ribbon's eigenvalue spectrum induces only quantitative differences in the light spreading regimes. Due to the universality of studied phenomena, our findings may be of interest in various disordered physical systems with flat spectral bands, ranging from photonics to ultracold matter systems and plasmonics.

Swimming bacteria exhibit a repertoire of motility patterns, in which persistent motion is interrupted by turning events. What are the statistical properties of such random walks? If some particular instances have long been studied, the general case where turning times do not follow a Poisson process has remained unsolved. We present a generic extension of the continuous time random walks formalism relying on operators and noncommutative calculus. The approach is first applied to a unimodal model of bacterial motion. We examine the existence of a minimum in velocity correlation function and discuss the maximum of diffusivity at an optimal value of rotational diffusion. The model is then extended to bimodal patterns and includes as particular cases all swimming strategies: run-and-tumble, run-stop, run-reverse and run-reverse-flick. We characterize their velocity correlation functions and investigate how bimodality affects diffusivity. Finally, the wider applicability of the method is illustrated by considering curved trajectories and Lévy walks. Our results are relevant for intermittent motion of living beings, be they swimming micro-organisms or crawling cells.

For the relativistic Kelvin-Helmholtz instability (KHI), which occurs at shear interfaces between two plasma streams, we report results on the polarized radiation over all observation directions and frequencies emitted by the plasma electrons from ab initio kinetic simulations. We find the polarization of the radiation to provide a clear signature for distinguishing the linear phase of the KHI from its other phases. During the linear phase, we predict the growth rate of the KHI radiation power to match the growth rate of the KHI to a high degree. Our predictions are based on a model of the vortex dynamics, which describes the electron motion in the vicinity of the shear interface between the two streams. Albeit the complex and turbulent dynamics happening in the shear region, we find excellent agreement between our model and large-scale particle-in-cell simulations. Our findings pave the way for identifying the KHI linear regime and for measuring its growth rate in astrophysical jets observable on earth as well as in laboratory plasmas.

Time delay in general leads to instability in some systems, while specific feedback with delay can control fluctuated motion in nonlinear deterministic systems to a stable state. In this paper, we consider a stochastic process, i.e., a random walk, and observe its diffusion phenomenon with time-delayed feedback. As a result, the diffusion coefficient decreases with increasing delay time. We analytically illustrate this suppression of diffusion by using stochastic delay differential equations and justify the feasibility of this suppression by applying time-delayed feedback to a molecular dynamics model.

Data-based and model-free accurate identification of intrinsic time delays and directional interactions is an extremely challenging problem in complex dynamical systems and their networks reconstruction. A model-free method with new scores is proposed to be generally capable of detecting single, multiple, and distributed time delays. The method is applicable not only to mutually interacting dynamical variables but also to self-interacting variables in a time-delayed feedback loop. Validation of the method is carried out using physical, biological, and ecological models and real data sets. Especially, applying the method to air pollution data and hospital admission records of cardiovascular diseases in Hong Kong reveals the major air pollutants as a cause of the diseases and, more importantly, it uncovers a hidden time delay (about 30–40 days) in the causal influence that previous studies failed to detect. The proposed method is expected to be universally applicable to ascertaining and quantifying subtle interactions (e.g., causation) in complex systems arising from a broad range of disciplines.

We study the stability of asymptotic states displayed by a complex neural network. We focus on the loss of stability of a stationary state of networks using recurrence quantifiers as tools to diagnose local and global stabilities as well as the multistability of a coupled neural network. Numerical simulations of a neural network composed of 1024 neurons in a small-world connection scheme are performed using the model of Braun et al. [Int. J. Bifurcation Chaos 08, 881 (1998)], which is a modified model from the Hodgkin-Huxley model [J. Phys. 117, 500 (1952)]. To validate the analyses, the results are compared with those produced by Kuramoto's order parameter [Chemical Oscillations, Waves, and Turbulence (Springer-Verlag, Berlin Heidelberg, 1984)]. We show that recurrence tools making use of just integrated signals provided by the networks, such as local field potential (LFP) (LFP signals) or mean field values bring new results on the understanding of neural behavior occurring before the synchronization states. In particular we show the occurrence of different stationary and nonstationarity asymptotic states.

The propagation of a beneficial mutation in a spatially extended population is usually studied using the phenomenological stochastic Fisher-Kolmogorov-Petrovsky-Piscounov (SFKPP) equation. We derive here an individual-based, stochastic model founded on the spatial Moran process where fluctuations are treated exactly. The mean-field approximation of this model leads to an equation that is different from the phenomenological FKPP equation. At small selection pressure, the front behavior can be mapped into a Brownian motion with drift, the properties of which can be derived from the microscopic parameters of the Moran model. Finally, we generalize the model to take into account dispersal kernels beyond migration to nearest neighbors. We show how the effective population size (which controls the noise amplitude) and the diffusion coefficient can both be computed from the dispersal kernel.

This paper presents a numerical investigation of the dynamics of pinch-off in liquid drops and jets during injection of a liquid through an orifice into another fluid. The current study is carried out by solving axisymmetric Navier-Stokes equations and the interface is captured using a coupled level-set and volume-of-fluid approach. The delicate interplay of inertia and viscous effects plays a crucial role in deciding the dynamics of the formation as well as breakup of liquid drops and jets. In the dripping regime, the growth and breakup rate of a drop are studied and quantified by corroborating with theoretical predictions. During the growth stage of the drops, a self-similar behavior of the drop profile is identified over a relatively short duration of time. The viscosity of the drop liquid shows substantial influence on the thinning behavior of a liquid neck and a transition is observed from an inertia dominated regime to an inertia-viscous regime beyond a critical minimum value of the neck radius. The phenomenon of interface overturning is fundamentally related to the magnitude of drop viscosity. The variation of overturning angle as a function of drop viscosity is computed and a critical value of Ohnesorge number is obtained beyond which overturning ceases. Increasing the inertia of drop liquid transforms the system from a periodically dripping regime to a quasiperiodic regime and finally it culminates into an elongated liquid jet. Another interesting transition from dripping to jetting regime is demonstrated by varying the viscosity of the ambient medium. The breakup of jets in Rayleigh mode is explored and the breakup length obtained from our computations shows excellent agreement with the theoretical predictions owing to Rayleigh's analysis. The ambient medium is entrained as the jet moves downstream with the creation of a vortical structure just outside the jet signifying increased participation of the ambient medium in the dynamics of jet breakup at higher inflow rates.

Gas-surface interaction is studied by the molecular dynamics method to investigate qualitatively characteristics of accommodation coefficients. A large number of trajectories of gas molecules colliding to and scattering from a surface are statistically analyzed to calculate the energy (thermal) accommodation coefficient (EAC) and the tangential momentum accommodation coefficient (TMAC). Considering experimental measurements of the accommodation coefficients, the incident velocities are stochastically sampled to represent a bulk condition. The accommodation coefficients for noble gases show qualitative coincidence with experimental values. To investigate characteristics of these accommodation coefficients in detail, the gas-surface interaction is parametrically studied by varying the molecular mass of gas, the gas-surface interaction strength, and the molecular size of gas, one by one. EAC increases with increasing every parameter, while TMAC increases with increasing the interaction strength, but decreases with increasing the molecular mass and the molecular size. Thus, contradictory results in experimentally measured TMAC for noble gases could result from the difference between the surface conditions employed in the measurements in the balance among the effective parameters of molecular mass, interaction strength, and molecular size, due to surface roughness and/or adsorbed molecules. The accommodation coefficients for a thermo-fluid dynamics field with a temperature difference between gas and surface and a bulk flow at the same time are also investigated.

In conventional thermodynamics, it is widely acknowledged that the realization of an isothermal process for a system requires a quasistatic controlling protocol. Here we propose and design a strategy to realize a finite-rate isothermal transition from an equilibrium state to another one at the same temperature, which is named the “shortcut to isothermality.” By using shortcuts to isothermality, we derive three nonequilibrium work relations, including an identity between the free-energy difference and the mean work due to the potential of the original system, a Jarzynski-like equality, and the inverse relationship between the dissipated work and the total driving time. We numerically test these three relations by considering the motion of a Brownian particle trapped in a harmonic potential and dragged by a time-dependent force.

We study critical and magnetic properties of a bilayer Ising system consisting of two triangular planes A and B, with the antiferromagnetic (AF) coupling $J_{\mathrm{A}}$ and the ferromagnetic (FM) one $J_{\mathrm{B}}$ for the respective layers, which are coupled by the interlayer interaction $J_{\mathrm{AB}}$ by using Monte Carlo simulations. When $J_{\mathrm{A}}$ and $J_{\mathrm{B}}$ are of the same order, the unfrustrated FM plane orders first at a high temperature $T_{c1}\sim J_{\mathrm{B}}$. The spontaneous FM order then exerts influence on the other frustrated AF plane as an effective magnetic field, which subsequently induces a ferrimagnetic order in this plane at low temperatures below $T_{c2}$. When short-range order is developed in the AF plane while the influence of the FM plane is still small, there appears a preemptive Berezinskii-Kosterlitz-Thouless–type pseudocritical crossover regime just above the ferrimagnetic phase transition point, where the short-distance behavior up to a rather large length scale exponentially diverging in $\propto J_{\mathrm{A}}/T$ is controlled by a line of Gaussian fixed points at $T=0$. In the crossover region, a continuous variation in the effective critical exponent $\frac{4}{9}\lesssim {\eta }^{\mathrm{eff}}\lesssim \frac{1}{2}$ is observed. The phase diagram by changing the ratio $J_{\mathrm{A}}/J_{\mathrm{B}}$ is also investigated.

For a fixed-energy Manna sandpile model belonging to a Manna class in one dimension $\left(d=1\right)$, we recently showed that the critical decay is different for random and regular initial conditions (ICs). Compared with previous results of natural IC for several models, we suggested for the Manna class that the critical decay depends on the characteristics of the three ICs. But the dependence on the random and regular ICs was shown only for a single model. In this work, we study the critical decay for the random and regular ICs for another model of the Manna class in $d=1$, a diffusive epidemic process. It is shown that the critical decay exponent agrees with the previous result for each IC, which verifies that IC dependence is a common feature of the Manna class. In addition, for the random and regular ICs, we measure the variance ${\sigma }^2\left(r\right)$ of total particle density in a region of size $r$ by increasing $r$ up to system size and investigate its temporal evolution toward the value ${\sigma }_q^2\left(r\right)$ of the quasisteady state at criticality. In $d=1,{\sigma }^2\left(r\right)$ scales as ${\sigma }^2\left(r\right)\sim r^{-\psi }$ with $\psi =1$ for random distributions and $1<\psi \le 2$ for hyperuniform ones. The temporal evolution shows that ${\sigma }^2\left(r\right)$ of the two ICs differently relax toward ${\sigma }_q^2\left(r\right)$ and the regular IC becomes a hyperuniform distribution of $\psi =2$ in the beginning of the evolution. We estimate $\psi =1.45\left(3\right)$ for both the quasisteady state and absorbing states, so the quasisteady state is also as hyperuniform as absorbing states. The hyperuniformity of the quasisteady state shows that the natural IC also should be hyperuniform as much as the quasisteady state, because the natural IC is obtained from particle configurations close to the quasisteady state. Consequently, the different $\psi$ of the three ICs suggest that ${\sigma }^2\left(r\right)$ can classify the characteristics of the three ICs in a unified way and the different degree of hyperuniformity of the ICs provides another explanation for the observed IC-dependent critical decay in a point of view of initial fluctuations and correlations.

Using numerical simulations we investigate the space-time properties of a system in which spirals emerge within coarsening domains, thus giving rise to nontrivial internal dynamics. Initially proposed in the context of population dynamics, the studied six-species model exhibits growing domains composed of three species in a rock-paper-scissors relationship. Through the investigation of different quantities, such as space-time correlations and the derived characteristic length, autocorrelation, density of empty sites, and interface width, we demonstrate that the nontrivial dynamics inside the domains affects the coarsening process as well as the properties of the interfaces separating different domains. Domain growth, aging, and interface fluctuations are shown to be governed by exponents whose values differ from those expected in systems with curvature driven coarsening.

Phase oscillator lattices subject to noise are one of the most fundamental systems in nonequilibrium physics. We have discovered a dynamical transition which has a significant impact on the synchronization dynamics in such lattices, as it leads to an explosive increase of the phase diffusion rate by orders of magnitude. Our analysis is based on the widely applicable Kuramoto-Sakaguchi model, with local couplings between oscillators. For one-dimensional lattices, we observe the universal evolution of the phase spread that is suggested by a connection to the theory of surface growth, as described by the Kardar-Parisi-Zhang (KPZ) model. Moreover, we are able to explain the dynamical transition both in one and two dimensions by connecting it to an apparent finite-time singularity in a related KPZ lattice model. Our findings have direct consequences for the frequency stability of coupled oscillator lattices.

We investigate the mechanisms of anomalous diffusion in cationic surfactant micelles using molecular dynamics simulations in the presence of explicit salt and solvent-mediated interactions. Simulations show that when the counterion density increases, saddle-shaped branched interfaces manifest. In experiments, branched structures exhibit lower viscosity as compared to linear and wormlike micelles. This has long been attributed to stress relaxation arising from the sliding motion of branches along the main chain. Our simulations reveal a mechanism of branch motion resulting from an enhanced counterion condensation at the branched interfaces and provide quantitative evidence of stress relaxation facilitated by branched sliding. Furthermore, depending on the surfactant and salt concentrations, which in turn determine the microstructure, we observe normal, subdiffusive, and superdiffusive motions of surfactants. Specifically, superdiffusive behavior is associated with branch sliding, breakage and recombination of micelle fragments, as well as constraint release in entangled systems.

We report a study of how a bend in a quasi-one-dimensional (q1D) channel containing a colloid suspension at equilibrium that exhibits single-file particle motion affects the hydrodynamic coupling between colloid particles. We observe both structural and dynamical responses as the bend angle becomes more acute. The structural response is an increasing depletion of particles in the vicinity of the bend and an increase in the nearest-neighbor separation in the pair correlation function for particles on opposite sides of the bend. The dynamical response monitored by the change in the self-diffusion $[D_{11}\left(x\right)]$ and coupling $[D_{12}\left(x\right)]$ terms of the pair diffusion tensor reveals that the pair separation dependence of $D_{12}$ mimics that of the pair correlation function just as in a straight q1D channel. We show that the observed behavior is a consequence of the boundary conditions imposed on the q1D channel: both the single-file motion and the hydrodynamic flow must follow the channel around the bend.

We combine experiments and theory to investigate the dynamics and orientational fluctuations of ferromagnetic microellipsoids that form a ribbonlike structure due to attractive dipolar forces. When assembled in the ribbon, the ellipsoids display orientational thermal fluctuations with an amplitude that can be controlled via application of an in-plane magnetic field. We use video microscopy to investigate the orientational dynamics in real time and space. Theoretical arguments are used to derive an analytical expression that describes how the distribution of the different angular configurations depends on the strength of the applied field. The experimental data are in good agreement with the developed model for all the range of field parameters explored. Understanding the role of fluctuations in chains composed of dipolar particles is important not only from a fundamental point of view, but it may also help understanding the stability of such structures against thermal noise, which is relevant in microfluidics and laboratory-on-a-chip applications.