Dynamical Symmetries of the 2D Newtonian Free Fall Problem Revisited

: Among the few exactly solvable problems in theoretical physics, the 2D (two-dimensional) Newtonian free fall problem in Euclidean space is perhaps the least known as compared to the harmonic oscillator or the Kepler–Coulomb problems. The aim of this article is to revisit this problem at the classical level as well as the quantum level, with a focus on its dynamical symmetries. We show how these dynamical symmetries arise as a special limit of the dynamical symmetries of the Kepler–Coulomb problem, and how a connection to the quartic anharmonic oscillator problem, a long-standing unsolved problem in quantum mechanics, can be established. To this end, we construct the Hilbert space of states with free boundary conditions as a space of square integrable functions that have a special functional integral representation. In this functional space, the free fall dynamical symmetry algebra is shown to be isomorphic to the so-called Klink’s algebra of the quantum quartic anharmonic oscillator problem. Furthermore, this connection entails a remarkable integral identity for the quantum quartic anharmonic oscillator eigenfunctions, which implies that these eigenfunctions are in fact zonal functions of an underlying symmetry group representation. Thus, an appropriate representation theory for the 2D Newtonian free fall quantum symmetry group may potentially open the way to exactly solving the difﬁcult quantization problem of the quartic anharmonic oscillator. Finally, the initial value problem of the acoustic Klein–Gordon equation for wave propagation in a sound duct with a varying circular section is solved as an illustration of the techniques developed here.


Introduction
In an elementary freshman physics course, the problem of particle motion under constant gravitational acceleration on the surface of the earth in Euclidean space, usually considered as one dimensional, is frequently dubbed as a free fall problem [1,2]. This is not the more general relativistic problem in which an inertial particle subject to no force moves along a space-time geodesic [3].
To the best knowledge of the author, up to now, no comprehensive report on the dynamical symmetries of this 2D Newtonian free fall problem has been found in the literature besides the seminal work of T. Iwai and S. G. Rew [4]. As the one-dimensional free fall problem, which has been thoroughly treated classically and quantum mechanically, displays limited interesting features, we introduce an extra space dimension to allow for new specific features to emerge. This is how a body of dynamical symmetries comes about, and how a connection to the quartic anharmonic oscillator problem is generated, which is one of the most challenging theoretical problems as it is generally thought of as a non-trivial quantum field theory in zero space dimension. These issues were not considered in [4]. This paper is divided into two main sections and a short section on an application. Section 2 is devoted to the classical physics of the 2D Newtonian free fall problem. For the convenience of the reader, some basic concepts on the symmetries of a system are recalled in the framework of Hamiltonian mechanics before a derivation of the dynamical symmetry algebra is undertaken along the lines of [4]. Then, we prove a theorem which states that the dynamical symmetry algebra of the free fall problem may be obtained as a special limit of the the dynamical symmetry algebra of the Kepler-Coulomb problem, which concerns the dynamics of a particle moving in the inverse separation distance potential field. Both problems are known to be super-integrable, and now they are shown to be connected at the classical level. Hence, one may access the free fall dynamical symmetry algebra from the well-known Kepler-Coulomb dynamical symmetry algebra. In addition, this connection goes even further when a higher order integral of motion in the free fall problem is shown to be an outgrowth of a Kepler-Coulomb integral of motion at zero coupling constant. Furthermore, a passage to parabolic coordinates reveals a connection between the free fall problem and the quartic anharmonic oscillator.
In Section 3, we apply the Schrödinger quantization to the system, which is in fact not the quantization scheme adopted by T. Iwai and S. G. Rew in [4]. For this, we introduced a Hilbert space of square integrable wave functions that verify free boundary conditions in R 2 and are suitable for obtaining an image representation of the free fall dynamical symmetry algebra as Klink's so-called one-variable algebra of the quartic anharmonic oscillator. Thus, the Schrödinger quantization does naturally lead to the quartic anharmonic oscillator besides the passage to parabolic coordinates. This is due to a particular integral representation of the free fall wave functions, which takes the form of an Airy-Fourier transform. Finally, by performing a passage to parabolic coordinates, one obtains an integral identity for the eigenfunctions of the quartic anharmonic oscillator problem, which turns out to be the characteristic integral equation for zonal functions in the representation of some underlying group. As the only natural group arising from the 2D free fall problem is its own dynamical symmetry group, we suspect that a comprehensive development of its representation theory would reveal that the eigenfunctions of the quartic anharmonic oscillator problem are just the corresponding zonal functions of this group. We defer this tantalizing investigation to a future work. In short, the present results on quantizing the 2D Newtonian free fall dynamical symmetry algebra appear as a necessary intermediate step in the search for an exact solution to the quantum quartic anharmonic oscillator problem.
The last section, Section 4, uses the functional techniques of Section 3 to solve the initial problem of sound wave propagation in a cylindrical duct with a particular varying circular section.
The paper ends with a short conclusion and perspectives on possible further research.

Generalities on Classical Symmetries, Canonical Structure, and Integrals of Motion
As generally admitted, the classical state of a non-moving physical system is described by a set of functions in a coordinate system. Under a coordinate transformation, these functions may take different forms. But if they remain invariant, the system is said to admit a symmetry. Now, for a physical system in motion, its states are specified by functions of time t and dynamical coordinates in phase space. If under a transformation of dynamical coordinates its state functions remain invariant, we say that the system admits a dynamical symmetry. The set of symmetry transformations may have a group structure, which is called dynamical symmetry group. Such a dynamical symmetry group (or equivalently, the algebra of its generators) is of highest interest in the search for solutions of its equation of motion (The term dynamical symmetry was coined by A. O. Barut [5]. Earlier, such symmetry was known as hidden or accidental [6]).
In this paper, we are concerned with the dynamical symmetries of a system in the framework of Hamilton's canonical formalism of mechanics, with the time-independent Hamiltonian function H. We now recall some useful main points of this framework.
-A system with n degrees of freedom is described by 2n canonical variables (q i , p i ) with i = 1, . . . , n in phase space, verifying the following fundamental Poisson bracket commutation relations: where the Poisson bracket between two functions in phase space φ and ψ is defined by the following: and the canonical equations of motion: the solutions of which depend on the 2n initial conditions (q i (t), p i (t))| t=0 = (q 0 i , p 0 i ). Consequently a dynamical symmetry exists if it originates from a canonical coordinate transformation in phase space, leaving both the Hamiltonian H and the fundamental Poisson bracket relations invariant. In differential geometry, this is called a symplectic structure. Non-trivial dynamical symmetry exists only for n > 1.
This can be equivalently expressed by requiring the fundamental Poisson brackets in the variables (Q 1 , . . . , Q n , P 1 , . . . , P n ) to be valid: The Poisson bracket of two functions (φ, ψ) is invariant under the contact transformations {φ, ψ} (Q,P) = {φ, ψ} (q,p) . The lower index (q, p) (resp. (Q, P)) refers to the variables in the partial derivatives of the respective Poisson brackets. A contact transformation (q 1 , . . . , q n , p 1 , . . . , p n ) → (Q 1 , . . . , Q n , P 1 , . . . , P n ), which preserves the equations of motion of the system is a canonical transformation if ∑ n i=1 P i dQ i is an integral invariant of the system (Jacobi's theorem). Consequently, a canonical transformation implements a dynamical symmetry.
An infinitesimal canonical transformation of parameter τ has the following form: where F τ is an arbitrary function of (q 1 , . . . , q n , p 1 , . . . , p n ) and ∆τ is a small increment of the parameter τ. F τ is called the the generating function of the contact transformation.
The variation ∆ f of an arbitrary function f in phase space under an infinitesimal contact transformation of generator F τ is given by ∆ f = { f , F τ }∆τ. F τ can be an integral of the motion, a function which takes a constant value and does not explicitly depend on time t, see [7]. A system of integrals of motion F i is said to be in involution when {F i , F j } = 0 for all (i, j = 1, . . . , m). If F and F are two integrals of the motion, their Poisson bracket is also an integral of the motion (Poisson's theorem). The set of all F j forms a Lie algebra with respect to the Poisson bracket. The determination of all canonical transformations for a dynamical system is at the core of finding the dynamical symmetries of this system. -For a system with an n degree of freedom, with a time-independent Hamiltonian H, total energy is conserved and H = E. If there exists n functionally independent integrals of motion F i , where i = 1, . . . , n in involution (or {F i , F j } = 0 for all (i, j)), the system is said to be Liouville integrable and its solution can be given up to quadratures. If there are further m functionally independent integrals of motion F m , where m = 1, . . . , (n − 1), the system is called super-integrable [8]. For m = 1 (resp. m = (n − 1)) it is called minimal super-integrable (resp. maximal super-integrable). These m extra integrals of motion usually build a Lie algebra. Maximal super-integrable systems are also known as exactly solvable systems, and their properties can be derived algebraically. Each integral of motion F i may originate from Noether's conservation law or from a coordinate variable separation. The Kepler-Coulomb (or inverse distance potential) and the isotropic harmonic oscillator problems are known to be super-integrable systems in two dimensions [9].

Statement of the Classical Two-Dimensional Free Fall Problem
As it is widely known, the motion of a particle under constant force in two-dimensional space occurs along parabolas with a symmetry axis parallel to the direction of the constant force. To simplify the writing, the particle is assumed to have a unit mass and moves in a two-dimensional configuration space (v, u) ∈ R 2 , with Ov being the horizontal axis and the gravitational constant set is also equal to one. Let (v,u) be the time derivatives of the Cartesian coordinates (v, u). From the Lagrangian L = 1 2 (u 2 +v 2 ) − u , one gets the conjugate momenta and the Hamiltonian: Note that H is the sum of the v-free motion part and of the u-free fall part.

Dynamical Symmetries
As in two dimensions, the maximal number of possible symmetries is three. Iwai and Rew [4] were the first to obtain dynamical symmetries by considering linear inhomogeneous transformations S(ρ, σ, τ) in the phase space dependent on three parameters (ρ, σ, τ): Proposition 1. S(ρ, σ, τ) is a canonical transformation.
Proof. By the simple substitution of the (V, U, P v , P u ) in terms of the (v, u, p v , p u ), as given by Equation (6), it appears that the Hamiltonian function form is invariant: Moreover, using the expression of the Poisson bracket and Equation (6), one can also check the following: Proof. Composition law. Let another canonical transform S(ρ , σ , τ ) with parameters (ρ , σ , τ ) act on the previous one S(ρ, σ, τ) and compute the resulting product S(ρ , σ , τ ) S(ρ, σ, τ), which is expressed by the set of new parameters (V, U, P v , P v ) : Then, the substitution of Equation (6) into Equation (9) yields the following equations: Hence, we conclude that the composition of two operations is as follows: Because of the extra term (στ − σ τ), the previous formula does not correspond to the additive structure of a group operation with respect to the parameters (ρ, σ, τ). But since the v variable is a cyclic variable, this term is physically irrelevant ; in this way, the group structure is physically restored.
Iwai and Rew proposed to represent each element of the group with fifth-order, upper-triangular matrices, as can be checked by the matrix multiplication rule: -The identity matrix is obviously the neutral group element.

Infinitesimal Iwai-Rew Canonical Linear Transforms and Integrals of Motion
The infinitesimal form of the Iwai-Rew canonical linear transformations is written in the form of Equation (8), with infinitesimal (∆ρ, ∆σ, ∆τ) as: Proof. By working out the Poisson brackets using Equation (12), they will appear to remain in the canonical form of Equation (8), while the form of the Hamiltonian in the variables (V, U, P v , P u ) is the same as the form in variables (v, u, p v , p u ).
We now determine the three generating functions (F ρ , F σ , F τ ) from the three subgroups of canonical transformations, respectively defined by the following equations: Then, partial derivatives of (F ρ , F σ , F τ ) with respect to (v, u, p v , p u ) can be deduced as follows: Therefore, we obtain the exact differentials (dF ρ , dF σ , dF τ ) in terms of (dv, du, dp v , dp u ) (Schwarz's theorem is trivially verified for all pairs of variables). Their integration yields the sought integrals of motion: Proposition 4. The free fall system in two dimensions is super-integrable since its has three integrals of the motion: Proof. Compute the three Poisson brackets and observe that they are zero.

Proposition 5.
The three generating functions (F ρ , F σ , F τ ) build a Lie algebra structure, the dynamical symmetry algebra A FF of the two-dimensional free fall with respect to the Poisson bracket: Proof. Work out the Poisson brackets and use the expressions of the (F ρ , F σ , F τ ) from Equation (15). They are identical to those of Iwai-Rew [4]. If we make the substitutions (15), we recover the Iwai-Rew commutation relations given by their Equation (2.21).

Remark 1.
Observe that the classical trajectory data (two initial position coordinates, two initial momentum coordinates) can be used to compute the values of H, F ρ , F σ , F τ and vice versa.

Proposition 6.
Origins of the integrals of motion (a) F ρ is due to translational invariance of H in the Ov direction; (b) F σ is due the separability of H in a π/4-rotated coordinate system; (c) F τ is due to the manifest separability of H in v and u variables.

Proof. (a) For this case, the proof is trivial because
. Substitution into the expressions of H and F σ yields the following equation: where H and F σ are both separable in the new coordinate system obtained by a π/4 rotation of the (v, u) coordinate system around the origin. In fact, they are the sum and difference of the one-dimensional free fall Hamiltonian in the ξ and ζ directions. Hence, {F σ , H} = 0.
In reference [10], it was claimed that at the quantum level, F σ is due to separability in translated parabolic coordinates. But so far, no proof has appeared in print.
(c) As H is the sum of a free-motion Hamiltonian in the v-direction, and the Hamiltonian is that of a one-dimensional free fall in the u-direction, it is then clear that F τ , which is the difference of these two Hamiltonians, should verify {F τ , H} = 0. Since the total energy is conserved as H = E, it may be simpler to write F τ = (p 2 v − E).

Remark 2.
We notice that v-parity is also a dynamical symmetry and a discrete canonical transformation since v → −v (which also implies p v → −p v ) leaves H as well as {v, p v } = 1 and all other Poisson brackets invariant.

The Free Fall Problem as a Special Limiting Case of the Kepler-Coulomb Problem
In this section, we show that the free fall problem can emerge from a special limit of the Kepler-Coulomb (KC) problem or the problem of the inverse distance potential, occurring either in gravitational or electrical interaction. This problem is known to be super-integrable and in two dimensions, and the set of its three integrals of motion makes up the components of the so-called Runge-Lenz vector [11]. Here, we adopt the notations of [12], with the particle mass set equal to one.
The idea is simple. The inverse distance potential κ/r arising from a source of strength κ is rotation-symmetric, where r is the distance from the source to the observation site. If the source recedes to infinity in the u-direction, the potential at the observation site tends toward zero. But one may compensate this potential decrease by taking a source strength κ(r), which increases with the distance. Thus, we may choose an increasing functional dependence so that in the limit of infinite source-observation separation, a linear potential appears.
In our Cartesian coordinate system, if the source of strength κ is placed at the coordinate origin, the Hamiltonian of the inverse distance problem is the following: and the three integrals of motion are given by [12]: where L w is the angular momentum around Ow, which is orthogonal to the plane Ovu. They verify the Poisson bracket relations of the dynamical symmetry algebra A KC Now, if the source is no longer at the origin O but situated on the Ou axis at a distance l from the origin, the Hamiltonian and the components of the Runge-Lenz vector have new expressions H KC and (L w , K v , K u ), which are deduced from the previous expressions in which u is replaced by (u + l) and κ by κ(l): and Now, as l → ∞, the asymptotic behaviour of the inverse distance potential is as follows: Hence, if the source strength increases as κ(l) = l 2 , then the inverse distance potential reaches the limiting form of a linear potential in u: In this limit, the inverse distance potential problem tends toward the free fall problem up to a negative infinite constant: and the components of the Runge-Lenz vector take the following asymptotic forms: Theorem 1. The dynamical symmetry algebra of the free fall problem A FF is a contraction of the dynamical symmetry algebra of the Kepler-Coulomb problem A KC .
Proof. We now rewrite the Poisson brackets of the Kepler-Coulomb problem when the potential source is at a large distance l from the origin, and then replace the generators (L w , K v , K u ) by their asymptotic expansions for l → ∞. Therefore: Extracting the leading order in l 2 on the left-hand side and on the right-hand side, we obtain {F ρ , F σ } = −I, as expected; see Equation (17).
Collecting terms of leading order in l 2 on both sides of this equation, we get {F τ , F ρ } = 0 since there is no term in l 2 on the right-hand side.
Equating terms of order d 2 on both sides of this equation yields precisely {F σ , F τ } = 2F ρ . Hence, we reproduce all the Poisson brackets of A FF .

A "Higher" Order Integral of the Motion
We now raise the question whether there exists a "higher" order integral of the motion as a construct of the dynamical symmetry algebra generators. What comes to mind is a weighted sum of squares of the (F ρ , F σ , F τ ), an object similar to the square of the angular momentum in a rotation algebra. Instead of a tedious systematic search, a more astute way of finding such an integral of the motion would start by observing that the symmetry algebra of the Kepler-Coulomb problem does have a non-trivial limit when the source strength is turned off, i.e., κ = 0. Then, the generators take the following form: They fulfil the same Poisson bracket relations as those with κ = 0: where H 0 is the Hamiltonian for free particle motion in two dimensions. Now, if we allow a linear potential to "grow" in the u-direction, rotational symmetry disappears. Hence, L w must be discarded as a possible integral of motion in the presence of a linear potential u. Thus, from the two remaining (K v (0), K u (0)), only one can survive under a modified form as a "higher" order integral of the motion, because of Poisson 's theorem in Section 2.1.
Let F 2 be this hypothetical "higher" integral of the motion. As the introduced linear potential is in the u direction, we may assume F 2 to be of the simple form: where h(v, u, p v , p u ) is an unknown function in phase space.

Proposition 7.
There exists a second order integral of the motion F 2 given by the following equation: Proof. The function h(v, u, p v , p u ) is determined by the condition F 2 , H = 0. Since from explicit computation one gets the following: it is obvious that one should require that h(v, u, p v , p u ) = − 1 2 v 2 .
Remark 3. The same search procedure for another quadratic integral of the motion does not work with K v (0) because the commutativity of the Poisson bracket with H leads to impossible conditions having to be satisfied.

Proposition 8.
The expression of F 2 in terms of the generators (F r , F s , F t ) is as follows: Proof. Substitute the expressions (H, F ρ , F σ , F τ ) into the expression of F 2 .
Corollary 1. The Poisson brackets of F 2 with the generators (F ρ , F σ , F τ ) are easily obtained: Proposition 9. The infinitesimal transformation generated by F 2 with parameter z is canonical.
Proof. The infinitesimal canonical transform generated by F 2 with the parameter z is given by the following equations: We can check that the Hamiltonian remains invariant if the terms of order (∆z) 2 are ignored: Next, we can verify that the six canonical Poisson brackets are preserved, but the details are not presented here. Moreover, it does possess the additive abelian group property with respect to z, i.e., S(z) · S(z ) = S(z + z ), as can be checked explicitly.

Passage to Parabolic Coordinates and the Physical Meaning of F 2
The issue here is to understand how F 2 arises as a dynamical symmetry. The comprehensive work of Miller et al. [13] on quantum separability has revealed that parabolic coordinates do play a central role. Following this indication, we make a passage to parabolic coordinates (x, y) from our Cartesian coordinates (v, u), as defined by u = 1 2 (x 2 − y 2 ), v = xy. The Lagrangian L then changes to a new expression: From the following definitions of conjugate momenta: we deduce a relation between the (v, u) and the (x, y) momenta: This allows for the acquisition of new expressions of H and F 2 in parabolic coordinates: As such, these expressions do not show any obvious x and y variable separation. The three integrals of motion (F ρ , F σ , F τ ) also do not display any obvious separation into an x part and a y part when re-expressed in the parabolic coordinates. Next, we recall that the initial conditions (v 0 , u 0 , p 0 v , p 0 v ) fully determine the four integrals of motion (H, F ρ , F σ , F τ ). Hence, F 2 has a fixed value because it is a construct of these four integrals of motion.

Proposition 10.
For H = E, F 2 takes a constant value in the range of energy values E + (E) of a confining quartic oscillator with an angular frequency square (−E).
Proof. As total energy is conserved, H = E implies that the following relation must be verified for all (x, y) = (0, 0): This means that for a given E, the sum of the Hamiltonians of a confining quartic xoscillator and a non-confining quartic y-oscillator must be equal to zero for all (x, y) = (0, 0). Since these one dimensional quartic oscillators are time-independent, their respective Hamiltonians have fixed values, i.e.: however, subjected to the condition E + (E) + E − (E) = 0. Note that for a given E, E + (E) takes all real values above the minimum value of the quartic x-oscillator Hamiltonian polynomial in phase space. On the other hand, F 2 may be rewritten as follows: Since E + (E) + E − (E) = 0, one gets F 2 = E + (E). This is due to this new aspect of separation of variables called the Stäckel separation of variables.

Third Order Integrals of the Motion
In [13], it is shown that third order super-integrable systems separable in parabolic coordinates admit second order integrals, while third order integrals are reducible, i.e., they are Poisson brackets of second order integrals of motion. Here, we have: However, these third order reducible integrals of motion do not generate one parameter infinitesimal canonical transformations because they do not leave H functionally invariant and are consequently uninteresting.

Quantization
The Hamiltonian form of the free fall problem lends itself nicely to its quantization. The classical canonical variables are put into one-to-one correspondence with their quantum counterparts, as self-adjoint operators in a Hilbert space of states H: They build a direct product of v and u Heisenberg algebras (hereh = 1 for ease of writing): where The quantum Hamiltonian, the dynamical symmetry algebra generators, and the quadratic integral of motion are then given by the following equation: They fulfil the quantum commutation relations deduced from their Poisson brackets classical counterparts: where I is the identity operator. The quadratic integral of motion F 2 is as follows:

Schrödinger Representation
It is convenient to work with the Schrödinger coordinate representation. The Hilbert space H of states |ψ ∈ H is represented by the square integrable functions ψ(v, u) = v, u|ψ , where |v, u for (v, u) ∈ R 2 is a continuous set of complete and total set in H. The canonical dynamical operators by differential operators in (v, u) are as follows: Consequently, the differential operator representation of the Hamiltonian, the dynamical symmetry generators, and the quadratic integral of motion are as follows: All Poisson brackets from the classical (H, F ρ , F σ, F τ, , F 2 ) now become quantum commutators between ( H, F ρ , F σ , F τ , F 2 ), as can be easily checked. Let A FF be the algebra generated by ( H, F ρ , F σ , F τ ). (56) Proof. Use the expressions in Equation (55) and substitute in (56).
As usual, the space of relevant wave functions L 2 (v, u) are generated by the eigenfunctions of the stationary Hamiltonian operator H. At this step, to determine this functional space, boundary conditions should be specified. As pointed out in the previous section on the classical mechanics of free fall motion, we are concerned with global motion along parabolic trajectories in R 2 , with concavity turned downward and not with the bouncing motion on a horizontal line v =constant or the billiard motion inside a two-dimensional box. Both have overly restrictive boundary conditions, rending the wave functions uninteresting. Therefore, only stationary wave functions with free boundary conditions on the v-axis and on the negative u-axis solutions of H ψ E (v, u) = E ψ E (v, u) are considered here.
Since H is separable in the v and u parts, elementary solutions are of the product form: where γ 3 = 2 and Ai(x) is the first Airy function which decreases at x → ∞; see [14]. Hence, the spectrum of H is real and continuous. An arbitrary stationary eigenfunction ψ E (u, v) with an eigenvalue E is given by an integral on k: where ψ E (k), yet to be determined, appears as a density amplitude representing the relative distribution of the v-motion with respect to the u-motion parts for a fixed E in ψ E (u, v).
Both ψ E (u, v) and ψ E (k) describe two aspects of the same quantum state |ψ E .
Proof. Let us compute the overlap integral ψ E |ψ E between two eigenstates of energies E and E : The 2π factor is just for convenience. Integration over dv yields 2π δ(k − k ). Then, after k -integration, one may consider the remaining u-integration.
This u-integral is readily given in [14] by equation 3.108 on page 57 as follows: Hence, we have a relation between the inner product in L 2 (v, u) and the inner product in L 2 (k): It expresses the completeness of the eigenfunctions ψ E (v, u), provided that: Hence, ψ E (k) ∈ L 2 (k).
Proof. By integrating both sides of Equation (58) with the following equation: we get the following integral on the right-hand side: Assuming Fubini's theorem hypothesis, we can exchange integrations and the dv 2πintegration would yield δ(k − k ). Then, after integration on dk, the right-hand side becomes: Ai γu + k 2 which, upon application of the integration formula (60), yields the following: A last integration on dE on both sides of this equation gives the final result: Hence, the integral mapping ψ E (v, u) → ψ E (k) is invertible, provided that a summation on E is performed or uses complete data. In this respect, it may be called the Fourier-Airy Transform. These two amplitudes describe the same physics in two different contexts.

Representation of the Dynamical Symmetry Algebra A FF by Klink's Algebra
Proposition 12. The action of the A FF generators on the wave function ψ E (v, u) is easily transferred to the wave function ψ E (k) according to the following: Proof. The action of H is E since ψ E (v, u) is an eigenfunction of H.
As F ρ and F τ are represented in L 2 (v, u) by v-derivatives acting under the integral sign on e ikv , we successively get k and k 2 acting on ψ E (k).
For the action of Since the Airy function vanishes for k = ±∞, we perform a partial integration in k to get the following: v where Ai'(x) is the derivative of Ai(x). We then observe that the second integral is as follows: Hence, the action of F σ on ψ E (v, u) is replaced by i ∂ ∂k ψ E (k) in L 2 (k).

Corollary 3.
As a consequence, the action of F 2 on ψ E (v, u) is translated into action on ψ E (k) as the action of the one-dimensional confining quartic anharmonic oscillator Hamiltonian in the k variable on the wave function ψ E (k): Proof. Use the expression of F 2 in terms of the following dynamical symmetry algebra generators ( F ρ , F σ , F τ ): and the previous proposition.

Corollary 4.
The mapping ψ E (v, u) → ψ E (k) induces an isomorphism between the two-variable algebra generated by ( F ρ , F σ , F τ ), and the one-variable algebra generated by (k, i d dk , k 2 ), which is known as Klink's algebra for quartic anharmonic oscillator and for which the quadratic integral of motion F 2 is the Schrödinger Hamiltonian of this quartic anharmonic oscillator [15].
Proof. The commutators of ( F ρ , F σ , F τ ) are isomorphic to those of (k, i d dk , k 2 ) by computational checking.
This relation shows how the one-dimensional confining quartic anharmonic oscillator is linked to the two-dimensional free fall. Note that this relation cannot be established if one had started with the one-dimensional free fall problem. In classical physics, the presence of quartic anharmonic oscillators in the free fall problem arises only when parabolic coordinates are introduced. This is not the case here. There is still another curious relation between the two systems discovered by Voros in [16].

Finite Quantum Canonical Transforms
Proposition 13. The finite quantum unitary transforms generated by ( F ρ , F σ , F τ ) are expressed as follows: They produce the quantum version of the classical Iwai-Rew transform (see Equation (8)).
Proof. The proof is trivial and involves the use of the Baker-Hausdorff-Campbell formula: where ( X, Y) are operators. Because of the commutation relations of A FF , the computation of this Baker-Hausdorff-Campbell on ( Q v , P v , Q u , P u ) yields only a few terms, whose coefficients are precisely those in Equation (8).

Proposition 14.
The action of ( U ρ , U σ , U τ ) on the wave function ψ E (v, u) may be represented by the following integral transforms: Proof. The proof is straightforward. It uses the action of each generator given by Equation (69) and then exponentiates it as action on ψ E (k). This leads to unitary factors for U ρ and U τ as well as a shift in the argument of ψ E (k) for U σ . This is an alternative form to the one obtained by Iwai-Rew [4].

Quantum Integrals of Motion and Consequences on the Schrödinger Wave Functions
In this subsection, we study the nature of the quantum integrals of motion. In particular, when a manifest separation of variables occurs in an operator (a1) F ρ exists since [ H, F ρ ] = 0 because ρ is a cyclic variable, as in the classical case.
(a2) F τ is due to the separation in the (v, u) Cartesian variables. Hence, commute.
(b2) However, there was an unproven claim of separation of variables in displaced parabolic coordinates [10] as a justification for the existence of ( F σ ). Hereafter, we provide an argument that may explain this claim. Consider the intermediate space (k, u), obtained after a partial v-Fourier transform from the (v, u) space. The Schrödinger equation in the (k, u) space for stationary states of the eigenvalue E admits solutions in the form of: This solution has a constant value C on the parabolas of equation u + 1 2 k 2 − E = C, where C is a shift of the parabolas. Consequently, the tangential derivative of any function of (k, u) along these parabolas is zero. This tangential derivative is just t · ∇, where ∇ = ∂ ∂k , ∂ ∂u is the gradient operator and t = (1, −k) is the tangent vector to the parabola. Thus, we have the following: Going back to the (v, u) space by the k-Fourier inverse transform, this tangential derivative reappears as a (v, u)−partial differential operator ∂ 2 ∂u ∂v − v , which is just − F σ . Therefore, the claim of reference [13] is only valid in the (k, u)-space.
(c) F 2 is due to a special form of separation of variables (Stäckel separation of variables), when one changes from Cartesian (v, u) to parabolic (x, y) coordinates by the following formulas: (u, v) = x 2 −y 2 2 , xy . After working out the expressions of the (x, y) partial derivatives, we end up with new expressions in (x, y) for ( H, F 2 ): These expressions can also be obtained from the quantization of classical expressions by replacing (p x , p y ) with −i ∂ ∂x , −i ∂ ∂y . The meaning of this Stäckel separation of variables is given in the following proposition.
Proof. Using the expression of H in parabolic coordinates given by Equation (82) and calling the corresponding wave function ψ E (x, y), we can transform the stationary Schrödinger equation Hψ E (x, y) = Eψ E (x, y) for all (x, y) = (0, 0) into the following condition: Since this condition is to be satisfied for all (x, y) = (0, 0), ψ E (x, y) should be a product of eigenfunctions of the two separate Schrödinger operators in x and in y, with opposite eigenvalues. For the confining quartic potential in x, it is known that this Schrödinger operator has discrete non-degenerate point spectrum with the eigenvalues E = n , with n ∈ N, which are bounded below. On the other hand, the Schrödinger operator with non-confining quartic potential in y has a continuous, real, non-degenerate, scattering-type spectrum. Thus, ψ E (v, u) ∼ ψ E (x, y) must be of the product form (up to a multiplicative constant): verifying the stationary Schrödinger equations: Then, after rewriting the expression of F 2 under the following form: we see that Hence, for H = E, we have necessarily −E as an eigenvalue of F 2 , where E = n is any eigenvalue of the confining quartic potential + 1 2 x 4 .

Integral Relation for Schrödinger Eigenfunctions of Quartic Oscillators
The passage to parabolic coordinates has a remarkable consequence on the eigenfunctions of quantum quartic oscillators, as stated in the following theorem. Theorem 3. The eigenfunctions of quantum confining and non-confining quartic oscillators verify the following integral relation: Proof. Recast the expression of ψ E (v, u) in parabolic coordinates to get the following: where µ E takes care of the fact that ψ E (v, u) ∼ ψ E (x, y) up to a multiplicative constant. Let F 2 operate on both sides of this equation. On the left-hand side, we get −Eψ E (y) according to Equation (86). On the right-hand side, we obtain: according to Equation (73). Hence, to have consistency between the two sides, one must require that ψ E (k) = ψ (+) E (k). Therefore, it follows that the eigenfunctions of the confining and non-confining quartic oscillators should verify the integral relation above.
Since the confining quartic oscillator is of dominant physical interest as a non-trivial field theory model in zero space dimension, an integral identity for its eigenfunctions can be derived from the previous identity via a "Wick-rotation" [17]. E (x) of the confining quartic oscillator fulfil the following integral identity: Proof. Since H replace y with y = ix in Equation (88). Then, as parity is a good quantum number for H (+) osc , the factor e −kxy in the integrand should be replaced by either sinh xx k or cosh xx k, according to the parity of the state of energy E = n for n ∈ N. This last integral identity was discovered long ago via the Weyl quantization of anharmonic oscillators [18] and recently rediscovered by [19].

An Application to Wave Propagation in Duct of Varying Section
In 2003, B J Forbes et al. [20] observed that the Webster equation describing the excess pressure p(u, t) in a fluid flowing in a duct with a slowly varying circular cross-sectional area S(u): 1 where u is the space coordinate along the axis of the duct, t the time, and c is the constant wave speed in the fluid, which can be turned into a Klein-Gordon equation with a potential V(u) for a wave function ψ(u, t): if ψ(u, t) = p(u, t) S(u) and Hence, if we require that the potential to be linear V(u) = (u − d), where d is an arbitrary distance, Equation (94) shows that the duct section should vary as S(u) =Ai(u − d), and consequently, Equation (93) is obtained from the two-dimensional free fall Schrödinger equation by "Wick-rotation" v = ict [17] (without a factor 1/2 in the partial derivatives). This Klein-Gordon Equation (93) admits a dynamical symmetry algebra which is the Wick rotated A FF . Note that the free fall F τ integral of motion now becomes the Klein-Gordon operator, and conversely, the free fall Hamiltonian operator is now an integral of the motion for the Klein-Gordon operator. The quadratic integral operator of the Klein-Gordon problem is now the Wick rotated free fall F 2 operator. In this case, one directly gets the integral identity for the Schrödinger eigenfunctions of the confining quartic anharmonic oscillator; see Equation (90).
The general solution of this Klein-Gordon wave equation φ d (t, u) is formally analogous to ψ E (v, u); it is given by the ω-integral for t > 0: with the initial conditions φ d (0, u) andφ d (0, u) = d dt φ d (t, u)| t=0 .

Theorem 4.
The angular frequency distribution f (ω) of φ d (t, u) is fully determined by the initial conditions.