1. Introduction
Unlike many energy-conserving oscillatory systems, the relativistic harmonic oscillator (RHO) represents a unique case. In many of these systems, the nonlinearity is a nonlinear function of the amplitude, x(t). Usually, the occurrence in the nonlinearity of the velocity,
, represents damping and energy non-conservation. The RHO is an energy-conserving system, in which the nonlinearity involves the velocity. Furthermore, in popular systems, as the amplitude of oscillation, a, tends to infinity, the period of oscillations, T, tends to zero and the maximal velocity tends to infinity. In the case of the RHO, the period tends to infinity and the velocity tends to a constant, the speed of light, c.
The study of relativistic oscillations has progressed in several avenues: The exploitation of relativistic oscillators for the purpose of the description of bound states of quarks in High-Energy Physics [1] [2] [3] [4]; The study of the Dirac oscillator theoretically [5] - [12] and experimentally [13] [14], and the Quantum-Mechanical RHO [15] - [23]. Oscillations in the weak-relativistic limit have been of interest in Plasma Physics [24] - [36]. The emergence of chaotic behavior in a parametrically forced relativistic van der Pol oscillator is discussed in Ref. [37]. The behavior of an ensemble of RHO’s is studied in Ref. [38]. Interest in the RHO, in both weak- and extreme-relativistic limits has won renewed attention in the study of media, through which light propagates at velocities that are much lower than its speed in vacuum. A recent example is that of the dynamics of an oscillator composed of ultra-cold Bose-condensed atoms [39] [40] [41], where the effective speed of light is lower than in vacuum by orders of magnitude, e.g., 143 mm/s [41]. The energy-momentum relation of the atoms is nearly identical to that of a massive relativistic particle, with an effective mass.
Over many years, the major effort in the analysis has been dedicated to either generating numerical solutions or to various approximation schemes in the weak relativistic limit by expansions in powers of β2 (
) for
[42] - [50]. Owing to the emergence of secular terns in asymptotic expansion methods, the validity in time of approximate solutions is limited. Typically, the lowest-order approximation (O(β2)) in such analyses has a 1% - 2% error relative to the numerical solution up to only β ≈ 0.2. One important result of Ref. [43] is the presentation the period of oscillations in the frame of reference of an outside observer in terms of complete elliptic integrals.
Qualitative progress has been attained through the exploitation of the transformation from the time variable in the observer’s frame of reference to that of the oscillating mass (proper time) [51] - [57]. The oscillation period in proper time was derived in terms of a complete elliptic integral in Ref. [51]. The equivalence of the equation of motion of the RHO with a specific form of the equation of motion of a Duffing equation was found in Refs. [52] [53] [55] [57]. This has led to the derivation of a closed-from solution for the equation of the RHO in terms of Jacobi elliptic functions [51] [56].
The purpose of his paper is to expose the detailed characteristics of the solution of the equation of motion of the RHO in the weak- and extreme-relativistic limits.
Preliminary information is presented in Section 2. Exploiting the known closed-form expressions of the period, T, and the oscillation amplitude, xmax, the time and position variables, t and x, are transformed in Section 3 into dimensionless variables, θ and η, respectively. η is (2π)-periodic in θ and bounded: |η| ≤ 1. This transformation relieves the analysis in the weak-relativistic limit of the restriction on the validity in time of approximate solutions, encountered in standard asymptotic expansion perturbative analyses of the weak-relativistic limit. In addition, it naturally leads to a novel approach to the analysis of the extreme-relativistic limit.
The weak-relativistic limit is discussed in Section 4 through a Normal Form expansion [58] - [66]. As the solution is (2π)-periodic in θ, there is no need to worry about the growth of secular errors in the approximate solution for
. A direct consequence is that the error in the approximate solution turns out to remain very small up to rather high values of β. For example, in an expansion through O(β2), the approximations for the dimensionless position, η, and for the velocity, η', agree with the full numerical solution within a fraction of 1% up to β ≈ 0.6 and β ≈ 0.45, respectively.
The expansion in powers of β2 fails in the extreme relativistic limit (β 1). In this limit, over most of the period of oscillations, the velocity
,oscillates between two constant values, ± vmax, (vmax = c·β) and is reduced to zero rapidly over a short time interval near the turning point. The analysis is presented in Section 5 through the introduction of a small parameter, μ, which measures the proximity to the extreme relativistic limit.
Near the turning point (θ = 0), the velocity varies rapidly over a range of O(μ2) in θ. Hence, for finding an approximate solution for η over this narrow range in the extreme-relativistic limit, a variable transformation is required,
, so that z is of O(1). The resulting expansion of the solution in powers of μ2, expected to provide a good approximation to the solution near the turning point (θ = O(μ2) (z = O(1))), generates an excellent approximation to the solution for 0 ≤ θ ≤ 2π, even in the lowest order in μ2, for small μ. Unique features of the Fourier expansion of the solution are discussed in Section 6. Some concluding comments are presented in Section 7.
2. Preliminaries
The dynamical equation of the RHO is:
, (1)
where m is the mass of the oscillation particle, ω0 is the period of the oscillator in the non-relativistic limit and c is the speed of light.
The total energy is given by:
. (2)
Solving Equation (2) for
yields the period of oscillations, T [43]:
. (3)
In Equation (3), E(k) and K(k) are the complete elliptic integrals [67]:
. (4)
The maximal values of the amplitude of oscillations and of the velocity are obtained by setting
and x = 0, respectively, in Equation (2), yielding:
. (5)
In the limit of small relativistic effects (
), the expansion of E, T (using expansions of elliptic functions presented in [67]) and xmax in powers of β2 is given by:
(6)
In the far relativistic limit (
) it turns out to be convenient to define a new small parameter:
. (7)
Equation (6) is then replaced by (the expansion of the elliptic integrals can be found in Ref. [67]):
(8)
Note that both T and xmax tend to infinity as E1/2.
Figure 1 shows the phase-space plot of the numerical solution of Equation (1) for a sequence of values of β. x and t are expressed in terms of dimensionless entities:
. (9)
Figure 1. Phase space plots of solution of Equation (1) in terms of dimensionless entities (Equation (9)). From inner to outer plot: β = 0.1, 0.5, 0.999.
At high values of β, the phase-space plot tends to a rectangle. This has been observed in the many numerical analyse cited above.
3. Scaled Equation
Rescaling x and t as follows:
, (10)
Equation (1) becomes
. (11)
η(θ) is (2π)-periodic in θ and bounded: |η(θ)| ≤ 1. Exploiting Equations (10) and (6), yields that the maximal values of η and η' are both finite for all 0 ≤ β2 ≤ 1:
(12)
Figures 2-4 show the dependence of the solution on β. For
, the phase-space plot is very close to a circle. As β is increased, it approaches a rectangle. Concurrently, η' evolves from a sine-like function of θ to a periodic step function and that of η – to a saw-tooth pattern.
For
, Equation (11) is reduced to the equation of a classical harmonic oscillator perturbed by a small nonlinear term:
Figure 2. Phase-space plots, numerical solution of Equation (11).
Figure 3. η(θ) vs. θ, numerical solution of Equation (11).
Figure 4. η'(θ) vs. θ, numerical solution of Equation (11). Horizontal lines are maximal values of η'.
. (13)
Hence, this calls for an expansion in powers of β2.
In the extreme-relativistic limit, Equation (11), using Equations (7) and (8), and expanding in powers of
:
. (14)
4. Weak-Relativistic Limit
A convenient way to analyze Equation (11) is to transform it to a first-order equation for a complex entity:
. (15)
Expressing η and η' in terms of z and z*, Equation (11) is transformed into:
(16)
The boundary conditions at the turning point are:
. (17)
One now writes z(θ) as:
(18)
and exploits the expansion of T and xmax (Equation (6)) to obtain an expansion of Equation (16) in powers of β2.
The zero-order term of this expansion yields:
. (19)
Thus, as expected, the zero-order term is (2π)-periodic in θ. The O(β2) term in the expansion of Equation (16) is:
. (20)
For the solution to be periodic in θ, the coefficient of the resonant term cos(θ + θ0) must vanish, yielding
. (21)
One now proceeds in the same fashion through higher orders, requiring that the higher-order corrections, z4 and z6 of Equation (18) are periodic in θ (i.e., the coefficients of resonant terms in the dynamical equations for z4 and z6 vanish) and that the initial conditions are obeyed by the zero-order term:
. (22)
The result is the following expression for η:
(23)
The fact that the period, T, is known explicitly means that the restriction of the validity of the approximation that emerges in well-known asymptotic expansion methods is avoided. In these methods, the period is approximated by a truncated expansion in powers of the small parameter, limiting the validity of the approximation to times of O(1) or O(1/β2), depending on the analysis method used. This limitation is avoided here. The approximation to the full numerical solution up the rather high values of β is excellent. For β = 0.1, the maximal relative error generated by the approximate solution of Equation (23) is, at most, O(.01%). In expansions through O(β2), O(β4) and O(β6), the approximation for η agrees with the full numerical solution of Equation (11) within a fraction of 1% up to β ≈ 0.6, 0.7 and 0.8, respectively, for all θ. The approximation for η', agrees with the full numerical solution within a fraction of 1% up to β ≈ 0.4, 0.55 and 0.65, respectively, for all θ. Even for β ≈ 0.8, the relative error of the approximate expression for the amplitude (through O(β6), Equation (23)) varies between 0 and close to 1% (see Figure 5) and for the velocity—between 1% and 5% (see Figure 6).
5. Extreme Relativistic Limit
In the extreme relativistic limit, the expansion in powers of β2 fails. Figures 7(a)-(c) show how the approximate solution (Equation (23)) begins to deviate
Figure 5. Error in approximate solution for η(θ) Equation (23) relative to numerical solution of Equation (11), β = 0.8.
Figure 6. Error in approximate velocity η'(θ) computed from Equation (23) relative to numerical solution of Equation (11), β = 0.8.
Figure 7. Phase-space plots showing evolution of deviation of approximate solution (through O(β6), Equation (23)) from numerical solution of Equation (11) for high β. Continuous: numerical solution; Dashed: Approximate solution.
from the numerical solution of Equation (11) as β grows. The errors in η and η' are shown in Figure 8 and Figure 9, respectively, in a computation through O(β6), for β = 0.99. The error in η is, at most about 10%. In η', it exceeds 60% near the turning point.
In this limit, the velocity,
, is constant at (±c β) most of the time. It vanishes rapidly over a short time gap near the turning point (See Figure 1). In the scaled coordinates, the dependence of η on θ is a straight line, except very close to the turning point (see Figure 3).
To guide oneself towards a choice of a scaled variable that is appropriate for the analysis of Equation (11) near the turning point, one uses
of Equation (7), which is a measures of the proximity to the extreme-relativistic limit. Observe that η' is now very close to zero. Using Equation (8) for T and xmax and replacing η' by zero in Equation (14), the O(μ2) approximation to the equation becomes:
Figure 8. Error in approximate solution for η(θ) (Equation (23)) relative to numerical solution of Equation (11), β = 0.99.
Figure 9. Error in approximate velocity η'(θ) computed from Equation (23) relative to numerical solution of Equation (11), β = 0.99.
. (24)
For the boundary conditions at the turning point (θ = 0), η(0) = 1 and η'(0) = 0, Equation (24) is solved by
. (25)
Exploiting the leading-order term in η'max (Equation (12)), one finds an approximation for θS, the value of θ, at which the slope η' from Equation (25) becomes equal to the constant (maximal) slope that prevails throughout most of the range 0 ≤ θ ≤ (π/2). Equation (25) yields θS = O(μ2):
. (26)
Thus, for a valid expansion in powers of μ2 near the turning point, one should rescale θ as:
, (27)
so that z is of O(1). One now expands η(θ) in powers of μ2:
. (28)
Using Equations (28) and (7), one expands Equation (11) in powers of μ2, and solves for the higher-order corrections in Equation (28). The boundary conditions for all F2n are:
. (29)
Through O(μ6), F2n are found to be given by:
, (30)
, (31)
(32)
Detailed inspection of F2n for 2 n = 2, 4, 6, reveals that they all have the same qualitative dependence on θ. Except very near the turning point (θ = 0) they all tend to a linear dependence on θ. They vanish at roughly the same rate as θ varies from O(μ2) down to zero.
The surprising result is that an expansion that was meant to be valid around the turning point (θ ≈ O(μ2)) yields an excellent approximation to the solution throughout 0 ≤ θ ≤ 2 π for small μ. For sufficiently small μ, an excellent approximation is obtained with just the O(μ2) term:
. (33)
Figure 10 and Figure 11 show a comparison between the approximation, Equation (33), and the numerical solution of Equation (11) for μ = 0.1 (β = 0.99995). θS is then ≈ 0.0078. To stress the quality of the approximation, Figure 12 shows
Figure 10. Comparison between the O(μ2) approximation to the velocity, η'(θ) (using Equation (33), dashed) and the numerical solution of Equation (11). μ = 0.1. Dotted line—constant limit of velocity beyond turning point.
Figure 11. Comparison between O(μ2) approximation to η(θ) (Equation (33), dashed) and numerical solution of Equation (11). μ = 0.1.
Figure 12. Comparison between the O(μ2) approximation to the position, η(θ) (Equation (33), dashed) and the numerical solution of Equation (11) near the turning point. μ = 0.1.
how small the deviation of the approximate solution from the numerical solution remains even for
. Figure 13 shows the relative error in η'. It remains small for all 0 ≤ θ ≤ 2 π.
Finally, for θ =O(1), namely, away from the turning point, F4[z = θ/μ2] adds a μ2 contribution. An excellent approximation is obtained by:
. (34)
6. Fourier Expansion
The evolution of the velocity profile, η'(θ), in the extreme relativistic limit (β 1, μ 0) into a periodic step function (see Figure 4) offers unique characteristics for the Fourier spectrum of η'(θ). The Fourier coefficients tend to the limiting spectrum of the step function, the latter oscillating between ±(2/π) with a cycle length of (2π). The Fourier series for this step function is:
. (35)
The purpose of this Section is to discuss two aspects. The first is a demonstration of the qualitative difference between the Fourier spectra in the weak- and extremely-relativistic limits. The second is an additional test of the quality of the approximate solution obtained in Section 5.
6.1. Fourier Spectra in Weak- and Extreme-Relativistic Limits
Figure 14 presents a comparison between the β dependence of the Fourier coefficients of η'(θ), obtained from the numerical solution of Equation (11), and the spectrum of Equation (35). There is a qualitative difference between the β—dependence of the Fourier coefficients for small β and the spectrum of the
Figure 13. Relative error between O(μ2) approximation to velocity, η'(θ) (using Equation (33)) and numerical solution of Equation (11). μ = 0.1.
Figure 14. Fourier coefficients of η'(θ), from numerical solution of Equation (11), and spectrum of step-function Equation (35).
asymptotic step-function profile. For small β, the spectrum is dominated by a small number of modes. For β 0, the spectrum tends to a single mode, sin(θ). The analysis in Section 4 exhibits the β—dependence of the modes for small β: sin((2 n +1)θ) modes appear in O(β2n). As β approaches 1, the Fourier coefficients tend to the limiting spectrum, which is O(1) in μ.
6.2. Test of Quality of Approximate Solution in Extreme-Relativistic Limit
In Figure 15, a comparison between the Fourier coefficients of η'(θ), computed
Figure 15. Fourier coefficients of η'(θ) for approximate solution, Equations (28) and (32) and from numerical solution of Equation (11). μ = 0.2.
for the approximate solution, Equations (28) and (32), with those obtained from the numerical solution of Equation (11) for μ = 0.2 (β = 0.992) is presented for 1 ≤ 2 n +1 ≤ 41. The relative error between the exact numerical and the approximate values is below 5% for 1 ≤ 2 n + 1 ≤ 13. For μ = 0.1, the error is smaller than 0.5% for all 1 ≤ 2 n +1 ≤ 41.
7. Concluding Comments
The dynamics of the Relativistic Harmonic Oscillator (RHO) have been analyzed in the weak- and extreme-relativistic limits. Exploiting the fact that the period of oscillations can be computed in closed form and to any desired level of numerical accuracy, and using the scaled variables of Equation (10), for β < 1, a Normal Forms expansion yields a high quality approximation to the solution that is by far better than the standard application of asymptotic expansion methods.
The same formulation allows for the analysis of the extreme-relativistic limit as a boundary-layer problem, leading to a high-quality approximation for the solution for all 0 ≤ θ ≤ (2π) for (1 − β) << 1.
Statements and Declarations
The author declares that no funds, grants, or other support were received during the preparation of this manuscript.
The author has no relevant financial or non-financial interests to disclose.
1) The author made all contributions to the conception or design of the work. No data acquisition was required. No new software was created;
2) The author drafted the work or revised it critically for important intellectual content;
3) The author approved the version to be published; and
4) The author agrees to be accountable for all aspects of the work in ensuring that questions related to the accuracy or integrity of any part of the work are appropriately investigated and resolved.
Data Availability Statement
Data sharing is not applicable to this article, as no datasets were generated or analyzed during the current study.