Papers

In the whole space of low spatial dimensions, namely $d \leq 2$, we study the minimizers of an energy functional with an attractive cubic nonlinearity and a repulsive quintic nonlinearity, which describes a quantum Bose gas with a two-body attraction and a three-body repulsion. We prove the existence of minimizers at fixed effective statistics parameter. In the limit of a large effective statistics parameter of the mass-critical nonlinearity (with respect to the kinetic energy), we derive an effective Thomas-Fermi-like model for the homogeneous Bose gas. We also consider other limit regimes depending on the mass-critical nonlinearity.

We consider a homogeneous Bose gas composed of $N$ identical bosons occupying an infinite two-dimensional space. This gas exhibits both an attractive two-body interaction and a repulsive three-body interaction. Via the intermediate Hartree theory, we rigorously derive the homogeneous cubic-quintic nonlinear Schroedinger functional as the mean-field limit of the model. Our investigation focuses on the system's behavior in relation to the two-body interaction.

We investigate the behavior of a homogeneous one-dimensional Bose gas composed of $N$ identical bosons experiencing attractive two-body and repulsive three-body interactions. Utilizing intermediate Hartree theory, we rigorously derive the homogeneous cubic-quintic nonlinear Schrödinger functional as the mean-field limit of the many-body Hamiltonian. Particular attention is given to the impact of the repulsive three-body term on ground state energy estimates and mass concentration properties

The Virasoro conjecture predicts a set of universal relations among all genera Gromov–Witten invariants of any smooth projective variety. The conjecture is well understood for semisimple theories, but remains largely open in the non-semisimple setting. We prove the genus-one Virasoro conjecture on the ambient state space of smooth Fano complete intersections in projective space. For most of these complete intersections, the big quantum cohomology is nowhere semisimple. We also generalize the wall-crossing formula for quasimap invariants with weighted markings to the equivariant twisted setting, allowing descendant insertions at light markings. Together with genus-one quantum Lefschetz for quasimaps with light markings, this wall-crossing formula provides the key bridge from the Gromov–Witten theory of the complete intersection to the semisimple equivariant twisted theory of the projective space.

Quantum resonances and black-hole quasinormal modes (QNMs) are described by closely related mathematics: radiative boundary conditions, analytic continuation of Green functions or resolvents, poles on nonphysical sheets, and non-self-adjoint spectral representations. This article does not introduce that equivalence but reconstructs its genealogy. Nuclear resonance theory grew from radioactive decay, reaction cross sections, metastable compound states, and complex energies, whereas black-hole perturbation theory grew from spacetime stability and causal response. We place the two developments on a common chronology, from Sommerfeld's radiation condition and early nuclear-decay theory through Regge–Wheeler perturbations, Vishveshwara's 1968 dissertation and 1970 papers, Press's 1971 quasi-normal terminology, and the independent Aguilar–Balslev–Combes theory of analytic dilation, then follow their later mathematical reunification in scattering-resonance theory. Alongside the chronology, we give a working dictionary relating no-incoming and horizon-ingoing/infinity-outgoing conditions, pole energy and complex frequency, width and damping rate, residues and excitation factors, and nonresonant background with branch-cut and prompt contributions. The aim is to make the cross-disciplinary equivalence operational without flattening its history, and to distinguish mathematical equivalence from historical genealogy.

Poisson algebras are the quasiclassical limits of associative algebra deformations of commutative associative algebras. This paper extends this process to the level of bialgebras. We derive Poisson bialgebras as the quasiclassical limits of antisymmetric infinitesimal bialgebra deformations of commutative and cocommutative antisymmetric infinitesimal bialgebras. It might be regarded as the “infinitesimal" version of the quantization process of Lie bialgebras in terms of Hopf algebras. Such deformations-quasiclassical limits process for a Poisson bialgebra is equivalently characterized in terms of the introduced notions of deformations of a matched pair of associative algebras as well as a standard Manin triple of associative algebras, whose corresponding quasiclassical limits are a matched pair of Poisson algebras and a standard Manin triple of Poisson algebras, respectively. We illustrate these equivalent deformations-quasiclassical limits processes via coherent derivations.

We develop simplicial-set models for dagger $(\infty,1)$-categories. We first establish a model structure on dagger simplicial categories whose weak equivalences are local weak equivalences that are essentially surjective up to coherently unitary equivalence. We then construct the dagger Joyal model structure on dagger simplicial sets and prove that the dagger rigidification–nerve adjunction is a Quillen equivalence. We next introduce the unitary core and use it to characterize dagger Joyal equivalences intrinsically. Finally, we compare the dagger Joyal model structure with the DCH anti-involutive Joyal model structure.

We introduce latent geometry in BFSS/BMN models. Before an emergent geometry evaporates at the geometric transition, large-$d$ Gaussianization encodes its information in the exceptional $N=2$ sector selected by the baby fuzzy sphere, a single active spin-$1/2$ Pauli block. Dimensional reduction, finite-$N$ Gaussianization and exact Molien–Weyl projection transport this encoding to mass-deformed BFSS$_2$. At $N=2$, the same stable BFSS$_2$ oscillator admits inequivalent unitary $\mathbf{H}^2_\theta=\mathbf{EAdS}^2_\theta$ and $\mathbf{dS}^2_\theta$ decodings, realizing spacetime transmutation and the emergence of Lorentzian time.

The bosonic Kitaev chain under open boundary conditions has attracted recent attention due to its realization in driven-dissipative systems and its intriguing non-Hermitian boundary physics. The model is known to be solvable via local squeezing transformations in the position-momentum representation. In this paper, we present an alternative, purely algebraic solution that relies entirely on the standard bosonic Bogoliubov transformation. For an $N$-site chain, we propose an asymmetric plane-wave ansatz with unequal left- and right-moving momenta to analytically solve the associated $2N\times 2N$ non-Hermitian “associated matrix". The left eigenvalue problem yields $N$ distinct eigenvalues, each of which is twofold degenerate. By carefully resolving these degeneracies using the bosonic commutation relations, we construct the $N$ physical Bogoliubov quasiparticle operators. The construction reveals non-uniquenesses that in special cases exactly mirror the freedom in the local squeezing transformations of the original approach. The diagonal form of the Hamiltonian is obtained explicitly and is shown to be equivalent to the original Hamiltonian. The proposed asymmetric plane-wave ansatz and degeneracy-resolution technique are not limited to the present model and can be generalized to other bosonic pairing systems, including those with inhomogeneous pairing or hopping terms.

We develop a geometric framework for implicit discrete Hamiltonian systems based on discrete Morse families, Lagrangian relations, and discrete analogues of Tulczyjew's triple. The main idea is to regard the Lagrangian submanifold defining the discrete dynamics, rather than an explicit symplectic evolution map, as the fundamental geometric object. This viewpoint naturally accommodates implicit, constrained, and degenerate discrete systems. Within this framework, we formulate a Type–II discrete Hamilton–Jacobi theory in terms of the propagation of Lagrangian submanifolds between consecutive discrete steps. When these submanifolds are locally represented by exact one-forms $dW_k$ and $dW_{k+1}$, the resulting equations provide a discrete Hamilton–Jacobi relation between consecutive generating functions. More generally, when the projection onto configuration space becomes singular and a single-valued generating function ceases to exist, we show that the evolution can be described by the composition of Type–II discrete dynamics with Morse families. This yields a generating family for the propagated Lagrangian submanifold without requiring the dynamics to be represented as a graph. As an application, we consider the propagation of optical wavefronts through fold caustics. A Type–II discrete Hamiltonian yields a symplectic ray integrator, while a Morse family represents the multivalued wavefront near the caustic. Their composition provides a discrete propagation rule for the complete Lagrangian manifold across the singularity. In this way, the same geometric construction simultaneously provides a discrete Hamiltonian integrator and a regular representation of multivalued Hamilton–Jacobi solutions.

In Einstein gravity, every null geodesic of a non-static Bianchi I spacetime satisfying the null energy condition is incomplete in at least one direction within the Bianchi I region. We prove a stronger directional rigidity theorem: a single null geodesic complete in both affine directions forces staticity if the null curvature is nonnegative along it and, at each time, in one transverse null direction. The spatial metric need not be diagonal, and the expansion rates may have either sign. The proof uses caustic-free null congruences whose conserved spatial covectors define fixed screen hyperplanes. The two curvature conditions admit averaged replacements along the geodesic and, when cosmic time is complete, along comoving worldlines. In a non-static metric, at most four oriented directions are optically saturated; every other complete null geodesic has strictly negative averaged null curvature. Periodic and finite-pulse examples attain this bound with NEC-violating sources and nonpositive total energy density. For finite pulses, we derive bounds on metric excursions from three directional curvature budgets at fixed cosmic duration, and on the mode mixing of a free, massless, minimally coupled test scalar field from the averaged null-energy deficit at fixed affine duration.

In an inhomogeneous dissipative medium, different rays experience different round-trip attenuation, so an instantaneous time mirror may retrace rays without producing a balanced point focus. We solve this problem for a microlocally selected back-propagating transverse-magnetic component generated by a short pulsating modulation of a thin shell surrounding an inhomogeneous Debye medium with Ohmic conductivity. An exact decomposition of the Maxwell–Debye field isolates the pulse-generated return, and a semiclassical mode reduction yields a scalar returned operator whose principal symbol separates the temporal pulse response, shell modulation, and accumulated material loss. Our analysis shows that geometric distortion and amplitude attenuation are both involved in wave propagation, but the geometric distortion cancels between the outgoing and returning branches, whereas physical attenuation remains. A normalized directional energy, computed from pre-modulation measurements using the known lossless backward propagator, contains the same leading attenuation factor as the returned symbol. Its reciprocal defines a smooth, nonnegative, physically admissible modulation that equalizes the selected source-side gain up to a relative $O(h)$ error without requiring a complete attenuation model. With such a modulation, the field generated by a point source has a focal-scale leading profile independent of directional attenuation, with a unique maximum at the source and coherent return time. In every fixed focal window, the maximum is displaced by at most $O(h^2)$ in space and time, and the spatial half-amplitude area is $O(h^2)$.

Strong topological insulators are classified by integer invariants which admit different real-space expressions. Inspired by recent work studying these invariants for spaces with curved boundaries, we revisit the problem of equivalence of the various expressions, including Fredholm index pairings and Kubo formulas involving half-space projections. Using Roe algebras and a variant of KK-theory suitable for non-separable $C^*$-algebras, we prove a general form of bulk-edge correspondence for the index pairing with a general position-space Dirac operator in the presence of arbitrary boundaries. For spaces coarsely equivalent to $\mathbb{R}^d$, we show that up to multiplicity, these pairings are equal to those obtained with the standard dual Dirac operator, with the multiplicity explicitly given as the topological degree of the symbol function. In addition, we prove a generalized Kubo formula which computes the general index pairing by a real-space formula.

Physics-informed neural networks (PINNs) have emerged as a meshless alternative to conventional numerical methods for solving partial differential equations (PDEs). However, their limited ability to capture sharp gradients can lead to substantial errors when resolving boundary and interior layers. Here, we introduce an energy-norm-enhanced PINN (ENPINN) that incorporates gradient information and variational structure into the loss function to improve the resolution of layer-dominated solutions. We first examine two related formulations: weak-loss PINNs (WLPINNs), which incorporate test functions into the conventional PINN residual, and gradient-enhanced PINNs (gPINNs), which augment the loss with spatial derivatives of the PDE residual. By analyzing these formulations, we identify their limitations in resolving steep solution gradients and motivate the systematic construction of ENPINN. We establish theoretically how the energy-norm error depends on the ENPINN loss and show that a suitably modified residual-derivative term is essential for accurately capturing boundary layers. We further establish the existence of neural-network approximations with arbitrarily small energy error and derive corresponding derivative bounds, providing a theoretical foundation for the proposed framework. The performance of ENPINN is assessed through systematic comparisons with existing PINN variants for convection-diffusion-reaction problems exhibiting steep gradients. Numerical experiments include a combustion model, a coupled multi-scale system, a two-dimensional Burgers equation with an interior layer, and a three-dimensional time-dependent problem.

We develop a representation-independent symmetric-Fock formulation of the nonlinear Boltzmann collision operator for Maxwell molecules. Starting from the standard velocity-space operator, we give a short derivation of Bobylev's Fourier identity and isolate its characteristic structure: two linearly related Fourier arguments followed by multiplication. We then prove, independently of collision geometry and of any coordinate realization, a canonical lift–fusion intertwining theorem: one-particle linear maps lift to the two incoming Fock legs and symmetric-algebra multiplication fuses them into a single outgoing state. Bargmann coordinates make this intrinsic construction especially transparent, but do not define it; Bobylev's collision maps subsequently select the Maxwell vertex from the already graded lift–fusion class. The resulting bilinear map obeys the exact grading \[ \widehat{\mathcal Q}(\mathscr H_p,\mathscr H_q)\subseteq\mathscr H_{p+q}, \] which gives an intrinsic origin for the triangular structures of Maxwell kinetics. Its one-vacuum restriction yields the complete Wang Chang–Uhlenbeck spectrum, while moment adaptation makes the nonequilibrium hierarchy strictly triangular and every finite compatible Fock jet an exact autonomous factor with global exponential relaxation. For the inhomogeneous problem, we extend realization covariance to the nonlinear two-leg vertex and evaluate the moving-frame connection independently through Hermite-function and Fourier realizations. The two routes give the same abstract Fock propagation and the same first Chapman–Enskog source; its traceless level-$2$ and vector level-$3$ sectors are relaxed by the exact Maxwell rates, giving $\Pr=2/3$. This provides an explicit computational demonstration that coordinate realizations may be chosen for economy without entering the structural Fock formulation.

Let $f:(0,\infty)\to(0,\infty)$ be a nontrivial normalized operator monotone function and set $s=f'(1)$. We introduce the alternative-mean trace functional $$ \altPhi_f(A,B) :=\Tr(A\nabla_s B) -\Tr\!\left( f(A^{-1}\sharp B)\,A\,f(A^{-1}\sharp B) \right). $$ on the positive definite cone. We prove that $\altPhi_f$ is a quantum divergence in the sense of Bhatia–Gaubert–Jain whose diagonal Hessian induces a positive multiple of the Bures–Wasserstein Riemannian metric. We also establish the sharp comparison $$ s(1-s)d_{\rm BW}(A,B)^2 \le \altPhi_f(A,B) \le (1-s+s^2)d_{\rm BW}(A,B)^2. $$

Let $(M,g)$ be an asymptotically flat manifold such that its tangent bundle $TM$ is diffeomorphically asymptotically Euclidean. We prove a positive mass theorem and a positive mass rigidity theorem for a natural class of asymptotically flat metrics $\widetilde g$ on $TM$ exhibiting slow asymptotic decay. These metrics are constructed by interpolating, along the horizontal distribution of $TM$, between the Euclidean metric and the base metric $g$ near infinity. The main difficulty is that the induced metrics on $TM$ decay below the standard threshold for the ADM mass in dimension $2n$; nevertheless, we show that their asymptotic mass is well-defined in a generalized sense and is determined by geometric data on the underlying manifold $M$.

The Cardy-Smirnov function of a four-marked planar domain encodes the conjectural scaling limit of the crossing probability for critical percolation. Cardy predicted a closed formula for this limit on the rectangle in the cases of both site and bond percolation, and Smirnov proved the formula together with conformal invariance for critical site percolation on the triangular lattice. Kleban and Zagier later showed that on the rectangle the function admits a modular interpretation: it is determined by a modular functional equation together with a mild analytic ansatz. We carry out the analogue on the $\pi/3$ parallelogram. We derive a closed conformal-map formula for the Cardy-Smirnov function as an incomplete beta integral, prove a closed modular formula realizing it as an integral of the modular form $\eta(\tau)^2\eta(3\tau)^2$, and establish a uniqueness theorem showing that a single functional equation, together with a $q$-expansion ansatz and a nondegeneracy condition, determines the function uniquely.

We study a null model of high-dimensional projection pursuit: we are given $M$ points sampled i.i.d. from a standard gaussian in $N$ dimensions, where $M,N\to\infty$ with $M/N\to\alpha\in(0,\infty)$. Our goal is to characterize the possible empirical distributions of these points' projections along a data-dependent direction $x$, which ranges over either the sphere $S_N=\sqrt{N}\mathbb{S}^{N-1}$ or cube $\Sigma_N=\{-1,+1\}^N$. We consider this problem in an algorithmic setting, where $x$ must be the output of an algorithm with dimension-free Lipschitz dependence on the input; this class of algorithms includes general gradient-based methods such as Langevin dynamics and approximate message passing (AMP). Our main result exactly characterizes the set of empirical distributions attainable by this class in terms of a one-dimensional stochastic control problem. As a consequence of our main result, we obtain exact algorithmic thresholds for optimizing the Hamiltonian of a spherical or Ising perceptron model with general bounded continuous activation. For the spherical problem, independent work of Montanari and Zhou (2024) characterized the empirical distributions attainable by a related two-stage AMP algorithm, also in terms of stochastic control. Our proof of hardness builds on the branching overlap gap property introduced in earlier work by the first two authors. Our main innovation is to develop stochastic control theory within the branching OGP framework, significantly expanding the settings in which it locates an exact algorithmic threshold. Notably, our methods apply even though the non-algorithmic problem of characterizing all feasible projections remains a major outstanding challenge. For the matching algorithmic result, we construct a new incremental AMP algorithm that acts on a Brownian-bridge revelation of the gaussian disorder and simulates the same family of controlled SDEs.

This article presents a systematic study of the critical cohomology of $\QM^{\xi}(\bbP^1, X)$, the moduli space of based quasimaps from $\bbP^1$ to a Nakajima quiver variety $X$. We develop two derived equivalent presentations of this moduli space as a global critical locus. Accordingly, we obtain two equivalent explicit realizations of the canonical DT sheaf of the moduli space. Using each critical presentation, we describe canonical actions of shifted Yangians and quantized Coulomb branches on the quasimap critical cohomologies. We produce, by geometric means, a canonical morphism from an appropriate Hall algebra to the Coulomb branch, and show that the two actions we describe are compatible via this morphism under an additional hypothesis. We conjecture the compatibility holds generally, and reduce the claim to a sheaf-theoretic statement that may be approached via Joyce conjecture.

This is the first publication, in which an inverse problem of the simultaneous reconstruction of an unknown (x,t) - dependent coefficient and an unknown initial condition in a general linear parabolic equation of the second order is considered. The input data depend on (n+1) variables and are, therefore, formally determined ones. Stability estimates and uniqueness theorem are obtained for this inverse problem. As a by-product, logarithmic stability estimates for initial conditions of parabolic equations and inequalities with time reversed data are obtained for the first time. All known publications about inverse problems for parabolic equations with formally determined input data and unknown coefficients assume that those coefficients depend either only on x or only on t. In addition, the initial condition is assumed to be known, except of two publications cited in the text. The inverse problem of this paper has two potential applications. The first one is in forecasting of public opinions in the framework of the Mean Field Games theory. The second one is in tracking spatiotemporal outbreaks of epidemics.

On the space of density matrices $D(H)$, Kadison's theorem establishes that all invertible convex transformations $K:D(H)\to D(H)$ are unitarity or antiunitary. For a finite Hilbert space, we propose a similar theorem which replaces the invertibility requirement with that of entropy conservation. Thus, we argue, the supposition of reversibility in quantum mechanics may be replaced with a principle of entropy conservation.

A parametrized quantum circuit reports its state geometry through a quantum Fisher information matrix (QFIM), often singular. A small Fisher value can reflect exact state-preserving redundancy, compression by the circuit chart, or weak intrinsic distinguishability, and these mechanisms call for different numerical treatments. We show that the circuit metric factors as $F=B^{*}MB$, where $B$ is the state-level circuit-to-orbit differential and $M$ is the intrinsic Fisher operator on the reachable orbit. The factorization identifies the exact kernel as $\ker B$, separates coordinate transfer from intrinsic geometry, and yields the condition for a circuit to realize an orbit-level quantum natural-gradient (QNG) direction. When the prescribed redundancy exhausts the Fisher kernel, the Moore–Penrose update is the minimum-norm horizontal lift of the quotient Riemannian gradient. At critical points with a locally diffeomorphic quotient-to-orbit map, chart singular values cancel from the linearized QNG operator while intrinsic anisotropy remains; in the trace-orthonormal full-control generator frame, excitation-gap anisotropy gives $\kappa_{\mathrm{QNG}}=\kappa_{\mathrm{Eucl}}$. Representation theory makes $M$ explicit on highest-weight, Slater, and fermionic-Gaussian orbits, and cominuscule fidelity flow becomes integrable, with conserved principal-defect ratios, cubic Lie-retracted convergence at $\eta=2$, and stability boundary $\eta=4$. Finite data impose a second boundary: an estimated QFIM and its confidence radius alone cannot distinguish an exact zero from a small physical mode, so the estimated spectrum alone cannot license hard projection. Under depolarization, inverse-Fisher scaling amplifies mean updates and fluctuations together and cannot restore update signal-to-noise. A redundant Slater/Givens circuit confirms exact transfer identities and illustrates finite-shot tradeoffs.

Dynamical localization requires all position moments of a quantum wavepacket to remain bounded in time, but for quasiperiodic Schrödinger operators such bounds are generally not uniform in phase. In the positive Lyapunov exponent regime, the best known phase-uniform estimates instead grow on a logarithmic scale. We prove that both the logarithmic scale and the dependence on the moment order are sharp for a class of one-frequency quasiperiodic Schrödinger operators with even potentials. Our main ingredient is a reflective version of semi-uniformly localized eigenfunctions, adapted to the two localization centers forced by a completely resonant phase, from which we obtain matching logarithmic lower bounds along sequences of times.

We present a systematic Eulerian study of passive scalar advection in rough two-dimensional Kraichnan flows, covering the full range of velocity roughness exponent $h\in(0,1)$. The advection–diffusion equation is integrated directly using a pseudo-spectral method at resolutions up to $2048^2$ grid points, with a frozen-noise Runge–Kutta scheme consistent with the white-in-time construction of the carrier flow. The simulations recover the duality between advected scalar and advecting flow—the smoother the carrier, the rougher the scalar—together with the predicted scaling laws for the second-order statistics. Once the scaling range is properly identified, the fourth-order flatness anomaly is measured across the whole range of $h$, in agreement with the Lagrangian estimates of Frisch et al. (1999) and with the perturbative predictions of Bernard et al. (1998) and Pumir et al. (1997). Benefiting from the fact that our Eulerian simulations give direct access to the full scalar field, we also examine the probability density functions of scalar increments and the corresponding higher-order statistics, which show systematic departures from Gaussianity and from log-normality, with the strongest deviations manifesting for intermediate values of $h$. A central outcome of this work is a systematic account of how the simulation parameters, in particular the molecular diffusivity, must be adjusted with $h$, providing practical guidelines for reliable simulations; we further show that the residual deviations from the theoretical scaling laws are quantitatively accounted for by the finite spectral representation of the carrier flow. Our analysis also serve as a numerical baseline for simulating more realistic extensions of the Kraichnan model, where the carrier flow is coupled with a Gaussian multiplicative chaos and theoretical developments are limited.

This paper studies global existence for the Cauchy problem of nonlinear Klein-Gordon equations in three space dimensions, strictly within the regularity regime of classical local existence. We prove it via higher-order and lower-order energy estimates. The proof relies on two key ingredients. The first is due to the work of Georgiev and Popivanov, which reduces quadratic nonlinearities to cubic terms plus ghost-energy-controllable terms, securing lower-order estimates. The second is a sharp Klainerman-Sobolev-type inequality without the scaling operator, established herein, which yields enough time decay for derivatives up to second order; integration by parts then controls derivative loss terms in the higher-order estimates.

Synchronization is conventionally understood as the emergence of a phase-locked collective state among interacting dynamical units. However, extending this notion to systems whose dynamical variables are defined on higher-dimensional simplices introduces fundamentally new constraints arising from the topology of the underlying simplicial complex. In the simplicial Kuramoto model, nontrivial topological cycles give rise to a higher dimensional harmonic subspace that is unaffected by the coupling and can therefore drift indefinitely, preventing a globally synchronized state. To resolve this we employ Hodge decomposition on simplicial Kuramoto dynamics and investigated the components. Despite this topological obstruction, we show that the simplicial Kuramoto model admits fixed-point states in the exact and coexact sectors of the Hodge decomposition above critical coupling strengths, which we derive analytically for both sectors. This decomposition provides a generalized notion of synchronization in which the non-harmonic components converge to fixed states. We further investigate the robustness of these fixed-point states to external perturbations. Excluding the drifting and non-interacting harmonic component, the perturbation response is governed by the spectrum of weighted Laplacian and recovers a Kirchhoff index dependence analogous to standard Kuramoto case. In contrast, the harmonic sector is non-dissipative, and consequently the fragility of the system increases with the dimension of this topological subspace. We characterize this effect numerically using triangulated tori with varying first Betti number and find a superlinear scaling relation between system fragility and the dimension of the harmonic sector.

We prove an all genus full descendant open mirror symmetry for footballs. The B-model is given by the Chekhov-Eynard-Orantin topological recursion on the mirror curve.

We introduce a new class of commutative associative algebras whose structural constants are linear operators acting in a linear topological space. An analogue of the theory of monogenic functions is developed for such algebras. It is shown that the components of monogenic functions satisfy second-order operator-differential equations. This approach provides a constructive method for obtaining solutions of various partial differential and integro-differential equations. Several illustrative examples are presented.

We introduce a new class of algebras with variable operator-valued structural constants and develop elements of the theory of monogenic functions in such algebras. Compatibility conditions for the structural operators are established. It is shown that the components of monogenic functions satisfy second-order operator-differential equations with variable coefficients, and converse results are proved. The proposed approach provides a constructive method for obtaining solutions of various differential, integro-differential, and functional-differential equations. Several illustrative examples are presented.