Papers

We investigate Hamiltonian curl-force systems with indefinite kinetic energy. We first reconsider Berry's polynomial Hamiltonian curl-force model, whose numerically observed closed trajectories motivated an integrability conjecture. A Painlevé analysis yields a non-principal resonance spectrum, so that the corresponding Laurent series cannot accommodate the required number of arbitrary constants of the general solution. The model therefore fails the standard Painlevé test. We then introduce a four-parameter curl-force family and identify the parameter locus on which this system is integrable. We construct a second Hamiltonian, compatible Poisson tensors, separated complex characteristic variables, and a Lax representation. More generally, the separated form yields polynomial integrable curl-force Hamiltonians of arbitrary degree. We also show that the same construction admits a higher time-derivative potentialisation whose free limit is the degenerate Pais-Uhlenbeck oscillator. Finally, we analyse zero-curl invariant reductions and elliptic periodic solutions, and exhibit an isolated periodic orbit outside the integrable regime. This illustrates that closed trajectories alone do not imply Liouville or Painlevé integrability.

Quantum cryptographic primitives beyond key distribution remain a less well understood area of research. In classical cryptography, one-way functions underpin nearly all standard cryptographic protocols, motivating the search for meaningful quantum analogues and for a clear understanding of the physical and computational mechanisms that could enforce one-wayness. In this article, we review quantum one-way functions and a range of closely related quantum-state primitives, including one-way state generators, pseudorandom quantum states, and efficiently indistinguishable pairs of states. We discuss both computational and information-theoretic notions of quantum one-wayness, emphasizing the different adversarial models and security assumptions that underlie these constructions. We compare and contrast the various proposed primitives, and clarify their conceptual relationships. Particular emphasis is placed on questions of physical realizability, experimental feasibility, and robustness to noise. Finally, we outline open problems and future directions toward the development of practical quantum cryptographic primitives beyond key distribution, and the emergence of a broader quantum-cryptographic ecosystem.

We solve cyclic finite-dimensional quantum histories for arbitrary time-dependent unitary steps, without assuming that one step has finite order. The propagation Hamiltonian is a unitary connection Laplacian on a cycle; its complete gauge invariant is the monodromy $M=U_{L-1}\cdots U_0$. Its spectrum is $\lambda_{a,k}=1-\cos((2\pi k-\theta_a)/L)$, where $e^{i\theta_a}\in\mathrm{spec}(M)$. Thus the exact history sector is isomorphic to $\mathrm{Fix}(M)$, while frustration, the gap above a nonempty zero-energy sector, the determinant, and the finite-temperature trace are obtained in closed form. Ordinary spectral data recover the multiset of monodromy phase cosines but not phase orientation; low energy certifies proximity to an exact relational history. We then define the predictive quotient of a sharp finite clock relative to an accessible operator system as the unique coarsest event alphabet preserving all conditional statistics on a history sector. A finite-error theorem shows that threshold clustering recovers this quotient when the minimum diamond separation of inequivalent event channels exceeds four times the estimation error, and proves an optimal record-count bound. With full matrix access and homogeneous step $U$, the minimal number of clock events is the projective order of $U$. We distinguish the normalizer of the clock algebra from transformations preserving the coherent history code and classify oriented exact sharp clock changes by $U(r)\times\mathbb{Z}_L$ on a rank-$r$ history sector; without orientation the cyclic factor becomes dihedral. Reversible changes of full-information clock fibers are necessarily unitary, so irreversible coarse-graining is not exact clock covariance. Minimal realizations of a complete history Gram kernel are uniquely unitarily equivalent, with a finite-data Procrustes bound. Independent finite-matrix code verifies the main results.

In this paper, we develop a unified Hamiltonian framework for dynamical systems with multiple independent dissipation channels based on q-contact geometry. As an application, we embed an elastoplastic damage model with two internal variables into the non-uniform 2-contact framework, explicitly construct the underlying geometric structure, derive the effective equations of motion, and validate the theoretical predictions through numerical simulations, which confirm the exponential dissipation law and the associated rescaled conserved quantity.

In this paper, we study the time-periodic problem for the Vlasov-Poisson-Boltzmann (VPB) system with a given time-periodic external force in the whole space $\mathbb{R}^3$. The force is allowed to be non-potential. Around the global Maxwellian, we prove the global existence of small solutions in a hybrid function space that combines the low-frequency Besov framework for the forced Boltzmann equation with a corresponding control of the self-consistent electric field. The main novelty lies in the treatment of the nonlinear Vlasov force $-\nabla_x\phi \cdot \nabla_vf + \frac{1}{2}(v \cdot \nabla_x\phi)f$ at low frequencies. Rather than treating it as a generic source term, we exploit the Poisson equation and macroscopic balance laws to recover the structural cancellation required for the VPB semi-group estimates, which combined with high-frequency energy estimates and weighted microscopic propagation, yields a closed global well-posedness theory. We further prove the asymptotic stability of small solutions driven by the same force. When the external force is time-periodic, Serrin's method yields a unique time-periodic solution with the same period, together with its stability. As a direct consequence, our result also gives the existence and stability of stationary solutions when the external force is time-independent.

The higher-order Craig-Bampton (HCB) method improves reduced-model accuracy through successive residual-mode enrichments. This study proposes an a priori estimator for the eigenvalue error remaining in the first-order HCB model (HCB-1) without using the full-order eigensolution. The estimator is defined as the relative difference between the HCB-1 and HCB-2 eigenvalues. Under symmetric Ritz approximation assumptions, ordered positive eigenvalues, nested augmented spaces, and preservation of retained eigenvalues by the System Equivalent Reduction Expansion Process (SEREP), the estimator is shown to be non-negative and to provide a lower bound on the relative HCB-1 eigenvalue error. A projection-based expression and an exact error identity are also derived, clarifying the conditions under which the estimator is accurate and showing that the lower-bound property alone does not guarantee a prescribed tolerance. Numerical results for a plate, an elbow pipe, and a reactor pressure vessel show estimation-to-error ratios of about 0.85-0.99, while a cost assessment reports the additional expense associated with the HCB-2 reference solution.

We introduce a notion of completeness for two-dimensional second-gradient elastic continua and propose a microstructural route toward its synthesis. A continuum is said to be complete if the Hessian of the stored energy with respect to the second-gradient variable is locally positive definite, so that every nonzero admissible increment of the placement second gradient is quadratically controlled in the highest-order part of the energy about each configuration. Starting from the Principle of Virtual Work, we derive the constitutive relations, equilibrium equations, and admissible boundary interactions for a broad class of fibrous second-gradient continua whose stored energies depend on fiber stretch, stretch gradients, and curvature. The theory is then applied to several continua motivated by pantographic microstructures. Classical pantographic sheets are shown to be incomplete, while bi-pantographic fabrics enlarge the class of components of the second gradient detected by the energy but remain incomplete. Finally, we formulate a tri-pantographic continuum associated with a proposed three-family architecture and prove that it is complete. The examples illustrate how microstructural architecture influences both the completeness properties of an effective continuum and the pointwise form of its higher-order boundary interactions.

Two hypergeometric families of mixed-type multiple orthogonal forms are constructed for rank-one $q\times p$ matrices of weights, with arbitrary $p$ and $q$: Jacobi and Laguerre I systems. Both normalized mixed forms are obtained explicitly for admissible near-diagonal multi-indices. The power-vector components are terminating generalized hypergeometric polynomials, while the hypergeometric-vector components are finite sums of such polynomials. In the Jacobi case, these sums are expressed as finite combinations of terminating Kampé de Fériet polynomials evaluated at $(x,1)$. For the mixed beta–Euler Laguerre I system, a finite triangular system relates the residues at finite poles to the terms generated by the Euler operator. Gamma-quotient Mellin formulas, Meijer $G$-representations, and Rodrigues formulas are derived for the complete mixed forms. On the step-line, the two biorthogonal systems satisfy dual recurrences governed by matrices with $p$ subdiagonals and $q$ superdiagonals. All recurrence coefficients are given by finite Gamma–Pochhammer expressions. Under normality and nonvanishing-pivot assumptions, Christoffel transformations and Gauss–Borel factorization yield bidiagonal factorizations of the recurrence matrices. The lower factors have closed Pochhammer formulas, while the upper factors are expressed through finite Christoffel tau-determinants or, equivalently, cross-ratios of shifted moment minors. Both Christoffel chains close explicitly in the mixed Piñeiro specialization.

A magnetic billiard describes a charged particle constrained to a planar domain: the particle moves along circular Larmor arcs in the interior and undergoes specular reflection at the boundary. The round disk has an explicit first integral that is polynomial in the velocity, and a central rigidity question asks whether any other smooth convex table can have such an integral. We answer this question negatively in the strong-field regime. Let $\Omega\subset\mathbb R^2$ be a bounded strictly convex domain with smooth boundary $\gamma$, let $r=|B|^{-1}$ be the Larmor radius, and assume $0<r<r_0(\gamma)/2$, where $r_0(\gamma)$ is the maximal embedded tubular radius. If the magnetic billiard admits a nonconstant first integral polynomial in the velocity variables, of any finite degree, then $\Omega$ is a disk. This removes the finite exceptional set of strong field strengths left by the earlier polynomial nonintegrability theory. The rigidity mechanism has two logically independent stages. First, the highest reflection mode gives a boundary winding identity. This determines the degree of the top coefficient and places all of its roots strictly inside the table, but it does not show that those roots coincide. Second, after the two leading reflection identities are continued to the normalization of the complexified boundary, their valuations at infinity exclude simultaneous poles of the coordinate functions. The remaining one-sided poles force the entire top coefficient to be one linear factor of multiplicity equal to the Fourier degree. Combining the location theorem with this root-collapse theorem produces a constant-angle relation between the boundary tangent and a radial direction, hence circularity. No real-analytic boundary hypothesis is imposed: analyticity follows from the algebraic strong-field parallel curves.

Deep Operator Networks (DeepONets; arXiv:1910.03193) typically encode an input function through point values on a fixed discretization. Building on the Topological DeepONet framework of Ismailov (arXiv:2603.11972), we replace point samples by continuous linear functionals drawn from the continuous dual of a Hausdorff locally convex space $({V},\{p_\alpha\}_{\alpha\in A})$, whose topology is generated by a point-separating family of seminorms rather than a single norm, and develop fixed and adaptive functional measurement systems. Measurements are combined with the coefficient-space Two-Step procedure of Lee and Shin (arXiv:2309.01020), while a training-only decoder and regularization stabilize the adaptive coordinates. We derive a discrete error decomposition separating measurement, output-basis, and neural-approximation errors, together with a Barron-rate refinement. The framework is evaluated on the antiderivative operator, a non-normable locally convex input space, heterogeneous Darcy flow, a controlled operator, and fixed-time and time-evolving Navier-Stokes vorticity operators. In the heterogeneous Darcy problem, the functional models retain nearly resolution-independent errors of 5.5-5.6% on unseen grids, while in the controlled problem adaptive measurements reduce the mean error below 1.2%. For the fixed-time Navier-Stokes problem, the Adaptive Topological DeepONet is the most accurate DeepONet-based model, attaining a mean relative $L^2$ error of 1.685% +/- 0.017% using 128 functional coordinates. A comparably sized Fourier neural operator (FNO; arXiv:2010.08895) achieves the lower error 0.832% +/- 0.172%, but requires the full 64x64 input field, twice the training time, and 10.7x greater peak GPU memory. The formulation provides compact, interpretable, and discretization-portable coordinates in the continuous dual $V'$, including for non-normable input spaces.

Starting from SU(N) on the lattice, we give a rigorous derivation of the Dyson–Schwinger equations in the continuum limit. We formulate the Dyson–Schwinger identities for the lattice Yang–Mills theory directly in terms of the link variables $U_\mu(m)\in SU(N)$, exploiting the invariance of the Haar measure under left group translations. This provides an exact lattice derivation of the corresponding master equation for the Wilson action, expressed through left-invariant Lie derivatives acting on individual links. Because the construction is carried out directly on the compact gauge group, it avoids the ambiguities associated with introducing Lie-algebra valued gauge potentials as primary integration variables at finite lattice spacing. For practical applications, in a second part we then break down the gauge degree of freedom by choosing Feynman gauge. We analyze the continuum-limit form of the resulting lattice identities and derive equations for the one- and two-point connected functions. Under a further simplifying reduction, these equations close to a tractable scalar system. Our results establish a direct bridge between exact lattice identities and the functional equations commonly used in continuum nonperturbative studies of Yang–Mills theory.

Accurate prediction of tritium behaviour in molten salt breeding blankets is essential for the design and safe operation of ARC-class fusion reactors. This work presents a fully open-source, component-scale multiphysics framework for modelling tritium transport in an ARC liquid immersion blanket. Neutron transport, thermal hydraulics, and hydrogen isotope transport are coupled using OpenMC, OpenFOAM, and FESTIM, leveraging dedicated tools enabling direct transfer of spatially resolved fields between solvers. Assuming a zero inlet concentration, steady-state simulations predict a total tritium inventory of approximately 243 mg, with the blanket reaching steady-state tritium throughput within approximately 30 min, which is of a similar order to previous system-level estimates. The results show that tritium transport is dominated by turbulence-enhanced diffusion, with strong localisation in flow stagnation regions and reduced accumulation in highly turbulent zones. Sensitivity analyses indicate that predicted inventories are governed primarily by the numerical stabilisation scheme, with only a modest dependence on the turbulent Schmidt number. The proposed workflow provides a transparent and extensible basis for high-fidelity analysis of tritium transport in ARC-class breeding blankets.

We are interested in modelling superfluid helium-4, the common isotope of helium, of atomic number 4. We present a mathematically solvable toy model of the phase transition between the normal liquid and superfluid phases and use it to show how an order parameter can be obtained for the superfluid.

We construct exact transparent boundary conditions (TBCs) for a time-continuous, spatially discrete Schrödinger equation that models a one-dimensional quantum lattice. Using a recently developed exact solution for the discrete system, we derive the Dirichlet-to-Neumann maps analytically via Laplace transforms. This yields a convolution-type boundary condition governed by Bessel functions. We rigorously demonstrate the consistency of this discrete formulation with its continuous counterpart in the continuum limit. Additionally, we present an efficient time-discretization scheme based on the trapezoidal rule for practical implementation. Numerical experiments using a Crank-Nicolson solver verify that our proposed TBCs eliminate spurious backscattering entirely and preserve reflectionless propagation of a Gaussian wave packet exiting the computational domain

Motivated by the study of the level lines of the $(2+1)$D Solid-On-Solid (SOS) model above a floor, near the corners of the box, we derive the limit law of an ensemble of $K$ curves from a $(1+1)$D SOS model, with an area tilt, and above a wedge-shaped floor in $\{-N,\ldots,N\}$. We show that there exist explicit critical points $\alpha_0=1>\alpha_1>\ldots>\alpha_K>0$ such that, for each $r \geq 1$, along the intervals $\pm(\alpha_{r} N,\alpha_{r-1} N)$, the bottom $K+1-r$ curves, rescaled by $(N^{2/3},N^{1/3})$, tend to the law of a Geometrically-Area-Tilted Ensemble of non-crossing Brownian paths (Brownian GATE), independently across those $2K$ intervals. All other curves, centered and rescaled by $(N,\sqrt{N})$, tend to a product of $K$ suitable Brownian bridges.

We consider a stochastic process on $\mathbb{Z}^d$ for $d \geq 1$. Given a translation invariant and ergodic starting configuration of finite clusters, each cluster $C$ performs a continuous time simple random walk with rate $|C|^{-\alpha}$. If it attempts to move to a vertex occupied by another cluster, it does not move, and instead the two clusters connect via a new edge. In all dimensions, we show that if $\alpha\ge 0$, there is almost surely no spontaneous creation of an infinite cluster within finite time. Moreover, for any $\alpha\le-1-2/d$ there is a finite-time blowup almost surely. In the regime $\alpha\in(-1,0)$ we show that the behavior greatly depends on the initial configuration. In addition, in dimension one, we establish the exact phase diagram.

Data assimilation (DA) uses Bayesian inference to update the state of a numerical forecast model with observed data. In this study, we propose a fundamentally different, unified approach to atmospheric data assimilation. We use latent video flow-matching to sample temporally consistent trajectories from a prior trained using ERA5 reanalysis (69 variables over an 8-day window). We also use posterior sampling to assimilate real observation sources, such as those from the NOAA Integrated Global Radiosonde Archive and the Integrated Surface Database. Because the prior generates a continuous trajectory, it naturally propagates information between observed and unobserved frames. Therefore, we can perform various DA tasks, such as filtering and smoothing, simply by changing the observed frames. Moreover, we generate full-state ensemble forecasts directly from sparse observations, achieving performance competitive with state-of-the-art observation-to-forecast models.

Random surfaces on a hard substrate often exhibit entropic repulsion, wherein the surface is propelled upwards to allow entropically preferable downward fluctuations. A particularly rich class of examples arises from the low-temperature 3D Ising model. A powerful approach to studying such surfaces is through their level curves, which form a family of non-intersecting random curves. In [CIW18, CIW19], an ensemble of Brownian lines with geometrically increasing area tilts was proposed as a putative limiting model in this case. This model falls outside the scope of techniques based on integrable or SDE structures, which have been key ingredients in the study of the Airy line ensemble. A particularly intriguing question about such line ensembles concerns their mixing properties when viewed as a Markov process, and in particular the rate of decay of correlations in time. For the Airy line ensemble, this decay is known to be inverse quadratic. The first quantitative bound on the decay of correlations in the area-tilted model, established in [CG25], was slower than polynomial in time. An earlier result [DLZ24] had established positivity of the spectral gap for the finite-line version of the ensemble, without quantitative bounds. This left open the important question of the true decay rate of correlations for the infinite ensemble. Settling this question for sufficiently large area-tilt strength, corresponding to sufficiently low temperature for the 3D Ising model, we prove exponential decay of correlations for the infinite ensemble and a uniform (in the number of lines) positive spectral gap for the finite ensemble. Our proof is based on establishing a precise form of separation of scales between curves of different indices, using a novel probabilistic approach involving embedding supercritical branching processes in the line ensemble.

The tube algebra, which carries the structure of a $C^*$ weak Hopf algebra, is a fundamental tool for characterizing topological excitations in topological phases. In this work, we generalize the tube algebra framework to codimension-2 defects in $(2+1)$D gapped phases described by Turaev–Viro–Barrett–Westbury TQFTs, with particular emphasis on Levin–Wen string-net models. We show that such codimension-2 defects are described by a defect tube algebra, which naturally carries the structure of a comodule algebra over the weak Hopf tube algebra associated with the topological excitations. Boundary and domain wall defects are then characterized by the representations of the corresponding boundary and domain wall defect tube algebras. In particular, a domain wall defect tube algebra can be regarded as a generalized Drinfeld double of the tube algebras associated with the two adjacent boundary defects. More generally, for a domain wall joining $N$ bulk phases, the associated domain wall defect tube algebra can be viewed as an $N$-tuple algebra. This perspective extends naturally to general $k$-defects, namely codimension-2 defects joining $k$ codimension-1 defects, for which the resulting generalized tube algebra carries the structure of a multicomodule algebra.

We associate height function models to a range of generalized two-dimensional XY models and prove that delocalisation of the height function model rules out exponential decay of the nematic order parameter in the primal spin model. The argument is based on a generalized loop representation which relates the variance of a height difference between two faces to the nematic order parameter. This generalizes the results by van Engelenburg and Lis (2023). The range of models also includes the much-studied generalized XY model as proposed by Korshunov and independently by Lee and Grinstein, where neighbouring spins interact ferromagnetically favouring parallel alignment and nematically favouring parallel or antiparallel alignment. Moreover, we recover the existence of a nematic region by a classical comparison argument of Pfister (1982).

We classify central extensions of the dg Lie algebra of derived global sections of the tangent sheaf on the punctured, formal d-disk for d > 2. (The d=2 case has appeared in previous work). We also prove a local, universal form of the Grothendieck-Riemann-Roch theorem for families of d-dimensional complex varieties.

Fast preparation of quantum error-correcting codes is essential for scalable quantum memories, but geometric locality and $U(1)$ charge conservation impose an unavoidable transport constraint. We combine exact complementary-channel geometry, charge-sector Haar analysis, and a gate-resolved connected-moment expansion to study one-dimensional covariant encoders under flagged erasure. Charge-Haar codes attain the universal adjacent-charge lower bound up to exponentially small corrections, yielding an exact $n^{-1/2}$ extensive-erasure law and a sharp half-erasure transition. For local number-conserving brickwork circuits, diffusion of the logical charge enforces an $\Omega(n^2)$ encoding-time lower bound; we also prove an $O(n^3)$ mixing bound for the classical component and reduce the remaining full-channel upper bound to a source-restricted low-support operator-spreading problem. These results identify diffusion as an operational limit on symmetry-constrained quantum coding and establish a route to its exact formation time.

We determine the sharp modulus of continuity, in trace distance, of the optimized Petz and sandwiched Rényi conditional entropies for every order $\alpha\in[\frac12,1)$. If two bipartite states are within trace distance $\delta$, then both conditional entropies differ by at most $\frac{1}{1-\alpha} \log[(1-\varepsilon)^{\alpha} +(D-1)^{1-\alpha}\varepsilon^{\alpha}]$, where $\varepsilon := \min\{\delta,1-1/D\}$ and $D$ is the effective dimension, given by the dimension of the first subsystem times the largest possible Schmidt rank. For every distance constraint $\delta\in[0,1]$, the bound is attained by an isotropic pair with a maximally entangled anchor. Taking $\alpha\uparrow1$ recovers the recent sharp continuity bound of quantum conditional entropy by Berta et al. [arXiv:2607.24687]. The proof linearizes the relevant concave Rényi functional at a comparison point dictated by the isotropic equality family. Schmidt-rank domination extends the equality geometry to an arbitrary anchor state, after which trace-distance duality and a noncommutative calibration estimate control the perturbation and anchor term without weakening the sharp constant. The latter estimate requires matrix analysis and is assisted by ChatGPT 5.6 Sol.

Extracting total correlations from a quantum system usually requires reconstructing its state, whereas many experiments access only a few measurement settings. A possible shortcut is to add the mutual informations obtained from complementary measurements; in dimensions above two, however, this procedure can count the same classical correlation twice. We establish that qubits are protected from such overcounting. For every two-qubit state, the correlations observed in two complementary local bases are bounded by the premeasurement quantum mutual information. The proof traces this protection to binary-entropy curvature on the Bloch ball and combines a qubit information-exclusion tradeoff with data processing under local dephasing. Consequently, two correlation tables give a tomography-free lower bound on total correlation. A score above one bit also certifies a quantitative one-way entanglement-distillation rate; when applied to the Choi state of a qubit channel, the same data lower bound its quantum capacity. The theorem therefore identifies both an operational use of complementarity and the trusted two-dimensional setting in which its correlation accounting is valid.

We study local minimizers of a sextic-potential Landau–de Gennes energy for nematic liquid crystals in the small-elastic-constant limit. Under the uniform energy and $L^\infty$ bounds, these minimizers converge to a locally energy-minimizing harmonic map $\mathbf{Q}_0$ into a biaxial vacuum manifold. The main result of this paper is that such a limiting map $\mathbf{Q}_0$ has no interior point singularities. The proof relies on a geometric identification of the lifted Frobenius metric on the universal cover $\mathbb{S}^3$ with a rescaled Berger metric. For a hypothetical tangent cone at a point singularity, its link is a nonconstant harmonic two-sphere into the Berger sphere. We construct a smooth variation field adapted to the Hopf direction and show that it induces a quantitative instability estimate, in contradiction with the shifted stability inequality inherited from local minimality. This replaces the usual round-sphere test fields by a target-specific construction and rules out interior point defects in the biaxial setting. The result is sharp in view of the well-known existence of point defects in the uniaxial theory.

We introduce a variational approach to a finite-temperature continuous-spin perceptron trained on a Gaussian mixture. The model allows for a broad class of concave utilities and log-concave separable prior measures on the spins. By combining the interpolation method with log-concavity and concentration estimates, we derive lower and upper minimax variational bounds for the limiting quenched pressure. Remarkably, the two bounds differ only in the order of optimization of two variational parameters, while all remaining extrema are controlled by the concave–convex structure of the variational potential. Whenever the two optimizations commute, the two bounds match and identify the solution of the model. The same potential yields the fixed-point equations as stationarity conditions and provides a unified route to the computation of the ground-state energy, training loss, and generalization error.

Although Arnoldi reduction of a generally non-Hermitian Hamiltonian yields an upper Hessenberg matrix rather than the tridiagonal form of Hermitian Lanczos theory, we show that a closed Toda sector survives in its diagonal and subdiagonal coefficients. For a fixed finite-dimensional Hamiltonian and a cyclic state vector deformed holomorphically, the Krylov Gram determinants are $\tau$ functions of the finite two-dimensional Toda lattice, whose Flaschka variables coincide exactly with these Arnoldi coefficients. The Toda dynamics therefore closes on this sector without determining the remaining upper Hessenberg entries. The subdiagonal part of the same sector also has a direct geometric meaning: the squared subdiagonal coefficients determine both the Fubini–Study metric and the Berry curvature of holomorphic Krylov subspaces, whereas the geometric quantities associated with subspaces lost at Arnoldi breakdown cease to be defined. Along a smooth real path in the cyclic region, the Arnoldi-frame connection further provides a Hermitian tridiagonal generator of exact isospectral transport. When added to the Arnoldi matrix, this generator cancels transitions between instantaneous eigenspaces and realizes counterdiabatic driving whenever the matrix is diagonalizable with a nondegenerate spectrum.

The main objective of this work is to investigate the traveling wave solution and dynamic characteristics of the (2 + 1)-dimensional seventh-order Caudrey-Dodd-Gibbon-KP (sCDG-KP) equation. Applying the ($\frac{G^{\prime}}{G^{\prime}+G+A}$) method, we examine the exact solution of the (2 + 1)-dimensional seventh-order Caudrey-Dodd-Gibbon-KP (sCDG-KP) equation by altering it into a reduced ODE via a suitable wave transformation. Graphical representations, such as 2D, 3D, and a heat map of the ascertained solution, are present to facilitate comprehension of the empirical relevance of the obtained solutions. As a result, we acquired a bright and anti-kink soliton solution. Next, we alter the ODE into a 2D system of equations to analyze the dynamical behavior of the reduced system via bifurcation analysis, phase portrait, and attractor analysis. During this process, we portray the graphical visualization of the bifurcation phase portrait, 2D phase portrait, 3D phase portrait, time series, chaotic attractor, sensitive analysis, fractal dimension, recurrence plot, and power spectrum of the dynamical system.

Reconstructing unknown electromagnetic sources from far-field radiation patterns is a fundamental inverse problem with broad applications in biomedical imaging, non-destructive testing, and telecommunications. In practical settings, however, collecting dense multi-frequency far-field measurements at the Nyquist sampling rate is often infeasible. Under-sampled or sparse data introduce non-radiating source components that sever the uniqueness of the solution, creating severe artifacts when standard inversion techniques are applied. To overcome this limitation, we present a two-stage reconstruction strategy exploiting the physical property that compactly supported, geometrically sparse sources exhibit a finite rate of innovations (FRI). In the first stage, we construct an associated wrap-around structured Hankel matrix. By leveraging the low-rank property of the matrix due to FRI of the unknown sources, we enrich the sub-sampled data. To that end, we convert missing multi-frequency far-field data recovery into a constrained matrix completion task solved via Annihilating Filter-based Low-rank Hankel Matrix Completion Approach (ALOHA). In the second stage, a Fourier inversion scheme reconstructs the current source density from the enriched dataset. Extensive numerical evaluations on electromagnetic source models show that our enrichment framework effectively eliminates under-sampling artifacts and resolves non-uniqueness challenges. The method delivers accurate and stable reconstructions under high sub-sampling rates (e.g., with $30$% to $50$% available samples) and strong noise conditions ($10$ dB SNR), outperforming standard $\ell_1$-compressed sensing baselines.

We study a two-dimensional semiclassical Schrödinger operator whose potential admits one reflection symmetry. Under a rational independence assumption on the harmonic frequencies and a nonvanishing condition \(a_{30}\ne0\), we show that the first two layers of the quantum Birkhoff normal form, together with the sign of \(a_{30}\) and the transverse data \(\{a_{1,2k}\}_{k\ge1}\), determine the Taylor series of the potential at the bottom of the well. The proof is constructive: the \(\hbar^2\)-layer gives a triangular recursion for odd-degree terms, while the classical layer recovers even-degree terms from their resonant projections.