text:
# Specific heat and magnetic measurements in Nd0.5Sr0.5MnO3 and R0.5Ca0.5MnO3 (R=Nd, Sm, Dy and Ho) samples

## I Introduction

The physical properties in charge ordering (CO) manganese perovskites are thought to arise from the strong competition involving a ferromagnetic double exchange interaction, an antiferromagnetic superexchange interaction, and the spin-phonon coupling,,,,,,. These interactions are determined by intrinsic parameters such as doping level, average cationic size, cationic disorder and oxygen stoichiometry. CO compounds are particularly interesting because spin, charge and orbital degrees of freedom are at play simultaneously and classical simplifications, that neglect some of these interactions, do not work. More detailed information on the physics of manganites can be found in a review paper by Myron B. Salamon and Marcelo Jaime.

L. Ghivelder et al. reported specific heat measurements in LaMnO3+δ samples and found that the specific heat at low temperature is very sensitive to small variations of $`\delta `$, similarly to results published by W. Schnelle et al. in a Nd0.67Sr0.33MnO3-δ sample. In this latter work, and also in a paper by J. E. Gordon et al., a Schottky-like anomaly was found at low temperatures. They associated this result to the magnetic ordering of Nd3+ ions and to the crystal-field splitting. F. Bartolomé et al. also found a Schottky-like anomaly in a closely related compound of NdCrO3. They proposed a crystal-field energy level scheme in agreement with neutron-scattering studies in the same sample. In two papers V. N. Smolyaninova et al., studied the low temperature specific heat in Pr1-xCaxMnO3 (0.3$`<`$x$`<`$0.5) and La1-xCaxMnO3 (x=0.47, 0.5 and 0.53). They fitted the data with an excess specific heat, C´(T), of non-magnetic origin associated with charge ordering. They also showed that a magnetic field, sufficiently high to induce a transition from the charge ordered state to the ferromagnetic metallic state, did not completely remove C´(T). However, no Schottky anomaly was found in any of these compounds.

We have already presented specific heat measurements with applied magnetic fields between 0 and 9 T and temperatures between 2 and 30 K for Nd0.5Sr0.5MnO3, Nd0.5Ca0.5MnO3, Sm0.5Ca0.5MnO3, Dy0.5Ca0.5MnO3 and Ho0.5Ca0.5MnO3 samples. All these compounds presented a Schottky-like anomaly at low temperatures,. Here, we report a general magnetic characterization and specific heat measurements in the full temperatures interval, between 2 and 300 K, for the same five samples. As far as we know, the high temperature interval of the specific heat measurements in these compounds have not been published yet.

## II Experimental methods

Polycrystalline samples of Nd0.5Sr0.5MnO3, Nd0.5Ca0.5MnO3 and Ho0.5Ca0.5MnO3 were prepared by the sol-gel method. Stoichiometric parts of Nd2O3 (Ho2O3) and MnCO3 were dissolved in HNO3 and mixed to an aqueous citric acid solution, to which SrCO3 or CaCO3 was added. The mixed metallic citrate solution presented the ratio citric acid/metal of 1/3 (in molar basis). Ethylene glycol was added to this solution, to obtain a citric acid/ethylene glycol ratio 60/40 (mass ratio). The resulting solution was neutralized to pH$``$7 with ethylenediamine. This solution was turned into a gel, and subsequently decomposed to a solid by heating at 400 oC. The resulting powder was heat-treated in vacuum at 900 oC for 24 hours, with several intermediary grindings, in order to prevent formation of impurity phases. This powder was pressed into pellets and sintered in air at 1050 oC for 12 hours.

Polycrystalline samples of Sm0.5Ca0.5MnO3 and Dy0.5Ca0.5MnO3 were prepared from stoichiometric amounts of Sm2O3 or Dy2O3, CaO, and MnO2 by standard solid-state reaction method. All the powders were mixed and ground for a long time in order to produce a homogeneous mixture. First, the mixture was heated at 927 oC for 24 hours and after that it was ground and heated at 1327 oC (72 hours) and 1527 oC (48 hours). X-ray diffraction measurements indicated high quality samples in all cases.

The magnetization measurements were done with a Quantum Design MPMS-5S SQUID magnetometer. Specific heat measurements were done with a Quantum Design PPMS calorimeter. The PPMS used the two relaxation time technique, and data was always collected during sample cooling. The intensity of the heat pulse was calculated to produce a variation in the temperature bath between 0.5 % (at low temperatures) and 2% (at high temperatures). Experimental errors during the specific heat and magnetization measurements were lower than 1 % for all temperatures and samples.

## III Results and Discussion

### A Magnetization measurements

Figure 1 shows the magnetization temperature dependence, measured with a 5 T applied magnetic field in field cooling conditions, in polycrystalline samples of Nd0.5Sr0.5MnO3, Nd0.5Ca0.5MnO3, Sm0.5Ca0.5MnO3, Dy0.5Ca0.5MnO3 and Ho0.5Ca0.5MnO3. The curves are plotted with a logarithmic scale in the y-axes to facilitate comparisons. Charge ordering transition temperatures (TCO) are indicated by arrows at 160, 250, 270, 280 and 271 K, respectively. These temperatures are associated to peaks in the magnetization curves, in agreement with previous reports , ,,. It is interesting to note that the relation between the charge ordering temperature and the antiferromagnetic ordering temperature (TN) changes from one sample to the other, ,,. In the first case they are approximately coincident, in the second and third cases the charge ordering temperatures are much higher, and in the fourth and fifth cases a long range antiferromagnetic transition is not observed.

The Nd0.5Sr0.5MnO3 sample presented a ferromagnetic transition at T$`{}_{C}{}^{}`$244 K and an antiferromagnetic transition at T$`{}_{N}{}^{}`$160 K. The Nd0.5Ca0.5MnO3 compound presented a strong magnetization maximum near TCO, but showed an unexpected minimum close to the antiferromagnetic transition temperature T$`{}_{N}{}^{}`$160 K. Usually an antiferromagnetic transition is accompanied by a maximum in the temperature dependence of magnetization. The antiferromagnetic transition in Sm0.5Ca0.5MnO3 presented a maximum at T$`{}_{N}{}^{}`$150 K. For temperatures lower than 10, 20 and 50 K the Nd0.5Sr0.5MnO3, Nd0.5Ca0.5MnO3 and Sm0.5Ca0.5MnO3 samples respectively, showed a sharp increase in the magnetization. This trend have been associated to a short range magnetic ordering of the intrinsic magnetic moment of Nd3+ ions. However, no long range ferromagnetic order of the Nd3+ ions was found in neutron diffraction measurements at these low temperatures, .

Differently from the three previous samples, the Dy0.5Ca0.5MnO3 and Ho0.5Ca0.5MnO3 compounds do not present a strong magnetization maximum at the charge ordering temperature. However, a clear inflection is observed at TCO for both samples, as revealed by the temperature derivative shown in the inset of fig. 1b. The existence of charge ordering in Dy0.5Ca0.5MnO3 and Ho0.5Ca0.5MnO3 was suggested by T. Terai et al.  after studies of magnetization and resistivity curves. As shown ahead our high temperature measurements of specific heat present peaks at around the same temperature interval of the suggested charge ordered transition.

### B Specific heat at high temperatures

Figure 2 shows specific heat measurements with a zero applied magnetic field from 2 to 300 K in the Nd0.5Sr0.5MnO3, Nd0.5Ca0.5MnO3, Sm0.5Ca0.5MnO3, Dy0.5Ca0.5MnO3 and Ho0.5Ca0.5MnO3 samples. In order to facilitate the visualization, the curves for Nd0.5Sr0.5MnO3 and Ho0.5Ca0.5MnO3 were displaced 20 J/mol K upside and downside in figure 2a, and the curve for Dy0.5Ca0.5MnO3 was displaced 20 J/mol K downside in figure 2b. Specific heat measurements give information about both lattice and magnetic excitations. At high temperatures the excitations from the lattice vibrations are dominant and decrease as the temperature decreases. The magnetic contribution can be obtained approximately by subtracting the lattice part from the experimental values.

Continuous lines in figure 2 represent the fitting of the thermal background, in the interval from 30 to 300 K, by the Einstein model given by:

$$C_{Einstein}=3nR\underset{i}{}a_i\left[\frac{x_i^2e^{x_i}}{\left(e^{x_i}1\right)^2}\right]$$
(1)
where $`x_i=T_i/T`$ . We used three optical phonons ($`i`$ =1, 2, 3) with energies $`\mathrm{T}_i`$ (in Kelvin) and relative occupations $`\mathrm{a}_i`$. The Einstein model for the specific heat considers the oscillation frequency (or energy) independently of the wave vector, which is a valid approximation for the optical part of the spectrum. The values of temperatures (energies) and relative occupations are shown in table 1. These values are similar to those reported, using the same model, by A. P. Ramirez et al. in a La0.37Ca0.63MnO3 sample and Raychaudhuri et al. in a Pr0.63Ca0.37MnO3 sample. However, we should point out that the values found are not unique because there are six free parameters during the fitting (one energy and one occupation coefficient for each oscillation mode). Besides, the thermal background determination is risky, particularly because it could have a ”tail” of any magnetic or charge ordering anomaly. Nonetheless, this seems to be the best possible trial to quantify the specific heat at high temperatures.

Figure 3 represents the differences between the experimental data and the fitted curves in figure 2. There is a maximum at 231 K for the Nd0.5Sr0.5MnO3 sample, which is correlated with the ferromagnetic transition at 250 K in the corresponding magnetization curve (figure 1). A second maximum, which could be partially associated to the antiferromagnetic and charge ordering transitions at 160 K, appears at 180 K. However, lattice parameters in this compound change rapidly between approximately 110 K and 250 K. The variation in lattice parameters changes the intensity of the interactions between the atoms and the oscillation frequency of the phonons mode, contributing to the specific heat in the second peak.

Figure 3 also shows that there is a maximum at 243 K for the Nd0.5Ca0.5MnO3 sample. This maximum correlates with the corresponding charge ordering temperature at 250 K. Differently from the Nd0.5Sr0.5MnO3 compound, there is no magnetic ordering in this high temperature interval for the Nd0.5Ca0.5MnO3 sample. However, in this latter case lattice parameters change very rapidly between approximately 200 and 250 K . An inflection point in the C vs. T curve (fig. 2a) appears at 141 K for the Nd0.5Ca0.5MnO3 sample. Besides, there is a maximum at 141 K in figure 3, but its height is relatively small compared to the maximum at 243 K. The Neel temperature corresponding to this compound is 160 K. Therefore, these results lead us to conclude that the specific heat variations due to the antiferromagnetic order are small compared to those induced in the charge ordering and ferromagnetic transitions.

For the Sm0.5Ca0.5MnO3, Dy0.5Ca0.5MnO3 and Ho0.5Ca0.5MnO3 samples, the maxima were found at 254 K, 281 K and 276 K, respectively. These experimental results and the resistivity measurements reported by Y. Tokura et al. and T. Terai et al. , are strong evidences that indicate the existence of CO transitions in the Sm0.5Ca0.5MnO3, Dy0.5Ca0.5MnO3 and Ho0.5Ca0.5MnO3 samples. However, electron diffraction studies would be needed to unambiguously classify this transition as CO. The other small peaks, comparable to the experimental error, were not found to be correlated to any magnetic transition.

Results in figure 3 allow us to calculate the variation in entropy ($`\mathrm{\Delta }`$S), within the limits of the model, associated to the charge ordering, ferromagnetic and antiferromagnetic transitions:

$$\mathrm{\Delta }S=\underset{Ti}{\overset{Tf}{}}\frac{\left(CC_{ph}\right)}{T}𝑑T$$
(2)
where $`\mathrm{T}_i`$ and $`\mathrm{T}_f`$ are two temperatures conveniently chosen to delimitate the interval of interest and $`\mathrm{C}_{ph}`$ is the specific heat due to the lattice oscillations.

For the Nd0.5Ca0.5MnO3 sample the entropy variation between 201 and 301 K, with H=0, was $`\mathrm{\Delta }S(`$T$`{}_{CO}{}^{})=`$ 2.0 J/(mol K). Raychaudhuri et al. reported an entropy variation, close to the charge ordering transition in the compound Pr0.63Ca0.37MnO3, of 1.8 J/(mol K) with zero applied magnetic field and 1.5 J/(mol K) with an 8 T magnetic field. On the other hand, Ramirez et al.  found $`\mathrm{\Delta }S(`$T$`{}_{CO}{}^{})=`$ 5 J/(mol K) in a La0.37Ca0.63MnO3 sample. All these results correspond to those expected for a charge ordering transition . Besides, the entropy variation (not related with phonons), calculated between 118 and 201 K for the Nd0.5Ca0.5MnO3 sample, was $`\mathrm{\Delta }S(`$T$`{}_{N}{}^{})=`$ 0.80 J/(mol K). We associate this smaller $`\mathrm{\Delta }\mathrm{S}`$ value to the antiferromagnetic order at T$`{}_{N}{}^{}=`$160 K.

The entropy variation between 133 and 274 K for the Nd0.5Sr0.5MnO3 sample was $`\mathrm{\Delta }S(`$T$`{}_{CO+FM}{}^{})=`$ 3.6 J/(mol K). Considering that the entropy variation associated to the charge ordering transition is the same as for the Nd0.5Ca0.5MnO3 sample, one finds that the entropy variation associated to the ferromagnetic transition is approximately 1.6 J/(mol K). Using a model proposed by J. E. Gordon et al.  we estimate that the change in entropy needed to a full ferromagnetic transition in the Nd0.5Sr0.5MnO3 sample is 12.45 J/(mol K). Therefore, the fact that the entropy variation found is 13 % of the theoretical value, suggests that only a small part of the spins order ferromagnetically. J. E. Gordon et al.  found also that the entropy variation associated to the ferromagnetic order in a Nd0.67Sr0.33MnO3 sample was approximately 10 % of the theoretical value. The entropy variation, not associated to lattice oscillations, between 200 K and 300 K, for the Sm0.5Ca0.5MnO3 and Dy0.5Ca0.5MnO3 samples were $`\mathrm{\Delta }S(`$T$`{}_{CO}{}^{})=`$ 1.15 J/(mol K) and 2.80 J/(mol K); between 216 K and 294 K for the Ho0.5Ca0.5MnO3 sample was $`\mathrm{\Delta }S(`$T$`{}_{CO}{}^{})=`$ 0.78 J/(mol K). The entropy variation for Sm0.5Ca0.5MnO3 and Dy0.5Ca0.5MnO3 are closer to the ones of the Nd0.5Ca0.5MnO3 sample. The $`\mathrm{\Delta }S(`$T$`{}_{CO}{}^{})`$ for the Ho0.5Ca0.5MnO3 sample is smaller than for the other samples, suggesting a different nature of the charge ordering transition.

Figure 4 shows the specific heat measurements, between 200 and 300 K, with a zero applied magnetic field (open symbols) and with H = 9 T (closed symbols) for the five measured samples. To allow comparisons the Nd0.5Ca0.5MnO3 curve was displaced 20 J/mol K upside; the Sm0.5Ca0.5MnO3, Dy0.5Ca0.5MnO3 and Ho0.5Ca0.5MnO3 curves were displaced 20, 40 and 60 J/mol K downside, respectively. The external magnetic field suppresses the peak around the ferromagnetic temperature in the Nd0.5Sr0.5MnO3 sample. Although a 9 T magnetic field is much smaller than a thermodynamic field given by kBTC / $`\mu _B`$, near the ferromagnetic phase transition the system is in an unstable condition, which is very sensitive to the external magnetic field. This fact explains the suppression of the peak in the presence of field in the specific heat curve and the decrease in the slope of the magnetization curve (not shown).

On the other hand, the external magnetic field seems to have no effect in the specific heat of the other four compounds. This result indicates that, although the sample magnetization increases with temperature close to the charge ordering transition, the CO phase do not appear to depend thermodynamically on the application of an external magnetic field. The increase of magnetization close to TCO is associated with an abrupt change in lattice parameters, affecting the distance among magnetic ions and hence their interactions. Our experiments show that an external magnetic field of 9 T is not strong enough to produce variations in the specific heat of the CO phase occurring in these samples.

## IV Conclusions

We have made a magnetic characterization of Nd0.5Sr0.5MnO3, Nd0.5Ca0.5MnO3, Sm0.5Ca0.5MnO3, Dy0.5Ca0.5MnO3 and Ho0.5Ca0.5MnO3 polycrystalline samples. Ferromagnetic, antiferromagnetic and charge ordering transitions in our samples agreed with previous reports. We also reported, to our knowledge for the first time, specific heat measurements with applied magnetic fields between 0 and 9 T and temperatures between 2 and 300 K in all these five samples. Each curve was successfully fitted at high temperatures by an Einstein model with three optical phonon modes. Close to the charge ordering and ferromagnetic transition temperatures the specific heat curves showed peaks superposed to the characteristic response of the lattice oscillations. The entropy variation corresponding to the charge ordering transition was higher than the one corresponding to the ferromagnetic transition. Near 160 K the specific heat curve showed an abrupt change in slope for the two compounds with Nd3+ ions, which were correlated to the corresponding antiferromagnetic transition. However, an external 9 T magnetic field seems to have no effect in the specific heat of the CO phase transition.

We thank Dr. P. N. Lisboa-Filho for the preparation of Nd0.5Sr0.5MnO3, Nd0.5Ca0.5MnO3 and Ho0.5Ca0.5MnO3 samples and the Brazilian science agencies FAPESP and CNPq for the financial support.


format: text

text:
# Spectroscopy of BL Lac Objects: new redshifts and mis-identified sources.

## 1. Introduction and Observations

At variance with other classes of AGN BL Lac Objects are characterized by quasi featureless optical spectra. This often hinders the determination of the redshift (and therefore of the distance) of these sources. In the Véron Cétty and Véron catalog (2001) of AGN there are $``$600 objects classified as BL Lacs and only for half of them the redshift is known. In order to contribute to the measurement of this fundamental parameter we are carrying out a systematic study of the optical spectra of objects with unknown or tentative redshift. We present here preliminary results from two campaigns at the European Southern Observatory (ESO) telescope performed in 2001 July and 2002 January.

The observations were gathered using the 3.6m telescope equipped with EFOSC2. Standard data reduction was performed using different packages in IRAF in order to obtain 1-dimensional wavelength and flux calibrated extracted spectra.

## 2. Results

### 2.1. New redshifts

PKS 0754+100 is a highly polarized and variable source that was discovered by Tapia et al. (1977) and classified as a blazar by Angel & Stockman (1980). A first tentative redshift proposed by Persic and Salucci (1986, z=0.66) was put in doubt after the detection of the host galaxy by Abraham et al. (1991) and by Falomo (1996), because at this redshift the host would be extremely luminous (M$`{}_{R}{}^{}`$-25).

Our spectrum (see fig. 1) shows clearly two emission lines at $`\lambda `$=4717.2 Å (EW=0.9 Å) and $`\lambda `$=6340.4 Å(EW=1.3 Å). Identification of these features with \[OII\](3727 Å) and \[OIII\] (5007 Å) yields a redshift z=0.266$`\pm `$0.001. This is consistent with the tentative redshift proposed by Falomo and Ulrich (2000, z$``$0.28). At z=0.266 the host galaxy luminosity becomes MR=-22.9. We note that at the same redshift there is a companion galaxy, 13 ′′ NE, corresponding to a projected distance of $``$ 70 kpc (Pesce et al. 1995).

1ES 0715-259 belongs to the Einstein Slew Survey (Elvis et al. 1992) and was classified as a BL Lac Object by Perlman et al. (1996). They assumed that a VLA radio source could be identified with the corresponding Slew Survey source if it fell within 80′′ of X-ray position. However in this case they reported an offset in the radio-X position of $`>`$200 ′′. The VLA radio map of 1ES 0715-259 (Perlman et al. 1996) shows a clear FR II morphology. The MMT spectrum obtained by Perlman et al.(1996) is featureless which motivated the BL Lacs classification. The value of the ratio Fx/Fr (logFx/Fr=-5.18) places 1ES 0715-259 near the borderline between RBL and XBL objects. There are no data of polarimetry for this source. HST (Urry et al. 2000) clearly resolved 1ES0715-259 into a point source surrounded by a small, rather elongated host galaxy.

We took an optical spectrum of the counterpart of the radio source. This shows clearly strong emission lines at $`\lambda `$=7114.5 (EW=85.3 Å), $`\lambda `$=7262.4 Å (EW= 14.3 Å) and $`\lambda `$=7332.0 Å (EW=46 Å) (see figure 1). Identification respectively with the Hβ, \[OIII\] 4959 Å and \[OIII\] 5007 Å yields a redshift z=0.464$`\pm `$0.001. The Hβ luminosity ( L$`_{H_\beta }`$=2.5 $`\times `$ 1043 erg s-1) is typical of radio loud quasar. On the basis of our optical spectroscopy and the radio properties we believe that this source is not a BL Lac but a radio loud quasar. The mis-identification by Perlman et al. (1996) is likely due to the very poor quality of their optical spectrum.

### 2.2. Lineless objects

For some relatively bright objects (V$``$ 15-16) despite the high signal to noise (S/N$``$200) the spectra appear lineless. We report here two examples: PG1553+ 113 and PKS1722+119 (see figure 2).

PG 1553+113 is a very bright optically selected source (Green et al. 1986) classified as a BL Lacs on the basis of the featureless spectrum (Miller & Green 1983) and the large optical variability ($`\mathrm{\Delta }`$m$``$ 2 mag, Miller et al. 1988). Falomo & Treves (1990) showed that the emission feature in the UV region identified as Lyman-$`\alpha `$ at z=0.36 by Miller & Green (1983) is spurious. Falomo & Treves (1990) did not find absorption features with an EW upper limit of 1 Å.

Our high signal to noise optical spectrum of the object confirms the absence of lines with a much lower of EW ($``$0.1 Å). The redshift of this source remains unknown. Assuming that it is hosted by a galaxy typical of BL Lacs (MR=-23.7 Re=10 kpc, Urry et al. 2000) from our spectrum we set a lower limit of z$`>`$0.3 consistent with the fact that the source is unresolved in HST images (Urry et al. 2000).

PKS1722+119 is a highly polarized X-ray selected BL Lac Object originally discovered by Uhuru. Griffith et al. (1989) claimed the observation of an absorption feature near 6000 Å which could be Na I (5892 Å) redshifted at z=0.018. This was disproved by Brissenden et al. (1990) who obtained a completely lineless spectrum, consistently with more recently results by Vèron Cètty and Vèron 1993 and by Falomo et al. 1993.

The spectrum presented here (see fig. 2) is lineless and we can set a lower limit of z$`>`$0.3. Also this source is unresolved by HST (Urry et al. 2000).

### 2.3. Mis-identified sources

From the spectra examined so far we found two cases UM493 and 1620+103 where we observe the typical stellar absorption lines of the Balmer series clearly indicating a wrong classification. Probably the mis-identifications (e.g.Vèron-Cètty & Vèron 2001, Hewitt & Burbidge 1993) were due to the poor informations on these sources, classification being based only on modest optical spectroscopy. Our high signal-to noise spectra clarify their nature and we can definitively classify these sources as stars of spectral type A.

## References

Abraham R.G., McHardy I.M., Crawford C.S. 1991, MNRAS, 252, 482

Angel J.R.P. & Stockman H.S. 1980, ARA&A, 18, 321

Brissenden R.J.V. et al. 1990, ApJ, 350, 578

Elvis M. et al. 1992, ApJS, 80, 257

Falomo R. & Treves A. 1990, PASP, 102, 1120

Falomo R. et al. 1993, ApJ, 411, L63

Falomo R. 1996, MNRAS, 283, 241

Falomo R. & Ulrich M-H 2000, A&A, 357, 91

Green R.F., Schmidt M. and Liebert J. 1986, ApJS, 61, 305

Griffiths R.E. et al. 1989, MNRAS, 240, 33

Hewitt A. & Burbidge G., 1993, ApJS, 87, 451

Miller H.R. and Green R.F. 1983, Bull. AAS, 15, 957

Miller H.R. et al. 1988, in proc. IUE Symposium, Greenbelt, ESA SP-281, p303

Perlman E.S. and Stocke J.T. et al. 1996, ApJS, 104, 251

Persic M. & Salucci P. 1986, Struc. and Evol. of AGN, ed Giuricin et al. , p657

Pesce,J.E., Falomo,R. Treves,A., 1995, AJ, 110, 1554

Tapia et al. 1977, ApJ, 215, 71

Urry C.M. et al. 2000, ApJ, 532, 816

Vèron-Cètty M.P. and Vèron P. 2001, A&A, 374, 92


text:
# Analysis of the Anomalous Brillet and Hall Experimental Result

## 1 INTRODUCTION

### 1.1 Background

In 1978, Brillet and Hall used a Fabry-Perot interferometer that rotated with respect to the lab to measure the isotropy of space (i.e., of the speed of light). As the apparatus turned, anisotropic speed of light would cause otherwise coherent light waves to change phase and partially interfere, thereby reducing the transmitted wave intensity. A servo adjusted the frequency (and therefore the wavelength and phase) to maximize the transmitted wave intensity. Thus, the change in frequency was a direct measure of the degree of anisotropy.

Brillet-Hall found a ”null” effect at the $`\mathrm{\Delta }f/f`$= 3X10-15 level, where $`f`$ is the frequency of the electromagnetic wave employed and $`\mathrm{\Delta }f`$ is the amplitude of the change in $`f`$ as the apparatus is rotated. However, to obtain this result they subtracted out an anomalous non-null signal of approximate amplitude 2X10-13 at twice the apparatus rotation rate. This signal had effectively constant phase in the laboratory (earth) frame, and though “persistent”, was deemed “spurious” without further comment.

In 1981, Aspden noted that this anomalous signal corresponded closely to what pre-relativity physicists would have expected solely from the earth surface velocity. According to logic employed by Michelson and Morley, $`\mathrm{\Delta }f/f`$ is related to $`1/2v^2/c^2`$, where $`v`$ is the speed of the test apparatus relative to the frame in which light speed is isotropic. Aspden suggested that perhaps in some little understood way, the earth surface speed might be somehow different in nature than the solar and galactic orbit speeds. His analysis implied that, if light speed were anisotropic, $`\mathrm{\Delta }f`$/$`f`$ would have an amplitude of (.262)($`1/2v^2/c^2)`$ . Setting this result equal to the measured signal of $``$2X10-13, one finds a $`v`$ value of approximately 362 m/sec. The earth surface velocity at the test site due to rotation about the earth’s axis is 355 m/sec.

The present author believes Aspden’s analysis to be flawed111 Aspden considered the change in angle of the output signal due to “motion of the mirror relative to the light reference frame” to be responsible for the non-null signal variation observed. However, the output signal in the Brillet-Hall experiment was monitored only for frequency change, and for an observer and apparatus fixed in the same frame, this does not vary with frame orientation, even in a Galilean analysis. (Wave speed and wavelength would change, but not frequency.) See Section 4.2 herein, ref. , and ref. .. Nevertheless, the detected signal is within a factor of two of what Michelson-Morley would have attributed to the earth surface velocity. (It turns out, as we will see, that $`1/2v^2/c^2`$ would be peak-to-peak in the Brillet-Hall test, and half of this would be the amplitude.) In the present article, this signal is evaluated using general relativity/differential geometry in which off diagonal time-space components occur in the rotating (earth) frame metric. For such a metric, time is not orthogonal to space, and this is referred to as non-time-orthogonality. As shown below, non-time-orthogonal (NTO) analysis, implies an anisotropy in the speed of light in rotating frames, which when applied to the earth, leads to a prediction of the observed $``$2X10-13 signal.

In this context, we note that as a result of studies on global positioning system (GPS) electromagnetic signal data, recognized world leading GPS expert Neil Ashby recently noted in Physics Today,

“ .. the principle of the constancy of c cannot be applied in a rotating reference frame ..” .

He has also stated,

“Now consider a process in which observers in the rotating frame attempt to use Einstein synchronization \[constancy of the speed of light\] ….. Simple minded use of Einstein synchronization in the rotating frame … thus leads to a significant error”.

It is emphasized preemptively that NTO analysis does not contravene the geometric foundation of relativity theory. In fact, it is a direct consequence of straightforward application of differential geometry to cases where time is not orthogonal to space. These include rotating frames and spacetime around a star possessing angular momentum (in which off diagonal space-time metric components also exist.) NTO methodology leaves unchanged all analyses of systems in which time is orthogonal to space. Such systems comprise free fall frames in gravitational orbits and the vast majority of all other relativistic systems. Hence, the null signals obtained in many tests for effects of the solar or galactic orbit speeds are not in conflict with NTO analysis.

Further, NTO analysis is in complete agreement with all known experiments,,,. Of these experiments, only the Brillet-Hall test has been capable of monitoring any effect due to non-time-orthogonality of the rotating earth fixed reference frame222 Arguments can, however, be made that the Sagnac effect (see ref. ) results from such non-time-orthogonality..

### 1.2 Analysis Approach

The steps in analysis presented herein are as follows.

1. Introductory overview of NTO analysis (Section 2.)
2. Determination of physical light speed (i.e., measured with physical instruments) in the circumferential and radial directions of a rotating frame via NTO analysis. (Section 3.1), and from them, the difference in time delay for the two directions (Sections 3.2, 3.3, and 3.4.)
3. A review of the Brillet and Hall apparatus, analysis of the result predicted by the NTO approach due to differences in time delay of light signals in the radial and circumferential direction, and a comparison with the actual results. (Section 4.) The signal predicted is seen to be of the same order of magnitude but off by somewhat less than a factor of two.
4. An analysis of effects on the light signal other than time delay predicted by the NTO approach. (Section 5.) These other effects are transverse to the light propagation direction, whereas the signal time delay is a parallel direction effect. The two combined are then shown to result in a predicted signal strikingly close to that measured by Brillet and Hall.
5. An explanation for why null signals are predicted by NTO analysis for solar and galactic orbital speeds but not for the earth surface speed due to rotation about its axis. (Section 6.)
6. Summary and Conclusions (Section 7.)

## 2 OVERVIEW OF NTO ANALYSIS

The most widely accepted transformation333 There is significant background behind this transformation that is discussed in greater depth in refs  and . We shall tentatively accept this transformation as valid and find in Section 5 that it results in a correct prediction for the experiment under consideration. from the non-rotating (lab, upper case symbols) frame to a rotating (lower case) frame is

$$\begin{array}{c}cT=ct\hfill \\ R=r\hfill \\ \mathrm{\Phi }=\varphi +\omega t\hfill \\ Z=z\hfill \end{array}$$
(1)
where $`\omega `$ is the angular velocity of the rotating frame and cylindrical spatial coordinates are used. The coordinate time $`t`$ for the rotating system equals the proper time of a standard clock located in the lab, which equals the time on a standard clock at the origin (axis of rotation) of the rotating frame.

Substituting the differential form of (1) into the line element in the lab frame

$$ds^2\text{ }=\text{ }\text{ }c^2dT^2+\text{ }dR^2\text{ }+\text{ }R^2d\mathrm{\Phi }^2\text{ }+\text{ }dZ^2$$
(2)
results in the line element for the rotating frame444 The first person to determine this metric for rotating frames may have been Paul Langevin , who used it successfully to analyze the Sagnac experiment, though he did not suggest it would lead to a non-null Michelson-Morley result for rotation. See further comments in this regard by the present author in refs , , and .

$$ds^2=c^2(1\frac{r^2\omega ^2}{c^2})dt^2+dr^2+r^2d\varphi ^2+\mathrm{\hspace{0.33em}2}r^2\omega d\varphi dt+dz^2=g_{\alpha \beta }dx^\alpha dx^\beta .$$
(3)
Note that the metric in (3) is not diagonal, since $`g_{\varphi t}0`$, and this implies that time is not orthogonal to space (i.e., an NTO frame.) This mathematical feature becomes physically significant due to the fact that in a rotating frame it is not possible to synchronize clocks continuously without producing an NTO coordinatization. (See refs.  and .)

Time on a standard clock at any fixed 3D location in the rotating frame, found by taking ds$`{}_{}{}^{2}=c^2d\tau `$ and dr = d$`\varphi `$ = dz = 0, is

$$d\tau =d\widehat{t}=\sqrt{g_{tt}}dt=\sqrt{1r^2\omega ^2/c^2}dt,$$
(4)
where the caret over dt indicates physical (proper) time (i.e., time measured with physical world standard clocks fixed in the rotating frame.) Obviously, time becomes imaginary for $`r>c/\omega `$, speed becomes faster than light, and physical clocks (or any material body) cannot exist at such locations.

## 3 LIGHT SPEEDS AND TRAVEL TIMES

### 3.1 Light Speeds in Two Directions

#### 3.1.1 Circumferential Direction Light Speed

For light ds2 = 0. Inserting this into (3), taking dr=dz=0, and using the quadratic equation formula, one obtains a light coordinate velocity (generalized coordinate spatial grid units per coordinate time unit) in the circumferential direction

$$v_{light,coord,circum}=\frac{d\varphi }{dt}=\omega \pm \frac{c}{r},$$
(5)
where the sign before the last term depends on the circumferential direction of travel of the light ray. The physical velocity (the value one would measure in experiment using standard meter sticks and clocks in units of meters per second) is found from this to be555 The texts and articles listed in ref.  are among those that discuss physical vector and tensor components (the values one measures in experiment) and the relationship between them and coordinate components (the mathematical values that depend on the generalized coordinate system being used.) This relationship is $`u^{\widehat{\mu }}=\sqrt{g_{\underset{¯}{\mu }\underset{¯}{\mu }}}u^\mu `$ where the caret over the index indicates a physical vector component and underlining implies no summation.

$$v_{light,phys,circum}=\frac{\sqrt{g_{\varphi \varphi }}d\varphi }{\sqrt{g_{tt}}dt}=\frac{r\omega \pm c}{\sqrt{1\frac{\omega ^2r^2}{c^2}}}=\frac{v\pm c}{\sqrt{1\frac{v^2}{c^2}}},$$
(6)
where $`v=\omega r`$ is the circumferential speed of a point fixed in the rotating frame as seen from the lab. Note that for rotation, the physical speed of light is not invariant or isotropic, and that this lack of invariance/isotropy depends on $`\omega `$, the angular velocity. If $`\omega `$=0, light speed is isotropic and invariant, the metric (3) is diagonal, and time is orthogonal to space.

#### 3.1.2 Radial Direction Light Speed

For a radially directed ray of light, d$`\varphi `$ = dz = 0, and ds = 0. Solving (3) for dr/dt one obtains

$$\frac{dr}{dt}=c\sqrt{1r^2\omega ^2/c^2}.$$
(7)
Since $`g_{rr}`$ = 1, the physical component (measured with standard meter sticks) for radial displacement $`d\widehat{r}`$ equals the coordinate radial displacement dr. The physical (measured) speed of light in the radial direction is therefore

$$v_{light,phys,radial}=\frac{d\widehat{r}}{d\widehat{t}}=\frac{\sqrt{g_{rr}}dr}{\sqrt{g_{tt}}dt}=\frac{dr}{\sqrt{1r^2\omega ^2/c^2}dt}=c.$$
(8)
For details on NTO analysis, the reader is referred to Klauber,.

### 3.2 Round Trip Travel Times

#### 3.2.1 Circumferential Round Trip Travel Time

In a given direction the outward and return round trip time over a path length $`l`$ is

$$t_{RT}=\frac{l}{v_+}+\frac{l}{v_{}},$$
(9)
where notation should be obvious. Inserting (6) into (9), and using absolute values for speed in (6), one obtains the circumferential round trip light time

$$t_{RT,circum}=\frac{l\sqrt{1v^2/c^2}}{cv}+\frac{l\sqrt{1v^2/c^2}}{c+v}.$$
(10)
This readily becomes

$$t_{RT,circum}=\frac{2lc\sqrt{1v^2/c^2}}{c^2v^2}=\frac{2l}{c\sqrt{1v^2/c^2}}\frac{2l}{c}\left(1+\frac{1}{2}\frac{v^2}{c^2}\right).$$
(11)
#### 3.2.2 Radial Round Trip Travel Time

From (8), we know that in rotating frames for the radial direction $`v_+`$ = $`v_{}=c`$. Using these values in (9), one finds the round trip (assuming $`g_{tt}`$ is effectively constant over short distances) time, as measured on standard clocks fixed in the rotating frame, to be

$$t_{RT,radial}=\frac{2l}{c}.$$
(12)
This is no different from that of inertial systems in special relativity.

### 3.3 Round Trip Time Difference

Consider the earth as a rotating frame with either a Michelson-Morley or Brillet-Hall experiment mounted on the earth’s surface. Each apparatus turns so that the perpendicular and parallel (to the earth surface speed around its axis) directions are alternately tested for the round trip light travel times.

The difference between the parallel and perpendicular directions for rotation is (11) minus (12),

$$\mathrm{\Delta }t_{RT,NTOframerotation}\frac{2l}{c}\left(\frac{1}{2}\frac{v^2}{c^2}\right)=(t_{RT,nomotion})\left(\frac{1}{2}\frac{v^2}{c^2}\right)=\frac{1}{2}\beta ^2(t_{RT,nomotion}).$$
(13)
### 3.4 Rotation vs Translation

For the NTO rotating frame analysis, the radial direction corresponds to the perpendicular direction for a translating system. Note, the travel times for this direction between Galilean translation (reviewed in Appendix B) and NTO relativistic rotation are not the same, and differ by a factor of $`1/2`$ before the $`\beta ^2`$. (Compare (12) to (46).)

Note also that the travel times between Galilean translation and NTO relativistic rotation in the parallel to velocity direction also differ by a factor of $`1/2`$ before $`\beta ^2`$ (compare (11) to (47).)

However, the difference between the parallel and perpendicular directions for both cases is (13), i.e., the same as (48). One can thus conclude that according to NTO analysis, the Michelson-Morley and Brillet-Hall experiments performed on a rotating frame (i.e., the earth) would detect signals of the same magnitude as predicted for Galilean translation.

## 4 BRILLET-HALL

### 4.1 Brillet-Hall’s Fabry-Perot Interferometer

A Fabry-Perot interferometer consists of two partially transmitting, partially reflecting mirrors separated by an “etalon” (of length $`d)`$ and facing one another as shown in Figure 1. Light entering the left side of the Fabry-Perot interferometer in Figure 1 emerges as a series of rays on the right side, and creates a fringe pattern. For reference, a derivation of the well known constructive interference relation (shown inside the box of Figure 1 and repeated below as (14)) is provided in Appendix A.

To help make a point later, we emphasize here that a wave front entering from the left is visualized as a collection of Huygen’s type wavelets. The particular ray components of such wavelets entering the left glass section at angle $`\alpha `$ produce constructive interference at an exit angle of $`\alpha `$. Other ray components of the same wave front, entering at other angles, produce various degrees of destructive interference. In essence then, a change in angle of incidence of the centerline of a macroscopic beam of light would have no effect on the angle $`\alpha `$ at which constructive interference occurs.

The angle $`\alpha `$ of peak intensity for light exiting the Fabry-Perot interferometer is found from

$$\mathrm{cos}\alpha =\frac{m\lambda }{2d},$$
(14)
where $`m`$, the mode number, is an integer, and different values for $`m`$ correspond to different angles of local peaks in intensity.

### 4.2 Fabry-Perot and the Brillet-Hall Experiment

The incoming light ray depicted in the Fabry-Perot interferometer of Figure 2 has an incident angle of zero. The upper part of the figure shows the right traveling incident ray (ray #1), its left traveling reflection off of the RHS (ray #1r), and the right traveling reflection of the ray #1r (ray #2). For simplicity, additional successive rays (#2r, #3, #3r, etc.) that would result from reflection of ray #2 and subsequent reflections are not shown. Isotropy of light speed is assumed in the upper part of Figure 2, and the wavelength has been adjusted (by varying the incident light frequency) so that ray #1 and ray #2 are in phase when they strike the RHS. The portions of #1 and #2 that exit the RHS will then combine to maximize intensity of the resultant beam. Higher order rays (#3, #4, etc. not shown) would then be in phase with #1 and #2 as well.

From (14), the condition for constructive interference with $`\alpha =0`$ and isotropic light speed is

$$\lambda =\frac{2d}{m}.$$
(15)
Note that a change in either the etalon length or the wavelength will reduce the constructive interference at $`\alpha `$ = 0 and thus, the intensity sensed.

For anisotropic light speed (lower part of Figure 2), the round trip time for a given wave maximum to reflect from the right mirror and return to the right mirror would be greater than for the isotropic case. Hence, the phase relation between waves #1 and #2 would vary as the apparatus turned, and thus so would intensity of the light entering the sensor. By judiciously changing the wavelength as the orientation changed, one could, in principle, maintain maximum constructive interference.

Brillet-Hall used this property of the Fabry-Perot interferometer to measure light speed anisotropy. They employed a He Ne generated laser that was successively reflected off two mirrors and into an Fabry-Perot interferometer, as shown in Figure 3. A sensor at fixed angle relative to the Fabry-Perot interferometer monitored the output radiation intensity, and the signal from it was fed into a servo unit. The entire apparatus shown in Figure 3 was mounted on a 95 cm by 40 cm by 12 cm granite table. The table was then turned about an axis perpendicular to the table surface at the rate of one revolution every 10 seconds. As the apparatus turned, the servo unit continually adjusted the frequency (and hence the wavelength) of the generated laser to maintain sensor intensity at a maximum. Variation in this frequency was therefore a direct measure of the anisotropy of light speed.

Thus, if light speed were anisotropic in rotating frames as NTO analysis suggests, the $`1/2\beta ^2`$ of (13) would represent the fractional change in time between maxima of the first and second right traveling waves as they exit the etalon on the RHS. That is, it represents the phase lag of the second wave behind the first upon arrival at the RHS. Changing the wavelength (done by changing the frequency) can reduce this phase lag to zero and maximize the detected intensity. Hence, the amount the servo changes frequency in order to maximize intensity indicates the magnitude of $`1/2\beta ^2`$ and therefore the magnitude of $`v`$ detectable (i.e., the degree of anisotropy.)

### 4.3 Time Delay Due to Anisotropy in the Brillet-Hall Test

We saw in Section 3.4 that the time delay variation imparted to wave #2 in Figure 2 as a result of changing orientation of the Fabry-Perot interferometer through 90o for both translating frames in a Galilean universe and NTO rotating frames is (13), i.e.,

$$\mathrm{\Delta }t_{RT,movingframe}\frac{1}{2}\beta ^2(t_{RT,nomotion}).$$
(16)
Note, however, that the Michelson-Morley apparatus has two perpendicular legs and results in (16) representing a single-sided amplitude variation in the fringing as the apparatus turns. This is a result of the lead time of one leg over the other becoming a lag time when the apparatus is turned 90o. In the Brillet-Hall experiment, on the other hand, the round trip time difference, although numerically equal to (16), is a peak-to-peak (double-sided) amplitude. The Brillet-Hall apparatus uses a Fabry-Perot interferometer as a single arm and the time delay varies between zero for the radial direction and the same value (given by (11)) for both the positive and negative circumferential directions.

### 4.4 Brillet-Hall: NTO Time Delay Prediction

As Brillet-Hall noted, with a single arm interferometer, any variation due to anisotropy would have a peak-to-peak magnitude (total change) represented by the magnitude of (16). That is, since an increase in time delay must be compensated for by an increase in wavelength (i.e., a decrease in frequency),

$$\frac{\mathrm{\Delta }t}{t}\frac{\mathrm{\Delta }\lambda }{\lambda }\frac{\mathrm{\Delta }f}{f}\frac{1}{2}\beta ^2\text{(total}\text{change, Brillet-Hall)}.$$
(17)
This contrasts with the Michelson-Morley experiment in which there were two perpendicular arms and the peak-to-peak signal would be twice the magnitude of (17). In both experiments, this detected signal would vary with 2X the apparatus rotation rate. In the Brillet-Hall experiment the amplitude of the fractional change in frequency $`\mathrm{\Delta }f/f`$ of the HeNe laser radiation source employed may thus be expected, to lowest order, to be half of the absolute value of (17), i.e

$$\frac{\mathrm{\Delta }f}{f}\frac{1}{4}\beta ^2\text{(amplitude, Brillet-Hall)}.$$
(18)
At the location of the Brillet-Hall test, the earth surface speed is 355 m/sec. So, for $`c=`$ 3X108 m/sec one finds $`1/2\beta ^2`$ = 7.00 X10-13 (peak-to-peak). The amplitude of the fractional time delay $`\mathrm{\Delta }t/t`$ (= $`|\mathrm{\Delta }f/f|`$ ) between successive exiting waves obtained as the apparatus rotates would be half of this or 3.5 X10-13. This signal should have constant phase angle relative to the lab (earth) frame.

### 4.5 Billet and Hall: Test Result

The mean value of the anomalous signal Brillet-Hall measured was $`\mathrm{\Delta }f/f2\times 10^{13}`$ at twice the table rotation rate and constant phase angle in the lab. Data varied between –42% and +52% of the mean.

The detected mean value is approximately 57% (2/3.5) of that predicted by NTO analysis of the signal time delay due to an anisotropy induced by the earth surface speed. Considering the experiment covered a range that extended to 10-15, one could consider this result alone to be quite significant. However, there are additional effects on the Brillet-Hall signal predicted by NTO analysis that modify the prediction and bring it significantly closer to the measured value. These are discussed in the following section.

## 5 FURTHER NTO EFFECTS ON BRILLET-HALL SIGNAL

### 5.1 Qualitative Considerations

We emphasize that the effect of the earth surface speed on the speed of light in (6) is due to non-time-orthogonality in a rotating frame and is not the result of an “ether wind”. However, it can aid in analysis and visualization of the effect described below if one temporarily ignores the second order contribution in (6), and thinks of the effect somewhat classically. Then the earth surface velocity v (magnitude $`v)`$ can be effectively considered as a simple vector addition to the usual velocity (magnitude $`c)`$ of a light ray.

In Figure 4 one sees the NTO effect of v in the Brillet-Hall experiment when it is transverse to the path of a beam of light. In the left side of the figure the portion of the incoming laser beam between the two mirrors, which we designate as “leg 2”, is perpendicular to the direction of v. The result is a deflection of the beam by an amount $`\mathrm{\Delta }_2`$ at the right hand mirror. Since the mirror is aligned at 45o to leg 2, the same $`\mathrm{\Delta }_2`$ deflection of the beam will occur to the beam as it enters, and leaves, the Fabry-Perot interferometer. Hence, the peak intensity fringe that would otherwise strike the sensor directly is moved to the right in the figure by the amount $`\mathrm{\Delta }_2`$.

Similar effects can be seen in the right hand side of the figure for leg 1 ($`\mathrm{\Delta }_1`$ deflection) and leg 3 ($`\mathrm{\Delta }_3`$ deflection) when the apparatus mounting table is turned by 90o relative to the earth surface. Note there is an additional effect of a change in angle of incidence for the ray as it traverses leg 3. However, as noted in section 4.1, this has no effect on the angle at which constructive interference occurs, and for our purposes, can therefore be ignored.

The important point to note is that we now have four separate effects on intensity detected by the sensor that are due to the earth surface speed. These are

1) Time delay (Section 4.4),

2) $`\mathrm{\Delta }_1`$ from transverse effect on leg 1,

3) $`\mathrm{\Delta }_2`$ from transverse effect on leg 2, and

4) $`\mathrm{\Delta }_3`$ from transverse effect on leg 3.

### 5.2 Quantification of Transverse Effects

Taking the orientation for v relative to the test table of the left side of Figure 4 as $`t`$ = 0, and using the dimensions for each leg shown in Figure 3, the lowest order value for $`\mathrm{\Delta }_2`$ is

$$\mathrm{\Delta }_2=d_2\left(\frac{v}{c}\right)\mathrm{cos}\omega _Tt,$$
(19)
where $`\omega _T`$ is the rate at which the table turns in radians/sec.

Similarly, we have for the other two legs,

$$\mathrm{\Delta }_1=d_1\left(\frac{v}{c}\right)\mathrm{cos}\left(\omega _Tt\frac{\pi }{2}\right),$$
(20)
and

$$\mathrm{\Delta }_3=d_3\left(\frac{v}{c}\right)\mathrm{cos}\left(\omega _Tt\frac{\pi }{2}\right).$$
(21)
Note each of these quantities represents the respective displacement of the entire fringe pattern at the sensor location to the right of the solid line (isotropic peak intensity path) entering the sensor in Figures 3 and 4.

### 5.3 Cumulative Result of NTO Effects

From Figure 1, it is apparent that a decrease in wavelength increases the peak intensity angle $`\alpha `$. Further, from Figure 2, one sees that the time delay due to anisotropy is similar in effect to a decrease in wavelength, and thus also increases the peak intensity angle. The effect of time delay on the fringe pattern is shown in the left side of Figure 5, i.e., it moves the fringes, including the peak intensity fringe, outward.

The total displacement due to transverse effects is

$$\mathrm{\Delta }_{Transv}=\mathrm{\Delta }_1+\mathrm{\Delta }_2+\mathrm{\Delta }_3.$$
(22)
The effect of positive $`\mathrm{\Delta }_{Transv}`$ on the interference fringes is depicted in the right side of Figure 5, i.e., it displaces the fringes to the right. Negative $`\mathrm{\Delta }_{Transv}`$ displaces the fringes to the left.

The peak intensity angle relation (14) can be expanded as

$$\frac{m\lambda }{2d}=\mathrm{cos}\alpha =\sqrt{1\mathrm{sin}^2\alpha }1\frac{1}{2}\alpha ^2,$$
(23)
since $`\alpha `$ is small. Consider, for simplicity, the ideal isotropic case in which a local peak intensity is at $`\alpha `$ = 0 and is directed precisely at the center of the sensor. (Imagine the sensor at the origin in the left side of Figure 5, light speed is isotropic, and the peak intensity is located at the origin.) From (23) this means

$$\lambda _0=\frac{2d}{m},$$
(24)
where $`\lambda _0`$ is the wavelength in the isotropic case that maximizes intensity at the center of the Fabry-Perot interferometer. (For the case where $`\alpha `$ is not initially zero, see Appendix C.) Thus, from (23) and (24), we have

$$\frac{\lambda }{\lambda _0}=\frac{\lambda _0+\mathrm{\Delta }\lambda }{\lambda _0}1\frac{1}{2}\alpha ^2\frac{\mathrm{\Delta }\lambda }{\lambda _0}\frac{1}{2}\alpha ^2.$$
(25)
Note that $`\mathrm{\Delta }\lambda `$ is always negative, and this makes sense physically. A decrease in $`\lambda `$ means an increase in $`|\alpha |`$.

Consider (25) with regard to the time delay only case of Figure 2, i.e., we introduce anisotropy into our heretofore isotropic picture, but ignore transverse effects. When the Fabry-Perot interferometer is aligned in the direction of v, the effect is similar to that of decreased $`\lambda `$. From our analysis of Section 4.4 and (25), we expect the relative change in wavelength to vary with an amplitude of $`\beta ^2/4`$ between zero and a negative value of $`\beta ^2/2`$ at twice the rate of test table rotation. From (17) and (25) then, the time delay effect on wavelength is

$$\frac{\mathrm{\Delta }t}{t}\frac{\mathrm{\Delta }\lambda }{\lambda _0}\frac{1}{4}\left(\frac{v}{c}\right)^2\left(\mathrm{cos}\left(2\omega _Tt\right)+1\right)\frac{1}{2}(\alpha _{TimeDelay})^2.$$
(26)
Trigonometry then gives us

$$\left(\alpha _{TimeDelay}\right)^2\frac{1}{2}\left(\frac{v}{c}\right)^2\left(\mathrm{cos}\left(2\omega _Tt\right)+1\right)=\left(\frac{v}{c}\right)^2\mathrm{cos}^2\left(\omega _Tt\right),$$
(27)
and the variation in $`\alpha _{TimeDelay}`$ is

$$\alpha _{TimeDelay}\pm \frac{v}{c}\mathrm{cos}\left(\omega _Tt\right).$$
(28)
The displacement variation at the sensor location due to time delay alone, where $`d_s`$ is the distance from the exit end of the Fabry-Perot etalon to the sensor (see Figure 3), is thus

$$\mathrm{\Delta }_{TimeDelay}(d+d_s)\alpha _{TimeDelay}(d+d_s)\frac{v}{c}\mathrm{cos}\left(\omega _Tt\right).$$
(29)
From Figures 2, 4 and 5, we can convince ourselves that $`\mathrm{\Delta }_{TimeDelay}`$ at the sensor location due to time delay is 180o out of phase with $`\mathrm{\Delta }_2`$. Hence, from (19), we see that the proper sign in (28) and (29) is the negative one.

The total displacement then is

$$\mathrm{\Delta }_{Tot}=\mathrm{\Delta }_{TimeDelay}+\mathrm{\Delta }_{transv}=\mathrm{\Delta }_{TimeDelay}+\mathrm{\Delta }_1+\mathrm{\Delta }_2+\mathrm{\Delta }_3$$
(30)
or

$$\mathrm{\Delta }_{Tot}=\frac{v}{c}\left((d+d_s)\mathrm{cos}\left(\omega _Tt\right)+d_1\mathrm{cos}\left(\omega _Tt\frac{\pi }{2}\right)+d_2\mathrm{cos}\left(\omega _Tt\right)d_3\mathrm{cos}\left(\omega _Tt\frac{\pi }{2}\right)\right).$$
(31)
This displacement must be corrected for by the servo to maximize intensity. That is, the servo must change the frequency, and hence the wavelength, to bring $`\mathrm{\Delta }_{Tot}`$ back to zero. In other words, the induced change in wavelength must be such that if acting alone (in an isotropic world) it would produce an angle $`\alpha _{Tot}`$ that would move the peak intensity fringe by $`\mathrm{\Delta }_{Tot}`$. This angle is then

$$\alpha _{Tot}\frac{\mathrm{\Delta }_{Tot}}{d+d_s}\frac{v}{c}\left(\left(\frac{d_2(d+d_s)}{(d+d_s)}\right)\mathrm{cos}\left(\omega _Tt\right)+\frac{d_1d_3}{(d+d_s)}\mathrm{cos}\left(\omega _Tt\frac{\pi }{2}\right)\right).$$
(32)
This can be represented as

$$\alpha _{Tot}=\frac{v}{c}A\mathrm{cos}\left(\omega _Tt\varphi \right)$$
(33)
where

$$A=\frac{\sqrt{(d_2dd_s)^2+(d_1d_3)^2}}{d+d_s}$$
(34)
and $`\varphi `$ is readily determinable. From (25), and a standard trigonometric relation, we have

$$\frac{\mathrm{\Delta }f}{f_0}\frac{\mathrm{\Delta }\lambda }{\lambda }=\frac{1}{2}\left(\alpha _{Tot}\right)^2=\frac{1}{2}\left(\frac{v}{c}\right)^2A^2\mathrm{cos}^2\left(\omega _Tt\varphi \right)=\frac{1}{4}\left(\frac{v}{c}\right)^2A^2\left(1+\mathrm{cos}\left(2\omega _Tt2\varphi \right)\right).$$
(35)
Thus, the fractional change in frequency induced by the servo varies at twice the table rotation rate at amplitude

$$\frac{\mathrm{\Delta }f}{f_0}|_{amplitude}=\frac{1}{4}\left(\frac{v}{c}\right)^2A^2.$$
(36)
### 5.4 Comparison of NTO Theory with Brillet-Hall Experiment

Brillet-Hall noted that the Fabry-Perot interferometer they used had $`d`$ = 30.5 cm, but it is unfortunate that the dimensions $`d_1`$, $`d_2`$, $`d_3`$, and $`d_s`$ were not specified. However, from the schematic of their Figure 1, and knowing the dimensions of the table illustrated therein, one can make reasonable estimates for these values. For example, they appear to have values in centimeters which may be reasonably close to ($`d_2`$ appears to be between 20 and 25 cm)

$$d_s=10d_1d_3=25d_2=22.5.$$

Then $`A^2`$ = .58 and

$$\frac{\mathrm{\Delta }f}{f_0}|_{amplitude}=.138\left(\frac{v}{c}\right)^2=2.03\times 10^{13}.$$
(37)
As noted in Section 4.5, the Brillet-Hall anomalous signal mean was approximately this value, and like the NTO prediction, had both constant phase (value not noted in the article) relative to the lab and variation at twice the table rotation rate.

## 6 GRAVITATIONAL ORBITS VS TRUE ROTATION

Why should we get a null signal for the solar and galactic orbital velocities, but a non-null signal for the earth surface speed from its own rotation? The answer is that bodies in gravitational orbits follow geodesics, i.e., they are in ”free fall”. That is, they are in locally inertial (time orthogonal) frames and therefore obey Lorentzian mechanics. An object in such a frame feels no force, provided it does not rotate about its own axis. There is no experimental means by which one could measure (without looking outside at the stars) one’s rate of revolution (and hence one’s speed in orbit).

On the other hand, objects (like the Brillet-Hall or Michelson-Morley apparatuses) fixed in ”true” rotational frames are held in place by non-gravitational forces and do not travel geodesic paths. Via experiment, one can actually determine one’s speed in such frames in the circumferential direction (without looking outside) in an absolute (or Machian) sense.

For example, in true rotation one can measure the motion of a Foucault pendulum’s arc over time and determine the rate of rotation $`\omega `$. One can also use a mass-spring system to determine the distance $`r`$ to the center of rotation, and thereby calculate the circumferential speed $`v=\omega r`$ in an absolute sense. In gravitational orbit this can not be done. If the planet does not rotate about its own axis, there is no movement of a pendulum’s arc, even thought the planet travels around its orbit. No centrifugal acceleration is felt, so one can not use a spring-mass system to determine where the center of revolution is. Hence, speed in gravitational orbit is indeterminate via experiment.

Thus, the earth surface speed is that of a true rotational system; orbital speed is not. This has at least two immediate ramifications. First, if for true rotation one can measure circumferential speed via one experiment (pendulum, spring, mass), then why not by another (Michelson-Morley or Brillet-Hall)? Second, the non-time-orthogonal factor (see off diagonal term in (3)) is a direct function of $`\omega `$ and $`r`$. Measuring those, we effectively measure the degree of non-time-orthogonality, and hence the degree of anisotropy.

Hence, the usual form of relativity should hold for gravitational orbits and we should expect null Michelson-Morley and Brillet-Hall results for orbital speeds, which is just what we measure. But the NTO form of relativity theory should be applied for the earth surface speed due to rotation about the earth’s axis. This issue is treated in detail in Klauber.

Note that the anisotropy in light speed is predicted to exist on the earth provided the earth is rotating, though it would be undetectable in practice unless the experimental apparatus turned with respect to the earth’s surface, thereby enabling speed comparisons in different directions (e.g., E-W then N-S) within the rotating earth frame itself. Although other tests have been performed with accuracy comparable to, or better than, that of Brillet-Hall, none of these turned the test apparatus relative to the earth surface.

## 7 SUMMARY AND CONCLUSIONS

Non-time-orthogonal analysis of relativistically rotating frames predicts a signal that is in accord with the anomalous non-null signal on the order of 10-13 found in the Brillet-Hall experiment. Contributions to the signal come from time delay and three separate transverse effects on wave propagation. Both the predicted signal and the measured signal have fixed phase in the lab frame and vary with twice the test apparatus rotation rate. Signal magnitude and phase depend on particular dimensions between components of the test apparatus, which are unknown, though estimable. Reasonable estimates for these dimensions result in a prediction that is well within the range of data found via experiment.

NTO analysis predicts a signal that is dependent on the surface velocity of the earth due to rotation about its axis, but none due to the velocity of the earth in gravitational orbit. This is in agreement with all relevant experiments. This result is due to off diagonal time-space terms that arise in the metric for true rotating frames (where angular velocity and centrifugal acceleration can be measured). No such off diagonal terms exist in the metric for bodies in the free fall frames of gravitational orbit.

This result contrasts with the null prediction of the usual relativistic analysis approach, utilizing local co-moving Lorentz frames. The Brillet-Hall measurement, GPS data, and the analysis herein raise significant questions relating to relativity theory and the nature of the spacetime continuum. (Can time truly be non-orthogonal to space, and if so, what ramifications are there for physical measurements?) The author hopes that someone will soon repeat the Brillet-Hall experiment (with apparatus rotating with respect to the earth surface) and provide definitive answers to these questions.

## APPENDIX A: THE FABRY-PEROT INTERFEROMETER

For reference, this appendix provides a derivation of the well known Fabry-Perot constructive interference condition for thin films666 This article does not delve into certain intricacies of geometrical optics in rotating frames, particularly with regard to reflection and half-reflection, and in this sense, derivations herein are approximate..

Figure 1 depicts the Fabry-Perot interferometer with two partially transmitting glass (silica) mirror sections mounted a distance $`d`$ (the “etalon” length) apart. A light ray enters from the left at angle $`\alpha `$ to the horizontal, is refracted through the first glass section, and then is partially reflected at the surface of the second glass section. The reflected ray is reflected once again at A, then partially refracted at B, and upon exiting the right side of the second glass section interferes with the original refracted ray. Two phase changes of 180o occur (at the first reflection and the second reflection, point B) so the net phase change due to reflection upon exiting the right hand glass section is zero. The reflected portion of the ray at B leads to a third ray, as well as subsequent rays, exiting the RHS. Each successive such ray is diminished in intensity from the prior one, and the interference behavior of all of these rays together parallels that of the first two shown in Figure 1.

Note that for constructive interference, we must have

$$2ABCD=m\lambda ,$$
(38)
where $`m`$ is an integer. The following trigonometry is straightforward.

$$AB=\frac{d}{\mathrm{cos}\alpha }$$
(39)
$$l=d\mathrm{tan}\alpha $$
(40)
$$CD=2l\mathrm{sin}\alpha =2d\mathrm{tan}\alpha \mathrm{sin}\alpha $$
(41)
Using the above values for AB and CD in (38), one finds

$$\frac{2d}{\mathrm{cos}\alpha }2d\mathrm{tan}\alpha \mathrm{sin}\alpha =m\lambda .$$
(42)
Multiplication by cos$`\alpha `$ yields

$$1\mathrm{sin}^2\alpha =\mathrm{cos}\alpha \frac{m\lambda }{2d},$$
(43)
or

$$\frac{m\lambda }{2d}=\mathrm{cos}\alpha .$$
(44)
## APPENDIX B: GALILEAN ROUND TRIP TIMES IN TRANSLATING FRAMES

The general relation for determining measured round trip travel times for different speeds in each direction is

$$t_{RT}=\frac{l}{v_+}+\frac{l}{v_{}}$$
(45)
where $`l`$ is the distance traveled in one direction measured by standard meter sticks, $`v_+`$ is the physical (measured) speed in the outward direction and $`v_{}`$ is the physical return speed.

Using (45), the pre-relativistic ether approach to light travel times in directions parallel and perpendicular to velocity of translation through the ether yielded, to second order,

$$t_{RT,perpend}\frac{2l}{c}\left(1+\frac{1}{2}\frac{v^2}{c^2}\right)$$
(46)
and

$$t_{RT,parallel}\frac{2l}{c}\left(1+\frac{v^2}{c^2}\right).$$
(47)
The difference between these is

$$\mathrm{\Delta }t_{\begin{array}{c}RT,Galilean\hfill \\ movingframe\hfill \end{array}}\frac{2l}{c}\left(\frac{1}{2}\frac{v^2}{c^2}\right)=(t_{RT,nomotion})\left(\frac{1}{2}\frac{v^2}{c^2}\right)=\frac{1}{2}\beta ^2(t_{RT,nomotion})$$
(48)
The fractional difference in round trip light travel time between the two paths equals $`1/2\beta ^2`$.

This value is reflected directly in the predicted (for classical ether) change in interference fringing in the Michelson-Morley experiment and in the change in frequency of the Fabry-Perot interferometer in the Brillet-Hall experiment. This, of course, presumes that motion in each case can be modeled as translational and Galilean.

It is significant that the Michelson-Morley experiment was not sensitive enough to detect any effect from the earth surface speed, though the Brillet-Hall experiment was.

## APPENDIX C: CASE OF INITIAL PEAK INTENSITY AT $`\alpha 0`$

Consider the case where the sensor is not located at the precise centerline of the Fabry-Perot interferometer, such that the local peak in intensity striking the sensor is not at $`\alpha `$ = 0. That is, $`\alpha _00`$, though still small. Then

$$\frac{m\lambda _0}{2d}=\mathrm{cos}\alpha _01\frac{1}{2}\alpha _0^2$$
(49)
yields the wavelength $`\lambda _o`$ for initial peak intensity at the sensor. As the apparatus turns, due to anisotropy, the servo changes the wavelength to maintain constant $`\alpha _o`$. The extremum value for this wavelength is $`\lambda _1`$, and this corresponds to what the peak intensity angle $`\alpha _1`$ would be for such a wavelength if the apparatus had not turned. Thus,

$$\frac{m\lambda _1}{2d}=\mathrm{cos}\alpha _11\frac{1}{2}\alpha _1^2.$$
(50)
From (49) and (50),

$$\frac{\mathrm{\Delta }\lambda }{\lambda _0}=\frac{\lambda _1\lambda _0}{\lambda _0}=\frac{\lambda _1}{\lambda _0}1=\frac{\mathrm{cos}\alpha _1}{\mathrm{cos}\alpha _0}1\frac{1\frac{1}{2}\alpha _1^2}{1\frac{1}{2}\alpha _0^2}1\frac{1}{2}\left(\alpha _1^2\alpha _0^2\right).$$
(51)
For

$$\alpha _1=\alpha _0+\mathrm{\Delta }\alpha ,$$
(52)
(51) becomes

$$\frac{\mathrm{\Delta }\lambda }{\lambda _0}\frac{1}{2}\left((\mathrm{\Delta }\alpha )^2+2(\mathrm{\Delta }\alpha )\alpha _0\right).$$
(53)
The first term on the right side will cause a second order variation in $`\mathrm{\Delta }\lambda `$ at twice the table rotation rate. (See (26) to (28).) The last term should cause a variation at the table rotation rate with a magnitude dependent on the initial condition $`\alpha _o`$.

In this context, it may be significant that Brillet-Hall detected another non-null signal of magnitude $``$1.2X10-12 at the table rotation rate. They attributed this signal to gravitational stretching of their apparatus, but in light of (53), one might wonder if the cause lies elsewhere. A repetition of their test would tell us a lot.