Gelfand Triplets, Ladder Operators and Coherent States

1. Introduction and Motivation

As well known, observables in the usual formalism of Quantum Mechanics are given by self adjoint operators on a Hilbert space. Although the modern characterization of observable has shifted from the self-adjoint property to the PT-invariance property, we should keep Hamiltonians self adjoint for two reasons: i.) The evolution law for a given Hamiltonian H, $e^{-itH}$, is well defined and is unitary and, therefore, preserves the probabilities defined by wave functions and ii.) the current density is preserved if the Hamiltonian is self-adjoint.

However, self adjoint and the vast majority of quantum observables are given by unbounded self adjoint operators on a separable infinite dimensional Hilbert space. This destroys the intuition derived from matrix calculus to deal with operators. The use of unbounded operators involve many subtleties that do not appear in the matrix calculus or on its direct generalization, the calculus with bounded operators. The simple idea of unbounded operator implies the existence of its domain, a subspace of the Hilbert space on which the operator is well defined, which cannot be the whole Hilbert space. The situation is quite different in some other aspects. Let A be an operator with dense1 domain $\mathcal{D}\left(A\right)\subset \mathcal{H}$, where $\mathcal{H}$ is a separable infinite dimensional Hilbert space. Contrarily to what happens for bounded operators (i.e. matrices on a finite dimensional Hilbert space), the relation

$$\langle \psi |A\phi \rangle =\langle A\psi |\phi \rangle ,\forall \varphi ,\phi \in \mathcal{D}\left(A\right),$$

(1)

does not imply the self-adjointness of A. Just that its adjoint $A^{†}$ extends2A, $A\prec A^{†}$. Self-adjointness of A implies not only $A\prec A^{†}$, but also that $\mathcal{D}\left(A\right)\equiv \mathcal{D}\left(A^{†}\right)$, so that $A\equiv A^{†}$.

Bounded operators are continuous linear actions on the Hilbert space, while unbounded operators although linear are not continuous on Hilbert space. However, for an unbounded self adjoint operator A, there is always a dense subdomain thereof, $\mathrm{\Phi }$, endowed with a topology finer3 than the Hilbert space topology, on where A is continuous. This result is the well known Gelfand-Maurin theorem [1,2]. Note that this discussion can take place on infinite dimensional Hilbert spaces only.

The study of quantum one dimensional systems is quite important for two particular reasons. One is that one dimensional systems provides simple models with full quantum properties. In addition, these models are very often solvable or quasi solvable, which makes their study particularly interesting. Among these one dimensional systems there are some particularly interesting in which the Hamiltonian admits a factorization in terms of some operators called ladder operators (plus some additional term) [3,4,5,6,7]. This factorization may have an important consequence as it allows to construct a sequence of Hamiltonians each one with a similar spectrum than the precedent. Here, Hamiltonians and the ladder operators which factorize them are unbounded.

Then, our question is. Can we find a domain of the Hilbert space and a topology on this domain for which Hamiltonian and ladder operators which factorize it are continuous? This may go beyond a mathematical curiosity, as continuous linear operators have properties that may serve to cancel some of the complicated subtleties of non-continuous operators. The answer pass through a mathematically concept widely used in quantum physics, such that the notion of Gelfand triplet also known as rigged Hilbert space (RHS). Let us define this notion.

A Gelfand triplet or rigged Hilbert space is a tern of linear spaces over the complex field [8,9,10,11]:

$$\mathrm{\Phi }\subset \mathcal{H}\subset {\mathrm{\Phi }}^{\times },$$

(2)

where, i.) $\mathcal{H}$ is an infinite dimensional separable Hilbert space; ii.) $\mathrm{\Phi }$ is a dense subspace of $\mathcal{H}$ endowed with its own topology which is finer than the topology that $\mathrm{\Phi }$ has inherited from4$\mathcal{H}$; iii.) ${\mathrm{\Phi }}^{\times }$ is the space of continuous antilinear5 functionals6 on $\mathrm{\Phi }$ endowed with a topology compatible7 with the dual pair $\{\mathrm{\Phi },{\mathrm{\Phi }}^{\times }\}$. Often, ${\mathrm{\Phi }}^{\times }$ is called the antidual space of $\mathrm{\Phi }$.

The space $\mathrm{\Phi }$ is a locally convex space [12], with a topology defined by a set of seminorms8 of which one has to be the Hilbert space norm defined on $\mathrm{\Phi }$ as a subspace of $\mathcal{H}$. In our particular context, the topology on $\mathrm{\Phi }$ will be defined by a countably infinite set of seminorms, so that the space $\mathrm{\Phi }$ be metrizable9.

A sequence $\left\{{\varphi }_n\right\}\in \mathrm{\Phi }$ converges to $\varphi \in \mathrm{\Phi }$ if and only if for any seminorm p giving the topology on $\mathrm{\Phi }$, it happens that $p\left({\varphi }_n\right)⟼p\left(\varphi \right)$. If $\mathrm{\Phi }$ is metrizable, a linear10 mapping $F:\mathrm{\Phi }⟼\mathbb{C}$ is continuous if and only if for any convergent sequence ${\varphi }_n⟼\varphi$ in $\mathrm{\Phi }$, we have that $F\left({\varphi }_n\right)⟼F\left(\varphi \right)$. Similarly, an operator A on $\mathrm{\Phi }$ is continuous if $A{\varphi }_n⟼A\varphi$.

Nevertheless if $\mathrm{\Phi }$ were a locally convex space, other more operative method to check the continuity of linear or antilinear mappings is in order. If $F:\mathrm{\Phi }⟼\mathbb{C}$ is linear, then it is continuous if and only if there exists a positive constant $C>0$ and a finite number of seminorms $\{p_1,p_2,\cdots ,p_n\}$ from those which define the topology on $\mathrm{\Phi }$ such that for all $\varphi \in \mathrm{\Phi }$, we have:

$$|F\left(\varphi \right)|\le C\{p_1\left(\varphi \right)+p_2\left(\varphi \right)+\cdots +p_n\left(\varphi \right)\}.$$

(3)

Analogously, if $A:\mathrm{\Phi }⟼\mathrm{\Phi }$ is linear, it is continuous if and only if, for any seminorm $p_i$ defining the topology on $\mathrm{\Phi }$ there is a positive constant $C>0$ and k seminorms out of those which define the topology on $\mathrm{\Phi }$ such that for all $\varphi \in \mathrm{\Phi }$ [13],

$$p_i\left(A\varphi \right)\le C\{p_{i_1}\left(\varphi \right)+p_{i_2}\left(\varphi \right)+\cdots +p_{i_k}\left(\varphi \right)\}.$$

(4)

In principle, the constant C, the seminorms $p_{i_j}$ and its number, k, depend on the particular $p_i$ chosen. Obviously, C and the choice of the seminorms is not unique.

This result may be extended to any linear operator $A\mathrm{\Phi }⟼\mathrm{\Psi }$, where $\mathrm{\Phi }$ and $\mathrm{\Psi }$ are different locally convex spaces. If $\left\{p_i\right\}$ and $\left\{q_j\right\}$ are the families of seminorms that produce the topologies on $\mathrm{\Phi }$ and $\mathrm{\Psi }$, respectively, then (4) is now written as

$$q_j\left(A\varphi \right)\le C_j\{p_{i_{j_1}}\left(\varphi \right)+p_{i_{j_2}}\left(\varphi \right)+\cdots +p_{i_{j_k}}\left(\varphi \right)\},$$

(5)

for all seminorm $q_j$ on $\mathrm{\Psi }$, where $C_j$ and the seminorms $p_{i_{j_s}}$ depend on $q_j$.

At this point, let us define the notion of weak topology on the antidual space ${\mathrm{\Phi }}^{\times }$. This is a locally convex topology, for which the seminorms are defined like this: For each $\varphi \in \mathrm{\Phi }$, we define a seminorm $p_{\varphi }$ on ${\mathrm{\Phi }}^{\times }$, such that for any $F\in {\mathrm{\Phi }}^{\times }$, one has,

$$p_{\varphi }\left(F\right):=|F\left(\varphi \right)|.$$

(6)

Finally, it is convenient to use $\langle F|\varphi \rangle$ instead of $F\left(\varphi \right)$ to denote the action of $F\in {\mathrm{\Phi }}^{\times }$ on $\varphi \in \mathrm{\Phi }$. As a matter of fact, the weak topology11 induces on $\mathcal{H}$ a topology weaker (has less open sets) than the Hilbert space topology. Note that any $\psi \in \mathcal{H}$ gives an $F_{\psi }\in {\mathrm{\Phi }}^{\times }$ defined as

$$\langle F_{\psi }|\varphi \rangle =F_{\psi }\left(\varphi \right):=\langle \psi |\varphi \rangle ,\forall \varphi \in \mathrm{\Phi },$$

(7)

where $\langle -|-\rangle$ denotes the scalar product on $\mathcal{H}$. The functional $F_{\psi }$ is obviously linear on $\mathrm{\Phi }$. Furthermore, it is also continuous since

$$|\langle F_{\psi }|\varphi \rangle |=|\langle \psi |\varphi \rangle |\le ||\psi ||||\varphi ||=C{p}_0\left(\varphi \right),\forall \varphi \in \mathrm{\Phi },$$

(8)

with $C=||\psi ||$. Being given $\psi \in \mathcal{H}$, the $F_{\psi }\in {\mathrm{\Phi }}^{\times }$ is unique. This mapping $i\left(\psi \right)=F_{\psi }$, for all $\psi \in \mathcal{H}$ is antilinear and continuous12. Note that the choice of antilinear functionals instead of linear functionals in order to comply with the Dirac notation becomes obvious after (6).

Now, one understands the structure in (2). These spaces are related through continuous canonical injections (always $i\left(\psi \right)=\psi$, whatever $\psi$ is), $\mathrm{\Phi }⟼\mathcal{H}⟼{\mathrm{\Phi }}^{\times }$.

1.1. The Schwartz Space

Let us give an example of Gelfand triplet, which is more than a simple example, since it will help us to construct other triplets for our purposes. The Schwartz space $\mathcal{S}$ is the linear space of all functions $f\left(x\right):\mathbb{R}⟼\mathbb{C}$, where $\mathbb{R}$ is the real line, such that:

i.) Any $f\left(x\right)\in \mathcal{S}$ is indefinitely differentiable at all points of the real line $\mathbb{R}$.

ii.) Any $f\left(x\right)\in \mathcal{S}$ converges to zero at the infinite faster than the inverse of any polynomial. This means that for any (complex, x real) polynomial $P\left(x\right)$ and any $f\left(x\right)\in \mathcal{S}$, one has that

$$\underset{x\to \pm \infty }{lim}{\displaystyle \frac{f\left(x\right)}{P\left(x\right)}}=0.$$

(9)

All functions $f\left(x\right)\in \mathcal{S}$ are in $L^2\left(\mathbb{R}\right)$. Furthermore, $\mathcal{S}$ is dense in $L^2\left(\mathbb{R}\right)$ with the topology of the latter.

We may endow $\mathcal{S}$ with a locally convex metrizable topology into three equivalent forms, although we choose here one particularly interesting for our purposes. This comes after an interesting characterization of the Schwartz space $\mathcal{S}$ in terms of the Hermite functions, $\left\{{\psi }_n\left(x\right)\right\}$ to be defined in the next Section. Hermite functions form an orthonormal basis (complete orthonormal set) for the Hilbert space $L^2\left(\mathbb{R}\right)$, so that any $f\left(x\right)\in L^2\left(\mathbb{R}\right)$ takes the form

$$f\left(x\right)=\sum _{n=0}^{\infty }{a}_n{\psi }_n{\left(x\right),||f\left(x\right)||}^2=\sum _{n=0}^{\infty }{|a_n|}^2<\infty .$$

(10)

Then, $f\left(x\right)\in L^2\left(\mathbb{R}\right)$ is in $\mathcal{S}$ if and only if for the sequnce $\left\{a_n\right\}$ in (10), one has that for all $k=0,1,2,\cdots$,[13]

$$\sum _{n=0}^{\infty }{(n+1)}^{2k}{|a_n|}^2<\infty .$$

(11)

Then, let us define the following set of norms $p_k(-)=||-{||}_k$ for all $f\left(x\right)\in \mathcal{S}$:

$$p_k^2\left(f\right):={||f\left(x\right)||}_k^2:=\sum _{n=0}^{\infty }{(n+1)}^{2k}{|a_n|}^2,k=0,1,2,\cdots .$$

(12)

These norms (which are also seminorms) define the topology on $\mathcal{S}$. Note that i.) Since $k=0$ is a possible value, the list includes the norm on $L^2\left(\mathbb{R}\right)$. This implies that the topology on $\mathcal{S}$ is finer than the Hilbert space topology on $L^2\left(\mathbb{R}\right)$; ii.) Since the number of norms (seminorms) is countably infinite, the space $\mathcal{S}$ is metrizable. We may add that $\mathcal{S}$ is Frèchet space having the property of nuclearity13. Thus,

$$\mathcal{S}\subset L^2\left(\mathbb{R}\right)\subset {\mathcal{S}}^{\times }$$

(13)

is a Gelfand triplet or RHS. The antidual is isomorphic algebraic and topologically to the space of tempered distributions, usually defined as continuous linear functionals on $\mathcal{S}$.

This paper is organized as follows: In the next Section, we discuss the well known case of the Harmonic oscillator on the Schwartz space and the continuity of all operators involved on it, with also a mention to usual coherent states. Section 3 introduces the notion of Laguerre-Gauss special functions as orthonormal bases of spaces of type $L^2\left({\mathbb{R}}^{+}\right)$ with mention of an optical model which serves as a physical motivation for this mathematical construction. Functions of each orthonormal bases are characterized by a fixed value of a constant $l=0,1,2,\cdots$. For each one of the values of l, we define some ladder operators and give their intertwining relations with a Hamiltonian derived from the optical model. In Section 4, we introduce for each value of $l=0,1,2,\cdots$ a Gelfand triplet. Then, the Hamiltonian and ladder operators are continuous as linear mappings on these triplets. This is the main Section of the present article, in which we define new set of Gelfand triplets in order to fit as continuous linear functionals those Laguerre-Gauss functions with negative value of l. On Section 5, we give a brief presentation of various coherent states constructed on our spaces spanned by the Laguerre-Gauss functions. We finish with some concluding remarks and two appendices.

2. The Harmonic Oscillator

Although the contents of this Section are well known, we consider that a pedagogical account of them will help to a much better understanding of the motivation, purpose and proofs for the results contained in the main body of the present article. We begin with the ubiquitous Harmonic oscillator.

As is well known, the Hamiltonian of one dimensional quantum harmonic oscillator is given by

$$H={\displaystyle \frac{p^2}{2m}}+{\displaystyle \frac{1}{2}}m{\omega }^2{q}^2,\omega :=\sqrt{k/m}.$$

(14)

This Hamiltonian has a pure non-degenerate discrete spectrum with infinite values given by $E_n=\hslash \omega (n+1/2)$, $n=0,1,2,\cdots$. They respective eigenfunctions are the normalized Hermite functions:

$${\psi }_n\left(x\right):={\displaystyle \frac{1}{2^{n}n!}}\left({\displaystyle \frac{m\omega }{\pi \hslash }}\right)e^{-{\displaystyle \frac{m\omega x^2}{2\hslash }}}{H}_n\left(\sqrt{{\displaystyle \frac{m\omega }{\hslash }}}x\right),n=0,1,2,\cdots ,$$

(15)

where $H_n\left(y\right)$ are the Hermite polynomials. The annihilation and creation operators are, respectively, given by

$$a:=\sqrt{{\displaystyle \frac{m\omega }{2\hslash }}}\left(q+{\displaystyle \frac{i}{m\omega }}p\right),a^{†}:=\sqrt{{\displaystyle \frac{m\omega }{2\hslash }}}\left(q-{\displaystyle \frac{i}{m\omega }}p\right).$$

(16)

These operators are usually called, the ladder operators. Note that q and p are the multiplication and differentiation operators, respectively, i.e., $q\psi \left(x\right)=x\psi \left(x\right)$ and $p\psi \left(x\right)=-i\hslash {\psi }^{\prime }\left(x\right)$, where the prime denotes derivation with respect to x. In terms of the ladder operators, the Hamiltonian (14) is written as

$$H=\hslash \omega \left(N+{\displaystyle \frac{1}{2}}\right),N:=a^{†}a.$$

(17)

The operator N is the number operator. When using the ladder operators is somehow convenient a change in the notation for the sake of simplicity. Henceforth, we shall use $|n\rangle :={\psi }_n\left(x\right)$, all $n=0,1,2,\cdots$. The action of the ladder operators on the normalized Hermite function is given in this latter notation as

$$a|n\rangle =\sqrt{n}|n-1\rangle ,a|0\rangle =0,a^{†}|n\rangle =\sqrt{n+1}|n+1\rangle .$$

(18)

Some properties:

• The Hilbert space on which Hamiltonian and ladder operator act is $L^2\left(\mathbb{R}\right)$.
• The Hermite functions form an orthonormal basis (also called complete orthonormal set) in $L^2\left(\mathbb{R}\right)$. Then, the subspace of $L^2\left(\mathbb{R}\right)$ of (finite) linear combinations of Hermite functions is dense in $L^2\left(\mathbb{R}\right)$.
• Hermite functions are Schwartz functions.
• The Hamiltonian and the ladder operators are unbounded operators. Hence, they do not act on the whole $L^2\left(\mathbb{R}\right)$, but just on subspaces thereof called the domains of the operators. In any case Hermite functions lie in the domains of all these three operators, so that these domains are always dense in $L^2\left(\mathbb{R}\right)$. All these domains contain the Schwartz space as a subspace.

As shown in the Introduction, (13) is a Gelfand triplet. Let us show that H, a and $a^{†}$ are continuous operators on the Schwartz space $\mathcal{S}$. Let $\psi \left(x\right)\in \mathcal{S}$ and write

$$\psi \left(x\right)=\sum _{n=0}^{\infty }{a}_n{\psi }_n\left(x\right)=\sum _{n=0}^{\infty }{a}_n|n\rangle .$$

(19)

Since $\psi \left(x\right)\in \mathcal{S}$, the coefficients $\left\{a_n\right\}$ must satisfy (12). Let us define $a\psi$, the action of the operator a on $\psi \left(x\right)$ as

$$a\psi :=\sum _{n=1}^{\infty }{a}_n\sqrt{n}|n-1\rangle .$$

(20)

In a unique operation, we show that a as in (19) is well defined and is continuous on on $\mathcal{S}$. For any norm $p_k(-)$ as in (12), $k=0,1,2,\cdots$, we have for any $\psi \in \mathcal{S}$ that ($m=n+1$)

$$\begin{array}{c}\hfill p_k\left(a\psi \right)=\sum _{n=0}^{\infty }{b}_n^2{(n+1)}^{2k}=\sum _{n=0}^{\infty }|a_{n+1}{|}^2{(2n+1)}^{2k+1}=\sum _{m=1}^{\infty }{|a_m|}^2{m}^{2k+1}\\ \hfill \le \sum _{m=0}^{\infty }{|a_m|}^2{(m+1)}^{2k+2}=p_{k+1}\left(\psi \right),\end{array}$$

(21)

where the convergence of these sums for all k show that (20) is well defined. Now, let us define the action of $a^{†}$ on $\psi \left(x\right)\in \mathcal{S}$ as

$$a^{†}\psi :=\sum _{n=0}^{\infty }{a}_n\sqrt{n+1}|n+1\rangle .$$

(22)

Thus, for $k=0,1,2,\cdots$,

$$p_k\left(a^{†}\psi \right)=\sum _{n=0}^{\infty }{b}_n^2{(n+1)}^{2k}=\sum _{n=0}^{\infty }|a_n{|}^2{(n+1)}^{2k+1}\le \sum _{n=0}^{\infty }{|a_n|}^2{(n+1)}^{2k+2}=p_{k+1}\left(\psi \right),$$

(23)

which shows both that (22) is well defined and that $a^{†}$ is continuous on $\mathcal{S}$. Once we have shown the continuity of a and $a^{†}$, the continuity of H is obvious after (16). Nevertheless, this continuity may be proven directly just by noting that for any $k=0,1,2,\cdots$,

$$p_k\left(H\psi \right)={\hslash }^2{\omega }^2\sum _{n=0}^{\infty }{(n+1/2)}^2|a_n{|}^2{(n+1)}^{2k}\le {\hslash }^2{\omega }^2\sum _{n=0}^{\infty }{|a_n|}^2{(n+1)}^{2k+2}={\hslash }^2{\omega }^2{p}_{k+1}\left(\psi \right).$$

(24)

A note about the coherent states. Let $\alpha \in \mathbb{C}$ arbitrary but fixed. Its coherent state is given by

$$|\alpha \rangle :=e^{-{\displaystyle \frac{1}{2}}{|\alpha |}^2}\sum _{n=0}^{\infty }{\displaystyle \frac{{\alpha }^n}{\sqrt{n!}}}|n\rangle .$$

(25)

Coherent states are eigenvectors of the annihilation operator, so that $a|\alpha \rangle =\alpha |\alpha \rangle$. Coherent states evolve classically. In the present case, $|\alpha \rangle \in \mathcal{S}$ for all $\alpha \in \mathbb{C}$. To prove it, we just have to show that for all $k=0,1,2,\cdots$,

$$\sum _{n=0}^{\infty }{\displaystyle \frac{{|\alpha |}^{2n}}{n!}}{(n+1)}^{2k}<\infty .$$

(26)

Then, we just need to show that the general term in the series (26) goes to zero at the infinity faster than $1/n$. To see that this is true, just note that

$$\underset{n\to \infty }{lim}{n}^2{\displaystyle \frac{{|\alpha |}^{2n}}{n!}}{(n+1)}^{2k}=0,$$

(27)

which is a simple exercise. Thus coherent states for the Harmonic oscillator are Schwartz functions. In the sequel, we shall consider more general types of coherent states.

3. On Laguerre-Gaussian Ladder Operators

The Laguerre-Gauss functions have recently been considered by some authors in the study of two dimensional systems. As an example, they appear as the radial part of common eigenfunctions of the angular momentum and number operators for the two dimensional harmonic oscillator written in cylindric coordinates [15]. In the present paper, we are considering another type of model: the paraxial wave equation for parabolic media, given by the following three dimensional partial differential equation:

$$\left[-\left({\displaystyle \frac{{\partial }^2}{\partial x^2}}+{\displaystyle \frac{{\partial }^2}{\partial y^2}}\right)+{\displaystyle \frac{{\mathrm{\Omega }}^2}{2}}{\rho }^2\right]U=i{\displaystyle \frac{\partial }{\partial z}}U,{\mathrm{\Omega }}^2{\rho }^2<<1,$$

(28)

where $\mathrm{\Omega }$ is a constant and $\rho =\sqrt{x^2+y^2}$. This model describes the z-propagation of electromagnetic waves through media with square refractive index $1-{\mathrm{\Omega }}^2{\rho }^2$. This model has been studied in [16,17,18,19,20] and its physical motivation is not particularly relevant to our purposes.

Which really concerns to us is the form of (28). It is like a two dimensional time dependent Schrödinger equation for the harmonic oscillator, where the variable z plays the role of time. In this sense, it may be looked as a two dimensional model in spatial coordinates. The plane $(x,y)$ may be also described by its polar coordinates $(\rho ,\theta )$. Then, (28) admits as solutions functions of the type (see [16]):

$$U_{l,p}=C_{l,p}(z,\theta ){\mathrm{\Phi }}_{|l|,p}\left(\rho \right)l\in \mathbb{Z},p=0,1,2,\cdots .$$

(29)

Under some conditions [16], these solutions are square integrable, a fact that we are assuming from now on. We are not interested in the explicit form of the functions $C_{l,p}(z,\theta )$, our concern goes to ${\mathrm{\Phi }}_{|l|,p}\left(\rho \right)$, which are of the form [16]:

$${\mathrm{\Phi }}_{|l|,p}\left(\rho \right)={\displaystyle \frac{2{(-1)}^p}{w_0}}\sqrt{{\displaystyle \frac{\mathrm{\Gamma }(p+1)}{\mathrm{\Gamma }(|l|+p+1)}}}{\left[{\displaystyle \frac{\sqrt{2}\rho }{w_0}}\right]}^{|l|}exp\left(-{\displaystyle \frac{{\rho }^2}{w_0^2}}\right)L_p^{(|l|)}\left({\displaystyle \frac{2{\rho }^2}{w_0^2}}\right),$$

(30)

where $w_0$ is a constant and $L_p^{(|l|)}\left(x\right)$ are the associated Laguerre polynomials of degree p and order $|l|$. These special functions are the mentioned Laguerre-Gauss functions. The most important properties of these functions that concerns to us are [15,16]:

• For each fixed value of l, ${\mathrm{\Phi }}_{|l|,p}\left(\rho \right)$ form an orthonormal basis for $L^2({\mathbb{R}}^{+},\rho d\rho )$. Note that ordinary Laguerre functions

  $$M_n^{\alpha }\left(y\right)=\sqrt{{\displaystyle \frac{\mathrm{\Gamma }(n+1)}{\mathrm{\Gamma }(n+\alpha +1)}}}{y}^{\alpha /2}{e}^{-y/2}{L}_n^{\alpha }\left(y\right),$$

  (31)

  form an orthonormal basis for $L^2\left({\mathbb{R}}^{+}\right)$ for any fixed $\alpha \in (-1,\infty )$.
• The Laguerre-Gauss functions satisfy the following differential equation:

  $$\left[-{\displaystyle \frac{w_0^2}{2}}\left({\displaystyle \frac{{\partial }^2}{\partial {\rho }^2}}+{\displaystyle \frac{1}{\rho }}{\displaystyle \frac{\partial }{\partial \rho }}-{\displaystyle \frac{l^2}{{\rho }^2}}\right)+2{\displaystyle \frac{{\rho }^2}{w_0^2}}-2{\beta }_{|l|,p}\right]{\mathrm{\Phi }}_{|l|,p}\left(\rho \right)=0,$$

  (32)

  with

  $${\beta }_{|l|,p}=|l|+2p+1.$$

  (33)

Next, we define the following functions:

$${\psi }_{|l|,p}\left(\rho \right):={\rho }^{1/2}{\mathrm{\Phi }}_{|l|,p}\left(\rho \right).$$

(34)

After the aforementioned properties of ${\mathrm{\Phi }}_{|l|,p}\left(\rho \right)$, it is straightforward to show that, for each fixed value of l, the sequence $\left\{{\psi }_{|l|,p}\left(\rho \right)\right\}$, $p=0,1,2,\cdots$, is an orthonormal basis in $L^2\left({\mathbb{R}}^{+}\right)$. Functions ${\psi }_{|l|,p}\left(\rho \right)$ are also called Laguerre-Gauss.

Now for a fixed value of the integer, $l\in \mathbb{Z}$ and $\rho \ge 0$, let us consider the following operator on $L^2\left({\mathbb{R}}^{+}\right)$ [21]:

$$H_l:={\displaystyle \frac{w_0^2}{2}}\left[-{\displaystyle \frac{d^2}{d{\rho }^2}}+{\displaystyle \frac{l^2-\frac{1}{4}}{{\rho }^2}}\right]+2{\displaystyle \frac{{\rho }^2}{w_0^2}}.$$

(35)

After (32) and (34), we obviously have for $l\in \mathbb{Z}$ and $p=0,1,2,\cdots$,

$$H_l{\psi }_{|l|,p}\left(\rho \right)=2{\beta }_{|l|,p}{\psi }_{|l|,p}\left(\rho \right).$$

(36)

In the sequel, it is convenient to choose $w_0^2=2$ in (35) for simplicity in the notation. This will not change any of our relevant results. Since for fixed l, the sequence ${\psi }_{|l|,p}\left(\rho \right)$, $p=0,1,2,\cdots$, forms an orthonormal basis in $L^2\left({\mathbb{R}}^{+}\right)$, and these basis functions are eigenfunctions of $H_l$, being obviously $H_l$ symmetric, then, $H_l$ is essentially self adjoint on any subspace of $L^2\left({\mathbb{R}}^{+}\right)$ containing this basis for l. In particular, each of the $H_l$, $l\in \mathbb{Z}$, is essentially self adjoint on the subspaces ${\mathrm{\Phi }}_l$ to be defined later.

Then, (35) defines a family of Hamiltonians on $L^2\left({\mathbb{R}}^{+}\right)$ indexed by $\mathbb{Z}$ (or rather by $|l|=0,1,2,\cdots$, since $H_l$ depends on $l^2$). In any case, Hamiltonians of this family, $l\in \mathbb{Z}$, can be factorized in terms of the following operators [7]:

$$A_l^{-}:=\left({\displaystyle \frac{d}{d\rho }}-{\displaystyle \frac{l+\frac{1}{2}}{\rho }}+\rho \right),A_l^{+}:=\left(-{\displaystyle \frac{d}{d\rho }}-{\displaystyle \frac{l+\frac{1}{2}}{\rho }}+\rho \right),$$

(37)

and

$$B_l^{-}:=\left({\displaystyle \frac{d}{d\rho }}+{\displaystyle \frac{l-\frac{1}{2}}{\rho }}+\rho \right),B_l^{+}:=\left(-{\displaystyle \frac{d}{d\rho }}+{\displaystyle \frac{l-\frac{1}{2}}{\rho }}+\rho \right).$$

(38)

These objects as defined on (37) and (38) are the Laguerre-Gauss ladder operators. This factorization of $H_l$, $l\in \mathbb{Z}$, can be performed into four different modes and this goes as follows:

$$H_l=B_l^{+}{B}_l^{-}-{\epsilon }_{l-2}=B_{l+1}^{-}{B}_{l+1}^{+}-{\epsilon }_l=A_l^{+}{A}_l^{-}+{\epsilon }_l=A_{l-1}^{-}{A}_{l-1}^{+}+{\epsilon }_{l-2},$$

(39)

where $\left\{{\epsilon }_l\right\}$ is a sequence of real numbers. Some other relations are the following:

$$\begin{array}{c}\hfill A_l^{-}{H}_l=(H_{l+1}+2)A_l^{-},H_l{A}_l^{+}=A_l^{+}(H_{l+1}+2),\\ \hfill B_l^{-}{H}_l=(H_{l-1}+2)B_l^{-},H_l{B}_l^{+}=B_l^{+}(H_{l-1}+2).\end{array}$$

(40)

4. Gelfand Triplets and Continuous Operators

Let $\phi \left(\rho \right)\in L^2\left({\mathbb{R}}^{+}\right)$ arbitrarily chosen. Then, for fixed l, we know that

$$\phi \left(\rho \right)=\sum _{p=0}^{\infty }{\alpha }_{l,p}{\psi }_{|l|,p}\left(\rho \right),$$

(41)

where $\left\{{\alpha }_{l,p}\right\}$ is a sequence of complex numbers such that

$$\sum _{p=0}^{\infty }{|{\alpha }_{l,p}|}^2<\infty .$$

(42)

The convergence of the series in (41) is in the norm of $L^2\left({\mathbb{R}}^{+}\right)$, as is well known.

For fixed $l\ge 0$, let us construct now a Gelfand triplet in the following form: Let us choose the subspace ${\mathrm{\Phi }}_l$ of $L^2\left({\mathbb{R}}^{+}\right)$ of all vectors $\phi \left(p\right)$ as in (41) such that

$$p_{s,|l|}^2\left(\phi \right):=\sum _{p=0}^{\infty }|{\alpha }_{l,p}{|}^2{(|l|+2p+1)}^{2s}<\infty ,s=0,1,2,3,\cdots .$$

(43)

Then, $p_{s,|l|}$, $s=0,1,2\cdots$ is a family of seminorms14 that produce the topology on ${\mathrm{\Phi }}_l$. In Appendix B, we shall prove some topological properties of this space for $l>0$.

Note that ${\mathrm{\Phi }}_l$ contains all the functions on the orthonormal basis ${\left\{{\psi }_{|l|,p}\left(\rho \right)\right\}}_{p=0}^{\infty }$, so that it is dense in $L^2\left({\mathbb{R}}^{+}\right)$. Since the choice $s=0$ in (43) gives the norm on $L^2\left({\mathbb{R}}^{+}\right)$, the topology induced on ${\mathrm{\Phi }}_l$ by the Hilbert space norm is coarser than the topology on ${\mathrm{\Phi }}_l$. Hence the canonical injection $j:{\mathrm{\Phi }}_l⟼L^2\left({\mathbb{R}}^{+}\right)$, $j\left(\phi \right)=\phi$, $\forall \phi \in {\mathrm{\Phi }}_l$ is continuous. Therefore, if ${\mathrm{\Phi }}_l^{\times }$ is the antidual space of ${\mathrm{\Phi }}_l$, we conclude that

$${\mathrm{\Phi }}_l\subset L^2\left({\mathbb{R}}^{+}\right)\subset {\mathrm{\Phi }}_l^{\times },$$

(44)

is a Gelfand triplet or rigged Hilbert space (RHS in the sequel).

Let us prove that the Hamiltonian $H_l$ leaves invariant ${\mathrm{\Phi }}_l$, which means that $H_l\psi \in {\mathrm{\Phi }}_l$ for any $\phi \in {\mathrm{\Phi }}_l$, and that it is a continuous operator on ${\mathrm{\Phi }}_l$. First of all, we need to define the action of $H_l$ on each of the $\phi \in {\mathrm{\Phi }}_l$, with the requirement that it extends (36). This definition should go as

$$H_l\sum _{p=0}^{\infty }{\alpha }_{l,p}{\psi }_{|l|,p}\left(\rho \right)=\sum _{p=0}^{\infty }{\alpha }_{l,p}(|l|+2p+1){\psi }_{|l|,p}\left(\rho \right).$$

(45)

To see that (45) is well defined, note that

$$p_{s,|l|}^2\left(H_l\phi \right)=\sum _{p=0}^{\infty }|{\alpha }_{l,p}{|}^2{(|l|+2p+1)}^{2(s+1)}=p_{s+1}^2\left(\phi \right).$$

(46)

This shows that the coefficients ${\alpha }_{l,p}(|l|+2p+1)$ of $H_l\phi$ verify relations (43) if $\phi \in {\mathrm{\Phi }}_l$, so that $H_l\phi \in {\mathrm{\Phi }}_l$ and $H_l$ leaves ${\mathrm{\Phi }}_l$ invariant for all $l\in \mathbb{Z}$. In addition, after (8) one concludes that $H_l$ is continuous on ${\mathrm{\Phi }}_l$. Note that none of the $H_l$ is a bounded (continuous) operator on $L^2\left({\mathbb{R}}^{+}\right)$.

4.1. The Ladder Operators

Although the index l could, in principle, run out the set of entire numbers, the Laguerre-Gauss basis functions (30) and (34) are just defined for $l\ge 0$. Laguerre-Gauss functions for negative values of l can be defined, although are highly singular at the origin. This means that such functions cannot be square integrable for the Lebesgue measure on the positive semiaxis, although some other measures of the type ${\rho }^{|l|}d\rho$ could be used.

The action of the ladder operators (37) and (38) on the basis functions ${\psi }_{|l|,p}\left(\rho \right)$ has been well studied [16,22]. For $l\ge 0$, we have that

$$\begin{array}{c}\hfill A_l^{-}{\psi }_{l,p}\left(\rho \right)=2\sqrt{p}{\psi }_{l+1,p-1}\left(\rho \right),A_{l-1}^{+}{\psi }_{l,p}\left(\rho \right)=2\sqrt{p+1}{\psi }_{l-1,p+1}\left(\rho \right)(l\ne 0)\\ \hfill B_l^{-}{\psi }_{l,p}\left(\rho \right)=2\sqrt{p+l}{\psi }_{l-1,p-1}\left(\rho \right)(l\ne 0),B_{l+1}^{+}{\psi }_{l,p}\left(\rho \right)=2\sqrt{p+l+1}{\psi }_{l+1,p}\left(\rho \right).\end{array}$$

(47)

Functions with subindex $p-1$ are identically equal to zero if $p=0$.

Next, we have the following result: Ladder operators are linear continuous mappings between the following respective spaces:

$$\begin{array}{ccc}\hfill A_l^{-}& :& {\mathrm{\Phi }}_l⟼{\mathrm{\Phi }}_{l+1},\hfill \end{array}$$

(48)

$$\begin{array}{ccc}\hfill A_{l-1}^{+}& :& {\mathrm{\Phi }}_l⟼{\mathrm{\Phi }}_{l-1},(l\ne 0),\hfill \end{array}$$

(49)

$$\begin{array}{ccc}\hfill B_l^{-}& :& {\mathrm{\Phi }}_l⟼{\mathrm{\Phi }}_{l-1},(l\ne 0),\hfill \end{array}$$

(50)

$$\begin{array}{ccc}\hfill B_{l+1}^{+}& :& {\mathrm{\Phi }}_l⟼{\mathrm{\Phi }}_{l+1}.\hfill \end{array}$$

(51)

The proof is simple. First, take $\phi \in {\mathrm{\Phi }}_l$, so that taking into account (41), we have

$$\left(A_l^{-}\psi \right)\left(\rho \right)=A_l^{-}\sum _{p=0}^{\infty }{\alpha }_{l,p}{\psi }_{l,p}\left(\rho \right)=\sum _{p=0}^{\infty }{\alpha }_{l,p}2\sqrt{p}{\psi }_{l+1,p-1}\left(\rho \right).$$

(52)

Then,

$$p_{s,|l+1|}^2\left(A_l^{-}\phi \right)=4\sum _{p=0}^{\infty }p|{\alpha }_{l,p}{|}^2{(|l+1|+2p+1)}^{2s},$$

(53)

with $s=0,1,2,\cdots$. Clearly, (henceforth, we omit the modulus for l since we are assuming that $l\ge 0$)

$$p\le (l+1+2p+1),$$

(54)

and

$$\begin{array}{c}\hfill l+1+2p+1=l+2p+1+1\le 2(l+2p+1),\end{array}$$

(55)

so that, (53) becomes

$$p_{s,l+1}^2\left(A_l^{-}\phi \right)\le 8\sum _{p=0}^{\infty }{|{\alpha }_{l,p}|}^2{(l+2p+1)}^{2(s+1)}=8p_{s+1,l}^2\left(\phi \right),\forall \phi \in {\mathrm{\Phi }}_l.$$

(56)

Equation (56) shows that

1.- The series in (53) converges for all values of $s=0,1,2,\cdots$ for $l+1$, which implies that $A_l^{-}\phi$ belongs to ${\mathrm{\Phi }}_{l+1}$ for all $\phi \in {\mathrm{\Phi }}_l$.

2.- In accordance to (5), the mapping $A_l^{-}:{\mathrm{\Phi }}_l⟼{\mathrm{\Phi }}_{l+1}$ is continuous with the topologies defined on both spaces.

The definition for the remainder ladder operators on ${\mathrm{\Phi }}_l$ is similar. In fact,

$$A_{l-1}^{+}\phi :=2\sum _{p=1}^{\infty }\sqrt{p+1}{\alpha }_{l,p}{\psi }_{l-1,p+1}\left(\rho \right),$$

(57)

$$B_l^{-}\phi :=2\sum _{p=1}^{\infty }\sqrt{p+l}{\alpha }_{l,p}{\psi }_{l-1,p-1}\left(\rho \right),$$

(58)

and

$$B_{l+1}^{+}\phi :=2\sum _{p=0}^{\infty }\sqrt{p+l+1}{\alpha }_{l,p}{\psi }_{l+1,p}\left(\rho \right).$$

(59)

The proof of the continuity of these operators is made in analogy with the proof we have given for $A_l^{-}$. Take for instance,

$$\begin{array}{c}\hfill p_{s,l+1}^2\left(B_{l+1}^{+}\phi \right)=4\sum _{p=0}^{\infty }(p+l+1){|{\alpha }_{l,p}|}^2{(l+1+2p+1)}^{2s}\\ \hfill \le 8\sum _{p=0}^{\infty }|{\alpha }_{l,p}{|}^2{(l+2p+1)}^{2s+3}\le 8\sum _{p=0}^{\infty }{|{\alpha }_{l,p}|}^2{(l+2p+1)}^{2(s+2)}=8p_{s,l}^2\left(\phi \right),\forall \phi \in {\mathrm{\Phi }}_l.\end{array}$$

(60)

This proves, first that $B_{l+1}^{+}$ maps ${\mathrm{\Phi }}_l$ into ${\mathrm{\Phi }}_{l+1}$ and, then, that this map is continuous with respect to the topologies of the initial and final spaces. Same kind of arguments goes for $A_{l-1}^{+}$ and $B_l^{-}$.

4.2. The Laguerre-Gauss Functions with Negative Value of l

Let us go back to the family of functions in (30), were for simplicity we choose $w_0^2=2$. This choice simplifies the notation and does not affect to the mathematical discussion under process. Then, the functions in (30) for fixed $l\ge 0$ and $p=0,1,2,\cdots$ form an orthonormal basis for $L^2({\mathbb{R}}^{+},\rho d\rho )$. Then, for any $\psi \left(\rho \right)\in L^2({\mathbb{R}}^{+},\rho d\rho )$, we have the following span:

$$\psi \left(\rho \right)=\sum _{p=0}^{\infty }{a}_{l,p}{\mathrm{\Phi }}_{l,p}\left(\rho \right),\sum _{p=0}^{\infty }{|a_{l,p}|}^2<\infty ,$$

(61)

where the first sum in (61) converges in the norm of $L^2({\mathbb{R}}^{+},\rho d\rho )$.

Now, let ${\mathrm{\Psi }}_l$, $l\ge 0$, the space of all $\psi \left(\rho \right)\in L^2({\mathbb{R}}^{+},\rho d\rho )$, such that for fixed $l\ge 0$, we have that

$$q_{s,l}^2\left(\psi \right):=\sum _{p=0}^{\infty }{|a_{l,p}|}^2{[(l+2p+1)!]}^{2s}<\infty ,s=0,1,2,\cdots .$$

(62)

No matter how small the spaces ${\mathrm{\Psi }}_l$ could be as subspaces of $L^2({\mathbb{R}}^{+},\rho d\rho )$, they are dense in $L^2({\mathbb{R}}^{+},\rho d\rho )$ with the Hilbert space topology, since all ${\mathrm{\Psi }}_l$ contains a orthonormal basis of $L^2({\mathbb{R}}^{+},\rho d\rho )$. In addition, the seminorm in (62) $q_{l,0}$ is nothing else than the Hilbert space norm, so that the canonical injection $i:{\mathrm{\Psi }}_l⟼L^2({\mathbb{R}}^{+},\rho d\rho )$, $i\left(\psi \right)=\psi$ is continuous for all values of $l=0,1,2,\cdots$. We have a new series of Gelfand triplets given by

$${\mathrm{\Psi }}_l\subset L^2({\mathbb{R}}^{+},\rho d\rho )\subset {\mathrm{\Psi }}_l^{\times },l=0,1,2,\cdots .$$

(63)

Let us go back to the beginning of Section 4. Also for $l\ge 0$ and a subspace, ${\tilde{\mathrm{\Phi }}}_l$, of vectors $\psi \left(\rho \right)\in L^2({\mathbb{R}}^{+},\rho d\rho )$, we may define the set of norms (43), wich for an “abuse de langage” shall denote as $p_{s,l}\left(\psi \right)$, exactly as in the case described there15. This gives a new Gelfand triplet:

$${\tilde{\mathrm{\Phi }}}_l\subset L^2({\mathbb{R}}^{+},\rho d\rho )\subset {\tilde{\mathrm{\Phi }}}_l^{\times }.$$

(64)

Furthermore, it is obvious that for $l=0,1,2,\cdots$, ${\mathrm{\Psi }}_l\subset {\tilde{\mathrm{\Phi }}}_l$. Let us call $j_l:{\mathrm{\Psi }}_l⟼{\tilde{\mathrm{\Phi }}}_l$ to the corresponding canonical injection. Then for any $\psi \in {\mathrm{\Psi }}_l$ as in (61), we have that

$$p_{s,l}\left(j_l\left(\psi \right)\right)=\sum _{p=0}^{\infty }|a_{l,p}{|}^2{(l+2p+1)}^{2s}\le \sum _{p=0}^{\infty }{|a_{l,p}|}^2{[(l+2p+1)!]}^{2s}=q_{s,l}\left(\psi \right),$$

(65)

so that $j_l$ is continuous. As a consequence, we have the following sequence of locally convex spaces for which all the canonical mappings are continuous:

$${\mathrm{\Psi }}_l\subset {\tilde{\mathrm{\Phi }}}_l\subset L^2({\mathbb{R}}^{+},\rho d\rho )\subset {\tilde{\mathrm{\Phi }}}_l^{\times }\subset {\mathrm{\Psi }}_l^{\times },l=0,1,2,\cdots .$$

(66)

The reason to adopt such an strong topology will be evident later.

Next, let us take the “l negative version of the functions (30)”, which is ($l\ge 0$)

$${\mathrm{\Phi }}_{-l,p}\left(\rho \right)={(-1)}^p\sqrt{{\displaystyle \frac{\mathrm{\Gamma }(p+1)}{\mathrm{\Gamma }(l+p+1)}}}{\rho }^{-l}exp(-{\rho }^2/2)L_p^l\left({\rho }^2\right).$$

(67)

Note that the associate Laguerre polynomial $L_m^{\alpha }\left(x\right)$ makes sense only if $\alpha \in (-1,\infty )$, so that we cannot use $\alpha =-l$. The objective is to show that ${\mathrm{\Phi }}_{-l,p}\left(\rho \right)\in {\mathrm{\Psi }}_l^{\times }$, $l=0,1,2,\cdots$.

First of all, let us prove that if $\psi \left(\rho \right)\in {\mathrm{\Psi }}_l$, then, the first series in (61) converges also uniformly. To this end, let us consider the following relation concerning Laguerre functions:

$$e^{-x}{x}^{k/2}{L}_n^k\left(x\right)={\displaystyle \frac{1}{n!}}{\int }_0^{\infty }{e}^{-t}{t}^{n+k/2}{J}_k\left(2\sqrt{tx}\right)dt,$$

(68)

where $J_k\left(t\right)$ is the Bessel function with index k. If k is an integer number, which is the case in our discussion, the Bessel function $J_k\left(t\right)$ is given by

$$J_k\left(y\right)={\displaystyle \frac{1}{\pi }}{\int }_0^{\pi }cos(k\tau -ysin\tau )d\tau ,$$

(69)

so that $|J_k\left(y\right)|\le 1$, k being an integer number. For us, $k=l$ and $n=p$. Then,

$$|e^{-x}{x}^{l/2}{L}_p^l\left(x\right)|\le {\displaystyle \frac{1}{p!}}{\int }_0^{\infty }{e}^{-t}{t}^{p+l/2}dt={\displaystyle \frac{\mathrm{\Gamma }(p+l/2+1)}{p!}}.$$

(70)

Taking $x={\rho }^2$, we have that

$$|e^{-{\rho }^2}{\rho }^l{L}_p^l\left({\rho }^2\right)|\le {\displaystyle \frac{\mathrm{\Gamma }(p+l/2+1)}{p!}}.$$

(71)

If $\psi \left(\rho \right)\in {\mathrm{\Psi }}_l$ as in (61), one has the following series of inequalities for each $l=0,1,2,\cdots$:

$$\begin{array}{c}\hfill ||\psi \left(\rho \right)||=\left|\sum _{p=0}^{\infty }{a}_{l,p}{\mathrm{\Phi }}_{l,p}\left(\rho \right)\right|\le \sum _{p=0}^{\infty }|a_{l,p}||{\mathrm{\Phi }}_{l,p}\left(\rho \right)|\le \sum _{p=0}^{\infty }|a_{l,p}|\sqrt{{\displaystyle \frac{\mathrm{\Gamma }(p+1)}{\mathrm{\Gamma }(l+p+1)}}}{\displaystyle \frac{\mathrm{\Gamma }(p+l/2+1)}{p!}}\\ \hfill \le \sum _{p=0}^{\infty }|a_{l,p}|(l+2p+1)!\le \sum _{p=0}^{\infty }|a_{l,p}|{[(l+2p+1)!]}^2{\displaystyle \frac{1}{(l+2p+1)!}}\\ \hfill \le \sqrt{{\sum }_{p=0}^{\infty }{|a_{l,p}|}^2{[(l+2p+1)!]}^4}\times \sqrt{{\sum }_{p=0}^{\infty }{\displaystyle \frac{1}{{[(l+2p+1)!]}^2}}}=C_l\left[q_{2,l}\left(\psi \right)\right],\end{array}$$

(72)

where $C_l:=\sqrt{{\sum }_{p=0}^{\infty }{\displaystyle \frac{1}{{[(l+p+1)!]}^2}}}$. The last inequality in (72) is the Schwartz inequality. Due to the Weierstrass M-test, the series ${\sum }_{p=0}^{\infty }|a_p||{\mathrm{\Phi }}_{l,p}\left(\rho \right)|$ converges uniformly for $\rho \in [0,\infty )$16.

Now, let $\psi \left(\rho \right)$ be an arbitrary function in ${\mathrm{\Psi }}_l$. The action of ${\mathrm{\Phi }}_{-l,q}\left(\rho \right)$ as in (66), $q=0,1,2,\cdots$, as a functional on ${\mathrm{\Psi }}_l$ should be for all $\psi \left(\rho \right)\in {\mathrm{\Psi }}_l$,

$$\begin{array}{c}\hfill {\int }_0^{\infty }\psi \left(\rho \right){\mathrm{\Phi }}_{-l,q}\left(\rho \right)\rho d\rho ={\int }_0^{\infty }\left(\sum _{p=0}^{\infty }{a}_{l,p}{\mathrm{\Phi }}_{l,p}\left(\rho \right)\right){\mathrm{\Phi }}_{-l,q}\left(\rho \right)\rho d\rho \\ \hfill =\sum _{p=0}^{\infty }{a}_{l,p}{\int }_0^{\infty }{\mathrm{\Phi }}_{l,p}\left(\rho \right){\mathrm{\Phi }}_{-l,q}\left(\rho \right)\rho d\rho ,\end{array}$$

(73)

where the last identity in (73) is due to the uniform convergence of the series. The last integral in (73) gives

$${\int }_0^{\infty }{\mathrm{\Phi }}_{l,p}\left(\rho \right){\mathrm{\Phi }}_{-l,q}\left(\rho \right)\rho d\rho =\sqrt{{\displaystyle \frac{\mathrm{\Gamma }(p+1)\mathrm{\Gamma }(q+1)}{\mathrm{\Gamma }(l+p+1)\mathrm{\Gamma }(l+q+1)}}}{\int }_0^{\infty }{e}^{-{\rho }^2}{L}_p^l\left({\rho }^2\right)L_q^l\left({\rho }^2\right)\rho d\rho .$$

(74)

Next, let us use the explicit form for the associated Laguerre polynomials

$$L_p^l\left({\rho }^2\right)={\displaystyle \frac{1}{p!}}\sum _{s=0}^p{\displaystyle \frac{p!}{s!}}\left(\begin{array}{c}p+l\\ p-s\end{array}\right){(-{\rho }^2)}^s,$$

(75)

so that

$$\begin{array}{c}\hfill \left|{\int }_0^{\infty }{e}^{-{\rho }^2}{L}_p^l\left({\rho }^2\right)L_q^l\left({\rho }^2\right)\rho d\rho \right|\\ \hfill \le \sum _{s=0}^p{\displaystyle \frac{1}{s!}}{\displaystyle \frac{(p+l)!}{(p-s)!(l+s)!}}\sum _{r=0}^{p^{\prime }}{\displaystyle \frac{1}{r!}}{\displaystyle \frac{(p^{\prime }+l)!}{(p^{\prime }-r)!(l+r)!}}{\int }_0^{\infty }{e}^{-{\rho }^2}{\rho }^{2(r+s)}\rho d\rho \end{array}$$

(76)

Let us make the change of variables ${\rho }^2=t$ in the integral. It gives

$${\displaystyle \frac{1}{2}}{\int }_0^{\infty }{e}^{-t}{t}^{r+s}dt={\displaystyle \frac{1}{2}}\mathrm{\Gamma }(r+s+1)={\displaystyle \frac{1}{2}}(r+s)!,$$

(77)

so that (75) is smaller or equal to

$$\sum _{s=0}^p\sum _{r=0}^{p^{\prime }}(p+l)!(p^{\prime }+l)!(s+r)!.$$

(78)

Now, we have two possibilities. Either $p^{\prime }>p$, so that $(s+r)!<\left(2p^{\prime }\right)!$ or $p^{\prime }\le p$, in which case $(s+r)!\le \left(2p\right)!$. Also note that $rs\le (r+s)!$, so that if $p^{\prime }>p$, then (78) is smaller or equal to

$$(p+l)!(p^{\prime }+l)!\left(2p^{\prime }\right)!\left(rs\right)\le (p+l)!\{(p^{\prime }+l)!{[\left(2p^{\prime }\right)!]}^2\}\le (l+2p+1)!\{(p^{\prime }+l)!{[\left(2p^{\prime }\right)!]}^2\}.$$

(79)

Note that the term between the brackets does not depend on p and only of l and $p^{\prime }$ which are fixed here. Then, if $F_{p^{\prime }}$ represents the functional on ${\mathrm{\Psi }}_l$ defined by ${\mathrm{\Phi }}_{-l,p^{\prime }}\left(\rho \right)$ with (73), we have for all $\psi \in {\mathrm{\Psi }}_l$:

$$\begin{array}{c}\hfill |F_q\left(\psi \right)|=\left|{\int }_0^{\infty }\psi \left(\rho \right){\mathrm{\Phi }}_{-l,q}\left(\rho \right)\rho d\rho \right|\le \sum _{p=0}^{\infty }|a_{l,p}|q!\left|{\int }_0^{\infty }{e}^{-{\rho }^2}{L}_p^l\left({\rho }^2\right)L_q^l\left({\rho }^2\right)\rho d\rho \right|\\ \hfill \le C_{q,l}\sum _{p=0}^{q-1}|a_{l,p}|(l+2p+1)!+\sum _{p=q}^{\infty }|a_{l,p}|q!\left|{\int }_0^{\infty }{e}^{-{\rho }^2}{L}_p^l\left({\rho }^2\right)L_q^l\left({\rho }^2\right)\rho d\rho \right|,\end{array}$$

(80)

with $C_{q,l}:=q!(q+l)!{[\left(2q\right)!]}^2$. We analyze both summands in the last row of (80) separately. For the first one, we have:

$$\begin{array}{c}\hfill C_{q,l}\sum _{p=0}^{q-1}|a_{l,p}|(l+2p+1)!=C_{p^{\prime },l}\sum _{p=0}^{q-1}|a_{l,p}|{((l+2p+1)!)}^2{\displaystyle \frac{1}{(l+2p+1)!}}\\ \hfill \le C_{q,l}\sqrt{{\sum }_{p=0}^{q-1}|a_{l,p}{{|}^2(l+2p+1)!)}^4}\sqrt{{\sum }_{p=0}^{q-1}{\displaystyle \frac{1}{{((l+2p+1)!)}^2}}}\\ \hfill \le C_{q,l}\sqrt{{\sum }_{p=0}^{\infty }|a_{l,p}{{|}^2(l+2p+1)!)}^4}\sqrt{{\sum }_{p=0}^{\infty }{\displaystyle \frac{1}{{((l+2p+1)!)}^2}}}=C_{q,l}{C}_l\left[q_{2,l}\left(\psi \right)\right],\end{array}$$

(81)

where the second inequality is nothing else that the Cauchy-Schwarz inequality and the meaning of $C_l$ is obvious. For the second term in the last row of (80) ($p\ge p^{\prime }$), we just have to replace ${[\left(2p^{\prime }\right)!]}^2$ in (80) by ${[\left(2p\right)!]}^2$, which is smaller or equal to ${[(l+2p+1)!]}^2$ and $\left(p^{\prime }\right)!$ by $p!\le (l+2p+1)!$. Then if $K_{p^{\prime },l}:=(p^{\prime }+l)!$ and following similar arguments as just before, we have that

$$\begin{array}{c}\hfill \sum _{p=q}^{\infty }|a_{l,p}|q!\left|{\int }_0^{\infty }{e}^{-{\rho }^2}{L}_p^l\left({\rho }^2\right)L_q^l\left({\rho }^2\right)\rho d\rho \right|\le K_{q,l}\sum _{p=q}^{\infty }|a_{l,p}|{[(l+2p+1)!]}^4\\ \hfill \le K_{q,l}\sum _{p=0}^{\infty }|a_{l,p}|{[(l+2p+1)!]}^4\le \sqrt{{\sum }_{p=0}^{\infty }{|a_{l,p}|}^2{[(l+2p+1)!]}^{10}}\sqrt{{\displaystyle \frac{1}{{[(l+2p+1)!]}^2}}}\\ \hfill \le {\tilde{K}}_{q,l}\left[q_{4,l}\left(\psi \right)\right],\end{array}$$

(82)

where ${\tilde{K}}_{q,l}$ is a constant. Let $K:=max\{C_{q,l}{C}_l,K_{q,l}\}$ and take into account that for all $\psi \in {\mathrm{\Psi }}_l$, $q_{2,l}\left(\psi \right)\le q_{4,l}\left(\psi \right)$, we finally arrive to

$$|F_q\left(\psi \right)|\le 2K\left[q_{4,l}\left(\psi \right)\right],\forall \psi \in {\mathrm{\Psi }}_l,$$

(83)

which proves the continuity of $F_q$ on ${\mathrm{\Psi }}_l$. This goes for all $q=0,1,2,\cdots$, which ends the proof of our statement.

The need of introducing such an strong topology is due to some difficulties producing by the terms in the first sum of (76) and which cannot be simplified so that the use of the topology (43) can be feasible. Since this new topology given by the norms (62) is stronger than (43), it would be interesting to check whether the properties of continuity of creation and annihilation operators are preserved under the new topology.

4.3. On the Continuity of Ladder Operators

In this subsection, we would make a brief analysis on the continuity of the ladder operators when they act between the spaces of the type ${\mathrm{\Psi }}_l$. It is noteworthy that equations (48-) also hold with the spaces ${\mathrm{\Psi }}_l$ with continuity. Just need to check this result with the new version of (48), since the other cases go similarly. The equivalent formula to (48) is

$$A_l^{-}:{\mathrm{\Psi }}_l⟼{\mathrm{\Psi }}_{l+1}$$

(84)

Now, take $\psi \left(\rho \right)$ as in (61) and assume that $\psi \left(\rho \right)\in {\mathrm{\Psi }}_l$. Then ($l\ge 0$) after (52), we have for $s=0,1,2,\cdots$

$$q_{s,l+1}^2\left(A_l^{-}\psi \right)=4\sum _{p=0}^{\infty }p{|{\alpha }_{l,p}|}^2{[(l+1+2p+1)!]}^{2s}.$$

(85)

Clearly,

$$(l+2p+2)!=(l+2p+2)\left[\right(l+2p+1)!]=2(l/2+p)\left[\right(l+2p+1)!].$$

(86)

Since

$${\displaystyle \frac{l}{2}}+p<l+2p+1<(l+2p+1)!\mathrm{and}p<(l+2p+1)!,$$

(87)

we have that

$$p{[(l+1+2p+1)!]}^{2s}\le 2^{2s}{[(l+2p+1)!]}^{4s+2},$$

(88)

so that

$$q_{s,l+1}^2\left(A_l^{-}\psi \right)\le 2^{2s+2}\sum _{p=0}^{\infty }{|{\alpha }_{l,p}|}^2{[(l+2p+1)!]}^{4s+2}=2^{2s+2}{q}_{2s+1,l}^2\left(\psi \right)\forall \psi \in {\mathrm{\Phi }}_l.$$

(89)

This shows, at the same time, that $A_l^{-}\psi \in {\mathrm{\Psi }}_l$, for all $\psi \in {\mathrm{\Psi }}_l$ and that the mapping, $A_l^{-}$ as in (83) is continuous. Similar proofs will go for the equivalent of mappings (-).

Note that the interest of the topology introduced here lately by (62) is due to a proper characterization of the functions (66) as continuous functionals. For the purpose of the continuity of the ladder operators is not only enough but it also looks more appropriate the use of the topology introduced by the norms (43). The topology given by (62) makes the spaces ${\mathrm{\Psi }}_l$ too small and the topology itself is too strong. This is not a serious inconvenience for the characterization of the ladder operators as continuous operators, although it seems an excessive use of the resources. We have just pointed out that the topology given by norms (43) seems not to be appropriate to prove the continuity of the functionals ${\mathrm{\Phi }}_{-l,q}\left(\rho \right)$.

4.4. The Laguerre-Gauss Modes

Along the present subsection, we intend to recall a discussion already presented in [16]. That paper deals fundamentally with the physical aspects of the problem. Here, we wish to make a series of comments from the mathematical point of view. Let us go back to (29), which is explicitly given in [16] as

$$U_{l,p}(\theta ,z,\rho )={\displaystyle \frac{e^{il\theta }}{2\pi }}exp\left[-iz{\beta }_{|l|,p}\right]{\mathrm{\Phi }}_{|l|,p}\left(\rho \right),$$

(90)

where ${\beta }_{|l|,p}$ is given by (33). For each fixed value of $z\in \mathbb{R}$, the functions $U_{l,p}(\theta ,z,\rho )$, $l\in \mathbb{Z}$ and $p=0,1,2,\cdots$ form an orthonormal basis on $L^2\left({\mathbb{R}}^2\right)$. These functions are called the Laguerre-Gauss (LG) modes. Thus, for each fixed z, the orthogonality (we omit the dependence on the variables for simplicity)

$${\int }_{-\pi }^{\pi }{\int }_0^{\infty }{U}_{l,p}{U}_{l^{\prime },p}\rho d\rho d\theta ={\delta }_{l,l^{\prime }}{\delta }_{p,p^{\prime }},$$

(91)

and completeness

$$\sum _{l\in \mathbb{Z}}\sum _{p=0}^{\infty }{U}_{l,p}(\theta ,\rho )U_{l,p}({\theta }^{\prime },{\rho }^{\prime })=\delta (\theta -{\theta }^{\prime })\delta (\rho -{\rho }^{\prime })$$

(92)

relations hold.

In [16] was defined the following operator for each fixed value of $l\in \mathbb{Z}$ (we have chosen $w_0^2=2$ in the formulas given in [16]):

$${\mathcal{L}}_l:={\rho }^{-1/2}{L}_l{\rho }^{1/2},\mathrm{with}{L}_l={\displaystyle \frac{1}{2}}\left[-{\displaystyle \frac{{\partial }^2}{\partial {\rho }^2}}+{\displaystyle \frac{(l+1/2)(l-1/2)}{{\rho }^2}}\right]+{\rho }^2$$

(93)

It is proven in [16] that the LG modes $U_{l,p}(\theta ,\rho )$, $p=0,1,2,\cdots$ are eigenfunctions of ${\mathcal{L}}_l$ for each fixed value of $l=0,1,2,\cdots$ with eigenvalue ${\displaystyle \frac{1}{2}}{\beta }_{|l|,p}$, i.e.,

$${\mathcal{L}}_l{U}_{l,p}={\displaystyle \frac{1}{2}}{\beta }_{|l|,p}{U}_{l,p},p=0,1,2,\cdots ,l\in \mathbb{Z}.$$

(94)

Let us call ${\mathcal{H}}_l$ the subspace of $L^2\left({\mathbb{R}}^2\right)$ by the functions $U_{l,p}(\theta ,\rho )$, $p=0,1,2,\cdots$ and $l\in \mathbb{Z}$ fixed, which form an orthonormal basis for ${\mathcal{H}}_l$. Obviously $L^2\left({\mathbb{R}}^2\right)={\oplus }_{l\in \mathbb{Z}}{\mathcal{H}}_l$. The functions $U_{l,p}(\theta ,\rho )$ differ from the functions ${\mathrm{\Phi }}_{|l|,p}\left(\rho \right)$ in (30) just for a phase and a multiplicative constant. Thus, they could be identified by an “abuse de langage” and hence ${\mathcal{H}}_l$ and $L^2({\mathbb{R}}^{+},\rho d\rho )$ for each value of l. Any ${\psi }_l(\theta ,\rho )\in {\mathcal{H}}_l$ with fixed l has the form:

$${\psi }_l(\theta ,\rho )=\sum _{p=0}^{\infty }{a}_{l,p}{U}_{l,p}(\theta ,\rho )\mathrm{with}\sum _{p=0}^{\infty }{|a_{l,p}|}^2<\infty .$$

(95)

Then, let us consider the subspace of functions ${\psi }_l$ as in (95) such that

$$\sum _{p=0}^{\infty }|a_{l,p}{|}^2{(|l|+2p+1)}^{2s}<\infty ,s=0,1,2,\cdots .$$

(96)

This is a subspace of ${\mathcal{H}}_l$ which may be identify with ${\tilde{\mathrm{\Phi }}}_l$ after the aforementioned “abuse de langage”, so that ${\tilde{\mathrm{\Phi }}}_l\subset {\mathcal{H}}_l\subset {\tilde{\mathrm{\Phi }}}_l^{\times }$ is a RHS for all values of l. If we define the action of ${\mathcal{L}}_l$ on each ${\psi }_l\in {\tilde{\mathrm{\Phi }}}_l$, for fixed l, as

$${\mathcal{L}}_l{\psi }_l(\theta ,\rho ):={\displaystyle \frac{1}{2}}\sum _{p=0}^{\infty }{a}_{l,p}(|l|+2p+1)U_{l,p}(\theta ,\rho ),$$

(97)

we immediately see using the arguments at the beginning of the present Section that (90) is well defined and that ${\mathcal{L}}_l$ is continuous on ${\tilde{\mathrm{\Phi }}}_l$.

Ladder operators, ${\mathcal{L}}_l^{\pm }$, have been also defined in [16] and this definition can be extended to ${\mathrm{\Phi }}_l$ as follows:

$${\mathcal{L}}_l^{\pm }{\psi }_l(\theta ,\rho )=\sum _{p=0}^{\infty }{a}_{l,p}{U}_{l,p\pm 1},\mathrm{with}{a}_{l,-1}=0.$$

(98)

Same arguments prove that ${\mathcal{L}}_l^{\pm }$ are continuous on ${\tilde{\mathrm{\Phi }}}_l$.

This construction may be extended trivially to $L^2\left({\mathbb{R}}^2\right)$. In fact, if $\psi (\theta ,\rho )\in L^2\left({\mathbb{R}}^2\right)$, we have

$$\psi (\theta ,\rho )=\sum _{l=-\infty }^{\infty }\sum _{p=0}^{\infty }{a}_{l,p}{U}_{l,p}(\theta ,\rho ),\mathrm{with}\sum _{l=-\infty }^{\infty }\sum _{p=0}^{\infty }{|a_{l,p}|}^2<\infty .$$

(99)

Then, select the subspace $\mathbf{\Phi }$ of all vectors, $\psi (\theta ,\rho )$, in $L^2\left({\mathbb{R}}^2\right)$ such that, for all $l\in \mathbb{Z}$ and $p=0,1,2,\cdots$,

$$\sum _{l=-\infty }^{\infty }\sum _{p=0}^{\infty }|a_{l,p}{|}^2{(|l|+2p+1)}^{2s}<\infty ,s=0,1,2,\cdots .$$

(100)

Since $\mathbf{\Phi }$ contains all the LG modes $\left\{U_{l,p}\right\}$, $l\in \mathbb{Z}$, $p=0,1,2,\cdots$, it is dense in $L^2\left({\mathbb{R}}^2\right)$. Since $s=0$ gives the Hilbert space topology, then the topology on $\mathbf{\Phi }$ given by the seminorms for which their squares are given in (100), the topology on $\mathbf{\Phi }$ is finer that the Hilbert space topology. Consequently,

$$\mathbf{\Phi }\subset L^2\left({\mathbb{R}}^2\right)\subset {\mathbf{\Phi }}^{\times }$$

(101)

is a RHS. On $\mathbf{\Phi }$ are continuous the following operators:

$$\left[\underset{l=-\infty }{\overset{\infty }{⨁}}{\mathcal{L}}_l\right]\psi (\theta ,\rho ):={\displaystyle \frac{1}{2}}\sum _{l=-\infty }^{\infty }\sum _{p=0}^{\infty }{a}_{l,p}(|l|+2p+1)U_{l,p}(\theta ,\rho ),$$

(102)

and

$$\left[\underset{l=-\infty }{\overset{\infty }{⨁}}{\mathcal{L}}_l^{\pm }\right]\psi (\theta ,\rho ):=\sum _{l=-\infty }^{\infty }\sum _{p=0}^{\infty }{a}_{l,p}{U}_{l,p\pm 1}(\theta ,\rho ),$$

(103)

with $a_{l,-1}=0$ for all value of l. Proofs are as in previous examples.

5. Coherent States and Continuous Generators for the Algebra $su(1,1)$

Let us give a brief analysis of different types of coherent states from the point of view of the above defined Gelfand triplets. To begin with, let us define the so called index free operators, which have been introduced in [16]. These are ladder operators defined on ${\mathrm{\Phi }}_l$ for fixed values of $l>0$. On the basis vectors of ${\mathrm{\Phi }}_l$ these operators act as follows [7]:

$$A^{+}{\psi }_{l,p}:=\sqrt{p+1}{\psi }_{l,p+1}\left(\rho \right),A^{-}{\psi }_{l,p}:=\sqrt{p}{\psi }_{l,p-1}\left(\rho \right),$$

(104)

so that their action on an arbitrary vector in ${\mathrm{\Phi }}_l$, $l>0$, of the form (41) with (43) takes the form:

$$A^{+}\phi \left(\rho \right):=\sum _{p=0}^{\infty }{a}_{l,p}\sqrt{p+1}{\psi }_{l,p+1}\left(\rho \right),A^{-}\phi \left(\rho \right):=\sum _{p=1}^{\infty }{a}_{l,p}\sqrt{p}{\psi }_{l,p-1}\left(\rho \right).$$

(105)

Following similar techniques than those given in Section 4.1, we easily show that both ladder operators $A^{\pm }$ are well defined and continuous on each of the ${\mathrm{\Phi }}_l$, $l=1,2,\cdots$.

Since, we are focusing our attention on coherent states, we are mainly interested on the operator $A^{-}$, since coherent states are eigenvectors of the annihilation operator with arbitrary eigenvalue. Recall that coherent states are the eigenstates of the annihilation operator, so that if ${\eta }^{\alpha }\left(\rho \right)$ is a coherent state $A^{-}{\eta }^{\alpha }\left(\rho \right)=\alpha {\eta }^{\alpha }\left(\rho \right)$ with $\alpha \in \mathbb{C}$. Clearly coherent states must satisfy the following relation:

$$A^{-}{\eta }^{\alpha }\left(\rho \right)=\sum _{p=0}^{\infty }{b}_{l,p}\sqrt{p}{\psi }_{l,p-1}\left(\rho \right)=\alpha \sum _{p=0}^{\infty }{b}_{l,p}{\psi }_{l,p}\left(\rho \right).$$

(106)

A straightforward calculation gives

$${\eta }^{\alpha }\left(\rho \right)\equiv {\eta }_l^{\alpha }\left(\rho \right)=\sum _{p=0}^{\infty }{\displaystyle \frac{{\alpha }^p}{\sqrt{p!}}}{e}^{-{|\alpha |}^2/2}{\psi }_{l,p}\left(\rho \right).$$

(107)

Thus, coherent states are defined for each value of $l=1,2,\cdots$. Furthermore,

Proposition.- For fixed $\alpha \in \mathbb{C}$ and $l=1,2,\cdots$, ${\eta }_l^{\alpha }\left(\rho \right)\in {\mathrm{\Phi }}_l$.

Proof.- We just need to show that

$$\sum _{p=0}^{\infty }{\displaystyle \frac{{|\alpha |}^{2p}}{p!}}{e}^{-{|\alpha |}^2}{(l+2p+1)}^{2s}<\infty ,s=0,1,2,\cdots .$$

(108)

Note that the term $e^{-{|\alpha |}^2}$ in (108) is a constant and, hence, irrelevant. To show this, it is sufficient to prove that for large values of p we have17

$${\displaystyle \frac{{|\alpha |}^{2p}}{p!}}{(l+2p+1)}^{2s}<{\displaystyle \frac{1}{p^2}}s=0,1,2,\cdots .$$

(109)

Note that, according to the Stirling formula,

$$p!\sim \sqrt{2\pi p}{\left({\displaystyle \frac{p}{e}}\right)}^p,$$

(110)

so that it is enough to show that for large values of p,

$${\displaystyle \frac{a^p}{\sqrt{2\pi p}{p}^p}}{(l+2p+1)}^{2s}{p}^2<1,a={|\alpha |}^{2}e.$$

(111)

Then, it is enough to show two inequalities valid for large values of p. The former is

$$a^p<p^{p/2}⟺p>a^2,$$

(112)

for any fixed value of the positive constant a and large p. The second one is

$${(l+2p+1)}^{2s}{p}^2<p^{p/2},$$

(113)

which is straightforward for p large. Then, (111) comes and, hence, (108). □

These coherent states are often referred as Glauber-Klauder-Sudarshan coherent states. Note that the coherent state ${\eta }_l^{\alpha }\left(\rho \right)\notin {\mathrm{\Psi }}_l$, since for $s=3$,

$$\sum _{p=0}^{\infty }{\displaystyle \frac{{|\alpha |}^{2p}}{{(p!)}^2}}{[(l+2p+1)!]}^6=\infty .$$

(114)

Next, let us consider a new set of ladder operators, ${\mathcal{A}}_l^{\pm }$, formally defined as

$${\mathcal{A}}_l^{\pm }:=\left({\rho }^2\mp \rho {\displaystyle \frac{d}{d\rho }}-\left(l+2p\mp {\displaystyle \frac{1}{2}}\right)\right).$$

(115)

Their action on the basis of functions ${\psi }_{l,p}\left(\rho \right)$ is given by

$${\mathcal{A}}_l^{+}{\psi }_{l,p-1}\left(\rho \right)=\sqrt{p(l+p)}{\psi }_{l,p}\left(\rho \right),{\mathcal{A}}_l^{-}{\psi }_{l,p}\left(\rho \right)=\sqrt{p(l+p)}{\psi }_{l,p-1}\left(\rho \right).$$

(116)

Their extension to any vector $\phi \left(\rho \right)\in {\mathrm{\Phi }}_l$ is straightforward. This extensions are well defined and continuous on ${\mathrm{\Phi }}_l$.

As in the previous case, we may here define coherent states for each positive value of l, ${\psi }_{l,\alpha }^{\mathrm{BG}}\left(\rho \right)$, as ${\mathcal{A}}_l^{-}{\psi }_{l,\alpha }^{\mathrm{BG}}\left(\rho \right)=\alpha {\psi }_{l,\alpha }^{\mathrm{BG}}\left(\rho \right)$, $\forall \alpha \in \mathbb{C}$. They are called the Barut-Girardello coherent states. They have the following form:

$${\psi }_{l,\alpha }^{\mathrm{BG}}\left(\rho \right)=\sum _{p=0}^{\infty }{\displaystyle \frac{{\alpha }^p}{\sqrt{p!(l+p)!}}}{\displaystyle \frac{{\alpha }^{l/2}}{\sqrt{I_l\left(2\alpha \right)}}}{\psi }_{l,p}\left(\rho \right),$$

(117)

where $I_l\left(x\right)$ is the modified Bessel function. Again, ${\psi }_{l,\alpha }^{\mathrm{BG}}\left(\rho \right)\in {\mathrm{\Phi }}_l$ and ${\psi }_{l,\alpha }^{\mathrm{BG}}\left(\rho \right)\notin {\mathrm{\Psi }}_l$.

Finally, we have the Perelomov coherent states. Their construction requieres some analysis. First of all, let us define a new operator ${\mathcal{A}}_0$ as

$$2{\mathcal{A}}_l^0{\psi }_{l,p}\left(\rho \right):=[{\mathcal{A}}_l^{-},{\mathcal{A}}_l^{+}]{\psi }_{l,p}\left(\rho \right)=(l+2p+1){\psi }_{l,p}\left(\rho \right).$$

(118)

Clearly, ${\mathcal{A}}_l^0$ may be extended to ${\mathrm{\Phi }}_l$ with continuity. Note that the operators ${\mathcal{A}}_l^{\pm }$ and ${\mathcal{A}}_l^0$ give a system of generators on each of the ${\mathrm{\Phi }}_l$ of the Lie algebra $su(1,1)$, since

$$[{\mathcal{A}}_l^0,{\mathcal{A}}_l^{\pm }]=\pm {\mathcal{A}}_l^{\pm },[{\mathcal{A}}_l^{-},{\mathcal{A}}_l^{+}]=2{\mathcal{A}}_l^0.$$

(119)

In consequence, each of the ${\mathrm{\Phi }}_l$ supports an irreducible representation of the algebra $su(1,1)$ given by continuous operators with the topology on ${\mathrm{\Phi }}_l$. Note that the operators on the corresponding enveloping algebras are also continuous on each of the ${\mathrm{\Phi }}_l$. Their canonical extension to the duals ${\mathrm{\Phi }}_l^{\times }$ are also continuous with any topology on ${\mathrm{\Phi }}_l^{\times }$ compatible with duality.

In order to construct the Perolomov coherent states, let us consider for each $l>0$ the ground state ${\psi }_{l,0}\left(\rho \right)$ and define

$${\xi }_l\left(\rho \right):=D(\lambda ,t){\psi }_{l,0}\left(\rho \right),$$

(120)

where $\lambda$ is a complex number and

$$D(\lambda ,t):=e^{t(\lambda {\mathcal{A}}_l^{+}-{\lambda }^{\ast }{\mathcal{A}}_l^{-})}=e^{\alpha {\mathcal{A}}_l^{+}}{e}^{\beta {\mathcal{A}}_l^0}{e}^{\gamma {\mathcal{A}}_l^{-}}=K_{+}{K}_0{K}_{-},$$

(121)

where the last identity in (121) is just the definition of $K_{\pm }$ and $K_0$ and the second one is a consequence of the Baker-Campbell-Hausdorff formula. In order to go further, let us derive with respect to the parameter t, so that

$$\begin{array}{c}\hfill {\displaystyle \frac{d}{dt}}D(\lambda ,t)=(\lambda {\mathcal{A}}_l^{+}-{\lambda }^{\ast }{\mathcal{A}}_l^{-})K_{+}{K}_0{K}_{-}\\ \hfill =K_{+}{K}_0{\gamma }^{\prime }{\mathcal{A}}_l^{-}{K}_{-}+{\alpha }^{\prime }{\mathcal{A}}_l^{+}{K}_{+}{K}_0{K}_{-}+K_{+}{\mathcal{A}}_l^0{\beta }^{\prime }{\mathcal{A}}_l^0{K}_0{K}_{-}.\end{array}$$

(122)

Then using the commutation relations and some algebra, we arrive to

$$\lambda ={\alpha }^{\prime }-\alpha {\beta }^{\prime }+{\alpha }^2{e}^{-\beta }{\gamma }^{\prime },{\beta }^{\prime }=2\alpha {\gamma }^{\prime },{\lambda }^{\ast }=-e^{-\beta }{\gamma }^{\prime },$$

(123)

where the star denotes conjugation. If we now define

$$\zeta :={\displaystyle \frac{\lambda }{|\lambda |}}tanh|\lambda |,$$

(124)

we have

$$\alpha =\zeta ,\beta ={log(1-|\zeta |}^2),\gamma =-{\zeta }^{\ast }.$$

(125)

Thus, we have all parameters in terms of $\zeta$. After some algebra, we can obtain the explicit expressions for the Perolomov coherent states for each value of $l>0$ as

$${\xi }_l^{\zeta }\left(\rho \right)=\sum _{p=0}^{\infty }{\displaystyle \frac{{\zeta }^p}{\sqrt{p!}}}{log(1-|\zeta |}^2{)}^{(l+1)}{\psi }_{l,p}\left(\rho \right).$$

(126)

Analogously, ${\xi }_l^{\zeta }\left(\rho \right)\in {\mathrm{\Phi }}_l$ and ${\xi }_l^{\zeta }\left(\rho \right)\notin {\mathrm{\Psi }}_l$ for $l=1,2,\cdots$. Nevertheless, after (65), it is clear that ${\eta }_l^{\alpha }\left(\rho \right),{\psi }_{l,\alpha }^{\mathrm{BG}}\left(\rho \right),{\xi }_l^{\zeta }\left(\rho \right)\in {\mathrm{\Psi }}_l^{\times }$, for all positive values of l.

5.1. Resolutions of the Identity

In this Subsection, we briefly analyze the meaning and some properties of the resolutions of the identity given by the coherent states (107). A standard calculation with $d\alpha =dxdy$ shows that for it fixed $l>0$, one has

$$\begin{array}{c}\hfill {\displaystyle \frac{1}{\pi }}{\int }_{\mathbb{C}}|{\eta }_l^{\alpha }\left(\rho \right)\rangle \langle {\eta }_l^{\alpha }\left(\rho \right)|d\alpha =\sum _{p,q=0}^{\infty }|{\psi }_{l,p}\left(\rho \right)\rangle \langle {\psi }_{l,q}\left(\rho \right)|{\displaystyle \frac{1}{\pi }}{\int }_{\mathbb{C}}{\displaystyle \frac{{\left({\alpha }^{\ast }\right)}^p}{\sqrt{p!}}}{\displaystyle \frac{{\left(\alpha \right)}^q}{\sqrt{q!}}}{e}^{-{|\alpha |}^2}d\alpha \\ \hfill =\sum _{p,q=0}^{\infty }|{\psi }_{l,p}\left(\rho \right)\rangle \langle {\psi }_{l,q}\left(\rho \right)|{\displaystyle \frac{p!}{\sqrt{p!q!}}}{\delta }_{p,q}=\sum _{p=0}^{\infty }|{\psi }_{l,p}\left(\rho \right)\rangle \langle {\psi }_{l,p}\left(\rho \right)|=I,\end{array}$$

(127)

where I is the identity on $L^2\left({\mathbb{R}}^{+}\right)$.

Note that the resolution of the identity given by the first integral in (127) not only represents an identity on the Hilbert space $L^2\left({\mathbb{R}}^{+}\right)$. It is, in addition, a representation of the identity mapping, $I_l$, from ${\mathrm{\Phi }}_l$ into its dual ${\mathrm{\Phi }}_l^{\times }$. In fact, let us take an arbitrarily fixed ${\phi }_l\left(\rho \right)\in {\mathrm{\Phi }}_l$ and consider

$${\displaystyle \frac{1}{\pi }}{\int }_{\mathbb{C}}|{\eta }_l^{\alpha }\left(\rho \right)\rangle \langle {\eta }_l^{\alpha }\left(\rho \right)|{\phi }_l\left(\rho \right)\rangle d\alpha =\sum _{p=0}^{\infty }|{\psi }_{l,p}\left(\rho \right)\rangle \langle {\psi }_{l,p}\left(\rho \right)|{\phi }_l\left(\rho \right)\rangle ={\phi }_l\left(\rho \right).$$

(128)

Thus, $I_l{\phi }_l={\phi }_l$. Let us apply the first term in (128) to an any ${\varphi }_l\left(\rho \right)\in {\mathrm{\Phi }}_l$. We have,

$${\displaystyle \frac{1}{\pi }}{\int }_{\mathbb{C}}\langle {\varphi }_l\left(\rho \right)|{\eta }_l^{\alpha }\left(\rho \right)\rangle \langle {\eta }_l^{\alpha }\left(\rho \right)|{\phi }_l\left(\rho \right)\rangle d\alpha =\sum _{p=0}^{\infty }\langle {\varphi }_l\left(\rho \right)|{\psi }_{l,p}\left(\rho \right)\rangle \langle {\psi }_{l,p}\left(\rho \right)|{\phi }_l\left(\rho \right)\rangle =\langle {\varphi }_l|{\phi }_l\rangle .$$

(129)

Therefore, $I_l{\phi }_l$ acts on each of the ${\varphi }_l\left(\rho \right)\in {\mathrm{\Phi }}_l$ as a linear functional, $F_{{\phi }_l}$, so that $F_{{\phi }_l}\left({\varphi }_l\left(\rho \right)\right)\equiv \langle {\varphi }_l|I_l{\phi }_l\rangle =\langle {\varphi }_l|{\phi }_l\rangle$. Thus,

$$|F_{{\phi }_l}\left({\varphi }_l\left(\rho \right)\right)|=|\langle {\varphi }_l|{\phi }_l\rangle |\le ||{\varphi }_l||||{\phi }_l||=C_l{p}_{l,0}\left({\varphi }_l\right),\forall {\varphi }_l\left(\rho \right)\in {\mathrm{\Phi }}_l,$$

(130)

where $C_l\equiv ||{\phi }_l||$. Note that ${\phi }_l\left(\rho \right)\in {\mathrm{\Phi }}_l$ has been chosen to be fixed, which proves that $F_{{\phi }_l}\in {\mathrm{\Phi }}_l^{\times }$ and this is true for each ${\phi }_l\left(\rho \right)\in {\mathrm{\Phi }}_l$. It is customary to identify $F_{{\phi }_l}$ with ${\phi }_l$ and this is the true meaning of the identity $I_l:{\mathrm{\Phi }}_l⟼{\mathrm{\Phi }}_l^{\times }$.

The weak topology on ${\mathrm{\Phi }}_l^{\times }$ is given by the set of seminorms, $q_{{\psi }_l}$, which are given by $q_{{\psi }_l}\left({\phi }_l\right):=|F_{{\psi }_l}\left({\phi }_l\right)|=|\langle {\psi }_l|{\phi }_l\rangle |$, $\forall {\phi }_l\in {\mathrm{\Phi }}_l$. Each vector ${\psi }_l\in {\mathrm{\Phi }}_l$ defines a unique seminorm $q_{{\psi }_l}$, the correspondence ${\psi }_l⟼q_{{\psi }_l}$ is one to one and vectors ${\psi }_l$ run out the space ${\mathrm{\Phi }}_l$. Then, for all ${\phi }_l\left(\rho \right)\in {\mathrm{\Phi }}_l$ one has

$$q_{{\psi }_l}\left(I_l{\phi }_l\right)=q_{{\psi }_l}\left(F_{{\phi }_l}\right)=|\langle {\psi }_l|{\phi }_l\rangle |=K_l||{\phi }_l||=K_k{p}_{l,0}\left({\phi }_l\right),$$

(131)

which shows the continuity of the canonical identity mapping $I_l:{\mathrm{\Phi }}_l⟼{\mathrm{\Phi }}_l^{\times }$. Note that in this argument, we may find also the proof that the identity mapping $I:L^2\left({\mathbb{R}}^{+}\right)⟼{\mathrm{\Phi }}_l^{\times }$ is also continuous. Needless to say that canonical identity maps are linear.

A resolution of the identity for (117) has been obtained in [23]. It has the form for each $l>0$

$${\int }_{\mathbb{C}}|{\psi }_{l,\alpha }^{\mathrm{BG}}\left(\rho \right)\rangle \langle {\psi }_{l,\alpha }^{\mathrm{BG}}\left(\rho \right)|d\sigma \left(\alpha \right)=\sum _{p=0}^{\infty }|{\psi }_{l,p}\left(\rho \right)\rangle \langle {\psi }_{l,p}\left(\rho \right)|=I.$$

(132)

Here, $d\sigma \left(\alpha \right)=\sigma \left(r\right)rdrd\theta$, $r:=|\alpha |$, is a measure on the complex plane. The function $\sigma \left(r\right)$ is a product of a constant times a power of r times a Bessel function of third kind with argument proportional to r. Properties of this resolution of the identity are similar as in the previous case.

We do not have a resolution of the identity for the coherent states of the form (126). Nevertheless, a resolution of the identity for the coherent states obtained for the representation of the discrete series of $su(1,1)$ exists [24]. This would requiere a discussion which does not fit within the context of the present article.

6. Concluding Remarks

In a recent series of articles, it has been shown that Gelfand triplets, also named as rigged Hilbert spaces, are the precise mathematical framework that includes several tools currently used in Quantum Mechanics. These are discrete and continuous bases, a representation of Lie algebras by continuous operators and discrete bases given by specific special functions.

A former and well know example considers the Heisenberg-Weyl Lie algebra. This algebra includes the momentum and poistion operators along the harmonic oscillator ladder operators, which are continuous on the Schwartz space (and not on Hilbert space). Here, we take as special functions the normalized Hermite functions [25]. This example has been the guide to analyze other systems in which Legendre, Laguerre, Zernike or Jacobi spacial functions take the place of Hermite functions [26,27]. All the mentioned special functions are very important in a wide range of quantum problems.

Recent works on Quantum Optics reveal to importance of the Laguerre-Gauss special functions. Following the ideas of previous works, we construct suitable Gelfand triplets in which ladder operators, Hamiltonians and other operators are continuous on spaces spanned by Laguerre-Gauss functions for each fixed positive values of the parameter l.

A refinement on the topology of test functions shows that Laguerre-Gauss functions with negative values of l, ${\mathrm{\Phi }}_{-l,q}\left(\rho \right)$, can be looked as functionals on the space of test functions labelled by $l=|-l|$. Thus, we have two types of Gelfand triplets, one in which the topology of the space of test functions is similar to the topology of the Schwartz space, as proven in Appendix II and another one with a more restrictive topology and, therefore, smaller.

We also show that different types of coherent states can be constructed with the help of the Laguerre-Gauss special functions and that belong to the test spaces with topology similar to the Schwartz space, although do not belong to the more restrictive subspace. In the derivation of the Perelomov coherent states, we include a representation of $su(1,1)$ by continuous operators.