The broad band spectral index of TeV blazars detected by Fermi LAT

Using γ-ray data detected by Fermi Large Area Telescope (LAT) and multi-wave band data for 35 TeV blazars sample, we have studied the possible correlations between different broad band spectral indices ( $\alpha_{\rm r.ir}$ , $\alpha_{\rm{r.o}}$ , $\alpha_{\rm r.x}$ , $\alpha_{\rm r.\gamma}$ , $\alpha_{\rm{ir.o}}$ , $\alpha_{\rm ir.x}$ , $\alpha_{\rm ir.\gamma}$ , $\alpha_{\rm o.x}$ , $\alpha_{\rm o.\gamma}$ , $\alpha_{\rm r.x}$ , $\alpha_{\rm x.\gamma}$ ) in all states (average/high/low). Our results are as follows: (1) For our TeV blazars sample, the strong positive correlations were found between $\alpha_{\rm r.ir}$ and $\alpha_{\rm{r.o}}$ , between $\alpha_{\rm r.ir}$ and $\alpha_{\rm r.x}$ , between $\alpha_{\rm r.ir}$ and $\alpha_{\rm r.\gamma}$ in all states (average/high/low); (2) For our TeV blazars sample, the strong anti-correlations were found between $\alpha_{\rm r.ir}$ and $\alpha_{\rm x.\gamma}$ , between $\alpha_{\rm{r.o}}$ and $\alpha_{\rm ir.\gamma}$ , between $\alpha_{\rm{r.o}}$ and $\alpha_{\rm o.\gamma}$ , between $\alpha_{\rm{r.o}}$ and $\alpha_{\rm x.\gamma}$ , between $\alpha_{\mathrm{ir.o}}$ and $\alpha_{\rm o.\gamma}$ , between $\alpha_{\rm r.x}$ and $\alpha_{\rm x.\gamma}$ , between $\alpha_{\rm ir.x}$ and $\alpha_{\rm x.\gamma}$ in all states (average/high/low). The results suggest that the synchrotron self-Compton radiation (SSC) is the main mechanism of high energy γ-ray emission and the inverse Compton scattering of circum-nuclear dust is likely to be a important complementary mechanism for TeV blazars. Our results also show that the possible correlations vary from state to state in the same pair of indices, Which suggest that there may exist differences in the emitting process and in the location of the emitting region for different states.     

The planetary spin and rotation period: a modern approach

Using a new approach, we have obtained a formula for calculating the rotation period and radius of planets. In the ordinary gravitomagnetism the gravitational spin (S) orbit (L) coupling, $\vec{L}\cdot\vec{S}\propto L^{2}$ , while our model predicts that $\vec{L}\cdot\vec{S}\propto\frac{m}{M}L^{2}$ , where M and m are the central and orbiting masses, respectively. Hence, planets during their evolution exchange L and S until they reach a final stability at which MS∝mL, or $S\propto\frac{m^{2}}{v}$ , where v is the orbital velocity of the planet. Rotational properties of our planetary system and exoplanets are in agreement with our predictions. The radius (R) and rotational period (D) of tidally locked planet at a distance a from its star, are related by, $D^{2}\propto\sqrt{\frac{M}{m^{3}}}R^{3}$ and that $R\propto\sqrt{\frac {m}{M}}a$ .     

Gaussian variational equations for osculating elements of an arbitrary separable reference orbit

If a satellite orbit is described by means of osculating Jacobi α's and β's of a separable problem, the paper shows that a perturbing forceF makes them vary according to $$\dot \alpha _\kappa = {\text{F}} \cdot \partial {\text{r/}}\partial \beta _k {\text{ }}\dot \beta _k = {\text{ - F}} \cdot \partial {\text{r/}}\partial \alpha _k ,{\text{ (}}k = 1,2,3).{\text{ (A1)}}$$ Herer is the position vector of the satellite andF is any perturbing force, conservative or non-conservative. There are two special cases of (A1) that have been previously derived rigorously. If the reference orbit is Keplerian, equations equivalent to (A1), withF arbitrary, were derived by Brouwer and Clemence (1961), by Danby (1962), and by Battin (1964). IfF=?gradV ₁(t), whereV ₁ may or may not depend explicitly on the time, Equations (A1) reduce to the well known forms (e.g. Garfinkel, 1966) $$\dot \alpha _\kappa = {\text{ - }}\partial V_1 {\text{/}}\partial \beta _k {\text{ }}\dot \beta _k = \partial V_1 {\text{/}}\partial \alpha _k ,{\text{ (}}k = 1,2,3).{\text{ (A2)}}$$ holding for all separable reference orbits. Equations (A1) can of course be guessed from Equations (A2), if one assumes that \(\dot \alpha _k (t)\) and \(\dot \beta _k (t)\) depend only onF(t) and thatF(t) can always be modeled instantaneously as a potential gradient. The main point of the present paper is the rigorous derivation of (A1), without resort to any such modeling procedure. Applications to the Keplerian and spheroidal reference orbits are indicated.     

The Chandra X-ray galaxy clusters at z<1.4: constraints on the evolution of L X −T−M g relations

We analyzed the luminosity-temperature-mass of gas (L _X ?T?M _g ) relations for a sample of 21 Chandra galaxy clusters. We used the standard approach (β?model) to evaluate these relations for our sample that differs from other catalogues since it considers galaxy clusters at higher redshifts (0.4<z<1.4). We assumed power-law relations in the form $L_{X} \sim(1 +z)^{A_{L_{X}T}} T^{\beta_{L_{X}T}}$ , $M_{g} \sim(1 + z)^{A_{M_{g}T}} T^{\beta_{M_{g}T}}$ , and $M_{g} \sim(1 + z)^{A_{M_{g}L_{X}}} L^{\beta_{M_{g}L_{X}}}$ . We obtained the following fitting parameters with 68 % confidence level: $A_{L_{X}T} = 1.50 \pm0.23$ , $\beta_{L_{X}T} = 2.55 \pm0.07$ ; $A_{M_{g}T} = -0.58 \pm0.13$ and $\beta_{M_{g}T} = 1.77 \pm0.16$ ; $A_{M_{g}L_{X}} \approx-1.86 \pm0.34$ and $\beta_{M_{g}L_{X}} = 0.73 \pm0.15$ , respectively. We found that the evolution of the M _g ?T relation is small, while the M _g ?L _X relation is strong for the cosmological parameters Ω _m =0.27 and Ω _Λ =0.73. In overall, the clusters at high-z have stronger dependencies between L _X ?T?M _g correlations, than those for clusters at low-z. For most of galaxy clusters (first of all, from MACS and RCS surveys) these results are obtained for the first time.     

On the Maximum Separation of Visual Binaries

In this paper, an efficient algorithm is established for computing the maximum (minimum) angular separation ρ _max(ρ _min), the corresponding apparent position angles ( $\theta|_{\rho_{\rm max}}$ , $\theta|_{\rho_{\rm min}}$ ) and the individual masses of visual binary systems. The algorithm uses Reed’s formulae (1984) for the masses, and a technique of one-dimensional unconstrained minimization, together with the solution of Kepler’s equation for $(\rho_{\rm max}, \theta|_{\rho_{\rm max}})$ and $(\rho_{\rm min}, \theta|_{\rho_{\rm min}})$ . Iterative schemes of quadratic coverage up to any positive integer order are developed for the solution of Kepler’s equation. A sample of 110 systems is selected from the Sixth Catalog of Orbits (Hartkopf et al. 2001). Numerical studies are included and some important results are as follows: (1) there is no dependence between ρ _max and the spectral type and (2) a minor modification of Giannuzzi’s (1989) formula for the upper limits of ρ _max functions of spectral type of the primary.     

Symmetric and regularized coordinates on the plane triple collision manifold

The planar problem of three bodies is described by means of Murnaghan's symmetric variables (the sidesa _j of the triangle and an ignorable angle), which directly allow for the elimination of the nodes. Then Lemaitre's regularized variables \(\alpha _j = \sqrt {(\alpha ^2 - \alpha _j )}\) , where \(\alpha ^2 = \tfrac{1}{2}(a_1 + a_2 + a_3 )\) , as well as their canonically conjugated momenta are introduced. By finally applying McGehee's scaling transformation \(\alpha _j = r^{{1 \mathord{\left/ {\vphantom {1 2}} \right. \kern-0em} 2}} \tilde \alpha _j\) , wherer ² is the moment of inertia a system of 7 differential equations (with 2 first integrals) for the 5-dimensional triple collision manifold \(T\) is obtained. Moreover, the zero angular momentum solutions form a 4-dimensional invariant submanifold \(N \subset T\) represented by 6 differential equations with polynomial right-hand sides. The manifold \(N\) is of the topological typeS ²×S ² with 12 points removed, and it contains all 5 restpoint (each one in 8 copies). The flow on \(T\) is gradient-like with a Lyapounov function stationary in the 40 restpoints. These variables are well suited for numerical studies of planar triple collision.     

Nonlinear newtonian theory and the effect of time dilatation

The fact that the energy density ρ_g of a static spherically symmetric gravitational field acts as a source of gravity, gives us a harmonic function \(f\left( \varphi \right) = e^{\varphi /c^2 } \) , which is determined by the nonlinear differential equation $$\nabla ^2 \varphi = 4\pi k\rho _g = - \frac{1}{{c^2 }}\left( {\nabla \varphi } \right)^2 $$ Furthermore, we formulate the infinitesimal time-interval between a couple of events measured by two different inertial observers, one in a position with potential φ-i.e., dt _φ and the other in a position with potential φ=0-i.e., dt ₀, as $${\text{d}}t_\varphi = f{\text{d}}t_0 .$$ When the principle of equivalence is satisfied, we obtain the well-known effect of time dilatation.     

Stability of the triangular Lagrange points beyond Gascheau��s value

We examine the stability of the triangular Lagrange points L ₄ and L ₅ for secondary masses larger than the Gascheau??s value ${\mu_{\rm G}= (1-\sqrt{23/27}/2)= 0.0385208\ldots}$ (also known as the Routh value) in the restricted, planar circular three-body problem. Above that limit the triangular Lagrange points are linearly unstable. Here we show that between??? _G and ${\mu \approx 0.039}$ , the L ₄ and L ₅ points are globally stable in the sense that a particle released at those points at zero velocity (in the corotating frame) remains in the vicinity of those points for an indefinite time. We also show that there exists a family of stable periodic orbits surrounding L ₄ or L ₅ for ${\mu \ge \mu_G}$ . We show that??? _G is actually the first value of a series ${\mu_0 (=\mu_G), \mu_1,\ldots, \mu_i,\ldots}$ corresponding to successive period doublings of the orbits, which exhibit ${1, 2, \ldots, 2^i,\ldots}$ cycles around L ₄ or L ₅. Those orbits follow a Feigenbaum cascade leading to disappearance into chaos at a value ${\mu_\infty = 0.0463004\ldots}$ which generalizes Gascheau??s work.     

Collisionless tearing modes in the presence of shear flow

Tearing modes in a plane collisionless current sheet with shear bulk flow are studied. An analytic expression for the growth rate is obtained for the case \(M^2 = (1 - \varepsilon {\text{ sech}}^m \bar z)\) , whereM is the Mach number,m the shear flow index, ε a positive constant less than unity, and \(\bar z\) the (normalized) co-ordinate normal to the current sheet. The growth rates are large and the unstable wave number domain is increased as compared to the case without flow. The relevance of these results to time-dependent reconnection processes in the Earth's magnetosphere is discussed.     

Perspectives on effectively constraining the location of a massive trans-Plutonian object with the New Horizons spacecraft: a sensitivity analysis

The radio tracking apparatus of the New Horizons spacecraft, currently traveling to the Pluto system where its arrival is scheduled for July 2015, should be able to reach an accuracy of 10 m (range) and 0.1  $\text{ mm } \text{ s }^{-1}$ mm s ? 1 (range-rate) over distances up to 50 au. This should allow to effectively constrain the location of a putative trans-Plutonian massive object, dubbed Planet X (PX) hereafter, whose existence has recently been postulated for a variety of reasons connected with, e.g., the architecture of the Kuiper belt and the cometary flux from the Oort cloud. Traditional scenarios involve a rock-ice planetoid with $m_\mathrm{X}\approx 0.7\,m_{\oplus }$ m X ≈ 0.7 m ⊕ at some 100–200 au, or a Jovian body with $m_\mathrm{X}\lesssim 5\,m_\mathrm{J}$ m X ? 5 m J at about 10,000–20,000 au; as a result of our preliminary sensitivity analysis, they should be detectable by New Horizons since they would impact its range at a km level or so over a time span 6 years long. Conversely, range residuals statistically compatible with zero having an amplitude of 10 m would imply that PX, if it exists, could not be located at less than about 4,500 au ( $m_\mathrm{X}=0.7\,m_{\oplus }$ m X = 0.7 m ⊕ ) or 60,000 au ( $m_\mathrm{X}=5\,m_\mathrm{J}$ m X = 5 m J ), thus making a direct detection quite demanding with the present-day technologies. As a consequence, it would be appropriate to rename such a remote body as Thelisto. Also fundamental physics would benefit from this analysis since certain subtle effects predicted by MOND for the deep Newtonian regions of our Solar System are just equivalent to those of a distant pointlike mass.     

A note on theH-functions of transfer problems in multiplying media

Some useful results and remodelled representations ofH-functions corresponding to the dispersion function $$T\left( z \right) = 1 - 2z^2 \sum\limits_1^n {\int_0^{\lambda r} {Y_r } \left( x \right){\text{d}}x/\left( {z^2 - x^2 } \right)} $$ are derived, suitable to the case of a multiplying medium characterized by $$\gamma _0 = \sum\limits_1^n {\int_0^{\lambda r} {Y_r } \left( x \right){\text{d}}x > \tfrac{1}{2} \Rightarrow \xi = 1 - 2\gamma _0< 0} $$      

Cosmic-ray antimatter: A primary origin hypothesis

We examine the possibility that the observed cosmic-ray protons are of primary extragalactic origin. The present \(\bar p\) data are consistent with a primary extragalactic component having \(\bar p\) /p?3.2±0.7 x 10^-4 independent of energy. Following the suggestion that most extragalactic cosmic rays are from active galaxies, we propose that most of the observed \(\bar p\) 's are alos from the same sites. This would imply the possibility of destroying the corresponding \(\bar \alpha \) 'sat the source, thus leading to a flux ratio \(\bar \alpha \) /α< \(\bar p\) /p. We further predict an estimate for \(\bar \alpha \) α～10^-5, within the range of future cosmic-ray detectors. the cosmological implications of this proposal are discussed.     

On pulsar-driven isothermal blast-wave models of supernova remanants

We investigate the ‘equilibrium’ and stability of spherically-symmetric self-similar isothermal blast waves with a continuous post-shock flow velocity expanding into medium whose density varies asr ^?ω ahead of the blast wave, and which are powered by a central source (a pulsar) whose power output varies with time ast ^ω?3. We show that:

1. for ω<0, no physically acceptable self-similar solution exists;
2. for ω>3, no solution exists since the mass swept up by the blast wave is infinite;
3. ? must exceed zero in order that the blast wave expand with time, but ?<2 in order that the central source injects a finite total energy into the blast wave;
4. for 3>ω_min(?)>ω>ω_max(?)>0, where $$\begin{gathered} \omega _{\min } (\varphi ){\text{ }} = {\text{ }}2[5{\text{ }} - {\text{ }}\varphi {\text{ }} + {\text{ }}(10{\text{ }} + {\text{ 4}}\varphi {\text{ }} - {\text{ 2}}\varphi ^2 )^{1/2} ]^2 [2{\text{ }} + {\text{ (10 }} + {\text{ 4}}\varphi {\text{ }} - {\text{ 2}}\varphi ^2 {\text{)}}^{{\text{1/2}}} ]^{ - 2} , \hfill \\ \omega _{\max } (\varphi ){\text{ }} = {\text{ }}2[5{\text{ }} - {\text{ }}\varphi {\text{ }} - {\text{ }}(10{\text{ }} + {\text{ 4}}\varphi {\text{ }} - {\text{ 2}}\varphi ^2 )^{1/2} ]^2 [2{\text{ }} - {\text{ (10 }} + {\text{ 4}}\varphi {\text{ }} - {\text{ 2}}\varphi ^2 {\text{)}}^{{\text{1/2}}} ]^{ - 2} , \hfill \\ \end{gathered} $$ two critical points exist in the flow velocity versus position plane. The physically acceptable solution must pass through the origin with zero flow speed and through the blast wave. It must also pass throughboth critical points if \(\varphi > \tfrac{5}{3}\) , while if \(\varphi< \tfrac{5}{3}\) it must by-pass both critical points. It is shown that such a solution exists but a proper connection at the lower critical point (for ?>5/3) (through whichall solutions pass with thesame slope) has not been established;
5. for 3>ω>ω_min(?) it is shown that the two critical points of (iv) disappear. However a new pair of critical points form. The physically acceptable solution passing with zero flow velocity through the origin and also passing through the blast wave mustby-pass both of the new critical points. It is shown that the solution does indeed do so;
6. for 3>ω_min(?)>ω_max(?)>ω it is shown that the dependence of the self-similar solution on either ω or ? is non-analytic and therefore, inferences drawn from any solutions obtained in ω>ω_max(?) (where the dependence of the solutionis analytic on ω and ?) are not valid when carried over into the domain 3>ω_min(?)>ω_max(?)>ω;
7. all of the physically acceptable self-similar solutions obtained in 3>ω>0 are unstable to short wavelength, small amplitude but nonself-similar radial velocity perturbations near the origin, with a growth which is a power law in time;
8. the physical self-similar solutions are globally unstable in a fully nonlinear sense to radial time-dependent flow patterns. In the limit of long times, the nonlinear growth is a power law in time for 5<ω+2?, logarithmic in time for 5>ω+2?, and the square of the logarithm in time for 5=ω+2?.

The results of (vii) and (viii) imply that the memory of the system to initial and boundary values does not decay as time progresses and so the system does not tend to a self-similar form. These results strongly suggest that the evolution of supernova remnants is not according to the self-similar form.     

On the connection between the Nekhoroshev theorem and Arnold diffusion

The analytical techniques of the Nekhoroshev theorem are used to provide estimates on the coefficient of Arnold diffusion along a particular resonance in the Hamiltonian model of Froeschlé et al. (Science 289:2108–2110, 2000). A resonant normal form is constructed by a computer program and the size of its remainder ||R _opt || at the optimal order of normalization is calculated as a function of the small parameter ${\epsilon}$ . We find that the diffusion coefficient scales as ${D \propto ||R_{opt}||^3}$ , while the size of the optimal remainder scales as ${||R_{opt}|| \propto {\rm exp}(1/\epsilon^{0.21})}$ in the range ${10^{-4} \leq \epsilon \leq 10^{-2}}$ . A comparison is made with the numerical results of Lega et al. (Physica D 182:179–187, 2003) in the same model.     

Apsidal motion in the close binary IT cassiopeiae

In 1982 and 1993, we carried out highly accurate photoelectric WBVR measurements for the close binary IT Cas. Based on these measurements and on the observations of other authors, we determined the apsidal motion $\left[ {\dot \omega _{obs} = {{(11\mathop .\limits^ \circ 0 \pm 2\mathop .\limits^ \circ 5)} \mathord{\left/ {\vphantom {{(11\mathop .\limits^ \circ 0 \pm 2\mathop .\limits^ \circ 5)} {100 years}}} \right. \kern-0em} {100 years}}} \right]$ . This value is in agreement with the theoretically calculated apsidal motion for these stars $\left[ {\dot \omega _{th} = {{(14^\circ \pm 3^\circ )} \mathord{\left/ {\vphantom {{(14^\circ \pm 3^\circ )} {100 years}}} \right. \kern-0em} {100 years}}} \right]$ .     

Triple approaches in the plane isosceles equal-mass three-body problem

We analyze flyby-type triple approaches in the plane isosceles equal-mass three-body problem and in its vicinity. At the initial time, the central body lies on a straight line between the other two bodies. Triple approaches are described by two parameters: virial coefficient k and parameter $\mu = \dot r/\sqrt {\dot r^2 + \dot R^2 }$ , where $\dot r$ is the relative velocity of the extreme bodies and $\dot R$ is the velocity of the central body relative to the center of mass of the extreme bodies. The evolution of the triple system is traceable until the first turn or escape of the central body. The ejection length increases with closeness of the triple approach (parameter k). The longest ejections and escapes occur when the extreme bodies move apart with a low velocity at the time of triple approach. We determined the domain of escapes; it corresponds to close triple approaches (k>0.8) and to μ in the range ?0.2<μ<0.7. For small deviations from the isosceles problem, the evolution does not differ qualitatively from the isosceles case. The domain of escapes decreases with increasing deviations. In general, the ejection length increases for wide approaches and decreases for close approaches.     

A new look at the origin of the 6.67 hr period X-ray pulsar 1E 161348-5055

The point X-ray source 1E 161348-5055 is observed to display pulsations with the period 6.67?hr and $|\dot{P}| \leq1.6 \times10^{-9}\,{\rm s\,s^{-1}}$ . It is associated with the supernova remnant RCW?103 and is widely believed to be a ～2000?yr old neutron star. Observations give no evidence for the star to be a member of a binary system. Nevertheless, it resembles an accretion-powered pulsar with the magnetospheric radius ～3000?km and the mass-accretion rate $\sim 10^{14}\,{\rm g\,s^{-1}}$ . This situation could be described in terms of accretion from a (residual) fossil disk established from the material falling back towards the star after its birth. However, current fall-back accretion scenarios encounter major difficulties explaining an extremely long spin period of the young neutron star. We show that the problems can be avoided if the accreting material is magnetized. The star in this case is surrounded by a fossil magnetic slab in which the material is confined by the magnetic field of the accretion flow itself. We find that the surface magnetic field of the neutron star within this scenario is ～10¹²?G and that a presence of $\gtrsim10^{-7}\,{\rm M_{\odot}}$ magnetic slab would be sufficient to explain the origin and current state of the pulsar.     

Non-minimal derivatively coupled quintessence in the Palatini formalism

The quintessence dark energy model with a kinetic coupling to gravity within the Palatini formalism is studied in this paper. Two different coupling forms: $\hat{R}\partial^{\mu}\phi\partial_{\mu}\phi$ and $\hat {R}_{\mu\nu}\partial^{\mu}\phi\partial^{\nu}\phi$ are analyzed, respectively. We find that both the model with the $\hat{R}\partial^{\mu}\phi\partial_{\mu}\phi$ coupling and the one with the $\hat{R}_{\mu\nu}\partial^{\mu}\phi\partial^{\nu}\phi$ coupling can realize the phantom divide line crossing from phantom to quintessence at late time for its effective equation-of-state. Furthermore, the former can behave like phantom. These features are different from those found in the $\hat {R}\phi^{2}$ coupling case.     

A Statistical Study on Photospheric Magnetic Nonpotentiality of Active Regions and Its Relationship with Flares During Solar Cycles 22?–?23

A statistical study is carried out on the photospheric magnetic nonpotentiality in solar active regions and its relationship with associated flares. We select 2173 photospheric vector magnetograms from 1106 active regions observed by the Solar Magnetic Field Telescope at Huairou Solar Observing Station, National Astronomical Observatories of China, in the period of 1988??C?2008, which covers most of the 22nd and 23rd solar cycles. We have computed the mean planar magnetic shear angle ( $\overline{\Delta\phi}$ ), mean shear angle of the vector magnetic field ( $\overline{\Delta\psi}$ ), mean absolute vertical current density ( $\overline{|J_{z}|}$ ), mean absolute current helicity density ( $\overline{|h_{\mathrm{c}}|}$ ), absolute twist parameter (|?? _av|), mean free magnetic energy density ( $\overline{\rho_{\mathrm{free}}}$ ), effective distance of the longitudinal magnetic field (d _E), and modified effective distance (d _Em) of each photospheric vector magnetogram. Parameters $\overline{|h_{\mathrm{c}}|}$ , $\overline{\rho_{\mathrm{free}}}$ , and d _Em show higher correlations with the evolution of the solar cycle. The Pearson linear correlation coefficients between these three parameters and the yearly mean sunspot number are all larger than 0.59. Parameters $\overline {\Delta\phi}$ , $\overline{\Delta\psi}$ , $\overline{|J_{z}|}$ , |?? _av|, and d _E show only weak correlations with the solar cycle, though the nonpotentiality and the complexity of active regions are greater in the activity maximum periods than in the minimum periods. All of the eight parameters show positive correlations with the flare productivity of active regions, and the combination of different nonpotentiality parameters may be effective in predicting the flaring probability of active regions.     

Etude variationnelle des decalages spectraux

In a static gravitational field the paths of light are curved, as noticed by H. Weyl. This property can bea priori stated for aV ₃ Riemannian manifold: through any two points ofV ₃ it is possible to draw two families of curves, the straight lines of Euclidean geometry and the photon trajectoriesz. We can perform a fibration of the Galilean space-time in an original way, by taking thez-trajectories of the photons as the base, the isochronic surfaces as fibres, and ‘the equal length time on az trajectory to reach a given point’ as the equivalence relation. The straight lines of Euclidean geometry can then carry the classical mechanics timet, and thez trajectories can carry the optics time t. These times are related by dt=F(x,t) dt. If we class the Universe as a pseudo-Riemannian manifold of normal hyperbolic typeC ^∞, the time t determined above can be taken as the time coordinate inV ₄. Under these conditions we have \(d\overline s ^2 \) =F ² \(d\overline s ^2 \) , where \(d\overline s ^2 \) is the metric of the Riemannian manifold, conforming to the metric ds ² and allowing t as the cosmic time. We can then use the results previously achieved by the author (Peton, 1979) and write: 1 +Z _G =F(A _s,t _s,)/F(Ao_s,t _o) wherez _G denotes the shift of the spectral lines due to the metric. In the case of relative motion betweenO andS, we have $${\text{1 + z' = (1 + }}z_{\text{G}} {\text{)(1 + }}\beta _{\text{r}} {\text{)(1 }} - {\text{ }}\beta ^2 {\text{)}}^{ - 1/2} $$ The Doppler-Fizeau effect therefore appears as a result of the application of the Fermat principle.