Hermite Functions and Fourier Series

: Using normalized Hermite functions, we construct bases in the space of square integrable functions on the unit circle ( L 2 ( C ) ) and in l 2 ( Z ) , which are related to each other by means of the Fourier transform and the discrete Fourier transform. These relations are unitary. The construction of orthonormal bases requires the use of the Gramm–Schmidt method. On both spaces, we have provided ladder operators with the same properties as the ladder operators for the one-dimensional quantum oscillator. These operators are linear combinations of some multiplication-and differentiation-like operators that, when applied to periodic functions, preserve periodicity. Finally, we have constructed riggings for both L 2 ( C ) and l 2 ( Z ) , so that all the mentioned operators are continuous.


Introduction
Hermite functions have been an important tool in the development of elementary quantum mechanics as solutions of the quantum non-relativistic harmonic oscillator [1]. From a mathematical point of view, Hermite functions serve as an orthonormal basis (complete orthonormal set) for the Hilbert space L 2 (R). They are products of Hermite polynomials times and a Gaussian, so they are functions which are strongly localized near the origin [2,3].
The Fourier transform is a unitary operator on L 2 (R). The Hermite functions are its eigenfunctions and allow a division of L 2 (R) into four eigenspaces related to the cyclic group C 4 . This division can be relevant in applications.
The one-dimensional Fourier transform and its inverse are automorphisms on the Hilbert space L 2 (R), which preserve the Hilbert space norm, after the Plancherel theorem [4]. This result can be extended to some other spaces of interest in physics such as the space of infinitely differentiable functions converging to zero at the infinite faster than the inverse of a polynomial and the space of tempered distributions. In both cases the Fourier transform and the inverse Fourier transform are automorphisms, which are continuous with the standard topologies defined in both spaces [4].
Loosely speaking, the Fourier series may be looked at as a particular case of a span of square integrable functions on a finite interval [a, b] in terms of square integrable functions on this interval and extended beyond by periodicity. One usually takes the interval [0, 2π] or [−π, π]. Then, one may look at the Fourier series as mappings from the space of C ∞ functions on the unit circle to discrete functions over the set of integers Z. This mapping is invertible, which means that a sequence of complex numbers {a n } n∈Z with the property that ∑ n∈Z |a n | 2 < ∞ uniquely fixes (almost elsewhere) a square integrable function on the unit circle [5]. These numbers are obtained using a discrete-time Fourier transform [6] over the original function, so that the Fourier series and discrete-time Fourier transforms may be considered as operations inverse of each other.
In the present article, we construct a complete sequence of periodic functions using the Hermite functions, which is a non-orthonormal basis on L 2 [−π, π] ≡ L 2 (C), where C is the unit circle. Then, after the Gram-Schmidt procedure we obtain an orthonormal basis formed by periodic functions. All functions on this orthonormal basis can be spanned into a Fourier series with coefficients obtained from the Hermite functions. Vice-versa, these coefficients are obtained via the discrete Fourier transform of the functions belonging to the orthonormal basis.
We proceed to an equivalent construction on the space l 2 (Z) of square summable complex sequences indexed by the set of integer numbers. Each term of the sequence is given by the value on an integer of a given normalized Hermite function. Fourier series and Fourier transform give the equivalence between both constructions, which induces a unitary mapping between both spaces.
We also construct multiplication and differentiation operators on a subspace of L 2 (C) that includes our complete sequence of periodic functions in a way that these operators preserve periodicity. They play a similar role to that played by multiplication and derivation operators in the description of the one-dimensional quantum harmonic oscillator. In particular, they provide ladder operators for the complete sequence of periodic functions analogous to the creation and annihilation operators for the quantum oscillator. We may also construct a locally convex space of test functions, S, including the chosen complete sequence of periodic functions, such that these operators are continuous on S. This space, being dense in L 2 (C) with continuous injection, defines a rigged Hilbert space S ⊂ L 2 (C) ⊂ S × , so that the operators can be continuously extended to the dual S × . A similar construction is also possible for l 2 (Z).
The rigging of Hilbert spaces is necessary for the continuity of the relevant operators that we introduce in this presentation and this includes the new ladder operators. The formalism of rigged Hilbert spaces was introduced by Gelfand [7]. Although this is not quite familiar to many theoretical physicists, it has acquired more and more importance in the field of mathematical physics [8][9][10][11][12][13] and even in mathematics [14][15][16][17][18][19][20].
In a series of previous articles [21][22][23][24][25], we have discussed the relations between discrete and continuous basis, algebras of operators, special functions and rigged Hilbert spaces (also called Gelfand triplets). We have shown that all these concepts properly convive in the framework of rigged Hilbert spaces and not on Hilbert spaces. The present paper fits also in part within the same context, where our special functions are now series constructed from Hermite functions. In addition, it can be shown that the structure of rigged Hilbert spaces is suitable for the representation of our operators as an algebra of continuous operators on the same domain, while on Hilbert space these operators are unbounded. Thus, we create a representation with some formal and mathematical advantages derived from this continuity on a common domain.
The organization of this paper is as follows: In Section 2, we study the behavior of the Hermite functions under the Fourier transform. The well-defined parity of the Hermite functions allows that the even/odd Hermite functions span the subspace of the even/odd functions of L 2 (R). In a similar way, since the Hermite functions are eigenvectors of the Fourier transform with eigenvalues the four roots of unity, L 2 (R) can be split in four subspaces of functions, each one characterised by one eigenvalue of the Fourier transform as an operator on L 2 (R). In Section 3, we introduce a set of functions defined in the unit circle in terms of the Hermite functions. In Sections 4 and 5, we present some results with the aim of given a unitary view that include different concepts such as Fourier transform, Fourier series, discrete Fourier transform and Hermite functions. Later in Sections 6 and 7, we define ladder operators and those others related to them and the equipations of Hilbert spaces, or Gelfand triplets, or rigged Hilbert spaces on which these operators are well defined as continuous operators. We finish the present article with concluding remarks plus two Appendices. In Appendix A, we show two results which are important in the development of the material in Section 5. In Appendix B, we construct orthonormal systems in L 2 (C) and l 2 (Z) via the Gramm-Schmidt process.

The Hermite Functions and the Fourier Transform
Let us consider the normalized Hermite functions in one-dimension, sometimes also called the Gauss-Hermite functions [26][27][28][29]. As is well known, they have the following form where H n (x) are the Hermite polynomials, x ∈ R and n ∈ N, with N the set of non-negative integers, 0, 1, 2, . . . . It is well-known that the Hermite functions form an orthonormal basis in the Hilbert space of square integrable functions on the line, L 2 (R). In other words, they verify the following orthonormality and completeness relations Observe that the completeness for the Hermite functions is given and labelled by the natural numbers, so that n ∈ N, although the usual definition of the Parseval identity is a sum on the set of integers. This is not a contradiction, so both the set of natural numbers N (including the zero) and the set of integers Z are both infinite countable. One may associate the even Hermite functions, ψ 2n , to the positive integers and the odd Hermite functions, ψ 2n+1 , to the set of negative integers. The Hermite functions (1) are eigenfunctions of the Fourier Transform (FT) in the following sense [21,30] FT[ψ n (x), x, y] ≡ ψ n (y) : or more generally, This is the key issue of the present Section. We define the action of the cyclic group of order two, C 2 ≡ {I, P}, on the space of the square integrable functions on the line L 2 (R), where I is the identity and P is the reflection operator, with P 2 = I. Using the elements of C 2 , we can construct the following projectors that split L 2 (R) into two mutually orthogonal subspaces: the spaces of even and odd functions L 2 E (R) and L 2 O (R), respectively, Even or odd Hermite functions ψ n (x) are labeled by an even or odd index n, respectively. This shows that This property, together with the orthogonality of the Hermite functions on L 2 (R), suggests the following notation where the analogy with (7) is obvious. Next, let us use the property of the orthogonality (2) of the Hermite functions as well as their defined parity (8), so as to obtain the following identities We easily find from the first and the third rows of (10) that Observe that the former and second expressions in (11) gives completeness in the space of the even/odd functions, respectively. The group C 2 is a subgroup of C 4 , which is the cyclic group of the four roots of unity {1, i, −1, −i}. This comes from two simple facts: the set of Hermite functions {ψ n (x)} is a basis in L 2 (R) and the Hermite functions are eigenfunctions of the Fourier Transform (3). From the properties of the imaginary unit we easily obtain that Consequently, in analogy with (9), relations (12) yield to this identity having the form of direct sums Clearly, the subspaces spanned by the Hermite functions {ψ 4n+k (x)} with k = 0, 1, 2, 3 are orthonormal to each other. However, we know that each set of Hermite functions given by one of the values k = 0, 1, 2, 3 cannot fulfil the completeness relation. They satisfy instead In order to prove relations (14), let us consider the following Fourier transform (see notation in (3)) Then, we recall that the second identity in (2) means that, when we endow the space of the tempered distributions S with the weak topology where S is the Schwartz space, we have The Fourier transform is weakly continuous on S , hence we have and The uniqueness of the weak limit in S gives Obviously, (18) yields to After the transformation y → −y, we use the parity property of Hermite functions given in (8), so that (19) gives Summing (19) and (20), we obtain Subtracting (20) from (19), we obtain that Let us go back to (8), which can be now rewritten as Finally, let us sum (21) with the first identity of (23). The result is This is the first identity in (14). The validity of the other three relations in (14) is proven analogously. The proof is now complete. Next, let us consider an arbitrary function f (x) in L 2 (R). It allows an expansion in terms of the Hermite functions ψ n (x) as The even part of f (x), f E (x), allows an expansion in terms of even Hermite functions. Then, and taking into account (9), For the odd part, In All the above results together show that any function f (x) in L 2 (R) can be split into four parts Then, from (25) and the parity of the Hermite functions, it comes that the even part of f (x) has the following form Next, using the continuity of the Fourier transform in L 2 (R) and taking into account (25), Hence, From (32) and (34), we obtain the first expression of (31). Similarly, we prove all the other relations in (31). The projectors producing this splitting are They verify the following identities All these projections are orthogonal and the corresponding splitting of the Hilbert space L 2 (R) is given by We conclude this part of the discussion here.

Periodic Functions and Hermite Functions
As is well known and as we have mentioned before, the complete set of Hermite functions, {ψ n (x)} n∈N , forms an orthonormal basis on L 2 (R). Next, and using the Hermite functions, we shall construct a countable set of periodic functions that will be a system of generators of the space of square integrable functions on the unit circle.
The space L 2 (C) is the space of Lebesgue square integrable functions f (φ) : C −→ C. The norm of the function f (φ) is given by the following relation We intend to introduce a space of functions with given properties. Consider the angular variable −π ≤ φ < π and define The functions C n (φ) are obviously periodic with period 2π A first result about the convergence of the series defining C n (φ) is given by the following proposition. Proposition 1. The series defining each of the C n (φ) are absolutely convergent. Furthermore, each C n (φ) is bound on the interval −π ≤ φ < π and the square integrable on this interval.

Proof. Let us write (39) as
The first term in the right hand side of (40) is bound by ψ 0 (φ) ≤ 2 , ∀φ ∈ [−π, π), as we may check by using the Cramér inequality [31] (page 787, formula 22.14.17). Since both series in (41) are similar, it is sufficient to analyze one of them, as the conclusions for the other one would be the same. Let us consider the series with k ≥ 1. Then, we have to study the convergence of the series From (1) we have that with φ + 2kπ ≥ 1 , ∀k = 1, 2, . . . . For any real value of x, H n (x) is a real polynomial of order n with the property that for |x| > 1 has the following upper bound which is a straightforward consequence of the know formula of the Hermite polynomials, which states that Consequently, since π 1/4 2 n/2 √ n! ≥ 1, ∀n ∈ N, we have from (43) and (44) an upper bound for ψ n (φ + 2kπ) Hence, since n is fixed and after (46), the sum (42) converges if the following series converges. This convergence is a simple exercise of analysis. The conclusion is that the series in (41) with k ≥ 1 is absolutely convergent and hence pointwise convergent. In particular, this means that the functions C n are all Lebesgue measurable, since they are the pointwise limit of measurable functions. Similarly, the same property is valid for the series with k < −1 in (41). This shows the boundedness of C n (φ) on −π ≤ φ < π for n = 0, 1, 2, . . . , which, along with its measurability in the Lebesgue sense, shows that the functions C n (φ) are square integrable in the considered interval. Proposition 1 has an important consequence. We have shown that, for each n ∈ N, the function C n (φ) has an upper bound. The constant function that equals to this upper bound in the interval [−π, π) is square integrable, so that by using the Lebesgue dominated convergence theorem [32] we have that We also have that Finally, we may mention another interesting fact: the functions {C n (φ)} n∈N span L 2 [−π, π). This means that the subspace of all finite linear combinations of these functions is dense in L 2 [−π, π). The proof will be given later.
One of the objectives of the present article is to find some relations between functions that are of use in Fourier analysis. We shall discuss this idea along the next Section.

Fourier Transform on the Circle
After the general formula for the Fourier Transform given in (4), we obtain that Furthermore, the Inverse Fourier Transform (IFT) gives the following relation For the special case a = 2π, Equations (50) and (51) give, respectively, We may extend the previous formulae for a finite sum of Hermite functions, so that where D m (ay) is the Dirichlet kernel [33], i.e., From (48), we conclude that where In particular, for a = 2π (58) looks like The infinite sum X is called the Dirac comb [34][35][36] It converges in the weak sense as a distribution on the space, D, of all Schwartz functions with compact support on R endowed with the strict inductive limit topology. For the properties of the Fourier transform we can write that FT[C n (x), y] = i n X(y) ψ n (y) .
As is well known, the Dirac Comb is autoconjugate under Fourier transform. This is evident because and both the Dirichlet kernel and the r.h.s. of the first expression go to the Dirac Comb for k → ∞. Moreover, the functions {C n (φ)} n∈N can be split first into even and odd and then into four subspaces. Indeed

A Discretized Fourier Transform
To begin with, let us compare the space L 2 (C), also denoted as L 2 [−π, π), with the space l 2 (Z). As is well known, an orthonormal basis on where the sum converges in the sense of the norm (38) (for continuous functions f (φ) the series also converge pointwise [37]). The properties of orthonormal bases in Hilbert spaces show that We call to the complex numbers f m (m ∈ Z) the components of f . The Hilbert space l 2 (Z) is a space of sequences of complex numbers A ≡ {a m } m∈Z such that A This is a Hilbert space with a scalar product given by where the star denotes complex conjugation. An orthonormal basis for l 2 (Z) is given by the sequences that have all their components equal to zero except for one which is equal to one. Let us call {B m } m∈Z on this basis, where each of the B m represents each one of these series. Any F ∈ l 2 (Z) with components F ≡ { f m } m∈Z may be written as We readily see that there exists a correspondence between L 2 (C) and l 2 (Z). This correspondence relates any f (φ) ∈ L 2 (C) as in (63) with F = F/2π as in (67) with the same sequence { f m } m∈Z . This correspondence, f (φ) −→ F, is clearly linear, one to one and so on. In addition, f (φ) = F , which shows that it is, in addition, unitary (in fact any pair of infinite dimensional separable Hilbert spaces are unitarily equivalent in the sense that one may construct one, in fact infinite, unitary mappings from one to the other).
Equation (63) gives the Fourier series span for f (φ) ∈ L 2 (C). From this point of view, we may say that the Fourier series is a unitary mapping, F , from L 2 (C) onto l 2 (Z). It admits an inverse, F −1 , from l 2 (Z) onto L 2 (C), which is also unitary and is sometimes called the discrete Fourier transform [6].
We intend to offer a homogeneous version of concepts that are often introduced as separated. They are the Fourier transform, Fourier series, and discrete Fourier transform on one side and the Hermite functions on the other.
The first step is to construct a set of sequences in l 2 (Z) using the Hermite functions ψ n (x). For each n ∈ N, let us define the following sequence indexed by the set of integer numbers Z χ n := {ψ n (m)} m∈Z .
Our next result has not been proven so far. Its proof requires the previous demonstration of a couple of results displayed in Appendix A (Propositions A1 and A2).
Since the functions C n (φ) are in L 2 [−π, π), they admit a span in terms of the orthonormal basis in L 2 [−π, π) as commented upon at the beginning of the present Section (see Equation (63)). Thus, we can write with The continuity of the functions C n (φ) on [−π, π) guarantees the pointwise convergence of (69) [37]. In addition, since all C n (φ) are periodic with period 2π, (69) is valid for all real numbers of φ. We recall that each of the Hermite functions ψ n (x) are eigenfunctions of the Fourier Transform with eigenvalue i n . Let us use this idea in (70) in order to find an explicit expression of the coefficients c m n in terms of the values of the Hermite functions at the integers. Using the definition (39) of the C n (φ) in (70), we obtain The second identity in (71) makes use of the Lebesgue-dominated convergence theorem [32] in order to interchange the integral and the series as anticipated in (15). We have also used the change of variable s = φ + 2kπ and e imφ = e im(φ+2kπ) = e ims . This shows that (69) and (70) can be written, respectively, as and After (72) and (73), we may obtain a relation between the functions C n (φ) and the series χ n defined in (68) as C n (φ) ≡ i n χ n .
Thus, (72) and (73) provide a one to one relation between the functions C n (φ), as defined in (39) and the sequences χ n given in (68). Observe that the definition (39) of C n (φ) allows us to write (72) as So far, we have discussed the relation between a system of generators in L 2 [−π, π) ≡ L 2 (C) given by {C n (φ)} n∈N and a set of series {χ n } n∈N in l 2 (Z). These systems of generators do not form orthonormal bases on L 2 (C) and l 2 (Z), respectively.

Relevant Operators Acting on L 2 (C) and on l 2 (Z)
We want to discuss some of the properties of the functions {C n (φ)} in relation with its behavior under different operators. Here, we introduce a set of operators on L 2 (C) and on l 2 (Z) similar to the operators acting on the quantum harmonic oscillator. In particular, we have creation and annihilation and number operators that act on the chosen basis as expected. In addition, we have some other operators which play the role of multiplication by the variable and differentiation, with the expected relation with the ladder operators. In the present context, the definitions of these operators has to be done on a particular form. As also happens in relation to the harmonic oscillator, these operators are not bound on the Hilbert spaces they act on and have different domains. Nevertheless, we may equip these Hilbert spaces with a dual pair of locally convex spaces, so that these operators be continuous. This construction will be done in the next Section.
The possibility of introducing multiplication and differentiation operators on the circle has been previously considered and noticed that serious inconsistencies emerge when we try to extend these operators to the circle on its most natural form [38][39][40][41][42][43][44][45]. In particular, the need for boundary conditions for the wave functions of the form ϕ(0) = e iak ϕ(0) produce an ambiguity on the definition of the derivation operator, which now has to depend on the parameter k. In the formalism we introduced in the sequel, which in part has been based on the Weil-Brezin-Zak transformation [39,46,47], we try to avoid some of these inconveniences [38][39][40].

Multiplication and Derivation Operators on L 2 (C)
Let us give some definitions such as multiplication and derivation operators within our context.

Multiplication Operator
For the multiplication operator Φ, we have to discard the apparently most natural definition that, for any f (φ), it would have been Φ f (φ) = φ f (φ). Since f (φ) is a function on the circle that should be extended by periodicity to the real line, this definition is not appropriate as it does not provide a periodic function. Let us define the operator Φ by means of its action on each of the functions of the sequence {C n (φ)} and, then, extend it by linearity. Therefore, Φ would be defined on a dense set of the closed subspace spanned by the vectors of the sequence {C n (φ)}. Our definition is which is indeed periodic with period 2π. We may extend this definition for any real a ≥ 1 as since the series in the r.h.s. in (77) are absolutely convergent.
We want to study some properties of the operator Φ. We begin with the following property valid for Hermite polynomials Using this property in (76), we obtain Note that ΦC 0 (φ) ∝ C 1 (φ).

Derivative Operator
Next, let us define the derivative operator D φ on the subspace of all linear combinations of the functions {C n (φ)} n∈N . Clearly, we just need to define the action D φ on each of these functions. To begin with, let us write From the properties of the Hermite functions, we may obtain two different expressions, although equivalent, for the derivative of a Hermite function, which are Let us use the first expression of (80) in (81). The result is where the last identity on (82) comes from the absolute convergence of the series involved and the definition of Φ. Next, using the second expression of (81), we have We may use either one of the equivalent relations (82) or (83) as the definition of the operator D φ .

Ladder Operators
Let us define ladder operators on the subspace of linear combinations of the elements of {C n (φ)} as follows From (82)-(84), we obviously have that Exactly as with the quantum harmonic oscillator, we may define the following number operator It is also obvious that This completes the analogy with the quantum harmonic oscillator. All these operators admit closed extensions and are unbound.

Operators on l 2 (Z)
In the previous two subsections we have been concerned with functions in and operators on subspaces of L 2 (C). Now, let us find the equivalent objects in l 2 (Z).
Using the canonical orthonormal basis {B m } m∈Z of l 2 (Z) defined in Section 5, we can write after (68) a being any real number greater than or equal to 1. Expressions (90) are well defined and belong to l 2 (Z). In fact, after (67) and some analysis, we conclude that In addition, we may also define a formal derivative on the space of finite linear combinations of the {χ n } as follows Choosing the first expression in (81) and taking into account (91), we have where the first identity in (93) makes sense because {B m } is an orthonormal basis. Analogously, using the second row in (82), we have Both Equations (93) and (94) show independently that D is well defined in the space of finite linear combinations of vectors {χ n } and, therefore, (92) makes sense. In addition, if we define on the same space the creation, B + , and annihilation, B − , operators as we have the following relations and Formulas (95)-(97) are equivalent to (84)-(86), respectively. Then, we may define the corresponding number operator as which on the considered systems of generators gives These expressions give a harmonic oscillator like equation Equations (99) and (100) are valid for all n = 0, 1, 2, . . . . Again, these operators are closable and unbound.

On the Continuity of the Relevant Operators
Let us go back to the sequence of functions {C n (φ)}. These functions are linearly independent after Proposition A2 (Appendix A). Thus, the functions in the sequence {C n (φ)} are linearly independent. Then, we may consider the linear space G spanned by them and introduce on it a new scalar product defined as C n |C m = δ nm , where δ nm is the Kronecker delta. This scalar product is now extended to the whole G by linearity to the right and anti-linearity to the left. The resulting pre-Hilbert space may then be completed so as to obtain a Hilbert space that we shall denote as H. Next, let us define the space S of all functions f (φ) = ∞ ∑ n=0 a n C n (φ) ∈ H , a n ∈ C with the following property f (φ) 2 p := ∞ ∑ n=0 a n 2 (n + 1) 2p < ∞ , p = 0, 1, 2, . . . .
Clearly, S is isomorphic algebraic and topologically to the Schwartz space S of all C ∞ (R) functions that go to zero at infinity faster than the inverse of any polynomial, see [2], as the topology on S is given by the countable set of norms (102). The triplet S ⊂ H ⊂ S × , where S × is the dual space of S endowed with the weak topology corresponding to the dual pair {S, S × } [48] is a rigged Hilbert space or Gelfand triplet. The operators Φ, D φ , and A ± φ are continuous linear operators on S and the same property holds for the algebra spanned by these operators. For instance, let us pick f (φ) ∈ S (101). Then, using (79), we have This proves that Φ S ⊂ S with continuity. From here, it is obvious that the same property holds for Φ m , with m = 0, 1, 2, . . . . The proof for the same property concerning D φ comes from (83) and for A ± φ from (84). From these results, it is clear that all the operators in the algebra spanned by Φ, D φ and A ± φ are continuous on S, including the number operator N φ . Then, observe that A ± φ are the formal adjoint of each other. From (94), we see that D φ is formally symmetric and that the formal adjoint of Φ is −Φ. Let B be an arbitrary densely defined operator on the Hilbert space and B † its adjoint. Assume that B † leaves S invariant, which means that B † f ∈ S for any f ∈ S. Then, using the duality formula one shows that B may be extended as a linear operator to the dual S × . In addition, if B † is continuous on S, so is B on S × with the weak topology on S × . Therefore, the algebra spanned by the operators Φ, D φ , and A ± φ may be extended to the dual S × and these extensions are continuous with the weak topology on the dual. However, the above discussion, notwithstanding its simplicity, relies on unnatural Hilbert metrics and is somehow artificial. On the other hand, the use of the natural scalar product on L 2 (C) may require of the introduction of a more artificial test space that we shall denote with D. In order to construct D, let us consider the linear space of all finite linear combinations of functions C n (φ) On this linear space, we define the following set of seminorms (indeed norms), p k (−) The resulting locally convex space needs not be complete, although it is always possible to complete it with respect to the locally convex topology generated by the semi-norms (106). We call D to this completion. Note that where the identity in (107) is a consequence of (A16) and (A17). Then, combining inequalities ∑ ∞ m=−∞ e −m 2 m 2n ≤ c(n + 1)! 2 n+1 first and then (A6) , we obtain that the last term in (107) is smaller or equal to This chain of inequalities shows that the canonical identity i : is a rigged Hilbert space, or Gelfand triplet. The dual space D × is endowed with any topology compatible with duality (strong, weak, McKey). Now, proving the continuity of the operators defined in Sections 6.1 and 6.2 on D is rather trivial. For instance, take A + φ as defined in (84). Obviously, Hence, so that which proves, both, that A + φ D ⊂ D and that A + φ are continuous on D. Similar proofs apply to the other operators in (76) and (77).
Analogous results can be obtained when dealing with the operators defined in Section 6.3.

Concluding Remarks
We investigated the role of Hermite functions in Harmonic analysis in connection with Fourier analysis. We showed that Hermite functions permit the construction of a complete set of periodic functions defined in the unit circle that span L 2 (C). Using the Gramm-Schmidt procedure, we readily obtain an orthonormal basis for L 2 (C) out of these functions.
At the same time, and using the normalized Hermite functions, we constructed a system of generators in l 2 (Z), the space of square summable complex sequences indexed by the integer numbers. We showed that the use of Fourier series and Fourier transform relates both systems of generators, in L 2 (C) and l 2 (Z) in a very natural way, defining a unitary transformation between these two spaces.
On the subspace of L 2 (C), including a complete set of periodic functions, we defined a multiplication and a derivation operator that preserve periodicity in both cases. These operators generate creation and annihilation operators for the defined complete set of periodic functions, which behave just as creation and annihilation operators for the harmonic oscillator. Similar operators with identical properties are defined for l 2 (Z). We have constructed rigged Hilbert spaces supporting these operators on which they are continuous operators.

Funding:
We acknowledge partial financial support to the Junta de Castilla y León, grants BU229P18 and VA057U16.
Institutional Review Board Statement: Not applicable.

Informed Consent Statement: Not applicable.
Data Availability Statement: Not applicable.

Conflicts of Interest:
The authors declare no conflict of interest.

Appendix A. Determinants of Hermite Polynomials
In Appendix A, we shall prove some results relevant for the development of Section 5.
where H n (k) is the Hermite polynomial H n (x) evaluated at the integer point k and m ∈ N.
Proof. Let us prove that the second determinant is different from zero. The proof of this property for the first determinant is similar. We proceed by induction on the dimension m of the determinant. Since we know explicit expressions for the low order Hermite polynomials, we easily conclude that the determinant on the right hand side of (A1) is different from zero for m = 0, 1, 2, 3, 4. Let us assume that this property is true for m = 0, 1, 2, . . . , n. Take m = n + 1 and assume that the To prove that the sequences are linearly independent, let us consider an arbitrary finite linear combination (by a standard definition in linear algebra, all linear combinations are finite) of the form Since {B m } m∈Z is an orthonormal basis in l 2 (Z) the change in the order in (A8) is legitimate as well as that the p ∑ n=0 a n ψ n (m) = 0 , expression that must be valid for m = 0, ±1, ±2, . . . . If we choose m = 0, 1, 2, . . . , n, (A9) gives a linear homogeneous system of p + 1 equations with p + 1 indeterminates, which are a 0 , a 1 , . . . , a p . The determinant of the coefficients of this system is given by This determinant is different from zero if and only if the second determinant in (A1) is non-vanishing, which is true due to Proposition A1. Therefore, we have that a 0 = a 1 = · · · = a p = 0 and the {χ n } n∈N are linearly independent. Next, let us show that the set of sequences {χ n } n∈N spans l 2 (Z). We know that this happens if and only if for all n ∈ N, we have [5] for any ϕ in a dense subspace of l 2 (Z). As a dense subspace, we choose the space of all finite linear combinations of elements of the basis {B m }. These vectors are of the form where p is an arbitrary non-negative integer and the coefficients a m are arbitrary. Let us assume that, for all χ n , the scalar product χ n |ϕ p = 0. We have that This second determinant is just the determinant in (A1) as all the exponentials cancel out. Then, and due to Proposition A1, the determinant is non-vanishing and consequently, the indeterminates a −p = · · · = a 0 = · · · = a p = 0, so that ϕ n ≡ 0 and the set {χ n } n∈N is complete in l 2 (Z).

Appendix B. Orthonormal Systems in L 2 (C) and l 2 (Z)
In order to construct an orthonormal basis in L 2 (C) after the set of functions {C n (φ)}, we first evaluate the scalar product on L 2 (C) of two of these functions Note that after (A16) Now, we are in the position of using the system {C n (φ)} to construct an orthonormal system in L 2 (C). To do it, we need the Gramm-Schmidt process [49] (see, for instance, Ref. [50] for a Gramm-Schmidt algorithm in order to use computer techniques). The resulting orthonormal system would be an orthonormal basis for L 2 (C) as a consequence of Proposition A1 (see Appendix A). As is well known from the set {C n (φ)}, we can obtain firstly a set of orthogonal functions {D n (φ)} where However, taking into account (A18), we can rewrite the relations (A19) as follows distinguishing between even and odd subindices of the functions, i.e., (A20) Let us denote as C n (φ) the normalized functions Due to the properties of the Gramm-Schmidt process, the set of functions { C n (φ)} n∈N forms an orthonormal system in L 2 (C). Analogously, and using Formulae (74) and (A20), we may construct new sequences, {ζ n }, in l 2 (Z) defined as as well as the normalized series χ n := ζ n ζ n |ζ n , n ∈ N . (A23) Series (A23) gives an orthonormal system in l 2 (Z). We denote the components of the χ n as χ n = { ψ n (m)} m∈Z . After relations (74), (A20) and (A22), it is easy to show that (72) and (73) yield, respectively, and Thus, we have a partial isometry U between L 2 (C) and l 2 (Z) given for any f (φ) ∈ L 2 (C) in the closed subspace spanned by the { C n (φ)} by After Proposition A1 (Appendix A), we conclude that U is a unitary mapping between L 2 (C) and l 2 (Z).