Register to reply 
Maxwell equations in Lorenz gauge 
Share this thread: 
#1
Apr1312, 08:06 PM

P: 148

In the Lorenz gauge, the Maxwell equations reduce to four inhomogenous wave equations, with the charge density acting as the source for V, and the current density for A.
For now, just take a static charge distribution  say, a point charge at the origin. It is well known that a static charge distribution leads to the electrostatic field; for our point charge at the origin, V(r, t=infinity) = 1/r. However, if the charge is static, V would be constantly increasing...wouldn't it? So at t=infinity, V would be infinity everywhere? Now, running a simulation (because I suck at PDEs), it appears that while this appears to be true, E = gradV does appear to stay constant...so perhaps this is a red herring, but I would feel much more confident if I could find an analytical solution (or a logical explanation). Any takers? 


#2
Apr1412, 01:53 AM

Sci Advisor
HW Helper
PF Gold
P: 2,603

If the charge distribution is static, then [itex]\dot{\rho}(\mathbf{x})=0[/itex]. The scalar potential is
[tex] V = \int \frac{\rho(\mathbf{x}')}{\mathbf{x}\mathbf{x}'} d\mathbf{x}',[/tex] so [itex] \dot{V}=0[/itex] as well. So [itex]V[/itex] is static. Note that when the potential is static, the inhomogenous wave equation reduces to the Poisson equation of electrostatics. For the point charge at the origin, [itex]V\sim 1/r[/itex] at all times. The only divergence is at the origin. 


#3
Apr1412, 08:10 AM

P: 148

Right  but that assumes that it was always static, and therefore you're assuming what I want to derive: that's the green function for the Laplacian, not the Wave equation.
Put another way, you're assuming it was static for all time; But if the point charge was turned on any time after t=infinity, the static approximation won't hold (except possibly at t=infinity, but it seems that at any large but finite time, V will be much larger than 1/r... So to restate the problem (all scalar factors set to 1) V,tt = V,xx + V,yy + V,zz + p p(t<0) = 0 p(t>=0) = point charge at origin I'm wondering if one trick might be putting a boundary condition so V(r=infinity) = 0; since it only falls off as 1/r that would (I think?) be an extra assumption with no apparent justification (other than the fact it might give the right answer ;) ). This sort of reminds me of how I understand the process of renormalization in QFT, which of course is exactly that sort of unjustifiable hack, but again, I could be way off base here. Thanks, Justin 


#4
Apr1412, 08:42 AM

Sci Advisor
Thanks
P: 2,337

Maxwell equations in Lorenz gauge
Take the retarded Green's function of the wave equation and use it to calculate the electromagnetic field for a static charge distribution. On the one hand that's a weird way to calculate the Coulomb potential, on the other it furthers the understanding of the matter quite well :).



#5
Apr1412, 08:50 AM

Sci Advisor
HW Helper
PF Gold
P: 2,603

[tex] V(\mathbf{x},t) = \int \frac{\delta(t't + \mathbf{x}\mathbf{x}'/c)}{\mathbf{x}\mathbf{x}'} \rho(\mathbf{x}',t') d\mathbf{x}' dt' .[/tex] This is derived in, for example, Ch 6 of Jackson. [tex]V(\mathbf{x},t) = \begin{cases} 0, & t<\mathbf{x}/c, \\ \frac{\rho_0}{\mathbf{x}}, & t\geq \mathbf{x}/c. \end{cases} [/tex] At a given distance from the origin, the potential is zero until late enough times that a signal traveling at light speed can reach the point. After that, we have the ordinary static potential. This is as expected. By causality, we can't measure the potential until the signal has reached us. Once it does, the charge is effectively static, so we should measure the static field. 


#6
Apr1412, 10:30 AM

Sci Advisor
Thanks
P: 2,337

This cannot be since then you have an electric field popping out of nothing instantaneously and instantly at any point in space. Of course this comes from the totally unphysical assumption of "jjustinn" concerning his charge distribution: You cannot pop up a point charge out of nothing at [itex]t=0[/itex]. This violates not only causality but also the necessary condition of charge conservation and thus gauge invariance. What you can have, is a charge sitting at the origin forever (which is of course also an unrealistic assumption, but it's not in contradiction with any physical laws). Then you get immediately the Coulomb field as it shoud be (here written in its form using nonrationalized Gaussian units).



#7
Apr1412, 10:50 AM

P: 148

I'm familiar with the derivation of the Wave equation green function, and that's another thing that's puzzled me  the Green function seems to say, as you point out, that the potential at point r goes from 0 to 1/r, then stays at 1/r for all time  which not only sounds reasonable (since the charge appearing at t=0 is jut as discontinuous), but agrees with observations. So perhaps my problems are stemming from my latticebased simulation, where the charge is acting as a constant source, adding to the V around it, and since there is no dissipation / damping, it never goes away. Or, it might also have something to do with the fact I'm only updating it in the plane, and the 2D Green function isn't a spreading delta function, but IIRC a log function of some sort. To give some background, I originally thought my simulation was broken, since V was monotonically increasing everywhere (while E stabilized to a constant value)  but I started reading this (http://psych.colorado.edu/~oreilly/p...eilly05_md.pdf), and he mentions adding a empiricallydetermined damping factor at the (artificial, Sommerfieldconditioned) boundary "to achieve the 1/r dependence"; before I'd seen this, I'd found that adding a damping term at/near the boundary had a similar effect, but rejected it outright as nonphysical. Anyway, thanks for taking the time to respond. I'll play with this some more and will undoubtedly be back with more dumb questions with crucial unstated assumptions ;) 


#8
Apr1412, 12:55 PM

P: 148

Ah hah! I think this explains my problem: (from http://ieeexplore.ieee.org/xpl/login...er%3D5492190):
So...this actually brings up anther question, for a new thread: shouldn't our point charge at origin generate the same waves in the XY plane (or any plane through the origin) as a 2D charge would, via radial symmetry? 


#9
Apr1412, 01:05 PM

P: 148

However, I did find something interesting when trying to make it more physical: by adding an equal/opposite negative charge, its outward monotonicallydecreasing negative wave cancels the monotonicallyincreasing wave from the positive source, making the average value stable, rather than increasing to infinity...but it's not too useful to have a simulation where the net charge has to be zero. This is what initially made me think of "renormalization" (used in the lay 'removing infinities' sense)  it wasn't until later I found that gradV actually stayed constant. In any case, would you care to expand on that last comment? If you are uncomfortable with a single charge coming out of nowhere, how about a positron/electron pair pulled out of the vacuum? Or since this is explicitly classical, how about if I separate the charges on a glass rod and a bit of fur at t=0, then fire the glass rod from a cannon into the distance at t= 1? 


Register to reply 
Related Discussions  
Invariance of maxwell's equations under Gauge transformation  Advanced Physics Homework  6  
Use of the Lorenz Gauge  Advanced Physics Homework  6  
Electromag  Lorenz gauge  Advanced Physics Homework  3  
Lorenz Equations  Differential Equations  0 