Processing math: 68%

Friday, October 31, 2014

An Interesting Force.

A plot of the force between electrons in parallel orbits given by,

Fa=neq24εorre{(2ae)2+(rre)2}3/2dr

is provided below.


It is both repulsive (positive) and attractive (negative) depending on the value of re.  A most interesting force.

Sunday, October 26, 2014

Ionization Redefined

If the force between an orbiting electron and the nucleus is repulsive what then is ionization energy, Ei?

Ei=mec2{ln(reoref)+1refrefreo(refreo)}q24εo1re+Gmamere+reneq24εorre{(2ae)2+(rre)2}3/2dr

The terms are,

mec2{ln(rporpf)+1rpfrpfrpo(rpfrpo)}  

work done against the centripetal force, while the circular velocity remains constant.


q24εo1re

work recovered from the repulsive force of the nucleus B orbit.

Gmamere

work done against gravitational attraction.

and

reneq24εorre{(2ae)2+(rre)2}3/2dr

=neq24εo14a2e=neq24εo12ae

work done against the attraction of neighboring electron B orbits.

The work against centripetal force explodes immediately.  The centripetal force is provided by drag at terminal velocity.  Work done against the centripetal force is essentially work against drag.  The electron gains velocity initially as a result of electrostatic attraction.  Upon achieving light speed, drag curls its path into a helix. The electron perform circular motion in the plane perpendicular to its direction of travel.

To resolve this issue, we see that this work is zero for incremental change in ref, ie. when

ref=reo+Δr    as  Δr0

this work is zero.  In other words, move slowly lo.

And we are left with,

Ei=q24εo1re+Gmamere+neq24εo12ae

Ei=Gmamere+q24εo{ne2ae1re}

which boldly is the new expression for ionization energy.  ne is the number of electrons in the stack minus one,  ae is the radius of an electron and  ma is the mass of the nucleus.

Stacked Orbits

Such stacked orbits means that the electrons are held more firmly in orbit and we see a decrease in reactivity across the periodic table as electrons are added to each consecutive elements.  If not for the Lorentz's force, each additional electron should be more readily removed from the atom and thus makes the atom more reactive.  The existence of B orbit and the resulting stacked orbits provide an explanation for the increase in stability across the periodic table.

The nucleus is not involved in this stability.


A Brand New Attractive Force

Consider two electrons in parallel orbits,


FLsin(θ)=Fa

The attractive force between the electrons is the Lorentz's force.  Fa is the force in the plane of the orbit that keeps the orbital radius constant due to the attraction of the neighboring orbit.

tan(θ)=Δr2ae

FL=neq24εo(2ae)2cos2(θ)=neq24εo(2ae)24a2e(2ae)2+(Δr)2=neq24εo1(2ae)2+(Δr)2

Therefore,

Fa=neq24εo(2ae)2(2ae)2(2ae)2+(Δr)2(Δr)(2ae)2+(Δr)2

Fa=neq24εo(Δr){(2ae)2+(Δr)2}3/2

Δr=rre

Fa=neq24εorre{(2ae)2+(rre)2}3/2

A brand new attractive force! And we have,

mec2re=Ac2+Gmamer2e+neq24εorre{(2ae)2+(rre)2}3/2Zq24εo(rerp)2TnTe4πτor2e

a new expression for the forces in action on an orbiting electron.  Under normal circumstances,

r=re    and    Fa=0

But this force will definitely play a part in ionization.

2ae is twice the radius of an electron, and ne the number of electrons in orbit minus one.  We are assuming that the orbiting electrons are stacked closest possible in parallel orbits.

Wanted: Attractive Force. Electrostatic Need Not Apply

There is no other way.  The electron is attracted to the nucleus, it is accelerated to light speed and begins curl in the direction perpendicular to its forward motion.  It transcribe a helical path.  It is the drag force that provides the centripetal force that rotates the electron in the direction perpendicular to its direction of travel.  The magnetic repulsive force in the vicinity of the nucleus and thermal gravity pushes the electron way, and the new centripetal force is given by,

mec2re=Ac2+Gmamer2eZq24εo(rerp)2TnTe4πτor2e

where Z is the atomic number, ma is the atomic mass, and rp is the proton spin radius.

The equation does not have the electrostatic term,

+Zq24πεor2e

as the charges are in constant motion.

The gravitational term,

GZ=Gmamer2e

plays a much more important role.

Ac2r2emec2re+(GmameZq24εor2e(rerp)2TnTe4πεo)=0 --- (*)

we see that the equation for re is still essentially a quadratic equation.  For heavy nucleus elements, it is expected that the gravitation term is high enough to give a positive y-intercept to equation (*).  A positive y-intercept will allow the parabola to intercept the x-axis twice and give two positive solutions to re.  The lower solution corresponds to the an orbit in the valence band and the higher solution is in the conduction band.  An electron can make the transition from the valence band to the conduction band, both valid re solutions to the force equation (*).  This behavior is expected of heavier elements but not all elements.

Increasing temperature can bring the y-intercept down, and eventually a possible zero valued solution to re.  This temperature would correspond to Tc, the temperature for the formation of plasma.

Decreasing temperature or increasing the nucleus mass, ie the GZ term, moves the curve upwards and moves the two roots closer until a double root is created.  At the double root, the kink point, the root and the minima all coincide.  The element is very conductive, very shinny as the electron readily crosses kink point and emit a packet of energy.

But it is commonly believed that GZ is small and does not play a part in this situation.  Unless there is another attractive force in the equation, I believe common opinion is wrong, gravity plays a much more important role and cannot be ignored.  Electrostatic attraction between the charges cannot account for a positive y-intercept even if it is admitted because it is already smaller than the magnetic repulsive force by a factor of π.

Let's look for another attractive force that hold the electron to its orbit of re.

Saturday, October 25, 2014

Why Half The Intensity, Albert?

If the electron and the nucleus are held apart by a repulsive force then the work done as a photon passes is,

Es=mpc2{ln(rporpf)+1rpfrpfrpo(rpfrpo)}=q24εo{1ref1reo}=ΔPEre

where  rpf and rpf are the initial and final radius of the photon in its helical path, reo and ref are the initial and final radius of the electron orbit around the nucleus.  The expression on the left was derived in the post "Light Dispersion?" from considering the work done along the centripetal force in moving the photon to a different orbital radius.  The expression on the right is the potential energy change in a field of magnetic repulsive force FL=iπqE, also the Lorentz's force between the electron and the nucleus.  The negative sign indicates that a gain in energy in one equals the lost of energy in the other.

Up till now, the only postulation as to why a photon and an orbiting electron repel each other is thermal gravity between two hot bodies. If we examine the likely magnetic interaction between an orbiting electron and a passing photon,


We find that when the electron orbits are also in the same anti-clockwise sense looking onward to the direction of the approaching photon, the particles should be attracted to each other based on the their interacting magnetic fields alone.  The electrons are pulled further from their orbits and the photons pulled down to a smaller radii.  The electron on reaching high orbit is not going to to be ionized.  The photon is still above the electron, the charge cannot escape the hold of the atom.  Given the 1re relationship for PEre, the energy change from moving a distance Δre is always higher downwards than upwards.  (When the photon has a smaller radius, the electrons are pulled down to smaller orbital radii. This can cause ionization.)

The orbits can of course be re-oriented to an clockwise sense, in which case the particles are like opposite current in parallel wires and be repelled from one another.  This was the scenario proposed previously where the particles repel each other resulting in diffraction, ionization, glare, etc. (In the case of a smaller photon radial path the electrons are pushed upwards to higher orbit. This will not cause ionization.)

Note: From the above expression by setting ref for ionization and rpo<rre, it is not possible to prove conclusively for all cases that the atom cannot be ionized by pushing the electron to higher orbit.

Since, both orientations are equally likely and all orientations in between, any photon at and beyond the threshold frequency (at the appropriate interaction radius r) on approach to an atom has only half the chance of ejecting a electron.  As a whole, only half the photons would have a photoelectric effect.  This conclusion is the same as we consider the relative size of the orbital radius and helical path radius.

If the interaction is purely electrostatic, we would expect the photon to eject the electron in both orientations and that all photons at the threshold frequency would cause a photoelectric effect.

We have assumed that the electron orbits re-orientate themselves slowly, not at the light speed that the photon is approaching.  We have also ignored the B field due to the spinning positive charge which is always perpendicular to the B field due to the negative charge as they both travel down the helical path.

It is possible that the photon-electron interaction is limited to photon radius larger than the electron orbit, ie rph>re.  The following graph shows why pushing the electron down is likely to cause ionization.


Given a 1re profile for PEre, an equal distance upwards and downwards, the electron pushed downwards gain more energy.

ΔPEΔr>ΔPEΔr+

 If this energy gain is greater than ionization energy, the electron will break the hold of the nucleus.

ΔPEΔr>ΔPEionionization

Pushing the electron out is not going to cause ionization unless Δre.

Friday, October 24, 2014

Happy 500!

An electron in B orbit around the positive nucleus experiences a repulsive force given by,

qq4πεor2eπqq4πεor2e=qq4εor2e

Two positive charges in parallel B orbit experience an attractive force given by,

Fa=qq4εor2p

Similarly, two negative charges in parallel B orbit experience an attractive force given by,

Fa=qq4εor2e

The charges are in orbit because of the initial speed it acquired due to electrostatic attraction between opposite charge.  They are held in orbit and kept from collapsing by a magnetic repulsive force.

What happen to two like charges then?  Do they coalesce?  When they are stopped in their B orbits, they will immediately repel apart. An instant disintegration!

This might be the rationale behind a particle collider.



Thursday, October 23, 2014

Passing It On...

From the post "Catastrophe! No Small Change",

UB=μo32(eπae)2(sin(2ϕ)sin(ϕ/2)dϕdt)2(θo01cos4(θ)h(ϕ,θ)d(cos(θ)))2

where,

h(ϕ,θ)=e2resin(ϕ/2)aecos(θ)

What interest us here is the term,

(sin(2ϕ)sin(ϕ/2)dϕdt)2

where ϕ is revolves along the orbit of the charge.  Imagine a cylinder defined by this motion at a height of  y perpendicular to the plane of the orbit.  The wall of such a cylinder cuts the B field established by the moving charge perpendicularly, at the highest point of the loop center at the circumference perpendicular to the tangent of the orbit.

Considering the B field alone, the force line does no work as the motion is in a perpendicular direction, and the charge in constant dϕdt around the orbit experience no change in kinetic energy.

But,

dUBdϕ=Addϕ{(sin(2ϕ)sin(ϕ/2)dϕdt)2(θo01cos4(θ)h(ϕ,θ)d(cos(θ)))2}

where A is the preceding constant and  θo is up to where the cylinder cuts the B field at y.

Consider

Ω=ddϕ{(sin(2ϕ)sin(ϕ/2)dϕdt)2(θo01cos4(θ)h(ϕ,θ)d(cos(θ)))2}

Consider the first term,

2sin(2ϕ)sin(ϕ/2)dϕdt{2cos(2ϕ)sin(ϕ/2)dϕdt+12sin(2ϕ)cos(ϕ/2)dϕdt}

={sin(4ϕ)sin2(ϕ/2)+12sin2(2ϕ)sin(ϕ)}(dϕdt)2

Consider the second term,

2θo01cos4(θ)h(ϕ,θ)d(cos(θ))ddϕθo01cos4(θ)h(ϕ,θ)d(cos(θ))

=2θo01cos4(θ)h(ϕ,θ)d(cos(θ))ddϕθo01cos4(θ)h(ϕ,θ)(cos(ϕ)sin3(θ)dϕ)

From the post "Catastrophe! No Small Change",

cos(ϕ)sin3(θ)dϕ=dcos(θ),    so,

=2θo01cos4(θ)h(ϕ,θ)d(cos(θ)).1cos4(θo)h(ϕ,θo)cos(ϕ)sin3(θo)

Therefore,

Ω=(sin(2ϕ)sin(ϕ/2)dϕdt)(θo01cos4(θ)h(ϕ,θ)d(cos(θ))).[dϕdt{2cos(2ϕ)sin(ϕ/2)+12sin(2ϕ)cos(ϕ/2)}(θo01cos4(θ)h(ϕ,θ)d(cos(θ)))+2cos4(θo)h(ϕ,θo)cos(ϕ)sin3(θo)(sin(2ϕ)sin(ϕ/2)dϕdt)]

And we have the change in energy with time as,

dUBdt=dUBdϕdϕdt=AΩdϕdt

Clearly  Ω0 for all ϕ.  This implies that the charge while in orbit constantly absorb and radiate energy.  A plot of,

 O=(sin(2ϕ)sin(ϕ/2))2  and sin(ϕ/2) below clearly shows the periodicity in UB of  2π.


And that the net change in UB over one period is zero as the graph returns to the same point (zero) at the end of the period.

Together with the previous post "Diffusion, Only If...",  we have a non zero expression for A(x,˙x),

Pt=12A(x,˙x)=dUBdt

This is the mechanism by which energy is propagated throughout a material.  The orbiting charges are constantly radiating and absorbing energy, with no net change over a period of the orbiting motion.

Disclaimer: Amateur At Work

I need to check my Maths...Hold your horses!  Still checking pressure lamp.

Diffusion, Only If...

Given,

L=kineticenergypotentialenergy=TV

Consider,

L˙x=Lxx˙x+Ltt˙x

as L=f(x,˙x,t)

Take derivative wrt t,

t{L˙x}=t{Lxx˙x+Ltt˙x}

Consider the second term,

t{Ltt˙x}=2Lt2t˙xLt1{˙xt}22˙xt2

Consider the first term,

t{Lxx˙x}=t{Lx}x˙x+Lxt{x˙x}

=2Ltxx˙x+Lx˙x{xt}

=2Ltxx˙x+Lx

For the Euler-Lagrange equation to be true,

t{L˙x}=Lx

to be true, we have,

2Ltxx˙x+2Lt2t˙xLt1{˙xt}22˙xt2=0

But,

2Ltx=2Lt˙x˙xx+2Lttt˙x˙xx

We have,

2Lt˙x+22Lt2t˙xLt1{˙xt}22˙xt2=0

Lx+22Lt2t˙xLt1{˙xt}22˙xt2=0

Lx=1¨x{x¨xLt22Lt2}

If  x=0, ie the forces in the fields do not change with time,

2Lt2=¨x2Lx

That going through a L gradient of Lx with acceleration ¨x, the second rate of change in L experienced by the body is given by the expression above.

Obviously when F=0,  ¨x=0  and  2Lt2=0 this implies,

Lt=A(x,˙x)

A is a function in x and ˙x not necessarily zero.  Consider the total energy of the system,

E=L+2V

Edt=Lt+2Pt

when ¨x=0,  ˙x is constant, T is a constant,

Et=0

 Tt=0

the change in L is solely due to the change in P

Pt=12Lt=12A(x,˙x)=14¨xLxdt ---(**)

An expression in ¨x is different from the value of ¨x.

d2Pdt2=14¨xLx=0

as we started with ¨x=0

This suggests that when a body is in circular motion (for example electrons, protons in orbit), ¨x in the direction of ˙x is zero, the centripetal force does no work and ˙x is a constant, there can still be a change in potential energy with time given by expression (**).  When

Pt>0    energy is absorbed,

Pt<0     energy is radiated

Is this the explanation why energy tend to diffuse throughout a material?  Because the orbiting particles are constantly absorbing and radiating energy?  Only if,

A(x,˙x)0.

In general this is possible when the system is periodic as in circular motion where over one period, the net gain in energy is zero.  Which brings us to ϕ from the post "Catastrophe! No Small Change"...until next time...

Taking My Time

Given any x,

xx=1

2xx2=0=2xxttx=2xtttxtx=2xt2txtx

So for all x,

2xt2txtx=0

ie,

av2=0

In fact,

2xx2=0

=x{xx}

=t{xx}tx

=t{xttx}tx

=t{xt}tx+t{tx}xt

=2xt2tx1{xt}2t{xt}xt

=2xt2tx1{xt}2xt2

=2xt2tx2xt2tx=0

non sequitur.

Wednesday, October 22, 2014

Orbit Small

From the post "Electron Orbit B Field II",

E=Eoexae

x is measured from the surface of the charge, where ae is the radius of the charge. So

Eoe0ae=Eo

Eo is on the surface of the charge.


It is then possible that  ae>re, as the measurement of x does not include ae.


Smooth Flow

From the post "Catastrophe No Small Change",

The graphform,

en(θ)=sin3(θ)cos2(θ)e4reaecos(θ)

defines the envelope of the radiated energy.

den(θ)dθ=tan(x)e4resec(x)ae(sin(x)+2tan(x)sec(x)4reaetan3(x))

A plot of this graph whose zeros are the maxima of the envelop, with reae varying from 1 to 2, is shown below.  There are three zeros, two points of inflection at 0 and π/2 and one maximum where the gradient of the derivative is negative.

For the maximum points,

sin(x)+2tan(x)sec(x)4reaetan3(x)=0

For the points of inflections,

sin(x)=0


Do not steal my money. No discontinuity.  We are dealing with radiated energy from varying magnetic field.  Previously when we encountered a discontinuity in the gradient of  re vs T we were dealing with electrostatic, centripetal force and thermal gravity, even gravity itself.

May be there is a similar bandgap, a magnetic bandgap that is not visible here.

Tuesday, October 21, 2014

No Need For Packets, All Eat Here

From the previous post "Catastrophe! No Small Change",

¯UB=μoω2e2ae8π3r3esin3(θo)cos2(θo)e4reaecos(θ0)

The exponential term is very sensitive to the relative values of  reae  and by reciprocity, rpap.  A illustrative plot of UB using varies ratio of reae is given below,  The actual graph plotted is (sin(x))^3*e^{-4*a/cos(x)}*1/(cos(x))^2 where  1<a<2, in steps of 0.05.


It seems that the charge's orbit is only slightly larger than its size in order to give the radiation profile obtained experimentally.

Monday, October 20, 2014

Catastrophe! No Small Change

If from the post "Electron Orbit B Field II",

Bdr=2πrBo=μoεoEAt

2πrBo=μoεo2πEocos(ϕ/2)dxdtx2aeθo01cos4(θ)h(ϕ,θ)d(cos(θ))

If we retain r and label it rp, the orbital radius of the nucleus then,

rpBo=μoεoEocos(ϕ/2)dxdtx2aeθo01cos4(θ)h(ϕ,θ)d(cos(θ))

rpBo=μoεo(e4πεor2e)cos(ϕ/2)dxdtx2aeθo01cos4(θ)h(ϕ,θ)d(cos(θ))

Since,

dxdtx2=recos(ϕ)dϕdt(2resin(ϕ/2))2

dxdtx2=4r3ecos(ϕ)dϕdtsin2(ϕ/2)

And so,

rpBo=μoεo(e4πεoaer2e)cos(ϕ/2)(4r3ecos(ϕ)dϕdtsin2(ϕ/2))θo01cos4(θ)h(ϕ,θ)d(cos(θ))

Bo=μoe4πaererpsin(2ϕ)sin(ϕ/2)dϕdtθo01cos4(θ)h(ϕ,θ)d(cos(θ))

Bo is generated by the electron in circular motion of radius re, so Bo is proportional to re.  Such a field decreases with distance from the plane containing the orbit of  the electron, so Bo is inversely proportional to rp.  Consider the energy contained in such a field.

UB=12B2oμo

When the nucleus has a positive charge or a net positive charge of one by reciprocity,

re=rp  and so,  (rerp)2=1

UB=12μo(μoe4πae)2(sin(2ϕ)sin(ϕ/2)dϕdt)2(θo01cos4(θ)h(ϕ,θ)d(cos(θ)))2

UB=μo32(eπae)2(sin(2ϕ)sin(ϕ/2)dϕdt)2(θo01cos4(θ)h(ϕ,θ)d(cos(θ)))2

When such a system is subjected to heat, both re and rp will oscillate and so does the magnetic field.  The value θo will also varies as rp changes.  When Bo collapses from a higher value to a lower value, energy is radiated.  This radiated energy varies with the value of θo, given by the expression above.

The average energy over one period, (power) is given by,

¯UB=1TT0UBdt=1T.T0μo32(eπae)2(sin(2ϕ)sin(ϕ/2)dϕdt)2(θo01cos4(θ)h(ϕ,θ)d(cos(θ)))2dt

¯UB=1Tμo32(eπae)2.ϕT0(sin(2ϕ)sin(ϕ/2))2(θo01cos4(θ)h(ϕ,θ)d(cos(θ)))2dϕdtdϕ

assuming that dϕdt=ω is a constant.

¯UB=ωTμo32(eπae)2ϕT0(sin(2ϕ)sin(ϕ/2))2(θo01cos4(θ)h(ϕ,θ)d(cos(θ)))2dϕ

This is the integral of two functions involving  ϕ since,

h(ϕ,θ)=e2resin(ϕ/2)aecos(θ)

Consider,

Ψ=(sin(2ϕ)sin(ϕ/2))2(θo01cos4(θ)h(ϕ,θ)d(cos(θ)))2dϕ

Ψ=(4sin(ϕ/2)cos(ϕ/2)cos(ϕ)sin(ϕ/2))2(θo01cos4(θ)h(ϕ,θ)d(cos(θ)))2dϕ

Ψ=(4sin2(ϕ/2)cos(ϕ))2(aere)2.(θo01cos3(θ)(h(ϕ,θ)recos(ϕ/2)aecos(θ))d(cos(θ)))2dϕ

Becasue,

h(ϕ,θ)=e2resin(ϕ/2)aecos(θ)

dh(ϕ,θ)dϕ=e2resin(ϕ/2)aecos(θ)(2reaecos(θ))cos(ϕ/2)12

dh(ϕ,θ)dϕ=h(ϕ,θ)recos(ϕ/2)aecos(θ)

We have,

Ψ=(4sin2(ϕ/2)cos(ϕ))2(aere)2(θo01cos3(θ)(dh(ϕ,θ)dϕ)d(cos(θ)))2dϕ

Ψ=16sin4(ϕ/2)cos2(ϕ)(aere)2(ddϕθo01cos3(θ)h(ϕ,θ)d(cos(θ)))2dϕ

From,

rpx=tan(θ)

x=2resin(ϕ/2)cos(ϕ/2)=resin(ϕ)

rp=resin(ϕ)tan(θ),butrp=re

cot(θ)=sin(ϕ)

sin(θ)dθ=cos(ϕ)sin3(θ)dϕ=dcos(θ) --- (*)

Despite of this relationship, ϕ  and  θ are still independent variable.  The integration over θ is taken for each value of  ϕ as ϕ revolves around the orbit defined by re.

Substitute (*) into Ψ,

Ψ=16(aere)2sin4(ϕ/2)cos2(ϕ)(ddϕθo01cos3(θ)h(ϕ,θ)cos(ϕ)sin3(θ)dϕ)2dϕ

Ψ=16(aere)2sin4(ϕ/2)cos2(ϕ)(1cos2(θ)h(ϕ,θ)sin3(θ))2dϕ

Ψ=16(aere)2sin6(θ)cos4(θ)sin4(ϕ/2)cos2(ϕ)(h(ϕ,θ))2dϕ

Ψ=4a3er3esin3(θ)cos2(θ).sin3(ϕ/2)cos(ϕ)(h(ϕ,θ))2(4resin(ϕ/2)ae1cos2(θ))cos(ϕ)sin3(θ)dϕ

Substitute into Ψ again,

Ψ=4a3er3esin3(θ)cos2(θ)sin3(ϕ/2)cos(ϕ)(h(ϕ,θ))2(4resin(ϕ/2)ae1cos2(θ))dcos(θ)

Ψ=4a3er3esin3(θ)cos2(θ)sin3(ϕ/2)cos(ϕ)(h(ϕ,θ))2(4resin(ϕ/2)aesin(θ)cos2(θ))dθ

Because,

h(ϕ,θ)2=e4resin(ϕ/2)aecos(θ)

d(h(ϕ,θ))2dθ=h(ϕ,θ)2(4resin(ϕ/2)aesin(θ)cos2(θ))

And so,

Ψ=4a3er3esin3(θ)cos2(θ)sin3(ϕ/2)cos(ϕ)d(h(ϕ,θ))2dθdθ

Ψ=4a3er3esin3(θ)cos2(θ)ddθsin3(ϕ/2)cos(ϕ)(h(ϕ,θ))2dθ

Ψ=4a3er3esin3(θ)cos2(θ)sin3(ϕ/2)cos(ϕ)(h(ϕ,θ))2

And the average energy over one period (power) in the B field is,

¯UB=ωTμo8(eπae)2a3er3e{sin3(θ)cos2(θ)}θo02{sin3(ϕ/2)cos(ϕ)(h(ϕ,θo))2}π0

¯UB=ω22πμo8(eπ)2aer3esin3(θo)cos2(θo)2{1.(h(π,θo))2}

¯UB=μoω2e2ae8π3r3esin3(θo)cos2(θo)(h(π,θo))2

¯UB=μoω2e2ae8π3r3esin3(θo)cos2(θo)e4reaecos(θ0)

A plot of (sin(x))^3*e^{-1.1/cos(x)}*1/(cos(x))^2  for  0<x<π/2 is given below,


There is no catastrophe.  When subjected to high perturbations,( high temperature etc), both orbits oscillate.  High θo corresponds to high orbit.  When the charge return from high values of θo, higher amount of energy are released as the potential energy difference between those at high orbits and that at the centers of the oscillation are higher.  So, high θo corresponds to higher radiated energy.

If this interpretation is correct then there is no need for discrete amount of energy either.

Sunday, October 19, 2014

Capillary Action And Camera!

Here's my take on capillary action,


The logic is naive; capillary action act against gravity, thermal gravity is anti gravity.  So, we set off to investigate temperature and capillary action.  The metal should be of high heat capacity such that the portion surrounding the insulating fluid does not cool to the temperature of the insulating fluid.

An insulating fluid is used to maintain the the temperature gradient between the metal and the fluid as the measurement progresses.

Increase Tm and plot Temperature of Metal vs Height of Fluid; Tm (K) vs h (m) with  d, diameter of metal tube bore and Tf temperature of fluid as parameters.  Viscosity plays a part as fluid in thermal contact with the metal wall experiences a net upward force due to thermal gravity and drags the fluid in the center of the bore upwards.

Maybe you can make a droplet of fluid flow upwards with a temperature gradient decreasing upwards on a hot surface with low drag.

Have a nice day.

And Mochi Sticks Too

All the preceding posts just to answer why the positive charges at the nucleus stick together.

They stick, because they are revolving in parallel B orbits. The attractive force between two orbiting charges in parallel orbits is,

FL=πq|v|B=iπqE

where q is the charge in one orbit, v the velocity of the charge, B the magnetic field established by the other orbiting charge and E the electrostatic field at distance r where B is, due to the other moving charge, and r is also the distance between the two orbits.  i rotates the force FL in the plane perpendicular to B in the anticlockwise direction.

FL is the result of two interacting B fields.  The maximum interaction occurs when the B field lines are parallel.   The velocity v of the charge placed in the established B field was made perpendicular to the already established B field lines such that the B field set up by this moving charge is parallel to the first.

In the derivation of FL in the post "Lorentz's without q", a moving negative charge was used that established a E field downwards.  A moving positive charge was then added to the situation to obtain FL.

If v in rotated, viv  to the same direction as the moving negative charge then FL is along  E as the expression shows.  This is consistent with the observed fact that parallel wires with opposite current repulse each other.

There is no relative velocity concerns here, as long as one charge is moving, it establishes a B field that will interact with another moving charge.

Have a nice day.

Stable But Not Too Stable

Continuing from the post "Tie Me Not, Not Too Tight",

2Bor2e=2μoεoωEoaef(ϕ)cot(θo){Ar2eBre+C}

where,

C=θo01cos4(θ)h(ϕ,θ)d(cos(θ))

B=4sin(ϕ/2)aeθo01cos5(θ)h(ϕ,θ)d(cos(θ))

A=2(sin(ϕ/2)ae)2θo01cos5(θ)h(ϕ,θ)d(cos(θ))

A>0  and  C>0  for  0<θ<π/2 but B depends on  ϕ/2.

B24AC=16(sin(ϕ/2)ae)2θo01cos5(θ)h(ϕ,θ)d(cos(θ))[θo01cos5(θ)h(ϕ,θ)d(cos(θ))12cot(θo)θo01cos4(θ)h(ϕ,θ)d(cos(θ))]

From the graph of Ar2eBre+C which determines the sign of the second derivative,


The second derivative is negative over a range of re, beyond this range it is positive.  Under normal circumstances at,

Bore=0

a small increase in re result in a negative change in Bore as

2Bor2e<0.

Since we started with zero, a negative change means that the derivative itself is negative,

Bore<0

and the system is stable.

However, if there is energy input that move re beyond this range, where the second derivative becomes positive, an increase in re will result in a positive Bore.  Bo increases and the repulsion between the spinning nucleus and the revolving electron increases causing re to increase further.  The electron breaks away and ionization occurs.

Similarly, if re is moved below the stable region, an decrease in re results in a positive Bore that decreases B correspondingly.  A reduction in repulsion between the spinning nucleus and the revolving electron decreases re further and the electron collapses into the nucleus.  We have plasma.

So, re has a range of stable values beyond which the particle breaks away at the increasing end or collapses at the decreasing end.

Energy input is required to move re.  This can be achieved by decreasing drag through increasing temperature, by pushing the orbiting electron inwards/outwards using photons at resonance, decreasing/increasing the B field between the nucleus and electron by applying an opposing/parallel B field, etc.

Orbit In Blue

The positive nucleus in orbit around the B orbit generates a similar B field that interacts with the electron in a reciprocal way.  The electron begins to spin because any transverse motion from its orbit around the nucleus will cut this B field and result in a Lorentz's force in the transverse plane (a plane in the direction perpendicular to the blue field lines and the relative motion).  This force acts perpendicular to the transverse velocity component and causes the electron path to curve.  Given a strong B field, the electron spins in the transverse plane.


Which give definition to the confining sphere around the electron.  The electron has two degrees of freedom on the sphere, one due to its orbit around the nucleus and the other its motion on the reciprocal orbit, Br.  The confining sphere may not be a sphere at all, but a confining volume.

Cart First Because The Cow Is Driving

If we were to make the nucleus stationary,


then the B orbit will be seen as revolving in and out of the confining spheres the surface of which demarcate the presence of the nucleus and the electron.

Believe It Or Not

Here is another illustration for  B orbit.


Movement is more relative motion between electron and the nucleus.  Although, the argument for a stationary nucleus because of heavier mass (inertia) is not valid because the forces involved are very high.

Saturday, October 18, 2014

Tie Me Not, Not Too Tight

For θo  such that Bo  is at the extrema, we look at

g(ϕ,θo)=cos(θo)sin(θo)θo01cos4(θ)h(ϕ,θ)d(cos(θ))

g(ϕ,θo)=cot(θo)θo01cos4(θ)h(ϕ,θ)d(cos(θ))

where,

h(ϕ,θ)=e2resin(ϕ/2)aecos(θ)

Consider,

gθo=θo{cot(θo)θo01cos4(θ)h(ϕ,θ)d(cos(θ))}

gθo=csc2(θo)θo01cos4(θ)h(ϕ,θ)d(cos(θ))+cot(θo)θoθo01cos4(θ)h(ϕ,θ)(sin(θ))d(θ)

gθo=csc2(θo)θo01cos4(θ)h(ϕ,θ)d(cos(θ))cos(θo)sin(θo)1cos4(θo)h(ϕ,θo)sin(θo)

gθo=csc2(θo)θo01cos4(θ)h(ϕ,θ)d(cos(θ))1cos3(θo)h(ϕ,θo)

Consider the extrema,

gθo=0

1sin2(θo)θo01cos4(θ)h(ϕ,θ)d(cos(θ))+1cos3(θo)h(ϕ,θo)=0

Since,

h(ϕ,θ)>0  and  g(ϕ,θo)>0 for 0<θo<π/2 and all ϕ

A plot of  (1/sin(x))^2 and (1/cos(x))^3 shows that,


which means by itself  g(ϕ,θo) has no extrema in the range  0<θo<π/2.  Since, its first derivative of  g(ϕ,θo) is negative  Bo  decreasse with  (θo).  The following diagram shows an orbiting electron driving a B field.

The electron approaches the nucleus because of electrostatic attract, before it collides with the nucleus, the electron reaches the terminal speed (light speed) and begins a helical path towards the nucleus.  Right over the nucleus, the relative motion of the revolving electron and the nucleus causes the nucleus to spin.  This creates a repulsive Lorentz's force that pushes the particles apart.  The electron and nucleus do not collide.  As a result of its spin, the nucleus interacts with the B field established by the revolving electron (initially part of its helical path) and revolve around a B orbit.  In turn, the revolving nucleus generates a B field that causes the electron to spin and revolve around a second  Br orbit (Refer to post "Orbit In Blue").  rp  do not effect  re in reciprocal.  rp is the result of a centripetal force provided for by the B field established by the electron in orbit re.  The nucleus will be accelerated to terminal speed and reaches a minimum rp value such that the centripetal force balances the Lorentz's force (B field) and drag at terminal speed.  The radius re affects rp; rp does not affects re.


From the post "Electron Orbit B Field II",

Bo=μoεoωr2eEoaef(ϕ)g(ϕ,θo)

We consider Bo at its maximum,

Bore=re{μoεoωr2eEoaef(ϕ)g(ϕ,θo)}

Bore=2μoεoωreEoaef(ϕ)g(ϕ,θo)+μoεoωr2eEoaef(ϕ)g(ϕ,θo)re

Bore=μoεoωreEoaef(ϕ)[2g(ϕ,θo)+reg(ϕ,θo)re]

g(ϕ,θo)re=cot(θo)θo01cos4(θ)h(ϕ,θ)red(cos(θ))

h(ϕ,θ)re=e2resin(ϕ/2)aecos(θ).2sin(ϕ/2)aecos(θ)=h(ϕ,θ)2sin(ϕ/2)aecos(θ)

Therefore,

g(ϕ,θo)re=cot(θo)θo01cos4(θ)h(ϕ,θ)2sin(ϕ/2)aecos(θ)d(cos(θ))

g(ϕ,θo)re=2sin(ϕ/2)aecot(θo)θo01cos5(θ)h(ϕ,θ)d(cos(θ))

So,

Bore=μoεoωreEoaef(ϕ)[2g(ϕ,θo)re2sin(ϕ/2)aecot(θo)θo01cos5(θ)h(ϕ,θ)d(cos(θ))]

For extrema,

Bore=0

When

aesin(ϕ/2)g(ϕ,θo)=recot(θo)θo01cos5(θ)h(ϕ,θ)d(cos(θ)) --- (*)

re=aesin(ϕ/2)θo01cos4(θ)h(ϕ,θ)d(cos(θ))θo01cos5(θ)h(ϕ,θ)d(cos(θ))

Consider the second derivative,

2Bor2e=2μoεoωEoaef(ϕ)g(ϕ,θo)+2μoεoωreEoaef(ϕ)g(ϕ,θo)re+2μoεoωreEoaef(ϕ)g(ϕ,θo)re+μoεoωr2eEoaef(ϕ)2g(ϕ,θo)r2e

2Bor2e=2μoεoωEoaef(ϕ)[g(ϕ,θo)+2reg(ϕ,θo)re+r2e22g(ϕ,θo)r2e]

2g(ϕ,θo)r2e=(2sin(ϕ/2)ae)2cot(θo)θo01cos5(θ)h(ϕ,θ)d(cos(θ))

So,

2Bor2e=2μoεoωEoaef(ϕ)[g(ϕ,θo)re4sin(ϕ/2)aecot(θo)θo01cos5(θ)h(ϕ,θ)d(cos(θ))+2r2e(sin(ϕ/2)ae)2cot(θo)θo01cos5(θ)h(ϕ,θ)d(cos(θ))]

2Bor2e=2μoεoωEoaef(ϕ)cot(θo)[θo01cos4(θ)h(ϕ,θ)d(cos(θ))re4sin(ϕ/2)aeθo01cos5(θ)h(ϕ,θ)d(cos(θ))+2r2e(sin(ϕ/2)ae)2θo01cos5(θ)h(ϕ,θ)d(cos(θ))]

2Bor2e=2μoεoωEoaef(ϕ)cot(θo){Ar2eBre+C}

where,

C=θo01cos4(θ)h(ϕ,θ)d(cos(θ))

B=4sin(ϕ/2)aeθo01cos5(θ)h(ϕ,θ)d(cos(θ))

A=2(sin(ϕ/2)ae)2θo01cos5(θ)h(ϕ,θ)d(cos(θ))

A>0  and  C>0  for  0<θ<π/2 but B depends on  ϕ/2.

B24AC=16(sin(ϕ/2)ae)2θo01cos5(θ)h(ϕ,θ)d(cos(θ))[θo01cos5(θ)h(ϕ,θ)d(cos(θ))12cot(θo)θo01cos4(θ)h(ϕ,θ)d(cos(θ))]

If B24AC>0 then 2Bor2e can be negative, since A>0 (part of the graph of Ar2eBre+C is below the x-axis).  However this depends on the value of  θo and cot(θo).  cot(θo) is a monotonously decreasing function, that means for larger values of  θo or rp,  Bo is maximum.

Applying the condition for extrema, from  (*)

aesin(ϕ/2)g(ϕ,θo)=recot(θo)θo01cos5(θ)h(ϕ,θ)d(cos(θ))

aeresin(ϕ/2)θo01cos4(θ)h(ϕ,θ)d(cos(θ))=θo01cos5(θ)h(ϕ,θ)d(cos(θ))

For a possible maxima in Bo, B24AC>0 implies,

\cfrac { a_{ e } }{ r_{ e }sin(\phi /2) } -\cfrac { 1 }{ 2 } cot(\theta _{ o })\gt0

cot(\theta _{ o })\lt\cfrac { 2a_{ e } }{ r_{ e } } \cfrac { 1 }{ sin(\phi /2) } 

\theta_o  must be larger than,

\theta_{ o }\gt cot^{ -1 }\left( \cfrac { 2a_{ e } }{ r_{ e } } \cfrac { 1 }{ sin(\phi /2) }  \right) 

where  a_e is the size of the electron and  r_e the electron orbit.  By sysmmetry, for the positive nucleus,

\theta_{ oe }\gt cot^{ -1 }\left( \cfrac { 2a_{ p } }{ r_{ p } } \cfrac { 1 }{ sin(\phi /2) }  \right) 

where  a_p  is the size of the nucleus and  r_p  the radius of its orbit and \theta_{oe} is the angle subtended by the electron's orbit and the moving nucleus as it revolves along it own orbit.

In order to have a maximum B_o which corresponds to a maximum centripetal force, the respective \thetas have minimum values.  A smaller \theta means a smaller radius that would require a higher centripetal force which B_o at maximum can not provide for.

The conclusion is, the orbit does not collapse.

From the fist derivative of B_o wrt r_e,

\cfrac { \partial \, B_{ o } }{ \partial \, r_{ e } } =\cfrac { \mu _{ o }\varepsilon _{ o }\omega r_{ e }E_{ o } }{ a_{ e } } f(\phi )\\\left[ 2g(\phi ,\theta _{ o })-r_{ e }\cfrac { 2sin(\phi /2) }{ a_{ e } } cot(\theta _{ o })\int _{ 0 }^{ \theta _{ o } }{ \cfrac { 1 }{ cos^{ 5 }(\theta ) } h(\phi ,\theta )\,  } { d\, (cos(\theta )) } \right]

cot(\theta _{ o }) appears in both g(\phi ,\theta _{ o }) and the negative term.  An increasing  r_e eventually results in a negative rate of change in B_o that decreases B_o.  A decrease in B_o reduces the repulsion between the spinning nucleus and the revolving electron. And r_e decreases and the electron does not fly off.  Similarly, a small r_e eventually results in a positive rate of change in B_o that increases B_o.  The repulsive force between the spinning nucleus and the revolving electron increases.  The electron is pushed away and does not collide into the nucleus.

Once again the conclusion is, the orbit does not collapse.

Friday, October 17, 2014

B Orbits Here To Stay

The expression,

 \int _{ 0 }^{ \theta _{ o } }{ \cfrac { 1 }{ cos^{ 4 }(\theta ) } e^{ -\cfrac { x }{ a_{ e }cos(\theta ) }  } } { d\, (cos(\theta )) }

is an interesting one.  And has solution,

z=cos(\theta)

\int _{ 0 }^{ \theta _{ o } }{ \cfrac { 1 }{ z^{ 4 } } e^{ -\cfrac { x }{ a_{ e }z }  } } { d\, z}=\cfrac { 1 }{ z^{ 4 } } \left\{ \cfrac { x }{ { a }_{ e } } { E }_{ i }(-\cfrac { x }{ { a }_{ e }z } )+xe^{ -\cfrac { x }{ { a }_{ e }z }  } \right\} +c

where  E_i(x)  is the polyexponential function,

 OR

\cfrac { a_{ e }(x^{ 2 }+2a_{ e }xz+2a_{ e }^{ 2 }z^{ 2 }) }{ x^{ 3 }z^{ 2 } } e^{ (-\cfrac { x }{ a_{ e }z } ) }

And so,

 \int _{ 0 }^{ \theta _{ o } }{ \cfrac { 1 }{ cos^{ 4 }(\theta ) } e^{ -\cfrac { x }{ a_{ e }cos(\theta ) }  } } { d\, (cos(\theta )) }=\left[\cfrac { a_{ e } }{ x^{ 3 }cos^{ 2 }(\theta ) } \left\{x^{ 2 }+2a_{ e }xcos(\theta )+2a_{ e }^{ 2 }cos^{ 2 }(\theta )\right\}e^{ (-\cfrac { x }{ a_{ e }cos(\theta ) } ) }\right]^{\theta_o}_{0}

Which has a DC term,

\cfrac { a_{ e } }{ x^{ 3 } } \left\{ x^{ 2 }+2a_{ e }x+2a_{ e }^{ 2 } \right\} e^{ (-\cfrac { x }{ a_{ e } } ) }

This term suggests that the B field is always present and the positively charged nucleuii are held permanently in the B orbits.

1/x^4*e^(-a/(2*x)) for 0.5≤a≤5,  1/x^4*e^(-2.5/(2*x)),  and  1/x^5*e^(-2.5/(2*x)) are plotted below.


Since,  x=cos(\theta), x\le1 the graph \cfrac{1}{x^5}\gt\cfrac{1}{x^4}.  And

\cfrac{1}{x^4}e^{a/(2x)} increases with decreasing a which in our original expression is x, the distance between the center of the loop and the orbiting negative charge.

Electron Orbit B Field II

This post is wrong.  Please refer to the post "Not To Be taken Too Seriously, Please" dated 14 May 2016. 

Continuing from the post "Electron Orbit B Field", we replace

E=\cfrac{e}{4\pi\varepsilon_o d^2}  with,

E=E_oe^{-\cfrac{x}{a_{e}}},

where a_e  is the radius of the electron, and

E_o=\cfrac{e}{4\pi\varepsilon_o a^2_e}

This exponential form of  E  was developed analogously, using exponential gravity as guide.


Total  E  over flat circular x-section when the charge is at a distance  x  from the loop,

E_A=\int_{0}^{r} {2\pi xtan(\theta). E_oe^{-\cfrac{x}{a_{e}cos(\theta)} } }d\,r

E_{ A }=2\pi E_{ o }\int _{ 0 }^{ r }{ xtan(\theta )e^{-\cfrac{x}{a_{e}cos(\theta)} } } d\, r

\cfrac { \partial \, E_{ A } }{ \partial t } =2\pi E_{ o }\cfrac { \partial \, }{ \partial t }\int _{ 0 }^{ r }{ xtan(\theta )e^{-\cfrac{x}{a_{e}cos(\theta)} } } d\, r

Cosinder the term,

\Pi=\cfrac { \partial \, }{ \partial t }\int _{ 0 }^{ r }{ xtan(\theta )e^{-\cfrac{x}{a_{e}cos(\theta)} } } d\, r

\Pi=\int _{ 0 }^{ r }{ \cfrac { d\, x }{ d\, t } tan(\theta )e^{-\cfrac{x}{a_{e}cos(\theta)} }} \\+xtan(\theta )e^{-\cfrac{x}{a_{e}cos(\theta)} }\cfrac { -1 }{ a_{ e }cos(\theta ) } \cfrac { d\, x }{ d\, t } \\+xsec^{ 2 }(\theta )\cfrac { d\, \theta  }{ d\, t } e^{-\cfrac{x}{a_{e}cos(\theta)} }d\, r

\Pi=\int _{ 0 }^{ r }{ e^{-\cfrac{x}{a_{e}cos(\theta)} }\left\{ { \cfrac { d\, x }{ d\, t } tan(\theta ) }-xtan(\theta )\cfrac { 1 }{ a_{ e }cos(\theta ) } \cfrac { d\, x }{ d\, t } +xsec^{ 2 }(\theta )\cfrac { d\, \theta  }{ d\, t }  \right\}  } d\, r

Since,

\cfrac { r }{ x } =tan(\theta )

-\cfrac { r }{ x^{ 2 } } \cfrac { d\, x }{ d\, t } =sec^{ 2 }(\theta )\cfrac { d\, \theta  }{ d\, t }

xsec^{ 2 }(\theta )\cfrac { d\, \theta  }{ d\, t }=-\cfrac { r }{ x } \cfrac { d\, x }{ d\, t }=-\cfrac { d\, x }{ d\, t }tan(\theta)

Substitute into the expression for  \Pi,

\Pi=\int _{ 0 }^{ r }{ e^{-\cfrac{x}{a_{e}cos(\theta)} }\left\{ { \cfrac { d\, x }{ d\, t } tan(\theta ) }-xtan(\theta )\cfrac { 1 }{ a_{ e }cos(\theta ) } \cfrac { d\, x }{ d\, t } -\cfrac { d\, x  }{ d\, t }tan(\theta )  \right\}  } d\, r

\Pi=-\cfrac { d\, x }{ d\, t } \cfrac { x }{ a_{ e } } \int _{ 0 }^{ r }{ tan(\theta )sec(\theta )e^{-\cfrac{x}{a_{e}cos(\theta)} } } d\, r

and

\cfrac { r }{ x } =tan(\theta )

\cfrac { 1 }{ x } =sec^{ 2 }(\theta )\cfrac { d\, \theta  }{ dr }

dr=xsec^{ 2 }(\theta ){ d\, \theta  }

\Pi=-\cfrac { d\, x }{ d\, t } \cfrac { x }{ a_{ e } } \int _{ 0 }^{ r }{ tan(\theta )sec(\theta )e^{-\cfrac{x}{a_{e}cos(\theta)} } }xsec^{ 2 }(\theta ){ d\, \theta  } 

Substitute  \Pi  back into  E_A,

\cfrac { \partial \, E_{ A } }{ \partial t } =-2\pi E_{ o }\cfrac { d\, x }{ d\, t } \cfrac { x }{ a_{ e } } \int _{ 0 }^{ \theta _{ o } }{ tan(\theta )sec(\theta )e^{-\cfrac{x}{a_{e}cos(\theta)} }} xsec^{ 2 }(\theta ){ d\, \theta  }

\cfrac { \partial \, E_{ A } }{ \partial t } =-2\pi E_{ o }\cfrac { d\, x }{ d\, t } \cfrac { x^{ 2 } }{ a_{ e } } \int _{ 0 }^{ \theta _{ o } }{ sec^{ 4 (\theta )}e^{-\cfrac{x}{a_{e}cos(\theta)} } } sin(\theta ){ d\, \theta  }

\cfrac { \partial \, E_{ A } }{ \partial t } =2\pi E_{ o }\cfrac { d\, x }{ d\, t } \cfrac { x^{ 2 } }{ a_{ e } } \int _{ 0 }^{ \theta _{ o } }{ \cfrac{1}{cos^{ 4 }(\theta )}e^{-\cfrac{x}{a_{e}cos(\theta)} } } { d\, (cos(\theta )) } --- (*)

This however, is not the final solution.  We have to consider the fact that the charge is in orbit of radius  r_e  and  the loop defining  B is centered at the circumference along a radial line.


We are concerned with change in  E perpendicular to the loop.

\sigma=\cfrac{\pi-\phi}{2}=\cfrac{\pi}{2}-\cfrac{\phi}{2}

E_{\bot}=Esin(\sigma)=Esin(\cfrac{\pi}{2}-\cfrac{\phi}{2})=Ecos(\phi/2)

Consider,

x^2=r^2_e+r^2_e-2r^2_ecos(\phi)=2r^2_e(1-cos(\phi)),  

x=r_e\sqrt{2(1-cos(\phi))}=2r_{ e }sin(\phi /2)

\cfrac { dx }{ dt }=r_{ e }cos(\phi/2)\cfrac{d\,\phi}{d\,t}

and,

x_\bot=xcos(\phi/2)

\cfrac { dx_\bot }{ dt }=\cfrac { dx }{ dt }cos(\phi/2)-xsin(\phi/2).\cfrac{1}{2}\cfrac { d\phi }{ dt }  

\cfrac { dx_{ \bot  } }{ dt } =\left( r_{ e }cos^{ 2 }(\phi /2)-xsin(\phi /2).\cfrac { 1 }{ 2 }  \right) \cfrac { d\phi  }{ dt } 

Substitute x into the above.

\cfrac { dx_{ \bot  } }{ dt } =r_{ e }cos(\phi )\cfrac { d\phi  }{ dt } 

Consider the term,

e^{ -\cfrac { x }{ a_{ e }cos(\theta ) }  }=e^{ -\cfrac { r_{ e }\sqrt { 2(1-cos(\phi )) }  }{ a_{ e }cos(\theta ) }  }=e^{ -\cfrac { 2r_{ e }sin(\phi /2) }{ a_{ e }cos(\theta ) }  }

Let

h(\phi,\theta)=e^{ -\cfrac { 2r_{ e }sin(\phi /2) }{ a_{ e }cos(\theta ) }  }

which has both  \phi  and  \theta  dependence.

From (*),

\cfrac { \partial \, E_{ A } }{ \partial t } =2\pi E_{ o }\cfrac { d\, x }{ d\, t } \cfrac { x^{ 2 } }{ a_{ e } } \int _{ 0 }^{ \theta _{ o } }{ \cfrac{1}{cos^{ 4 }(\theta )}h(\phi,\theta) } { d\, (cos(\theta )) }

Adjusted for circular orbit,

\cfrac { \partial \, E_{ A } }{ \partial t } =2\pi E_{ o }cos(\phi/2 )\cfrac { d\, x_{ \bot  } }{ d\, t } \cfrac { x^{ 2 }_{ \bot  } }{ a_{ e } } \int _{ 0 }^{ \theta _{ o } }{ \cfrac { 1 }{ cos^{ 4 }(\theta ) } h(\phi,\theta) } { d\, (cos(\theta )) }

The direction of  E  is adjusted by  cos(\phi/2).  The the distance term in the exponential remains as  x as we see from the above diagram.

Assuming that  B  is constant around the loop,

\oint{B}d\,r=2\pi rB_o

2\pi rB_o=\mu _{ o }\varepsilon _{ o }2\pi E_{ o }cos(\phi /2)\cfrac { d\, x_{ \bot  } }{ d\, t } \cfrac { x^{ 2 }_{ \bot  } }{ a_{ e } } \int _{ 0 }^{ \theta _{ o } }{ \cfrac { 1 }{ cos^{ 4 }(\theta ) } h(\phi,\theta)} { d\, (cos(\theta )) }

This is wrong. r\ne x_\bot tan(\theta _{ o }) refer to post "Not To Be taken Too Seriously, Please" dated 14 May 2016 

x_\bot tan(\theta _{ o })B_{ o }=\cfrac {\mu _{ o }\varepsilon _{ o } E_{ o }  }{ a_{ e } }cos(\phi/2 )\cfrac { d\, x_{ \bot  } }{ d\, t } x^{ 2 }_{ \bot  } \int _{ 0 }^{ \theta _{ o } }{ \cfrac { 1 }{ cos^{ 4 }(\theta ) } h(\phi,\theta)} { d\, (cos(\theta )) }

B_{ o }=\cfrac { \mu _{ o }\varepsilon _{ o }E_{ o } }{ a_{ e } }cos(\phi/2 )\cfrac { d\, x_{ \bot  } }{ d\, t } x_{ \bot  } \cfrac {cos(\theta_o) }{ sin(\theta _{ o }) } \int _{ 0 }^{ \theta _{ o } }{ \cfrac { 1 }{ cos^{ 4 }(\theta ) } h(\phi,\theta) } { d\, (cos(\theta )) }

Consider,

\cfrac { d\, x_{ \bot  } }{ d\, t } x_{ \bot  }=2r^{ 2 }_{ e }cos(\phi /2)sin(\phi/2)cos(\phi )\cfrac { d\phi  }{ dt }  

Let

f(\phi)=cos(\phi /2).2cos(\phi /2)sin(\phi/2)cos(\phi )=\cfrac{1}{2}sin(2\phi)cos(\phi /2)

such that,

cos(\phi /2)\cfrac { d\, x_{ \bot  } }{ d\, t } x_{ \bot  } =r^2_ef(\phi)\cfrac { d\phi  }{ dt }

And we let,

g(\phi,\theta_o)= \cfrac { cos(\theta _{ o }) }{ sin(\theta _{ o }) } \int _{ 0 }^{ \theta _{ o } }{ \cfrac { 1 }{ cos^{ 4 }(\theta ) } h(\phi,\theta)} { d\, (cos(\theta )) }

and

\omega=\cfrac{d\,\phi}{d\,t}

B_o  becomes,

B_{ o }=\cfrac {\mu _{ o }\varepsilon _{ o }E_{ o } }{ a_{ e } } cos(\phi /2)\cfrac { d\, x_{ \bot  } }{ d\, t }x_{ \bot  }  \cfrac {cos(\theta_o) }{ sin(\theta _{ o }) } \int _{ 0 }^{ \theta _{ o } }{ \cfrac { 1 }{ cos^{ 4 }(\theta ) } h(\phi,\theta) } { d\, (cos(\theta )) }

B_{ o }=\cfrac { \mu _{ o }\varepsilon _{ o }\omega r^{ 2 }_{ e }E_{ o } }{ a_{ e } } f(\phi )g(\phi ,\theta _{ o })

where the term  g(\phi,\theta_o)  share the variable  \phi  with  f(\phi)  through  h(\phi,\theta_o)  and the two are not independent.

A plot of  f(\phi)  is shown below.  f(\phi) is periodic and bounded.


This shows that the used of the exponential form is appropriate for  E.

Thursday, October 16, 2014

Electron Orbit B Field..BLEWUP

What is the  B  field established by an electron in orbit.

From Maxwell,

\oint { B } \,dr=\varepsilon _{ o }\mu _{ o }\cfrac { \partial \,  }{ \partial t } \left\{ \oint { E } \,dA \right\} 


E=\cfrac{e}{4\pi\varepsilon_o d^2}


Total  E  over flat circular x-section when the charge is at a distance  x  from the loop,

E_A=\int_{0}^{r} {2\pi xtan(\theta).\cfrac{e}{4\pi\varepsilon_o}\left(\cfrac{cos(\theta)}{x}\right)^2}d\,r

\cfrac { \partial \, E_{ A } }{ \partial t } =\cfrac { e }{ 2\varepsilon _{ o } } \cfrac { \partial \,  }{ \partial t } \left\{ \int _{ 0 }^{ r }{ sin(\theta )cos(\theta ).\cfrac { 1 }{ x }  } d\, r \right\}

\cfrac { \partial \, E_{ A } }{ \partial t } =\cfrac { e }{ 4\varepsilon _{ o } } \cfrac { \partial \,  }{ \partial t } \left\{ \int _{ 0 }^{ r }{ sin(2\theta )\cfrac { 1 }{ x }  } d\, r \right\}

\cfrac { \partial \, E_{ A } }{ \partial t } =\cfrac { e }{ 4\varepsilon _{ o } } \left\{ \int _{ 0 }^{ r }{ 2cos(2\theta )\cfrac { 1 }{ x } \cfrac { d\, \theta  }{ d\, t }  } -sin(2\theta )\cfrac { 1 }{ x^{ 2 } } \cfrac { d\, x }{ d\, t } d\, r \right\}

\cfrac { \partial \, E_{ A } }{ \partial t } =\cfrac { e }{ 2\varepsilon _{ o } } \left\{ \int _{ 0 }^{ r }{ (2cos^{ 2 }(\theta )-1)\cfrac { 1 }{ x } \cfrac { d\, \theta  }{ d\, t }  } -sin(\theta )cos(\theta )\cfrac { 1 }{ x^{ 2 } } \cfrac { d\, x }{ d\, t } d\, r \right\} 

From geometry,

\cfrac { r }{ x } =tan(\theta )

-\cfrac { r }{ x^{ 2 } } \cfrac { d\, x }{ d\, t } =sec^{ 2 }(\theta )\cfrac { d\, \theta  }{ d\, t }

\cfrac { d\, \theta  }{ d\, t } =-\cfrac { rcos^{ 2 }(\theta ) }{ x^{ 2 } } \cfrac { d\, x }{ d\, t }

and

\cfrac { 1 }{ x } =sec^{ 2 }(\theta )\cfrac { d\, \theta  }{ dr }

dr=xsec^{ 2 }(\theta ){ d\, \theta  }

So we have,

\cfrac { \partial \, E_{ A } }{ \partial t } =\cfrac { e }{ 2\varepsilon _{ o } } \left\{ \int _{ 0 }^{ r }{ -(2cos^{ 2 }(\theta )-1)\cfrac { 1 }{ x } \cfrac { rcos^{ 2 }(\theta ) }{ x^{ 2 } } \cfrac { d\, x }{ d\, t }  } -sin(\theta )cos(\theta )\cfrac { 1 }{ x^{ 2 } } \cfrac { d\, x }{ d\, t } d\, r \right\}

\cfrac { \partial \, E_{ A } }{ \partial t } =\cfrac { e }{ 2\varepsilon _{ o } } \cfrac { 1 }{ x^{ 2 } } \cfrac { d\, x }{ d\, t } \left\{ \int _{ 0 }^{ r }{ -(2cos^{ 2 }(\theta )-1)tan(\theta )cos^{ 2 }(\theta )-sin(\theta )cos(\theta ) } d\, r \right\}

\cfrac { \partial \, E_{ A } }{ \partial t } =-\cfrac { e }{ 2\varepsilon _{ o } } \cfrac { 1 }{ x^{ 2 } } \cfrac { d\, x }{ d\, t } \left\{ \int _{ 0 }^{ r }{ 2cos^{ 2 }(\theta )sin(\theta )cos(\theta ) } d\, r \right\}

\cfrac { \partial \, E_{ A } }{ \partial t } =-\cfrac { e }{ 2\varepsilon _{ o } } \cfrac { 1 }{ x } \cfrac { d\, x }{ d\, t } \left\{ \int _{ 0 }^{ \theta _{ o } }{ cos^{ 2 }(\theta )sin(2\theta ) } sec^{ 2 }(\theta )d\, \theta  \right\}

\cfrac { \partial \, E_{ A } }{ \partial t } =-\cfrac { e }{ 2\varepsilon _{ o } } \cfrac { 1 }{ x } \cfrac { d\, x }{ d\, t } \left\{ \int _{ 0 }^{ \theta _{ o } }{ sin(2\theta ) } d\, \theta  \right\}

\cfrac { \partial \, E_{ A } }{ \partial t } =\cfrac { e }{ 4\varepsilon _{ o } } \cfrac { 1 }{ x } \cfrac { d\, x }{ d\, t } \left\{ cos(2\theta _{ o })-1 \right\} --- (*)

We have also to consider the fact that the charge is in orbit of radius  r_e  and  the loop defining  B is centered at the circumference along a radial line.


\sigma=\cfrac{\pi-\phi}{2}=\cfrac{\pi}{2}-\cfrac{\phi}{2}

sin(\sigma)=sin(\cfrac{\pi}{2}-\cfrac{\phi}{2})=cos({\phi}/{2})

E_\bot=Esin(\sigma)=Ecos(\phi/2)

x_\bot=xsin(\sigma)=xcos(\phi/2)

We are concerned with the change in E perpendicular to the smaller orbit in the loop.   From (*)

\cfrac { \partial \, E_{ A } }{ \partial t } =\cfrac { e }{ 2\varepsilon _{ o } } \cfrac { d\, x }{ d\, t } \cfrac { 1 }{ x } \left\{ cos(2\theta _{ o })-1 \right\} 

 So the equation (*) becomes,

\cfrac { \partial \, E_{ A } }{ \partial t } =\cfrac { e }{ 2\varepsilon _{ o } } sin(\sigma )\cfrac { d\, x_{ \bot  } }{ d\, t } \cfrac { 1 }{ x }  \left\{ cos(2\theta _{ o })-1 \right\} 

 \cfrac { \partial \, E_{ A } }{ \partial t } =\cfrac { e }{ 2\varepsilon _{ o } }cos(\phi/2)\cfrac { d\, x_{ \bot  } }{ d\, t }\cfrac { 1 }{ x } \left\{ cos(2\theta _{ o })-1 \right\}   --- (1)

The perpendicular E field component is given by E_\bot=Ecos(\phi/2), not by changing  x  to x_\bot.

Consider,

x^2=r^2_e+r^2_e-2r^2_ecos(\phi)=2r^2_e(1-cos(\phi)),  

x=r_e\sqrt{2(1-cos(\phi))}=2r_{ e }sin(\phi /2)

\cfrac { dx }{ dt }=r_{ e }cos(\phi/2)\cfrac{d\,\phi}{d\,t}

and,

x_\bot=xcos(\phi/2)

\cfrac { dx_\bot }{ dt }=\cfrac { dx }{ dt }cos(\phi/2)-xsin(\phi/2).\cfrac{1}{2}\cfrac { d\phi }{ dt }  

\cfrac { dx_{ \bot  } }{ dt } =\left( r_{ e }cos^{ 2 }(\phi /2)-xsin(\phi /2).\cfrac { 1 }{ 2 }  \right) \cfrac { d\phi  }{ dt } 

Substitute x into the above.

\cfrac { dx_{ \bot  } }{ dt } =r_{ e }cos(\phi )\cfrac { d\phi  }{ dt } 

\cfrac { dx_{ \bot  } }{ dt } \cfrac { 1 }{ x } =r_{ e }\cfrac { 1 }{ 2r_{ e }sin(\phi /2) } cos(\phi )\cfrac { d\phi  }{ dt } 

\cfrac { dx_{ \bot  } }{ dt } \cfrac { 1 }{ x } =\cfrac { 1 }{ 2 } \cfrac { cos(\phi ) }{ sin(\phi /2) } \cfrac { d\phi  }{ dt } 

From (1),

 \cfrac { \partial \, E_{ A } }{ \partial t } =\cfrac { e }{ 2\varepsilon _{ o } }cos(\phi/2)\cfrac { d\, x_{ \bot  } }{ d\, t }\cfrac { 1 }{ x }  \left\{ cos(2\theta _{ o })-1 \right\} 

\cfrac { \partial \, E_{ A } }{ \partial t } =\cfrac { e }{ 4\varepsilon _{ o } } cos(\phi /2)\cfrac { 1 }{ 2 } \cfrac { cos(\phi ) }{ sin(\phi /2) } \cfrac { d\phi  }{ dt }  \left\{ cos(2\theta _{ o })-1 \right\} 

We know at this point that B is going to blowup when integrated over  \phi.  This suggests that  \cfrac{1}{x^2} is not approprite here.

Assuming that  B  is constant around the loop,

\oint{B}d\,r=2\pi r_pB_o

where r_p is the radius of the B orbit.

2\pi r_pB_o=\mu _{ o }\varepsilon _{ o }\cfrac { \partial \, E_{ A } }{ \partial t } =\cfrac { \mu _{ o }e }{ 4 } cos(\phi /2)\left\{ \cfrac { 1 }{ sin(\phi /2) } -2sin(\phi /2) \right\} \cfrac { d\phi  }{ dt } \left\{ cos(2\theta _{ o })-1 \right\}

B_{ o }=\cfrac { \mu _{ o }e }{ 8\pi r_{ p } } cos(\phi /2)\left\{ \cfrac { 1 }{ sin(\phi /2) } -2sin(\phi /2) \right\} \cfrac { d\phi  }{ dt } \left\{ cos(2\theta _{ o })-1 \right\}

This is a time varying B  field, changing as the electron travels at  \omega=\cfrac{d\phi}{d\,t} around its orbit.  We can average B^2 over one period to obtain the energy in needed to establish the orbit.

 \cfrac { 1 }{ T } \int _{ 0 }^{ T }{ B^{ 2 }_{ o } } d\, t=\cfrac { 1 }{ T } \int _{ 0 }^{ T }{ \left( \cfrac { \mu _{ o }e }{ 8\pi r_{ p } } sin^{ 2 }(\theta _{ o }) \right) ^{ 2 }\left( cos(\phi /2)\left\{ \cfrac { 1 }{ sin(\phi /2) } -2sin(\phi /2) \right\}  \right) ^{ 2 }\cfrac { d\phi  }{ dt } \cfrac { d\phi  }{ dt }  }

=\cfrac { 1 }{ T } \left( \cfrac { \mu _{ o }e }{ 8\pi r_{ p } } sin^{ 2 }(\theta _{ o }) \right) ^{ 2 }\int _{ 0 }^{ 2\pi }{ \left( cos(\phi /2)\left\{ \cfrac { 1 }{ sin(\phi /2) } -2sin(\phi /2) \right\}  \right) ^{ 2 }\omega  } d\phi

as \omega=\cfrac{d\phi}{d\,t}.

Consider,

\int _{ 0 }^{ 2\pi}{ \left( cos(\phi /2)\left\{ \cfrac { 1 }{ sin(\phi /2) } -2sin(\phi /2) \right\}  \right) ^{ 2 }\omega  } d\phi 

KABOOM!!