next up previous contents
Next: The Weyl-Titchmarsh m-function Up: Resonance in the Dirac Previous: Numerical integration with graphical

Zero mass and its subsequent enhancement

Taking up the options offered by Equations 3 and 4 occasioned by seting the particle's mass to zero, the solution matrix for the massless Dirac harmonic oscillator would use the angle

$\displaystyle \varphi(x,0)$ = $\displaystyle \int_{x_0}^x\left(\frac{1}{2}t^2-E\right)dt$  
  = $\displaystyle \left[\frac{1}{6}t^3-Et\right]_0^x$  
  = $\displaystyle \frac{1}{6}x^3-Ex$ (9)

to give the solution
$\displaystyle {\bf Z}(x)$ = $\displaystyle \left( \begin{array}{cc}
\cos(\frac{1}{6}x^3-Ex) & - \sin(\frac{1...
...in(\frac{1}{6}x^3-Ex) & \cos(\frac{1}{6}x^3-Ex)
\end{array} \right) {\bf Z}(0),$ (10)

which is a pure rotation in the phase plane of rapidly increasing angular velocity. As a rotation, it will always describe solutions of constant amplitude.

When it comes time to include the mass in the solution, we need to work out the supplementary coefficient for ${\bf V}$ which is $m_0{\bf j}$ transformed by ${\bf U}$. Here ${\bf j}$ is the symmetric off-diagonal quaternion. Since it anticommutes with the antisymmetric off-diagonal ${\bf i}$ and given that ${\bf U}= \exp(\varphi {\bf i})$ we get altogether

$\displaystyle {\bf j}\ e^{\varphi {\bf i}}$ = $\displaystyle e^{-\varphi {\bf i}}\ {\bf j},$ (11)

and accordingly
$\displaystyle {\bf U}^{-1}m_0{\bf j}{\bf U}$ = $\displaystyle m_0 \left( \begin{array}{cc}
-\sin(\frac{1}{3}x^3-2Ex) & \cos(\fr...
...(\frac{1}{3}x^3-2Ex) & \sin(\frac{1}{3}x^3-2Ex)
\end{array} \right) {\bf Z}(0),$ (12)

which is a rotating mirror. When it rotates rapidly it will confuse the solution point, leaving it nearly constant. But when rotating slowly, as it does near the classical turning points, it can provoke considerable change in a solution matrix, as it also does when it is counterrotating.

When the matrices ${\bf U}$ and ${\bf V}$ are combined to get the complete solution ${\bf Z}= {\bf U}{\bf V}$ it is helpful to observe that one interpretation of a matrix product is that the second factor forms linear combinations of the columns of the first factor. Doing so, the solutions inherent in the matrix ${\bf V}$ are just rewritten in the basis defined by ${\bf U}$, but otherwise interpreted as any other soluton would be.


  
Figure 2: The coefficient matrix in the one dimensional Dirac equation can be split into a mass term and an energy term. Solving the energy first allows the mass to be included later on.
\begin{figure}
\centering
\begin{picture}
(300,380)
\put(0,-10){\epsfxsize=300pt \epsffile{factor2.0.eps}}
\end{picture}
\end{figure}

Figure 2 shows how little the mass eventually matters in comparison to the energy for a consistently increasing potential such as the parabolic potential well of the harmonic oscillator. On the other hand, near the classical turning points and through the forbidden zone, there are appreciable discrepancies which are more clearly evident in the diagrams of the factored potential than they were in their original form.

In principle, there is no reason why the factorization implementing the law of exponents for matrices should not be performed in the opposite order, regarding the mass term as fndamental. It would simply generate hyperbolic functions, with respect to which the energy terms would drive a rotating mirror which could be expected to raise and lower amplitudes passing from resonances to antiresonances.

After all, the final result is independent of the order of factorization.


  
Figure 3: The auxiliary factor shown in the phase plane for a series of energies between a resonance and an antiresonance. The tight centers of the spirals represent asymptotic values at infinity, where the mass has little influence on the potential. Due to the symmetry of the potential, the phase diagram at resonance and antiresonance is also symmetric.
\begin{figure}
\centering
\begin{picture}
(280,360)
\put(0,-10){\epsfxsize=280pt \epsffile{peak1.eps}}
\end{picture}
\end{figure}

In Figure 3 a series of six phase planes is shown, straddling the interval of energies between one resonance and the next. The point to be observed is the changing relationship between the amplitude of the interior wave function and the amplitude of the exterior wave function. The phase planes are shown separately to make the progression of ratios clearly evident. But then, for purposes of comparison, all six planes have been superimposed in Figure 4, where the magnitude of the variation is easier to discern.


  
Figure 4: For ready comparison, all the trajectories comprising Figure 3 have been superposed, to get a design which would do many a peacock proud.
\begin{figure}
\centering
\begin{picture}
(260,260)
\put(0,-5){\epsfxsize=260pt \epsffile{peak3.eps}}
\end{picture}
\end{figure}


next up previous contents
Next: The Weyl-Titchmarsh m-function Up: Resonance in the Dirac Previous: Numerical integration with graphical
Microcomputadoras
2001-01-09