1. Introduction
A nonrelativistic quantum particle with unit mass in
spacetime dimensions with coordinates
is given by the natural Lagrangian
. The wave function is expressed in terms of the propagator,
which, following Feynman’s intuitive proposal [
1], is obtained as,
where the (symbolic) integration is over all paths
that link the spacetime point
to
and where:
is the classical action calculated along
[
1,
2,
3].
The rigorous definition and calculation of (
2) are beyond our scope here. However, the
semiclassical approximation leads to the van Vleck–Pauli formula [
2,
3,
4,
5],
where
is the classical action calculated along the (supposedly unique (This condition is satisfied away from caustics [
2,
3,
6]. Moreover, (
5) and (
8) are valid only for
and for
, respectively, as discussed in
Section 4.)) classical path
from
and
. This expression involves data of the classical motion only. We note here also the van Vleck determinant
in the prefactor [
4,
5].
Equation (
4) is exact for a quadratic-in-the-position potentials in
dimension
that we consider henceforth.
For
, i.e., for a free nonrelativistic particle of unit mass in 1 + 1 dimensions with coordinates
X and
T, the result is [
1,
2,
3],
A harmonic oscillator with dissipation is in turn described by the Caldirola–Kanai (CK) Lagrangian and the equation of motion, respectively [
7,
8]. For constant damping and a harmonic frequency, we have,
with
and
. A lengthy calculation then yields the exact propagator [
2,
3,
9,
10,
11]:
where an irrelevant phase factor was dropped.
Inomata and his collaborators [
12,
13,
14,
15] generalized (
9) to a time-dependent frequency by redefining time,
, which allowed them to transform the time-dependent problem to one with a constant frequency (see
Section 2). Then, they followed by what they called a “
time-dependent conformal transformation”
such that:
which allowed them to derive the propagator from the free expression (
5). When spelled out, (
10) boils down to a generalized version, (
22), of the correspondence found by Niederer [
16].
It is legitimate to wonder:
in what sense are these transformations “conformal” ? In
Section 3, we explain that, in fact,
both mappings can be interpreted in the Eisenhart–Duval (E-D) framework as conformal transformations between two appropriate Bargmann spaces [
17,
18,
19,
20,
21]. Moreover, the change of variables
is a special case of the one put forward by Arnold [
22,
23] and is shown to be convenient to study time-dependent systems explicitly.
A bonus is the extension to the arbitrary time-dependent frequency
of the Maslov phase correction [
2,
4,
5,
6,
19,
24,
25,
26,
27,
28] even when no explicit solutions are available (see
Section 4).
In
Section 5.2, we illustrate our theory by the time-dependent Mathieu profile
, the direct analytic treatment of which is complicated.
2. The Junker–Inomata Derivation of the Propagator
Starting with a general quadratic Lagrangian in 1 + 1 spacetime dimensions with coordinates
and
t, Junker and Inomata derived the equation of motion [
12]:
which describes a nonrelativistic particle of unit mass with dissipation
. The driving force
can be eliminated by subtracting a particular solution
of (
11),
in terms of which (
11) becomes homogeneous,
This equation can be obtained from the time-dependent generalization of (
6),
The friction can be eliminated by setting
, which yields a harmonic oscillator with no friction, but with a shifted frequency [
29,
30,
31],
For and , for example, we obtain the usual harmonic oscillator with a constant shifted frequency,
The frequency is in general time-dependent, though
; therefore, (
14) is a
Sturm–Liouville equation that can be solved analytically only in exceptional cases.
Junker and Inomata [
12] followed another, more subtle path. Equation (
12) is a linear equation with time-dependent coefficients, the solution of which can be searched for within the ansatz (A similar transcription was used also by Rezende [
28].):
where
A,
B,
and are constants and
and
functions to be found. Inserting (
15) into (
12), putting the coefficients of the exponentials to zero, separating the real and imaginary parts, and absorbing a new integration constant into
provide us with the coupled system for
and
,
Manifestly,
. Inserting
into (
16) then yields the
Ermakov–Milne–Pinney (EMP) equation [
32,
33,
34] with time-dependent coefficients,
We note for later use that eliminating
would yield instead:
Conversely, the constancy of the r.h.s. here can be verified using Equation (
17). Equivalently, starting with the Junker–Inomata condition (
10),
To sum up, the strategy to follow is [
12,
35,
36]:
to solve first the EMP Equation (
18) for
;
Then, the trajectory is given by (
15).
Junker and Inomata showed, moreover, that substituting into (
13) the new coordinates:
allows us to present the Caldirola–Kanai action as (Surface terms do not change the classical equations of motion and multiply the propagator by an unobservable phase factor, and are therefore dropped.),
where we recognize the
action of a free particle of unit mass. One checks also directly that
satisfy the free equation, as they should. The conditions (
10) are readily verified.
The coordinates
X and
T describe a free particle; therefore, the propagator is (
5) (as anticipated by our notation). The clue of Junker and Inomata [
12] is that, conversely, trading
X and
T in (
5) for
x and
t allows deriving the propagator for the CK oscillator (see also [
11],
Section 5.1) (The extension of (
24) from
to all
t [
2,
3,
6,
11] is discussed in
Section 4.),
where we used the shorthands
, etc.
This remarkable formula says that in terms of “redefined time”,
, the problem is essentially one with a constant frequency. Equation (
24) is still implicit, though, as it requires solving first the coupled system (
17), which we can do only in particular cases.
When
where
, Equation (
12) describes a time-dependent oscillator with constant friction,
Then, setting
, Equation (
17) provide us with the EMP equation for
R, cf. (
18),
If, in addition, the
frequency is
constant then Equation (
26) is solved algebraically by:
Thus,
is a linear combination of
and
The spacetime coordinate transformation of
in (
22) simplifies to the friction-generalized form of that of Niederer [
16],
for which the general expression (
24) reduces to (
9) when
;
When the oscillator is turned off,
, but
, we have motion in a dissipative medium. The coordinate transformation propagator (
22) and (
24) become:
and:
respectively. A driving force
(e.g., terrestrial gravitation) could be added and then removed by
.
Further examples can be found in [
13,
14,
15]. An explicitly time-dependent example is presented in
Section 5.2.
3. The Eisenhart–Duval Lift
Further insight can be gained by “Eisenhart–Duval (E-D) lifting” the system to one higher dimension to what is called a “Bargmann space” [
17,
18,
19,
20,
21]. The latter is a
-dimensional manifold endowed with a Lorentz metric, the general form of which is:
which carries a covariantly constant null Killing vector
. Then:
Theorem 1 ([
18,
20]).
Factoring out the foliation generated by yields a nonrelativistic spacetime in dimensions. Moreover, the null geodesics of the Bargmann metric project to ordinary spacetime, consistent with Newton’s equations. Conversely, if is a solution of the nonrelativistic equations of motion, then its null lifts to Bargmann space are:where is an arbitrary initial value.
Let us consider, for example, a particle of unit mass with the Lagrangian of:
where
is a positive metric on a curved configuration space
Q with local coordinates
,
. The coefficients
and
may depend on time
t, and
is some (possibly time-dependent) scalar potential. The associated equations of motion are:
where the
are the Christoffel symbols of the metric
. For
,
and
for
, resp. for
, we obtain a (possible time-dependent) 1d oscillator without, resp. with, friction, Equation (
7) [
7,
8,
9,
29,
30,
31].
Equation (
34) can also be obtained by projecting a null-geodesic of
-dimensional Bargmann spacetime with coordinates
, whose metric is:
For
, we recover (
12).
Choosing
would describe motion with a time-dependent mass
. The friction can be removed by the conformal rescaling
, and the null geodesics of the rescaled metric describe, consistent with (
14), an oscillator with no friction, but with a time-dependent frequency,
[
37].
The friction term
in (
34) can be removed also by introducing a new time parameter
, defined by
[
21]. For
, for example, putting
eliminates the friction, but it does this at the price of obtaining a manifestly time-dependent frequency [
38,
39]:
3.1. The Junker–Inomata Ansatz as a Conformal Transformation
The approach outlined in
Section 2 admits a Bargmannian interpretation. For simplicity, we only consider the frictionless case
.
Theorem 2. The Junker–Inomata method of converting the time-dependent system into one with a constant frequency by switching from “real” to “fake time”, induces a conformal transformation between the Bargmann metrics: Proof. Putting
allows us to present the constant-frequency
(
19) as:
Then, with the notation
=
dζ/
dτ, we find,
Let us now recall that the null lift to the Bargmann space of a spacetime curve is obtained by subtracting the classical action as the vertical coordinate,
Setting here
and dropping surface terms yield, using the same procedure for the time-dependent-frequency case,
up to surface terms. Then, inserting all our formulae into (
38a) and (
38b) yields (
39), as stated. In Junker–Inomata language (
10),
. □
Our investigation has so far concerned classical aspects. Now, we consider what happens quantum mechanically. Restricting our attention at
space dimensions as before (In
, conformal invariance requires adding a scalar curvature term to the Laplacian.), we posit that the E-D lift
of a wave function
is equivariant,
Then, the massless Klein–Gordon equation for
associated with the
d Barmann metric implies the Schrödinger equation in 1+1 d,
where
is the Laplace–Beltrami operator associated with the metric. In
, it is of course
.
A conformal diffeomorphism
with conformal factor
,
projects to a spacetime transformation
. It is implemented on a wave function lifted to the Bargmann space as:
In
Section 4.2, these formulae are applied to the Niederer map (
73).
3.2. The Arnold Map
The general damped harmonic oscillator with time-dependent driving force
in 1 + 1 dimensions, (
11),
can be solved by an
Arnold transformation [
22,
23], which “straightens the trajectories” [
21,
29,
30,
31,
40]. To this end, one introduces new coordinates,
where
and
are solutions of the associated homogeneous Equation (
46) with
and
is a particular solution of the full Equation (
46). It is worth noting that (
47) allows checking, independently, the Junker–Inomata criterion in (
10). The initial conditions are chosen as,
Then, in the new coordinates, the motion becomes free [
22,
23],
Equation (
46) can be obtained by projecting a null geodesic of the Bargmann metric:
Completing (
47) by:
lifts the Arnold map to Bargmann spaces,
(In the Junker–Inomata setting (
10),
and
.),
The oscillator metric (
50) is thus carried conformally to the free one, generalizing earlier results [
18,
19,
41]. For the damped harmonic oscillator with
and
,
is a particular solution. When
, for example,
are two independent solutions of the homogeneous equation with initial conditions (
48) and provide us with:
In the undamped case,
; thus,
, and (
56) reduces to that of Niederer [
16] lifted to the Bargmann space [
19,
20],
The Junker–Inomata construction in
Section 2 can be viewed as a particular case of the Arnold transformation. We chose
and the two independent solutions:
The initial conditions (
48) at
imply
Then, spelling out (
51),
completes the lift of (
22) to Bargmann spaces. In conclusion, the one-dimensional damped harmonic oscillator is described by the conformally flat Bargmann metric,
The metric (
60) is manifestly conformally flat; therefore, its geodesics are those of the free metric,
. Then, using (
47) with (
58) yields:
The bracketed quantity here describes a constant-frequency oscillator with “time” . The original position, x, obtains a time-dependent “conformal” scale factor.
4. The Maslov Correction
As mentioned before, the semiclassical formula (
9) is correct only in the first oscillator half-period,
. Its extension for all
t involves the
Maslov correction. In the constant-frequency case with no friction, for example, assuming that
is not an integer, we have [
2,
3,
6],
where the integer:
is called the
Maslov index (
is the integer part of
x.).
ℓ counts the completed half-periods and is related also to the
Morse index, which counts the negative modes of
[
4,
5].
Now, we generalize (
62) to the
time-dependent frequency:
Theorem 3.In terms of and τ introduced in Section 2, Outside caustics, i.e., for , the propagator for the harmonic oscillator with the time-dependent frequency and friction is: -4.6cm0cm
Proof. In terms of the redefined coordinates:
cf. (
37), and using the notation
=
d/
dτ, the time-dependent oscillator Equation (
12) is taken into:
Thus, the problem is reduced to one with a
time-independent frequency,
in (
19) (We record for the sake of later investigations that (turning off
) (
Section 4) can be presented as:
where
is the
Schwarzian derivative of
[
42]). □
Let us now recall Formula (19) of Junker and Inomata in [
12], which tells us how propagators behave under the coordinate transformation
:
Here,
is the propagator of an oscillator with a time-dependent frequency and friction,
and
, respectively—the one we are trying to find.
is in turn the Maslov-extended propagator of an oscillator with no friction and a constant frequency, as in (
62). Then, the propagator for the
harmonic oscillator with a time-dependent frequency and friction, Equation (
64), is obtained using (
67).
Notice that (
64) is
regular at the points
where
. However, at caustics,
,
diverges, and we have instead (
66).
Henceforth, we limit our investigations to .
4.1. Properties of the Niederer Map
More insight is gained from the perspective of the generalized Niederer map (
22). We first study their properties in some detail. For simplicity, we chose, in the rest of this section,
and
and
.
We start with the observation that the Niederer map (
22) becomes singular where the cosine vanishes, i.e., where:
because
is an increasing function by (
21). Moreover, each interval:
is mapped by (
22) onto the full range
. Therefore, the inverse mapping is
multivalued, labeled by integers
k,
where
with
the principal determination, i.e., in
.
Then,
and
imply that:
Therefore, the intervals and are joined at and the form a partition of the time axis,
Returning to (
64) (which is (
62) with
), we then observe that, whereas the propagator is regular at
, it diverges at caustics,
cf. (
65). Thus,
, and:
Thus,
maps the full
T-line into
with
an internal point. Conversely,
is an
internal point of
. The intervals
cover again the time axis,
By (
61) the classical trajectories are regular at
. Moreover, for arbitrary initial velocities,
implying that after a half-period
, all classical motions are focused at the same point. The two entangled sets of intervals are shown in
Figure 1.
The Niederer map (
57) “E-D lifts” to the Bargmann space.
Theorem 4. The E-D lift of the inverse of the Niederer map (57), which we shall denote by (), is: Proof. These formulae follow at once by inverting (
57), at once with the cast
. Alternatively, it could also be proven as for Theorem 2.
For each integer
k (
78) maps the real line
into the “open strip” [
19]
with
defined in (
71). Their union covers the entire Bargmann manifold of the oscillator. □
Now, we pull back the free dynamics by the multivalued inverse (
78). We put
for simplicity. The free motion with initial condition
,
E-D lifts by (
78) to:
consistent with
, as can be checked directly. Note that the
s coordinate oscillates with a doubled frequency.
At (where the Niederer maps are joined), we have Thus, the pull backs of the Bargmann lifts of free motions are glued to smooth curves;
Similarly, at t caustics
, we infer from (
80) that for all initial velocities
a and for all
ℓ Thus, the lifts are again smooth at
, and after each half-period, all motions are focused above the initial position
.
4.2. The Propagator by the Niederer Map
Now, we turn to quantum dynamics. Our starting point is the free propagator (
5), which (as mentioned before) is valid only for
. Its extension to all
T involves the
sign of
[
19].
Let us explain this subtle point in some detail. First of all, we notice that the usual expression (
5) involves a square root, which is double-valued, obliging us to
choose one of its branches. Which one we choose is irrelevant: it is a mere gauge choice. However, once we do choose one, we must stick to our choice. Take, for example, the one for which
, then the prefactor in (
5) is:
Let us now consider what happens when
changes sign. Then, the prefactor becomes multiplied by
so it becomes,
for the same choice of the square root,
In conclusion, the formula valid for all
T is,
where:
is the free action calculated along the classical trajectory. Let us underline that (
82) already involves a “Maslov jump”
, which, for a free particle, happens at
. For
, we have
.
Accordingly, the wave function
of a free particle is, by (
1),
Now, we pull back the free dynamics using the multivalued inverse Niederer map. It is sufficient to consider the constant-frequency case
and to denote time by
t. Let
t belong to the range of
in (
73),
Then, applying the general formulae in
Section 3.1 yields [
19],
However, the second exponential in the middle line combines with the integrand in the braces in the last line to yield
the action calculated along the classical oscillator trajectory,
Thus, using the equivariance, we end up with,
Now, we recover the Maslov jump, which comes from the first line here. For simplicity, we consider again and denote .
Firstly, we observe that the conformal factor
has a constant sign in the domain
and changes sign at the end points. In fact,
The cosine enters into the van Vleck factor, while the phase combines with
. Recall now that
divides
into two pieces,
cf.
Figure 1. However,
is precisely where the tangent changes sign: this term contributes to the phase in
and
in
. Combining the two shifts, we end up with the phase:
which is the Maslov jump at
.
Intuitively, the multivalued “exports” to the oscillator at the phase jump of the free propagator at . Crossing from to shifts the index ℓ by one.
6. Conclusions
The Junker–Inomata–Arnold approach yields (in principle) the exact propagator for any quadratic system by switching from a
time-dependent to a
constant frequency and redefined time,
The propagator (
64)–(
66) is then derived from the result known for the constant frequency. A straightforward consequence is the Maslov jump for arbitrary time-dependent frequency
: everything depends only on the product
.
By switching from
t to
, the Sturm–Liouville-type difficulty is not eliminated, but only transferred to that of finding
following the procedure outlined in
Section 2. We have to first solve EMP Equation (
18) for
(which is nonlinear and has time-dependent coefficients) and then integrate
; see (
21). Although this is as difficult to solve as solving the Sturm–Liouville equation, it provides us with theoretical insights.
When no analytic solution is available, we can resort to numerical calculations.
The Junker–Inomata approach of
Section 2 is interpreted as a Bargmann-conformal transformation between time-dependent and constant frequency metrics; see Equation (
39).
Alternatively, the damped oscillator can be converted to a free system by the generalized Niederer map (
22), whose Eisenhart–Duval lift (
47)–(
51) carries the conformally flat oscillator metric (
60) to the flat Minkowski space.
Two sets of points play a distinguished role in our investigations: the
in (
71) and the
in (
75). The
divides the time axis into domains
of the (generalized) Niederer map (
22). Both classical motions and quantum propagators are
regular at
, where these intervals are joined. The
are in turn the caustic points where all
classical trajectories are focused, and the
quantum propagator becomes
singular.
While the “Maslov phase jump” at caustics is well established when the frequency is constant,
, its extension to the time-dependent case
is more subtle. In fact, the proofs we are aware of [
25,
26,
27,
28] use sophisticated mathematics, or a lengthy direct calculation of the propagator [
44]. A bonus from the Junker–Inomata transcription (
10) we followed here is to provide us with a straightforward extension valid to an arbitrary
. Caustics arise when (
65) holds, and then, the phase jump is given by (
88).
The subtle point mentioned above comes from the standard (but somewhat sloppy) expression (
5), which requires choosing a branch of the double-valued square root function. Once this is done, the sign change of
induces a phase jump
. Our “innocent-looking” factor
is in fact the Maslov jump for a free particle at
(obscured when one considers the propagator for
only). Moreover, it then becomes the key tool for the oscillator: intuitively, the multivalued inverse Niederer map repeats, again and again, the same jump. The details are discussed in
Section 4.
The transformation (
10) is related to the
nonrelativistic “Schrödinger” conformal symmetries of a free nonrelativistic particle [
45,
46,
47], later extended to the oscillator [
16] and an inverse-square potential [
48]. These results can in fact be derived using a time-dependent conformal transformation of the type (
10) [
19,
42].
The above results are readily generalized to higher dimensions. For example, the oscillator frequency can be time-dependent, uniform electric and magnetic fields, and a curl-free “Aharonov–Bohm” potential (a vortex line [
49]) can also be added [
41]. Further generalization involves a Dirac monopole [
50].
Alternative ways to relate free and harmonically trapped motions are studied, e.g., in [
51,
52,
53,
54]. Motions with the Mathieu profile were considered also in [
55].