On Apr 17, 1:12 pm, black head <
larryhar...@softhome.net> wrote:
> My question is:
>
> Can a center of momentum frame always be found for just the mechanical
> momentum of an electromagnetic system?
You means a rest frame? Take the 4x4 stress/energy/momentum tensor and
do a Lorentz transform that renders the {10,20,30,01,02,03} components
0. It should work if those components are small enough relative to the
{00} component.
If the field dynamics is described by a Lagrangian function L(B, E) of
the B and E fields, then the D and H fields can be defined as the
partial derivatives D = dL/dE, H = -dL/dB. The usual expression used
for L is L = 1/2 epsilon_0 (E^2 - B^2 c^2).
For a Lorentz-invariant and gauge-invariant Lagrangian, the most
general form is a function
L = L(I, J)
where
I = (E^2 - B^2 c^2), J = B.E.
Then, defining derivatives epsilon = dL/dI, theta = dL/dJ (the
coefficients, themselves, being functions of I and J), you get the
relations D = epsilon E + theta B, H = B/mu - theta E, where mu = 1/
(epsilon c^2).
Regardless of what L is, the stress/energy/momentum tensor have
components:
T^0_0 = D.E - L
(T^1_0, T^2_0, T^3_0) = (E x H)^i
(T^0_1, T^0_2, T^0_3) = (B x D)_j
(T^i_j: i,j = 1, 2, 3) = (D^i E_j + B^i H_j - delta^i_j (B.H + L))
This has a decomposition into the form
L = epsilon I + theta J - Delta/4, for a suitably defined function
"Delta"
T^0_0 = n U + Delta/4
T^0_j = -n P K
T^i_j = n m c^2 (XX - YY - U/(mc^2) ZZ)^i_j + delta^i_j Delta/4
where (X, Y, Z) for an orthonormal triad, with dyad notation used
(e.g. the ^i_j component of XX is X^i X_j); and with functions n, U
and m given by solving
n epsilon_0 = epsilon
B x D = -n/c root(U^2 - (mc^2)^2) Z
E x H = nc root(U^2 - (mc^2)^2) Z
D.E + B.H = 2n U.
E + i c B = (X + i Y)/root(epsilon_0) * e^{i phi}
I + i c J = mc^2/epsilon_0 e^{2i phi}
for a suitably defined phase phi.
For the frame (X,Y), one can define it by the conditions:
E = (x X cos phi - y Y sin phi)/root(epsilon_0)
B = root(mu_0) (y Y cos phi - x X sin phi)
provided x and y are chosen appropriately (x^2 = U + mc^2 and y^2 = U
- mc^2 should work).
I think that covers the basis. Check the calculations.
If you take components relative to the frame (X,Y,Z) you should be
able to readily define the Lorentz transformation that yields 0 for B
x D and E x H.
I think the relevant velocity u is the "Lorentz 1/2" of the velocity v
that emerges from
D x B = E x H/c^2 = m v/root(1 - (v/c)^2).
So it would be a solution to
v = 2 u/root(1 - (u/c)^2).