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
. 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
. 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
. A dynamical symmetry group isomorphic to
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 4
D 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 asymptotic state, while in the 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 4
D oscillator with a non-stationary frequency under the action of an external force:
where
denotes the 4
D configuration space and the operator
is the Hamiltonian of the 4
D quantum harmonic oscillator (QHO), which is represented in the following form:
Mind that Equations (
1) and (
2) are written in units of
, where
m is the mass of the oscillator and
ℏ is the Plank constant.
We suppose that the problem is described by two
(at
) and
(at
) asymptotic subspaces. In particular, let the frequency
and the external force
satisfy the following conditions:
where
and
denote constant frequencies in the
and
asymptotic subspaces, respectively.
Based on the conditions (
3), it follows that the quantum system in the asymptotic subspaces
and
is described by pure states of orthonormal 4
D 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., and we have the case of a quantum parametric oscillator (QPO), and
When an external force acts on the non-stationary oscillator, i.e., .
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):
The solution of Equation (
4) will be sought in the form:
where
specifies the length scale at time
t, and
denotes the corresponding time scale.
It is assumed that the solution
satisfies the initial condition:
As a consequence, it is obvious that
Now, let us consider the equation of an oscillator that is affected by an external force:
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):
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:
where
denotes the derivative with respect to time, and
is the wave function of a 4
D isotropic harmonic oscillator (IHO) with constant frequency
. In addition,
is the displacement of the quantum oscillator in the 4
D space, and
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:
where
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
first proposed by Kustaanheimo and Stiefel (KS transformation) [
4] (see also [
8]):
where
, in addition,
and
.
The length element in these coordinates would be:
Rewriting Equation (
10) in the new coordinates and replacing
, we can find the following equation for the 4
D IHO:
Recall that is the wave function of the IHO in the new coordinates; in addition, the following notation is used for the set of new coordinates; .
Since Equation (
14) has an autonomous form, its solution can be represented as:
where
E denotes the total energy of the 4
D IHO in
the new space–time. Substituting (
13) into Equation (
14), we obtain:
where the angular operator
has the form:
Representing the wave function
in the form:
we can obtain two new equations from Equation (
14). In particular, the angular part of the wave function satisfies the following equation:
where
is some constant, which will be defined below. The radial part of the wave function satisfies the following second-order ODE:
Concerning Equation (
16), as we know, this is the equation of a 3
D spherical rotator, which has the following solution:
where
In addition, the following notations are made:
Mind that
is a hyper-geometric function, and its parameters are defined as follows:
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]):
where
denotes the associated Laguerre polynomial; in addition, recall that for an isolated hydrogen atom, the energy is
.
Thus, by combining (
9), (
13), (
15) and (
18), we can write the wave function of the reacting 4
D isotropic quantum oscillator in the explicit form:
where
denotes a set of quantum numbers describing the quantum state. However, given that
, the wave function (
20) can be rewritten as:
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
for the total wave function of the 4
D IHO. It is easy to check that this condition is equivalent to the case when we put the quantum number
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
with frequency
to the asymptotic subspace
, where the oscillator frequency is equal to
:
where
is some constant to be defined below; in addition,
denotes the set of quantum numbers characterizing the hydrogen atom,
and
. In the limit
or in the
subspace, the total wave function (
22) must match with the asymptotic wave function
:
Note that the constant
can be found from the normalization condition for the wave function of an isolated hydrogen atom. In particular, if this is completed in the 3
D Euclidean space
, then we obtain:
However, we complete all the construction of the problem in the 4
D space, so the normalization of the 3
D wave function (
22) must also be performed in the 4
D space. This new constant
will obviously be different from
represented by formula (
23). Since the coefficient
does not depend on the quantum number
K, in what follows, we will replace
with
.
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., . 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]:
where
and
are unknown functions that we define below based on the requirement of reducing Equations (
1) and (
2) to an autonomous form; in addition,
denotes the displacement of the oscillator center of mass in 4
D Euclidean space. Substituting (
24) into Equations (
1) and (
2), we obtain:
It is easy to see that Equation (
25) can be reduced to the case of a parametric oscillator with
. In particular, this can be completed if we require that the following inhomogeneous equation holds for the classical oscillator (see also (
7)):
and also for the following equalities:
where
is the classical Lagrangian of the 1
D oscillator.
Thus, we were able to determine the unknown functions
and
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:
where
denotes the scaled coordinate shift; in addition,
is the radius;
denotes a set of coordinates, and finally,
is the action of the 1
D classical oscillator. If we now set
, 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 4
D 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 and for the QPO. Particular attention will be paid to the case when the quantum number , 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
and
asymptotic subspaces, the QPO will be described by purely quantum states:
where
and
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
and
. In particular, if we assume that in the
asymptotic subspace the QPO was in a pure quantum state
, then the evolving total wave function of the oscillator
can be represented in the form:
Calculating the expansion coefficient of the sum
, we obtain:
where the infinitesimal volume in the new coordinates is defined as follows:
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):
Recall that the calculation shows that the value of the determinant is
The transition
S-matrix element can be obtained from the expression (
31) in the limit
, the square of the modulus of which will denote the transition probability between
and
quantum states:
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:
where
and
.
Note that if we assume that the wave function
describes an isolated hydrogen atom, then the wave function
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
. Taking into account that the elementary atomic–molecular process occurs in the 4
D Euclidean space
, the transition
S-matrix element can be written as:
where
and
denote the sets of quantum numbers of the hydrogen atom, respectively, in
and
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:
Substituting the expression for the hydrogen wave function (
22) into (
35), we find:
Now, using a simple formula [
11]:
we can calculate the integral (
36) for the normalization constant:
Given (
37), the wave function of the reacting hydrogen atom
normalized in 4
D space can be written as:
Note that the normalized wave function for an isolated hydrogen atom
can be obtained from (
38) in the
limit, when
and
. We normalize the asymptotic wave function
in a similar way.
Now, given that the following more general integral can be computed exactly [
11]:
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:
where
and
The integral (
40) can be simplified by writing it in the form:
where
Continuing the calculations taking into account (
39), we find:
With the help of (
39), we can explicitly write the first two elements of the
S-matrix:
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
. 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
limit, when the frequency
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
.
In the
subspace, i.e., in the
limit, Equation (
4) is simplified and, accordingly, its solution can be formally represented as a combination of the incident and reflected waves [
7]:
where
and
are some complex constants that can be found after solving Equation (
4) over the entire time interval
Given (
43) for the
function, it is easy to find the following expression:
where
and
Continuing the calculations and taking into account (
44), we find that:
and accordingly, the
S-matrix elements will be:
where
.
From (
46), it becomes obvious that in the limit of
, 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:
Accordingly, the transition probability per unit time will be determined as follows:
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
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 3
D sphere of unit radius.
The point becomes much more complicated when we consider an arbitrarily moving non-stationary 4
D 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 4
D space and the Kepler problem in an ordinary 3
D space would also be broken [
5,
12,
13]. However, as we showed in this paper using the example of a non-stationary 4
D-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 4
D 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
).
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.