Papers

By constructing Lax operators in the loop algebra of the $\mathbb{Z}_2^2$-graded extension of the Lie superalgebra $\mathfrak{osp}(1|2)$, we derive a $\mathbb{Z}_2^2$-graded extension of the Camassa-Holm equation. The resulting equation is an integrable nonlinear PDE for a system of four $\mathbb{Z}_2^2$-graded commutative functions, each associated with a distinct $\mathbb{Z}_2^2$-degree. We further show that the $\mathbb{Z}_2^2$-Camassa-Holm equation admits a bi-Hamiltonian structure. As a consequence, it possesses infinitely many conserved quantities, including one with non-trivial $\mathbb{Z}_2^2$-degree, which are mutually in involution with respect to the $\mathbb{Z}_2^2$-graded Poisson brackets.

The controlled dynamical generation of radial dam break flows (DBFs) is demonstrated in two-dimensional droplet environments. These ultracold mixtures feature competing mean-field and quantum fluctuation effects and are initialized on a disk-shaped uniform density. Depending on the value of the latter, reflecting also the dominance of attractive or repulsive interactions, three dynamical response regimes are identified. At small densities a moving wavetrain of bright ring solitons is observed (dubbed ring DBF), at intermediate densities a ring kink-dispersive shock wave, while at large densities the flows encompass composite ring kink-rarefaction waves moving towards the center. The emergence of these patterns, whose formation would be prohibited in two-dimensional Kerr media due to wave-collapse, is rooted in the competition between attraction and repulsion inherent to droplet environments, modeled by the extended Gross-Pitaevskii equation. The characterization of the ensuing nonlinear waveforms is further corroborated by reduced models based on the cubic-quintic nonlinear Schrödinger equation. The different wave patterns emanating from dam break flows can be experimentally realized using current state-of-the-art ultracold atom experiments.

The most energetic coherent structure of a hovering full-scale quadcopter is associated with a self-organizing process in which the individual rotor vortices synchronize their frequencies while undergoing merging events, yielding a globally correlated structure. We identify and characterize this phenomenon by applying spectral modal and conditional analyses to assimilated three-dimensional velocity data acquired via Shake-The-Box Lagrangian particle tracking. The dataset captures a high-Re, turbulent flow further complicated by time-varying rotor speeds stemming from active flight control, low-frequency vehicle drift, finite spatio-temporal resolution, and measurement uncertainty. Most coherent structures recover established single-rotor features such as tip vortices and their subharmonic pairing. The globally synchronized vortex merging manifests as a spectral peak at an incommensurate frequency below the rotor band, which cannot be explained by single-rotor aerodynamics, subharmonic instabilities, or band-to-band triadic interactions. Instead, conditional averaging provides evidence of the aforementioned intermittent, distinctly non-subharmonic vortex-merging process involving all four rotor wakes. Establishing whether or not this phenomenon is observed across different flight conditions and configurations remains speculative; however, the consistent characterization of the globally synchronized vortex merging using complementary frequency- and time-domain analyses despite experimental complexities, in particular rotor speed variations, demonstrates its robustness.

Buckling of thin cylindrical shells under axial compression, a classical example of a subcritical instability, is highly sensitive to small imperfections, with minute geometric variations causing large changes in buckling threshold. To uncover the origin of this sensitivity, we use numerical continuation and bifurcation analysis while systematically varying the depth and size of a single localized Gaussian defect. We show that the instabilities of the imperfect shell originate from localized equilibria already present in the perfect shell. By breaking translation symmetry, the defect pins these equilibria and changes how they connect to the imperfect base state. Small changes in defect geometry can thereby switch the bifurcation that triggers buckling, producing non-monotonic and discontinuous changes in buckling threshold and abrupt changes in buckling mode. Imperfection sensitivity is therefore not simply sensitivity to imperfection magnitude, but sensitivity of the underlying bifurcation structure to imperfection geometry.

It is often assumed that general stationary solutions to the wave kinetic equation are described solely by conserved fluxes and temperature. The families of solutions parameterized by these quantities are referred to as warm cascades, and are derivable from differential models of wave kinetics. These models rely on the assumption of local interactions, and we demonstrate by a time scale argument that the nonlocal effect of thermal damping on the cascade may result in dependence on an additional parameter, the UV cutoff, calling into question the validity of warm cascade solutions. This is shown for systems with both three-wave and four-wave interactions. We then consider the depletion of the cascade by this damping, and present a numerical demonstration of this nonlocal and nonuniversal behavior.

What can and cannot be predicted by a model when only limited information is available? We answer this question by constructing a ladder of information limits that organizes prediction tasks according to their information requirements: from memoryless and memory-augmented trajectory forecasts to event prediction, stationary statistics, and generative laws. The analysis, which is independent of the model structure or architecture, identifies the minimum attainable error along with the information limitations that give rise to it. The approach distinguishes information hidden in unresolved variables from information recovered through time-delayed observations and connects these contributions to Mori–Zwanzig memory. The ladder also reveals why Lyapunov growth alone cannot characterize reduced-order prediction error, why stationary statistics may remain predictable even when individual trajectories become unpredictable, and why stochasticity can represent (but cannot recover) missing information. Numerical studies of the Kuramoto–Sivashinsky equation and the Lorenz system illustrate these results across the ladder. This unified perspective provides a principled language for understanding the fundamental limitations of ROMs for a given prediction task and for clarifying whether improved predictions require richer observations, greater precision in the input, or additional memory.

Learning in recurrent neural networks can fundamentally reshape their underlying dynamics, transforming initially chaotic activity into stable task-dependent behavior. We develop a non-equilibrium dynamical mean-field theory(DMFT) to describe this transition during learning. We show that a slow feedback-driven learning process generates an evolving effective feedback strength that drives the network through a transition from chaotic to stable dynamics defined by a bifurcation of the DMFT solution. By deriving the two-time correlation function throughout learning, we identify a critical feedback strength and a corresponding learning rate dependent critical time separating these regimes. The transition arises from the progressive deformation of an effective dynamical landscape by the growing learned feedback structure. Starting from the untrained state, the theory predicts the time evolution of the network output during training and shows quantitative agreement with numerical simulations.

Many examples of integrable sigma models admit families of integrable deformations that involve coupling the undeformed “seed” theory to auxiliary fields with algebraic equations of motion. In this work, we apply this paradigm to the case where the seed theory is the lambda model, which interpolates between the Wess-Zumino-Witten model and the non-Abelian T-dual of the principal chiral model. We present an auxiliary field deformation of the lambda model and exhibit a Lax connection whose flatness is equivalent to the remaining equations of motion after imposing the algebraic field equations. We also carry out a Hamiltonian analysis of the deformed theory, verifying that it obeys the assumptions of theorems in arXiv:2605.18213 which guarantee both (i) that the model enjoys the Maillet structure that ensures involution of its conserved charges, and (ii) that the theory possesses a classical Yangian. Thus, as in other cases, the auxiliary field deformation of the lambda model changes its classical dynamics while leaving the integrable structure essentially unchanged.

Crises and population responses can form coupled dynamical systems, with crisis conditions shaping protective behaviours and collective responses altering the crisis trajectory. Existing models rarely capture both evolving crisis conditions and the reinforcement-dependent spread of competing behaviours. We propose a complex contagion with feedback (CCF) model that couples competing complex contagion with crisis dynamics through bidirectional feedback. Agents stochastically switch between competing states according to social reinforcement and adoption complexity, which is influenced by crisis conditions and information campaigns, while population-level behavioural change feeds back into the crisis. We formulate and analyse the model for a well-mixed population, characterising its equilibria and associated stability conditions. As a case study for model validation, we integrate CCF model with a census-calibrated agent-based model of COVID-19 transmission in Australia to represent social distancing adoption and discontinuation. We compare the CCF model with an existing opinion dynamics model and find that, despite using fewer parameters, it achieves comparable performance in reproducing recurrent infection waves. The resulting framework is parsimonious, analytically tractable, and can be adapted to different types of crises.

Phase-field trajectories governed by the Allen–Cahn equation determine the effective force $G=F'/\eps^2$, but not the free-energy potential $F$ and interfacial scale $\eps$ separately. We introduce surface-tension-calibrated Allen–Cahn identification (STAC), a two-step method that resolves this ambiguity using an independently supplied surface tension. First, space–time weak moments produce a linear inverse problem without differentiating the measured field. A Bernstein representation with positive coefficients enforces two stable pure phases and a single central energy maximum. Second, the surface tension fixes the remaining scale algebraically. We prove scale equivalence of the forward model, uniqueness of the calibrated pair within this class, and finite perturbation bounds that remain valid when the potential vanishes at the pure phases. Experiments with resolved synthetic phase fields separate the effects of the structural constraint, ridge regularization, and constitutive approximation error. Across 80 shared noisy records, unconstrained degree-five regression produces six potentials with an incorrect well structure. Positivity combined with ridge regularization eliminates these violations and reduces the mean force error from 8.70% to 5.66%. Lower-degree models are nevertheless more accurate in some tests, and an admissible law outside the positive Bernstein family shows how the constraint can introduce approximation bias. STAC thus recovers a structurally admissible, physically scaled potential while making explicit which information comes from the observed dynamics and which comes from the independent surface-tension datum.

We develop a geometric framework in which classical thermodynamics is reformulated as a Hamiltonian dynamical system. An explicit canonical mapping $(q,p,t,H)\leftrightarrow(V,-P,S,-T)$ identifies Maxwell's relations as the characteristic equations of the thermodynamic Poincaré–Cartan one-form, in direct parallel with Hamilton's equations of motion. Treating thermodynamic potentials as action functionals, a variational principle recovers both the Maxwell relations and the thermodynamic constraints (adiabaticity, isothermality) as conserved first integrals. For adiabatic processes in an ideal gas, this yields explicit second-order ordinary differential equations for $V(T)$ and $P(T)$, each governed by a temperature-dependent Lagrangian. Assembling the component Lagrangians into a multi-parameter Lagrangian 1-form $\mathcal{L}$, we prove that the closure condition $d\mathcal{L}=0$ holds on the solution manifold, establishing that path-independence of thermodynamic state functions is a geometric consequence of multi-time integrability rather than an independent postulate.

We obtain exact solutions of the nonlinear Dirac equation in $1+1$ dimensions of the form $\psi(x,t) = e^{-i\omega t} \psi(x)$ for scalar-scalar (SS) plus vector-vector (VV) interaction with interaction Lagrangian given by $$L_{I} = \frac{g^2}{\kappa+1}[(\bar{\psi} \psi)^{\kappa+1} +\frac{1}{p} (\bar{\psi} \gamma_\mu \psi \bar{\psi} \gamma^{\mu} \psi)^{\kappa+1}]$$ where $p > 0$ but arbitrary otherwise. We look for solutions with $0 < \omega < m$ where $\omega, m$ are frequency and mass, respectively. We find solutions for all values of $\omega$ in this range. We compute the charge $Q$ and the energy $E$ for each of the solitary wave solutions and explore the region in the ($p, \kappa$) parameter space in which solitary wave bound states exist (i.e., for which $E/Q < m$). We show that for all the cases while both $E$ and $Q$ depend on the coupling constant $g$, their ratio $E/Q$ is independent of $g$. We further find that in case $p\kappa \le 1$, the charge density for all the solitary waves have only single hump while for $p \kappa > 1$ there is a transition from double to single hump and we determine it as a function of $\omega/m$. We notice that for all $p$ there is a transition at $\kappa=2$ in the behavior of $E/Q$ as a function of $\omega$ which we speculate is related to the onset of instability of the solutions at $\kappa=2$. We obtain the nonrelativistic reduction of the two-parameter family to a non-relativistic modified nonlinear Schrödinger equation (NLSE) and discuss the stability of the single hump olitary waves in the domain of validity of the modified NLSE.

Chaos emerges when stretching is repeatedly recycled by reinjection. We uncover a new route to chaos in which the decisive variable is the temporal order of non-normal tangent maps: periodic and chaotic states can share essentially the same one-step stretching statistics while their ordered products acquire opposite Lyapunov growth. We introduce the ordered-product growth rate $h_L$ over $L$ successive tangent maps, which reveals how states indistinguishable at one step separate under geometric composition and identifies the finite composition scale at which chaos emerges. We use this mechanism to establish a new form of global chaos control: minute phase actions reorient the successive non-normal amplification directions so that their geometric composition becomes contracting, suppressing chaos at fixed dissipation without reducing local amplification or targeting a preselected orbit.

Whether human or large language model (LLM), an agent in a discussion reads only a few of the others' contributions, bounded by cognition, context, or cost. LLM collectives can settle on a wrong consensus even when a majority starts out correct; we ask how far that reading bound alone decides the outcome. We model the bound with one number, the message capacity, which sets how many of the others' messages an agent reads, and generate the communication network from it. Over 31,824 randomized queries, we found that an 8-billion-parameter model's judgment of a claim effectively reduces to a logistic function of a weighted sum of its inbox, the update rule of a stochastic binary neuron with divisively normalized weights. From these weights and the network's degree statistics alone, the wrong consensus should become unreachable from any start once agents read, on average, fewer than 6.4 of their 31 sources. In 1,414 episodes with assigned starts the prediction failed: the correct side won in fewer than 50% of episodes from every start, and in only 28-45% when 75% of agents started correct. The failure traces to the field, the threshold that a claim's wording sets for the agent's answer before any message is read: the experimental claims' fields lay below the calibration mean, and with each claim's own field the same weights reproduce the outcomes. Reversing the wording showed that the threshold follows what a claim asserts, not whether it is true. On a second 8B model the pipeline predicts claim-dependent bistability; transition points appeared where computed, and an eight-claim calibration matched in 15 of 16 conditions. At 70B the assertion bias is not detected. Thus a collective's fate is largely set by two single-agent measurements: the threshold a claim's wording sets, and the message capacity that sets the transition point.

Labor strikes do not necessarily lead to complete institutional disruption, particularly when non-striking sectors maintain part of normal activities. We propose a minimal dynamical model to investigate how such residual activity affects the persistence of labor conflicts. The model couples strike participation, institutional response, and overall functionality through two competing mechanisms: a direct suppressive effect of residual activity on strike participation and an indirect effect arising from the reduced institutional pressure to resolve the conflict. Analytical calculations reveal a counterintuitive regime in which increasing residual activity increases, rather than decreases, stationary strike participation. We characterize this behavior through a strike-response susceptibility and derive an exact compensation threshold separating the conventional regime from a resilience-paradox regime. At the threshold, the two mechanisms exactly balance, making stationary strike participation independent of residual activity. Numerical results confirm the analytical predictions and illustrate the underlying mechanism. Our results show that institutional resilience can have an unintended dynamical consequence: by absorbing part of the disruption generated by a strike, continued operation may reduce the urgency of institutional response and thereby favor conflict persistence. This provides a simple example of how individually stabilizing processes can collectively generate counterintuitive outcomes in complex social systems.

We study the Arnold-Beltrami-Childress (ABC) flow in perturbative regimes near its three elementary integrable coordinate axes. Starting from a degenerate family of periodic streamlines of the limiting integrable system, we use first-order averaging to prove the existence of two isolated periodic solutions of the perturbed flow. The three perturbative regimes are related by the cyclic symmetry of the ABC vector field. In each case, the corresponding two-dimensional averaged system possesses two simple equilibria, one elliptic and the other hyperbolic, depending on the sign of the relevant parameter ratio. Poincaré sections illustrate how the degenerate periodic structure of the integrable limit breaks under perturbation and how the isolated periodic orbits predicted by averaging emerge within the surrounding dynamics. These periodic orbits provide the natural link between the perturbative analysis and the integrability problem. The eigenvalues of the linearized averaged system determine the leading-order behavior of the nontrivial characteristic multipliers of the bifurcating periodic solutions. Combining this relation with the Poincaré-Llibre-Valls criterion, we prove that, in a neighborhood of each of these periodic orbits, there exists no nonconstant first integral $H\in C^1$ that is regular along the orbit. The analytical results are further illustrated and confirmed by direct shooting, monodromy-matrix computations, and continuation of the elliptic and hyperbolic branches, demonstrating that the first-order averaging approximation remains quantitatively accurate over a substantially wider parameter range than the strict asymptotic regime.

We address three-dimensional (3D) solitons maintained by spin-orbit coupling (SOC) in the binary self-interacting Bose-Einstein condensate (BEC) held in the 3D optical lattice (OL). The analysis reveals that the SOC-OL interplay results in the formation of stable full-vortex (FV) soli-tons, built as sets of four density peaks residing in neighboring wells of the lattice potential, with the superimposed global vortical phase, and site-centered semi-vortices (SV), in which the vorticity is present in only one component. The full-vortex solitons, with their specific phase textures, do not exist in a uniform BEC with SOC. Full-vortex and semi-vortex states in the binary self-attractive BEC are stable despite the possibility of the supercritical collapse in the 3D system. Such states also exist, as gap solitons, in the self-repulsive 3D system. In terms of the chemical potential and number of particles, the stability regions of the 3D full-vortex solitons and semi-vortices expand with the increase of the SOC strength and OL depth. The results open the route to the creation of 3D complexes of vorticity-carrying condensates that can be realized with existing experimental techniques.

The hulusi is a wind instrument that was invented in Yunnan Province, China, and has become tremendously popular in recent years. It consists of a mouthpiece, a gourd, and three bamboo tubes, all with free reeds made of copper. The main bamboo tube in the middle has seven finger holes. In this instrument, the pipe length, not the free reed's eigenfrequency, determines the instrument's pitch, unlike, for example, with the Western accordion or the blues harp. In this study, a machine learning model implemented in the COMSAR framework (https://github.com/ifsm) was used to investigate the timbre characteristics of the hulusi to cluster different instruments and pitches. The measured hulusi pitches C, B, A, G, and F were analyzed according to seven psychoacoustic features, among which only the spectral centroid, sharpness, and fractal correlation dimension are shown to form pitch clusters. These timbre features were used to train Kohonen self-organizing maps (SOMs) for clustering. Brightness and sharpness analysis revealed that the highest pitches were less bright and less sharp than mid- and low-range pitches were. Furthermore, the fractal correlation dimension, which mainly determines the chaoticity of the initial transients, was the best-clustering timbre feature for the hulusi, with the highest pitches showing the least chaoticity. This result is supported by defining a cluster quality index for the SOMs.

We revisit the Nakamura conjecture, which relates the Tomimatsu-Sato solutions of stationary axisymmetric gravity to finite Toda molecules. While the conjecture has been established partially, its general rotating sector remains an open problem. We show that the Toda determinants underlying the conjecture possess a natural weight grading. In particular, the two functions entering the Ernst potential have weights n^2 and n^2-1, and this grading extends systematically to shifted determinants labelled by partitions. In coordinates adapted to the Toda generators, each differentiation corresponds to adding one box to the associated Young diagram and increases the weight by one. The same integer n^2 also appears in the zero-order term of the Nakamura bilinear operator, revealing a compatibility between the differential equation and the determinant grading. The partition structure further explains the previously unresolved behavior of second derivatives. Repeated differentiation in one direction produces an internal sector and an external sector requiring only a one-step extension of the Wronskian hierarchy; the latter is reduced by a local three-term Pluecker relation. Thus weight grading, Young-diagram growth, Wronskian enlargement, and Pluecker reduction emerge as parts of a single determinant structure. The unit weight relation n^2 = (n^2-1) + 1 also singles out the elementary Toda seed as a natural third object, suggesting a possible route toward a genuine trilinear formulation. Although no trilinear closure is assumed here, the present construction reduces the remaining general-n Nakamura problem to definite determinant-minor identities and provides a structural framework in which such a formulation can be investigated.

In nonlinear prediction, two state representations with the same local uncertainty volume can have radically different predictive value. We consider differentiable finite-time dynamics together with a quadratic terminal requirement that specifies which terminal state differences matter. Pulling this requirement back through the tangent map defines a requirement-induced predictive geometry on the present tangent space. This geometry is fixed by the dynamics and the requirement; under the fixed local uncertainty-volume constraint used here, the capacity parameter fixes only the scale of the representation that realizes it. We show that the unique optimal covariance has the inverse predictive shape and, equivalently, that optimality is characterized by isotropization of propagated uncertainty in the terminal requirement geometry. The gain over a baseline-isotropic representation is the arithmetic-to-geometric mean ratio of the predictive spectrum, making it a coordinate-invariant measure of predictive anisotropy. In Euclidean coordinates this yields a general finite-time spectral law, and for two-dimensional area-preserving dynamics an exact identity between predictive gain and the condition number of the optimal geometry. Stronger directional stretching therefore increases potential predictive gain only together with the anisotropy required to realize it; bounding that anisotropy imposes a finite gain ceiling. Standard-map and double-pendulum calculations demonstrate these relations in nonlinear Hamiltonian systems. The framework provides a basis for adaptive sensing, state estimation, and model-predictive computation in which representation geometry is matched to predicted task consequences.

Recent experiments demonstrating superconducting-like characteristics in materials above their equilibrium transition temperatures following strong sub-picosecond optical excitation have stimulated theoretical interest in phase-disordered superconductors subject to a periodically modulated superfluid stiffness. Previous analyses of the single-layer optical conductivity have shown that such a drive can enhance the low-frequency $1/\omega$ response, even when the period-averaged stiffness remains unchanged. We extend these studies to bilayer systems, such as the YBa$_2$Cu$_3$O$_{6+x}$ (YBCO) cuprate, by considering drives with different symmetries. Specifically, we investigate a layer-symmetric drive, in which the two planes stiffen and soften simultaneously, and a layer-antisymmetric drive, in which one plane stiffens while the other softens. Both drives enhance correlations of the layer-averaged phase and generate comparable superconducting-like optical conductivities. Only the symmetric drive, however, strengthens relative-phase locking and the associated counterflow response to a magnetic field parallel to the planes. Thus, an enhanced $1/\omega$ optical conductivity does not by itself imply enhanced magnetic screening. In view of applications to driven YBCO, where superconducting-like in-plane optics and transient diamagnetism have both been reported, our result provides a symmetry diagnostic: within a phase-synchronization mechanism, the experimentally active low-frequency stiffness modulation must contain a substantial layer-symmetric component.

Repeated decisions exhibit temporal structures such as persistence, direction-dependent switching, recurrent alternation, and abrupt transitions. We examine the generative sufficiency of a two-dimensional fast-slow dynamical system. The system combines a cubic fast equation with linear slow feedback and is analyzed through its equilibrium geometry, trace-determinant structure, equilibrium-fold loci, candidate Hopf boundaries, and singular critical manifold. An explicit observation operator projects continuous trajectories to a scalar signal and applies a binary readout, separating latent state-space dynamics from observable behavior. The analysis establishes a unique-equilibrium regime and, for suitable parameters, a three-equilibrium wedge with a central saddle. The outer equilibria are attracting only where their traces are negative. The analysis also identifies the simple-zero condition required for ordinary saddle-nodes, trace-zero positive-determinant spectral boundaries compatible with oscillatory instability, and the attracting and repelling branches of the critical manifold. Prescribed nonautonomous sweeps numerically illustrate direction-dependent switching (T1), transient episodic recurrent switching (T2), and an abrupt localized regime shift (T3). The attracting-equilibrium regime associated with prolonged state retention (T4) is characterized analytically. The resulting correspondence is intended as a test of generative sufficiency at the level of observable temporal organization, rather than as an identification or empirical validation of latent cognitive mechanisms.

We present a direct and nonrecursive construction that maps an arbitrary planar polynomial dissipative flow $\dot{\alpha}=h(\alpha,\alpha^*)$ to an open quantum system in Gorini–Kossakowski–Sudarshan–Lindblad form. Given the homogeneous components of the classical drift, the corresponding Hamiltonian and collapse-operator blocks are obtained algebraically and independently at each degree. No recursive cancellation of lower-order terms is required: in the large-amplitude limit $|\alpha|\sim S\to\infty$, ordering-generated lower-degree terms are parametrically suppressed, while mean-field closure yields the prescribed leading $O(S)$ drift for semiclassically localized states. This provides a simple and systematic route from dissipative classical dynamics to explicit open-quantum-system realizations, without claiming a unique microscopic quantization of the classical flow. We demonstrate the construction for stable fixed points, a Hopf bifurcation, and a bistable-ring flow. The corresponding Liouvillian spectra recover the classical relaxation exponents, the radial Floquet exponent and neutral phase direction of a limit cycle, and the separation between local relaxation and inter-ring switching in the bistable case. These examples show that the construction can be used not only to analyze a given quantum model, but also to design open quantum systems with prescribed classical dynamical structures in the semiclassical limit.

Ultrashort laser pulses can induce order out of chaos. A femtosecond pulse melts a metal surface for about a hundred picoseconds; under photoexcitation with no net polarization axis, the resolidifying film nonetheless organizes into regular nanoscale patterns. Here we derive from first principles the nonlinear hydrodynamics responsible for it and show it accounts for the full phase diagram of morphologies from measurable laser and material parameters. Starting from the Navier-Stokes equations in the lubrication limit, the molten film obeys a damped Kuramoto-Sivashinsky equation that reduces to a Swift-Hohenberg formalism near threshold, with every coefficient fixed by the laser fluence, melt thickness, optical skin depth, and thermocapillary response. The ratio of melt depth to optical skin depth sets the pattern wavelength, while the inter-pulse delay, via post-pulse damping, gates the route from flat surface to ordered pattern to chaos. We validate the phase diagram (hexagonal cavity arrays, bistable mixed states, labyrinths) numerically and confront it with pulse-resolved experiments on two metals, Ni and Fe3Cr. The agreement between theory and experiment quantitatively links the microscopic laser-material response to the macroscopic pattern selected after many pulses. This order is self-limiting: the same corrugation the instability builds scatters the modes that build it, and simulations, with every scale fixed by material parameters, place the coherence breakdown near a roughness of 30 nm after a few tens of pulses, a finite window confirmed by pulse-by-pulse measurements of orientational order on the experimental surfaces; at stronger coupling, disorder does not destroy the pattern but saturates and pins it instead. This yields a closed, predictive description of laser-driven self-organization, from single-pulse melt dynamics to the arrest of nanoscale order.

We study the quantised version of the Kadomtsev-Petviashvili equation in the formulation proposed in [4]. We analyse two aspects. One is the conservation of the first higher charge $H_3$ (besides the total mass $H_0$, the momentum $H_1$ and the Hamiltonian $H_2$). We find the expression of $H_3$ and prove commutativity, provided we regularise a specific ultraviolet divergence. The result then does not depend on the choice of a regulator, under fairly general assumptions. We then study the Bethe ansatz in finite volume and obtain a general solution to the two-body problem. We have used AI supervised by us to aid the proofs.

In this work, we construct analytical kink-antikink $(K\bar{K})$ configurations in the non-integrable $\varphi^4$ model by mapping exact solutions of the integrable sine-Gordon system via a field deformation. This procedure yields two distinct classes of configurations, parametrized by the initial half-separation and velocity. We compare these profiles with the standard additive superposition Ansatz in terms of vacuum structure and equation-of-motion residuals, deriving explicit closed-form expressions for the corresponding integrated squared residuals. For small separations, the mapped soliton-antisoliton configurations exhibit appreciably smaller residuals than the naive superposition Ansatz, which in turn performs better than the mapped two-soliton configurations. Although the mapped fields are not exact solutions of the $\varphi^4$ equation of motion, they provide mathematically consistent, topologically sound, and physically motivated initial data for numerical studies of kink-antikink scattering and resonance phenomena.

Connection formulas for the asymptotic solutions of a system of n> 1 coupled Painleve'-II equations with symmetry-breaking parameters are written explicitly. An asymptotically exact WKB approach to these formulas relies on the quantum-mechanical independent crossing approximation for an explicitly time-dependent Schroedinger equation.

Hamiltonian and Lagrangian mechanics are powerful frameworks for analyzing physical systems. Previous work has extended these formalisms to ecological systems, such as the predator-prey Lotka-Volterra (LV) system. In this Article, we derive a canonical Lagrangian for the two-dimensional LV model directly from its Hamiltonian representation. We find that the two-dimensional LV system admits a standard canonical Lagrangian formulation with one degree of freedom and a non-quadratic kinetic structure. This formulation admits a mechanical interpretation in terms of a particle moving in a potential well, where the non-standard kinetic structure produces a position-dependent damping term that can instead act as "revving." The derivation provides a direct connection between predator-prey dynamics and a canonical formulation of mechanical dynamics. As a verification of the construction, we apply Noether's procedure to the explicitly time-independent derived Lagrangian and reveal that the well-known Hamiltonian of the LV system is the corresponding conserved quantity. We also uncover a subtle redundancy associated with the choice of canonical momentum and its identification with the original population variables.

We investigate local conserved quantities in non-Hermitian bosonic lattice systems whose underlying Yang–Baxter structure remains unclear. Focusing on a one-dimensional bosonic chain with quartic interactions of creation operators, we provide a new representation of the local conserved quantities previously constructed by Sanatani and Shiraishi. Using Fourier transformation and trigonometric identities, we systematically derive these conserved quantities and show directly that they form a mutually commuting family. We further extend the construction to interactions extending beyond a single site, cubic interactions, asymmetric hopping, and an on-site potential. Our results provide a unified framework for constructing and characterizing local conserved quantities in this class of non-Hermitian many-body systems.

We establish rogue wave patterns associated with the fourth Painlevé equation $(\mathrm{P}_{\mathrm{IV}})$ in the multi-component nonlinear Schrödinger and Hirota equations. The generalized Hermite and generalized Okamoto polynomials arise in representations of rational solutions of $\mathrm{P}_{\mathrm{IV}}$, and we show that their roots determine two classes of rogue wave patterns when one of the internal parameters of rogue wave solutions is large. Specifically, the generalized Hermite polynomials arise from rogue wave solutions represented by Schur-polynomial determinants with consecutive indices, whereas the generalized Okamoto polynomials arise from analogous determinants with index jumps of three. Numerical examples for both equations agree with the predictions.