Source author record

Hidetsugu Sakaguchi

Hidetsugu Sakaguchi appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

50works
22topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

50 published item(s)

preprint2022arXiv

One- and two-dimensional solitons in spin-orbit-coupled Bose-Einstein condensates with fractional kinetic energy

We address effects of spin-orbit coupling (SOC), phenomenologically added to a two-component Bose-Einstein condensate composed of particles moving by Levy flights, in one- and two-dimensional (1D and 2D) settings. The corresponding system of coupled Gross-Pitaevskii equations includes fractional kinetic-energy operators, characterized by the Levy index, α< 2 (the normal kinetic energy corresponds to α= 2). The SOC terms, with strength λ, produce strong effects in the 2D case: they create families of stable solitons of the semi-vortex (SV) and mixed-mode (MM) types in the interval of 1 < α< 2, where the supercritical collapse does not admit the existence of stable solitons in the absence of the SOC. At λ--> 0, amplitudes of these solitons vanish as (λ)^{1/(α- 1)}.

preprint2021arXiv

Dragging spin-orbit-coupled solitons by a moving optical lattice

It is known that the interplay of the spin-orbit-coupling (SOC) and mean-field self-attraction creates stable two-dimensional (2D) solitons (ground states) in spinor Bose-Einstein condensates. However, SOC destroys the system's Galilean invariance, therefore moving solitons exist only in a narrow interval of velocities, outside of which the solitons suffer delocalization. We demonstrate that the application of a relatively weak moving optical lattice (OL), with the 2D or quasi-1D structure, makes it possible to greatly expand the velocity interval for stable motion of the solitons. The stability domain in the system's parameter space is identified by means of numerical methods. In particular, the quasi-1D OL produces a stronger stabilizing effect than its full 2D counterpart. Some features of the domain are explained analytically.

preprint2020arXiv

Singular mean-field states: A brief review of recent results

This article provides a focused review of recent findings which demonstrate, in some cases quite counter-intuitively, the existence of bound states with a singularity of the density pattern at the center, while the states are physically meaningful because their total norm converges. One model of this type is based on the 2D Gross-Pitaevskii equation (GPE) which combines the attractive potential ~ 1/r^2 and the quartic self-repulsive nonlinearity, induced by the Lee-Huang-Yang effect (quantum fluctuations around the mean-field state). The GPE demonstrates suppression of the 2D quantum collapse, driven by the attractive potential, and emergence of a stable ground state (GS), whose density features an integrable singularity ~1/r^{4/3} at r --> 0. Modes with embedded angular momentum exist too, and they have their stability regions. A counter-intuitive peculiarity of the model is that the GS exists even if the sign of the potential is reversed from attraction to repulsion, provided that its strength is small enough. This peculiarity finds a relevant explanation. The other model outlined in the review includes 1D, 2D, and 3D GPEs, with the septimal (seventh-order), quintic, and cubic self-repulsive terms, respectively. These equations give rise to stable singular solitons, which represent the GS for each dimension D, with the density singularity ~1/r^{2/(4-D). Such states may be considered as a result of screening of a "bare" delta-functional attractive potential by the respective nonlinearity.

preprint2020arXiv

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, namely, represented by septimal, quintic, and usual cubic terms in the framework of the one-, two-, and three-dimensional (1D, 2D, and 3D) nonlinear Schroedinger 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 delta-functional potential by the defocusing nonlinearity. The strength ("bare charge") of the attractive potential is infinite in 1D, finite in 2D, and vanishingly small in 3D. Analytical asymptotics of the singular solitons at small and large distances are found, entire shapes of the solitons being produced in a numerical form. Complete stability of the singular modes is accurately predicted by the anti-Vakhitov-Kolokolov criterion (under the assumption that it applies to the model), as verified by means of numerical methods. In 2D, the NLSE with a quintic self-focusing term admits singular soliton solutions with intrinsic vorticity too, but they are fully unstable. We also mention that dissipative singular solitons can be produced by the model with a complex coefficient in front of the nonlinear term.

preprint2020arXiv

Star shaped patterns caused by colloidal aggregation during the spreading process of a droplet

This research found that when an acidic solution with a low surface tension spread on the surface of a glycerol solution mixed with milk, a star shaped pattern was spontaneously formed on the surface in the horizontal plane during the spreading process. We investigated the emergence of the star shaped pattern owing to an interfacial instability in experiments using glycerol solutions with several viscosities and 2-methoxyethanol aqueous solutions, which are acidic solutions, with several concentrations. This result demonstrated that the star shaped pattern emerged in the high concentration of 2-methoxyethanol. We proposed a phenomenological model, based on our experimental results, which explains three points as follows; the spreading of the 2-methoxyethanol aqueous solution on the surface of the glycerol solution, the colloidal aggregation of the milk protein colloids caused by the denaturation that occurs when mixed with 2-methoxyethanol, and the accumulation of the aggregates toward the dent regions of the moving interface by a sweeping effect. The model reproduces the formation of the star shaped pattern which was similar to the experimental one. Furthermore, the model provided a phase diagram against the concentration of the 2-methoxyethanol solution and the viscosity of the glycerol solution as control parameters in our experiments. The phase diagram was close to that obtained from our experiments. The results suggest that the above three points are important for the formation of the star shaped pattern.

preprint2020arXiv

Symmetry breaking in a two-component system with repulsive interactions and linear coupling

We extend the well-known theoretical treatment of the spontaneous symmetry breaking (SSB) in two-component systems, combining linear coupling and self-attractive nonlinearity, to a system in which the linear coupling competes with repulsive interactions. First, we address one- and two-dimensional (1D and 2D) ground-state (GS) solutions and 2D vortex states with topological charges S=1 and 2, maintained by a confining harmonic-oscillator potential. The system can be implemented in BEC and optics. By means of the Thomas-Fermi approximation and numerical solution of the underlying coupled Gross-Pitaevskii equations, we demonstrate that SSB takes place, in the GSs and vortices alike, when the cross-component repulsion is stronger that the self-repulsion in each component. The SSB transition is categorized as a supercritical bifurcation, which gives rise to states featuring broken symmetry in an inner area, and intact symmetry in a surrounding layer. Unlike stable GSs and vortices with S=1, the states with S=2 are unstable against splitting. We also address SSB for 1D gap solitons in the system including a lattice potential. In this case, SSB takes place under the opposite condition, i.e., the cross-component repulsion must be weaker than the self-repulsion, and SSB is exhibited by antisymmetric solitons.

preprint2020arXiv

The Gross-Pitaevskii-Poisson model for an ultracold plasma: density waves and solitons

We introduce 1D and 2D models of a degenerate bosonic gas composed of ions with positive and negative charges (cations and anions). The system may exist in the mean-field condensate state, enabling the competition of the Coulomb coupling, contact repulsion, and kinetic energy of the particles, provided that their effective mass is reduced by means of a lattice potential. The model combines the Gross-Pitaevskii (GP) equations for the two-component wave function of the cations and anions, coupled to the Poisson equation for the electrostatic potential mediating the Coulomb coupling. The contact interaction in the GP system can be derived, in the Thomas-Fermi approximation, from a system of three GP equations, which includes the wave function of heavy neutral atoms. In the system with fully repulsive contact interactions, we construct stable spatially periodic patterns (density waves, DWs). The transition to DWs is identified by analysis of the modulational instability of a uniformly mixed neutral state. The DW pattern, which represents the system's ground state (GS), is predicted by a variational approximation. In 2D, a stable pattern is produced too, with a quasi-1D shape. The 1D system with contact self-attraction in each component produces bright solitons of three types: neutral ones, with fully mixed components; dipoles, with the components separated by the inter-species contact repulsion; and quadrupoles, with a layer of one component sandwiched between side lobes formed by the other. The transition from the neutral solitons to dipoles is accurately modeled analytically. A chart of the GSs of the different types (neutral solitons, dipoles, or quadrupoles) is produced. Different soliton species do not coexist as stable states. Collisions between traveling solitons are elastic for dipole-dipole pairs, while dipole-antidipole ones merge into stable quadrupoles via multiple collisions.

preprint2019arXiv

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. A solitary wave appears also in our model similarly to previous models. In some parameter range, the moving direction of the solitary wave sometimes changes rather suddenly, which is like the sudden change of moving direction of a flock of birds. We study the average reversal time of traveling direction numerically and compare the results with a mean-field theory. The one-dimensional model is generalized to a two-dimensional model. Flip motion of a bandlike soliton is observed in the two-dimensional model.

preprint2016arXiv

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 quasi-periodic motions by direct numerical simulations of this coupled oscillator system. A low-dimensional equation for order parameters is derived using the Ott-Antonsen ansatz. The applicability of the ansatz is checked by the comparison of numerical results of the coupled oscillator system and the reduced low-dimensional equation. We investigate further several interesting phenomena in which mutual interactions between the fast and slow oscillators play an essential role. Fast oscillations appear intermittently as a result of excitatory interactions with slow oscillators in a certain parameter range. Slow oscillators experience an oscillator-death phenomenon owing to their interaction with fast oscillators. This oscillator death is explained as a result of saddle-node bifurcation in a simple phase equation obtained using the temporal average of the fast oscillations. Finally we show macroscopic synchronization of the order 1:m between the slow and fast oscillators.

preprint2016arXiv

Matter-wave soliton interferometer based on a nonlinear splitter

We elaborate a model of the interferometer which, unlike previously studied ones, uses a local (delta-functional) nonlinear repulsive potential, embedded into a harmonic-oscillator trapping potential, as the splitter for the incident soliton. An estimate demonstrates that this setting may be implemented by means of the localized Feshbach resonance controlled by a focused laser beam. The same system may be realized as a nonlinear waveguide in optics. Subsequent analysis produces an exact solution for scattering of a plane wave in the linear medium on the delta-functional nonlinear repulsive potential, and an approximate solution for splitting of the incident soliton when the ambient medium is nonlinear. The most essential result, obtained by means of systematic simulations, is that the use of the nonlinear splitter provides the sensitivity of the soliton-based interferometer to the target, inserted into one of its arms, which is much higher than the sensitivity provided by the usual linear splitter.

preprint2016arXiv

One- and two-dimensional solitons in PT-symmetric systems emulating the spin-orbit coupling

We introduce a two-dimensional (2D) system, which can be implemented in dual-core planar optical couplers with the Kerr nonlinearity in its cores, making it possible to blend effects of the PT symmetry, represented by the balanced linear gain and loss in the two cores, and spin-orbit coupling (SOC), emulated by a spatially biased coupling between the cores. Families of 1D and 2D solitons and their stability boundaries are identified. In the 1D setting, the addition of the SOC terms leads, at first, to shrinkage of the stability area for PT-symmetric solitons, which is followed by its rapid expansion. 2D solitons have their stability region too, in spite of the simultaneous action of two major destabilizing factors, viz., the collapse driven by the Kerr nonlinearity, and a trend towards spontaneous breakup of the gain-loss balance. In the limit of the SOC terms dominating over the intrinsic diffraction, the 1D system gives rise to a new model for gap solitons, which admits exact analytical solutions.

preprint2016arXiv

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

We present an analysis of two-dimensional (2D) matter-wave solitons, governed by the pseudo-spinor 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 semi-vortex (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. Existence domains and stability boundaries are identified for the soliton families. For physically relevant parameters of the SOC system, the number of atoms in the 2D solitons is limited by $\sim 1.5\times 10^{4}$. The results are obtained by means of combined analytical and numerical methods.

preprint2016arXiv

Widening Disparity and its Suppression in a Stochastic Replicator Model

Winner-take-all phenomena are observed in various competitive systems. We find similar phenomena in replicator models with randomly fluctuating growth rates. The disparity between winners and losers increases indefinitely, even if all elements are statistically equivalent. A lognormal distribution describes well the nonstationary time evolution. If a nonlinear load corresponding to progressive taxation is introduced, a stationary distribution is obtained and disparity widening is suppressed.

preprint2015arXiv

Aftershocks and Omori's law in a modified Carlson-Langer model with nonlinear visco-elasticity

A modified Carlson-Langer model for earthquakes is proposed, which includes nonlinear visco-elasticity. Several aftershocks are generated after the main shock owing to the damping of the additional visco-elastic 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.

preprint2015arXiv

Anomalous Fluctuations in Autoregressive Models with Long-Term Memory

An autoregressive model with a power-law type memory kernel is studied as a stochastic process that exhibits a self-affine-fractal-like behavior for a small time scale. We find numerically that the root-mean-square displacement for the time interval increases with a power law for small time but saturates at sufficiently large time. The exponent changes with the power exponent of the memory kernel.

preprint2015arXiv

Numerical Study of Crystal Size Distribution in Polynuclear Growth

The crystal size distribution in polynuclear growth is numerically studied using a coupled map lattice model. The width of the size distribution depends on c/D, where c is the growth rate at interface sites and $D$ is the diffusion constant. When c/D is sufficiently small, the width W increases linearly with c/D and saturates at large c/D. Monodisperse square and cubic crystals are obtained respectively on square and cubic lattices when c/D is sufficiently small for a small kinetic parameter b. The linear dependence of W on c/D in a parameter range of small c/D is explained by the eigenfunction for the first eigenvalue in a two-dimensional model and a mean-field model. For the mean-field model, the slope of the linear dependence is evaluated theoretically.

preprint2014arXiv

Domain Wall Dynamics in Ginzburg-Landau-Type Equations with Conservative Quantities

In the Ginzburg-Landau equation, there are domain walls connecting two metastable states. The dynamics of domain walls has been intensively studied, but there remain still unsolved but crucial problems even for a single domain. We study the domain wall dynamics in three different Ginzburg-Landau-type equations satisfying conservation laws. In a modified $ϕ^4$ model satisfying the law of energy conservation and the Lorentz invariance, the motion of a domain wall is accelerated and the velocity approaches its maximum. In a one-dimensional model of eutectic growth, the order parameter is conserved and a domain wall connecting a metastable uniform state and a spatially periodic pattern appears. We try to find a selection rule for the wavelength of a spatially periodic pattern. In a model equation for martensitic transformation, a domain wall connecting a uniform metastable state and a zigzag structure appears which propagates at a high velocity.

preprint2014arXiv

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, which is consistent with observations in mammalian lungs. We find furthermore that there is an optimal branching number when the Peclet number is large. Under the doubly optimal conditions with respect to the size ratio and branching number, the optimal value of r is close to (1/2)^{1/3}.

preprint2014arXiv

The 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 the cubic self-attraction cannot produce stable solitons in free space because of the occurrence of the collapse in the same setting. By means of the 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: semi-vortices (SVs, with a vortex in one component and a fundamental soliton in the other), and mixed modes (MMs, with topological charges 0 and +1/-1 mixed in both components). These two- dimensional composite modes can be created using the trapping harmonic-oscillator (HO) potential, but remain stable in the free space, if the trap is gradually removed. The SVs and MMs realize the ground state of the system, provided that the self-attraction in the two components is, respectively, stronger or weaker than the cross-attraction between them. The SVs and MMs which are not the ground state are subject to a drift instability. In the free space (in the absence of the HO trap), modes of both types degenerate into unstable Townes solitons when their norms attain the respective critical values, while there is no lower existence threshold for the stable modes. Moving free-space stable solitons are also found in the present non-Galilean-invariant system, up to a critical velocity. Collisions between two moving solitons lead to their merger into a single one.

preprint2013arXiv

Dynamical quorum sensing and clustering dynamics in a population of spatially distributed active rotators

A model of clustering dynamics is proposed for a population of spatially distributed active rotators. A transition from excitable to oscillatory dynamics is induced by the increase of the local density of active rotators. It is interpreted as dynamical quorum sensing. In the oscillation regime, phase waves propagate without decay, which generates an effectively long-range interaction in the clustering dynamics. The clustering process becomes facilitated and only one dominant cluster appears rapidly as a result of the dynamical quorum sensing. An exact localized solution is found to a simplified model equation, and the competitive dynamics between two localized states is studied numerically.

preprint2013arXiv

Intermittency and Slip-Size Distribution in the Block-Spring Model of Earthquakes

The Carlson-Langer model is a deterministic model of earthquakes. There were many investigations of this model, but its complicated spatio-temporal dynamics is not yet completely understood. We again study the model equation numerically, and obtain two new results. Firstly, we try to understand the complicated dynamics from the viewpoint of spatio-temporal intermittency. We study the intermittent dynamics of a two-block system using a one-dimensional map, which is useful for understanding an intermittent spatio-temporal chaos in larger systems. Large slips occur intermittently after repeated small slips when the parameter of the velocity weakening friction is large. Small slips reduce the spatial irregularity of block displacement, and large slips occur on a relatively flat configuration of blocks. Secondly, we study the effect of spatial heterogeneity such as asperity using percolation clusters. We show the numerical results of the slip-size distribution in the block-spring model on percolation clusters, and find that the b parameter for the magnitude distribution is close to 2/3 when the velocity weakening parameter is intermediate.

preprint2013arXiv

Phase Model with Feedback Control for Power Grids

A phase model with feedback control is studied as a dynamical model of power grids. As an example, we study a model network corresponding to the power grid in the Kyushu region. The standard frequency is maintained by the mutual synchronization and the feedback control. Electric failures are induced by an overload. We propose a local feedback method in which the strength of feedback control is proportional to the magnitude of generators. We find that the electric failures do not occur until the utilization ratio is close to 1 under this feedback control. We also find that the temporal response for the time-varying input power is suppressed under this feedback control. We explain the mechanisms using the corresponding global feedback method.

preprint2013arXiv

Shock Waves in Falling Coupled Harmonic Oscillators

Shock waves propagate in falling coupled harmonic oscillators. The bottom end of coupled harmonic oscillators does not fall downwards until a shock wave reaches the bottom end. The exact solution can be expressed by the Fourier series expansion, and an approximate solution can be expressed by the integral of the Airy function. The width of the shock wave increases slowly in accordance with a power law.

preprint2013arXiv

Simple Model of Splitting Instability in Swollen Membranes

A simple one-dimensional mechanical model is proposed for splitting instability in swollen membranes. The splitting instability occurs by ring constriction. The bifurcation can be both subcritical and supercritical, depending on the slenderness of spheroidal membranes. The lower critical point is estimated theoretically for the simplest circular system. Necking instability occurs in a spheroidal membrane when the twist angle in the boundary conditions is increased.

preprint2013arXiv

Suppression of the quantum collapse in binary bosonic gases

Attraction of the quantum particle to the center in the 3D space with potential V/r^2 gives rise to the quantum collapse, i.e., nonexistence of the ground state (GS) when the attraction strength exceeds a critical value (V = 1/8, in the present notation). Recently, we have demonstrated that the quantum collapse is suppressed, and the GS is restored, if repulsive interactions between particles in the quantum gas are taken into account, in the mean-field approximation. This setting can be realized in a gas of dipolar molecules attracted to the central charge, with dipole-dipole interactions taken into regard too. Here we analyze this problem for a binary gas. GSs supported by the repulsive interactions are constructed in a numerical form, as well as by means of analytical approximations for both miscible and immiscible binary systems. In particular, the Thomas-Fermi (TF) approximation is relevant if V is large enough. It is found that the GS of the miscible binary gas, both balanced and imbalanced, features a weak phase transition at another critical value, V = 1/2. The transition is characterized by an analyticity-breaking change in the structure of the wave functions at small r. To illustrate the generic character of the present phenomenology, we also consider the binary system with the attraction between the species (rather than repulsion), in the case when the central potential pulls a single component only.

preprint2013arXiv

Vortex lattice solutions to the Gross-Pitaevskii equation with spin-orbit coupling in optical lattices

Effective spin-orbit coupling can be created in cold atom systems using atom-light interaction. We study the BECs in an optical lattice using the Gross-Pitaevskii equation with spin-orbit coupling. Bloch states for the linear equation are numerically obtained, and compared with stationary solutions to the Gross-Pitaevskii equation with nonlinear terms. Various vortex lattice states are found when the spin-orbit coupling is strong.

preprint2013arXiv

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). They exist above a threshold value of the total energy (norm) and below another critical value, which determines a stability boundary. The existence threshold vanishes is a part of the parameter space, depending on the mismatch parameter, which is explained by means of the comparison with the 2D counterpart of the system. Above the stability boundary, the vortex features shape oscillations, periodically breaking its axisymmetric form and restoring it. Collisions between vortices moving in the longitudinal direction are studied too. The collision is strongly inelastic at relatively small values of the velocities, breaking the two colliding vortices into three, with the same vorticity. The results suggest a possibility of the creation of stable 3D optical solitons with the intrinsic vorticity.

preprint2012arXiv

Cascade Failure in a Phase Model of Power Grids

We propose a phase model to study cascade failure in power grids composed of generators and loads. If the power demand is below a critical value, the model system of power grids maintains the standard frequency by feedback control. On the other hand, if the power demand exceeds the critical value, an electric failure occurs via step out (loss of synchronization) or voltage collapse. The two failures are incorporated as two removal rules of generator nodes and load nodes. We perform direct numerical simulation of the phase model on a scale-free network and compare the results with a mean-field approximation.

preprint2012arXiv

Melting Transition of Vortex Lattice in Point Vortex Systems

Point vortices take a triangular lattice structure in a rotating system as a minimum energy state. We perform a numerical simulation of point vortex systems using initial conditions indicating that the triangular lattice is randomly perturbed. The total energy increases with the magnitude of the perturbation. When the energy is increased, the vortex lattice becomes irregular and a layered structure appears. When the energy is further increased, the layered structure disappears and a liquidlike state appears. We interpret the melting transition with a mean-field approximation for layered structures.

preprint2012arXiv

Reflection of Channel-Guided Solitons at Junctions in Two-Dimensional Nonlinear Schroedinger Equation

Solitons confined in channels are studied in the two-dimensional nonlinear Schrödinger equation. We study the dynamics of two channel-guided solitons near the junction where two channels are merged. The two solitons merge into one soliton, when there is no phase shift. If a phase difference is given to the two solitons, the Josephson oscillation is induced. The Josephson oscillation is amplified near the junction. The two solitons are reflected when the initial velocity is below a critical value.

preprint2012arXiv

Stabilizing single- and two-color vortex beams in quadratic media by a trapping potential

We consider two-dimensional (2D) localized modes in the second-harmonic-generating χ^{(2)} system with the harmonic-oscillator (HO) trapping potential. In addition to its realization in optics, the system describes the mean-field dynamics of mixed atomic-molecular Bose-Einstein condensates (BECs). The existence and stability of various modes is determined by their total power, N, topological charge, m/2 [m is the intrinsic vorticity of the second-harmonic (SH) field], and χ^{(2)} mismatch, q. The analysis is carried out in a numerical form and, in parallel, by means of the variational approximation (VA), which produces results that agree well with numerical findings. Below a certain power threshold, N\leq N_{c}^{(m)}(q), all trapped modes are of the single-color type, represented by the SH component only, while the fundamental-frequency (FF) one is absent. In contrast with the usual situation, where such modes are always unstable, we demonstrate\ that they are stable, for m=0,1,2 (the mode with m=1 may be formally considered as a semi-vortex with topological charge m/2=1/2), at N\leq N_{c}^{(m)}(q), and unstable above this threshold. On the other hand, N_{c}^{(m)}(q)\equiv 0 at q\geq q_{\max} (in our notation, q_{\max}=1$), hence the single-color modes are unstable in the latter case. At N=N_{c}^{(m)}, the modes with m=0 and m=2 undergo a pitchfork bifurcation, which gives rise to two-color states, which remain completely stable for m=0. The two-color vortices with m=2 (topological charge 1) have an upper stability border, N=N_{c2}(q). Above the border, they exhibit periodic splittings and recombinations, while keeping their vorticity. The semi-vortex does not bifurcate; at N=N_{c}^{(m=1)}, it exhibits quasi-chaotic oscillations and a rotating "groove" resembling a screw-edge dislocation induced by the semi-integer vorticity.

preprint2012arXiv

Stable two-dimensional solitons supported by radially inhomogeneous self-focusing nonlinearity

We demonstrate that modulation of the local strength of the cubic self-focusing (SF) nonlinearity in the two-dimensional (2D) geometry, in the form of a circle with contrast $Δg$ of the SF coefficient relative to the ambient medium with a weaker nonlinearity, stabilizes a family of fundamental solitons against the critical collapse. The result is obtained in an analytical form, using the variational approximation (VA) and Vakhitov-Kolokolov (VK) stability criterion, and corroborated by numerical computations. For the small contrast, the stability interval of the soliton's norm scales as $ΔN\sim Δg$ (the replacement of the circle by an annulus leads to a reduction of the stability region by perturbations breaking the axial symmetry). To further illustrate this mechanism, we demonstrate, in an exact form, the stabilization of 1D solitons against the critical collapse under the action of a locally enhanced quintic SF nonlinearity.

preprint2011arXiv

Locomotive and reptation motion induced by internal force and friction

We propose a simple mechanical model of locomotion induced by internal force and friction. We first construct a system of two elements as an analog of the bipedal motion. The internal force does not induce a directional motion by itself because of the action-reaction law, but a directional motion becomes possible by the control of the frictional force. The efficiency of these model systems is studied using an analogy to the heat engine. As a modified version of the two-elements model, we construct a model which exhibits a bipedal motion similar to kinesin's motion of molecular motor. Next, we propose a linear chain model and a ladder model as an extension of the original two-element model,. We find a transition from a straight to a snake-like motion in a ladder model by changing the strength of the internal force.

preprint2011arXiv

Mechanism of Stepped Leaders in a Simple Discharge Model

We construct a one-dimensional model for the stepped leader in the filamental discharge by simplifying an electric-circuit model of discharge. We find that the leader of the discharge moves stepwise by direct numerical simulations, and then we try to understand the mechanism of the stepwise motion by reducing the spatially extended system to the dynamics of the tip position of the discharge.

preprint2010arXiv

Branching Patterns and Stepped Leaders in an Electric-Circuit Model for Creeping Discharge

We construct a two-dimensional electric circuit model for creeping discharge. Two types of discharge, surface corona and surface leader, are modeled by a two-step function of conductance. Branched patterns of surface leaders surrounded by the surface corona appear in numerical simulation. The fractal dimension of branched discharge patterns is calculated by changing voltage and capacitance. We find that surface leaders often grow stepwise in time, as is observed in lightning leaders of thunder.

preprint2010arXiv

Reflection and Splitting of Channel-Guided Solitons in Two-dimensional Nonlinear Schrödinger Equation

Solitons confined in a channel are studied in the two-dimensional nonlinear Schroedinger equation. When a channel branches into two channels, a soliton is split into two solitons, if the initial kinetic energy exceeds a critical value. The branching point works as a pulse splitter. If it is below the critical value, the soliton is reflected. The critical kinetic energy for splitting is evaluated by a variational method. The variational method can be applied in the design of channel systems with reduced reflection.

preprint2010arXiv

Solitons in combined linear and nonlinear lattice potentials

We study ordinary solitons and gap solitons (GSs) in the effectively one-dimensional Gross-Pitaevskii equation, with a combination of linear and nonlinear lattice potentials. The main points of the analysis are effects of the (in)commensurability between the lattices, the development of analytical methods, viz., the variational approximation (VA) for narrow ordinary solitons, and various forms of the averaging method for broad solitons of both types, and also the study of mobility of the solitons. Under the direct commensurability (equal periods of the lattices, the family of ordinary solitons is similar to its counterpart in the free space. The situation is different in the case of the subharmonic commensurability, with L_{lin}=(1/2)L_{nonlin}, or incommensurability. In those cases, there is an existence threshold for the solitons, and the scaling relation between their amplitude and width is different from that in the free space. GS families demonstrate a bistability, unless the direct commensurability takes place. Specific scaling relations are found for them too. Ordinary solitons can be readily set in motion by kicking. GSs are mobile too, featuring inelastic collisions. The analytical approximations are shown to be quite accurate, predicting correct scaling relations for the soliton families in different cases. The stability of the ordinary solitons is fully determined by the VK (Vakhitov-Kolokolov) criterion, while the stability of GS families follows an inverted ("anti-VK") criterion, which is explained by means of the averaging approximation.

preprint2010arXiv

Suppression of the quantum-mechanical collapse by repulsive interactions in a quantum gas

The quantum-mechanical collapse (alias fall onto the center of particles attracted by potential -1/r^2), or "quantum anomaly", is a well-known issue in the quantum theory. We demonstrate that the mean-field repulsive nonlinearity prevents the collapse and thus puts forward a solution to the quantum-anomaly problem different from that previously developed in the framework of the linear quantum-field theory. This solution may be realized in the 3D or 2D gas of dipolar bosons attracted by a central charge, and in the 2D gas of magnetic dipoles attracted by a current filament. In the 3D setting, the dipole-dipole interactions are also taken into regard, in the mean-field approximation. In lieu of the collapse, the cubic nonlinearity creates a 3D ground state (GS), which does not exist in the respective linear Schroedinger equation (SE). The addition of the harmonic trap gives rise to a tristability, in the case when the SE still does not lead to the collapse. In the 2D setting, the cubic nonlinearity is not strong enough to prevent the collapse; however, the quintic term does it, creating the GS, as well as its counterparts carrying the angular momentum (vorticity). Counter-intuitively, such self-trapped 2D modes exist even in the case of a weakly repulsive potential 1/r^2. In the presence of the harmonic trap, the 2D quintic model with a weakly repulsive central potential 1/r^2 gives rise to three confined modes, the middle one being unstable, spontaneously developing into a breather. In both the 3D and 2D cases, the GS wave functions are found in a numerical form, and also in the form of an analytical approximation, which is asymptotically exact in the limit of the large norm.

preprint2009arXiv

Cooperative jump motions of jammed particles in a one-dimensional periodic potential

Cooperative jump motions are studied for mutually interacting particles in a one-dimensional periodic potential. The diffusion constant for the cooperative motion in systems including a small number of particles is numerically calculated and it is compared with theoretical estimates. We find that the size distribution of the cooperative jump motions obeys an exponential law in a large system.

preprint2009arXiv

Kinetic Friction by a Small Number of Intervening Inelastic Particles between Rough Surfaces

We investigate a mechanism of the appearance of kinetic friction in granular materials. We consider a small number of intervening inelastic particles between two rough surfaces as one of the simplest dynamical models to study granular friction. The resistance force applied to the upper surface is numerically calculated. We find that the resistance force F(t) is scaled as F'(vt) for a small pulling velocity v. The time average F_0=<F(t)> in the limit v->0 is not zero owing to the mutual collisions between the intervening particles. The nonzero F_0 implies the appearance of kinetic friction in this simple dynamical system.

preprint2009arXiv

Ratio control in a cascade model of cell differentiation

We propose a kind of reaction-diffusion equations for cell differentiation, which exhibits the Turing instability. If the diffusivity of some variables is set to be infinity, we get coupled competitive reaction-diffusion equations with a global feedback term. The size ratio of each cell type is controlled by a system parameter in the model. Finally, we extend the model to a cascade model of cell differentiation. A hierarchical spatial structure appears as a result of the cell differentiation. The size ratio of each cell type is also controlled by the system parameter.

preprint2009arXiv

Two-dimensional dissipative gap solitons

We introduce a model which integrates the complex Ginzburg-Landau (CGL) equation in two dimensions (2D) with the linear-cubic-quintic combination of loss and gain terms, self-defocusing nonlinearity, and a periodic potential. In this system, stable 2D \textit{dissipative gap solitons} (DGSs) are constructed, both fundamental and vortical ones. The soliton families belong to the first finite bandgap of the system's linear spectrum. The solutions are obtained in a numerical form and also by means of an analytical approximation, which combines the variational description of the shape of the fundamental and vortical solitons and the balance equation for their total power. The analytical results agree with numerical findings. The model may be implemented as a laser medium in a bulk self-defocusing optical waveguide equipped with a transverse 2D grating, the predicted DGSs representing spatial solitons in this setting.

preprint2008arXiv

Synchronization in coupled phase oscillators

We make a short review about the synchronization in coupled phase oscillator models. Next, we study the common-noise-induced synchronization among active rotators. At an intermediate noise strength, the noise-induced synchronization takes place most effectively, which is analogous to the stochastic resonance. Finally, we study the synchronization of coupled phase oscillators with nonvariational interaction on scale-free networks. We find a sharp transition and a weak hysteresis in the nonvariational systems. The sharp transition is found also in the mean-field approximation.

preprint2005arXiv

Higher-order vortex solitons, multipoles, and supervortices on a square optical lattice

We predict new generic types of vorticity-carrying soliton complexes in a class of physical systems including an attractive Bose-Einstein condensate in a square optical lattice (OL) and photonic lattices in photorefractive media. The patterns include ring-shaped higher-order vortex solitons and supervortices. Stability diagrams for these patterns, based on direct simulations, are presented. The vortex ring solitons are stable if the phase difference Δϕbetween adjacent solitons in the ring is larger than π/2, while the supervortices are stable in the opposite case, Δϕ<π/2. A qualitative explanation to the stability is given.