Papers

We present SymQuPS, a SymPy-based algebra system in Python mainly aimed toward the phase space representation of quantum mechanics within the Cahill-Glauber formalism (including the Glauber-Sudarshan $P$, Wigner, and Husimi $Q$ representations). By extension, the package serves as an algebraic venue for canonical quantization. A key feature is the phase space representation of an arbitrary Lindblad master equation, which gives the phase space equation of motion of the quantum system. We describe the core functionalities of the package, consisting of $s$-ordered operators, the star products, and the phase space representation. Some examples of use are given to illustrate the application of the package, and the package's performance in typical use cases is discussed.

We develop the kinematical and operator-algebraic foundations of a groupoid-based relational quantum field theory (RQFT) on curved spacetimes. Indeed, the usual group-based quantum reference frame (QRF) formalism is not directly suited to generic curved Lorentzian backgrounds as global symmetry groups are typically absent or too small. We formulate a notion of QRF for a continuous groupoid. This yields a groupoid relativization map and relational observables. We show that a localization limit recovers the ordinary non-relational description. We construct canonical sharp groupoid QRFs, which form the groupoidal counterpart of the ideal group QRFs based on $L^2(G)$. We further prove that the groupoid QRF construction reduces to the standard operational QRF formalism for locally compact groups. The action groupoid QRFs are torsor QRFs only for specific classes of fields of positive operator-valued measures (POVMs), and the torsor relativization map only applies to constant operator fields of system observables. We review the foundations of RQFT in Minkowski spacetime. We prove new results that further link RQFT to Wightman QFT. We show that covariant POVMs are $\mu$-continuous with respect to quasi-invariant $\sigma$-finite positive Borel measures $\mu$. Thus, relational quantum fields can be understood as the smearing of pointwise-defined kernels with respect to the QRF's statistics. We develop RQFT in curved spacetime, where we argue that the correct replacement for the Poincaré group is the Poincaré groupoid of the spacetime. We also indicate how the framework extends further to internal gauge symmetry and relational gauge-covariant quantum field theory via Atiyah groupoids, providing a first step towards formulating a relational quantum Yang-Mills field theory.

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

We introduce cleft extensions for Hopf algebroids. We prove the equivalence between cleft extensions, $\sigma$-twisted crossed products, and Hopf-Galois extensions with the normal basis property, thereby generalizing the theory of cleft extensions for Hopf algebroids developed by Böhm and Brzeziński, and fitting in with the general theory of Galois and biGalois extensions over Hopf algebroids developed by the authors. We investigate the Ehresmann Hopf algebroid associated with a cleft extension and show that it is isomorphic to a generalized version of the Connes-Moscovici Hopf algebroid. A special case of the Connes-Moscovici Hopf algebroid, namely the case where the coinvariants of the cleft extension coincide with the base of the Hopf algebroid, is a Drinfeld twist of a Hopf algebroid by a two-cocycle, generalizing work of Böhm, Han and Majid.

We extend de Broglie–Bohm (dBB) pilot-wave theory by introducing a complex source field into the Schrödinger equation and examining its effects on quantum equilibrium and nonequilibrium dynamics. In dBB theory, the physical particle distribution P, rather than the Born density $|{\psi}|^2$, carries the ensemble probability, so a source term may modify the pilot wave without violating conservation of total particle probability. We derive the source-modified Hamilton–Jacobi and continuity equations, the transport equation for the nonequilibrium ratio $f=P/|{\psi}|^2$, and an exact entropy-production formula. Three applications follow. First, suitably designed sources can drive exponential relaxation toward quantum equilibrium. Second, the entropy-production rate admits a Prigogine-type bilinear form, providing a basis for a subquantum thermodynamics with entropy-producing and entropy-extracting regimes. Third, a tuned source can exactly cancel the Bohmian quantum potential, yielding classical particle trajectories while the guiding wave retains nontrivial structure and the nonequilibrium ratio remains conserved. The resulting framework provides a unified setting for studying source-driven quantum nonequilibrium, entropy exchange, and classicalization, and motivates further investigation of the physical and ontological status of the source field.

The Carrier–Greenspan hodograph transformation linearises the nonlinear shallow-water equations on a beach of constant slope. We extend the transformation to a moving bottom and identify the family of bottom motions that preserves the Carrier–Greenspan invariants. The bottom gradient must be uniform in space, that is, $h(x,t)=\rho(t)x+\mathcal{B}(t)$, and constant $\rho$ returns the classical slope. For this family a horizontally accelerating frame removes the forcing. The bottom motion then disappears from the hodograph equation, and two quadratures recover the physical variables. If the gradient varies in time, the bed displacement grows without bound offshore and the far field cannot remain at rest. The domain must therefore be finite, although the pivot may be placed anywhere. A bottom-tilting wave maker realises this family. Closed forms follow for the hinge response, for the crest-to-trough steepness of the generated wave, and for the leading nonlinear correction. The correction depends on the point to which the steepness is referred, and so differs between a leading-elevation and a leading-depression wave. Because the hinge is at the beach toe, these closed forms provide the incident wave required by the classical plane-beach solution. The run-up can therefore be predicted from the plate motion alone. The steepness is predicted without a fitted constant and agrees with the reference computations. Sixty-one runs over a level bed, reported here for the first time, carry the test to deeper water and longer plate motions. The measured run-up follows the predicted run-up, with an offset that the inviscid theory does not account for.

Graphons are measurable functions used to describe the asymptotic behavior of convergent graph families. Originally motivated by problems in combinatorics and graph theory, graphons have found numerous applications in the modeling and analysis of dynamical processes on networks. In this work, we use graphons to formulate the Ising model on convergent graph sequences, which include many network topologies common in applications. We derive the mean-field limit for the resulting model and obtain exact results for phase transitions in such systems. Specifically, we show that the critical temperatures of the Ising model on graphons are determined by the eigenvalues of the Hilbert-Schmidt operator associated with the graph limit. For many important network topologies, these eigenvalues can be computed explicitly. We illustrate our results with three representative random network models: Erdős-Rényi, small-world, and power-law. In the small-world case, we demonstrate phase transitions to both ferromagnetic and antiferromagnetic phases, as well as coexistence of local minima of the free energy. The latter gives rise to multistability, as confirmed by Monte Carlo simulations. The results of this work demonstrate that the Ising model on graphons combines the analytical tractability of exactly solvable mean-field models with the ability to accommodate a broad range of network topologies. We expect that the use of graphons in spin models will lead to new insights into the statistical physics of interacting systems on complex networks.

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

Channel–state duality identifies completely positive maps with bipartite states. Its generalisation to von Neumann algebras is known to produce states not on $\mathcal{M}\bar{\otimes}\mathcal{N}$ but on the tensor product of $\mathcal{M}$ with the opposite algebra of $\mathcal{N}$, or equivalently with a commutant; returning to $\mathcal{N}$ itself requires a transposition, and it has been observed that this step typically fails. We determine exactly when it does not fail. Our main result is that a Choi–Jamiołkowski type isomorphism which is natural in the first algebra exists if and only if the second algebra is $*$-anti-isomorphic to itself. Naturality cannot be dropped: whenever $\mathcal{M}$ is anti-isomorphic to itself such an isomorphism exists for trivial reasons, since the order structure of a predual is a Jordan invariant and cannot distinguish an algebra from its opposite. Consequently no such correspondence exists for the type III factors constructed by Connes. We show in addition that the domain admits no description in terms of properties of individual maps: the space of all normal maps with the operator norm is too large, and neither complete boundedness nor complete positivity cuts it down to the right size.

Everybody learns in school that an out-of-equilibrium system coupled to a heat bath at a fixed temperature evolves to thermodynamic equilibrium as time goes on, and the free energy will only decrease on its way there. At least since Holley's 1971 work, mathematicians know this too, in the modest context of classical Ising Glauber dynamics. But does the free energy also asymptotically decrease to the free energy of an equilibrium state? A typical school kid would say "yes, of course", but since the free energy is only lower semicontinuous, this question is less straightforward than it initially seems. To the best of our knowledge, apart from the uniqueness regime, where one can make use of classical functional inequalities, this question has not previously been resolved rigorously. We prove the asymptotically complete dissipation of the free energy for the classical Ising Glauber dynamics by introducing a coarse modified log-Sobolev inequality, which holds at every positive temperature, in particular in the phase-coexistence regime.

Bipartite qubit entanglement distillation is the process of converting noisy ebits into pure ebits using only local operations and classical communication. This is a core operation for quantum repeaters, enabling such crucial tasks as long-distance quantum communication and distributed quantum computing. In this work, we introduce a method for searching for entanglement distillation protocols and, using this technique, distil qubit Werner states at a higher rate than could be achieved using previously discovered protocols. In particular, we demonstrate the advantage of our new distillation strategy by improving the best-known lower bound for the two-way-assisted quantum capacity of the qubit depolarising channel across a wide range of channel parameters, making progress in one of the long-standing problems of quantum information theory.

We prove that every PPT linear map has finite entanglement breaking index, thereby establishing the eventual entanglement breaking property of PPT channels in full generality. Furthermore, by utilizing completely positive maps with low entanglement dimensionality, we show that a large family of PPT maps, which strictly contains the class of 2-superpositive maps, has entanglement breaking index bounded above by 3, uniformly in the dimension. In particular, these results provide evidence that the PPT-cubed conjecture may hold in full generality.

Quasi-hydrodynamics describes systems with quasi-conserved degrees of freedom, namely observables that relax on timescales that are finite but parametrically longer than microscopic relaxation times. Examples include kinetic chemistry and linear viscoelasticity. Here, we develop a rigorous effective-field-theory framework for linear quasi-hydrodynamics from kinetic-type theories. Starting from any linearized, causal kinetic-like theory endowed with slow degrees of freedom, we show that the exact dynamics of conserved and quasi-conserved observables admits a systematic expansion in the fast relaxation timescale. At zeroth order, the resulting equations form a causal, symmetric-hyperbolic theory belonging to the appropriate transient-hydrodynamic universality class, establishing Israel-Stewart-like dynamics as the universal description of slow relaxation modes. Higher-order corrections can be computed systematically and inherit universal symmetry, Onsager, positivity, and causality constraints from the underlying microscopic theory.

We describe a mechanism of affine symmetry emergence in the maximal density regime of the prototypical model of self-organized criticality when the inverse square of the mesh of the underlying lattice is much larger than the number of random perturbation points distributed according to a prescribed probability measure supported in the interior of the ambient convex domain. Moreover, an appropriate scaling limit of the toppling function (aka odometer), which counts the number of operations per site, is a solution to a non-linear partial differential equation well known in the context of optimal transport and differential geometry, making it possible to accurately estimate the deviation of the density from its maximal value in any macroscopic window. The mechanism for the affine symmetry emergence is due to the novel empirical fact, supported in addition by inductive arguments that have recently being upgraded to a rigorous proof, that the scaling limit of the toppling function is the unique concave solution of the Monge-Ampère equation with Dirichlet boundary condition on the convex domain with the potential given by the probability measure used above as the infinite-perturbation profile.

Dyson universality typically manifests itself in local eigenvalue statistics. Here we show that its signature survives a nonlinear inverse-spectral reconstruction and reappears in the matrix geometry of the reconstructed operator. Using a dressing transformation, we map each unfolded spectrum to a deformation $f(x)$ of a fixed harmonic oscillator and represent it in the common oscillator basis by $F_{mn}=\bra m|f|n\ket$. We resolve the matrix-element weight into shells of fixed distance $d=|m-n|$, corresponding to the energy-transfer channels of the reference oscillator, and characterize the resulting distribution by distance-shell moments. Independently calibrated on Gaussian $\beta$-ensembles, these moments vary smoothly with $\beta$ and distinguish the GOE, GUE, and GSE. With this calibration fixed, applying the same diagnostic to the nontrivial zeros of the Riemann zeta function places the reconstructed operators in the GUE sector. Thus, the GUE character of the zeros is recovered not through direct statistics of the input levels, but from the distance-resolved geometry of the reconstructed operator.

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

It has been known for a long time that the equation of motion for dissipative linear dynamical systems with constant coefficients cannot be derived from the classical principle of stationary action because the term proportional to velocity in the equation of motion leads to the time derivative of order one half in the Lagrangian. Thus, approaches utilising fractional calculus have been used; however, they suffer from deficiencies both from mathematical and physical points of view. We here present our version of such fractional calculus based approach that provides correct Euler-Lagrange and, ultimately, the Hamilton equations, energy change of the moving body, and an attempt for a geometric interpretation of how energy dissipates.

This article contrasts the concepts of a proper condensate and the Onsager-Penrose criterion for Bose-Einstein condensation in the form of a debate between two proponents of the respective concepts. The prologue briefly introduces the two criteria. The epilogue contains remarks on properties of the underlying resolvent algebra that are physically relevant in this context.

Recently we started building a bridge between two cases of Langlands duality. The latter is one of the most influential topics in mathematical research. It has many different appearances and influential subtopics. Yet there is a topic that until now seems unrelated to the Langlands program. That is the topic of invariant differential operators. That is strange since both items are deeply rooted in Harish-Chandra's representation theory of semisimple Lie groups. We started with the case of the group ~$SL(2n)$. In the present paper we deal with the group $SL(2n+1)$. The two mentioned groups are similar, but their representation theories is rather different.

The existence of global classical solutions for reaction-diffusion systems arising from chemical reaction networks remains a major open problem in the deterministic setting, especially for reactions with high polynomial growth. We prove that a suitably chosen, physically motivated transport noise yields unique strong solutions for complex balanced chemical reaction networks that are global in time with arbitrarily high probability. These solutions possess paths in $C^{\theta}_t C^{\infty}_x$ for all $\theta<1/2$ and are, in particular, classical in space. Furthermore, we show that a suitable transport noise can enhance the dissipation of spatial fluctuations at an arbitrarily prescribed exponential rate. Our proofs rely on a combination of scaling-limit arguments, maximal $L^p(L^q)$-regularity, and entropy-entropy dissipation estimates.

We construct three derived endofunctors on the derived category of graded integrable representations of the current algebra of a flat degeneration of $D(2|1,\alpha)$, and prove that their classes in the Grothendieck group are the genus two Macdonald operators. Computing the graded $\operatorname{Ext}$ pairing explicitly, we deduce the self-adjointness of the genus two Macdonald operators and the orthogonality of genus $2$ Macdonald polynomials with respect to the certain skew-bilinear form.

The data processing inequality is a fundamental property that describes the loss of information through noisy channels. A more refined description is characterized by the strong data processing inequality (SDPI). In classical information theory, the tensorization of strong data processing inequality holds for a whole family of $f$-divergences. However, its quantum counterpart is less known. The tensorization of SDPI was shown only for some special cases previously, and the general understanding about the tensorization property of SDPI for quantum divergences remains open. In this work, we report two negative results: the tensorization property fails for certain quantum chi-square divergences, and it also does not hold for the quantum relative entropy.

We study Toeplitz determinants $\det T_n(e^f)$ for $f$ whose Fourier coefficients satisfy $f_k=O(|k|^{-1})$. This regime extends beyond $H^{1/2}$ and includes symbols with Fisher-Hartwig singularities. We develop an operator-theoretic approach based on the Baker-Campbell-Hausdorff formula that separates the quadratic term \[ \sum_{k=1}^{\infty}\min(k,n)f_kf_{-k} \] from the higher-order terms in the expansion of $\log\det T_n(e^{tf})$. We show that this quadratic term accounts for the possible growth with $n$, while every fixed higher-order coefficient remains bounded. For symbols with bounded positive and negative Fourier parts, our estimates yield two-sided bounds for the determinant after removal of the quadratic contribution. For a broader admissible class, including Fisher-Hartwig-type symbols, we obtain uniform higher-order coefficient bounds and a central limit theorem for the associated CUE linear statistics. We also obtain bounds on mixed exponential moments for CUE-derived random fields beyond the characteristic polynomial.

The most well-known beyond-Fourier heat conduction model, the Maxwell–Cattaneo–Vernotte equation is solved analytically in the presence of heat transfer boundary condition. In contrast to the corresponding Fourier problem, this boundary condition renders the underlying differential operator non-selfadjoint. With a suitable scalar product, the adjoint eigenvalue problem is established. The resulting left and right eigenfunctions constitute a biorthogonal system, allowing the expansion coefficients to be determined from arbitrary square-integrable initial conditions. This enables a convenient infinite-sum analytical solution, which is presented and thoroughly investigated for various values of the model parameters (including near-Fourier and highly hyperbolic regimes) and for two practically important initial conditions (equilibrium and flash pulse initiated). The spectral structure is analyzed in detail, including the occurrence of real, imaginary, and complex-conjugate eigenvalue roots, their asymptotic distribution, and their dependence on the dimensionless relaxation time and Biot number. We find good agreement with corresponding finite-difference numerical solutions. The comparison also illustrates the effects of spectral truncation (i.e., the Gibbs phenomenon) and numerical dissipation near propagating thermal-wave fronts. Completeness of the eigenfunction set is numerically demonstrated. The initial condition induced by the flash pulse is derived analytically. The results provide an analytical benchmark for non-selfadjoint hyperbolic heat conduction, unveil how boundary heat transfer and relaxation time influence the transition between Fourier-like diffusion and thermal-wave propagation, and open the possibility to find and investigate beyond-Fourier heat conduction via heat transfer in experiments and practical applications.

We prove uniqueness in the Lorentzian Calderón problem in a semiglobal setting, comparing a potentially large perturbation of the Minkowski metric against a small one. Our proof uses a reduced amount of data, solving a formally determined version of the problem. In contrast to traditional approaches, it employs neither control-theoretic arguments nor a reduction to a geometric inverse problem via microlocal methods or high-frequency solutions. Instead, it relies on distorted plane waves and weighted $L^2$-estimates.

The main goal of tomography is the reconstruction of density function out of its line integrals (integral measurements of this density along x-ray lines). Such construction is possible and known as inverse x-ray/Radon transformations. We are interested in generalization of this problem to the case of particles. For this we need to limit the speed of these particles, otherwise the problem is reduced to the previous one. The question we discuss in this paper is how low can this speed be depending on the parameters of the studied potential. We study this problem using simple theoretical examples and Machine Learning pipeline to restore the potential.

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

We study transparent boundary conditions (TBCs) for the time-dependent heat equation in a branching, quasi-one-dimensional domain modeled as a metric star graph. By combining the classical concept of TBCs with the theory of partial differential equations (PDEs) on networks, we derive an exact vertex condition that ensures unobstructed thermal flow across junctions. Specifically,we derive a sum rule for the diffusion coefficients that eliminates thermal backflow at the vertex.We demonstrate the validity of this analytical model numerically using a Crank-Nicolson finite-difference method, which confirms smooth, unobstructed heat propagation through the network.These results provide a practical mathematical framework for tunable control and optimization of thermal diffusion in low-dimensional structures. In particular, combining the well-known concept of the TBCs and theory for PDE on metric graphs, we propose a mathematical model providing a control tool for thermal diffusion in networks.

We generalize the classical vector-valued tangential and normal trace theory on Lipschitz domains to the setting of non-Lipschitz $H^1$-extension domains, which includes domains with fractal boundaries such as the Koch snowflake. We define and study these operators based on the surjectivity of the trace operator for elements of $H^1(\Omega)$ and the Hilbert structure of the associated trace space. The normal trace operator is defined on $\boldsymbol{H}({\rm div},\Omega)$ in $\mathbb{R}^n$. A generalized Stokes formula allows us to introduce the tangential trace operator on $\boldsymbol{H}({\bf curl},\Omega)$ and $\boldsymbol{H}^1(\Omega)$ in two and three dimensions. Following the approach of Buffa, Costabel and Sheen (2002) for Lipschitz domains, we define two abstract tangential boundary spaces as images of these trace operators, establish their Hilbert structure, and construct an abstract rotation operator linking them, in place of the geometric rotation based on the normal vector. We also extend Costabel's (1991) coercive bilinear form approach for Maxwell's equations to non-Lipschitz domains, yielding a theory that supports the treatment of the Hodge-Dirac operator and Green's formulas and enables the solution of boundary-value problems for the $\operatorname{\textbf{curl}}\operatorname{\textbf{curl}}+\textbf{I}$ operator within the $H^1$-extension domains framework.

We construct a discrete Smorodinsky–Winternitz II superintegrable system on a triangular region of the two-dimensional square lattice. The model is built from a pair of commuting finite-difference number operators with finite spectrum together with an associated ladder-operator structure. We show that it is maximally superintegrable and that its symmetry algebra admits a Hahn-algebra presentation. The spectral problem is solved exactly in terms of bivariate orthogonal polynomials of mixed Krawtchouk and dual Hahn type associated with the factorized $A_2$-Leonard pair. Finally, we show that the continuum limit recovers the continuous Smorodinsky–Winternitz II system together with its Hermite–Laguerre eigenfunctions. We further explain how the Hahn presentation of the discrete symmetry algebra becomes singular in this limit, while the limiting algebraic structure is naturally described by the Laguerre–Heun algebra associated with Cartesian and parabolic separation of variables.