Sunday, May 9, 2010

An interesting feature of the quasiclassical atom

As everyone who's taken a college level physics class should know, the problem of the classical atom in the days of Rutherford, with (for hydrogen, say) the electron orbiting the proton like the moon orbits the earth, is radiative decay. This system has a nonzero, nonstationary, electric dipole moment, and the radiation intensity due to this is proportional to the square of the magnitude of the electric dipole acceleration. This causes energy loss via radiation reaction, and the classical atom decays to a point in a small fraction of a second.

Something I've come to realize just recently however is that if the electron has an intrinsic magnetic moment (which of course it does, unbeknownst to Rutherford, or Bohr in 1913), this radiative decay is no longer so obviously necessarily the case. This is because a moving magnetic dipole necessarily acquires an electric dipole moment, as the cross product of the magnetic moment vector with its velocity, and it's easy to arrange things so that this opposes the electric dipole moment part that's due to the separation of charge. For example, if the electron is in a circular orbit and its intrinsic magnetic moment is oriented perpendicular to the orbital plane, the electric dipole moments from the two sources (separation of charge and moving magnetic dipole) are either parallel or antiparallel, depending on the orbital direction.

So then, for a perfectly circular orbit, since the magnitude of the motion-acquired electric dipole is proportional to the velocity, and the velocity increases as the orbital radius shrinks, whilst the dipole moment due to the separation of charge shrinks, there is a unique radius where the total electric dipole moment vanishes. If the quasiclassical atom were dynamically stable otherwise at this radius, it would not need to decay radiatively.

When this realization came to me the other day, the next thing I wanted to do was determine what this radius would be and whether it would turn out to be the Bohr radius. After all, it will certainly involve Planck's constant via the intrinsic magnetic moment being the Bohr magneton (which is e * hbar/ (m_e * c) in Gaussian units).

Turns out, the radius where the electric dipole moment vanishes is not the Bohr radius. The formula for it is easy to get, but it isn't much like the Bohr radius formula, although it also involves the charge and mass and h-bar but in the wrong powers and also has the speed of light which is absent in the Bohr radius formula. In magnitude however the radius is about a couple of orders of magnitude smaller than the Bohr radius. So it's between the atomic scale and the nuclear scale. Isn't it interesting that the scale where the two dipole moments can cancel is near the atomic scale? I think so. It's worth remembering here that in the classical atom without spin there is nothing obvious to determine scale at all.

I suppose, there's no reason for this particular radius to be dynamically preferrred, in the sense of it being a local energy minima. The dipole moment can momentarily vanish yet its acceleration not. So it has no direct particular significance. I think it's interesting though that acquired electric dipole moment is about the right magnitude so that if it jiggled around at the orbital frequency and just so, it might cancel the charge-separation dipole acceleration. After all, it's not the total electric dipole moment that has to vanish but rather its acceleration.

This is something I'm going to pursue further. I suspect pretty strongly that nobody has ever investigated it. When people originally tried to make the classical atom work, they had no knowledge of the electron spin. Then, within a year or so of the spin becoming known, both the Schroedinger and Heisenberg formulations of quantum theory were established and completely supplanted the classical atom idea (for good reason). A couple of years later the existence of the spin and the spin orbit coupling were even fully described by Dirac's relativistic quantum theory. Physicists have never looked back, but fortunately they have me willing to go back and make sure they didn't miss something important.

Here is how to get an equation of nonradiativity of electric dipole radiation for a classical atom where the particles are allowed to have intrinsic spin and associated intrinsic magnetic moment.

Let d represent the total electric dipole moment of the atom. Since the electric dipole radiation intensity is proportional to |d^2d/dt^2|^2, and assuming d is an ordinary vector, the nonradiativity requirement can be expressed simply as d^2d/dt^2 = 0.

For simplicity, I will treat hydrogen, but the approach generalizes to more complex atoms easily enough.

The total electric dipole moment of the atom consists of a part due to separation of the point charges, and, since the electron is also intrinsically a magnetic dipole, a part v x m/c due to the translational motion of the electron, where "x" indicates the vector cross product. So the total electric dipole moment is, for hydrogen

d = -e r + v x m/c

where e is the fundamental charge, r is the displacement of the electron from the proton, v = dr/dt is the electron velocity, c is the speed of light (approximating that the proton is stationary), and m is the intrinsic magnetic moment of the electron (which is known to be the Bohr magneton). The condition d^2d/dt^2 = 0 is then easily found to be equivalent to

e c a = (da/dt) x m + v x (d^2m/dt^2) + 2a x (dm/dt)

where a = dv/dt = d^2r/dt^2

It's also possible to write an equation for magnetic dipole nonradiativity, and that one obtains something like quantization of the orbital angular momentum, if the nuclear magnetic moment is not neglected. However the composite nature of a nucleus even consisting of a single proton cannot be disregarded here, so I was forced to consider Positronium. In Positronium I get that the magnetic dipole radiation vanishes if the orbital angular momentum is (3/4)*h-bar, if it is assumed the electron spin component parallel to the orbital angular momentum is h-bar/2. It would be more interesting if the factor of 3/4 were not present, of course, but another problem is that the magnetic dipole radiation intensity is negligible compared to the electric dipole radiation intensity, in the atomic scale (although not at the nuclear scale, where the magnetic dipole radiation intensity becomes larger than the electric). So naturally it's more interesting if the electric dipole radiation could be shown to vanish for a nonzero atomic radius. I simply never until recently recognized it could even plausibly vanish, because I was not realizing the total electric dipole moment has the additional contribution due to the intrinsic magnetic moment. I knew about the acquired electric dipole moment of a moving magnetic dipole, but I failed to realize it must contribute to the radiation intensity in this way.

So now I am trying to determine if there are any interesting consequences of this equation, that might also correspond to known physics.

Monday, May 3, 2010

Two faces of one coin?

As a model of atomic spin-orbit interaction in atomic hydrogen, consider two continuous-loop superconducting wire coil magnets, one larger than the other, in free-fall in outer space. Once set up, a current will circulate indefinitely without any applied voltage in a superconducting wire loop, creating a permanent magnet. Think of the small coil's magnetic dipole moment as representing the electron intrinsic magnetic moment due to its spin, and the large coil's as representing the electron orbital magnetic moment. The amount of current and number of windings here are arbitrary and not necessarily the same in both coils.

Now suppose the smaller coil is located within the empty space in the center of the larger, and that the solenoidal axes of the two coils are not perfectly aligned. Then, basic electromagnetics tells us the two coils will mutually precess. That is, the axis of each coil instantaneously precesses around the axis of the other. For the smaller coil inside the larger, this is very basic electromagnetic theory, since (assuming the inner coil is sufficiently small) the magnetic field due to the larger coil is for practical purposes here a constant value aligned with the large coil axis. The smaller coil, being an ideal current-loop magnetic dipole, feels a torque in the large coil's magnetic field, and according to basic mechanics the instantaneous change in the small coil angular momentum must equal the torque exactly. (The angular momentum of the coil is that of the circulating electrons forming the current.) However, since the torque at all times is perpendicular to the angular momentum of the small coil, it can't change the angular momentum magnitude, only its direction, and so the only allowed motion is the precession of the small coil around the axis of the larger.

For similar but slightly more complicated reasons, the larger coil also precesses around the axis of the smaller. The added complication to precession of the outer coil is that the magnetic field due to the smaller coil is not everywhere the same direction relative to the motion of the the charge carriers in the larger coil. The net force must be found by averaging over small current elements of the larger coil, but the end result is similar in that the larger coil also precesses instantaneously around the instantaneous axis of the other.

That the two coils are instantaneously precessing around each others' axes may seem difficult to envision, but visualization can be made easier by invoking the principle of conservation of angular momentum. If angular momentum is conserved, and there are no other torques acting on the coils than those already noted, the total angular momentum of the two coils must be a constant. Since the total angular momentum here is simply the vector sum of the angular momenta of the two coils, the only possible motions are those where the two coils precess around the total angular momentum direction with equal angular velocity. The situation is illustrated on Figure 1. The vector angular momenta of the large and small coils are represented as arrows labeled L and s, and their vector sum is J = L + s.

Figure 1

We should also consider what happens to the total magnetic moment of the system consisting of the two coils in free space. If the total magnetic moment is not stationary, energy will be lost as radiation into free space, and the precession motion will eventually cease as the angular momentum vectors come into alignment due to the radiation reactive forces. However, classical electrodynamics also tells us that the magnetic dipole moment of each coil is proportional to the angular momentum of each coil. Further, the constant of proprtionality is the same for each coil provided that the charge-to-mass ratio of the current carriers is the same for each coil. For our superconducting coils the charge carriers are electrons in each case, so it is concluded that stationarity of the total angular momentum implies in this case stationarity of the total magnetic moment as well.

If, instead of simply assuming that the total angular momentum is a constant of the motion, we had gone through the machinations of actually calculating the magnetic moments and fields and torques and rates of changes of the angular momenta of each of the coils, we would have found that indeed that angular momentum is conserved according to electrodynamics. Thus we have arrived at a pair of completely unsurprising conclusions. The total angular momentum and total magnetic moment of the isolated system of the two coils, that is not being acted on by any external mechanical or electromagnetic forces, are both constants of the motion.

Next it's worthwhile to consider the case where the charge carriers in one coil have a different charge-to-mass ratio than those in the other, and determine whether angular momentum and magnetic moment are still conserved in spite of this. So, suppose that the charge carriers in one of the coils are muons instead of electrons. The muonic charge is the same as the electron's, but the muon mass is much larger than the electron mass. If we suppose the current is kept constant in spite of switching the carriers from electrons to muons, then the magnetic field generated is the same, but the angular momentum of the muonic coil is much larger than it was for electrons. The torque on the muonic coil due to the electronic coil is unchanged, however, so the angular velocity of the precession must decrease in the muonic coil. The net effect is that although the total angular momentum remains approximately constant, in this case the total magnetic moment is no longer stationary. The moving total magnetic dipole moment radiates energy, and so the back-force of radiation reaction will cause the two coils eventually to align with each other. Although the total mechanical angular momentum of the system is no longer a constant of the motion, angular momentum is nonetheless conserved because the electromagnetic radiation also carries angular momentum, and the total of the mechanical and field angular momenta is constant.

To this point we have been considering only classical electrodynamics. However, when it was first recognized that electrons must possess an "intrinsic", or non-orbital, angular momentum and an associated intrinsic magnetic dipole moment, it was also immediately recognized that the constant of proportionality, or gyromagnetic ratio, between intrinsic angular momentum and magnetic moment must be twice the classically-expected value. This doubling of the gyromagnetic ratio is conventionally expressed by including in the expression for the magnetic moment in terms of the intrinsic angular momentum (or "spin"), an additional "g-factor" of two. This factor of two has never been derived from classical electrodynamic electron models, if it is assumed that the charge-to-mass ratio of the electron is a constant throughout its volume.

Now let us consider what happens when one of the two coils is allowed to have a nonunity g-factor. For an electron-like g-factor of two, the magnetic moment is doubled compared to the classical g-factor of unity, for a fixed value of angular momentum. If we suppose the g-factor of the small coil is changed to two, while keeping the large coil classical g-factor of unity (and returning again to electron currents in both coils), then it's clear that the rate of precession of the small coil is doubled. The magnetic field at the small coil due to the large coil is unchanged, but the torque on the small coil has doubled due to its magnetic moment doubling, whilst the small coil angular momentum is unchanged. This is a net doubling of the rate of precession of the small coil. However, the rate of precession of the large coil also doubles, since the magnetic field of the small coil that drives its precession has doubled. So, the total angular momentum is a constant of the motion, if there are no radiative effects. However, evaluating the total magnetic moment, it's plain that it cannot be stationary if the total angular momentum is stationary (and both coils are precessing), since the magnetic moments are no longer both in equal proportion to the coils' angular momenta. Thus, seen from afar, the magnetic moment of the system is precessing, and this type of motion must radiate electromagnetic energy and momentum, and decay due to radiation reactive force.

Finally, there is an important additional effect that needs to be added to our system of two electronic-current-carrying coils, if it is to be a useful toy representation of atomic spin-orbit coupling, and that is the Thomas precession. The Thomas precession (also sometimes called Wigner rotation) is a relativistic effect wherein a reference frame that is accelerating and translating relative to another must also rotate compared to the other. The effect of the Thomas precession is to approximately halve the rate of precession of the intrinsic spin of a particle with a g-factor of two such as an electron. If the small coil with the g-factor of two also undergoes Thomas precession, while the large coil does not, plainly the angular momentum will no longer be a constant of the motion, even neglecting radiation reaction, if there is no change in the rate of the precession of the large coil due to the Thomas precession acting on the small coil. But, the nature of the Thomas precession is such that it only acts on the electron intrinsic spin and not on the orbital angular momentum vector, represented here by the angular momentum of the current-carrying electrons of the large coil. Also, the Thomas precession has only reduced the rate of precession of the small coil, not its magnetic field as experienced by the large coil. So the rate of precession of the large coil cannot change as would be necessary to maintain total angular momentum constancy, and thus it's clear that the total mechanical angular momentum is not a constant of the motion. But what about the field angular momentum? Can it restore angular momentum constancy? It turns out it cannot. Calculating the total magnetic moment, we find that the factor of one-half due to the Thomas precession (the celebrated Thomas factor that resolved the spin-orbit coupling anomaly of the 1920s) has compensated the g-factor of two that previously caused the total magnetic moment to move in spite of stationarity of total angular momentum. Present the Thomas precession, we now have the situation of a stationary magnetic moment in spite of moving total angular momentum, so that there is no need for the motion of the total angular momentum to lead to radiative decay, nor require an external torque.

The precession of the total angular momentum, even absent an externally-applied magnetic field, due to the spin-orbit coupling is a feature of quantum theory, but it has been seen here to be simply the classical electrodynamic consequence of the electron possessing a g-factor of two. This raises a further question. Will this same sort of behavior, of motion of the total angular momentum but stationarity of total magnetic moment, occur for other g-factor values than two? It turns out, g=2 is the unique g factor where the total magnetic moment can be stationary in the presence of spin-orbit interaction. (The proof can be found in my paper, "Regarding L. H. Thomas's paper of 1927 ....," linked at right.) For any other g-factor, including the classical value of unity, in the presence of Thomas precession, neither the total angular momentum nor the total magnetic moment is a constant of the motion. This suggests to me that the electron g-factor being two may be related to or even a consequence of the Thomas precession.