Papers

A one-dimensional SQUID (superconducting quantum interference device) array/metamaterial is investigated numerically with respect to its localization properties due to nonlinearity in the absence of dissipation and periodic driving. The system possesses a conserved Hamiltonian function representing its energy, and supports localized modes of the discrete breather type even in the presence of a moderately high dc flux bias. The appearance of discrete breathers in that system has been largely overlooked in literature. We find a new type of discrete breather that is asynchronous, meaning that the frequency of oscillation of the SQUID at the central breather site is different than that of the SQUIDs at the other sites of the metamaterial. Nonlinear localization is investigated by initializing the system with a single-site excitation of given amplitude (initial amplitude) for a fixed value of the coupling coefficient, while parameters such as the dc flux bias, the single-site initial excitation amplitude, and/or the SQUID can vary independently. Using the energetic participation ratio as a measure of the degree of localization, the existence of asynchronous highly localized modes and transitions between delocalized (extended) and localized modes re identified.

We study a two–dimensional discrete nonlinear Schrödinger equation with saturable nonlinearity on the lattice $\mathbb Z^2$. Using a variational approach based on the Nehari manifold, we establish the existence of nontrivial periodic ground states on finite lattices and establish the existence of exponentially localized ground states in $\ell^2(\mathbb Z^2)$. A principal result is the rigorous passage from periodic to localized states: we show that, up to lattice translations, periodic ground states converge strongly in $\ell^2(\mathbb Z^2)$ to a localized ground state as the lattice periods tend to infinity. The analysis combines variational methods, spectral properties of the discrete Laplacian, and concentration–compactness techniques adapted to the two–dimensional discrete setting. We further derive qualitative properties of the resulting solutions, including positivity and exponential localization, and establish a conditional orbital stability result within the Grillakis–Shatah–Strauss framework. Numerical computations illustrate the theoretical results and confirm the predicted convergence and localization behavior.

We develop a systematic adiabatic perturbation theory for the integrable $F=1$ spinor nonlinear Schrödinger equation under nonvanishing boundary conditions, formulated entirely within the framework of the associated Riemann–Hilbert problem. In this setting, localized nonlinear excitations are characterized by discrete spectral data consisting of a complex eigenvalue and an associated polarization vector. For a general class of small perturbations preserving the background, we derive the perturbation-induced evolution of the scattering data directly at the level of the Riemann–Hilbert problem. In the one-soliton sector, this yields a closed finite-dimensional dynamical system governing the slow evolution of the effective soliton parameters, including the spectral variables, the soliton center and phase, the residue amplitude, and the internal polarization state. The latter evolves according to a constrained dynamical equation with no scalar analogue. For localized perturbations, the modulation equations are expressed in explicit integral form in terms of the one-soliton eigenfunctions, providing a fully computable description of the dynamics. In the limit of vanishing boundary conditions, the resulting system reduces to the perturbation theory obtained by E. V. Doktorov, et al, Phys. Rev. A 77 (2008), no. 4, 043617.

We investigate the dynamics and spectral stability of traveling kink and antikink solutions in a dissipative sine–Gordon equation with two distinct fourth–order regularization mechanisms: a mixed space–time (inertial) term and a purely spatial (elliptic) term. The model includes damping, bias forcing, and higher–order dissipative effects, and is motivated by refined descriptions of fluxon dynamics in long Josephson junctions. Using a collective–coordinate reduction, we derive a Melnikov–type condition for speed selection, yielding explicit predictions for asymptotic propagation speeds, which are validated by direct numerical simulations of the full partial differential equation. Spectral stability is analyzed using Evans function techniques adapted to the singular slow–fast structure induced by the regularization. By formulating the linearized problem on a consistent–splitting domain, we show that no additional point spectrum bifurcates from the origin. Near the edges of the essential spectrum, a square–root transformation is used to resolve branch singularities and establish analyticity in a lifted spectral variable. Numerical Evans function computations near $\lambda=0$ and near the essential spectrum edges confirm the analytical results, indicating absence of unstable eigenvalues for both kink and antikink solutions.

Environmental variability is a defining feature of natural ecosystems, yet many theories of ecological stability assume static environments. Here, we develop an analytical framework for local stability in complex ecological communities subjected to periodic environmental forcing. We show that, in slowly varying environments and under a large-system approximation, the maximal Lyapunov exponent can be approximated by the time-averaged rightmost spectral edge of the instantaneous interaction matrix, yielding explicit stability criteria for large ecological communities. Within the random-matrix interaction ensembles considered here, our minimal theoretical framework predicts a topology-dependent ordering of local perturbation stability under periodic temporal modulation in high-dimensional ecological networks, which is supported by numerical simulations. Beyond the adiabatic regime, rapid environmental oscillations dynamically stabilize otherwise unstable ecosystems, revealing a high-frequency rescue effect that is absent from static theories. These results extend ecological stability theory beyond autonomous systems and provide a general framework for understanding resilience in fluctuating environments.

We derive a one-dimensional reduction for nonlinear dynamics on simplicial complexes containing both pairwise and triangular (higher-order) interactions. The effective state is defined using a mixed weight determined by the pairwise and triangular degrees of each node. The resulting reduced equation retains two structural coefficients, associated separately with the pairwise and higher-order coupling channels. A fluctuation expansion identifies the closure assumptions underlying the reduction and shows how deviations of individual node states from the effective state contribute to the approximation error. We numerically validate the proposed framework on Gene-regulatory dynamics, the double-well system, and SIS spreading. The states of the reduced model are compared with full-network simulations through coupling-parameter sweeps, steady-state branch calculations, and progressive node-removal experiments on synthetic and real-world networks. The reduced model successfully reproduces the principal transitions and steady-state branches in all three dynamical systems considered. Agreement is strongest for relatively homogeneous networks and deteriorates when structural heterogeneity produces a broader distribution of node states. The closure diagnostics account for this loss of accuracy and indicate when a single effective state is no longer sufficient. The reduction therefore provides a tractable description of resilience in systems with coexisting pairwise and higher-order interactions.

We introduce a class of nonautonomous nonlinear oscillator equations of mixed Liénard type that arises from a parametric deformation of the commutative factorization procedure applied to second-order ordinary differential equations with periodic solutions. Their solutions, in particular the isochronous waveforms, are obtained in closed form through a Riccati reduction scheme for the power-law choice of the factorization function, $\phi(x)=kx^{q}$, $k\in \mathbb{R}$, and $q\in\mathbb{N}$, corresponding to the so-called modified Emden oscillators, and do not depend on the arbitrary deformation parameter $A_1$. The variational structure of the equation is characterised by a Lagrangian supplemented with a generalised Rayleigh dissipation function that contains a non-standard cubic term in $\dot{x}$. We also show that multiplying the equation of motion by the Jacobi multiplier $M(x)=x^{-2A_1}$ a position-dependent-mass (PDM) form is obtained with related friction and restoring force.

Self-organized spatial patterns are central to nonequilibrium physics and cell biology, yet locating instabilities in multi-component, reaction-diffusion networks remains challenging because standard eigenvalue analyses scale with the number of biochemical states and rely on reaction kinetics often poorly constrained by experiments. Exploiting the common mass-conserving structure of protein reaction kinetics and the fact that nonlinear feedback is typically confined to membrane reactions while lateral membrane diffusion is negligible, we develop a geometric classification that predicts stationary, and approximately also oscillatory, pattern-forming instabilities from reactive equilibria. The developed criteria reduce the stability analysis from the full component space to the space of conserved species. On this reduced space, slope matrices—describing the change of equilibrium cytosolic densities with respect to total species densities—govern onset. Based on densities in chemical equilibrium, this approach eliminates the requirement of comprehensive kinetic knowledge frequently lacking in experimental systems. Thus, a broad range of instabilities can be understood as mass-redistribution instabilities—self-amplifying mass redistribution caused by shifting local equilibria. We apply these criteria to models for the Escherichia coli Min system and the Caenorhabditis elegans polarity system and show that the reduction extends to the mixed-dimensional dynamics in systems coupling bulk cytosolic dynamics with membrane dynamics on the boundary. Together, these results provide an interpretable, broadly applicable, and experimentally accessible framework for diagnosing and designing pattern formation in multicomponent nonequilibrium systems on the basis of conservation laws.

Maxwell's wave equation in the presence of a cubic nonlinearity and a periodically time-varying refractive index (a photonic time crystal) is reduced, for spatially monochromatic waves, to a nonlinear Mathieu equation. Near the principal momentum gap this equation admits an autonomous two-dimensional reduction whose complete Hamiltonian phase portrait can be obtained analytically. We derive the two homoclinic separatrices corresponding to temporally localised momentum gap solitons, identify the nonlinear centres and the critical Hamiltonian value $H_c$, and calculate the point of maximum linear parametric gain. We then consider spatially localised pulses and show how the nucleation of multiple spatiotemporal gap solitons can produce a broad supercontinuum in momentum space; for stronger seeds, transient extreme nonlinear localisation can accompany an abrupt additional broadening of this momentum spectrum. These results establish a direct connection between Floquet amplification, nonlinear saturation, homoclinic dynamics, and momentum space spectral broadening in nonlinear photonic time crystals.

We identify a general mechanism for the finite-time thermalization of macroscopic observables, such as coarse-grained charge densities, in terms of elementary forms of the quantum dynamics of initial states: (1) aperiodicity, which provides a computable measure of (2) a dynamical partially ergodic exploration of the Hilbert space. Specifically, this mechanism predicts the equilibration of all (concentrated) macroscopic observables, in almost all states in an initial ensemble and almost all times within finite and longer intervals, given only the observable-independent information that the return probability of the ensemble of initial states is small over a finite time range. As a special case, it also accesses standard results on equilibration over infinitely long times in terms of (stronger versions of) the effective dimension of initial state delocalization in the energy eigenbasis. Our results incorporate macroscopic thermalization into the domain of operational quantum statistical mechanics, recently developed to provide finitely computable criteria for microscopic thermalization. We discuss an overall characterization of this approach as establishing connections between (1) the decay of a (theoretically or experimentally) computable probe indicating memorylessness, (2) a fundamental invariant mechanism in terms of the alignment of observables or states in the Hilbert space, and (3) predicting different natural forms of (classical and) quantum thermalization, most of which rigorously recover conventional eigenstate-based descriptions of infinite-time thermalization as a special case but provide stronger accessible predictions over finite observation times in the thermodynamic limit.

Many traditional models of electoral dynamics emphasize opinion change, social influence and persuasion mechanisms. In contrast, contemporary polarized societies often exhibit relatively stable partisan blocs coexisting with an electorally relevant mobile population. We propose a minimal dynamical model in which electoral competition is governed not by direct persuasion between opposing blocs, but by the dynamics and allocation of a mobile electoral interface. The electorate is partitioned into two stable partisan blocs and a mobile fraction, while direct transitions between the polarized blocs are strongly suppressed. Mobile voters are allocated between the competing blocs through a Fermi-like probabilistic rule governed by rejection asymmetries. Analytical calculations show that the electoral susceptibility to political shocks is proportional to the stationary size of the mobile electoral interface, identifying electoral mobility as the key quantity controlling the macroscopic response of polarized systems to external perturbations. The model naturally predicts a continuum of mobility regimes, ranging from frozen polarization to highly responsive electoral states characterized by a broad mobile interface. These results suggest that, in strongly polarized elections, aggregate electoral changes may be governed primarily by fluctuations at the mobile electoral interface rather than by large-scale ideological conversion.

Using patch-matrix-product methods, we study two reversible three-state interaction-round-a-face cellular automata introduced by Klobas and Prosen [J. Phys. A 55, 094003 (2022)] - the species-preserving and species-flipping rules - and a coherent quantum deformation interpolating between them. Within an explicit translationally invariant ansatz, the two classical rules share a one-parameter family of quasi-local conservation laws with an auxiliary-dimension-three realization. Every regular member is also the first centered logarithmic derivative of an analytic auxiliary-dimension-two family passing through the identity observable. On the closed span of this family and the three elementary local charges, the only nonzero Euler velocities are $\pm\sqrt{3/23}$. Consequently, $\sqrt{3/23}$ is a rigorous lower bound on the maximal Euler speed in any larger conserved sector. For the species-flipping rule, this value agrees with the reported extrapolated value $0.361$ of Klobas and Prosen, strongly supporting completeness of the known sound-active charges. The species-preserving rule admits, in addition, a staggered generating family whose logarithmic derivatives form an infinite tower of strictly local charges. In the quantum deformation, the same algebraic structures yield exact low-bond-dimension invariant states (scar candidates), exponentially long-lived quasi-local quasimodes, and diagonal quasi-local charges that produce a nonzero Mazur bound after projection away from the three elementary local charges.

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.

The study of magnetic chaos has traditionally focused on macroscopic variables under external driving. Here we demonstrate a new type of magnetic chaos, termed moiré magnetic chaos, associated with mesoscopic magnetic domain variables in twisted bilayer CrI3 without external driving. The domains are stabilized by a characteristic interlayer exchange frustration, which supplies the multiple dynamical degrees of freedom required for autonomous chaos. Through micromagnetic simulations, we show that relaxation toward moiré magnetic textures is extremely sensitive to minute local perturbations of the initial state, characterized by substantial finite-time Lyapunov exponents and a final-state sensitivity that persists over five decades of perturbation amplitude. Statistical analysis further reveals that the resulting domain configurations are stochastic and pairwise uncorrelated. Our results identify a form of microscopic, undriven chaos in twisted magnets that extends nonlinear magnetism beyond the conventional driven regime.

Nonequilibrium membrane pattern formation is studied using meshless membrane simulation. Bound proteins are considered to have two states that generate different membrane spontaneous curvatures. Protein binding and unbinding occur cyclically owing to chemical potential differences, as an off-lattice active Potts model. It is found that this cyclic binding/unbinding can induce blinking domains, with oscillating size: convex domains of the proteins with a higher spontaneous curvature grow, and subsequently, the proteins change to the other state with a lower spontaneous curvature, resulting in domain shrinkage. These processes repeat. In thermal equilibrium, hexagonal convex domains are formed by the competition between bending and surface tension energies, so that they are stably formed only under positive surface tension. However, blinking domains can form even in tensionless membranes.

In this paper, we travel through an unknown territory of dromions, to unearth and show the new properties/attributes of dromions which have never been brought to the fore since they were first discovered by Boiti etal [1]. These new attributes are brought out by investigating a generalized (2+1) dimensional Nonlinear Schrodinger (NLS) equation by exploiting Truncated Painleve approach. The signatures that can be attributed to dromions include existence of firewall and reflection at the boundary, uneven distribution of energy among different bound states, amplitude dependence on adjacent dromions, etc. We have corroborated the main analytical results with numerical simulations also. We categorically state that these properties are universal and can be extracted in any (2+1) dimensional nonlinear partial differential equation. We do believe that these properties which shed more light on the behaviour of dromions may have wider repercussions in nonlinear optics, Bose-Einstein condensates and plasma physics.

In this article, the nonlinear Helmholtz equation with non-Kerr nonlinearity, such as self steepening and self frequency shift, is considered. A traveling wave transformation is applied, and the extended nonlinear Helmholtz equation is reduced to a Hamiltonian dynamical system. Then, the reduced Hamiltonian system is analyzed by classification of equilibrium points, phase space analysis, and the construction of exact wave solutions. The relationship between the reduced dynamical coefficients and the original physical parameters is further established through a parameter space analysis. It is shown that self steepening directly modifies the reduced dynamics, whereas self frequency shift acts through the compatibility condition for the real traveling wave reduction. Together, these non-Kerr effects reshape the phase space geometry and traveling wave structure. Localized and periodic traveling waves are obtained, with their existence determined by the balance among dispersion, nonparaxiality, Kerr nonlinearity, and non-Kerr effects. Furthermore, a periodically forced version of the reduced system is examined to study the transition from regular to irregular dynamics. It has been observed that external forcing can induce complex oscillatory behavior. Bifurcation analysis, time series evolution, phase space analysis, largest Lyapunov exponent, and Poincaré section demonstrate the emergence of quasiperiodic and chaotic responses under sufficiently strong forcing. All analytical branches are verified through full-equation residual evaluation, while a few selected branches are additionally examined through direct numerical propagation and robustness tests under complex Gaussian perturbations. The results show that self steepening directly renormalizes the effective nonlinear dynamics, whereas self frequency shift restricts the admissible real-envelope traveling wave manifold.

Low-dimensional reductions provide a useful framework for studying high-dimensional dynamics on complex networks, but most existing approaches are restricted to pairwise interactions. Here, we develop a one-dimensional reduction for dynamical systems on networks with purely higher-order interactions. The reduction is formulated through an effective higher-order interaction strength ($\beta_{\Delta}$), associated with the triangular interactions of the underlying network and the dynamical system's effective state. We present a theoretical framework for the dimension-reduction approach and validate it across three dynamical models with exclusively higher-order interactions. We find that the reduction accuracy is mainly determined by the homogeneity of node states, i.e., the deviations in state values become very small. Numerical results on synthetic and real networks show that the reduced model captures the effective steady states and transitions of the full system with good accuracy.

Abrupt transitions in complex systems are often preceded by early warning signals. However, most indicators rely on the notion of critical slowing down and do not generally extend to rate-induced tipping where transitions can occur without local loss of stability. This is problematic in stochastic, nonautonomous systems where internal variability and time-varying variables interact to shape tipping onset. We use Koopman operator theory to develop a unified early warning framework for both bifurcation and rate-induced tipping in stochastic systems. Our approach builds on residual Koopman mode decomposition that measures discrepancies between dynamics and their finite-dimensional approximation, and extends it to the control setting by augmenting the observable space with time-varying control variables. In idealized examples, the resulting indicators recover expected signatures near bifurcation points and improve detection in rate-induced regimes where classical indicators fail. We further show that learned embeddings through deep learning outperform prescribed dictionaries, especially in a high-dimensional setting. Applied to simulations of the Atlantic Meridional Overturning Circulation, our Koopman-based indicators distinguish tipping from non-tipping trajectories and reveal interpretable spectral signatures prior to critical transition.

In our previous paper, we described the proposal and validation of a new physical implementation-type simulator (spontaneous symmetry breaking machine (SSBM)) for a combinatorial optimization problem that utilizes a phenomenon dual to the model in which a pseudo-spin many-body system is created by spontaneous symmetry breaking. Furthermore, we performed numerical simulations with different initial fluctuations on a large-scale benchmark problem (K_{2000}) and reported that we identified a condition under which pseudo-spin patterns gradually converge during the solution search process to ultimately form a single pattern and that the cut value of this converged solution reached 99.7% of the known best. In this paper, we discuss the conditions under which the SSBM can exhibit behavior analogous to that of a continuous-variable-type Ising machine, and report that by adding an effect similar to physical annealing to these conditions, the SSBM can be made to search for the best known cuts in K_{2000}.

Microbial populations generally evolve in fluctuating environments under time-varying conditions. These are often described by binary switching models, sometimes seen as coarse-grained feast-famine cycles, in which resource availability switches abruptly between abundant and scarce conditions. However, experimental studies suggest that feast-famine environments actually exhibit more complex temporal dynamics. Here, we study how two strains, one growing slightly slower than the other, compete for the same resources in fluctuating environments comprising a finite number of intermediate states, each having its own carrying capacity. Environmental switching between these states and their carrying capacities represents gradual changes in nutrient availability. This class of multi-state stochastic switching models can be interpreted as a coarse-grained description of feast–famine cycles and allows us to investigate strain competition under the gradual recovery and depletion of resources. By computational and analytical means, we characterise the population dynamics in these multi-state fluctuating environments. In particular, we study how the switching rates and distribution of carrying capacities affect the population-size statistics, fixation probability, and mean fixation time. By comparing these results with their counterparts in binary environments, we clarify how the frequency and amplitude of environmental fluctuations influence population dynamics in coarse-grained feast-famine cycles.

In this paper, we investigate the integrability of Lotka-Volterra (replicator) systems arising from interaction matrices generated from corresponding graph structures, continuing work started by Visomirski and Griffin [J. Phys. A., 58:015701, 2025] and Evripidou et al. [J. Phys. A., 55:325201, 2022] (among others). In particular, we define a new family of graphs, the spoked graphs, and show that all dynamics generated from this family are integrable. In reference to Evripidou et al. (2022), we define a new anti-cloning operator and show that its action on balanced tournament graphs (with odd vertex count) generates new graphs whose dynamics are integrable. Interestingly, we provide numerical evidence that this anti-cloning operation leads to chaotic behaviour when applied to other graph families (e.g., the directed cycles that generate the classically integrable Volterra lattice). This work completes a taxonomy of all integrable dynamics generated by directed graphs with up to six vertices started by Visomirski and Griffin (2025), and suggests several future directions of study on this topic.

Memory is a key determinant when cultural items compete for attention and, consequently, for success, as in the case of songs on a music chart. For modeling, one adds memory to Lotka-Volterra models, the reference for Markovian competitive processes. Here we treat memory in terms of an exponential moving average, finding that it leads to an extended region of winnerless chaos characterized by log-normal popularity statistics in trailing top-k charts. Importantly, the observed log-normal behavior collapses to a power-law distribution when the feedback dynamics is fast on the scale of the charting period. This result is in agreement with the observed statistics of real-world music charts (e.g., Billboard and Spotify). In a chaotic state, the size of the largest Lyapunov exponent is a measure of how unpredictable the system is. We find that the largest Lyapunov exponent varies strongly as a function of parameters in the phase where winnerless chaos is stable. Interestingly, the sets of parameters obtained by comparing simulations with real-world cultural-item dynamics extracted from Google Books and Google Trends, movies, Reddit, Wikipedia, Twitter, and scientific publications, are located close to the points in parameter space where the largest Lyapunov exponent reaches its local maximum, namely, close to the point of maximal unpredictability. This result suggests that cultural competitive processes are maximally chaotic when memory is a key determinant.

We present quasistationary ring-shaped soliton necklaces in a two-component envelope propagating in a medium with competing cubic-quintic nonlinearity. Metastable propagation of soliton necklaces results from a balance of repulsion between adjacent out-of-phase solitons in one component and confinement by the boundary of a flattop region in the other. Numerical simulations demonstrate metastable propagation over about a hundred diffraction lengths, even with random noise added to the input envelopes. The maximum number of solitons in metastable necklaces can be controlled either by changing the size of individual solitons or by adjusting the width of the flattop region hosting the necklace.

We introduce Semantic Lenia, an artificial life framework that transforms Large Language Model (LLM) inference from a static optimization problem into a continuous, closed-loop dynamical system. By establishing a non-linear homeostatic feedback loop to dynamically balance semantic attraction and syntactic repulsion, we demonstrate the emergence of “Homeostatic Solitons”-metastable semantic structures that actively resist repetitive crystallization. Our exhaustive parameter sweeps map a critical “Habitable Ridge” where applied steering forces balance the model's intrinsic syntactic inertia. This approach successfully maintains generative trajectories in a numerically sensitive critical regime, triggering profound abductive leaps without structural collapse, and reveals a capacity-dependent scaling trend in the syntactic inertia across different model sizes.

We study data-driven prediction of coarse-grained dynamics in multiscale PDE systems. Adopting a closure-free operator-learning viewpoint, we apply a linear coarse-graining map and learn a surrogate evolution operator for the resolved field directly from filtered high-fidelity trajectories. Motivated by the Mori-Zwanzig formalism, we propose a spatiotemporal neural operator mapping a resolved history slab on $\Omega\times[-T_{\mathrm{in}},0]$ to a resolved future slab on $\Omega\times[0,T_{\mathrm{out}}]$. Spatial mixing uses Fourier convolution, while temporal mixing uses a causal kernel operator with position-attention weights on time lags. This causal temporal operator encodes finite-memory effects in the resolved dynamics while preserving the directionality of the history-to-future map. To improve rollout robustness and suppress nonconservative artifacts, we embed a flux-form inductive bias by parameterizing the windowed update in explicit divergence form. We also provide a data-driven guideline for selecting the memory length $T_{\mathrm{in}}$ via the decorrelation time of a closure-injection diagnostic computed from filtered trajectories. We validate on the coarse-grained viscous Burgers' equation, the Kuramoto-Sivashinsky equation, and two-dimensional turbulent flows, obtaining stable autoregressive rollouts with improved long-horizon accuracy and statistical fidelity.

We investigate the collision dynamics of two bright sphalerons in a (1+1)-dimensional deformed $\phi^6$ scalar field theory with a symmetric potential possessing false vacua. Two one-parameter realizations of the model, referred to as the barrier and well models, are considered and their static and linear instability properties are first reviewed. We then study head-on collisions of boosted sphalerons over a broad range of initial velocities and deformation parameters. The scattering dynamics exhibit a rich variety of final states, including the production of kink-antikink pairs, long-lived oscillons in true and false vacuum, multiple oscillons propagating in false-vacuum regions, and radiative decay. A particularly remarkable outcome is the emergence of a long-lived bubble of the false broken vacuum bounded by a kink-antikink pair, which repeatedly collapses and re-expands before eventually decaying into an oscillon. These results demonstrate that deformed $\phi^6$ theories with false vacua exhibit considerably richer sphaleron dynamics than previously known and provide new insight into the role of unstable localized configurations in nonlinear field theories.

In this paper we present Open Multi-Processing (OpenMP) Fortran 90/95 programs to solve the Gross-Pitaevskii equation for a rotating dipolar Bose-Einstein condensate (BEC) in two and three dimensions, which is a new version of our previous published programs for a dipolar Bose-Einstein condensate without rotation. After the recent experimental study of a rotating dipolar BEC [L. Klaus et al., Nature Phys. 18, 1453 (2022)], the present programs will be useful tools for related theoretical investigation. The algorithm used is the split-step semi-implicit Crank-Nicolson scheme for imaginary- and real-time propagation to obtain stationary states and BEC dynamics, respectively, as in the previous version [L. E. Young-S. et al., Comput. Phys. Commun. 286 (2023) 108669].

In this work, we derive selection rules imposed by Liouville integrability conditions for monomial charge densities with arbitrary powers. For a certain monomial Hamiltonian system, the selection rules reduce to a negative Pell equation, and its solutions generate an infinite set of integrals of motion that are mutually in involution. Furthermore, we study the correspondence between combinatorial polynomial sequences and Liouville integrable Hamiltonian field theories in 1+1 dimensions. We show that the Motzkin system coincides with the dispersionless limit of the Levi system, while the binomial system is equivalent to the dispersionless derivative nonlinear Schrödinger equation. Additionally, we show that the binomial Hamiltonian model admits a reduction to the inviscid Burgers equation and its higher-order charges generate generalized Burgers-type conservation laws.

In recent years, integrable hierarchies have been used to great advantage for the enumeration of combinatorial maps. They have led to recurrence formulas with respect to the size and genus of the maps, e.g. for triangulations, bipartite quadrangulations and bipartite maps, and for constellations. These formulas are not only remarkably simple but also provide the fastest way of calculating these numbers of maps. With the exception of Louf's work on constellations, it has however remained a challenge to obtain recurrence formulas that control the degrees of the faces of the maps. Here we show how to achieve this for bipartite maps with bounded face degrees. By combining equations from the KP hierarchy and from the Virasoro constraints, a differentially algebraic system is obtained. It couples the generating functions of bipartite maps with bounded root face degrees while controlling the numbers of edges, black vertices, white vertices and number of faces of each degree (and in particular the genus). Finally, this system of ODEs is shown to give recurrence formulas that allows to calculate all the corresponding numbers of maps.