Nonequilibrium dynamics of coupled oscillators under the shear-velocity boundary condition.

Deterministic and stochastic coupled oscillators with inertia are studied on the rectangular lattice under the shear-velocity boundary condition. Our coupled oscillator model exhibits various nontrivial phenomena and there are various relationships with wide research areas such as the coupled limit-cycle oscillators, the dislocation theory, a block-spring model of earthquakes, and the nonequilibrium molecular dynamics. We show numerically several unique nonequilibrium properties of the coupled oscillators.

Noise-induced excitation wave and its size distribution in coupled FitzHugh-Nagumo equations on a square lattice.

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.

Suppression and frequency control of repetitive spiking in the FitzHugh-Nagumo model.

The FitzHugh-Nagumo equation is a simple model equation that exhibits spiking. The output signal of a neuron is represented in the spiking frequency or firing rate. We consider a few control methods for spiking from the viewpoint of nonlinear dynamics.

Vortex motion and nonlinear response in coupled noisy phase oscillator lattices under shear stress.

Vortex motion in coupled phase oscillator lattices is analogous to the dislocation motion in crystals. A single vortex exhibits a glide motion by force at the boundaries. Thermal fluctuations induce the glide motion even below the critical point corresponding to the Peierls stress.

Vortex motions in coupled phase oscillator lattices with inertia under shear stress.

We propose a coupled phase oscillator model with inertia and study the vortex motion in the model when the external force is applied at the boundaries. The vortex exhibits a glide motion when the external force is larger than a critical value. We find a transition from the pair annihilation to passing for the collision of the vortex and antivortex when the external force is changed.

Slow decay of infection in the inhomogeneous susceptible-infected-recovered model.

The susceptible-infected-recovered (SIR) model with spatially inhomogeneous infection rate is studied with numerical simulations in one, two, and three dimensions, considering the case that the infection spreads inhomogeneously in densely populated regions or hot spots. We find that the total population of infection decays very slowly in the inhomogeneous systems in some cases, in contrast to the exponential decay of the infected population I(t) in the SIR model of the ordinary differential equation. The slow decay of the infected population suggests that the infection is locally maintained for long and it is difficult for the disease to disappear completely.

Singular solitons.

We demonstrate that the commonly known concept which treats solitons as nonsingular solutions produced by the interplay of nonlinear self-attraction and linear dispersion may be extended to include modes with a relatively weak singularity at the central point, which keeps their integral norm convergent. Such states are generated by self-repulsion, which should be strong enough, represented by septimal, quintic, and usual cubic terms in the framework of the one-, two-, and three-dimensional (1D, 2D, and 3D) nonlinear Schrödinger equations (NLSEs), respectively. Although such solutions seem counterintuitive, we demonstrate that they admit a straightforward interpretation as a result of screening of an additionally introduced attractive δ-functional potential by the defocusing nonlinearity.

Flip motion of solitary wave in an Ising-type Vicsek model.

An Ising-type Vicsek model is proposed for collective motion and sudden direction change in a population of self-propelled particles. Particles move on a linear lattice with velocity +1 or -1 in the one-dimensional model. The probability of the velocity of a particle at the next step is determined by the number difference of the right- and left-moving particles at the present lattice site and its nearest-neighboring sites.

Interactions of solitons with positive and negative masses: Shuttle motion and coacceleration.

We consider a possibility to realize self-accelerating motion of interacting states with effective positive and negative masses in the form of pairs of solitons in two-component BEC loaded in an optical-lattice (OL) potential. A crucial role is played by the fact that gap solitons may feature a negative dynamical mass, keeping their mobility in the OL. First, the respective system of coupled Gross-Pitaevskii equations (GPE) is reduced to a system of equations for envelopes of the lattice wave functions.

Stationary and oscillatory bound states of dissipative solitons created by third-order dispersion.

We consider the model of fiber-laser cavities near the zero-dispersion point, based on the complex Ginzburg-Landau equation with the cubic-quintic nonlinearity and third-order dispersion (TOD) term. It is known that this model supports stable dissipative solitons. We demonstrate that the same model gives rise to several specific families of robust bound states of solitons.

Transmission coefficient from generalized Cantor-like potentials and its multifractality.

We study the scattering problem at generalized Cantor-like potentials characterized by the expansion rate a and duplication number N, and derive an exact formula of transmittance. It was found that the transmittance is expressed with Chebyshev polynomials of the second kind, and the multifractality of the reflectance varies depending on a and N.

Rayleigh-Taylor instability and vortex rings in coupled Gross-Pitaevskii equations.

The Rayleigh-Taylor instability is a gravitational instability in two fluids where the heavier fluid is set over the lighter fluid. The instability occurs both in classical fluids and quantum fluids. We numerically study the Rayleigh-Taylor instability using coupled Gross-Pitaevskii equations for two-component Bose-Einstein condensates.

Scaling laws of reflection coefficients of quantum waves at a Cantor-like potential.

We reconsider a one-dimensional scattering problem in the Schrödinger equation with a Cantor-like potential. The reflection coefficient obeys a scaling law for sufficiently large wave number k. The scaling law is expressed with a universal function characterized by a multifractal.

Vortex solitons in two-dimensional spin-orbit coupled Bose-Einstein condensates: Effects of the Rashba-Dresselhaus coupling and Zeeman splitting.

We present an analysis of two-dimensional (2D) matter-wave solitons, governed by the pseudospinor system of Gross-Pitaevskii equations with self- and cross attraction, which includes the spin-orbit coupling (SOC) in the general Rashba-Dresselhaus form, and, separately, the Rashba coupling and the Zeeman splitting. Families of semivortex (SV) and mixed-mode (MM) solitons are constructed, which exist and are stable in free space, as the SOC terms prevent the onset of the critical collapse and create the otherwise missing ground states in the form of the solitons. The Dresselhaus SOC produces a destructive effect on the vortex solitons, while the Zeeman term tends to convert the MM states into the SV ones, which eventually suffer delocalization.

Breakup and deformation of a droplet falling in a miscible solution.

When a droplet with a higher density falls in a miscible solution, the droplet deforms and breaks up. The instability of a vortex ring, formed by droplet deformation during the falling process, causes the breakup. To determine the origin of the instability, the wavelengths and thicknesses of the vortex rings are investigated at the time when the instability occurs.

Cooperative dynamics in coupled systems of fast and slow phase oscillators.

We propose a coupled system of fast and slow phase oscillators. We observe two-step transitions to quasiperiodic motions by direct numerical simulations of this coupled oscillator system. A low-dimensional equation for order parameters is derived using the Ott-Antonsen ansatz.

Pattern formation in a sandpile of ternary granular mixtures.

Pattern formation in a sandpile is investigated by pouring a ternary mixture of grains into a vertical narrow cell. Size segregation in avalanches causes the formation of patterns. Four kinds of patterns emerge: stratification, segregation, upper stratification-lower segregation, and upper segregation-lower stratification.

Aftershocks and Omori's law in a modified Carlson-Langer model with nonlinear viscoelasticity.

A modified Carlson-Langer model for earthquakes is proposed, which includes nonlinear viscoelasticity. Several aftershocks are generated after the main shock owing to the damping of the additional viscoelastic force. Both the Gutenberg-Richter law and Omori's law are reproduced in a numerical simulation of the modified Carlson-Langer model on a critical percolation cluster of a square lattice.

Discrete and continuum composite solitons in Bose-Einstein condensates with the Rashba spin-orbit coupling in one and two dimensions.

We introduce one- and two-dimensional (1D and 2D) continuum and discrete models for the two-component BEC, with the spin-orbit (SO) coupling of the Rashba type between the components, and attractive cubic interactions, assuming that the condensate is fragmented into a quasidiscrete state by a deep optical-lattice potential. In 1D, it is demonstrated, in analytical and numerical forms, that the ground states of both the discrete system and its continuum counterpart switch from striped bright solitons, featuring deep short-wave modulations of its profile, to smooth solitons, as the strength ratio of the inter- and intracomponent attraction, γ, changes from γ<1 to γ>1. At the borderline, γ=1, there is a continuous branch of stable solitons, which share a common value of the energy and interpolate between the striped and smooth ones.

Reaction-diffusion-advection equation in binary tree networks and optimal size ratio.

A simple reaction-diffusion-advection equation is proposed in a dichotomous tree network to discuss an optimal network. An optimal size ratio r is evaluated by the principle of maximization of total reaction rate. In the case of reaction-limited conditions, the optimal ratio can be larger than (1/2)(1/3) for a fixed value of branching number n, which is consistent with observations in mammalian lungs.

Competitive aggregation dynamics using phase wave signals.

Coupled equations of the phase equation and the equation of cell concentration n are proposed for competitive aggregation dynamics of slime mold in two dimensions. Phase waves are used as tactic signals of aggregation in this model. Several aggregation clusters are formed initially, and target patterns appear around the localized aggregation clusters.

Creation of two-dimensional composite solitons in spin-orbit-coupled self-attractive Bose-Einstein condensates in free space.

It is commonly known that two-dimensional mean-field models of optical and matter waves with cubic self-attraction cannot produce stable solitons in free space because of the occurrence of collapse in the same setting. By means of numerical analysis and variational approximation, we demonstrate that the two-component model of the Bose-Einstein condensate with the spin-orbit Rashba coupling and cubic attractive interactions gives rise to solitary-vortex complexes of two types: semivortices (SVs, with a vortex in one component and a fundamental soliton in the other), and mixed modes (MMs, with topological charges 0 and ±1 mixed in both components). These two-dimensional composite modes can be created using the trapping harmonic-oscillator (HO) potential, but remain stable in free space, if the trap is gradually removed.

Modulational instability and breathing motion in the two-dimensional nonlinear Schrödinger equation with a one-dimensional harmonic potential.

Modulational instability and breathing motion are studied in the two-dimensional nonlinear Schrödinger (NLS) equation trapped by the one-dimensional harmonic potential. The trapping potential is uniform in the y direction and the wave function is confined in the x direction. A breathing motion appears when the initial condition is close to a stationary solution which is uniform in the y direction.

Vortical light bullets in second-harmonic-generating media supported by a trapping potential.

We introduce a three-dimensional (3D) model of optical media with the quadratic (χ((2))) nonlinearity and an effective 2D isotropic harmonic-oscillator (HO) potential. While it is well known that 3D χ((2)) solitons with embedded vorticity ("vortical light bullets") are unstable in the free space, we demonstrate that they have a broad stability region in the present model, being supported by the HO potential against the splitting instability. The shape of the vortical solitons may be accurately predicted by the variational approximation (VA).