Are OPERA neutrinos faster than light because of noninertial reference frames?
INAF – Astronomical Observatory of Padova, Italy
email: claudio.germana@gmail.com
Received:
26
December
2011
Accepted:
15
January
2012
Context. Recent results from the OPERA experiment reported a neutrino beam traveling faster than light. The challenging experiment measured the neutrino time of flight (TOF) over a baseline from the CERN to the Gran Sasso site, concluding that the neutrino beam arrives ~60 ns earlier than a light ray would do. Because the result, if confirmed, has an enormous impact on science, it might be worth doublechecking the time definitions with respect to the noninertial system in which the neutrino travel time was measured. An observer with a clock measuring the proper time τ free of noninertial effects is the one located at the solar system barycenter (SSB).
Aims. Potential problems in the OPERA data analysis connected with the definition of the reference frame and time synchronization are emphasized. We aim to investigate the synchronization of noninertial clocks on Earth by relating this time to the proper time of an inertial observer at SSB.
Methods. The Tempo2 software was used to timestamp events observed on the geoid with respect to the SSB inertial observer time.
Results. Neutrino results from OPERA might carry the fingerprint of noninertial effects because they are timed by terrestrial clocks. The CERNGran Sasso clock synchronization is accomplished by applying corrections that depend on special and general relativistic time dilation effects at the clocks, depending on the position of the clocks in the solar system gravitational well. As a consequence, TOF distributions are centered on values shorter by tens of nanoseconds than expected, integrating over a period from April to December, longer if otherwise. It is worth remarking that the OPERA runs have always been carried out from April/May to November.
Conclusions. If the analysis by Tempo2 holds for the OPERA experiment, the excellent measurement by the OPERA collaboration will turn into a proof of the general relativity theory in a weak field approximation. The analysis presented here is falsifiable because it predicts that performing the experiment from January to March/April, the neutrino beam will be detected to arrive ~50 ns later than light.
Key words: instrumentation: detectors / neutrinos / methods: numerical / reference systems / time / gravitation
© ESO, 2012
1. Introduction
The OPERA experiment (The OPERA collaboration 2011) recently reported a neutrino beam traveling faster than light. The experiment measures the distributions of neutrino time emission/detection over a baseline from the CERN to the Gran Sasso (CNGS) site. Data are collected within runs lasting for several months^{1}. In these data, the neutrino beam time of flight (TOF) turns out to be ~60 ns shorter than that calculated by taking the speed of light in vacuum. If confirmed, the result would have fundamental implications for modern physics. Basically, it reviews the fundamental postulate by Albert Einstein: “Nothing can travel faster than light” (Einstein 1905).
Although the implications are intriguing, one should consider possible systematic errors that are as yet undisclosed. The OPERA team pointed out that the research for possible sources of systematic errors is still in progress.
The scientific community has produced a plethora of works to interpret the puzzling neutrino speed result. These works were based on some systematic errors caused by the Earth’s motion (Monderen 2011) or errors in synchronization of clocks by either a third clock moving into the gravitational field of Earth (Contaldi 2011) or GPS satellites (van Elburg 2011), which were later reviewed (Besida 2011, Armando V.D.B. Assis 2011^{2}). However, several works were based on superluminal neutrino and its implications, therefore considering the OPERA neutrino speed result free of systematic errors (e.g. Tamburini & Laveder 2011; Laveder & Tamburini 2011).
In this letter we focus on the question with respect to what a measure is made. The CNGS baseline has been measured through GPS benchmarks to better than 20 cm in the European Terrestrial Reference Frame (ETRF; Colosimo et al. 2011^{3}). However, from the point of view of a light ray traveling over CNGS (or a beam of particles), the TOF might need to be corrected for effects caused by the noninertial observer ETRF, because the Earth rotates, orbits around the Sun and exhibits other motions (precession, nutation, polar motion). Nanosecondprecision measurements of OPERA are demanded to be in pulsar timing. To study the actual pulsar clock behavior, the time tagged by clocks running on Earth has to be corrected to that by an inertial observer.
A reference frame that approximates an inertial one is the Barycentric Celestial Reference System (BCRS), which is nonrotating and located at the solar system barycenter (SSB). To turn the problem into the inertial frame BCRS, we use the Tempo2 software (Hobbs et al. 2006; Edwards et al. 2006), a software conceived to correct for noninertial effects the time of arrival of photons coming from pulsars. This uses the most accurate planetary ephemerides to date (Standish 1998^{4}), in agreement with the IAU 2000 resolutions.
2. From ETRF to BCRS: Tempo2
Tempo2 is the software that bestmodels light rays traveling from a pulsar to the observatory, with an accuracy of ~1 ns. The software corrects the time of arrival of photons, as measured at the observatory, for general relativistic effects on both the photonpath and clocks running in a gravitational field, against which time is measured. In this framework Tempo2 is used to study how the inertial observer at SSB would see a signal traveling over CNGS.
After transforming the ETRF coordinates of the sites into the International Terrestrial Reference System ITRS^{5}, they are inserted into Tempo2. The software calculates the ephemerides of the vector components r(Σ,t_{UTC}), pointing from the SSB to the site Σ on the geoid at terrestrial UTC time t_{UTC}. The software prints the components of the vector r′(Γ,t_{UTC}) from the SSB to the geocenter Γ and those of the vector s(Σ,t_{UTC}) from the geocenter to the site. The vector r′(Γ,t_{UTC}) is constructed in the barycentric frame BCRS, while s(Σ,t) refers to the Geocentric Celestial Reference System (GCRS), which is nonrotating and located at geocenter. However, s(Σ,t_{UTC}) would have been equivalent if it had been constructed in the frame BCRS (Edwards et al. 2006). In s(Σ,t_{UTC}) the software accounts for Earth rotation, precessionnutation, polar motion and irregularities in both polar motion and rotation.
The vector r(Σ,t_{UTC}) in BCRS reads (1)where (i = x,y,z) and u_{x}, u_{y}, u_{z} are unit vectors in the barycentric frame BCRS.
If the light ray is emitted at the CERN at the spacetime coordinate r(CN,t_{UTC}) and is detected at Gran Sasso at , the two events are linked by the equation (2)where c is the speed of light, τ(Σ,t_{UTC}) is the barycentric coordinate time (TCB) as a function of the site Σ and UTC time. The TCB is the proper time experienced by the inertial observer at the SSB (Standish 1998). The TCB time is linked to UTC time by the relation (3)where ls are leap seconds to tie the UTC time to terrestrial time (TT) (~66 s; Seidelmann & Fukushima 1992), Δ_{E}(Σ,t_{UTC}) is the Einstein delay, which accounts for relativistic effects on clocks running in a gravitational field.
Expressing Eq. (2) in UTC we obtain (4)from which the TOF () of the signal in UTC proper time can be deduced.
The term is the difference of the Einstein delay at the two clocks and at , t_{UTC}. In the OPERA experiment clocks were not corrected for Δ_{E}, therefore they tick at UTC time, which is not a proper time.
3. Results
The geocentric coordinates of both the CERN and the Gran Sasso sites reported in Colosimo et al. (2011) were inserted into Tempo2. These coordinates refer to the origin of the OPERA detector reference frame at Gran Sasso and to the target focal point at the CERN, 730534.61 m away. An additional baseline is that between the beam current transformer (BCT) and the focal point at the CERN (743.391 m; The OPERA collaboration 2011). Because the geocentric coordinates of the BCT were not found in the literature, we assume the CNGS baseline to be 730534.61 m. An additional baseline of 743.391 m is equal to ~20 arcsec on the Earth’s surface and would induce a negligible timing correction in the relativistic term v·s/c^{2} (see Sect. 3.1).
Fig. 1
Top: timing residuals after subtracting the expected TOF from that of Eq. (4), neglecting the Einstein delay term. The simulation is performed for the year 2010, expressed in Modified Julian Day (MJD = JD2 400 000.5). The bin time of the temporal series is 10 min. Bottom: zoomin over five days. 

Open with DEXTER 
3.1. Tying noninertial clocks on Earth to the clock at BCRS
We now study the TOF from Eq. (4) neglecting the term , like in OPERA. Figure 1 (top) shows the TOF of the signal after subtracting the value one expects, 730534.61 m/c = 0.002436801 s. There is an excess of up to ~240 ns modulated by both an annual component and the Earth’s rotation (Fig. 1 bottom). Figure 1 shows that the neutrino TOF is never constant and depends on the epoch at which the measurement is performed.
To understand what happens, we need to analyze the Einstein delay difference . We print the Einstein delay for the clock at the CERN to tie UTC to TCB time (Eq. (3)). The amount of the correction at MJD = 55 197.0 (January 1, 2010) is up to ~16 s and grows with a slope of ~1.5 × 10^{8} s/s since MJD = 43 144.0003725. The linear drift of ~1.5 × 10^{8} s/s takes into account the linear term of the Einstein timedilation integral, and another term accounting for the gravitational redshift due to Earth’s potential (Irwin & Fukushima 1999). The Einstein integral accounts for the special relativistic timedilation (v^{2}/c^{2}, v velocity of the geocenter) and the gravitational redshift at geocenter (U/c^{2}, U gravitational potential where the geocenter moves). Because v ~ 30 km s^{1}, the special timedilation is v^{2}/c^{2} ~ 10^{8} s/s.
Figure 2 shows the behavior of the Einstein delay correction after fitting a firstorder polynomial over five days and subtracting it^{6}. The overall linear behavior hides a diurnal modulation of some microseconds owing to Earth’s rotation. The modulation originates from the special relativistic timedilation on the geoid with respect to the geocenter. The correction reads v·s/c^{2}, with v the velocity vector of the geocenter and s the site position vector pointing from the geocenter to the site. Because v ~ 30 km s^{1} and s ~ 6300 km (radius of Earth), the corrections amount to up to ~2 μs modulated by Earth’s rotation.
The nonzero slope in Fig. 2 after removing the linear drift signifies that the corrections hide other terms. The Einstein delay integral, referring to the motion of the geocenter, has several periodic terms (Seidelmann & Fukushima 1992).
Fig. 2
Einstein delay correction to tie the clock at CERN from UTC to TCB time after removing the fit firstorder polynomial. 

Open with DEXTER 
The inertial observer ties the UTC time of the clock at Gran Sasso to the TCB as well. The corrections to the clock at the CERN are subtracted from the Gran Sasso clock corrections to give the Einstein delay difference . Figure 3 (top) shows the result. The behavior is equal to that seen in Fig. 1. Hence, the cause for the residual on each single TOF measurement shown in Fig. 1 is that the UTC time is not corrected for relativistic effects through the Einstein delay difference , as Eq. (4) would instead require. Correcting the UTC time for the Einstein delay implies null residuals in Fig. 1, as one would expect.
Because the UTC time in OPERA is not corrected for , the two clocks at the CERN and Gran Sasso are not properly synchronized. To keep them synchronized, the noninertial observer at ETRF has to apply the corrections in Fig. 3 (top) to the UTC time at Gran Sasso.
Fig. 3
Top: corrections that the noninertial observer ETRF should apply to keep the clock at Gran Sasso synchronized with that at the CERN. Bottom: the difference behavior of the velocity components parallel to the position vector s at the two sites for the year 2010. 

Open with DEXTER 
Like Fig. 1 (bottom), the corrections to keep clocks synchronized oscillate with a diurnal behavior. Consequently, events that are simultaneous in TCB time are simultaneous in UTC time as well only when the correction in Fig. 3 (top) is zero. This happens whenever the special relativistic correction on the geoid with respect to the geocenter, v·s/c^{2}, is equal at both sites.
To verify that the behavior in Fig. 3 (top) is indeed caused by differences in the term v·s/c^{2} between the sites, Fig. 3 (bottom) shows the difference of the dot product v·j between the CERN and Gran Sasso (j unit vector along s). The difference of the velocity component parallel to s is up to 3 km s^{1}, modulated by both diurnal and annual components. Its shape reflects that of the relativistic corrections in Fig. 3 (top). A difference in velocity of up to 3 km s^{1} corresponds to a difference in timing corrections (δv × s)/c^{2} ~ 210 ns. The annual component comes from the modulation over the seasons of the angle between v and s (see also Monderen 2011).
Fig. 4
Distribution of timing corrections. Top: distribution calculated over a typical OPERA run (April to November 2010). It would imply an overall misssynchronization of ~–80 ns. Middle: distribution over an entire year. Bottom: distribution from January 1 to March 31. It would imply an overall misssynchronization of ~+50 ns. 

Open with DEXTER 
3.2. Distribution of timing corrections
To check whether the overall effect of the corrections in Fig. 3 completely vanishes, we can study their distribution.
We divide the yaxis of Fig. 3 (top) into 10 ns long channels. For each channel we count the number of events, integrating over a typical OPERA run. Figure 4 (top) shows the distribution over an OPERA run from April 29 to November 22, 2010 (from MJD = 55 314.0 to MJD = 55 522.0). It peaks at both the −240 ns and +160 ns channel, therefore it might induce an overall error on the clock synchronization of ~–80 ns^{7}.
To investigate farther out, we study the distribution for measures taken over an entire year, from January 1 to December 31, 2010. Figure 4 (middle) shows that the distribution is symmetric with respect to channel zero, inducing a null effect on the clock synchronization.
The difference between the two distributions in Fig. 4 (topmiddle) might be explained in terms of the variation of the angle between v and s over the seasons (see also Fig. 3 bottom). The OPERA runs have always been carried out from April/May to November, i.e., they missed the period from January to March that would be needed to complete one angle cycle. Therefore, if the measures are taken only in a window of the entire Earth orbit, it might be possible that the clock synchronization is affected by this systematic error and the noninertial observer on Earth would measure a TOF shorter by ~80 ns. Because he is unaware of noninertial effects, the noninertial observer should take care to carefully synchronize the clocks. Alternatively, the measurements could be integrated over one angle cycle: In this way, from the perihelion to the aphelion and back, the effect vanishes.
If we add the ~+80 ns correction to the discrepancy claimed in the OPERA experiment (−60 ns), the result might change: we obtain a neutrino beam arriving later than light.
4. Conclusions
The analysis by Tempo2 suggests that clocks on Earth should be tied to the TCB time to ensure their synchronization. The TCB time is the proper time experienced by an inertial observer located at the SSB. The UTC time ensures synchronization only when the relativistic correction v·s/c^{2} (dubbed differential special timedilation; Edwards et al. 2006) is equal at both clocks. Differences in this correction of up to ~240 ns are seen and are modulated by both the Earth’s rotation and revolution.
The overall effect of a nonsynchronization of clocks may vanish over one Earth orbit. Integrating the nonsynchronized clock over months might induce an overall effect on the order of several tens of nanoseconds.
If the analysis described in this letter holds for the OPERA experiment, the excellent measurement by The OPERA collaboration (2011) would turn into a proof of relativistic theories (Einstein 1905, 1916). The analysis might cast doubts on interpretations dealing with superluminal neutrino that are, or are not, accounting for relativistic effects (e.g. Tamburini & Laveder 2011; Laveder & Tamburini 2011).
To check whether or not the present study holds, it might be interesting to see results from OPERA runs performed in a different period of the year than usual. Figure 4 (bottom) shows the distribution of the corrections for noninertial effects from January–March. Again, it is not symmetric with respect to channel zero, but shows an opposite result to that of the distribution on the top. From this result the neutrino beam should arrive ~50 ns later than light.
For full accuracy it is suggested to use coordinates in ITRS. ETRF coordinates were transformed into ITRS through the online converter at http://www.epncb.oma.be/_dataproducts/coord_trans/index.php
The ~–80 ns overall effect holds for the last OPERA results release as well (The OPERA collaboration 2011), because the run lasted from October 21 to November 7, 2011.
Acknowledgments
I would like to thank Andrej Čadež (Faculty of Mathematics and Physics, University of Ljubljana) and Massimo Calvani (INAF – Astronomical Observatory of Padova) for stimulating discussions. I also thank Rodolfo Angeloni (Departamento de Astronomía y Astrofísica, Pontificia Universidad Católica de Chile) for his kind support.
References
 Besida, O. 2011 [arXiv:1110.2909v1] (In the text)
 Contaldi, C. R. 2011 [arXiv:1109.6160v2] (In the text)
 Edwards, R. T., Hobbs, G. B., & Manchester, R. N. 2006, MNRAS, 372, 1549 [NASA ADS] [CrossRef] (In the text)
 Einstein, A. 1905, Ann. Phys., 322, 891 [NASA ADS] [CrossRef] (In the text)
 Einstein, A. 1916, Ann. Phys., 354, 769 [NASA ADS] [CrossRef] (In the text)
 Hobbs, G. B., Edwards, R. T., & Manchester, R. N. 2006, MNRAS, 369, 655 [NASA ADS] [CrossRef] (In the text)
 Irwin, A. W., & Fukushima, T. 1999, A&A, 348, 642 [NASA ADS] (In the text)
 Laveder, M., & Tamburini, F. 2011 [arXiv:1111.4441v1] (In the text)
 Monderen, D. 2011 [arXiv:1110.3581v1] (In the text)
 Seidelmann, P. K., & Fukushima, T. 1992, A&A, 265, 833 [NASA ADS] (In the text)
 Standish, E. M. 1998, A&A, 336, 381 [NASA ADS] (In the text)
 Tamburini, F., & Laveder, M. 2011 [arXiv:1109.5445v7] (In the text)
 The OPERA collaboration 2011 [arXiv:1109.4897v2] (In the text)
 van Elburg, R. A. J. 2011 [arXiv:1110.2685v4] (In the text)
All Figures
Fig. 1
Top: timing residuals after subtracting the expected TOF from that of Eq. (4), neglecting the Einstein delay term. The simulation is performed for the year 2010, expressed in Modified Julian Day (MJD = JD2 400 000.5). The bin time of the temporal series is 10 min. Bottom: zoomin over five days. 

Open with DEXTER  
In the text 
Fig. 2
Einstein delay correction to tie the clock at CERN from UTC to TCB time after removing the fit firstorder polynomial. 

Open with DEXTER  
In the text 
Fig. 3
Top: corrections that the noninertial observer ETRF should apply to keep the clock at Gran Sasso synchronized with that at the CERN. Bottom: the difference behavior of the velocity components parallel to the position vector s at the two sites for the year 2010. 

Open with DEXTER  
In the text 
Fig. 4
Distribution of timing corrections. Top: distribution calculated over a typical OPERA run (April to November 2010). It would imply an overall misssynchronization of ~–80 ns. Middle: distribution over an entire year. Bottom: distribution from January 1 to March 31. It would imply an overall misssynchronization of ~+50 ns. 

Open with DEXTER  
In the text 