Next Article in Journal
On the Controllability of a System Modeling Cell Dynamics Related to Leukemia
Next Article in Special Issue
Special Issue: “Symmetries in Quantum Mechanics and Statistical Physics”
Previous Article in Journal
Characteristic Length Scale during the Time Evolution of a Turbulent Bose-Einstein Condensate
Previous Article in Special Issue
On the Supersymmetry of the Klein–Gordon Oscillator
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Time-Dependent Conformal Transformations and the Propagator for Quadratic Systems †

1
School of Physics and Astronomy, Sun Yat-sen University, Zhuhai 519082, China
2
Institut Denis Poisson, Tours University-Orléans University, UMR 7013, F-37200 Tours, France
*
Author to whom correspondence should be addressed.
This paper celebrates the 90th birthday of Akira Inomata.
Symmetry 2021, 13(10), 1866; https://doi.org/10.3390/sym13101866
Submission received: 7 July 2021 / Revised: 13 August 2021 / Accepted: 17 August 2021 / Published: 3 October 2021
(This article belongs to the Special Issue Symmetries in Quantum Mechanics and Statistical Physics)

Abstract

:
The method proposed by Inomata and his collaborators allows us to transform a damped Caldirola–Kanai oscillator with a time-dependent frequency to one with a constant frequency and no friction by redefining the time variable, obtained by solving an Ermakov–Milne–Pinney equation. Their mapping “Eisenhart–Duval” lifts as a conformal transformation between two appropriate Bargmann spaces. The quantum propagator is calculated also by bringing the quadratic system to free form by another time-dependent Bargmann-conformal transformation, which generalizes the one introduced before by Niederer and is related to the mapping proposed by Arnold. Our approach allows us to extend the Maslov phase correction to an arbitrary time-dependent frequency. The method is illustrated by the Mathieu profile.
PACS:
03.65.-w quantum mechanics; 03.65.Sq semiclassical theories and applications; 04.20.-q classical general relativity

1. Introduction

A nonrelativistic quantum particle with unit mass in d + 1 spacetime dimensions with coordinates x , t is given by the natural Lagrangian L = 1 2 x ˙ 2 V ( x , t ) . The wave function is expressed in terms of the propagator,
ψ ( x , t ) = K ( x , t | x , t ) ψ ( x , t ) d x
which, following Feynman’s intuitive proposal [1], is obtained as,
K ( x , t | x , t ) = exp i A ( γ ) D ,
where the (symbolic) integration is over all paths γ ( t ) =   x ( t ) , t that link the spacetime point ( x , t ) to ( x , t ) and where:
A ( γ ) = t t L γ ( t ) , γ ˙ ( t ) , t d t
is the classical action calculated along γ ( t ) [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],
K ( x , t | x , t ) = i 2 π 2 A ¯ x x 1 / 2 exp i A ¯ ( x , t | x , t ) ,
where A ¯ ( x , t | x , t ) = t t L ( γ ¯ ( t ) , γ ¯ ˙ ( t ) , t ) d t 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 0 < T T and for 0 < t t < π , respectively, as discussed in Section 4.)) classical path γ ¯ ( τ ) from ( x , t ) and ( x , t ) . This expression involves data of the classical motion only. We note here also the van Vleck determinant 2 A ¯ x x in the prefactor [4,5].
Equation (4) is exact for a quadratic-in-the-position potentials in 1 + 1 dimension V ( x , t ) = 1 2 ω 2 ( t ) x 2 that we consider henceforth.
For ω 0 , 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],
K f r e e ( X , T | X , T ) = 1 2 π i ( T T ) 1 / 2 exp i ( X X ) 2 2 ( T T ) .
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,
L C K = 1 2 e λ 0 t d x d t 2 ω 0 2 x 2 ,
d 2 x d t 2 + λ 0 d x d t + ω 0 2 x = 0
with λ 0 = const . > 0 and ω 0 = const . . A lengthy calculation then yields the exact propagator [2,3,9,10,11]:
K C K ( x , t | x , t ) = Ω 0 e λ 0 2 ( t + t ) 2 π i sin Ω 0 ( t t ) 1 2 × exp i Ω 0 2 sin Ω 0 ( t t ) ( x 2 e λ 0 t + x 2 e λ 0 t ) cos Ω 0 ( t t ) 2 x x e λ 0 t + t 2 ,
Ω 0 2 = ω 0 2 1 4 λ 0 2 ,
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, t τ , 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 ( x , t ) ( X , T ) such that:
x = f ( T ) X ( T ) exp 1 2 λ 0 T , t = g ( T ) , where f 2 ( T ) = d g d T ,
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 x , t X , T 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 ω ( t ) 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 ω 2 ( t ) = a 2 q cos 2 t , a , b const . , 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 x ˜ and t, Junker and Inomata derived the equation of motion [12]:
x ˜ ¨ + λ ˙ ( t ) x ˜ ˙ + ω 2 ( t ) x ˜ = F ( t ) ,
which describes a nonrelativistic particle of unit mass with dissipation λ ( t ) . The driving force F ( t ) can be eliminated by subtracting a particular solution h ( t ) of (11), x ( t ) = x ˜ ( t ) h ( t ) , in terms of which (11) becomes homogeneous,
x ¨ + λ ˙ ( t ) x ˙ + ω 2 ( t ) x = 0 .
This equation can be obtained from the time-dependent generalization of (6),
L C K = 1 2 e λ ( t ) [ x ˙ 2 ω 2 ( t ) x 2 ] .
The friction can be eliminated by setting x ( t ) = y ( t ) e λ ( t ) / 2 , which yields a harmonic oscillator with no friction, but with a shifted frequency [29,30,31],
y ¨ + Ω 2 ( t ) y = 0 where Ω 2 ( t ) = ω 2 ( t ) λ ˙ 2 ( t ) 4 λ ¨ ( t ) 2 .
For λ ( t ) = λ 0 t and ω = ω 0 = const . , for example, we obtain the usual harmonic oscillator with a constant shifted frequency, Ω 2 = ω 0 2 λ 0 2 / 4 = const .
The frequency is in general time-dependent, though Ω = Ω ( t ) ; 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].):
x ( t ) = ρ ( t ) A e i ω ¯ τ ( t ) + B e i ω ¯ τ ( t ) ,
where A, B, and ω ¯ are constants and ρ ( t ) and τ ( t ) 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 A , B provide us with the coupled system for ρ ( t ) and τ ( t ) ,
ρ ¨ + λ ˙ ρ ˙ + ( ω 2 ( t ) ω ¯ 2 τ ˙ 2 ) ρ = 0 ,
τ ˙ ( t ) ρ 2 ( t ) e λ ( t ) = 1 .
Manifestly, τ ˙ > 0 . Inserting τ ˙ into (16) then yields the Ermakov–Milne–Pinney (EMP) equation [32,33,34] with time-dependent coefficients,
ρ ¨ + λ ˙ ρ ˙ + ω 2 ( t ) ρ = e 2 λ ( t ) ω ¯ 2 ρ 3 .
We note for later use that eliminating ρ would yield instead:
ω ¯ 2 = 1 τ ˙ 2 ω 2 ( t ) 1 2 τ τ ˙ + 3 4 τ ¨ τ ˙ 2 λ ¨ 2 λ ˙ 2 4 .
Conversely, the constancy of the r.h.s. here can be verified using Equation (17). Equivalently, starting with the Junker–Inomata condition (10),
ω 2 ( t ) = f ¨ f 2 f ˙ 2 f 2 + λ ˙ 2 4 + λ ¨ 2 .
To sum up, the strategy to follow is [12,35,36]:
  • to solve first the EMP Equation (18) for ρ ;
  • to integrate (17),
    τ ( t ) = t e λ ( u ) ρ 2 ( u ) d u .
Then, the trajectory is given by (15).
Junker and Inomata showed, moreover, that substituting into (13) the new coordinates:
T = tan [ ω ¯ τ ( t ) ] ω ¯ , X = x e λ ( t ) 2 τ ˙ ( t ) 1 2 sec [ ω ¯ τ ( t ) ] ,
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.),
A C K = t t L C K d t = T T 1 2 d X d T 2 d T ,
where we recognize the action of a free particle of unit mass. One checks also directly that X , T 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 0 < ω ¯ ( τ τ ) < π to all t [2,3,6,11] is discussed in Section 4.),
K o s c ( x , t | x , t ) = ω ¯ e λ + λ 2 ( τ ˙ τ ˙ ) 1 2 2 π i sin [ ω ¯ ( τ τ ) ] 1 2 × exp i ω ¯ 2 sin [ ω ¯ ( τ τ ) ] ( x 2 e λ τ ˙ + x 2 e λ τ ˙ ) cos [ ω ¯ ( τ τ ) ] 2 x x e λ + λ 2 ( τ ˙ τ ˙ ) 1 2 ,
where we used the shorthands λ = λ ( t ) , τ = τ ( t ) , 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 λ ( t ) = λ 0 t where λ 0 = const . 0 , Equation (12) describes a time-dependent oscillator with constant friction,
    x ¨ + λ 0 x ˙ + ω 2 ( t ) x = 0 .
    Then, setting R ( t ) = ρ ( t ) e λ 0 t / 2 , Equation (17) provide us with the EMP equation for R, cf. (18),
    R ¨ + Ω 2 ( t ) R ω ¯ 2 R 3 = 0 , where Ω 2 ( t ) = ω 2 ( t ) λ 0 2 4 ;
  • If, in addition, the frequency is constant  ω ( t ) = ω 0 = const . , then Equation (26) is solved algebraically by:
    ω ¯ 2 = ω 0 2 λ 0 2 / 4 , R = 1 ρ ( t ) = e λ 0 t / 2 , τ ( t ) = t .
    Thus, x ( t ) is a linear combination of e 1 2 λ 0 t sin ω ¯ t and e 1 2 λ 0 t cos ω ¯ t . The spacetime coordinate transformation of ( x , t ) ( X , T ) in (22) simplifies to the friction-generalized form of that of Niederer [16],
    T = tan ( ω ¯ t ) ω ¯ , X = x exp 1 2 λ 0 t sec ( ω ¯ t ) ,
    for which the general expression (24) reduces to (9) when λ 0 = 0 ;
  • When the oscillator is turned off, ω 0 = 0 , but λ 0 > 0 , we have motion in a dissipative medium. The coordinate transformation propagator (22) and (24) become:
    X = 2 x 1 + exp λ 0 t , T = 2 λ 0 1 exp ( λ 0 t ) 1 + exp ( λ 0 t )
    and:
    K d i s s ( x , t | x , t ) = λ 0 2 π i [ exp ( λ 0 t ) exp ( λ 0 t ) ] 1 2 × exp i λ 0 2 ( x x ) 2 exp ( λ 0 t ) exp ( λ 0 t ) ,
    respectively. A driving force F 0 (e.g., terrestrial gravitation) could be added and then removed by x x + ( F 0 / λ 0 ) t .
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 d + 1 + 1 -dimensional manifold endowed with a Lorentz metric, the general form of which is:
g μ ν d x μ d x ν = g i j ( x , t ) d x i d x j + 2 d t d s 2 V ( x , t ) d t 2 ,
which carries a covariantly constant null Killing vector s . Then:
Theorem 1
([18,20]). Factoring out the foliation generated by s yields a nonrelativistic spacetime in d + 1 dimensions. Moreover, the null geodesics of the Bargmann metric g μ ν project to ordinary spacetime, consistent with Newton’s equations. Conversely, if ( γ ( t ) , t ) is a solution of the nonrelativistic equations of motion, then its null lifts to Bargmann space are:
γ ( t ) , t , s ( t ) , s ( t ) = s 0 A ( γ ) = s 0 t L ( γ ( r ) , r ) d r
where s 0 is an arbitrary initial value.
Let us consider, for example, a particle of unit mass with the Lagrangian of:
L = 1 2 α ( t ) g i j ( x k ) x ˙ i x ˙ j β ( t ) V ( x i , t ) ,
where g i j ( x k ) d x i d x j is a positive metric on a curved configuration space Q with local coordinates x i , i = 1 , , d . The coefficients α ( t ) and β ( t ) may depend on time t, and V ( x i , t ) is some (possibly time-dependent) scalar potential. The associated equations of motion are:
d 2 x i d t 2 + Γ j k i d x j d t d x k d t α ˙ α d x i d t = α β g i j j V ,
where the Γ j k i are the Christoffel symbols of the metric g i j . For d = 1 , g i j = δ i j and V = 1 2 ω 2 ( t ) x 2 for α = β = 1 , resp. for α = β 1 = e λ ( t ) , 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 d + 1 + 1 -dimensional Bargmann spacetime with coordinates ( x μ ) = ( x i , t , s ) , whose metric is:
g μ ν d x μ d x ν = 1 α g i j d x i d x j + 2 d t d s 2 β V d t 2 .
For α = β 1 = e λ ( t ) , we recover (12).
Choosing λ ( t ) = ln m ( t ) would describe motion with a time-dependent mass m ( t ) . The friction can be removed by the conformal rescaling x y = m x , and the null geodesics of the rescaled metric describe, consistent with (14), an oscillator with no friction, but with a time-dependent frequency, Ω 2 = ω 2 m ¨ 2 m + ( m ˙ 2 m ) 2 [37].
The friction term ( α ˙ / α ) x ˙ i in (34) can be removed also by introducing a new time parameter t ˜ , defined by d t ˜ = α d t [21]. For λ ( t ) = λ 0 t , for example, putting t ˜ = e λ 0 t / λ 0 eliminates the friction, but it does this at the price of obtaining a manifestly time-dependent frequency [38,39]:
d 2 x d t ˜ 2 + Ω ˜ 2 ( t ˜ ) x = 0 , Ω ˜ 2 ( t ˜ ) = ω 2 t ˜ 2 λ 0 2 .

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 λ = 0 .
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”,
t τ ( t ) , ξ = τ ˙ x
induces a conformal transformation between the Bargmann metrics:
d x 2 + 2 d t d s ω 2 ( t ) x 2 d t 2 frequency ω 2 ( t )
d ξ 2 + 2 d τ d σ ω ¯ 2 ξ 2 d τ 2 , frequency ω ¯ = const .
d ξ 2 + 2 d τ d σ ω ¯ 2 ξ 2 d τ 2 = τ ˙ ( t ) d x 2 + 2 d t d s ω 2 ( t ) x 2 d t 2 .
Proof. 
Putting μ = ln τ ˙ allows us to present the constant-frequency ω ¯ (19) as:
ω ¯ 2 = τ ˙ 2 ω 2 ( t ) 1 2 μ ¨ + 1 4 μ ˙ 2 .
Then, with the notation ζ = dζ/dτ, we find,
ζ 2 = τ ˙ 1 x ˙ 2 + 1 4 μ ˙ 2 x 2 1 2 μ ¨ x 2 + d d t 1 2 μ ˙ x 2 .
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,
d σ = L ( ξ , ξ , τ ) d τ = 1 2 ( ξ 2 ω ¯ 2 ξ 2 ) d τ
Setting here ξ = x ( τ ˙ 1 / 2 and dropping surface terms yield, using the same procedure for the time-dependent-frequency case,
d σ = d s = 1 2 x ˙ 2 ω 2 ( t ) x 2 d t
up to surface terms. Then, inserting all our formulae into (38a) and (38b) yields (39), as stated. In Junker–Inomata language (10), f ( t ) = τ ˙ 1 / 2 sec ( ω ¯ τ ) , g ( t ) = ( ω ¯ ) 1 tan ( ω ¯ τ ) . □
Our investigation has so far concerned classical aspects. Now, we consider what happens quantum mechanically. Restricting our attention at d = 1 space dimensions as before (In d > 2 , conformal invariance requires adding a scalar curvature term to the Laplacian.), we posit that the E-D lift ψ ˜ of a wave function ψ is equivariant,
ψ ˜ ( x , t , s ) = e i s ψ ( x , t ) s ψ ˜ = i ψ ˜ .
Then, the massless Klein–Gordon equation for ψ ˜ associated with the 1 + 1 + 1 = 3 d Barmann metric implies the Schrödinger equation in 1+1 d,
Δ g ψ ˜ = 0 i t ψ = 2 2 Δ x + V ( x , t ) ψ
where Δ g is the Laplace–Beltrami operator associated with the metric. In d = 1 , it is of course Δ x = x 2 .
A conformal diffeomorphism ( X , T , S ) f ˜ ( X , T , S ) = ( x , t , s ) with conformal factor σ f 2 , f ˜ * g μ ν = σ f 2 g μ ν , projects to a spacetime transformation ( X , T ) f ( X , T ) = ( x , t ) . It is implemented on a wave function lifted to the Bargmann space as:
ψ ˜ ( x , t , s ) = σ f 1 / 2 ψ ˜ ( X , T , S )
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 F ( t ) in 1 + 1 dimensions, (11),
x ¨ + λ ˙ x ˙ + ω 2 t x = F ( t ) ,
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,
T = u 1 u 2 , X = x u p u 2 ,
where u 1 and u 2 are solutions of the associated homogeneous Equation (46) with F 0 and u p 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,
u 1 t 0 = u ˙ 2 t 0 = 0 , u ˙ 1 t 0 = u 2 t 0 = 1 , u p t 0 = u ˙ p t 0 = 0 .
Then, in the new coordinates, the motion becomes free [22,23],
X ( T ) = a T + b , a , b = const .
Equation (46) can be obtained by projecting a null geodesic of the Bargmann metric:
g μ ν d x μ d x ν = e λ ( t ) d x 2 + 2 d t d s 2 e λ ( t ) 1 2 ω ( t ) 2 x 2 F ( t ) x d t 2 .
Completing (47) by:
S = s + e λ u 2 1 1 2 u ˙ 2 x 2 + u ˙ p x + g t where g ˙ = 1 2 e λ u ˙ p 2 ω 2 u p 2 + 2 F u p
lifts the Arnold map to Bargmann spaces, ( x , t , s ) ( X , T , S ) (In the Junker–Inomata setting (10), f = u 2 e λ / 2 and g ( t ) = u 1 / u 2 .),
g μ ν d x μ d x ν = e λ ( t ) u 2 2 ( t ) d X 2 + 2 d T d S .
The oscillator metric (50) is thus carried conformally to the free one, generalizing earlier results [18,19,41]. For the damped harmonic oscillator with λ ( t ) = λ 0 t and F ( t ) 0 , u p 0 is a particular solution. When ω = ω 0 = const . , for example,
u 1 = e λ 0 t / 2 sin Ω 0 t Ω 0 , u 2 = e λ 0 t / 2 cos Ω 0 t + λ 0 2 Ω 0 sin Ω 0 t , Ω 0 2 = ω 0 2 λ 0 2 / 4
are two independent solutions of the homogeneous equation with initial conditions (48) and provide us with:
T = sin Ω 0 t Ω 0 ( cos Ω 0 t + λ 0 2 Ω 0 sin Ω 0 t ) ,
X = e λ 0 t / 2 x cos Ω 0 t + λ 0 2 Ω 0 sin Ω 0 t ,
S = s 1 2 e λ 0 t x 2 ω 0 2 Ω 0 sin Ω 0 t cos Ω 0 t + λ 0 2 Ω 0 sin Ω 0 t .
In the undamped case, λ 0 = 0 ; thus, Ω 0 = ω 0 , and (56) reduces to that of Niederer [16] lifted to the Bargmann space [19,20],
T = tan ω 0 t ω 0 , X = x cos ω 0 t , S = s 1 2 x 2 ω 0 tan ω 0 t .
The Junker–Inomata construction in Section 2 can be viewed as a particular case of the Arnold transformation. We chose u p 0 and the two independent solutions:
u 1 = e λ / 2 τ ˙ 1 / 2 sin ω ¯ τ ω ¯ , u 2 = e λ / 2 τ ˙ 1 / 2 cos ω ¯ τ .
The initial conditions (48) at t 0 = 0 imply τ ( 0 ) = ρ ˙ ( 0 ) = 0 , ρ ( 0 ) = τ ˙ ( 0 ) = 1 . Then, spelling out (51),
S = s 1 2 e λ ω ¯ τ ˙ tan ω ¯ τ + 1 2 λ ˙ + 1 2 τ ¨ τ ˙ x 2
completes the lift of (22) to Bargmann spaces. In conclusion, the one-dimensional damped harmonic oscillator is described by the conformally flat Bargmann metric,
g μ ν d x μ d x ν = cos 2 ω ¯ τ τ ˙ d X 2 + 2 d T d S .
The metric (60) is manifestly conformally flat; therefore, its geodesics are those of the free metric, X ( T ) = a T + b . Then, using (47) with (58) yields:
x ( t ) = e λ ( t ) / 2 τ ˙ 1 / 2 ( t ) a sin [ ω ¯ τ ( t ) ] ω ¯ + b cos [ ω ¯ τ ( t ) ] .
The bracketed quantity here describes a constant-frequency oscillator with “time” τ ( t ) . 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, 0 < t t < π / Ω 0 . Its extension for all t involves the Maslov correction. In the constant-frequency case with no friction, for example, assuming that Ω 0 ( t t ) / π is not an integer, we have [2,3,6],
K e x t ( x , t | x , t ) = Ω 0 2 π | sin Ω 0 ( t t ) | 1 2 × e i π 4 ( 1 + 2 ) exp i Ω 0 2 sin Ω 0 ( t t ) ( x 2 + x 2 ) cos Ω 0 ( t t ) 2 x x ,
where the integer:
= Ent Ω 0 ( t t ) π
is called the Maslov index ( Ent [ x ] 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 2 A / x x [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
    K e x t ( x , t | x , t ) = ω ¯ e λ + λ 2 ( τ ˙ τ ˙ ) 1 2 2 π | sin ω ¯ ( τ τ ) | 1 / 2 exp i π 2 1 2 + Ent ω ¯ ( τ τ ) π × exp i ω ¯ 2 sin ω ¯ ( τ τ ) [ ( x 2 e λ τ ˙ + x 2 e λ τ ˙ ) cos [ ω ¯ ( τ τ ) ] 2 x x e λ + λ 2 ( τ ˙ τ ˙ ) 1 2 ] ;
  • At caustics, i.e., for:
    ω ¯ ( τ τ ) = π , = 0 , ± 1 ,
    we have instead [3,6],
    K e x t x , x , | τ τ = π ω ¯ = e λ + λ 2 ( τ ˙ τ ˙ ) 1 2 1 / 2 × exp i π 2 δ x exp ( λ / 2 ) τ ˙ 1 / 2 ( 1 ) k x exp ( λ / 2 ) τ ˙ 1 / 2 .
Proof. 
In terms of the redefined coordinates:
τ = τ ( t ) and ξ = x exp λ ( t ) 2 τ ˙ 1 / 2 ( t ) ,
cf. (37), and using the notation { · } = d/dτ, the time-dependent oscillator Equation (12) is taken into:
ξ + ω ¯ 2 ξ = 0 , where ω ¯ 2 = 1 τ ˙ 2 ( ω 2 ( t ) 1 2 τ τ ˙ + 3 4 ( τ ̈ τ ˙ ) 2 λ ̈ 2 λ ˙ 2 4 ) .
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:
ω 2 ( t ) τ ˙ 2 ω ¯ 2 = 1 2 S ( τ )
where S ( τ ) = τ τ ˙ 3 2 τ ¨ τ ˙ 2 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 ( ξ , τ ) ( x , t ) :
K 2 ( x , t | x , t ) = ξ x ξ x 1 2 K 1 ( ξ , τ | ξ , τ ) .
Here, K 2 = K e x t is the propagator of an oscillator with a time-dependent frequency and friction, ω ( t ) and λ ( t ) , respectively—the one we are trying to find. K 1 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 r k J k where sin = ± 1 . However, at caustics, τ τ = ( π / ω ¯ ) , K e x t diverges, and we have instead (66).
Henceforth, we limit our investigations to λ = 0 .

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, x = t = 0 and x x and t t .
We start with the observation that the Niederer map (22) becomes singular where the cosine vanishes, i.e., where:
cos [ ω ¯ τ ( r k ) ] = 0 , i . e . τ ( r k ) = ( k + 1 2 ) π ω ¯ , k = 0 , ± 1 ,
r k < r k + 1 because τ ( t ) is an increasing function by (21). Moreover, each interval:
I k = r k , r k + 1 , k = 0 , ± 1 ,
is mapped by (22) onto the full range < T < . Therefore, the inverse mapping is multivalued, labeled by integers k,
N k : T t = arctan k ω ¯ T ω ¯ , X x = X 1 + ω ¯ 2 T 2 ,
where arctan k ( · ) = arctan 0 ( · ) + k π with arctan 0 ( · ) the principal determination, i.e., in ( π / 2 , π / 2 ) .
Then, lim t r k tan t = and lim t r k + tan t = imply that:
lim T N k ( T ) = r k + 1 = lim T N k + 1 ( T ) .
Therefore, the intervals I k and I k + 1 are joined at r k + 1 , and the I k form a partition of the time axis, < t < = k I k .
Returning to (64) (which is (62) with Ω 0 ω ¯ , t τ ), we then observe that, whereas the propagator is regular at r k , it diverges at caustics,
sin [ ω ¯ τ ( t ) ] = 0 i . e . , τ ( t ) = π ω ¯ , = 0 , ± 1 , ,
cf. (65). Thus, t t + 1 , and:
N k ( ) = r k , N k ( T = 0 ) = t k + 1 , N k ( + ) = r k + 1 .
Thus, N k maps the full T-line into I k with t k an internal point. Conversely, r k is an internal point of J k . The intervals J = t , t + 1 cover again the time axis, J = < t < .
By (61) the classical trajectories are regular at t = r k . Moreover, for arbitrary initial velocities,
τ ˙ ( t + 1 ) x ( t + 1 ) = τ ˙ ( t ) x ( t )
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 N ˜ k : ( X , T , S ) ( x , t , s ) ( t I k ), is:
t = arctan k ω ¯ T ω ¯ , x = X 1 + ω ¯ 2 T 2 , s = S + X 2 2 ω ¯ 2 T 1 + ω ¯ 2 T 2 .
Proof. 
These formulae follow at once by inverting (57), at once with the cast ω 0 ω ¯ , t τ . Alternatively, it could also be proven as for Theorem 2.
For each integer k (78) maps the real line < T < into the “open strip” [19] r k , r k + 1 × R 2 I k × R 2 with r k 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 ω ¯ = 1 for simplicity. The free motion with initial condition X ( 0 ) = 0 ,
X ( T ) = a T , S ( T ) = S 0 a 2 2 T ,
E-D lifts by (78) to:
x ( t ) = a sin t s ( t ) = S 0 a 2 4 sin 2 t ,
consistent with s ( t ) = s 0 A ¯ o s c , as can be checked directly. Note that the s coordinate oscillates with a doubled frequency.
  • At t = r k = ( 1 2 + k ) π (where the Niederer maps are joined), we have lim t r k x ( t ) = ( 1 ) k + 1 a , lim t r k s ( t ) = S 0 . Thus, the pull backs of the Bargmann lifts of free motions are glued to smooth curves;
  • Similarly, at t caustics t = t = π , we infer from (80) that for all initial velocities a and for all   lim t t x ( t ) = 0 , lim t t s ( t ) = S 0 . Thus, the lifts are again smooth at t , and after each half-period, all motions are focused above the initial position ( x ( 0 ) = 0 , s ( 0 ) = S 0 ) .

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 0 < T T . Its extension to all T involves the sign of ( T T ) [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 i = e i π / 4 , then the prefactor in (5) is:
1 2 π i ( T T ) 1 / 2 = e i π / 4 1 2 π | T T | 1 / 2 .
Let us now consider what happens when T T changes sign. Then, the prefactor becomes multiplied by 1 so it becomes, for the same choice of the square root,
e i π / 2 e i π / 4 1 2 π | T T | 1 / 2 = e + i π / 4 1 2 π | T T | 1 / 2 .
In conclusion, the formula valid for all T is,
K f r e e ( X , T | X , T ) = e i π 4 sign ( T T ) 1 2 π | T T | 1 / 2 exp i A ¯ f r e e ,
where:
A ¯ f r e e = ( X X ) 2 2 ( T T )
is the free action calculated along the classical trajectory. Let us underline that (82) already involves a “Maslov jump” e i π / 2 , which, for a free particle, happens at T = 0 . For T T = 0 , we have K f r e e = δ ( X X ) .
Accordingly, the wave function Ψ Ψ f r e e of a free particle is, by (1),
Ψ X , T = e i π 4 sign T T 1 2 π | T T | 1 / 2 R exp i A ¯ f r e e Ψ X , T d X .
Now, we pull back the free dynamics using the multivalued inverse Niederer map. It is sufficient to consider the constant-frequency case ω ¯ = const . and to denote time by t. Let t belong to the range of N k in (73), t I k = [ r k , r k + 1 ] = N k { < T < } . Then, applying the general formulae in Section 3.1 yields [19],
ψ ˜ ( x , t , s ) = cos 1 / 2 [ ω ¯ ( t t ) ] Ψ ˜ ( X , T , S ) = e i π 4 sign tan ω ¯ ( t t ) ω ¯ × cos 1 / 2 ω ¯ ( t t ) exp i s exp i ( 1 2 ω ¯ x 2 tan [ ω ¯ ( t t ) ] | ω ¯ | 2 π | tan [ ω ¯ ( t t ) ] | R exp i ω ¯ | x cos [ ω ¯ ( t t ) ] x | 2 2 tan [ ω ¯ ( t t ) ] ψ ( x , t ) d x .
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,
A ¯ o s c = ω ¯ 2 sin ω ¯ ( t t ) ( x 2 + x 2 ) cos ω ¯ ( t t ) 2 x x .
Thus, using the equivariance, we end up with,
ψ o s c x , t = cos 1 / 2 [ ω ¯ ( t t ) ] exp i π 4 sign tan [ ω ¯ ( t t ) ] ω ¯ × | ω ¯ | 2 π | tan [ ω ¯ ( t t ) ] | R exp i A ¯ o s c ψ o s c x , t d x .
Now, we recover the Maslov jump, which comes from the first line here. For simplicity, we consider again t = 0 , x = 0 and denote t = t , x = x .
Firstly, we observe that the conformal factor cos ω ¯ t has a constant sign in the domain I k and changes sign at the end points. In fact,
cos ω ¯ t = 1 k + 1 | cos ω ¯ t | cos 1 / 2 ω ¯ t = e i π 2 ( k + 1 ) | cos ω ¯ t | 1 / 2 .
The cosine enters into the van Vleck factor, while the phase combines with exp i π 4 sign ( tan ω ¯ t ω ¯ ) . Recall now that t k + 1 = N k T = 0 divides I k into two pieces, I k = [ r k , t k + 1 ] [ t k + 1 , r k + 1 ] , cf. Figure 1. However, t k + 1 is precisely where the tangent changes sign: this term contributes to the phase in [ r k , t k + 1 ] π / 4 and + π / 4 in [ t k + 1 , r k + 1 ] . Combining the two shifts, we end up with the phase:
π 4 1 + 2 for r k < t < t k + 1 π 4 1 + 2 ( + 1 ) for t k + 1 < t < r k + 1 where = Ent ω ¯ τ π = k + 1
which is the Maslov jump at t .
Intuitively, the multivalued N k “exports” to the oscillator at t + 1 the phase jump of the free propagator at T = 0 . Crossing from J to J + 1 shifts the index by one.

5. Probability Density and Phase of the Propagator: A Pictorial View

5.1. For a Constant Frequency

We assume first that the frequency is constant. We split the propagator K ( x , t ) K ( x , t | 0 , 0 ) in (62) as,
K ( x , t ) = | K ( x , t ) | P ( t ) , P ( t ) = e i ( p h a s e ) .
The probability density,
| K ( x , t ) | 2 = Ω 0 2 π | sin Ω 0 t |
viewed as a surface above the x t plane, diverges at t = t = π , = 0 , ± 1 , .
Representing the phase of the propagator would require four dimensions, though. However, recall that that the dominant contribution to the path integral should come from where the phase is stationary [1], i.e., from the neighborhood of classical paths x ¯ ( t ) , distinguished by the vanishing of the first variation, δ A x ¯ = 0 . Therefore, we shall study the evolution of the phase along classical paths x ¯ ( t ) for which (61) yields, for = ω ¯ = 1 and a R , b = 0 ,
x ¯ a ( t ) = a sin t and P a ( t ) = exp i π 4 1 a 2 π sin 2 t i π 2 ,
as depicted in Figure 2.
An intuitive understanding comes by noting that when t π = t , then different initial velocities a yield classical paths x ¯ a ( t ) with different end points, and thus contribute to different propagators. However, approaching from the left -times a half period, t ( π ) , all classical paths become focused at the same end point ( x = 0 for our choice) and for all a,
P a ( t π ) = e i π 4 ( 1 + 2 ) P .
which is precisely the Maslov phase. Thus, all classical paths contribute equally, by P , and to the same propagator. Comparing with the right-limit,
P a ( t π + ) = e i π 4 ( 1 + 2 ( + 1 ) = P + 1 = e i π 2 P .
the Maslov jump is recovered. Choosing instead y 0 , there will be no classical path from ( 0 , 0 ) to ( y , π ) , and thus no contribution to the path integral.
To conclude this section, we just mention with that the extended Feynman method [6] with the cast ω ¯ = constant frequency and τ = “fake time” would lead also to (64) and (66) with the integer counting the number of negative eigenvalues (Morse index) of the Hessian [2,4,5,24].

5.2. A time-Dependent Example: The Mathieu Equation

The combined Junker–Inomata–Arnold method allows us to go beyond the constant-frequency case, as illustrated here for no friction or driving force, λ = F 0 , but with explicitly time-dependent frequency. For Ω 2 ( t ) = a 2 q cos 2 t , for example, (14) becomes the Mathieu equation,
x ¨ + ( a 2 q cos 2 t ) x = 0 .
This equation can be solved either analytically using Mathieu functions [43], or numerically, providing us for a = 2 and q = 1 (for which odd Mathieu functions are real) with the dotted curve (in red), shown in Figure 3.
Alternatively, we can use the Junker–Inomata–Arnold transformation (47) [22,23,40]. We first achieve ω ¯ = 1 by a redefinition, τ τ = ω ¯ τ . Inserting Ansatz (15) into (94) yields the pair of coupled Equations (16) and (17). We chose u p = 0 and two independent solutions u 1 ( t ) and u 2 ( t ) , (58), with initial conditions (48) with t 0 = 0 , i.e., τ ( 0 ) = ρ ˙ ( 0 ) = 0 , ρ ( 0 ) = τ ˙ ( 0 ) = 1 , which fix the integration constant, C = ρ 2 0 τ ˙ ( 0 ) = 1 . Then, consistent with the general theory outlined above, the Arnold map (47) lifted to the Bargmann space becomes (22), completed with (59) with λ = 0 .
Equation (17) is solved by following the strategy outlined in Section 2. Carrying out those steps numerically provides us with Figure 3.
From the general formula (24), we deduce, for our choice x = x , t = t , x = t = 0 , that the probability density (The wave function is multiplied by the square root of the conformal factor, cf. (39).).
| K ( x , t ) | 2 = τ ˙ 2 π | sin τ ( t ) | ,
happens, not depending on the position, and can therefore be plotted as in Figure 4.
The propagator K and hence the probability density (95) diverge at t , which are roughly t 1 1.92 , t 2 4.80 , t 3 7.83 . The classical motions are regular at the caustics, x ¯ ( t ) ρ ( t ) 0 ; see Section 4. The domains I k = [ r k 1 , r k ] of the inverse Niederer map are shown in Figure 4. Approximately, r 1 1.52 , r 2 4.49 , r 3 6.75 , r 4 8.44 . The evolution of the phase factor along the classical path is depicted in Figure 5.

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,
ω ( t ) ω ¯ = const . and t fake 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 ω ( t ) : 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 τ = τ ( t ) following the procedure outlined in Section 2. We have to first solve EMP Equation (18) for ρ ( t ) (which is nonlinear and has time-dependent coefficients) and then integrate ρ 2 ; 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 r k in (71) and the t in (75). The r k divides the time axis into domains I k of the (generalized) Niederer map (22). Both classical motions and quantum propagators are regular at r k , where these intervals are joined. The t 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, ω = ω 0 = const . , its extension to the time-dependent case ω = ω ( t ) 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 ω ( t ) . 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 T T induces a phase jump π / 2 . Our “innocent-looking” factor is in fact the Maslov jump for a free particle at T = 0 (obscured when one considers the propagator for T > 0 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].

Author Contributions

All authors have contributed equally and substantially to this paper, including its conception, development, preparation, writing and editing. All authors have read and agreed to publish the manuscript in its present form.

Funding

This research was funded by the National Natural Science Foundation of China (grant number 11975320).

Institutional Review Board Statement

Exclude this statement since the study did not involve humans or animals.

Informed Consent Statement

Exclude this statement since the study did not involve humans.

Data Availability Statement

Exclude this statement since the study did not report any data.

Acknowledgments

This work was partially supported by the National Natural Science Foundation of China (Grant No. 11975320).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Feynman, R.P.; Hibbs, A.R. Quantum Mechanics and Path Integrals; McGraw-Hill: New York, NY, USA, 1965. [Google Scholar]
  2. Schulman, L. Techniques and Applications of Path Integration; Wiley: New York, NY, USA, 1981. [Google Scholar]
  3. Khandekar, D.C.; Lawande, S.V.; Bhagwat, K.V. Path-Integral Methods and Their Applications, 1st ed.; World Scientific: Singapore, 1993. [Google Scholar]
  4. DeWitt-Morette, C. The Semiclassical Expansion. Ann. Phys. 1976, 97, 367–399. [Google Scholar] [CrossRef]
  5. Levit, S.; Smilansky, U. A New Approach to Gaussian Path Integrals and the Evaluation of the Semiclassical Propagator. Ann. Phys. 1977, 103, 198. [Google Scholar] [CrossRef]
  6. Horvathy, P.A. Extended Feynman Formula for Harmonic Oscillator. Int. J. Theor. Phys. 1979, 18, 245. [Google Scholar] [CrossRef]
  7. Caldirola, P. Forze non-conservative nella meccanica quantistica. Nuovo Cim. 1941, 18, 393. [Google Scholar] [CrossRef]
  8. Kanai, E. On the Quantization of the Dissipative Systems. Prog. Theor. Phys. 1948, 3, 440. [Google Scholar] [CrossRef]
  9. Dekker, H. Classical and quantum mechanics of the damped harmonic oscillator. Phys. Rep. 1981, 80, 1–112. [Google Scholar] [CrossRef]
  10. Khandekar, D.C.; Lawande, S.V. Feynman Path Integrals: Some Exact Results and Applications. Phys. Rep. 1986, 137, 115–229. [Google Scholar] [CrossRef]
  11. Um, C.I.; Yeon, K.H.; George, T.F. The Quantum damped harmonic oscillator. Phys. Rep. 2002, 362, 63–192. [Google Scholar] [CrossRef]
  12. Junker, G.; Inomata, A. Transformation of the free propagator to the quadratic propagator. Phys. Lett. A 1985, 110, 195–198. [Google Scholar] [CrossRef]
  13. Cai, P.Y.; Inomata, A.; Wang, P. Jackiw Transformation in Path Integrals. Phys. Lett. A 1982, 91, 331–334. [Google Scholar] [CrossRef]
  14. Cai, P.Y.; Cai, J.M.; Inomata, A. A. A time-dependent conformal transformation in Feynman’s path integral. In Path integrals from meV to MeV; World Scientific: Teaneck, NJ, USA, 1989; pp. 279–290. [Google Scholar]
  15. Inomata, A. Time-dependent conformal transformation in quantum mechanics. In Proceedings of the ISATQP-Shanxi 1992, Taiyuan, China, 12–16 June 1992; Liang, J.Q., Wang, M., Qiao, S.N., Su, D.C., Eds.; Science Press: Beijing, China, 1993; pp. 75–82. [Google Scholar]
  16. Niederer, U. The maximal kinematical invariance group of the harmonic oscillator. Helv. Phys. Acta 1973, 46, 191–200. [Google Scholar]
  17. Eisenhart, L.P. Dynamical trajectories and geodesics. Ann. Math. 1928, 30, 591–606. [Google Scholar] [CrossRef]
  18. Duval, C.; Burdet, G.; Kunzle, H.P.; Perrin, M. Bargmann Structures and Newton-cartan Theory. Phys. Rev. D 1985, 31, 1841–1853. [Google Scholar] [CrossRef] [PubMed]
  19. Burdet, G.; Duval, C.; Perrin, M. Time Dependent Quantum Systems and Chronoprojective Geometry. Lett. Math. Phys. 1985, 10, 255–262. [Google Scholar] [CrossRef]
  20. Duval, C.; Gibbons, G.W.; Horvathy, P. Celestial mechanics, conformal structures and gravitational waves. Phys. Rev. D 1991, 43, 3907–3922. [Google Scholar] [CrossRef] [Green Version]
  21. Cariglia, M.; Duval, C.; Gibbons, G.W.; Horvathy, P.A. Eisenhart lifts and symmetries of time-dependent systems. Ann. Phys. 2016, 373, 631–654. [Google Scholar] [CrossRef] [Green Version]
  22. Arnold, V.I. Supplementary Chapters to the Theory of Ordinary Differential Equations; Nauka: Moscow, Russia, 1978. [Google Scholar]
  23. Arnold, V.I. Geometrical Methods in the Theory of Ordinary Differential Equations; Springer: New York, NY, USA, 1983. (In English) [Google Scholar]
  24. Maslov, V.P.; Bouslaev, V.C.; Arnol’d, V.I. Théorie des Perturbations et Méthodes Asymptotiques; Dunod: Paris, France, 1972. [Google Scholar]
  25. Arnold, V.I. Characteristic class entering in quantization conditions. Funktsional’Nyi Anal. Ego Prilozheniya 1967, 1, 1–14. [Google Scholar] [CrossRef]
  26. Souriau, J.M. Construction explicite de l’indice de Maslov. Applications. Lect. Notes Phys. 1976, 50, 117–148. [Google Scholar] [CrossRef]
  27. Burdet, G.; Perrin, M.; Perroud, M. Generating functions for the affine symplectic group. Comm. Math. Phys. 1978, 58, 241–254. [Google Scholar] [CrossRef]
  28. Rezende, J. Quantum Systems with Time Dependent Harmonic Part and the Morse Index. J. Math. Phys. 1984, 25, 32643269. [Google Scholar] [CrossRef]
  29. Aldaya, V.; Cossío, F.; Guerrero, J.; López-Ruiz, F.F. The quantum Arnold transformation. J. Phys. A 2011, 44, 065302. [Google Scholar] [CrossRef]
  30. Guerrero, J.; López-Ruiz, F.F.; Aldaya, V.; Cossío, F.C. Symmetries of the quantum damped harmonic oscillator. J. Phys. A 2012, 45, 475303. [Google Scholar] [CrossRef] [Green Version]
  31. Guerrero, J.; Aldaya, V.; López-Ruiz, F.F.; Cossío, F. Unfolding the quantum Arnold transformation. Int. J. Geom. Meth. Mod. Phys. 2012, 9, 1260011. [Google Scholar] [CrossRef]
  32. Ermakov, V.P. Second order differential equations. Conditions of complete integrability. Univ. Izv. Kiev Series III 1880, 9, 1. [Google Scholar] [CrossRef] [Green Version]
  33. Milne, W.E. The numerical determination of characteristic numbers. Phys. Rev. 1930, 35, 863. [Google Scholar] [CrossRef] [Green Version]
  34. Pinney, E. The nonlinear differential equation y′′+p(x)y+cy3=0. Proc. Am. Math. Soc. 1959, 1, 68. [Google Scholar]
  35. Galajinsky, A. Geometry of the isotropic oscillator driven by the conformal mode. Eur. Phys. J. C 2018, 78, 72. [Google Scholar] [CrossRef] [Green Version]
  36. Cariglia, M.; Galajinsky, A.; Gibbons, G.W.; Horvathy, P.A. Cosmological aspects of the Eisenhart–Duval lift. Eur. Phys. J. C 2018, 78, 314. [Google Scholar] [CrossRef] [Green Version]
  37. Cheng, B. Exact propagator for the harmonic oscillator with time dependent mass. Phys. Lett. A 1985, 113, 293. [Google Scholar] [CrossRef]
  38. Ilderton, A. Screw-symmetric gravitational waves: A double copy of the vortex. Phys. Lett. B 2018, 782, 22–27. [Google Scholar] [CrossRef]
  39. Zhang, P.M.; Cariglia, M.; Elbistan, M.; Horvathy, P.A. Scaling and conformal symmetries for plane gravitational waves. J. Math. Phys. 2020, 61, 022502. [Google Scholar] [CrossRef] [Green Version]
  40. López-Ruiz, F.F.; Guerrero, J. Generalizations of the Ermakov system through the Quantum Arnold Transformation. J. Phys. Conf. Ser. 2014, 538, 012015. [Google Scholar] [CrossRef] [Green Version]
  41. Duval, C.; Horváthy, P.A.; Palla, L. Conformal properties of Chern-Simons vortices in external fields. Phys. Rev. 1994, D50, 6658. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Gibbons, G.W. Dark Energy and the Schwarzian Derivative. arXiv 2014, arXiv:1403.5431. [Google Scholar]
  43. Weisstein, E.W. Mathieu Function. 2003. Available online: https://mathworld.wolfram.com/MathieuFunction.html (accessed on 1 August 2021).
  44. Cheng, B. Exact propagator for the one-dimensional time-dependent quadratic Lagrangian. Lett. Math. Phys. 1987, 14, 7–13. [Google Scholar] [CrossRef]
  45. Jackiw, R. Introducing scale symmetry. Phys. Today 1972, 25N1, 23–27. [Google Scholar] [CrossRef]
  46. Niederer, U. The maximal kinematical invariance group of the free Schrodinger equation. Helv. Phys. Acta 1972, 45, 802–810. [Google Scholar] [CrossRef]
  47. Hagen, C.R. Scale and conformal transformations in galilean-covariant field theory. Phys. Rev. D 1972, 5, 377–388. [Google Scholar] [CrossRef]
  48. De Alfaro, V.; Fubini, S.; Furlan, G. Conformal Invariance in Quantum Mechanics. Nuovo Cim. A 1976, 34, 569. [Google Scholar] [CrossRef] [Green Version]
  49. Jackiw, R. Dynamical Symmetry of the Magnetic Vortex. Ann. Phys. 1990, 201, 83–116. [Google Scholar] [CrossRef]
  50. Jackiw, R. Dynamical Symmetry of the Magnetic Monopole. Ann. Phys. 1980, 129, 183. [Google Scholar] [CrossRef]
  51. Andrzejewski, K.; Prencel, S. Memory effect, conformal symmetry and gravitational plane waves. Phys. Lett. B 2018, 782, 421–426. [Google Scholar] [CrossRef]
  52. Andrzejewski, K.; Prencel, S. Niederer’s transformation, time-dependent oscillators and polarized gravitational waves. Class. Quantum Gravity 2019, 36, 155008. [Google Scholar] [CrossRef] [Green Version]
  53. Inzunza, L.; Plyushchay, M.S.; Wipf, A. Conformal bridge between asymptotic freedom and confinement. Phys. Rev. D 2020, 101, 105019. [Google Scholar] [CrossRef]
  54. Dhasmana, S.; Sen, A.; Silagadze, Z.K. Equivalence of a harmonic oscillator to a free particle and Eisenhart lift. arXiv 2021, arXiv:2106.09523. [Google Scholar]
  55. Guha, P.; Garai, S. Integrable modulation, curl forces and parametric Kapitza equation with trapping and escaping. arXiv 2021, arXiv:2104.06319. [Google Scholar]
Figure 1. The generalized Niederer map (22) maps each interval I k = ( r k , r k + 1 ) onto the entire real line < T < . Its inverse mapping is therefore multivalued, labeled by an integer k. The classical motions and the propagator are both regular at the separation points r k . All classical trajectories are focused at the caustic points t , where the propagator diverges.
Figure 1. The generalized Niederer map (22) maps each interval I k = ( r k , r k + 1 ) onto the entire real line < T < . Its inverse mapping is therefore multivalued, labeled by an integer k. The classical motions and the propagator are both regular at the separation points r k . All classical trajectories are focused at the caustic points t , where the propagator diverges.
Symmetry 13 01866 g001
Figure 2. The phase factor P ( t ) of the propagator in (89) lies on the unit circle of the complex plane plotted vertically along a classical path γ ¯ ( t ) . The orientation is positive if it is clockwise when seen from t = + . In the time interval J labeled by the Maslov index = Ent [ t / π ] , the factor P ( t ) precesses around P = exp [ i π 4 ( 1 + 2 ) ] with double frequency w.r.t. the classical path, γ ¯ ( t ) . Arriving at a caustic, the phase jumps by ( π / 2 ) (red becoming purple) and then continues until the next caustic when it jumps again (and becomes magenta), and so on.
Figure 2. The phase factor P ( t ) of the propagator in (89) lies on the unit circle of the complex plane plotted vertically along a classical path γ ¯ ( t ) . The orientation is positive if it is clockwise when seen from t = + . In the time interval J labeled by the Maslov index = Ent [ t / π ] , the factor P ( t ) precesses around P = exp [ i π 4 ( 1 + 2 ) ] with double frequency w.r.t. the classical path, γ ¯ ( t ) . Arriving at a caustic, the phase jumps by ( π / 2 ) (red becoming purple) and then continues until the next caustic when it jumps again (and becomes magenta), and so on.
Symmetry 13 01866 g002
Figure 3. The analytic solution of the Mathieu equation with a = 2 , q = 1 for x ( t ) (dotted in red) lies on the black curve obtained by (15) from combining the numerically obtained ρ ( t ) (in green) and τ ( t ) (in blue), which are solutions of the pair (18)–(21). The black curve is also obtained by pulling back the free solution (49) by the inverse Niederer map (73).
Figure 3. The analytic solution of the Mathieu equation with a = 2 , q = 1 for x ( t ) (dotted in red) lies on the black curve obtained by (15) from combining the numerically obtained ρ ( t ) (in green) and τ ( t ) (in blue), which are solutions of the pair (18)–(21). The black curve is also obtained by pulling back the free solution (49) by the inverse Niederer map (73).
Symmetry 13 01866 g003
Figure 4. The probability density | K ( x , t ) | 2 (95) does not depend on x and is regular in each interval J between the adjacent points t (75), where it diverges. The r k that determines the domains I k of the generalized Niederer map (22) lies between the t and conversely.
Figure 4. The probability density | K ( x , t ) | 2 (95) does not depend on x and is regular in each interval J between the adjacent points t (75), where it diverges. The r k that determines the domains I k of the generalized Niederer map (22) lies between the t and conversely.
Symmetry 13 01866 g004
Figure 5. For 0 < t < t 1 , the Mathieu phase factor P ( t ) plotted along a classical path γ ¯ ( t ) = ( x ¯ ( t ) , t ) precesses around e i π / 4 . Arriving at the caustic point τ ( t 1 ) = π , its phase jumps by ( π / 2 ) , then oscillates around e 3 i π / 4 until τ ( t 2 ) = 2 π , then jumps again, and so on.
Figure 5. For 0 < t < t 1 , the Mathieu phase factor P ( t ) plotted along a classical path γ ¯ ( t ) = ( x ¯ ( t ) , t ) precesses around e i π / 4 . Arriving at the caustic point τ ( t 1 ) = π , its phase jumps by ( π / 2 ) , then oscillates around e 3 i π / 4 until τ ( t 2 ) = 2 π , then jumps again, and so on.
Symmetry 13 01866 g005
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zhao, Q.; Zhang, P.; Horvathy, P.A. Time-Dependent Conformal Transformations and the Propagator for Quadratic Systems. Symmetry 2021, 13, 1866. https://doi.org/10.3390/sym13101866

AMA Style

Zhao Q, Zhang P, Horvathy PA. Time-Dependent Conformal Transformations and the Propagator for Quadratic Systems. Symmetry. 2021; 13(10):1866. https://doi.org/10.3390/sym13101866

Chicago/Turabian Style

Zhao, Qiliang, Pengming Zhang, and Peter A. Horvathy. 2021. "Time-Dependent Conformal Transformations and the Propagator for Quadratic Systems" Symmetry 13, no. 10: 1866. https://doi.org/10.3390/sym13101866

APA Style

Zhao, Q., Zhang, P., & Horvathy, P. A. (2021). Time-Dependent Conformal Transformations and the Propagator for Quadratic Systems. Symmetry, 13(10), 1866. https://doi.org/10.3390/sym13101866

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