Next Article in Journal
Application of Manifold Corrections in Tidal Evolution of Exoplanetary Systems
Next Article in Special Issue
Relativistic Fermion and Boson Fields: Bose-Einstein Condensate as a Time Crystal
Previous Article in Journal
Improvement of Unconstrained Optimization Methods Based on Symmetry Involved in Neutrosophy
Previous Article in Special Issue
Quantifying Complementarity via Robustness of Asymmetry
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Time-Dependent 4D Quantum Harmonic Oscillator and Reacting Hydrogen Atom

by
Ashot S. Gevorkyan
1,2,* and
Aleksander V. Bogdanov
3,4
1
Institute for Informatics and Automation Problems, NAS of RA, 1, P. Sevak St., Yerevan 0014, Armenia
2
Institute of Chemical Physics, NAS of RA, 5/2, P. Sevak St., Yerevan 0014, Armenia
3
Faculty of Applied Mathematics and Control Processes, St. Petersburg State University, 7/9 Universitetskaya Nab., 199034 St. Petersburg, Russia
4
Center for Advanced Digital Technologies, St. Petersburg State Marine Technical University, Lotsmanskaya D. 3, 190121 St. Petersburg, Russia
*
Author to whom correspondence should be addressed.
Symmetry 2023, 15(1), 252; https://doi.org/10.3390/sym15010252
Submission received: 12 December 2022 / Revised: 9 January 2023 / Accepted: 12 January 2023 / Published: 16 January 2023
(This article belongs to the Special Issue Symmetry and Asymmetry in Quantum Mechanics)

Abstract

:
With the help of low-dimensional reference equations (ordinary differential equations) and the corresponding coordinate transformations, the non-stationary 4D quantum oscillator in an external field is reduced to an autonomous form. The latter, in particular, reflects the existence of a new type of dynamical symmetry that reduces the equation of motion of a non-stationary oscillator to an autonomous form that does not change with time. By imposing an additional constraint on the wave function of the isotropic oscillator, we have obtained the total wave functions of the reacting hydrogen atom in two different cases: (a) when the non-stationary frequency has two asymptotic values and there is no external field; and (b) when, in addition to the non-stationary frequency, an external force acts on the hydrogen atom. The transition S-matrix elements of various elementary atomic–molecular processes are constructed. The probabilities of quantum transitions of the hydrogen atom to others, including new bound states, are studied in detail, taking into account the influence of external forces.

1. Introduction

The hydrogen atom still remains one of the most important systems for research due to its wide application in practice, and it is also a reference problem for testing new theoretical concepts of non-relativistic quantum mechanics. In addition, it is well known that the hydrogen atom is perhaps the most important and rich in hidden symmetries system in quantum physics. Already from Pauli’s investigation, the Lie algebra of the symmetry group can be identified as S O ( 4 ) . However, a reasonable question arises: where does this “accidental” symmetry come from? Fock addressed this question thoroughly in his celebrated work “On the theory of the hydrogen atom” [1] and explained the degeneracy of the energy levels of the Kepler problem (or hydrogen atom) in terms of the dynamical symmetry group S O ( 4 ) . In particular, he showed that if one makes a stereographic projection in momentum space and rescales momenta with energy eigenvalues, then the problem is equivalent to a geodesic flow on a sphere with unit radius S 3 . A dynamical symmetry group isomorphic to O ( 4 ) was used by Schwinger [2] to derive the Green’s function of the Coulomb field in the momentum representation, but a more complete discussion of the group theory application to the hydrogen atom was given by Bander and Itzykson [3].
As a generalization of Levi–Civita’s parabolic coordinates, Kustaanheimo and Stiefel [4] introduced a coordinate transformation as applied to the corresponding classical problem. Boiteux used these coordinate transformations and proved that the bound states of the hydrogen atom correspond to the energy states of a 4D isotropic oscillator with one additional constraint [5].
The main goal of this work is to construct a mathematically rigorous and consistent representation of the reacting hydrogen atom (RHA) using the explicit and hidden symmetries of a 4D isotropic harmonic oscillator. Recall that RHA is understood to mean a hydrogen atom that is free in the ( i n ) asymptotic state, while in the ( o u t ) asymptotic state, as a result of an elementary atomic–molecular process, it transits to another quantum state, including the possibility of forming a bound state with another atom or molecule (chemical reaction).
To achieve this goal from a mathematical point of view, it is necessary to solve the problem of separating variables in a second-order partial differential equation (PDE) of a fairly general form, which in itself is a very complex and important independent mathematical problem. Note that as Miller showed in his monograph [6], the problem of separation of variables in PDE is closely related to one of the most powerful methods of modern mathematics and mathematical physics, namely the theory of Lie algebras.
In this paper, at the first stage, it is shown that with the help of the so-called reference equations, i.e., a system of four ordinary differential equations of small dimension, the original non-stationary Schrödinger equation for an isotropic 4D oscillator is reduced to an autonomous form with a constant frequency. At the second stage, using coordinate transformations and imposing an additional constraint on the wave function, the initial equation for the oscillator is reduced to the equation of a hydrogen atom in an external field. In the third stage, we construct the S-matrix elements of the transitions and the probabilities of the corresponding elementary processes involving the hydrogen atom.
In conclusion, the question of the further development of the theory of the reacting hydrogen atom is considered in detail, taking into account its presence in a medium with a finite temperature. The latter can be extremely useful for its wide use in solving many applied problems in physics and chemistry.

2. Statement of the Problem

Let us consider the problem of a 4D oscillator with a non-stationary frequency under the action of an external force:
i Ψ t = H ^ ( x , t ) Ψ , x = ( x 1 , , x 4 ) R 4 ,
where R 4 denotes the 4D configuration space and the operator H ^ ( x , t ) is the Hamiltonian of the 4D quantum harmonic oscillator (QHO), which is represented in the following form:
H ^ ( x , t ) = l = 1 4 { 1 2 2 x l 2 + 1 2 Ω 2 ( t ) x l 2 F ( t ) x l } , x l , t ( , + ) .
Mind that Equations (1) and (2) are written in units of m = = 1 , where m is the mass of the oscillator and is the Plank constant.
We suppose that the problem is described by two ( i n ) , (at t ) and ( o u t ) , (at t + ) asymptotic subspaces. In particular, let the frequency Ω ( t ) and the external force F ( t ) satisfy the following conditions:
lim t ± Ω ( t ) = Ω ± , lim t ± F ( t ) = 0 ,
where Ω and Ω + denote constant frequencies in the ( i n ) and ( o u t ) asymptotic subspaces, respectively.
Based on the conditions (3), it follows that the quantum system in the asymptotic subspaces ( i n ) and ( o u t ) is described by pure states of orthonormal 4D bases of an isotropic oscillator.
Below, we will consider two typical cases:
  • When there is no external force acting on the non-stationary oscillator, i.e., F ( t ) 0 and we have the case of a quantum parametric oscillator (QPO), and
  • When an external force acts on the non-stationary oscillator, i.e., F ( t ) 0 .

3. Reduction of the 4D QHO Equation to the Autonomous Form

Let us consider the classical oscillator problem, which will play a key role in the further constructions. Without taking into account external influence, it satisfies the following homogeneous second-order ordinary differential equation (ODE):
ξ ¨ + Ω 2 ( t ) ξ = 0 , ξ ˙ = d ξ / d t .
The solution of Equation (4) will be sought in the form:
ξ ( t ) = σ ( t ) e i γ ( t ) , σ ( t ) = | ξ ( t ) | , γ ( t ) = t d t σ 2 ( t ) ,
where σ ( t ) specifies the length scale at time t, and τ = γ ( t ) / Ω denotes the corresponding time scale.
It is assumed that the solution ξ ( t ) satisfies the initial condition:
ξ ( t ) e i Ω t , a t t .
As a consequence, it is obvious that lim t σ ( t ) = 1 .
Now, let us consider the equation of an oscillator that is affected by an external force:
η ¨ + Ω 2 ( t ) η F ( t ) = 0 , η ( ) = η ˙ ( ) = 0 .
Using Green’s function method (see for example [7]), we can construct the solution of the inhomogeneous Equation (7) by representing it in terms of the complex solution of the homogeneous Equation (4):
η ( t ) = 1 Ω ξ * ( t ) d ( t ) + ξ ( t ) d * ( t ) , d ( t ) = i Ω t ξ ( t ) F ( t ) d t .

3.1. The Parametric Quantum Oscillator

Given the representation (5) and the condition (6) (see also [7]), the solution of Equations (1) and (2) can be written as:
Ψ ( x , t ) = 1 σ 2 l = 1 4 exp i 2 σ ˙ σ x l 2 Ψ ( y , τ ) ,
where σ ˙ = d σ / d t denotes the derivative with respect to time, and Ψ ( y , τ ) is the wave function of a 4D isotropic harmonic oscillator (IHO) with constant frequency Ω . In addition, y = ( y 1 , , y 4 ) is the displacement of the quantum oscillator in the 4D space, and y l = x l / σ ( t ) denotes a new space coordinate.
Substituting the expression (9) into Equations (1) and (2) and taking into account (4), we obtain the following Schrödinger equation:
i τ Ψ = 1 2 l = 1 4 2 y l 2 + Ω 2 y l 2 Ψ ,
where { y , τ } denotes a new space–time.
It is easy to see that (10) is an autonomous equation that does not change its form during the evolution. In other words, with the help of an ordinary differential Equation (4), which we will call a low-dimensional reference equation, the master Equations (1) and (2) are reduced to an autonomous form. The latter obviously reflects the presence of a hidden dynamical symmetry in the quantum system under consideration. The sequence of mathematical operations that makes it possible to reduce the master equation to an autonomous form will be referred to below as the reference equation method (REM).
Let us consider the new coordinate system ( y 1 , y 2 , y 3 , y 4 ) ( r , θ , ϕ , ψ ) first proposed by Kustaanheimo and Stiefel (KS transformation) [4] (see also [8]):
y 1 = r cos θ 2 cos ϕ + ψ 2 , y 2 = r sin θ 2 cos ϕ ψ 2 , y 3 = r cos θ 2 sin ϕ + ψ 2 , y 4 = r sin θ 2 sin ϕ ψ 2 ,
where r = l = 1 4 y l 2 , in addition, ϕ [ 0 , 2 π ] and θ , ψ [ 0 , + π ] .
The length element in these coordinates would be:
d s 2 = d r 2 + r 2 4 d θ 2 + d ϕ 2 + d ψ 2 + 2 cos θ d ϕ d ψ .
Rewriting Equation (10) in the new coordinates and replacing r 2 = ϱ , we can find the following equation for the 4D IHO:
i τ Ψ ˘ = 2 { 1 ϱ ϱ ϱ 2 ϱ + 1 ϱ [ 1 sin θ θ sin θ θ + 1 sin 2 θ 2 ϕ 2 + 2 ψ 2 2 cos θ 2 ϕ ψ ] 1 4 Ω 2 ϱ } Ψ ˘ .
Recall that Ψ ˘ ( { ϱ } , τ ) = Ψ ( y , τ ) is the wave function of the IHO in the new coordinates; in addition, the following notation is used for the set of new coordinates; { ϱ } = ( ϱ , θ , ϕ , ψ ) R 4 .
Since Equation (14) has an autonomous form, its solution can be represented as:
Ψ ˘ ( { ϱ } , τ ) = e i E τ Λ ( { ϱ } ) ,
where E denotes the total energy of the 4D IHO in the new space–time. Substituting (13) into Equation (14), we obtain:
1 ϱ ϱ ϱ 2 ϱ + 1 ϱ G ^ ( θ , ϕ , ψ ) + 1 2 E 1 4 Ω 2 ϱ Λ ( { ϱ } ) = 0 .
where the angular operator G ˜ ( θ , ϕ , ψ ) has the form:
G ˜ ( θ , ϕ , ψ ) = 1 sin θ θ sin θ θ + 1 sin 2 θ 2 ϕ 2 + 2 ψ 2 2 cos θ 2 ϕ ψ .
Representing the wave function Λ ( { ϱ } ) in the form:
Λ ( { ϱ } ) = Λ ( 0 ) ( ϱ ) Φ ( θ , ϕ , ψ ) ,
we can obtain two new equations from Equation (14). In particular, the angular part of the wave function satisfies the following equation:
1 sin θ θ sin θ θ + 1 sin 2 θ 2 ϕ 2 + 2 ψ 2 2 cos θ 2 ϕ ψ + λ Φ = 0 ,
where λ is some constant, which will be defined below. The radial part of the wave function satisfies the following second-order ODE:
1 ϱ 2 d d ϱ ϱ 2 d d ϱ + Ω 2 4 + E 2 ϱ λ ϱ 2 Λ ( 0 ) ( ϱ ) = 0 .
Concerning Equation (16), as we know, this is the equation of a 3D spherical rotator, which has the following solution:
Φ ( θ , ϕ , ψ ) = Θ J K M ( θ ) e i ( K ϕ + M ψ ) ,
where
Θ J K M ( θ ) = z | K M | / 2 ( 1 z ) | K + M | / 2 F ( α , β , γ ; z ) , z = 1 2 ( 1 cos θ ) ,
In addition, the following notations are made:
λ = J ( J + 1 ) , J = 0 , 1 , , M , K = J , J 1 , , 1 J , J .
Mind that F ( α , β , γ ; z ) is a hyper-geometric function, and its parameters are defined as follows:
α = J + | K M | / 2 , β = J + | K + M | / 2 + | K M | / 2 , γ = 1 + | K M | .
As for the equation for the radial wave function (17), it is not difficult to see that it describes the wave function of the hydrogen atom and has an exact solution (see for example [9]):
Λ n J ( 0 ) ( ϱ ) = e ρ / 2 ρ J L n J 1 2 J + 1 ( ρ ) , ρ = 2 Ω ϱ , n = 1 , 2 ,
where L n m ( x ) denotes the associated Laguerre polynomial; in addition, recall that for an isolated hydrogen atom, the energy is E = Ω = 1 / n 2 .
Thus, by combining (9), (13), (15) and (18), we can write the wave function of the reacting 4D isotropic quantum oscillator in the explicit form:
Ψ ˘ M ( { ϱ } , t ) = 1 σ 2 l = 1 4 exp i 2 σ ˙ σ y l 2 ( { ϱ } ) Λ n J ( 0 ) ( ϱ ) Θ J K M ( θ ) e i ( K ϕ + M ψ E τ ) .
where M = ( n , J , K , M ) denotes a set of quantum numbers describing the quantum state. However, given that l = 1 4 y l 2 = ϱ , the wave function (20) can be rewritten as:
Ψ ˘ M ( { ϱ } , t ) = 1 σ 2 e i E τ ( t ) i σ ˙ σ ϱ / 2 Λ n J ( 0 ) ( ϱ ) Θ J K M ( θ ) e i ( K ϕ + M ψ ) .
Returning to the problem of the hydrogen atom, we note that the wave function of the reacting hydrogen atom is easily found by requiring the following constraint condition Ψ ˘ / ψ = 0 for the total wave function of the 4D IHO. It is easy to check that this condition is equivalent to the case when we put the quantum number M = 0 in the wave function for the isotropic harmonic oscillator. In other words, we can write the exact expression for the wave function IHO evolving from the asymptotic subspace ( i n ) with frequency Ω to the asymptotic subspace ( o u t ) , where the oscillator frequency is equal to Ω + :
Ψ ˘ M 0 ( { ϱ 0 } , t ) = C n J K σ 2 ( t ) e i E τ ( t ) i σ ˙ σ ϱ / 2 Λ n J ( 0 ) ( ϱ ) Θ J K ( θ ) e i K ϕ ,
where C n J K is some constant to be defined below; in addition, M 0 = ( n , J , K ) denotes the set of quantum numbers characterizing the hydrogen atom, Θ J K ( θ ) = Θ J K M ( θ ) | M = 0 and { ϱ 0 } = ( ϱ , θ , ϕ ) . In the limit t or in the ( i n ) subspace, the total wave function (22) must match with the asymptotic wave function Φ ˘ M 0 ( { ϱ 0 } ) :
Ψ ˘ M 0 ( { ϱ 0 } , t ) t Φ ˘ M 0 ( { ϱ 0 } ) = C n J K Λ n J ( 0 ) ( ϱ ) Θ J K ( θ ) e i K ϕ .
Note that the constant C n J K can be found from the normalization condition for the wave function of an isolated hydrogen atom. In particular, if this is completed in the 3D Euclidean space R 3 ( ϱ , θ , ϕ ) , then we obtain:
C n J K = 0 0 π 0 2 π Λ n J ( 0 ) ( ϱ ) Θ J K ( θ ) 2 ϱ 2 sin θ d ϱ d θ d ϕ = ( n J 1 ) ! 2 n ( n + J ) ! · 2 n 3 2 .
However, we complete all the construction of the problem in the 4D space, so the normalization of the 3D wave function (22) must also be performed in the 4D space. This new constant C n J K will obviously be different from C n J K represented by formula (23). Since the coefficient C n J K does not depend on the quantum number K, in what follows, we will replace C n J K with C n J .

3.2. Non-Stationary QHO under the Influence of an External Force

Let now consider the case of the non-stationary QHO subjected to an external force, i.e., F ( t ) 0 . Our goal now is to reduce the problem to the QPO case.
The wave function of a non-stationary QHO, taking into account the external force action, can be represented as [10]:
Ψ ¯ ( x , t ) = l = 1 4 exp i [ η ˙ ( t ) x ¯ l + S l ( t ) ] Ψ ¯ ( x ¯ , τ ) , x ¯ l = x l η ( t ) ,
where η ( t ) and S l ( t ) are unknown functions that we define below based on the requirement of reducing Equations (1) and (2) to an autonomous form; in addition, x ¯ = x ¯ ( x ¯ 1 , , x ¯ 4 ) R 4 denotes the displacement of the oscillator center of mass in 4D Euclidean space. Substituting (24) into Equations (1) and (2), we obtain:
i Ψ ¯ t = l = 1 4 { 1 2 2 x l ¯ 2 + 1 2 Ω 2 ( t ) x ¯ l 2 + η ¨ + Ω 2 ( t ) η F ( t ) x ¯ l + S ˙ l ( t ) 1 2 η ˙ 2 + 1 2 Ω 2 ( t ) η 2 F ( t ) η } Ψ ¯ .
It is easy to see that Equation (25) can be reduced to the case of a parametric oscillator with F ( t ) = 0 . In particular, this can be completed if we require that the following inhomogeneous equation holds for the classical oscillator (see also (7)):
η ¨ + Ω ( t ) η = F ( t ) ,
and also for the following equalities:
L l ( t ) = S ˙ l ( t ) = 1 2 η ˙ 2 1 2 Ω 2 ( t ) η 2 + F ( t ) η ,
where L l ( t ) is the classical Lagrangian of the 1D oscillator.
Thus, we were able to determine the unknown functions η l ( t ) and S l ( t ) and, what is very important, reduce the problem to the case of the QPO. Combining (9), (24) and (27) for the wave function of a non-stationary QHO that is subjected to an external force, we can write the following expression:
Ψ ¯ ˘ M ( { ϱ ¯ } , t ) = 1 σ 2 ( t ) e i E τ i σ ˙ σ ϱ ¯ / 2 Λ n J ( 0 ) ( ϱ ¯ ) l = 1 4 e i η ˙ σ y ¯ l ( { ϱ ¯ } ) + S l ( t ) Θ J K M ( θ ) e i ( K ϕ + M ψ ) ,
where y ¯ l = x ¯ l / σ denotes the scaled coordinate shift; in addition, ϱ ¯ = l = 1 4 y ¯ l 2 is the radius; { ϱ ¯ } = ( ϱ ¯ , θ , ϕ , ψ ) R 4 denotes a set of coordinates, and finally, S l ( t ) = t L l ( t ) d t is the action of the 1D classical oscillator. If we now set M = 0 , then the wave function (28) will describe the reacting hydrogen atom, taking into account the external fields’ influence.
In the end, we note that the wave function (28) must be normalized in a 4D space before using it in some other constructions, such as calculating the matrix elements of quantum transitions, for example.

4. Transition Probabilities between (in) and (out) Asymptotic Quantum States

In this section, we construct the transition probabilities between two asymptotic states ( i n ) and ( o u t ) for the QPO. Particular attention will be paid to the case when the quantum number M = 0 , i.e., the problem of the hydrogen atom under the influence of external fields, which leads to a change in the frequency.
Taking into account the condition (3) for the external force, we can claim that in ( i n ) and ( o u t ) asymptotic subspaces, the QPO will be described by purely quantum states:
Φ ˘ M ( { ϱ } ) = Λ n J ( 0 ) ( ϱ ) Θ J K M ( θ ) e i ( K ϕ + M ψ ) ,
where Λ n J ( 0 ) ( ϱ ) = e ρ / 2 ( ρ ) J L n J 1 2 J + 1 ( ρ ) and ρ = 2 Ω ϱ denote the radial parts of the wave functions of the hydrogen atom in the corresponding asymptotic subspaces. In other words, our task will be to calculate the probability of transitions between asymptotic states Φ ˘ M ( { ϱ } ) and Φ ˘ M + ( { ϱ } ) . In particular, if we assume that in the ( i n ) asymptotic subspace the QPO was in a pure quantum state Φ ˘ M ( { ϱ } ) , then the evolving total wave function of the oscillator Ψ ˘ M ( { ϱ } , t ) can be represented in the form:
Ψ ˘ M ( { ϱ } , t ) = M S M M ( t ) Φ ˘ M + ( { ϱ } ) .
Calculating the expansion coefficient of the sum S M M ( t ) , we obtain:
S M ; M ( t ) = R 4 Ψ ˘ M ( { ϱ } , t ) Φ ˘ M + ( { ϱ } ) ¯ d 4 V ( { y } ) ,
where the infinitesimal volume in the new coordinates is defined as follows:
d 4 V ( { y } ) = ( x 1 , x 2 , x 3 , x 4 ) ( ϱ , ϑ , ϕ , ψ ) d ϱ d ϑ d ϕ d ψ = σ 4 ( t ) ( y 1 , y 2 , y 3 , y 4 ) ( ϱ , ϑ , ϕ , ψ ) d ϱ d ϑ d ϕ d ψ ,
in addition, in expression (31), the symbol { · } ¯ denotes the complex conjugate function. As for the determinant of coordinate transformations in the integral (31), we can calculate it explicitly using expressions (11):
( · ) ( · ) = 1 2 ϱ cos θ 2 cos ϕ + ψ 2 ϱ 2 sin θ 2 cos ϕ + ψ 2 ϱ 2 cos θ 2 sin ϕ + ψ 2 ϱ 2 cos θ 2 sin ϕ + ψ 2 1 2 ϱ cos θ 2 sin ϕ ψ 2 ϱ 2 sin θ 2 cos ϕ ψ 2 ϱ 2 cos θ 2 cos ϕ ψ 2 ϱ 2 cos θ 2 cos ϕ ψ 2 1 2 ϱ cos θ 2 sin ϕ + ψ 2 ϱ 2 sin θ 2 sin ϕ + ψ 2 ϱ 2 cos θ 2 cos ϕ + ψ 2 ϱ 2 cos θ 2 cos ϕ + ψ 2 1 2 ϱ sin θ 2 sin ϕ ψ 2 ϱ 2 cos θ 2 sin ϕ ψ 2 ϱ 2 sin θ 2 cos ϕ ψ 2 ϱ 2 sin θ 2 cos ϕ ψ 2 .
Recall that the calculation shows that the value of the determinant is ( · ) ( · ) = ϱ sin θ 16 .
The transition S-matrix element can be obtained from the expression (31) in the limit t + , the square of the modulus of which will denote the transition probability between Φ ˘ M ( { ϱ } ) and Φ ˘ M + ( { ϱ } ) quantum states:
W M ; M = lim t + | S M ; M ( t ) | 2 .
For definiteness, below, we are going to study quantum transitions involving the hydrogen atom in more detail. We are going to consider two asymptotic states described by the following wave functions:
Φ ˘ M 0 ( { ϱ 0 } ) = e ρ / 2 ρ J L n J 1 2 J + 1 ( ρ ) , Φ ˘ M 0 + ( { ϱ 0 } ) = e ν ϱ / 2 ( ν ρ ) J L n J 1 2 J + 1 ( ν ρ ) ,
where ρ = ρ = 2 Ω ϱ and ν = Ω + / Ω .
Note that if we assume that the wave function Φ ˘ M 0 ( { ϱ 0 } ) describes an isolated hydrogen atom, then the wave function Φ ˘ M 0 + ( { ϱ 0 } ) that includes the frequency Ω + Ω will describe the conditionally bonded hydrogen atom. Our task will be to calculate the probability of transition between these two quantum states W M 0 M 0 . Taking into account that the elementary atomic–molecular process occurs in the 4D Euclidean space R 4 ( ϱ , θ , ϕ , ψ ) , the transition S-matrix element can be written as:
S M 0 ; M 0 ( t ) = σ 4 ( t ) 16 0 0 π 0 π 0 2 π Ψ ˘ M 0 ( { ϱ 0 } , t ) Φ ˘ M 0 + ( { ϱ 0 } ) ¯ ϱ sin θ d ϱ d θ d ϕ d ψ ,
where M 0 = ( n , J , K ) and M 0 = ( n , J , K ) denote the sets of quantum numbers of the hydrogen atom, respectively, in ( i n ) and ( o u t ) asymptotic states.
Before proceeding to the calculation of the S-matrix element by the formula (34), it is necessary to normalize both the total and the asymptotic wave functions. In particular, for the full wave function, this normalization constant will have the form:
C n J K ( t ) = 4 σ 2 ( t ) 0 0 π 0 π 0 2 π Ψ ˘ M 0 ( { ϱ 0 } , t ) 2 ϱ sin θ d ϱ d θ d ϕ d ψ 1 / 2 .
Substituting the expression for the hydrogen wave function (22) into (35), we find:
C n J ( t ) = 4 σ 2 ( t ) 0 ρ 2 J + 1 e ρ L n J 1 2 J + 1 ( ρ ) 2 d ϱ 1 / 2 .
Now, using a simple formula [11]:
0 z λ e z L m λ ( z ) L n λ ( z ) d z = Γ ( n + λ + 1 ) n ! δ m n , Re ( λ ) > 1 ,
we can calculate the integral (36) for the normalization constant:
C n J ( t ) = 4 σ 2 ( t ) ( n J 1 ) ! Γ ( n + J + 1 ) .
Given (37), the wave function of the reacting hydrogen atom Ψ ̲ ˘ M 0 ( { ϱ 0 } , t ) normalized in 4D space can be written as:
Ψ ˘ ̲ M 0 ( { ϱ 0 } , t ) = C n J ( t ) Ψ ˘ M 0 ( { ϱ 0 } , t ) .
Note that the normalized wave function for an isolated hydrogen atom Φ ̲ ˘ M 0 ( { ϱ 0 } ) can be obtained from (38) in the t limit, when σ ( t ) 1 and C n J ( t ) C n J 0 = 4 ( n J 1 ) ! Γ ( n + J + 1 ) . We normalize the asymptotic wave function Φ ̲ ˘ M 0 + ( { ϱ 0 } ) in a similar way.
Now, given that the following more general integral can be computed exactly [11]:
0 z α 1 e p z L m λ ( a z ) L n β ( b z ) d z = Γ ( α ) ( λ + 1 ) m ( β + 1 ) n m ! n ! p α j = 0 m ( m ) j ( α ) j ( λ + 1 ) j j ! a p j × k = 0 n ( n ) k ( j + α ) k ( β + 1 ) k k ! b p k , Re ( α ) > 0 , Re ( p ) > 0 , m , n N ,
we can proceed to the calculation of the S-matrix element (34).
Taking into account that the deformation of the wave function of the reacting hydrogen occurs only along the radial coordinate and substituting (22), (37), (38) and (33) into the expression (34), we can find:
S M 0 ; M 0 ( t ) = δ J J δ K K 0 e i σ ˙ σ ϱ / 2 e ρ / 2 ρ J L n J 1 2 J + 1 ( ρ ) e ν ϱ / 2 ( ν ρ ) J L n J 1 2 J + 1 ( ν ρ ) ρ d ρ = B n J K 0 B n J K 0 δ J J δ K K ν J 0 e q ( t ) ρ ρ J + J + 1 L n J 1 2 J + 1 ( ρ ) L n J 1 2 J + 1 ( ν ρ ) d ρ ,
where q ( t ) = i σ ˙ ( t ) σ ( t ) + 1 + ν / 2 and B n J 0 = C n J 0 / 4 . The integral (40) can be simplified by writing it in the form:
S n J ; n J ( t ) = A n J ; n J 0 ρ 2 J + 1 e q ( t ) ρ L n J 1 2 J + 1 ( ρ ) L n J 1 2 J + 1 ( ν ρ ) d ρ ,
where A n J ; n J = ν J ( n J 1 ) ! ( n J 1 ) ! Γ ( n + J + 1 ) Γ ( n + J + 1 ) .
Continuing the calculations taking into account (39), we find:
S n J ; n J ( t ) = A n J ; n J 1 q 2 ( J + 1 ) ( t ) Γ ( 2 J + 2 ) ( 2 J + 1 ) n J 1 ( 2 J + 1 ) n J 1 ( n J 1 ) ! ( n J 1 ) ! × j = 0 n J 1 ( J n + 1 ) j ( 2 J + 2 ) j ( 2 J + 1 ) j j ! q j ( t ) k = 0 n J 1 ( J n + 1 ) k ( j + 2 J + 2 ) k ( 2 J + 2 ) k k ! ν q ( t ) k
With the help of (39), we can explicitly write the first two elements of the S-matrix:
S 10 ; 10 ( t ) = 4 [ 1 + ν + i σ ˙ ( t ) σ ( t ) ] 2 , S 21 ; 21 ( t ) = 16 ν [ 1 + ν + i σ ˙ ( t ) σ ( t ) ] 4 .
As we can easily see, there are strict selection rules for quantum transitions. Only those transitions are possible that can be represented in the form S n ( n 1 ) ; n ( n 1 ) ; n N . Returning to the transition probability issue, we note that in this case, it is very important to study the S-matrix elements’ behaviour in the t + limit, when the frequency Ω ( t ) tends to a constant value Ω + . To carry out the necessary research, we need to turn to the classical oscillator problem (4), since the modulus of its solution determines the function σ ( t ) = | ξ ( t ) | .
In the ( o u t ) subspace, i.e., in the t + limit, Equation (4) is simplified and, accordingly, its solution can be formally represented as a combination of the incident and reflected waves [7]:
ξ ( t ) = c 1 e i Ω + t c 2 e i Ω + t , | c 1 | 2 | c 2 | 2 = Ω Ω + ,
where c 1 and c 2 are some complex constants that can be found after solving Equation (4) over the entire time interval < t < + .
Given (43) for the σ ( t ) function, it is easy to find the following expression:
σ ( t ) = c ¯ 1 2 + c ¯ 2 2 + 2 c ¯ 1 c ¯ 2 cos ( δ 1 δ 2 2 Ω + t ) ,
where c ¯ j = | c j | and δ j = arg ( c j ) , j = 1 , 2 .
Continuing the calculations and taking into account (44), we find that:
σ ˙ ( t ) σ ( t ) = 2 c ¯ 1 c ¯ 2 Ω + sin ( δ 1 δ 2 2 Ω + t ) ,
and accordingly, the S-matrix elements will be:
S 10 ; 10 ( t ) = 4 χ 2 ( t ) , S 21 ; 21 ( t ) = 16 ν χ 4 ( t ) ,
where χ ( t ) = [ 1 + ν i 2 c ¯ 1 c ¯ 2 Ω + sin ( δ 1 δ 2 2 Ω + t ) ] .
From (46), it becomes obvious that in the limit of t + , the elements of the S-matrix do not converge to constant limits but perform periodic oscillations. Based on this, it is necessary to additionally integrate the elements of the S-matrix over the oscillation period:
S ¯ n J ; n J = Ω + π 0 π / Ω + S n J ; n J ( t ) d t ,
Accordingly, the transition probability per unit time will be determined as follows:
W n J ; n J = | S ¯ n J ; n J | 2 .
The analysis shows that quantum transitions can be calculated in a similar way also for the model of a reacting hydrogen atom when it is affected by an external force. However, in this case, there is no special selection law, and quantum transitions between any states are possible.

5. Conclusions

In addition to the obvious relevance of the results obtained for the theory of atomic–molecular collisions and reactions, where the importance of exact results cannot be overestimated, the developed approach is also of great importance for the theory of fundamental processes. The fact is that from the moment of the first attempts to construct a field theory to obtain fundamental solutions, the main approach was to reduce the dimension of the problem. However, after Fock’s fundamental work [1], it became obvious to everyone that the problem can also be simplified by increasing its dimension and bringing it to a form in which the symmetry of the problem is realized in a natural way and most completely. Recall that Fock used the dynamic symmetry group S O ( 4 ) to explain the degenerate energy levels of the hydrogen atom. In addition, which is very important, he proved the equivalence of the considered problem to the propagation of a geodesic flow on the surface of 3D sphere of unit radius.
The point becomes much more complicated when we consider an arbitrarily moving non-stationary 4D quantum oscillator, which, logically, should break the hidden symmetries of a quantum system. Moreover, one would expect that the connection between the problem of an isotropic harmonic oscillator in a real 4D space and the Kepler problem in an ordinary 3D space would also be broken [5,12,13]. However, as we showed in this paper using the example of a non-stationary 4D-isotropic oscillator, by the low-dimensional reference equation method, the initial Schrödinger equation reduces to an autonomous form. However, as we have shown in this article, using the low-dimensional reference Equation (4), the original Schrödinger Equations (1) and (2) for the non-stationary 4D isotropic oscillator is reduced to an autonomous form (10). In other words, we have proved the existence of a new type of dynamical symmetry that preserves the invariant form of the master equation with respect to time in which the reference equation plays a key role. Note that the Schrödinger equation also reduces to an autonomous form for the case when an external force acts on the oscillator. After reducing the equation to an autonomous form, the coordinate transformation KS (11) is used, which makes it possible to reduce the problem to the Kepler problem and solve it exactly both for the wave function of the oscillator (see expressions (21) and (28)) and for the hydrogen atom (see (22) and (28), substituting in it M = 0 ).
An important result is the calculation of the explicit expressions of transition S-matrix elements for a reacting hydrogen atom in the absence of an external field (41) and (42). This problem is basic for calculating rearrangement processes in the theory of atomic–molecular collisions, so that in the near future, we can expect a revision of many results for three- and many-particle reactive collisions and an increase in the efficiency of their calculations. On the other hand, in view of the special role of the hydrogen atom in the theory of fundamental processes, it seems to us that the performed calculations make it possible to shed new light on some approaches to the formulation of the theory of elementary particles, especially in terms of the influence of hidden symmetries on the selection rules.
Finally, note that the method can be easily generalized to the case when the oscillator is immersed in a thermostat with a finite temperature (see [14] and also [15]). This, in particular, will make it possible to construct the thermodynamics of the reacting hydrogen atom in various media from the first principles of quantum mechanics, which will be very instructive both from the point of view of the foundations of quantum physics and for solving many important applied problems of physics and chemistry.

Author Contributions

Conceptualization, methodology and investigation A.S.G.; Investigation and funding acquisition A.V.B. All authors have read and agreed to the published version of the manuscript.

Funding

Gevorkyan A.S. is grateful to grant N 21T-1B059 of the Science Committee of Armenia, which partially funded this work. The research was partially funded by the Ministry of Science and Higher Education of the Russian Federation as part of the World-Class Research Center program: Advanced Digital Technologies (contract No. 075-15-2022-312 dated 20 April 2022).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The author declares no conflict of interest.

Abbreviations

RHOreacting harmonic oscillator
PDEpartial differential equation
QHOquantum harmonic oscillator
IHOisotropic harmonic oscillator
REMreference equation method
ODEordinary differential equation
QPOquantum parametric oscillator

References

  1. Fock, V. Theory of the hydrogen atom. Z. Phys. 1935, 98, 145–154. (In German) [Google Scholar] [CrossRef]
  2. Schwinger, J. Coulomb Green’s Function. Math. Phys. 1964, 5, 1606–1608. [Google Scholar] [CrossRef]
  3. Bander, M.; Itzykson, C. Group Theory and the Hydrogen Atom (I). Rev. Mod. Phys. 1966, 38, 330. [Google Scholar] [CrossRef] [Green Version]
  4. Kustaanheimo, P.; Stiefel, E. Perturbation Theory of Kepler Motion. Based on Spinor Regularization. J. Reine Angew. Math. 1965, 218, 204–219. [Google Scholar] [CrossRef]
  5. Boiteux, M. The three-dimensional hydrogen atom as a restricted four-dimensional harmonic oscillator. Physica 1973, 65, 381–395. [Google Scholar] [CrossRef]
  6. Miller, W., Jr. Symmetry and Separation of Variables; Encyclopedia of Mathematics and Its Applications: Volume 4; Cambridge University Press: Cambridge, UK, 1977; 281p. [Google Scholar]
  7. Baź, A.N.; Zeĺdovich, Y.B.; Perelomov, A.M. Scattering Reactions and Decays in Nonrelativistic Quantum Mechanics; Nauka: Moscow, Russia, 1971. (In Russian) [Google Scholar]
  8. Hurt, N.E. Geometric Quantization in Action; Mathematics and its Applications 8; D. Reidel Publishing Company: Dordrecht, The Netherlands, 1983; 332p. [Google Scholar]
  9. Landau, L.D.; Lifshitz, E.M. Quantum Mechanics: Non-Relativistic Theory; Course of Theoretical Physics; Pergamon Press Ltd.: Oxford, UK, 1965; 632p, Volume 3. [Google Scholar]
  10. Husimi, K. Miscellanea in Elementary Quantum Mechanics, II. Prog. Theor. Phys. 1953, 9, 381. [Google Scholar] [CrossRef]
  11. Generalized Laguerre Polynomials. (1998–2023 Wolfram Research, Inc.). Available online: https://functions.wolfram.com/ (accessed on 20 April 2022).
  12. Bergmann, D.; Frishman, Y. A Relation between the Hydrogen Atom and Multidimensional Harmonic Oscillators. J. Math. Phys. 1965, 6, 1855. [Google Scholar] [CrossRef]
  13. Komarov, L.I.; Van Hoang, L. Generalized Kustaanheimo-Stiefel transformations. Theor. Math. Phys. 1994, 99, 437–440. [Google Scholar] [CrossRef]
  14. Gevorkyan, A.S.; Bogdanov, A.V.; Mareev, V.V. Hidden Dynamical Symmetry and Quantum Thermodynamics from the First Principles: Quantized Small Environment. Symmetry 2021, 13, 1546. [Google Scholar] [CrossRef]
  15. Gevorkyan, A.S.; Bogdanov, A.V.; Mareev, V.V.; Movsesyan, K.A. Theoretical and Numerical Study of Self-Organizing Processes in a Closed System Classical Oscillator and Random Environment. Mathematics 2022, 10, 3868. [Google Scholar] [CrossRef]
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Gevorkyan, A.S.; Bogdanov, A.V. Time-Dependent 4D Quantum Harmonic Oscillator and Reacting Hydrogen Atom. Symmetry 2023, 15, 252. https://doi.org/10.3390/sym15010252

AMA Style

Gevorkyan AS, Bogdanov AV. Time-Dependent 4D Quantum Harmonic Oscillator and Reacting Hydrogen Atom. Symmetry. 2023; 15(1):252. https://doi.org/10.3390/sym15010252

Chicago/Turabian Style

Gevorkyan, Ashot S., and Aleksander V. Bogdanov. 2023. "Time-Dependent 4D Quantum Harmonic Oscillator and Reacting Hydrogen Atom" Symmetry 15, no. 1: 252. https://doi.org/10.3390/sym15010252

APA Style

Gevorkyan, A. S., & Bogdanov, A. V. (2023). Time-Dependent 4D Quantum Harmonic Oscillator and Reacting Hydrogen Atom. Symmetry, 15(1), 252. https://doi.org/10.3390/sym15010252

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop