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.

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

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

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

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

The derivative ${\displaystyle ~{\dot {x}}={\frac {{\rm {d}}x}{{\rm {d}}z}}}$ may have sense of velocity of particle with coordinate ${\displaystyle ~x~}$; then ${\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 ${\displaystyle ~D~}$ may have sense of coefficient of the speed-proportional friction.

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

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

In the application in laser physics, ${\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 ${\displaystyle ~\exp(x)~}$ and may show pulsation at oscillation of ${\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.

Rigorously, the oscillation is periodic only at ${\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 ${\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 ${\displaystyle ~x=x(t)~}$ is almost periodic.

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

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

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

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

The ratio ${\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 ${\displaystyle \delta ={\frac {x_{\max }-x_{\min }}{1}}}$ as ${\displaystyle \delta =\ln {\frac {\sin(\gamma )}{\gamma }}}$; and the energy ${\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 ${\displaystyle E}$ need not be the physical energy of the system; in these cases, this dimensionless quantity may be called quasienergy.

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

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

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

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

${\displaystyle k_{\text{fit}}(\gamma )={\frac {2\pi }{T_{\text{fit}}(\gamma )}}=}$

${\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 ${\displaystyle 22\times 10^{-9}}$.

At small (but still positive) values of ${\displaystyle ~u~}$ and ${\displaystyle ~v~}$, the pulsation decays slowly, and this decay can be described analytically. In the first approximation, the parameters ${\displaystyle ~u~}$ and ${\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]

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.