Toda oscillator

This article may be too technical for most readers to understand. Please help improve it to make it understandable to non-experts, without removing the technical details. (July 2012) (Learn how and when to remove this message)

|

In physics, the Toda oscillator is a special kind of nonlinear oscillator. It represents a chain of particles with exponential potential interaction between neighbors.[1] These concepts are named after Morikazu Toda. The Toda oscillator is used as a simple model to understand the phenomenon of self-pulsation, which is a quasi-periodic pulsation of the output intensity of a solid-state laser in the transient regime.

## Definition

[edit]

The Toda oscillator is a dynamical system of any origin, which can be described with dependent coordinate  x {\displaystyle ~x~} and independent coordinate  z {\displaystyle ~z~} , characterized in that the evolution along independent coordinate  z {\displaystyle ~z~} can be approximated with equation

d 2 x d z 2 + D ( x ) d x d z + Φ ′ ( x ) = 0 , {\displaystyle {\frac {{\rm {d^{2}}}x}{{\rm {d}}z^{2}}}+D(x){\frac {{\rm {d}}x}{{\rm {d}}z}}+\Phi '(x)=0,}

where  D ( x ) = u e x + v {\displaystyle ~D(x)=ue^{x}+v~} ,  Φ ( x ) = e x − x − 1 {\displaystyle ~\Phi (x)=e^{x}-x-1~} and prime denotes the derivative.

## Physical meaning

[edit]

The independent coordinate  z {\displaystyle ~z~} has sense of time. Indeed, it may be proportional to time  t {\displaystyle ~t~} with some relation like  z = t / t 0 {\displaystyle ~z=t/t_{0}~} , where  t 0 {\displaystyle ~t_{0}~} is constant.

The derivative x ˙ = d x d z {\displaystyle ~{\dot {x}}={\frac {{\rm {d}}x}{{\rm {d}}z}}} may have sense of velocity of particle with coordinate  x {\displaystyle ~x~} ; then  x ¨ = d 2 x d z 2 {\displaystyle ~{\ddot {x}}={\frac {{\rm {d}}^{2}x}{{\rm {d}}z^{2}}}~} can be interpreted as acceleration; and the mass of such a particle is equal to unity.

The dissipative function  D {\displaystyle ~D~} may have sense of coefficient of the speed-proportional friction.

Usually, both parameters  u {\displaystyle ~u~} and  v {\displaystyle ~v~} are supposed to be positive; then this speed-proportional friction coefficient grows exponentially at large positive values of coordinate  x {\displaystyle ~x~} .

The potential  Φ ( x ) = e x − x − 1 {\displaystyle ~\Phi (x)=e^{x}-x-1~} is a fixed function, which also shows exponential growth at large positive values of coordinate  x {\displaystyle ~x~} .

In the application in laser physics,  x {\displaystyle ~x~} may have a sense of logarithm of number of photons in the laser cavity, related to its steady-state value. Then, the output power of such a laser is proportional to  exp ⁡ ( x ) {\displaystyle ~\exp(x)~} and may show pulsation at oscillation of  x {\displaystyle ~x~} .

Both analogies, with a unity mass particle and logarithm of number of photons, are useful in the analysis of behavior of the Toda oscillator.

## Energy

[edit]

Rigorously, the oscillation is periodic only at  u = v = 0 {\displaystyle ~u=v=0~} . Indeed, in the realization of the Toda oscillator as a self-pulsing laser, these parameters may have values of order of  10 − 4 {\displaystyle ~10^{-4}~} ; during several pulses, the amplitude of pulsation does not change much. In this case, we can speak about the period of pulsation, since the function  x = x ( t ) {\displaystyle ~x=x(t)~} is almost periodic.

In the case  u = v = 0 {\displaystyle ~u=v=0~} , the energy of the oscillator  E = 1 2 ( d x d z ) 2 + Φ ( x ) {\displaystyle ~E={\frac {1}{2}}\left({\frac {{\rm {d}}x}{{\rm {d}}z}}\right)^{2}+\Phi (x)~} does not depend on  z {\displaystyle ~z~} , and can be treated as a constant of motion. Then, during one period of pulsation, the relation between  x {\displaystyle ~x~} and  z {\displaystyle ~z~} can be expressed analytically: [2][3]

z = ± ∫ x min x max d a 2 E − Φ ( a ) {\displaystyle z=\pm \int _{x_{\min }}^{x_{\max }}\\!\\!{\frac {{\rm {d}}a}{{\sqrt {2}}{\sqrt {E-\Phi (a)}}}}}

where  x min {\displaystyle ~x_{\min }~} and  x max {\displaystyle ~x_{\max }~} are minimal and maximal values of  x {\displaystyle ~x~} ; this solution is written for the case when  x ˙ ( 0 ) = 0 {\displaystyle {\dot {x}}(0)=0} .

however, other solutions may be obtained using the principle of translational invariance.

The ratio  x max / x min = 2 γ {\displaystyle ~x_{\max }/x_{\min }=2\gamma ~} is a convenient parameter to characterize the amplitude of pulsation. Using this, we can express the median value  δ = x max − x min 1 {\displaystyle \delta ={\frac {x_{\max }-x_{\min }}{1}}} as  δ = ln ⁡ sin ⁡ ( γ ) γ {\displaystyle \delta =\ln {\frac {\sin(\gamma )}{\gamma }}} ; and the energy  E = E ( γ ) = γ tanh ⁡ ( γ ) + ln ⁡ sinh ⁡ γ γ − 1 {\displaystyle E=E(\gamma )={\frac {\gamma }{\tanh(\gamma )}}+\ln {\frac {\sinh \gamma }{\gamma }}-1} is also an elementary function of  γ {\displaystyle ~\gamma ~} .

In application, the quantity  E {\displaystyle E} need not be the physical energy of the system; in these cases, this dimensionless quantity may be called quasienergy.

## Period of pulsation

[edit]

The period of pulsation is an increasing function of the amplitude  γ {\displaystyle ~\gamma ~} .

When  γ ≪ 1 {\displaystyle ~\gamma \ll 1~} , the period  T ( γ ) = 2 π ( 1 + γ 2 24 + O ( γ 4 ) ) {\displaystyle ~T(\gamma )=2\pi \left(1+{\frac {\gamma ^{2}}{24}}+O(\gamma ^{4})\right)~}

When  γ ≫ 1 {\displaystyle ~\gamma \gg 1~} , the period  T ( γ ) = 4 γ 1 / 2 ( 1 + O ( 1 / γ ) ) {\displaystyle ~T(\gamma )=4\gamma ^{1/2}\left(1+O(1/\gamma )\right)~}

In the whole range  γ > 0 {\displaystyle ~\gamma >0~} , the period  T ( γ ) {\displaystyle ~{T(\gamma )}~} and frequency  k ( γ ) = 2 π T ( γ ) {\displaystyle ~k(\gamma )={\frac {2\pi }{T(\gamma )}}~} can be approximated by

k fit ( γ ) = 2 π T fit ( γ ) = {\displaystyle k_{\text{fit}}(\gamma )={\frac {2\pi }{T_{\text{fit}}(\gamma )}}=}
( 10630 + 674 γ + 695.2419 γ 2 + 191.4489 γ 3 + 16.86221 γ 4 + 4.082607 γ 5 + γ 6 10630 + 674 γ + 2467 γ 2 + 303.2428 γ 3 + 164.6842 γ 4 + 36.6434 γ 5 + 3.9596 γ 6 + 0.8983 γ 7 + 16 π 4 γ 8 ) 1 / 4 {\displaystyle \left({\frac {10630+674\gamma +695.2419\gamma ^{2}+191.4489\gamma ^{3}+16.86221\gamma ^{4}+4.082607\gamma ^{5}+\gamma ^{6}}{10630+674\gamma +2467\gamma ^{2}+303.2428\gamma ^{3}+164.6842\gamma ^{4}+36.6434\gamma ^{5}+3.9596\gamma ^{6}+0.8983\gamma ^{7}+{\frac {16}{\pi ^{4}}}\gamma ^{8}}}\right)^{1/4}}

to at least 8 significant figures. The relative error of this approximation does not exceed  22 × 10 − 9 {\displaystyle 22\times 10^{-9}} .

## Decay of pulsation

[edit]

At small (but still positive) values of  u {\displaystyle ~u~} and  v {\displaystyle ~v~} , the pulsation decays slowly, and this decay can be described analytically. In the first approximation, the parameters  u {\displaystyle ~u~} and  v {\displaystyle ~v~} give additive contributions to the decay; the decay rate, as well as the amplitude and phase of the nonlinear oscillation, can be approximated with elementary functions in a manner similar to the period above. In describing the behavior of the idealized Toda oscillator, the error of such approximations is smaller than the differences between the ideal and its experimental realization as a self-pulsing laser at the optical bench. However, a self-pulsing laser shows qualitatively very similar behavior.[3]

## Continuous limit

[edit]

The Toda chain equations of motion, in the continuous limit in which the distance between neighbors goes to zero, become the Korteweg–de Vries equation (KdV) equation.[1] Here the index labeling the particle in the chain becomes the new spatial coordinate.

In contrast, the Toda field theory is achieved by introducing a new spatial coordinate which is independent of the chain index label. This is done in a relativistically invariant way, so that time and space are treated on equal grounds.[4] This means that the Toda field theory is not a continuous limit of the Toda chain.

## References

[edit]

1. ^ a b Toda, M. (1975). "Studies of a non-linear lattice". Physics Reports. 18 (1): 1. Bibcode:1975PhR....18....1T. doi:10.1016/0370-1573(75)90018-6.
2. ^ Oppo, G.L.; Politi, A. (1985). "Toda potential in laser equations". Zeitschrift für Physik B. 59 (1): 111–115. Bibcode:1985ZPhyB..59..111O. doi:10.1007/BF01325388. S2CID 119657810.
3. ^ a b Kouznetsov, D.; Bisson, J.-F.; Li, J.; Ueda, K. (2007). "Self-pulsing laser as Toda oscillator: Approximation through elementary functions". Journal of Physics A. 40 (9): 1–18. Bibcode:2007JPhA...40.2107K. CiteSeerX 10.1.1.535.5379. doi:10.1088/1751-8113/40/9/016. S2CID 53330023.
4. ^ Kashaev, R.-M.; Reshetikhin, N. (1997). "Affine Toda field theory as a 3-dimensional integrable system". Communications in Mathematical Physics. 188 (2): 251–266. arXiv:hep-th/9507065. Bibcode:1997CMaPh.188..251K. doi:10.1007/s002200050164. S2CID 17196702.