Cosmology and the massive photon frequency shift in the Standard-Model Extension

The total red shift $z$ might be recast as a combination of the expansion red shift and a static shift due to the energy-momentum tensor non-conservation of a photon propagating through Electro-Magnetic (EM) fields. If massive, the photon may be described by the de Broglie-Proca (dBP) theory which satisfies the Lorentz(-Poincar\'e) Symmetry (LoSy) but not gauge-invariance. The latter is regained in the Standard-Model Extension (SME), associated with LoSy Violation (LSV) that naturally dresses photons of a mass. The non-conservation stems from the vacuum expectation value of the vector and tensor LSV fields. The final colour (red or blue) and size of the static shift depend on the orientations and strength of the LSV and EM multiple fields encountered along the path of the photon. Turning to cosmology, for a zero $\Omega_{\Lambda}$ energy density, the discrepancy between luminosity and red shift distances of SNeIa disappears thanks to the recasting of $z$. Massive photons induce an effective dark energy acting `optically' but not dynamically.


I. INTRODUCTION
Astronomy is almost entirely built up from information coming from Electro-Magnetic (EM) signals, interpreted with Maxwellian linear and massless electromagnetism. This latter is possibly an approximation of a broader theory, as Newtonian gravitation is of the Einsteinian one. Observations [1][2][3] have led to the proposal that the universe contains up to 96% dark matter and dark energy (both entities thus far remaining theoretically unexplained and experimentally undetected) and holds to general relativity, as the correct theory of gravitation. Others, unconvinced by these ad hoc ingredients, propose new theories of gravitation, though general rel-ativity scores high marks in all tests so far.
Faced with this dichotomy and the respective pitfalls, we turn to the signals. Could a different interpretation of light lead to a third option? After all, modern physics -relativity and quantum mechanics -started with a new interpretation of light. Ultimately, physics is in any case obliged to examine its foundations through either the formulation of new types of matter and energy or embracing new conceptions of gravitation or electromagnetism.
The photon is the only massless free particle of the successful Standard-Model (SM), lately challenged by the neutrino mass, the light mass of the Higgs boson and the absence of candidate particles for the dark universe. So let us focus on massive photons, leaving non-linear electromagnetism for future work.
In 1922 de Broglie proposed a massive photon [9] and, through the group velocity dispersion of the latter, estimated the mass upper limit as 10 −53 kg [10,11]. The origin of his conception of the photon mass resides in his the-ory of fusion, which couples two or more free Dirac fields to produce a scalar or vector Klein-Gordon equation, typically an electron-positron or neutrino-anti-neutrino duet [12][13][14][15][16][17][18]. In [18] the modified Maxwell equations were first written in a non-covariant form. Thus, de Broglie, and accidentally his disciple Proca [19][20][21][22], laid down the first massive electromagnetism, compliant with the Lorentz(-Poincaré) Symmetry (LoSy), though not gauge-invariant. Contrary to common belief, the Proca Lagrangian describes generic composite particles, among which Proca states that the photons are massless 'doublets de masse nulle' [22]. For a comprehensive monograph see [23].
Energy-momentum non-conservation already manifests itself in the Maxwellian theory for a photon crossing space-time dependent electromagnetic (EM) background fields (host galaxy, intergalactic and Milky Way). Instead, we show here that the de Broglie-Proca (dBP) massive photon undergoes energy-momentum non-conservation also in the case of constant EM fields, since it couples to the non-constant EM background potential.
Motivated by the above challenges, the SM Extension (SME), based on the LoSy Violation (LSV), was put forward [24,25]. The LSV breaks the invariance of the laws of physics for all positions and inertial velocities of the observers, and thus naturally accommodates anisotropies. The SM is LoSy invariant up to the energy scales at which the LSV occurred in the early universe. The SME allows the testing of the low-energy manifestations of the LSV in the present universe.
The foundations of the analysis of light propagation in the SME were laid down in [26,27] . Going beyond the SM, a massive photon emerges by scrutiny of the group velocity dispersion or by recasting the SME Lagrangian. The effective, and frame-dependent, mass is proportional to the value of the LSV parameters and, in contrast to the dBP formalism, the SME massive photon is gaugeinvariant [26,27] and is compatible with the LSV upper limits [28,29]. Some of the benefits and drawbacks of attributing a mass to photons have been analysed recently [30][31][32][33][34].
Moreover, in specific conditions sub-and super-luminal velocities, imaginary and complex frequencies, and birefringence appear. Evidently in any condition, LSV anisotropy and inhomogeneity are present.
Also in the SME, the photon energy-momentum tensor is not conserved in vacuum, and most remarkably, includes terms depending on constant LSV and EM fields [27]. Furthermore, we have determined a non-expansion related frequency shift (towards the red or the blue) z LSV [35].
Quantified by the Hubble(-Humason)-Lemaître constant H 0 and detected at the end of the '20s, the cosmological expansion stretches wavelengths from astrophysical sources, including Type Ia supernovae (SNeIa). But since the end of the '90s, these latter, considering their faded luminosity with respect to the expected energy emitted, prompted the conception of dark energy, a mysterious fluid acting as the source of the cosmic accelerated expansion [36,37].
If the cosmic flow expanded at a constant rate, an SNIa red shift would be directly proportional to its distance, and thus to its brightness. Instead, the accelerating universe filled with dark energy implies that space expanded less quickly in the past than it does now. The light from distant objects such as SNeIa is less stretched during its journey to us, given how slowly space expanded during much of the time. An SNIa at a given distance (computed through its brightness, since SNeIa are believed to be standard candles) appears less red shifted than it would in a universe without dark energy. To put it another way, the red shift is not a measurement of distance but of expansion. When a distant source is observed, the lower expansion induces a lower red shift.
The apparent acceleration of the younger Universe could be simply accounted for by a non-zero cosmological constant (originally proposed by Einstein to refute an expanding universe) on the left (geometrical) side of the field equations, or else by an extra term on the right (energy-matter) side acting as a source term for gravity, and corresponding to a vacuum energy density.
Somewhat similarly, in this paper we also refer to an extension, namely that of the SM. Indeed, we identify the vacuum energy density in the framework of the SME as being due to the LSV field, which manifests itself through the z LSV shift.
The above-mentioned third option may be checked by adding algebraically the z LSV shift to the cosmological expansion shift z C in the analysis of the data from SNeIa. We go through Baryonic Acoustic Oscillation (BAO), Cosmic Microwave Background (CMB) and gravitational lensing, often reputed to back the existence of dark energy, to see if any counter-evidence would threaten our hypothesis. We also deal with time dilation in SNeIa.
Attributing different origins to the nature of dark energy does not fully describe the scope of the debate that arose since the appearance of [36,37]. One school of criticism focuses on the reliability of SNeIa, BAO, CMB and gravitational lensing data, and challenges the counterarguments defending the evidence for dark energy based on data, e.g., [44][45][46][47][48].
Another school of thought addresses alternatives to dark energy, of which we mention a few: photons oscillating into axions to produce a dimming effect for SNeIa [49,50], the different spectrum of primordial fluctuations generated by inflation [51], timescape cosmology and the inhomogeneous Universe [52], SNIa luminosity depending on red shift [53,54], stronger evidence of acceleration dipole approximately aligned with the CMB dipole than evidence of a monopole dark energy acceleration [55,56].
Last but not least, the identification of SNeIa as standard candles is at risk of being breached, e.g., [57].
We wish to emphasise that we do not take a position on data reliability. Our main point is to recast the red shift, nowadays uniquely explained by expansion; our recasting may be tailored to whatever set of data on which consensus will be reached. Moreover, a lesser z LSV shift value could coexist with dark energy. Generally speaking, the Concordance Lambda-Cold Dark Matter (ΛCDM) model may not rest on very firm ground, if Planck satellite CMB data, commonly reputedly its pillar, are reread and significantly support a closed universe [58].
In conclusion, our work is part of a large and scientifically healthy afflatus that revisits modern cosmology to verify its soundness.
Some wording on the relations between LSV, SME, extensions of the SM and photon mass are in order. LSV is not a sufficient condition to induce an effective mass. Generally speaking, extension of the SM such as the Minimally Supersymmetrised SM (MSSM), the Massive Neutrino Model (MNM), the Grand Unification Theory (GUT) do not induce an effective mass either. Thereby, it is the SME which is based on LSV that may cause an effective mass. Let us be more specific.
The SME-LSV factors are represented by a k AF α [metre −1 ] 4-vector when the handedness of the Charge conjugation-Parity-Time reversal (CPT) symmetry is odd and by a k ανρσ F [dimensionless] tensor when even. The k AF α vector coming from the Carroll-Field-Jackiw Lagrangian [59] induces always a mass, while the k ανρσ F tensor only in a supersymmetrised context after photino integration [26,27].
The SME is built up by means of the Effective Field Theory approach and has the SM and General Relativity (GR) as possible limiting cases. In this paper, we work within a scenario such that the LSV does not leave its imprints in the GR sector. We instead focus on a special situation where the LSV is present in the photonic sector, through the k AF α and the k ανρσ F terms. The spacetime metric, spin connection and curvature are unaffected by the LSV; we actually maintain Minkowski space-time with the anisotropies parametrised by k AF and k F . This means that the effects of the LSV on the red shift stem exclusively from the presence of the violating parameters in the photonic sector. This is the viewpoint we adopt here. It is, however, a relevant matter, for future work to consider the presence of LSV in the gravitational sector [60][61][62][63][64][65] of SME and thereby reassess its effect on the red shift.
The rest of this paper is structured as follows. Section II shows the energy-momentum non-conservation of the dBP photon and prepares for the more difficult analysis in the SME dealt with in Sect. III, devoted to light propagation in this context; Sect. IV focuses on the reinterpretation of dark energy through the z LSV shift; Sections V and VI discuss the results and their implications.

II. ENERGY NON-CONSERVATION IN MAXWELL AND DE BROGLIE-PROCA THEORIES
We recall here a basic feature of the dBP theory concerning the non-conservation of the photon energymomentum and pave the way to the analysis in the SME. We imagine the photon crossing an electromagnetic background.
The dBP equations of motion derived from the Lagrangian [19] correspond to the divergence of the electric field and the curl of the magnetic field. They are given by is the speed of light and = 1.055 × 10 −34 [kg m 2 s −1 ]. T stands for the total EM quantities due to background and photon contributions; j β is an external current, if it exists, and µ 0 = 1.257 is the magnetic permeability. Splitting the EM tensor field and the EM 4-potential in the background (capital letters) and photon (small letters) contributions, we have gives Equation (3) tells us that the dBP photon interacts with the background through the potential even when the background field is constant. Indeed, if a field is constant, its associated potential is not Conversely, this is not the case for the Maxwell photon, Eq. (5), which interacts only with a non-constant field The energy-momentum density tensor θ α τ [Jm −3 ] for the dBP photon is obtained after a lengthy computation and is given by For E, B, the electric and magnetic fields of the background, e, b, the electric and magnetic fields of the photon, Φ, A, scalar and vector potentials of the background, and φ, a, scalar and vector potentials of the photon, we make explicit first the energy density [Jm −3 ] is the electric permittivity, and second we make explicit the generalised Poynting density vector [Jm −3 ] The energy-momentum density tensor variation ∂ α θ α τ [Jm −4 ] is given by (9) which, after multiplication by c, in the explicit form be- In conclusion, the energy-momentum density tensor of the dBP photon is not conserved. In addition to the Maxwellian terms, the mass couples with the background potential time-derivative.

A. Non-conservation
In contrast to the LSV vector, the LSV tensor does not violate CPT. The frequency LSV shift that we shall be dealing with here is an observable of CPT violation, since it depends on the vector and tensor formulations. Incidentally, we shall see that the leading term is k AF 0 which violates CPT.
Indicating by the symbol * the dual field, the photon energy-momentum density tensor θ α τ [Jm −3 ] is [27,35] We render explicit the energy density [Jm −3 ] where k AF 0 represents the time component of the breaking vector and the k ανκλ F tensor is decomposed as χ ij and ζ m,n are symmetric in the indexes, ξ il does not have symmetry properties; χ ij and ζ m,n have 6 components each, while ξ il has 9, totalling the 21 components of k ανκλ F . We now render explicit the generalised Poynting density vector [Jm −3 ] The energy-momentum density tensor variation ∂ α θ α τ [Jm −4 ] is given by After multiplication by c, the energy-momentum density tensor variation in the explicit form becomes The energy-momentum density tensor variation is due to the following contributions: • Maxwellian, LSV-independent terms that we have seen in Sect. II.
There are three massive contributions (though not all the components are mass dependent).
• EM background independent terms implying that the energy-momentum density flux is not conserved in the absence of EM fields, if the LSV fields are space-time dependent. This is really a distinctive feature of the SME.
• An LSV and EM space-time dependent term.
• A constant term coming solely from the CPTodd handedness represented by the Carroll-Field-Jackiw (CFJ) electrodynamics [59]. Its action entails a non-constant 4-potential, for a constant EM background and a constant k AF . Indeed, there is an explicit x α coordinate dependence at the level of the Lagrangian, exactly as in the dBP theory.
The breaking tensor k F appears either under a derivative or coupled to a derivative of the EM background.  [29]. 2 Astrophysical observations [28]. Such estimates are close to the Heisenberg limit on the smallest measurable energy or mass or length for a given time, t, set equal to the age of the universe. A full table of values was shown in [35].
The above contributions determine the energymomentum non-conservation, according to the Noether theorem (symmetry breaking). Put another way, the photon exchanges energy with the LSV and EM fields.
Comparing the non-conservation in dBP and SME cases, the former is relativistic and respects LoSy, while the latter is gauge-invariant. The advantage of the latter lies in tracing the origin of the mass in the LSV vacuum energy. In both cases a modification of the SM is necessary.

Origin of the non-conservation
An estimate of the energy change that light would undergo was given [35]. The wave-particle correspondence, even for a single photon [66], leads us to consider that the light-wave energy non-conservation is translated into photon energy variation and thereby into a red or a blue shift. The energy variations, if there are losses, would translate into frequency damping.
According to [67], the LoSy breaking 4-vector, k AF , and the rank-4 tensor, k F , correspond, respectively, to the vacuum condensation of a vector and a tensor field in the context of string models. They describe part of the vacuum structure, which appears in the form of spacetime anisotropies. Therefore, their presence on the righthand side of Eq. (17) reveals that vacuum effects are responsible for the energy variation of light waves, which, in turn, correspond to a photon frequency shift. In plain terms, vacuum anisotropies are the really responsible for changing the frequency of light emitted by astrophysical structures.

B. Sizing the LSV frequency shift
In Eq. (17), we may neglect the tensorial perturbation in the frame of Super-Symmetry (SuSy) since it is less likely to condense than the CFJ vectorial perturbation [35,68,69]. Independently of SuSy, the k F term can be neglected since it is quadratic in the field strength and in frequency. The CFJ term instead contains a single derivative and is linear in the frequency. Thus, for optical frequencies of the SNeIA, the k AF is the dominant contribution.
Excluding the space-time dependent components of the (inter-) galactic magnetic and LSV fields, only one contribution, the last term in Eqs. (17,18), survives. Finally, excluding external currents and large scale electric fields, an estimate was given along the line of sight observersource for this term [35], where * F 0i are the magnetic components of the dual EM background tensor field and f i0 are the electric components of the EM photon tensor field. There is no relation between the k AF 0 and the photon mass 1 , though there are obviously massive terms in Eq. (17). 1 For simplicity of notation we replace k AF by V ; further, the 4- [26] or (9) in [27] confirms the CFJ dispersion relation [59]; that is, For a massive photon, the rest mass is computed in the rest frame of the photon, k = 0, rendering Eq. (21) as Equation (21) has two solutions: one non-massive for ω = 0 (since gauge symmetry is not broken in CFJ formalism), the other for ω = c| V |. In conclusion, the time component of the CFJ perturbation vector does not contribute to the photon mass.
Assuming that the energy variation of a light-wave corresponds to a frequency shift for a photon, we converted the energy variation given by Eq. (19) into ∆ν, ν being the photon frequency. For ν = 486 THz [35] |∆ν| 486 THz where t LB is the look-back time and ̺ is an arbitrary attenuation parameter that takes into account that the magnetic fields (host galaxy, intergalactic, Milky Way) estimated at B = 5 × 10 −10 − 5 × 10 −9 T each, passed through by the photon, probably have different orientations and thus compensate for each other. We have not considered any relevant magnetic field at the source. For a numerical estimate of Eq. (22), let us consider a source at z = 0.5, t LB = 1.57 × 10 17 s, f i0 = 3.79 × 10 −9 V s m −2 , an average B = 2.75 × 10 −9 T. We wish to determine what the maximum level of z LSV could be.
This leads to a value of the order of 10% of the total z, for ̺ ≈ 8 × 10 −26 and ̺ ≈ 8 × 10 −8 , respectively. These estimates are arbitrary and just show that they can recover the largest frequency shifts in agreement with observational data from SNeIa. Nevertheless, with different assumptions on ̺ and k AF 0 , we can get considerably smaller percentages of the rate z LSV /z that are still in agreement with data.
Indeed, it is important to note that the single z LSV shift from a single SNIa may be small or large, red or blue, depending on the amplitude and orientations of the LSV (vector or tensor) and of the EM fields (host galaxy, intergalactic medium, Milky Way), and obviously the distance of the source. In any case, the final z LSV is the result of accumulated shifts, both red and blue, encountered along the path.
We have set upper limits starting from the physics; in the following, we shall set upper limits from the cosmological data in answer to the question of what range of values of z LSV can be accommodated by cosmology. In following Sections, we shall refer to a generic z LSV given by all the terms in Eq. (17).

C. Recasting z
The existence of a photon frequency shift not due to the relative motion of the source and the observer belongs to the realm of physics. The role of such shifts in cosmology can attain three levels: replace totally or accompany to a certain degree the expansion shift. In the case of total replacement, we would return to the conception of a static universe, which still passes some tests but fails many others [2,3]. At the other extreme, such a shift could be marginal for cosmology, but nevertheless be of relevance for fundamental physics.
Here, we explore an intermediate option for which the static shift is superposed on the expansion shift.
After recalling that the definition of z = ∆ν/ν o , where ∆ν = ν e − ν o is the difference between the observed ν o and emitted ν e frequencies, or else z = ∆λ/λ e for the wavelengths, we pose the following conjecture: expansion causes λ e to stretch to λ c ; that is, The wavelength λ c could be further stretched or shrunk for the LSV shift to where z o is the spectroscopically or photometrically observed z. The second order is non-negligible for larger z C . We model the behaviour of z LSV with distance in three 2 different ways, Table II, according to whether the frequency variation is proportional to • the instantaneous frequency and to the distance, • the emitted frequency and the distance, or • just the distance.
The z C shift stems from the expansion of the universe, whereas the z LSV would occur also along a static distance r. The k i parameters can take either positive (frequency increase) or negative (frequency decrease) values.

IV. THE LSV FREQUENCY SHIFT AND DARK ENERGY
A. Supernovae: luminosity and red shift distances A greater than expected SNIa luminosity distance d L for a given red shift led to the proposal of dark energy in order to reach consistency with the data [36,37]. We pursue the same consistency by using Eq. (23) instead.
We do not intend in this paper to propose a fully fledged alternative cosmology, but solely to take the first steps towards the third option presented in Section I. In the following we do not intend to state strict conclusions for cosmology; instead, we limit ourselves to exploring whether existing data can accommodate our reinterpretation; namely, the recasting of z through the additional non-expansion related shift. For this exploration, we pick a popular cosmology simulator [70,71] and draw from it our initial considerations.
The values that z C should assume for a fixed d L are computed for two values (0.72 and 0) of Ω Λ energy density and three values (0.7, 0.67, 0.74) of h in Table III.
We find that a small percentage correction of the red shift (from 0.01% up to 10%, a range we have proven feasible in Sect. III) allows us to recover the same luminosity distances. The explored range of z values for SNeIa is (0-2). The most distant SNIa is at z = 1.914 [72]. Inverting Eq. (23), we get If the z LSV shift is blue, and thus negative, the photon gains energy, which implies that the cosmological z C is higher than the observed z. If red, and therefore positive, the photon dissipates energy along its path, which implies that the cosmological z C is smaller than the observed z.
From a first glance at Table III, it appears that in most cases adding a negative blue z LSV shift to the cosmological z C is sufficient for agreement between d L and z. Instead, the red shifts appear for h = 0.67 and small values of z. It is also evident that for higher h values a larger z recasting is necessary, and that the shift is bluer. Patently, these are macroscopic and averaged indications. The next step necessitates picking each SNIa from a catalogue and evaluating the associated z LSV .
Preliminary simulations show that z LSV turns into red or blue shifts, depending also on other cosmological parameters, e.g., the value of matter density Ω m or of the curvature Ω k , and on the behaviour of z LSV with fre-quency and distance, Table II.
At the bottom of Table III, we have added two lines for z = 4 and z = 11 to gauge the behaviour of z LSV at large z: the most distant superluminous SN is at z = 3.8893 [73], and the most distant galaxy is at z = 11.09 [74].
The data in Table III can be commented on by momentarily adopting the usual terminology of accelerated expansion. The acceleration is not observable in our immediate neighbourhood; it would start to be noticeable further out and increase up to a maximum at z ≈ 0.5. After this threshold value, the acceleration would decrease and for approximately z ≥ 4, the so-called turning point, it would change sign corresponding to the deceleration phase associated with the older universe.
How to reinterpret the above through a recast z? Far out, for approximately z ≥ 4, z LSV becomes definitely (and more and more) red. Since a red z LSV shift corresponds to dissipation, we feel comforted by this easy interpretation for remote distances.
Instead, the major difficulty is explaining why z LSV starts from a small (blue or red) value in our immediate neighbourhood and increases up to a (blue) maximum for z ≈ 0.5. This is maybe due to our location in a local void, or other manifestations of inhomogeneities or anisotropies. But before embarking in this sort of argument, the controversy on the data supporting dark energy [44][45][46][47][48][49][50][51][52][53][54][55][56][57]75] should be first of all settled.

Time dilation and Supernovae
Since the SME induces an effective mass to the photon [26,27], we are dealing with massive photons. Blondin et al. [76] compared several spectra from the more distant SNeIa with those of nearer ones and found that the more distant explosions took longer to unfold. They pointed out that time dilation in SNeIa cannot be attributed to 'tired light'. Let us discuss not only this statement in the context of this paper, but also our assumptions on reinterpreting z by examining three related issues.
Firstly, just attributing mass to a photon does not determine per se a decaying frequency, a feature of the 'tired light'. Indeed, massive photon propagation obeys the dBP equation which, in the absence of an EM background, is written as [23] We hold to the observed zo := z, Eq. (23), and show the values that the cosmological shift z C and possibly z LSV should assume for obtaining a given luminosity distance d L , in the first column. We assume different values for H 0 = h × 100 km/s per Mpc and for Ω Λ densities (the curvature and radiation densities are set to zero, Ω k = Ω rad = 0). In the second column for matter density Ωm = 0.28, energy density Ω Λ = 0.72 and h = 0.7, we pose z LSV = 0 and thereby z = z C ; in the third column for Ωm = 0.28 but Ω Λ = 0 and h = 0.7, we show the values of z C that determine the same d L ; in the fourth, eighth and twelfth columns, the percentage variation of z C ; in the fifth, ninth and thirteenth columns, from Eq. (23), z LSV = z − z C 1 + z C ; in the sixth, tenth and fourteenth columns, the rate z LSV z ; in the seventh column for Ωm = 0.28 but Ω Λ = 0 and h = 0.67, the values of z C which determine the same d L ; in the eleventh column for Ωm = 0.28 but Ω Λ = 0 and h = 0.74, the values of z C which determine the same d L . Red or blue shifts correspond to positive and negative values of z LSV , respectively. The most distant SNIa is at z = 1.914 [72], the most distant superluminous SN is at z = 3.8893 [73], and the most distant galaxy is at z = 11.09 [74]. The numerical values are derived from a Cosmology Simulator [70,71].  [77,78]. For violating photon energy conservation, apart an external current term, there must be an EM background field, Eq. (9). Secondly, the energy variation computed here occurs regardlessly of the expansion. The associated shift has been named z LSV and evaluated to be at most 10% of the total z. It is legitimate to ask whether there is any impact on time dilation due to this reinterpretation of z.
In an expanding universe, the ratio of the observed frequency to that emitted by a distant object or the observed rate of any time variation in the intensity of the emitted radiation will be proportional to Blondin et al. [76] estimated (1 + z) b as b = 0.97 ± 0.10, b = 1 being the value predicted by the Friedmann-Lemaître-Robertson-Walker cosmology. A variation of 10% for b is tantamount to a variation of the same order for z, as the Taylor expansion shows Therefore, supposing the total z as the sum of an expansion related z C and a static z LSV is compatible with the findings in [76].
It appears that the SNIa spectra often shift from higher, f h ≈ 8 × 10 14 Hz, to lower, f l ≈ 4 × 10 14 Hz, optical frequencies during the burst [76]. Photons emitted at the end will therefore take more time to reach the observer than those at start. This delay would occur even when the source is static, mimicking time dilation. For a source at distance d, the difference in arrival times is to be added to the burst duration. In SI units, the difference in arrival times is [82] (28) Inserting the official upper limit on photon mass 10 −54 kg, d = 1.4 × 10 9 light-years (approximately corresponding to z = 0.1), ∆t is in the order of 10 −21 s. In an expanding universe, Eq. (28) becomes [79,81] where Equation (30) reduces the outcome of Eq. (28) by a factor 10 for common values of the Ω densities. Hence, we have shown that although massive photons determine an effect similar to time dilation for not-moving sources or add a static contribution to moving sources, such an effect is utterly marginal for SNeIa.
In conclusion, our assumptions and results are compatible with current literature on SNeIa.

B. BAO, CMB and gravitational lensing
In Sect. I, we gave a brief account on the debates on i) the nature of dark energy, ii) the data reliability from SNeIa, BAO, CMB and gravitational lensing, and iii) the alternatives to dark energy. The latter two debates are not dominant as there is a prevailing tendency to converge towards the ΛCDM Concordance Model.
Nevertheless, there is an increasing interest in discussing the foundations of common assumptions. A critical analysis of the supposed confirmation of dark energy by the Baryonic Atomic Oscillations (BAO) and by the Cosmic Microwave Background (CMB) has recently appeared [53,54] in which it is shown that that a nonaccelerated universe is nicely able to fit low and high red shift data from SNeIa, BAO and the CMB on the main assumption of SN luminosity dependency on red shift. Now, Kang et al. [75] have found a significant (a 99.5%) correlation between SN luminosity and stellar population age. It is worth recalling that the CMB can be interpreted differently, possibly without dark energy [89] and in any case outside a canonical flat universe [58].
We have interpreted z LSV as due to vacuum energy, an effective dark energy of LSV origin, although not producing acceleration of the expansion. We do not dispute the results from BAO and the CMB proving the existence of dark energy per se. Our interpretation would be falsified, if there were a strict proof of acceleration (m/s 2 ) with respect to an inertial frame.
Further, it is worth considering that BAO and CMB data can accommodate our recasting of z due to the uncertainties of spectroscopic and especially photometric measurements. In particular for the CMB, the photons were all emitted at the same time and same distance from us. Thus, any change in the definition of the red shift should not bear huge consequences.
For BAO, the angular size of the sound horizon is measured for the large-scale structure at different red shifts, which allows us to measure H(z), which depends on Ω m , Ω Λ and other variables. Since d A (angular distance) is proportional to d L /(1 + z) 2 and the angular size of an object is inversely proportional to d A , it would not be surprising that the same z LSV mimicking the dark energy for SNeIa, acts similarly for the BAO peak.
Finally, we tackle gravitational lensing. The classical method used to probe dark energy is to measure the ratio of average tangential shear as a function of the red shift behind the clusters [90]. This ratio depends on the relation between distance and red shift, which is modified by dark energy, or by z LSV .

V. DISCUSSION OF THE MAIN RESULTS
We have assumed the total red shift z as a combination of the expansion red shift z C and of a static, red or blue shift z LSV (r), r being the travelled distance. The latter shift is due to the energy non-conservation of a photon propagating through EM background fields (host galaxy, intergalactic and Milky Way). Beyond the Maxwellian contributions, if the photon is massive, such propagation may be described within the framework of de Broglie-Proca (dBP) formalism or within the framework of the Lorentz(-Poincaré) Symmetry Violation (LSV) associated with the Standard-Model Extension (SME). In the latter case, the non-conservation stems from the vacuum expectation value of the vector and tensor LSV fields.
Our understanding is that z LSV is a manifestation of an effective dark energy caused by the expectation values of the vacuum under LSV. Indeed, we suggest that dark energy, i.e. LSV vacuum energy, is not causing an accelerated expansion but a frequency shift.
The single z LSV shift from a single SNIa may be small or large, red or blue, depending on the orientations of the LSV (vector or tensor) and of the EM fields (host galaxy, intergalactic medium, Milky Way), as well as the distance of the source. In any case, the colour of z LSV is the final output of a series of shifts, both red and blue, encountered along the path.
If the z LSV shift is blue, thus negative, the photon gains energy; it implies that the real red shift is larger than the measured z. If red, thus positive, z LSV corresponds to dissipation along the photon path; it implies that the real red shift is smaller than the measured z.
Recasting z, on average, we observe a blue static shift for z ≤ 2, but a red one in our local Universe for smaller values of the Hubble(-Humason)-Lemâitre parameter (67 − 74 km/s per Mpc), and always red for z > 4.
The peculiarity of our approach is that a single mechanism could explain all the positions of the SNeIa in the (µ, z) plan, µ being the distance modulus, including the outliers. The experimental and observational limits on LSV and magnetic fields are fully compatible with our findings.
In future work, leaving aside massive photons, we will show that a frequency shift is also produced by a generalised non-linear electromagnetism, encompassing the formulations of Born and Infeld, and Euler, (Kockel) and Heisenberg. The calculated limits on z LSV will be applicable to z NL for the non-linear electromagnetism.
Apart from an additional shift occurring also in Maxwellian electromagnetism in the presence of a spacetime dependent EM field, we will show that departing from Maxwellian electromagnetism in three different directions (classically massive dBP electromagnetism or non-linear electromagnetism or the SME) leads to a common conclusion: the red shift is not only due to expansion. If and by how much such a departure is relevant for cosmology is the real question.
Nineteenth century Maxwellian electromagnetism and the more modern Einsteinian gravitation have been well tested. This has not impeded the proposition of alternative formulations of electromagnetism and gravity. The lack of experimental proof on the dark universe and the successes of general relativity prompt us to revisit astrophysical observations, largely based on light signals, with non-Maxwellian electromagnetism, opening the door to radically new interpretations.
In future explorations, we will have to deal with the comparison of z LSV with the error on z. The analysis of the error on spectroscopic and photometric measurements is destined to become a pivotal issue for cosmology [96][97][98].
The z LSV shift considered here is at most 10% of the z C expansion shift relative to H 0 = 70 km/s per Mpc and is thus below 2.3 × 10 −19 ∆ν/ν per m (units for a static red shift). It would be desirable to test frequency invariance in vacuo [99][100][101][102] with a ground or space based interferometer.