The para--Racah Hamiltonian as an exactly solvable discretization of the P\"oschl--Teller box

The para--Racah Hamiltonian is the Jacobi matrix $J$ of the para--Racah polynomials, taken at the persymmetric value of their isospectral deformation parameter: a mirror--symmetric tridiagonal matrix whose spectrum is a quadratic bi--lattice, the squares $(s+a{)}^2$ and $(s+c{)}^2$ of two interlaced arithmetic progressions of unit step, and which carries out fractional revival when the grid parameters $a$ and $c$ satisfy a Diophantine condition. As the number of sites grows, $J$ contracts to the trigonometric P\"oschl--Teller Hamiltonian on a finite interval. We establish this twice: on the recurrence relation, and on the polynomials themselves, for which it is the confluence Racah~$\to$~Jacobi. Both grid parameters survive: $c-a$ becomes a scalar centrifugal term at the interior point fixed by the mirror --- the exact analogue of the para--Krawtchouk origin --- while $a+c-1$ sets two confining singular walls at the ends, which turn the spectrum into the quadratic pair $(t+a{)}^2,(t+c{)}^2$ with Jacobi eigenfunctions $P_t^{(a+c-1,\mp (c-a))}$, the two sublattices converging to two families that differ only by the sign of the second Jacobi parameter. The mirror symmetry of $J$ becomes the reflection about the midpoint; it grades the two families and carries the fractional revival to the limit. Algebraically, the bispectral pair of the discrete model, which generates the Racah algebra ${\mathfrak{r}}_1$, becomes a pair of generators of the Jacobi algebra ${\mathfrak{j}}_1$ --- one rung down the Askey--Wilson hierarchy --- the dynamical algebra of the P\"oschl--Teller Hamiltonian.