How are Sheffer polynomials related to Lie theory? Sheffer polynomials $\{P_n(x)\}$ have generating function $P(x,t) = \sum_{n=0}^{\infty}P_n(x)t^n=A(t)e^{xu(t)}$.
This form reminds me of the Lie group–Lie algebra correspondence.  Is there any connection to Lie theory, obvious or not?
(I've deleted my reference to orthogonal polynomials, which was incorrect and irrelevant here.)
Here are some articles on Sheffer polynomials:

*

*Daniel Galiffa, Tanya Riston. An elementary approach to characterizing Sheffer A-type 0 orthogonal polynomial sequences. 2015.

*Dongsu Kim, Jiang Zeng. A Combinatorial Formula for the Linearization Coefficients of General Sheffer Polynomials. 2001.
 A: First, a general set of Sheffer polynomials is not orthogonal with respect to some weight function; for example, the prototypical sequence $p_n(x) = x^n$, which belongs to both special sub-groups of Sheffer polynomials—the Appell and binomial Sheffer polynomials—with the e.g.f. $e^{xt}$, is not an orthogonal set—neither is the celebrated Bernoulli Appell Sheffer sequence. (Usually the Sheffer polynomials are defined by $A(t) e^{xB(t)} = \sum_{n\geq 0} p_n(x) \frac{t^n}{n!}= e^{tp_\cdot(x)}$. The normalization $P_n(x)=n!p_n(x)$ often serves to give integer coefficients.)
However, there are connections to constructs associated to Lie theory:
1) The binomial Sheffer sequences have e.g.f.s of the form $e^{xh(t)}$, where $h(0)=0$ and $h'(0) \neq 0$. A dual sequence, its umbral compositional inverse, has the e.g.f. $e^{x\;h^{(-1)}(t)}$ defined by the compositional inverse function. As formulated by Charles Graves as early as 1853, with $g(z) = \frac{1}{\partial_z h(z)}$, the action on a function $f(z)$ analytic at $z$ of the exponential map of the infinitesimal generator (infinigen)
$ g(z)\frac{\partial }{\partial z} =  \frac{1}{h'(z)}\frac{\partial }{\partial z} =  \frac{\partial }{\partial h(z)} = \frac{\partial }{\partial \omega} $
is
(Eqn. 1)
$e^{t\;g(z)\frac{\partial }{\partial z}}\; f(z) = e^{t \;\frac{\partial }{\partial \omega}}\; f(h^{(-1)}(\omega))= f[h^{(-1)}(\omega + t)] =f[h^{(-1)}(h(z)+t)]$.
Then
$$e^{t\;g(z)\partial_z}\;z \rvert_{z=0}= h^{(-1)}(t),$$
so the e.g.f. for the associated binomial Sheffer sequence is
$$e^{t\;p.(x)}= e^{x\;h^{(-1)}(t)}= e^{t\;g(z)\partial_z} e^{zx} \;\rvert_{z=0}.$$
Edit Oct. 16, 2022: (Start)
Then
$$p_n(x) = (g(z)\partial_z)^n\; e^{zx} \;\rvert_{z=0}$$
$$= e^{-x} (g(y-1)\partial_y)^n e^{yx} \; |_{y=1}=  e^{-x} \left[x \;g\left(\frac{u}{x}-1\right)\partial_u\right]^n e^{u} \; |_{u=x}.$$
Specializing to $g(z)=(1+z)^{m+1}$, allows a general connection between the normal ordering of iterates of the Scherk-Witt-Lie vectors / infinigens $z^{m+1}\partial_z$ and well-known binomial Sheffer sequences. Then
$p_n^{(m)}(x) =((1+z)^{m+1}\partial_z)^{n} \;e^{zx}\;|_{z=0} = e^{-x}\; (y^{m+1}\partial_y)^n \;e^{yx} \; |_{y=1}$
$=  e^{-x}\; (x \;(\frac{u}{x})^{m+1}\partial_u)^n\; e^{u} \; |_{u=x} =  e^{-x}\; x^{-mn} \;(u^{m+1}\partial_u)^n\; e^{u} \; |_{u=x}$
$ =  e^{-x}\; x^{-mn} \;(x^{m+1}\partial_x)\;^n e^{x} =  e^{-x}\; p_n^{(m)}(:x\partial_x:)\; e^{x}$
with, by definition, $(:x\partial_x:)^n := x^n\partial_x^n$, implying
(Eqn. 2a)
$(x^{m+1} \partial_x)^n = x^{mn} p_n^{(m)}(:x\partial_x:) = x^{mn}\; m^n \;St1_n^{r}(\frac{x\partial_x}{m})$
where the e.g.f. for the $m$-th family of binomial Sheffer polynomials $p_n^{(m)}(x)$ is
(Eqn. 2b)
$e^{tp.^{(m)}(x)} = e^{h^{(-1)}(t)x}= \exp[((1-mt)^{-\frac{1}{m}}-1)x]$
and $St1_n^{r}(x)$ are the reversed unsigned Stirling polynomials of the first kind of OEIS A094638 (see also A008275, A048994, and A130534).
For the special linear Lie algebra $sl_2$:
For $m =-1$, the infingen $g(z)\partial_z = g(1+z)\partial_z= \partial_z$ has associated $h(z)=h^{(-1)}(z)=z$ giving the translation $e^{t\partial_z}\;f(z) = f(z+t)$, consistent with Eqn. (2) with $p^{(-1)}_n(z) = z^n$, the iconic Appell Sheffer sequence, which is also a binomial Sheffer sequence.
For $m=0$, the infingen $g(z)\partial_z = z\partial_z$ has associated $h(z)=\ln(z)$ and $h^{(-1)}(\omega)=e^\omega$ giving, from Eqn. (1), the scaling $e^{tz\partial_z}f(z) = f[\exp(\ln(z)+t)] = f(e^tz)$. Eqn. (2) associated with $g(z)=(1+z)$, $h(z) =\ln(1+z)$, and $h^{(-1)}(\omega) =e^{\omega}-1$ gives
$(z\partial_z)^n = p^{(0)}_n(:z\partial_z:) = St2_n(:z\partial_z:)$,
the Bell / Touchard / Scherk / Stirling polynomials of the second kind with the e.g.f. $\exp[(e^t-1)x]$ (cf. A008277 and A048993). This is easy to corroborate using $(z\partial_z)^n z^k = k^nz^k$.
For $m=1$, the infingen $g(z)\partial_z = z^2\partial_z$ has associated $h(x)=-\frac{1}{z}$ and $h^{(-1)}(\omega)=-\frac{1}{\omega}$ giving, from Eqn. (1), the vertical shearing $e^{tz^2\partial_z}f(z) = f\left [-\frac{1}{-\frac{1}{z}+t}\right] = f(\frac{z}{1-tz})$. Eqn. (2) associated with $g(z)=(1+z)^2$, $h(z) =\frac{z}{1+z}$, and $h^{(-1)}(\omega) =\frac{\omega}{1-\omega}$ gives
$(z^2\partial_z)^n =z^n p^{(1)}_n(:z\partial_z:) = z^n Lah_n(:z\partial_z:)$,
the shifted Lah polynomials, or shifted, normalized, unsigned Laguerre polynomials of order -1, $n! \; Lag_n^{(-1)}(-z)$, with the e.g.f. $\exp[z \frac{t}{1-t}]$ (cf. A008297 and A111596). This is corroborated by noting
$(z^2\partial_z)^n = z\;(z\partial_zz)^n\;z^{-1} = z\;z^n\partial_z^n z^n \;z^{-1} = z^n\; z \;(\partial_z^nz^n)\; z^{-1} = z^n\; n!\; z\; Lag_n(:z\partial_z:)\;z^{-1}$
from the Rodriguez formula for the associated Laguerre polynomial sequences.
(End)
2) Each Sheffer sequence has a pair of ladder ops—the raising/creation and lowering/destruction/annihilation ops defined by $L \; p_n(x) = n \; p_{n-1}(x)$ and $R \; p_n(x) = p_{n+1}(x)$—satisfying the Graves–Lie bracket of vector fields, the commutator, relation
$$[L,R] = 1,$$
from which the Graves–Pincherle derivative
$$[f(L),R] = f'(L)$$
follows.
There is also the commutator
$$[(g(z)\partial_z)^n,u(z)] = n \; (g(z)\partial_z)^{n-1}$$
for the powers of the Lie derivative / infinitesimal generator for the Scherk-Comtet partition polynomials of A139605, forming an underlying calculus for the binomial Sheffer sequences.
3) A general Sheffer sequence (a semidirect product of the Appell and binomial Sheffer sequences) is associated with a generalized Lie derivative via
$$e^{t(q(z)+g(z)\partial_z)}  \; e^{xz} |_{z=0} = A(t) e^{xf^{(-1)}(t)}$$
with $q(z) = \partial_t \ln(A(t)) \;|_{t=f(z)}$ and $g(z) = 1/f'(z)$.
I have an extensive set of posts at my web blog and links to numerous OEIS entries and MO-Q&As and MSE-Q&As related to this topic. All the symmetric polynomials/functions--complete, elementary, power, Faber--are related to Sheffer Appell sequences with their related ladder ops. The Faa di Bruno/Bell and cycle index polynomials of the symmetric groups, a.k.a the refined Stirling partition polynomials of the second and first kinds, and other compositional partition polynomials are all Appell Sheffer sequences in a distinguished indeterminate. The closely related sets of Lagrange inversion partition polynomials, including the refined Euler characteristic partition polynomials of the associahedra, all have the raising op $g(z)\partial_z$, as shown above--one set, OEIS A134264, is an Appell sequence and is related to free probability, inversion of Laurent series, and characterization of Kac-Schwarz operators related to Heisenberg-Virasoro groups. Multiplicative inversion is intimately bound with the refined Euler characteristic partition polynomials of the permutahedra and the operational (differential/matrix) calculus of Appell Sheffer polynomials. The Bernoulli Appell polynomials are of course related to the BCHD theorem, exponential mappings of the Lie commutator, topology, Todd operator and class of Hirzebruch, the Hurwitz zeta function, and more. Formal group laws (and generalized local Lie groups) are related to the linearization coefficients of products of binomial Sheffer polynomials. The list goes on.
The orthogonal Sheffer sequences such as the associated Laguerre polynomials have a  lot of underlying group theory and vector fields associated with them. Vilenkin with "Special Functions and the Theory of Group Representations", Talman, Miller, Gilmore, Feinsilver, and many others have written books on the underlying group theory. (None covers every aspect.)
One group of researchers has also written extensively on this topic. See, e.g., "One-parameter groups and combinatorial physics" by
Duchamp, Penson, Solomon, Horzela, and Blasiak (although quite negligent of reffing the OEIS--tribal instincts...sigh...What can you do?) as well as Wolfdieter Lang.
4) (Added Oct. 16, 2022) The raising operator of an Appell Sheffer sequence gives the infinitesimal generator, related to the Riemann zeta function, for Heaviside fractional differ-integral operators. See, e.g., my MO-Q "Riemann zeta function at positive integers and an Appell sequence of polynomials related to fractional calculus" and MO-A to "What's the matrix of logarithm of derivative operator (lnD)? What is the role of this operator in various math fields?".
