Daily Papers

We introduce spatiotemporal solitons of the two-dimensional complex Ginzburg-Landau equation (2D CCQGLE) with cubic and quintic nonlinearities in which asymmetry between space-time variables is included. The 2D CCQGLE is solved by a powerful Fourier spectral method, i.e., a Fourier spatial discretization and an explicit scheme for time differencing. Varying the system's parameters, and using different initial conditions, numerical simulations reveal 2D solitons in the form of stationary, pulsating and exploding solitons which possess very distinctive properties. For certain regions of parameters, we have also found stable coherent structures in the form of spinning (vortex) solitons which exist as a result of a competition between focusing nonlinearities and spreading while propagating through medium.

By using ZEUS cluster at Embry-Riddle Aeronautical University we perform extensive numerical simulations based on a two-dimensional Fourier spectral method Fourier spatial discretization and an explicit scheme for time differencing) to find the range of existence of the spatiotemporal solitons of the two-dimensional complex Ginzburg-Landau equation with cubic and quintic nonlinearities. We start from the parameters used by Akhmediev {\it et. al.} and slowly vary them one by one to determine the regimes where solitons exist as stable/unstable structures. We present eight classes of dissipative solitons from which six are known (stationary, pulsating, vortex spinning, filament, exploding, creeping) and two are novel (creeping-vortex propellers and spinning "bean-shaped" solitons). By running lengthy simulations for the different parameters of the equation, we find ranges of existence of stable structures (stationary, pulsating, circular vortex spinning, organized exploding), and unstable structures (elliptic vortex spinning that leads to filament, disorganized exploding, creeping). Moreover, by varying even the two initial conditions together with vorticity, we find a richer behavior in the form of creeping-vortex propellers, and spinning "bean-shaped" solitons. Each class differentiates from the other by distinctive features of their energy evolution, shape of initial conditions, as well as domain of existence of parameters.

We investigate soliton formation and interactions in the non-integrable Gardner-Ostrovsky equation (alias the rotation-modified Gardner equation) with anomalous (positive) dispersion. This equation is relevant to the description of wave dynamics in plasmas and other dispersive media. We show that a special class of solitons, characterized by zero total mass and non-monotonic asymptotics, can emerge from pulse-like initial perturbations. Depending on the initial conditions, these solitons may form regular amplitude-ordered trains, irregular nonstationary ensembles of interacting localized structures, or stationary propagating multi-soliton complexes. We further demonstrate that soliton interactions in the Gardner-Ostrovsky equation are inherently inelastic, leading to the emergence of a dominant "soliton-champion" in closed systems. For example, under periodic boundary conditions, only the largest-amplitude soliton ultimately survives, while all smaller solitons are gradually eliminated through successive interactions. Despite their inelastic interactions, Gardner-Ostrovsky solitons exhibit remarkable robustness.

Numerical simulations of the complex cubic-quintic Ginzburg-Landau equation (CCQGLE), a canonical equation governing the weakly nonlinear behavior of dissipative systems in a wide variety of disciplines, reveal five entirely novel classes of pulse or solitary waves solutions, viz. pulsating, creeping, snaking, erupting, and chaotical solitons. Here, we develop a theoretical framework for analyzing the full spatio-temporal structure of one class of dissipative solution (snaking soliton) of the CCQGLE using the variational approximation technique and the dynamical systems theory. The qualitative behavior of the snaking soliton is investigated using the numerical simulations of (a) the full nonlinear complex partial differential equation and (b) a system of three ordinary differential equations resulting from the variational approximation.

We present an open-source Python implementation of an idealized high-order pseudo-spectral solver for the one-dimensional nonlinear Schr\"odinger equation (NLSE). The solver combines Fourier spectral spatial discretization with an adaptive eighth-order Dormand-Prince time integration scheme to achieve machine-precision conservation of mass and near-perfect preservation of momentum and energy for smooth solutions. The implementation accurately reproduces fundamental NLSE phenomena including soliton collisions with analytically predicted phase shifts, Akhmediev breather dynamics, and the development of modulation instability from noisy initial conditions. Four canonical test cases validate the numerical scheme: single soliton propagation, two-soliton elastic collision, breather evolution, and noise-seeded modulation instability. The solver employs a 2/3 dealiasing rule with exponential filtering to prevent aliasing errors from the cubic nonlinearity. Statistical analysis using Shannon, R\'enyi, and Tsallis entropies quantifies the spatio-temporal complexity of solutions, while phase space representations reveal the underlying coherence structure. The implementation prioritizes code transparency and educational accessibility over computational performance, providing a valuable pedagogical tool for exploring nonlinear wave dynamics. Complete source code, documentation, and example configurations are freely available, enabling reproducible computational experiments across diverse physical contexts where the NLSE governs wave evolution, including nonlinear optics, Bose-Einstein condensates, and ocean surface waves.

The Korteweg-de Vries (KdV) equation serves as a foundational model in nonlinear wave physics, describing the balance between dispersive spreading and nonlinear steepening that gives rise to solitons. This article introduces sangkuriang, an open-source Python library for solving this equation using Fourier pseudo-spectral spatial discretization coupled with adaptive high-order time integration. The implementation leverages just-in-time (JIT) compilation for computational efficiency while maintaining accessibility for instructional purposes. Validation encompasses progressively complex scenarios including isolated soliton propagation, symmetric two-wave configurations, overtaking collisions between waves of differing amplitudes, and three-body interactions. Conservation of the classical invariants is monitored throughout, with deviations remaining small across all test cases. Measured soliton velocities conform closely to theoretical predictions based on the amplitude-velocity relationship characteristic of integrable systems. Complementary diagnostics drawn from information theory and recurrence analysis confirm that computed solutions preserve the regular phase-space structure expected for completely integrable dynamics. The solver outputs data in standard scientific formats compatible with common analysis tools and generates visualizations of spatiotemporal wave evolution. By combining numerical accuracy with practical accessibility on modest computational resources, sangkuriang offers a platform suitable for both classroom demonstrations of nonlinear wave phenomena and exploratory research into soliton dynamics.

We present a model for chi^{(2)} waveguides accounting for three modes, two of which make an avoided crossing at the second harmonic wavelength. We introduce two linearly coupled pure modes and adjust the coupling to replicate the waveguide dispersion near the avoided crossing. Analysis of the nonlinear system reveals continuous wave (CW) solutions across much of the parameter-space and prevalence of its modulational instability. We also predict the existence of the avoided-crossing solitons, and study peculiarities of their dynamics and spectral properties, which include formation of a pedestal in the pulse tails and associated pronounced spectral peaks. Mapping these solitons onto the linear dispersion diagrams, we make connections between their existence and CW existence and stability. We also simulate the two-color soliton generation from a single frequency pump pulse to back up its formation and stability properties.

Using a variational formulation for partial differential equations (PDEs) combined with numerical simulations on ordinary differential equations (ODEs), we find two categories (pulses and snakes) of dissipative solitons, and analyze the dependence of both their shape and stability on the physical parameters of the cubic-quintic Ginzburg-Landau equation (CGLE). In contrast to the regular solitary waves investigated in numerous integrable and non-integrable systems over the last three decades, these dissipative solitons are not stationary in time. Rather, they are spatially confined pulse-type structures whose envelopes exhibit complicated temporal dynamics. Numerical simulations reveal very interesting bifurcations sequences as the parameters of the CGLE are varied. Our predictions on the variation of the soliton amplitude, width, position, speed and phase of the solutions using the variational formulation agree with simulation results.

We present an approach to the bright soliton solution of the NLS equation from the standpoint of introducing a constant potential term in the equation. We discuss a `nongauge' bright soliton for which both the envelope and the phase depend only on the traveling variable. We also construct a family of generalized NLS equations with solitonic sech^p solutions in the traveling variable and find an exact equivalence with other nonlinear equations, such as the Korteveg-de Vries and Benjamin-Bona-Mahony equations when p=2

The coupled cubic nonlinear Schrödinger (CNLS) equations are used to study modulational instabilities of a pair of nonlinearly interacting two-dimensional waves in deep water. It has been shown that the full dynamics of these interacting waves gives rise to localized large-amplitude wavepackets (wave focusing). In this short letter we attempt to verify this result numerically using a Fourier spectral method for the CNLS equations.

We provide conditions for existence of hyperbolic, unbounded periodic and elliptic solutions in terms of Weierstrass wp functions of both third and fifth-order KdV--BBM (Korteweg-de Vries--Benjamin, Bona \& Mahony) regularized long wave equation. An analysis for the initial value problem is developed together with a local and global well-posedness theory for the third-order KdV--BBM equation. Traveling wave reduction is used together with zero boundary conditions to yield solitons and periodic unbounded solutions, while for nonzero boundary conditions we find solutions in terms of Weierstrass elliptic wp functions. For the fifth-order KdV--BBM equation we show that a parameter γ=1{12}, for which the equation has a Hamiltonian, represents a restriction for which there are constraint curves that never intersect a region of unbounded solitary waves, which in turn shows that only dark or bright solitons and no unbounded solutions exist. Motivated by the lack of a Hamiltonian structure for γneq1{12} we develop H^k bounds, and we show for the non Hamiltonian system that dark and bright solitons coexist together with unbounded periodic solutions. For nonzero boundary conditions, due to the complexity of the nonlinear algebraic system of coefficients of the elliptic equation we construct Weierstrass solutions for a particular set of parameters only.

Achieving robust high-dimensional solitons both in optics and Bose-Einstein condensates remains a great challenge, due to unpredictable high-order perturbations. Here we propose and demonstrate a novel topological spatiotemporal mode-locked fiber laser empowered by mode-division multiplexing, for precisely controlling of the spatiotemporal solitons. By introducing time delays between the individual transverse mode paths, a topological spatiotemporal lattice is formed in the multimode cavity, which in turn excites periodic spatiotemporal solitons. These solitons are stable against the perturbations, attributed to periodic potentials induced by the nonreciprocal dissipative coupling interactions. Moreover, under the thermodynamic limit, the solitons condense into nested-type spatiotemporal solitons, a phenomenon analogous to the Bose-Einstein condensate of matter waves. Our results indicate that topological spatiotemporal mode-locked lasers not only constitute a controllable source of robust spatiotemporal light bullets, but also provide a valuable platform for exploring nonlinear high-dimensional soliton dynamics and non-Hermitian topological physics.

In this manuscript we set up the direct and inverse scattering problems for step-like traveling-wave solutions of the nonlinear Schrödinger equation. Specifically, we consider initial data u(x,0) satisfying u(x,0)to u_0^ell(x) as xto-infty and u(x,0)to u_0^r(x) as xto+infty, where u_0^ell(x) and u_0^r(x) are elliptic traveling waves. Under suitable assumptions on the initial data we formulate the direct scattering problem and establish analytic properties of the scattering data. We then formulate the inverse problem as a Riemann--Hilbert problem and prove its solvability. Finally, we observe that this Riemann--Hilbert formulation is a special case of the one arising for full soliton-gas initial data.

External potentials break the translational invariance of nonlinear Schrödinger (NLS) equations and generally convert a continuous family of standing wave solutions into isolated pinned states. We study this mechanism for the prototypical case of a kink in the NLS equation with competing focusing cubic and defocusing quintic nonlinearities. The Lyapunov--Schmidt reduction method is used to prove that the location s_0 of the pinned kinks is selected near roots of an effective potential V_{rm eff}(s). We then determine how the translational zero eigenvalue of the Jacobian operator splits under the perturbation and prove that the sign of V_{rm eff}'(s_0) distinguishes eigendirections of positive and negative energy. We further show that if the Jacobian operator admits a negative eigenvalue, then the stability problem contains a real unstable eigenvalue. The derivation of the asymptotic expansion for the unstable eigenvalue is nonstandard because the zero eigenvalue is embedded in the continuous spectrum of the linearized operator. Finally, we show that the persistence of kinks lacks parity symmetry even if the external potential is odd. Numerical continuation, spectral computations with domains adapted to the slowly decaying eigenfunctions, and direct dynamical simulations of the time evolution corroborate the analytical predictions for different representative examples.

The Rayleigh-Plateau instability predicts that infinitesimal perturbations of a liquid jet grow only for kR < 1, with kR = 1 marking the classical stability boundary. Here, we show that femtosecond-laser-induced cavitation generates a localized, impulsive recoil perturbation on a flowing liquid jet, enabling a two-dimensional map of jet breakup across the (kR, epsilon/R) plane. The perturbation's spacing (wavelength) and strength (amplitude) are set independently by the laser repetition rate and pulse energy, respectively. At low amplitudes, tuning the repetition rate recovers the classical Rayleigh-Plateau dispersion relation and its kR = 1 cutoff, while increasing the wavelength enables controlled generation of monodisperse droplet trains, bidisperse droplet populations, and ultimately isolated single droplets. Along the amplitude axis, low pulse energies produce a seed that amplifies only the imposed wavenumber, whereas above a critical energy Ec, the seed broadens to excite many wavenumbers simultaneously. At still higher amplitudes, we demonstrate deterministic breakup for kR > 1, extending to kR = 1.9, with laser-synchronized droplet production rates up to 0.2 MHz. These results establish localized femtosecond-laser-induced cavitation as a route to controlled Rayleigh-Plateau breakup both within and beyond the classical linear-stability boundary.

In this paper we review the physical relevance of a Korteweg-de Vries (KdV) equation with higher-order dispersion terms which is used in the applied sciences and engineering. We also present exact traveling wave solutions to this generalized KdV equation using an elliptic function method which can be readily applied to any scalar evolution or wave equation with polynomial terms involving only odd derivatives. We show that the generalized KdV equation still supports hump-shaped solitary waves as well as cnoidal wave solutions provided that the coefficients satisfy specific algebraic constraints. Analytical solutions in closed form serve as benchmarks for numerical solvers or comparison with experimental data. They often correspond to homoclinic orbits in the phase space and serve as separatrices of stable and unstable regions. Some of the solutions presented in this paper correct, complement, and illustrate results previously reported in the literature, while others are novel.

Many machine learning algorithms rely on the Fourier transform, which is also the inner core of many analog computing devices and, specifically, of optical computing. The square of the modulus of the Fourier spectrum naturally maps onto pairwise interactions, as in the renowned Ising model, making it applicable to combinatorial optimization. Here we show how the nonlinear Fourier transform, related to integrable nonlinear partial differential equations, significantly widens the computational capability and also enables a wide variety of computing gates. We study specifically the case of the nonlinear Schrödinger equation; we introduce an input composed by ``tokens'', i.e., batches of information encoded in a multi-step potential, and we show how the nonlinear spectral computing arising in the wave propagation not only generalizes combinatorial optimization to higher order problems, but also realizes boolean gates. In the weakly nonlinear regime, nonlinearity generates higher-order interactions as perturbative corrections to the pairwise terms. In the solitonic regime, the discrete nonlinear spectrum performs digital computation through soliton generation. By using four input tokens, we demonstrate the programmability of all 256 Boolean functions of three variables, including 3-SAT instances. Our results show that integrable nonlinear propagation provides an unprecedented single physical platform for higher-order and logical computation, opening the way to a rich variety of new algorithms for analog computing, machine learning, and classical and quantum information processing.

Continuous cellular automata are rocketing in popularity, yet developing a theoretical understanding of their behaviour remains a challenge. In the case of Lenia, a few fundamental open problems include determining what exactly constitutes a soliton, what is the overall structure of the parameter space, and where do the solitons occur in it. In this abstract, we present a new method to automatically classify Lenia systems into four qualitatively different dynamical classes. This allows us to detect moving solitons, and to provide an interactive visualization of Lenia's parameter space structure on our website https://lenia-explorer.vercel.app/. The results shed new light on the above-mentioned questions and lead to several observations: the existence of new soliton families for parameters where they were not believed to exist, or the universality of the phase space structure across various kernels.

We study a one dimensional nonlinear Klein-Gordon equation with a potential, a localized quadratic nonlinearity, and a non-localized quartic nonlinearity. We assume that the linearized operator has a single discrete eigenvalue below the continuous spectrum, corresponding to an ''internal mode''. This eigenvalue generates localized, time-periodic solutions for the linear equation. Assuming the natural Fermi Golden Rule, we give a global description of the dynamics of the amplitude of the internal mode through the radiation damping mechanism, and prove scattering for the radiation. Our approach is based on the distorted Fourier transform and a collection of refined dispersive decay and smoothing estimates. The corresponding cubic problem is closely related to the question of global-in-space dispersive asymptotics for perturbations of kink solutions in the classical ϕ^4 model, which remains a challenging open question. From the perspective of nonlinear estimates, a quartic (non-localized) interaction is essentially sharp for our analysis, and closing the estimates with such a low degree nonlinear term requires a careful and delicate treatment. As a direct application of our analysis, we prove a result on the asymptotic stability of kinks for classes of scalar field models that can be seen as perturbations of the ϕ^4 model.

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, \alpha < 2 (the normal kinetic energy corresponds to \alpha = 2). The SOC terms, with strength \lambda, 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 < \alpha < 2, where the supercritical collapse does not admit the existence of stable solitons in the absence of the SOC. At \lambda --> 0, amplitudes of these solitons vanish as (\lambda)^{1/(\alpha - 1)}.

We establish the nonlinear stability of solitary waves (solitons) and periodic traveling wave solutions (cnoidal waves) for a Korteweg-de Vries (KdV) equation which includes a fifth order dispersive term. The traveling wave solutions which yield solitons for zero boundary conditions and wave-trains of cnoidal waves for nonzero boundary conditions are analyzed using stability theorems, which rely on the positivity properties of the Fourier transforms. We show that all families of solutions considered here are (orbitally) stable.

This paper derives the second-order positive Burgers' equation from the standard Burgers' hierarchy to explore its complete integrability and exact analytical solutions. We construct this higher-order nonlinear evolution equation by systematically applying the recursion operator to the classical Burgers' equation. Expanding on this structural framework, we then derive the explicit third-order equation and use Complete Bell Polynomials to generalize the n-th order hierarchy. Through the Cole-Hopf transformation, the second-order nonlinear equation rigorously maps to the linear third-order dispersion equation. We establish complete integrability by explicitly formulating the scalar Lax pair (zero-curvature representation), which allows us to directly derive the associated differential and algebraic Darboux transformations. To systematically classify explicit stationary, time-dependent traveling wave, and self-similar profiles, we apply several techniques to the linear domain: separation of variables, the Hirota perturbation method, traveling wave reduction via the Complete Discrimination System for Polynomial Method (CDSPM), and Lie similarity reduction. Crucially, our analysis demonstrates the exact finite truncation of the Hirota perturbation series. We conclude by outlining how this equation impacts the theoretical understanding of fluid dynamics, nonlinear transport phenomena, and higher-order wave propagation.

Generation of high-fidelity amplified pulse bursts with a regular interpulse interval yields, in the spectral domain, an equidistant pattern of narrowband spectral modes, similar to frequency combs produced by cw mode-locked lasers, but with greatly increased pulse energy. Despite their great potential for nonlinear spectroscopy, material processing, etc., such long frequency-stable bursts are difficult to generate and amplify because of prominent temporal intensity modulation even after strong dispersive pulse stretching. This study presents a burst generation method based on a master-oscillator regenerative-amplifier system that allows for chirped-pulse amplification (CPA) with high scalability in pulse number. A gradual smoothing of temporal intensity profiles at an increasing number of pulses is discovered, demonstrating an unexpected recovery of the CPA performance at terahertz (THz) intraburst repetition rates. In consequence, a self-referenced stable burst spectral peak structure with megahertz (MHz) peak width is generated, without risk of amplifier damage caused by interference of chirped pulses. This result eliminates limitations in burst amplification and paves the way for advancements in ultrashort-pulse burst technology, particularly for its use in nonlinear optical applications.

The basic framework of the superfluid vortex model for pulsar glitches, though, is well accepted; there is a lack of consensus on the possible trigger mechanism responsible for the simultaneous release of a large number (sim 10^{17}) of superfluid vortices from the inner crust. Here, we propose a simple trigger mechanism to explain such catastrophic events of vortex unpinning. We treat a superfluid vortex line as a classical massive straight string with well-defined string tension stretching along the rotation axis of pulsars. The crustquake-induced lattice vibration of the inner crust can act as a driving force for the transverse oscillation of the string. Such forced oscillation near resonance causes the bending of the vortex lines, disturbing their equilibrium configuration and resulting in the unpinning of vortices. We consider unpinning from the inner crust's so-called {\it strong (nuclear)} pinning region, where the vortices are likely pinned to the nuclear sites. We also comment on vortex unpinning from the interstitial pinning region of the inner crust. We sense that unifying crustquake with the superfluid vortex model can naturally explain the cause of large-scale vortex unpinning and generation of large-size pulsar glitches.

Repeated impacts of a droplet on a vertically forced bath can produce a growing Faraday-wave response under forcing conditions for which the corresponding initially undisturbed bath remains linearly stable. We investigate this mechanism using axisymmetric two-phase direct numerical simulations over a wide range of Reynolds and Weber numbers. The resulting phenomenological regime map identifies quasi-periodic, decaying, faint, and unstable-bath responses. Along the parameter path considered here, the transition to the unstable-bath regime occurs over a narrow range of Weber numbers and varies only weakly with Reynolds number. In the unstable regime, the bath develops a coherent subharmonic response at half the forcing frequency, while disturbances in the decaying regime remain broadband and weaken between successive rebounds. Bath-only simulations show that a finite-amplitude perturbation representative of droplet impact can seed the confined radial mode susceptible to parametric forcing, initiating a growing subharmonic response below the linear instability threshold of the initially flat bath. A reduced capillary-dominated model provides a scaling estimate for the impact-equivalent perturbation amplitude, and bath-only simulations initialized with these perturbations reproduce the bounded-to-growing transition observed in the coupled simulations. Floquet analysis of the finite cylindrical bath identifies the first axisymmetric Neumann mode as having the lowest subharmonic instability threshold. A time-resolved Bessel decomposition independently shows that the growing bath response is dominated by the same radial mode, while the droplet-induced subharmonic response projects predominantly onto this structure. These results identify a finite-amplitude route by which localized droplet impacts can access Faraday-wave dynamics below the linear onset threshold of an initially undisturbed bath.

Traveling wave solutions to Kawahara equation (KE), transmission line (TL), and Korteweg-de Vries (KdV) equation are found by using an elliptic function method which is more general than the tanh-method. The method works by assuming that a polynomial ansatz satisfies a Weierstrass equation, and has two advantages: first, it reduces the number of terms in the ansatz by an order of two, and second, it uses Weierstrass functions which satisfy an elliptic equation for the dependent variable instead of the hyperbolic tangent functions which only satisfy the Riccati equation with constant coefficients. When the polynomial ansatz in the traveling wave variable is of first order, the equation reduces to the KdV equation with only a cubic dispersion term, while for the KE which includes a fifth order dispersion term the polynomial ansatz must necessary be of quadratic type. By solving the elliptic equation with coefficients that depend on the boundary conditions, velocity of the traveling waves, nonlinear strength, and dispersion coefficients, in the case of KdV equation we find the well-known solitary waves (solitons) for zero boundary conditions, as well as wave-trains of cnoidal waves for nonzero boundary conditions. Both solutions are either compressive (bright) or rarefactive (dark), and either propagate to the left or right with arbitrary velocity.

Singularity Theory is used to comprehensively investigate the bifurcations of the steady-states of the traveling wave ODEs of the cubic-quintic Ginzburg-Landau equation (CGLE). These correspond to plane waves of the PDE. In addition to the most general situation, we also derive the degeneracy conditions on the eight coefficients of the CGLE under which the equation for the steady states assumes each of the possible quartic (the quartic fold and an unnamed form), cubic (the pitchfork and the winged cusp), and quadratic (four possible cases) normal forms for singularities of codimension up to three. Since the actual governing equations are employed, all results are globally valid, and not just of local applicability. In each case, the recognition problem for the unfolded singularity is treated. The transition varieties, i.e. the hysteresis, isola, and double limit curves are presented for each normal form. For both the most general case, as well as for various combinations of coefficients relevant to the particular cases, the bifurcations curves are mapped out in the various regions of parameter space delimited by these varieties. The multiplicities and interactions of the plane wave solutions are then comprehensively deduced from the bifurcation plots in each regime, and include features such as regimes of hysteresis among co-existing states, domains featuring more than one interval of hysteresis, and isola behavior featuring dynamics unrelated to the primary solution branch in limited ranges of parameter space.

Harnessing integrability at arbitrary interaction strength, we build long-lived and spatially localized soliton-like excitations in the 1D Bose gas with contact repulsion (Lieb-Liniger model). We show how these density profiles, built from selective superpositions of multi-hole eigenstates, display long lifetimes tunable through parameters like dispersion width and number of holes. We further illustrate a scaling protocol that in principle enables unbounded increase in lifetimes in the thermodynamic limit. Our results showcase classical-like excitations originating from unitary quantum dynamics, away from the semi-classical regime, and in regions where interactions cannot be treated perturbatively.

In this paper, we propose an alternating numerical scheme for constructing quasi-periodic solutions to a class of nonlinear partial differential equations (PDEs) with bounded perturbations, including the one-dimensional nonlinear Schrödinger (NLS) equations and nonlinear wave (NLW) equations under periodic and Dirichlet boundary conditions. A key observation concerning the resonant Q-equations allows us to establish directly a local diffeomorphism between the amplitudes and the drifted frequencies. This renders the proof considerably cleaner, since it avoids the Birkhoff normal form and the associated coordinate transformations required in classical KAM theory for PDEs.Unlike in the finite-dimensional setting, the lattice boxes must also include the spatial index, so that the multiscale induction has to be carried out in one additional dimension. This introduces two difficulties. First, the measure estimate for the excluded set of ``bad'' frequencies must account for the spatial index. For both the NLS and NLW equations, the linear frequencies and their differences are integers or close to integers; this allows the spatial dependence to be controlled through the modulus of the temporal index, so that the excluded measure remains under control. Second, a singular set may contain sites with the same temporal index but widely separated spatial indices, which obstructs the ``inversion implies localization'' process. For the one-dimensional NLS and NLW equations, however, such sites occur only in pairs with opposite spatial indices. They can therefore be separated by small boxes, after which the global inverse is obtained from the local inverses via the resolvent identity. Numerical experiments illustrate the scheme for nonlinear traveling and standing waves, exhibiting both the computed quasi-periodic solutions and the convergence of the iteration.

In this work, we consider the generalized Benjamin-Bona-Mahony equation $partial_t u+partial_x u+partial_x( |u|^pu)-partial_t partial_x^{2}u=0, quad(t,x) in R times R, with p>4. This equation has the traveling wave solutions \phi_{c}(x-ct), for any frequency c>1. It has been proved by Souganidis and Strauss Strauss-1990 that, there exists a number c_{0}(p)>1, such that solitary waves \phi_{c}(x-ct) with 1<c<c_{0}(p) is orbitally unstable, while for c>c_{0}(p), \phi_{c}(x-ct) is orbitally stable. The linear exponential instability in the former case was further proved by Pego and Weinstein Pego-1991-eigenvalue. In this paper, we prove the orbital instability in the critical case c=c_{0}(p)$.

In this paper, we use a traveling wave reduction or a so-called spatial approximation to comprehensively investigate periodic and solitary wave solutions of the modified Benjamin, Bona & Mahony equation (BBM) to include both dissipative and dispersive effects of viscous boundary layers. Under certain circumstances that depend on the traveling wave velocity, classes of periodic and solitary wave like solutions are obtained in terms of Jacobi elliptic functions. An ad-hoc theory based on the dissipative term is presented, in which we have found a set of solutions in terms of an implicit function. Using dynamical systems theory we prove that the solutions of BBMv experience a transcritical bifurcation for a certain velocity of the traveling wave. Finally, we present qualitative numerical results.

Core collapse supernovae are among the most energetic astrophysical events in the Universe. Despite huge efforts on understanding the main ingredients triggering such explosions, we still lack of compelling evidences for the precise mechanism driving those phenomena. They are expected to produce gravitational waves due to asymmetric mass motions in the collapsing core, and emit in the meanwhile neutrinos as a result of the interactions in their high-density environment. The combination of these two cosmic messengers can provide a unique probe to study the inner engine of these processes and unveil the explosion mechanism. Among the possible detectable signature, standing accretion shock instabilities (SASI) are particularly relevant in this context as they establish a direct connection between gravitational wave emission and the outcoming neutrino flux. In this work, Hilbert-Huang transform is applied to a selected sample of 3D numerical simulations, with the aim of identifying SASI contribution and extract its instantaneous frequency. The performance of the method is evaluated in the context of Einstein Telescope.

We use a simple method that leads to the integrals involved in obtaining the traveling wave solutions of wave equations with one and two exponential nonlinearities. When the constant term in the integrand is zero, implicit solutions in terms of hypergeometric functions are obtained while when that term is nonzero all the basic traveling wave solutions of Liouville, Tzitzeica and their variants, as well as sine/sinh-Gordon equations with important applications in the phenomenology of nonlinear physics and dynamical systems are found through a detailed study of the corresponding elliptic equations

In this paper, Liouville soliton surfaces based on some soliton solutions of the Liouville equation are constructed and displayed graphically, including some of those corresponding to Darboux-transformed counterparts. We find that the Liouville soliton surfaces are centroaffine surfaces of Tzitzeica type and their centroaffine invariant can be expressed in terms of the Hamiltonian. The traveling wave solutions to Liouville equation from which these soliton surfaces stem are also obtained through a modified variation of parameters method which is shown to lead to elliptic functions solution method.

In this paper, we study the long time dynamics of small perturbations of static domain walls in weighted Sobolev spaces for the one dimensional Schrödinger map with constant Dzyaloshinskii-Moriya interaction. The perturbation equation contains derivative nonlinearities, and its linearized operator has a threshold resonance. The results are twofold. (i) We prove the asymptotic stability of the domain wall modulo the translation and rotation symmetries. More precisely, the corresponding modulation parameters converge as time tends to infinity, while the radiation decays at the sharp rate t^{-1/2}. (ii) We prove that the radiation exhibits modified scattering, with an explicit logarithmic phase correction generated by the long range cubic interaction. The proof develops geometric reduction to overcome the loss of derivatives and exploits structural cancellations in the nonlinear interactions.

We prove that all excited states -sin(nx), ngeq2, of the De Gregorio equation on the torus are nonlinearly unstable in the sense of Hadamard. More precisely, perturbations arbitrarily small in H^2 leave a fixed H^1 neighborhood of the state. For every ngeq2, the linearized operator has a real eigenvalue λ_nin(0.56,0.6). To prove this, we transform the eigenvalue problem into a Fuchsian equation on the complex disk and obtain a matching condition for the eigenvalue. This condition is verified rigorously using interval arithmetic. The nonlinear instability follows from a bootstrap argument in the spirit of Guo, Hallstrom, and Spirn.

The transition of an impulsively excited kink oscillation of a solar coronal loop to an oscillation with a stationary amplitude, i.e., the damping pattern, is determined using the low-dimensional self-oscillation model. In the model, the decayless kink oscillations are sustained by the interaction of the oscillating loop with an external quasi-steady flow. The analytical solution is based on the assumption that the combined effect of the effective dissipation, for example, by resonant absorption, and interaction with an external flow, is weak. The effect is characterised by a dimensionless coupling parameter. The damping pattern is found to depend upon the initial amplitude and the coupling parameter. The approximate expression shows a good agreement with a numerical solution of the self-oscillation equation. The plausibility of the established damping pattern is demonstrated by an observational example. Notably, the damping pattern is not exponential, and the characteristic decay time is different from the time determined by the traditionally used exponential damping fit. Implications of this finding for seismology of the solar coronal plasmas are discussed. In particular, it is suggested that a very rapid, in less than the oscillation period, decay of the oscillation to the stationary level, achieved for larger values of the coupling parameter, can explain the relative rareness of the kink oscillation events.

In this work, we compute multi-field core and halo properties in wave Dark Matter models. We focus on the case where Dark Matter consists of two light (real) scalars, interacting gravitationally. As in the single-field Ultra Light Dark Matter (ULDM) case, the scalar field behaves as a coherent BEC with a definite ground state (at fixed total mass), often referred to in the literature as a gravitational soliton. We establish an efficient algorithm to find the ground and excited states of such two-field systems. We then use simulations to investigate the gravitational collapse and virialization, starting from different initial conditions, into solitons and surrounding halo. As in the single-field case, a virialized halo forms with a gravitational soliton (ground state) at the center. We find some evidence for an empirical relation between the soliton mass and energy and those of the host halo. We use this to then find a numerical relation between the properties of the two. Finally, we use this to address the issue of alleviating some of the tensions that single-field ULDM has with observational data, in particular, the issue of how a galaxy's core and radius are related. We find that if galaxies of different masses have similar percentages of the two species, then the core-radius scaling tension is not addressed. However, more general possibilities occur if the relative abundance of species in each halo correlates with the total mass of the galaxy. If this is the case, the model predicts several other phenomenological signatures.