Could second generation of weak bosons explain the reduction of proton charge radius?

The discovery by Pohl et al (2010) was that the charge radius of proton deduced from the muonic version of hydrogen atom - is .842 fm and about 4 per cent smaller than .875 fm than the charge radius deduced from hydrogen atom is in complete conflict with the cherished belief that atomic physics belongs to the museum of science (for details see the Wikipedia article). The title of the article Quantum electrodynamics-a chink in the armour? of the article published in Nature expresses well the possible implications, which might actually go well extend beyond QED.

Quite recently (2016) new more precise data has emerged from Pohl et al (see this). Now the reduction of charge radius of muonic variant of deuterium is measured. The charge radius is reduced from 2.1424 fm to 2.1256 fm and the reduction is .012 fm, which is about .8 per cent (see this). The charge radius of proton deduced from it is reported to be consistent with the charge radius deduced from deuterium. The anomaly seems therefore to be real. Deuterium data provide a further challenge for various models.

The finding is a problem of QED or to the standard view about what proton is. Lamb shift is the effect distinguishing between the states hydrogen atom having otherwise the same energy but different angular momentum. The effect is due to the quantum fluctuations of the electromagnetic field. The energy shift factorizes to a product of two expressions. The first one describes the effect of these zero point fluctuations on the position of electron or muon and the second one characterizes the average of nuclear charge density as "seen" by electron or muon. The latter one should be same as in the case of ordinary hydrogen atom but it is not. Does this mean that the presence of muon reduces the charge radius of proton as determined from muon wave function? This of course looks implausible since the radius of proton is so small. Note that the compression of the muon's wave function has the same effect.

Before continuing it is good to recall that QED and quantum field theories in general have difficulties with the description of bound states: something which has not received too much attention. For instance, van der Waals force at molecular scales is a problem. A possible TGD based explanation and a possible solution of difficulties proposed for two decades ago is that for bound states the two charged particles (say nucleus and electron or two atoms) correspond to two 3-D surfaces glued by flux tubes rather than being idealized to points of Minkowski space. This would make the non-relativistic description based on Schrödinger amplitude natural and replace the description based on Bethe-Salpeter equation having horrible mathematical properties.

The basic idea of the original model of the anomaly (see this) is that muon has some probability to end up to the magnetic flux tubes assignable to proton. In this state it would not contribute to the ordinary Schrödinger amplitude. The effect of this would be reduction of |Ψ|2 near origin and apparent reduction of the charge radius of proton. The weakness of the model is that it cannot make quantitative prediction for the size of the effect. Even the sign is questionable. Only S-wave binding energy is affected considerably but does the binding energy really increase by the interaction of muon with the quarks at magnetic flux tubes? Is the average of the charge density seen by muon in S wave state larger, in other words does it spend more time near proton or do the quarks spend more time at the flux tubes?

In the following a new model for the anomaly will be discussed.

  1. The model is inspired by data about breaking of universality of weak interactions in neutral B decays possibly manifesting itself also in the anomaly in the magnetic moment of muon. Also the different values of the charge radius deduced from hydrogen atom and muonium could reflect the breaking of universality. In the original model the breaking of universality is only effective.
  2. TGD indeed predicts a dynamical U(3) gauge symmetry whose 8+1 gauge bosons correspond to pairs of fermion and anti-fermion at opposite throats of wormhole contact. Throats are characterized by genus g=0,1,2, so that bosons are superpositions of states labelled by (g1,g2). Fermions correspond to single wormhole throat carrying fermion number and behave as U(3) triplet labelled by g.

    The charged gauge bosons with different genera for wormhole throats are expected to be very massive. The 3 neutral gauge bosons with same genus at both throats are superpositions of states (g,g) are expected to be lighter. Their charge matrices are orthogonal and necessarily break the universality of electroweak interactions. For the lowest boson family - ordinary gauge bosons - the charge matrix is proportional to unit matrix. The exchange of second generation bosons Z01 and γ1 would give rise to Yukawa potential increasing the binding energies of S-wave states. Therefore Lamb shift defined as difference between energies of S and P waves is increased and the charge radius deduced from Lamb shift becomes smaller.

  3. The model thus predicts a correct sign for the effect but the size of the effect from naive estimate assuming only γ2 and α21== α for M=2.9 TeV is almost by an order of magnitude too small. The values of the gauge couplings α2 and αZ,2 are free parameters as also the mixing angles between states (g,g). The effect is also proportional to the ratio (mμ/M(boson)2. It turns out that the inclusion of Z01 contribution and assumption α1 and αZ,1 are near color coupling strength αs gives a correct prediction.
Motivations for the breaking of electroweak universality

The anomaly of charge radius could be explained also as breaking of the universality of weak interactions. Also other anomalies challenging the universality exists. The decays of neutral B-meson to lepton pairs should be same apart from corrections coming from different lepton masses by universality but this does not seem to be the case (see this). There is also anomaly in muon's magnetic moment discussed briefly here. This leads to ask whether the breaking of universality could be due to the failure of universality of electroweak interactions.

The proposal for the explanation of the muon's anomalous magnetic moment and anomaly in the decays of B-meson is inspired by a recent very special di-electron event and involves higher generations of weak bosons predicted by TGD leading to a breaking of lepton universality. Both Tommaso Dorigo (see this) and Lubos Motl (see this) tell about a spectacular 2.9 TeV di-electron event not observed in previous LHC runs. Single event of this kind is of course most probably just a fluctuation but human mind is such that it tries to see something deeper in it - even if practically all trials of this kind are chasing of mirages.

Since the decay is leptonic, the typical question is whether the dreamed for state could be an exotic Z boson. This is also the reaction in TGD framework. The first question to ask is whether weak bosons assignable to Mersenne prime M89 have scaled up copies assignable to Gaussian Mersenne M79. The scaling factor for mass would be 2(89-79)/2= 32. When applied to Z mass equal to about .09 TeV one obtains 2.88 TeV, not far from 2.9 TeV. Eureka!? Looks like a direct scaled up version of Z!? W should have similar variant around 2.6 TeV.

TGD indeed predicts exotic weak bosons and also gluons.

  1. TGD based explanation of family replication phenomenon in terms of genus-generation correspondence forces to ask whether gauge bosons identifiable as pairs of fermion and antifermion at opposite throats of wormhole contact could have bosonic counterpart for family replication. Dynamical SU(3) assignable to three lowest fermion generations labelled by the genus of partonic 2-surface (wormhole throat) means that fermions are combinatorially SU(3) triplets. Could 2.9 TeV state - if it would exist - correspond to this kind of state in the tensor product of triplet and antitriplet? The mass of the state should depend besides p-adic mass scale also on the structure of SU(3) state so that the mass would be different. This difference should be very small.
  2. Dynamical SU(3) could be broken so that wormhole contacts with different genera for the throats would be more massive than those with the same genera. This would give SU(3) singlet and two neutral states, which are analogs of η' and η and π0 in Gell-Mann's quark model. The masses of the analogs of η and π0 and the the analog of η', which I have identified as standard weak boson would have different masses. But how large is the mass difference?
  3. These 3 states are expected top have identical mass for the same p-adic mass scale, if the mass comes mostly from the analog of hadronic string tension assignable to magnetic flux tube. connecting the two wormhole contacts associates with any elementary particle in TGD framework (this is forced by the condition that the flux tube carrying monopole flux is closed and makes a very flattened square shaped structure with the long sides of the square at different space-time sheets). p-Adic thermodynamics would give a very small contribution genus dependent contribution to mass if p-adic temperature is T=1/2 as one must assume for gauge bosons (T=1 for fermions). Hence 2.95 TeV state could indeed correspond to this kind of state.
The sign of the effect is predicted correctly and the order of magnitude come out correctly

Could the exchange of massive MG,79 photon and Z0 give rise to additional electromagnetic interaction inducing the breaking of Universality? The first observation is that the binding energy of S-wave state increases but there is practically no change in the energy of P wave state. Hence the effective charge radius rp as deduced from the parameterization of binding energy different terms of proton charge radius indeed decreases.

Also the order of magnitude for the effect must come out correctly.

  1. The additional contribution in the effective Coulomb potential is Yukawa potential. In S-wave state this would give a contribution to the binding energy in a good approximation given by the expectation value of the Yukawa potential, which can be parameterized as

    V(r)= g2 e-Mr/r ,&g2 = 4π kα .

    The expectation differs from zero significantly only in S-wave state characterized by principal quantum number n. Since the exponent function goes exponentially to zero in the p-adic length scale associated with 2.9 TeV mass, which is roughly by a factor 32 times shorter than intermediate boson mass scale, hydrogen atom wave function is constant in excellent approximation in the effective integration volume. This gives for the energy shift

    Δ E= g2| Ψ(0)|2 × I ,

    Ψ(0) 2 =[22/n2]×(1/a03) ,

    a0= 1/(mα) ,

    I= ∫ (e-Mr/r) r2drdΩ =4π/3M2.

    For the energy shift and its ratio to ground state energy

    En= α2/2n2× m

    one obtains the expression

    Δ En= 64π2 α/n2 α3 (m/M)2 × m ,

    Δ En/En= (27/3) π2α2 k2(m/M)2 .

    For k=1 and M=2.9 one has Δ En/En ≈ 3× 10-11 for muon.

Consider next Lamb shift.

  1. Lamb shift as difference of energies between S and P wave states (see this) is approximately given by

    Δn (Lamb)/En= 13α3/2n .

    For n=2 this gives Δ2 (Lamb)/E2= 4.9× 10-7.

  2. The parameterization for the Lamb shift reads as

    Δ E(rp) =a - brp2 +crp3 = 209.968(5) - 5.2248 × r2p + 0.0347 × r3p meV ,

    where the charge radius rp=.8750 is expressed in femtometers and energy in meVs.

  3. The reduction of rp by 3.3 per cent allows to estimate the reduction of Lamb shift (attractive additional potential reduces it). The relative change of the Lamb shift is

    x=[Δ E(rp))-Δ E(rp(exp))]/Δ E(rp)

    = [- 5.2248 × (r2p- r2p(exp)) + 0.0347 × ( r3p-r3p(exp))]/[209.968(5) - 5.2248 × r2p + 0.0347 × r3p(th)] .

    The estimate gives x= 1.2× 10-3.

This value can be compared with the prediction. For n=2 ratio of Δ En/Δ En(Lamb) is

x=Δ En/Δ En (Lamb)= k2 × [29π2/3×13α] × (m/M)2 .

For M=2.9 TeV the numerical estimate gives x≈ (1/3)×k2 × 10-4. The value of x deduced from experimental data is x≈ 1.2× 10-3. There is discrepancy of one order of magnitude. For k≈ 5 a correct order of magnitude is obtained. There are thus good hopes that the model works.

The contribution of Z01 exchange is neglected in the above estimate. Is it present and can it explain the discrepancy?

  1. In the case of deuterium the weak isospins of proton and deuterium are opposite so that their contributions to the Z01 vector potential cancel. If Z01 contribution for proton can be neglected, one has Δ rp=Δ rd.

    One however has Δ rp≈ 2.75 Δ rd. Hence Z01 contribution to Δ rp should satisfy Δ rp(Z01)≈ 1.75×Δ rp1). This requires αZ,11, which is true also for the ordinary gauge bosons. The weak isospins of electron and proton are opposite so that the atom is weak isospin singlet in Abelian sense, and one has I3pI3μ= -1/4 and attractive interaction. The condition relating rp and rZ suggests

    αZ,11≈ 286=4+13 .

    In standard model one has αZ/α= 1/[sin2W)cos2W)] =5.6 for sin2W)=.23 . One has upper bound αZ,11 ≥ 4 saturated for sin2W,1) =1/2. Weinberg angle can be expressed as

    sin2W,1)= (1/2)[1 - (1-4( α1Z,1)1/2] .

    αZ,11≈ 28/6 gives sin2W,1) = (1/2)[1 -(1/7)1/2] ≈ .31.

    The contribution to the axial part of the potential depending on spin need not cancel and could give a spin dependent contribution for both proton and deuteron.

  2. If the scale of α1 and αZ,1 is that of αs and if the factor 2.75 emerges in the proposed manner, one has k2≈ 2.75× 10= 27.5 rather near to the rough estimate k2≈ 27 from data for proton.

    Note however than there are mixing angles involved corresponding to the diagonal hermitian family charge matrix Q= (a,b,c) satisfying a2+b2+c2=1 and the condition a+b+c=0 expressing the orthogonality with the electromagnetic charge matrix (1,1,1)/31/2 expressing electroweak universality for ordinary electroweak bosons. For instance, one could have (a,b,c)= (0,1,-1)/21/2 for the second generation and (a,b,c)= (2,-1,-1)/61/2 for the third generation. In this case the above estimate would would be scaled down: α1→ 2α1/3≈ 1/20.5.

To sum up, the proposed model is successful at quantitative level allowing to understand the different changes for charge radius for proton and deuteron and estimate the values of electroweak couplings of the second generation of weak bosons apart from the uncertainty due to the family charge matrix. Muon's magnetic moment anomaly and decays of neutral B allow to test the model and perhaps fix the remaining two mixing angles.

See the article Could second generation of weak bosons explain the reduction of proton charge radius?

For background see the chapters New Physics Predicted by TGD: Part I and New Physics Predicted by TGD: Part II.

See the chapter New Particle Physics Predicted by TGD.