We consider the soliton solutions of a recently proposed coupled Sasa-Satsuma-mKdV equation using the Kadomtsev-Petviashvili reduction method. The system consists of a complex-valued component coupled with a real-valued one. Under zero or nonzero boundary conditions, we derive four distinct classes of soliton solutions: bright-bright, dark-dark, bright-dark, and dark-bright. These solutions are derived from the vector Hirota equation, for which the bright, dark, and bright-dark soliton solutions are provided in the Appendix. We perform asymptotic analysis of soliton collisions for each class of solutions, in which inelastic collisions are observed between bright-bright solitons. In the dark-dark case, we identify soliton profiles similar to the Sasa-Satsuma equation, including double-hole, Mexican hat, and anti-Mexican hat solutions; this study further explores the collisions between these structures and hyperbolic tangent shaped kink solitons. Regarding the bright-dark case, beyond the expected soliton-kink interactions, we report and analyze a notable collision occurring between kink solitons.

General equations describing shear displacements in incompressible hyperelastic materials, holding for an arbitrary form of strain energy density function, are presented and applied to the description of nonlinear Love-type waves propagating on an interface between materials with different mechanical properties. The model is valid for a broad class of hyper-viscoelastic materials. For a cubic Yeoh model, shear wave equations contain cubic and quintic differential polynomial terms, including viscoelasticity contributions in terms of dispersion terms that include mixed derivatives $u_{xxt}$ of the material displacement. Full (2+1)-dimensional numerical simulations of waves propagating in the bulk of a two-layered solid are undertaken and analyzed with respect to the source position and mechanical properties of the layers. Interfacial nonlinear Love waves and free upper surface shear waves are tracked; it is demonstrated that in the fully nonlinear case, the variable wave speed of interface and surface waves generally satisfies the linear Love wave existence condition $c_1 < \abs{v} < c_2$, while tending to the larger material wave speed $c_1$ or $c_2$ for large times.

In this paper, we present a unified approach to constructing continuous and discrete $\mathrm{PGL}(3)$-invariant integrable systems, formulated in terms of the common dependent variables $z_1,z_2$, from linear spectral problems and their factorisation. Starting from third-order spectral problems, we first provide explicit forms of the differential and difference invariants, generalising the Schwarzian derivative and cross-ratio to the rank-$3$ setting. The factorisation induces dualities among linear spectral problems, underlying the exact discretisation and multi-dimensional consistency of the associated Boussinesq systems. Then, we derive both continuous and discrete $\mathrm{PGL}(3)$-invariant Boussinesq systems, representing natural rank-$3$ generalisations of the Schwarzian KdV and cross-ratio equations. A geometric lifting-decoupling mechanism is developed to explain the reduction of these systems to the $\mathrm{PGL}(2)$-invariant Boussinesq equations. Finally, we derive a ${\mathrm{PGL}}(3)$-invariant system of generating PDEs together with its Lagrangian structure, in which the lattice parameters serve as independent variables, providing the generating PDE system for the Boussinesq hierarchy.

We investigate the long-time asymptotic behavior of a class of solutions to the defocusing Manakov system under nonzero boundary conditions. These solutions are characterized by a $3 \times 3$ matrix Riemann Hilbert problem. We find that they exhibit interesting asymptotic behavior within a narrow transition zone in the $x$-$t$ plane. We determine the leading-order asymptotic term and the error bound in this region, and we demonstrate that the leading term can be expressed in terms of the Hastings-McLeod solution of the Painlevé II equation. The proof is rigorously established by applying the Deift-Zhou nonlinear steepest descent method to the associated Riemann Hilbert problem.

In complex systems, events occur at irregular intervals that inherently encode the underlying dynamics of the system. Analyzing the temporal clustering of these events reveals critical insights into the non-random patterns and the temporal evolution. Existing techniques can effectively quantify the overall clustering tendency of events using global statistical measures. However, these macroscopic approaches leave a critical gap, as they do not attempt to investigate the dynamics of individual clusters. Analyzing individual clusters is essential, as it helps comprehend the local interactions that actively drive the system dynamics, which may be obscured by global averaging, while simultaneously revealing the time scales involved. To address these limitations, we propose a complex network-based framework for analyzing clustering of events occurring at irregular intervals. The framework establishes connections using arrival times, transforming the time series into a network. Network properties are then used to quantify the clustering. Further, a community detection algorithm is used to identify individual clusters in time series. We illustrate the method by applying it to standard arrival processes, such as the Poisson process and the Markov-modulated Poisson process. To further demonstrate its scope, we apply the method to two diverse systems: the time series of droplet arrivals in turbulent flows and the R-R intervals in electrocardiogram (ECG) signals.

Gottschalk's surjunctivity conjecture states that for all group universes and finite alphabets, every equivariant and continuous selfmap of the full shift, known as cellular automaton, cannot be a strict embedding. Not all surjective cellular automata are injective. However, if the surjectivity condition is replaced by a certain strengthened property called post-surjectivity then all post-surjective cellular automata must be bijective whenever the universe is a sofic group. A group universe is said to be post-injunctive if every post-surjective cellular automaton with finite alphabet over this group universe must be bijective. Gromov's injectivity lemma states each injective cellular automaton over a subshift can be extended to an injective cellular automaton over every subshift which is close enough to the initial subshift. In this paper, we obtain analogous results where injectivity is replaced by other fundamental dynamical properties namely post-surjectivity and pre-injectivity. We also study various stable properties of the class of post-injunctive groups in parallel to properties of surjunctive groups. Among the results, we show that semidirect extensions of post-injunctive groups with residually finite kernels must be post-injunctive.

The focus of this paper is the study of the moduli space of representations of fundamental groupoids of surfaces $\Sigma$ with boundaries with values in $G:=GL_n(\mathbb C)$. In absence of marked points on the boundary, this moduli space is realized in many equivalent ways: as the moduli space of linear local systems on $\Sigma$, as the moduli space of representations of the fundamental groupoid $\Pi_1 (\Sigma)$, as the space of monodromy data and as character variety. By adding marked points to the boundary of $\Sigma$ in order to capture irregular singularities, the Betti moduli space has been generalized in several ways by many authors. Although it is clear that these different approaches describe essentially the same spaces of mathematical objects, exactly how they fit together has not yet been established. Motivated by the broader programme of establishing an explicit and conceptually coherent relationship between the existing approaches to the study of the decorated Betti moduli space, in this paper, we develop a categorical framework that allows for a systematic definition of the \dfn{decorated Betti moduli spaces} space, in the presence of higher order poles, designed to specialize to the different points of view encountered in the literature.

Turbulent flows exhibit robust universal features -- including logarithmic mean velocity profiles, scale-invariant energy spectra, anisotropy constraints and strongly non-local transport -- yet a unifying dynamical principle underlying these phenomena remains elusive. We show here that turbulence can be organized around an emergent oscillatory degree of freedom governing the Reynolds stress. Starting from the exact non-local representation of the stress in terms of a propagator, we demonstrate that the spectral structure of the response contains a dominant complex-conjugate pair of poles, implying an effective oscillator coupled to the mean flow. In wall-bounded turbulence, the near-wall Airy structure selects and stabilizes this mode through non-local feedback, yielding the logarithmic velocity profile and fixing the asymptotic von Kármán constant, $\kappa \simeq 0.39$. In homogeneous turbulence, the same dynamical picture closes the inertial-range energy balance and yields the Kolmogorov constant as $C_k=2/[3(1-2^{-2/3})]\simeq 1.80$ at leading order. The resulting formulation leads to a closed tensorial set of mean-field equations in three spatial dimensions, significantly cheaper than direct numerical simulation yet rich enough to support geometry-dependent reduced dynamics interpretable as distributed networks of interacting oscillators. The associated phase field admits a geometric description connected with Berry phase, anisotropy evolution on the Lumley triangle, and an effective gauge-covariant structure of phase transport. These results suggest that turbulence is governed not by an algebraic closure, but by a dynamical and geometric organization of the mean stress.

Topological Bloch oscillations are a hallmark of quantum transport phenomenon in which wavepackets undergo oscillatory motion driven by the interplay between an external force and topological edge states and serve as a powerful dynamical probe for the geometric properties of topological bands. Spin-orbit coupling (SOC) has also emerged as a crucial ingredient for manipulating quantum states in materials, with the corresponding gauge fields arising from the Rashba and Dresselhaus interactions. In this work, we investigate the propagation of spinor wavepackets in a honeycomb Zeeman lattice governed by the Gross-Pitaevskii equation. By tuning the relative strengths of Rashba and Dresselhaus SOC, we engineer a non-Abelian gauge field that drives anomalous topological Bloch oscillations (ATBOs). Unlike conventional topological Bloch oscillation (TBOs), these ATBOs exhibit asymmetric motion, including a freezing effect in one half of the oscillation cycle, which can be tuned by the SOC parameters and external forces. Our findings establish SOC-based non-Abelian gauge fields as a powerful mechanism controlling topological quantum dynamics, with implications for spintronic devices and quantum data processing.

Contrary to accepted turbulence folklore, which holds that no mathematical relation exists between the Navier-Stokes equations (NSEs) and the multifractal model (MFM) of Parisi and Frisch, we develop a theory that reconciles the MFM with Leray's weak solutions of Navier-Stokes analysis. From a combination of Euler invariant scaling and the NSEs we also derive the Paladin-Vulpiani inverse scale $L\eta_{h,pav}^{-1} = Re^{1/(1+h)}$ which acts as a mediator between the two theories. This is achieved by considering $L^{2m}$-norms of the velocity gradient to find a correspondence between $m$ and the local scaling exponent $h$ in the multifractal model. The parameter $m$ acts as if it were the sliding focus control on a telescope which allows us to zoom in and out on different structures. The range $1 \leqslant m \leqslant \infty$ is equivalent to $-2/3 \leqslant h_{min} \leqslant 1/3$, which lies precisely in the region where Bandak et al. (2022, 2024) have suggested that thermal noise makes the NSEs inadequate and generates spontaneous stochasticity. The implications of this are discussed.

We investigate how the defining statistical features of three-dimensional turbulence respond to systematic reductions of the Fourier-space triadic interaction network. Using direct numerical simulations of both fractally and homogeneously decimated Navier-Stokes dynamics, we show that progressive thinning of the set of active modes leads to a systematic suppression of intermittency and, most strikingly, to the vanishing of the mean dissipation rate in the large-Reynolds-number limit. Structure-function exponents collapse onto their dimensional values, the multifractal singularity spectrum contracts, and the analyticity width extracted from the exponential spectral tail increases monotonically with decimation-each indicating a substantial regularization of the velocity field. Together, these results provide direct evidence that anomalous dissipation in incompressible turbulence is not a generic property of the Navier-Stokes equations, but instead requires the full combinatorial richness of their triadic nonlinear interactions.

Most real-world networks exhibit a significant degree of modularity. Understanding the effects of such topology on dynamical processes is pivotal for advances in social and natural sciences. In this work we consider the dynamics of Kuramoto oscillators on modular networks and propose a simple coarse-graining procedure where modules are replaced by effective single oscillators. The method is inspired by EEG measurements, where very large groups of neurons under each electrode are interpreted as single nodes in a correlation network. We expose the interplay between intra-module and inter-module coupling strengths in keeping the coarse-graining process meaningful. We show that, when modules are well synchronized, the phase transition from asynchronous to synchronous motion in networks with 2 and 3 modules is very well described by their respective reduced systems, regardless of the network structure connecting the modules. Applications of the method to real networks with small modularity coefficients reveals that the approximation is also very accurate if oscillators in each module are identical. The method reproduces global synchronization patterns despite the low synchronizability of some modules, possibly allowing for the inference of the mean synchrony of each module when individual dynamics are not known.