CyberSec Research

We obtain and investigate theoretically a broad family of stable and unstable time-periodic orbits-oscillating Turing rolls (OTR)-in the Lugiato-Lefever model of optical cavities. Using the dynamical systems tools developed in fluid dynamics, we access the OTR solution branches in parameter space and elucidate their bifurcation structure. By tracking these exact invariant solutions deeply into the chaotic region of the modulation instability, we connect the main dynamical regimes of the Lugiato-Lefever model: continuous waves, Turing rolls, solitons, and breathers, which completes the classical phase diagram of the optical cavity. We then demonstrate that the OTR periodic orbits play a fundamental role as elementary building blocks in the regime of the intracavity field transition from stable Turing rolls to fully developed turbulent regimes. Depending on the cavity size, we observe that the chaotic intracavity field driven by modulation instability displays either spatiotemporal or purely temporal intermittancy between chaotic dynamics and different families of the OTR solutions, exhibiting locally the distinctive wave patterns and large amplitude peaks. This opens avenues for a theoretical description of optical turbulence within the dynamical systems framework.

We are surrounded by spatio-temporal patterns resulting from the interaction of the numerous basic units constituting natural or human-made systems. In presence of diffusive-like coupling, Turing theory has been largely applied to explain the formation of such self-organized motifs both on continuous domains or networked systems, where reactions occur in the nodes and the available links are used for species to diffuse. In many relevant applications, those links are not static, as very often assumed, but evolve in time and more importantly they adapt their weights to the states of the nodes. In this work, we make one step forward and we provide a general theory to prove the validity of Turing idea in the case of adaptive symmetric networks with positive weights. The conditions for the emergence of Turing instability rely on the spectral property of the Laplace matrix and the model parameters, thus strengthening the interplay between dynamics and network topology. A rich variety of patterns are presented by using two prototype models of nonlinear dynamical systems, the Brusselator and the FitzHugh-Nagumo model. Because many empirical networks adapt to changes in the system states, our results pave the way for a thorough understanding of self-organization in real-world systems.

Ventricular arrhythmias, like ventricular tachycardia (VT) and ventricular fibrillation (VF), precipitate sudden cardiac death (SCD), which is the leading cause of mortality in the industrialised world. Thus, the elimination of VT and VF is a problem of paramount importance, which is studied experimentally, theoretically, and numerically. Numerical studies use partial-differential-equation models, for cardiac tissue, which admit solutions with spiral- or broken-spiral-wave solutions that are the mathematical counterparts of VT and VF. In silico investigations of such mathematical models of cardiac tissue allow us not only to explore the properties of such spiral-wave turbulence, but also to develop mathematical analogues of low-amplitude defibrillation by the application of currents that can eliminate spiral waves. We develop an efficient deep-neural-network U-Net-based method for the control of spiral-wave turbulence in mathematical models of cardiac tissue. Specifically, we use the simple, two-variable Aliev-Panfilov and the ionically realistic TP06 mathematical models to show that the lower the correlation length {\xi} for spiral-turbulence patterns, the easier it is to eliminate them by the application of control currents on a mesh electrode. We then use spiral-turbulence patterns from the TP06 model to train a U-Net to predict the sodium current, which is most prominent along thin lines that track the propagating front of a spiral wave. We apply currents, in the vicinities of the predicted sodium-current lines to eliminate spiral waves efficiently. The amplitudes of these currents are adjusted automatically, so that they are small when {\xi} is large and vice versa. We show that our U-Net-aided elimination of spiral-wave turbulence is superior to earlier methods.

We present a COMSOL Multiphysics implementation of a continuum model for directed cell migration, a key mechanism underlying tissue self-organization and morphogenesis. The model is formulated as a partial integro-differential equation (PIDE), combining random motility with non-local, density-dependent guidance cues to capture phenomena such as cell sorting and aggregation. Our framework supports simulations in one, two, and three dimensions, with both zero-flux and periodic boundary conditions, and can be reformulated in a Lagrangian setting to efficiently handle tissue growth and domain deformation. We demonstrate that COMSOL Multiphysics enables a flexible and accessible implementation of PIDEs, providing a generalizable platform for studying collective cell behavior and pattern formation in complex biological contexts.

Pattern formation often occurs in confined systems, yet how boundaries shape patterning dynamics is unclear. We develop techniques to analyze confinement effects in nonlocal advection-diffusion equations, which generically capture the collective dynamics of active self-attracting particles. We identify a sequence of size-controlled transitions that generate characteristic slow modes, leading to exponential increase of patterning timescales. Experimental measurements of multicellular dynamics confirm our predictions.

While real-valued solutions of the Korteweg--de Vries (KdV) equation have been studied extensively over the past 50 years, much less attention has been devoted to solution behaviour in the complex plane. Here we consider the analytic continuation of real solutions of KdV and investigate the role that complex-plane singularities play in early-time solutions on the real line. We apply techniques of exponential asymptotics to derive the small-time behaviour for dispersive waves that propagate in one direction, and demonstrate how the amplitude, wavelength and speed of these waves depend on the strength and location of double-pole singularities of the initial condition in the complex plane. Using matched asymptotic expansions in the limit $t\rightarrow 0^+$, we show how complex singularities of the time-dependent solution of the KdV equation emerge from these double-pole singularities. Generically, their speed as they move from their initial position is of $\mathcal{O}(t^{-2/3})$, while the direction in which these singularities propagate initially is dictated by a Painlev\'{e} II (P$_{\mathrm{II}}$) problem with decreasing tritronqu\'{e}e solutions. The well-known $N$-soliton solutions of KdV correspond to rational solutions of P$_{\mathrm{II}}$ with a finite number of singularities; otherwise, we postulate that infinitely many complex-plane singularities of KdV solutions are born at each double-pole singularity of the initial condition. We also provide asymptotic results for some non-generic cases in which singularities propagate more slowly than in the generic case. Our study makes progress towards the goal of providing a complete description of KdV solutions in the complex plane and, in turn, of relating this behaviour to the solution on the real line.

We present an analytical and numerical investigation into the phenomenon of filter-induced modulation instability in passive hybrid optical resonators exhibiting quadratic and cubic nonlinearity. We show that asymmetric spectral losses, with respect to the continuous-wave solution frequency, can trigger sideband amplification in the normal dispersion regime. We calculate the parametric gain and demonstrate the associated pattern formation process. We furthermore show how this parametric process can be exploited to generate optical frequency combs with tunable repetition rate.

In the present work we introduce and explore a technique for the efficient removal of vortices from an atomic Bose-Einstein condensate, through the application and subsequent removal of a one-dimensional optical lattice. We showcase a prototypical experimental realization of the technique that motivates a detailed theoretical study of vortex removal mechanisms. Through simulations of the condensate dynamics during application of the optical lattice, we also discover a vortex removal mechanism that arises in narrow, optical-lattice-induced atomic density channels for which the channel width is on the order of the nominal vortex core size and healing length. This mechanism involves the density profile typically associated with a vortex core spatially separating from the phase singularity associated with the vortex. By analyzing numerical experiments covering a wide range of variations of the optical lattice amplitude and fringe periodicity, we identify the existence of an optimal set of parameters that enables the efficient removal of all vortices from the condensate. This analysis paves the way for further studies aimed at understanding vortex dynamics in narrow channels, and adds to an experimental toolkit for working with vortices and controlling the dynamical states of condensates.

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.

We study collisions between a ferrodark soliton (FDS) and an antiFDS ($Z_2$ kinks in the spin order) in the easy-plane phase of spin-1 Bose-Einstein condensates (BEC). For a type-I pair (type-I FDS-antiFDS pair) at low incoming velocities, the pair annihilates followed by the formation of an extremely long-lived dissipative breather on a stable background, a spatially localized wave packet with out-of-phase oscillating magnetization and number densities. Periodic emissions of spin and density waves cause breather energy dissipation and we find that the breather energy decays logarithmically in time. When the incoming velocity is larger than a critical velocity at which a stationary FDS-anti FDS pair forms, a pair with finite separating velocity is reproduced. When approaching the critical velocity from below, we find that the lifetime of the stationary type-I pair shows a power-law divergence, resembling a critical behavior. In contrast, a type-II pair (type-II FDS-antiFDS pair)never annihilates and only exhibits reflection. For collisions of a mixed type FDS-antiFDS pair, as $Z_2$ kinks in the spin order, reflection occurs in the topological structure of the magnetization while the mass superfluid density profiles pass through each other, manifesting spin-mass separation.

In this paper, we study a non-integrable discrete lattice model which is a variant of an integrable discretization of the standard Hopf equation. Interestingly, a direct numerical simulation of the Riemann problem associated with such a discrete lattice shows the emergence of both the dispersive shock wave (DSW) and rarefaction wave (RW). We propose two quasi-continuum models which are represented by partial differential equations (PDEs) in order to both analytically and numerically capture the features of the DSW and RW of the lattice. Accordingly, we apply the DSW fitting method to gain important insights and provide theoretical predictions on various edge features of the DSW including the edge speed and wavenumber. Meanwhile, we analytically compute the self-similar solutions of the quasi-continuum models, which serve as the approximation of the RW of the lattice. We then conduct comparisons between these numerical and analytical results to examine the performance of the approximation of the quasi-continuum models to the discrete lattice.

We show that, in the thin-wall regime, $Q$-ball--anti-$Q$-ball collisions reveal chaotic behaviour. This is explained by the resonant energy transfer mechanism triggered by the internal modes hosted by the $Q$-balls and by the existence of {\it ephemeral} states, that is unstable, sometimes even short-lived, field configurations that appear as intermediate states. The most important examples of such states are the {\it bubble} of the false broken vacuum, which as intermediate states govern the $QQ^*$ annihilation, and the {\it charged oscillons}. The usually short-lived bubble can be dynamically temporarily stabilized, which explains their importance in the dynamics of $Q$-balls. This happens due to the excitation of massless Goldstone modes, which, exerting pressure on the bubble boundaries or being trapped as bound modes, prevent the bubble from collapsing.

We investigate the nonlinear Schr\"odinger equation on a three-edge star graph, where each edge contains a linear localized inhomogeneity in the form of a Dirac delta linear potential. Such systems are of significant interest in studying wave propagation in networked structures, with applications in, e.g., Josephson junctions. By reducing the system to a set of finite-dimensional coupled ordinary differential equations, we derive explicit conditions for the occurrence of a symmetry-breaking bifurcation in a symmetric family of solutions. This bifurcation is shown to be of the transcritical type, and we provide a precise estimate of the bifurcation point as the propagation constant, which is directly related to the solution norm, is varied. In addition to the symmetric states, we explore non-positive definite states that bifurcate from the linear solutions of the system. These states exhibit distinct characteristics and are crucial in understanding solutions of the nonlinear system. Furthermore, we analyze the typical dynamics of unstable solutions, showing their behavior and evolution over time. Our results contribute to a deeper understanding of symmetry-breaking phenomena in nonlinear systems on metric graphs and provide insights into the stability and dynamics of such solutions.

Parametric instabilities are a known feature of periodically driven dynamic systems; at particular frequencies and amplitudes of the driving modulation, the system's quasi-periodic response undergoes a frequency lock-in, leading to a periodically unstable response. Here, we demonstrate an analogous phenomenon in a purely static context. We show that the buckling patterns of an elastic beam resting on a modulated Winkler foundation display the same kind of frequency lock-in observed in dynamic systems. Through simulations and experiments, we reveal that compressed elastic strips with modulated height alternate between predictable quasi-periodic and periodic buckling modes. Our findings uncover previously unexplored analogies between structural and dynamic instabilities, highlighting how even simple elastic structures can give rise to rich and intriguing behaviors.

Filters typically play a crucial role in generating solitons by strictly confirming the working wavelength in a mode-locked laser. However, we have found that the broad smoothing filter with saturable gain also helps to regulate the pulses' operating wavelength. The effective gain, formed by the original gain spectrum and filters, shifts its center based on pulse energy, thereby affecting the pulse spectrum. A virtual mode-locked cavity is established to verify this approach in simulation. The results indicate that the spectrum changes with pulse energy, leading to a wider spectrum output. Our method paves the way for further investigations of the pulse spectrum and provides a practical approach to generating pulses at specific wavelengths.

Two-component Bose-Einstein condensates in the miscible phase can support polarization solitary waves, known as magnetic solitons. By calculating the interaction potential between two magnetic solitons, we elucidate the mechanisms and conditions for the formation of bound states -- or molecules -- and support these predictions with dynamical simulations. We analytically determine the dissociation energy of bound states consisting of two oppositely polarized solitons and find good agreement with full numerical simulations. Collisions between bound states -- either with other bound states or with individual solitons -- produce intriguing dynamics. Notably, collisions between a pair of bound states exhibit a dipole-like behavior. We anticipate that such bound states, along with their rich collision dynamics, are within reach of current experimental capabilities.

Higher-order topological insulators (HOTIs) are unique topological materials supporting edge states with the dimensionality at least by two lower than the dimensionality of the underlying structure. HOTIs were observed on lattices with different symmetries, but only in geometries, where truncation of HOTI produces a finite structure with the same order of discrete rotational symmetry as that of the unit cell, thereby setting the geometry of insulator edge. Here we experimentally demonstrate a new type of two-dimensional (2D) HOTI based on the Kekule-patterned lattice, whose order of discrete rotational symmetry differs from that of the unit cells of the constituent honeycomb lattice, with hybrid boundaries that help to produce all three possible corners that support effectively 0D corner states of topological origin, especially the one associated with spectral charge 5/6. We also show that linear corner states give rise to rich families of stable hybrid nonlinear corner states bifurcating from them in the presence of focusing nonlinearity of the material. Such new types of nonlinear corner states are observed in hybrid HOTI inscribed in transparent nonlinear dielectric using fs-laser writing technique. Our results complete the class of HOTIs and open the way to observation of topological states with new internal structure and symmetry.

In this paper, we analyse the dynamics of a pattern-forming system close to simultaneous Turing and Turing--Hopf instabilities, which have a 1:1 spatial resonance, that is, they have the same critical wave number. For this, we consider a system of coupled Swift--Hohenberg equations with dispersive terms and general, smooth nonlinearities. Close to the onset of instability, we derive a system of two coupled complex Ginzburg--Landau equations with a singular advection term as amplitude equations and justify the approximation by providing error estimates. We then construct space-time periodic solutions to the amplitude equations, as well as fast-travelling front solutions, which connect different space-time periodic states. This yields the existence of solutions to the pattern-forming system on a finite, but long time interval, which model the spatial transition between different patterns. The construction is based on geometric singular perturbation theory exploiting the fast travelling speed of the fronts. Finally, we construct global, spatially periodic solutions to the pattern-forming system by using centre manifold reduction, normal form theory and a variant of singular perturbation theory to handle fast oscillatory higher-order terms.