Stellar aberration (derivation from Lorentz transformation)

Stellar aberration is an astronomical phenomenon "which produces an apparent motion of celestial objects". It can be proven mathematically that stellar aberration is due to the change of the astronomer's inertial frame of reference. The formula is derived with the use of Lorentz transformation of the star's coordinates.

As the astronomer John Herschel has already explained in 1844, the stellar aberration does not depend on the relative velocity of the star towards Earth.^[1] Otherwise eclipsing binary stars would appear to be separated, in stark contrast to observation: both stars are rotating with high speed —and ever changing and different velocity vectors— around each other, but they appear as one spot all the time.

Stellar aberration is only due to the change of the astronomer's inertial frame of reference

In the year 1926 the astrophysicist Robert Emden published the article Aberration und Relativitätstheorie in the journal Naturwissenschaften.^[2] In this article he states that the direction of a light ray isn't influenced by the motion of the star or by the motion of Earth.^{[Notes 1]} At that time, the opponents of the special theory of relativity reasoned that the theory must be flawed, because it would state that the stellar aberration would depend on the relative velocity of the star — which would be in contradiction to observation — and R. Emden's article explains that the special theory of relativity does not predict this. Today, the special theory of relativity isn't contested anymore but there are still articles that suggest that the aberration would depend on the relative velocity of the star.^[3]

Although a (relativistic) velocity-addition formula can be used to explain stellar aberration, (see Aberration of light), another (relativistic) explanation using only the Lorentz transformation is also possible, as will be demonstrated. This derivation only uses the star's coordinates at the time of emission, and therefore has the formal advantage there is no place for the relative velocity of the star towards the astronomer and therefore it is evident that the observed position doesn't depend on the star's velocity — provided that the resultant change of position is much smaller than the distance between star and Earth.^{[Notes 2]} The observed position of the star wouldn't depend on Earth's motion either, if the astronomer could use the same inertial frame of reference all the time. But of course that is technically impossible,^{[Notes 3]} the astronomer uses his current rest frame and these current rest frames are different at different times as Earth orbits around the Sun. It is mathematically convenient to declare the position of star in a rest frame of the Sun (more exactly: the center of mass of the Solar System) as the "real" position and that the difference to this "real" position derives form the "aberration".^{[Notes 4]}

Sample calculation

path of light signal in reference frame S

path of light signal in reference frame S'

S and S' are (quasi-)inertial frames of reference and the reference frame S' is in uniform motion with v_x = 0,5c relative to S, such that in the future the star comes nearer to the origin of the coordinate system (and consequently further afar in the past). The x-,y- and z-axes of both systems ought to be parallel and at time t=t'=0 the origins of both systems ought to coincide. Therefore one gets according to Lorentz transformation: $\scriptstyle \beta=v/c$, $\scriptstyle \gamma=1/\sqrt{1-\beta^{2}}$: y'=y ; z=z'; $\scriptstyle x'=\gamma\cdot(x-\beta\cdot ct)$

Suppose now that the star emitted a light signal at time $\scriptstyle c\,t_1=-5\,Ly$ at location $\scriptstyle x_1=4\,Ly\;,\;y_1= 3\,Ly\;,\;z_1=0$ (S coordinates) and that this light signal is received by an astronomer at time $\scriptstyle c\,t_2=0$ at location $\scriptstyle x_2=0\;,\;y_2=0\;,\;z_2=0$.

In S the star's position and the x-axis form an angle $\delta$ with $\scriptstyle \tan \delta =3/4 \; \rightarrow \;\delta =36.87^\circ$ But in S' $\scriptstyle x_1'=\gamma \cdot (x_1-\beta \cdot c\,t_1)=1,1547\cdot (4 Lj - 0,5\cdot(-5 Lj))=7,50555 Lj; \quad y_1'=y_1=3 Lj$ and therefore the angle $\delta'$ between the star's position and x'-axis^{[Notes 5]} is $\scriptstyle \tan \delta'=y_1'/x_1'=3/7,50555\;\rightarrow\;\delta' = 21,79^\circ$

Calculation with help of the formula in aberration of light#Explanation gives the same result: $\scriptstyle \tan (\delta'/2) = \tan (\delta /2) \cdot \sqrt {(1-0,5)/(1+0,5)} = \tan (36,87^\circ/2) \cdot \sqrt{1/3} = 0,19245\;$ $\scriptstyle\;\rightarrow\;\delta'/2= 10,89^\circ\;\rightarrow\;\delta'=21,79^\circ$.

Two-dimensional problem

Derivation of the formula for motion along the x-axis

For the derivation it is assumed, that the light signal only travels through space regions where the gravitation field is negligible. Hence is suffices to use special relativity and the path of the light signal is a straight line in any inertial frame of reference.

Observation in the rest frame S of the center of mass of our Solar System

The rest frame of the center of mass (barycenter) is a very good^{[Notes 6]} quasi-inertial frame of reference for periods of time in the order of thousands of years, since our solar system needs about 230 million years (galactic year) to move completely around the center of the Milky Way. The space coordinates of this frame of reference form a Cartesian coordinate system.

In the reference frame S $(x_1|y_1|0|c\,t_1)$ with $c\,t_1<0$ and $(x_1|y_1) \ne (0|0)$ are the (space-time) coordinates at which the star emits a light signal and $(0|0|0|0)$ are the coordinates at which the astronomer receives the light signal.

In reference frame S the light signal starts at $t_1<0$ and stops at time $t_2=0$ and therefore the light signal did cover the distance $d = c\cdot(0-t_1) = -c\,t_1$ .

In S the path of the light signal is a straight line and it forms an angle $\delta$ with the x-axis with: $\sin \delta = \frac{y_1}{d}=\frac{y_1}{-c\,t_1}$  and $\cos\delta=\frac{x_1}{-c\,t_1}$

Observation in the inertial frame of reference S' which is in uniform motion (relative to S) along the x-axis

The origin of the reference frame S' is in uniform motion relative to S with $(v|0|0)$, i.e. moves along the x-axis, and the x-,y- und z-axes of S' and S are parallel to each other and at time $\scriptstyle t=t'=0$ the origins of S and S' coincide. Let $\beta = \frac{v}{c} , \gamma = \frac{1}{\sqrt{1-\beta^2}}$

S' now is an equally good quasi-inertial frame of reference as S: the space coordinates form a Cartesian coordinate system and the path of the light signal is a straight line.

According to the Lorentz transformation one gets: $x_1' = \gamma \cdot(x_1-\beta \cdot c\,t_1) , y_1'=y_1 , z_1'=0 , c\,t_1' = \gamma \cdot (c\,t_1-\beta \cdot x_1)$

In reference frame S' the light signal starts at $t_1'<0$ and stops at time $t_2'=0$ and therefore the light signal did cover the distance $d\,' = c\cdot(0-t_1') = -c\,t_1'$ .

In S' the path of the light signal is a straight line, too. It forms an angle $\delta'$ with the x'-axis and one gets:

$\sin \delta' = \frac{y_1'}{d\,'}=\frac{y_1}{-c\,t_1'} =\frac{y_1}{-\gamma \cdot (c\,t_1-\beta \cdot x_1)} =\frac{y_1}{-c\,t_1 \cdot \gamma\cdot \left(1-\beta\cdot\frac{x_1}{c\,t_1}\right)}=\frac{\frac{y_1}{-c\,t_1}}{\gamma\cdot\left(1+\beta\cdot\frac{x_1}{-c\,t_1}\right)} =\frac{\sin \delta}{\gamma\cdot(1+\beta\cdot\cos\delta)}$

$\cos\delta'=\frac{x_1'}{-c\,t_1'} =\frac{\gamma \cdot(x_1-\beta \cdot c\,t_1)}{-\gamma \cdot (c\,t_1-\beta \cdot x_1)} =\frac{-c\,t_1 \cdot \gamma \cdot \left ( \frac{x_1}{-c\,t_1}+\beta \right)}{-c\,t_1 \cdot \gamma \cdot \left( 1 + \beta \cdot \frac{x_1}{-c\, t_1} \right) } =\frac{\frac{x_1}{-c\,t_1}+\beta}{1 + \beta \cdot \frac{x_1}{-c\, t_1}} =\frac{\cos \delta +\beta}{1 + \beta \cdot \cos \delta}$

Hence: $\tan \delta' = \frac{\sin \delta'}{\cos \delta'}=\frac{\sin \delta}{\gamma\cdot(\cos \delta +\beta)}$

These are the same formulas as in aberration of light#Explanation.

Approximate formula for motion along the x-axis in case of v/c <<1

$\triangle\delta\;=\;\delta' - \delta\;$ is the change of the angle δ. As β<<1 this change is also very small.

Case I: $\;\delta\ne\pm 90^\circ\; \rightarrow \;\cos \delta \ne 0$

As Δδ<<1 one gets: $\frac{\tan \delta' - \tan \delta}{\triangle\delta} \approx\frac{d}{d\delta}\tan \delta = \frac{1}{(\cos \delta)^2}\; \rightarrow \; \triangle\delta\;=\;(\cos \delta)^2 \cdot (\tan \delta' - \tan \delta)$

As β<<1 one gets: $\tan\delta'=\frac {\sin \delta}{\gamma(\cos \delta +\beta)} \approx \frac {\sin \delta}{\cos \delta +\beta} = \frac{\sin \delta}{\cos \delta \cdot \left(1+\frac{\beta}{\cos \delta}\right)} \approx \tan \delta \cdot \left(1-\frac{\beta}{\cos \delta}\right)$

Therefore $\triangle\delta\;=\;(\cos \delta)^2 \cdot (\tan \delta' - \tan \delta) = (\cos \delta)^2 \cdot \tan \delta \left(1-\frac{\beta}{\cos \delta}\;-1\right) = - \cos \delta \cdot \tan \delta \cdot \beta = -\beta \cdot \sin \delta$

Case IIa: $\;\delta\;=90^\circ\;$, hence: $\tan \delta' =\frac{\sin 90^\circ}{\gamma\cdot(\cos 90^\circ + \beta)}= \frac{1}{\gamma \cdot \beta}\;\approx \frac {1}{\beta} \quad \rightarrow \quad \cot \delta' = \beta \quad \rightarrow \quad \delta' = \arccot \beta \approx \frac{\pi}{2} - \beta$

and therefore: $\triangle\delta = \delta' - \delta = -\beta \quad \left(= -\beta \cdot \sin(90^\circ)\right)$

Case IIb: $\;\delta\;=-90^\circ\;$, and hence: $\tan \delta' = \frac{-1}{\gamma \cdot \beta}\;\approx -\frac {1}{\beta} \quad \rightarrow \quad \cot \delta' = -\beta \quad \rightarrow \quad \delta' = \arccot \beta \approx -\frac{\pi}{2} + \beta \quad$    ^{[Notes 7]}

Hence: $\triangle\delta = \delta' - \delta = +\beta\quad \left(= -\beta \cdot \sin(-90^\circ)\right)$

Conclusion:The change of the angle Δδ = δ'-δ in the case of β = v/c << 1 can be described by the approximate formula $\triangle\delta = -\frac{v}{c} \cdot \sin \delta\;$ resp. in the degree measure $\triangle\delta = -\frac{v}{c} \cdot \sin \delta \cdot \frac{180^\circ}{\pi}$

Symmetric form of the (exact) formula for motion along the x-axis

With help of tangent half-angle formula $\scriptstyle \tan (\alpha/2)=\sin\alpha/(1+\cos\alpha)$  one can prove the symmetric form: $\tan \frac{\delta'}{2}=\sqrt{\frac{1-\beta}{1+\beta}}\cdot\tan\frac{\delta}{2}$   (derivation found in a SR-textbook^[4])

$\tan \frac{\delta'}{2}=\frac{\sin \delta'}{1+\cos \delta'} =\frac{\frac{\sin \delta}{\gamma\cdot(1+\beta\cdot\cos\delta)}}{\frac{1 + \beta \cdot \cos \delta}{1 + \beta \cdot \cos \delta}+\frac{\cos \delta +\beta}{1 + \beta \cdot \cos \delta}} =\frac{\sin \delta}{\gamma\cdot(1 + \beta \cdot \cos \delta + \cos \delta +\beta)}=\frac{\sin \delta}{\gamma\cdot(1 + \beta)\cdot(1+ \cos \delta)} =\frac{\tan \frac{\delta}{2}}{\gamma\cdot(1+ \beta)}$

And as $\frac{1}{\gamma\cdot(1+\beta)} =\frac{\sqrt{1-\beta^2}}{1+\beta} =\frac{\sqrt{(1+\beta)(1-\beta)}}{1+\beta} =\sqrt{\frac{1-\beta}{1+\beta}}$ the symmetric form follows.

Formula for motion along the y-axis

Let $\theta$ be the angle between the light ray (=path of the light signal which is a straight line) and the y-axis whereby θ is positive if the y-axis would have to rotate counter-clockwise to coincide with the light ray. Then the derivation of the formula of angle θ' for the motion along the y-axis is the same as the derivation of the formula of angle δ' for the motion along the x-axis.

Hence: $\tan\theta\,' = \frac{\sin\theta}{\gamma\cdot(\cos\theta+\beta)}$

The symmetric form is: $\tan \frac{\theta\,'}{2}=\sqrt{\frac{1-\beta}{1+\beta}}\cdot\tan\frac{\theta}{2}$ and the approximate formula is: $\triangle\theta = -\beta\cdot\sin\theta$

As $\theta=\delta-90^\circ$ and as $\scriptstyle\sin(\delta-90^{\circ})=-\sin(90^\circ-\delta)=-\cos\delta\;,\;\cos(\delta-90^\circ)=\cos(90^\circ-\delta)=\sin\delta$ and $\scriptstyle\tan(\delta-90^\circ)=-\tan(90^\circ-\delta)=-\cot\delta$  one gets:

$\tan(\delta'-90^\circ) = \frac{\sin(\delta-90^\circ)}{\gamma\cdot\left(\cos(\delta - 90 ^\circ)+\beta\right)}$ and therefore:  $-\cot\delta' = \frac{-\cos\delta}{\gamma\cdot(\sin\delta+\beta)}$ and hence $\cot\delta' = \frac{\cos\delta}{\gamma\cdot(\sin\delta+\beta)}$

Since $\cot\delta'=\frac{1}{\tan\delta'}$  one also gets: $\tan\delta' = \frac{\gamma\cdot(\sin\delta+\beta)}{\cos\delta}$

And the symmetric form is: $\tan \frac{\delta'-90^\circ}{2}=\sqrt{\frac{1-\beta}{1+\beta}}\cdot\tan\frac{\delta-90^\circ}{2}$

And as $\scriptstyle\triangle\theta=(\delta'-90^\circ)-(\delta-90^\circ)=\triangle\delta$  the approximate formula is: $\triangle\delta = \triangle\theta =-\beta\cdot\sin(\delta-90^\circ)=+\beta\cdot\cos\delta$   

Formula for motion along a ray lying in the x-y-plane with direction vector (cos α | sin α)

With the same reasoning as above one gets the formula: $\tan(\delta'-\alpha) = \frac{\sin(\delta-\alpha)}{\gamma\cdot\left(\cos(\delta-\alpha)+\beta\right)}$

The symmetric form is: $\tan \frac{\delta'-\alpha}{2}=\sqrt{\frac{1-\beta}{1+\beta}}\cdot\tan\frac{\delta-\alpha}{2}$

The approximate formula is: $\triangle\delta = - \beta\cdot\sin(\delta-\alpha)$  and in the degree measure: $\triangle\delta = -\frac{v}{c}\cdot\sin(\delta-\alpha)\cdot\frac{180^\circ}{\pi}$

Three-dimensional problem

Heliocentric ecliptic coordinates

Earth-centered Ecliptic coordinates

Application: Aberration in astronomy

The stellar aberration is only an effect of the change of the reference frame. The astronomer orbits (with Earth) around the Sun in one year and furthermore rotates around the axis of Earth in one day. His current rest frame S' therefore has different velocities relative to the rest frame S of the barycenter of the Solar System at different times. Hence the astronomer observes, that the position of the star changes. The formula is derived under the condition that the change of the position of the star and of Earth is negligible in the period of observation ^{[Notes 8]} That is correct for almost all stars: The amplitude of the parallax of a star in a distance of ≥ n parsec is ≤ 1/n ".

Stellar aberration due to the orbit of Earth (around the Sun)

The mean orbital speed of Earth is $v_e = \frac{2 \cdot \pi \cdot 1\,AE}{365,25\,d} = 29,78\,km/s$, and therefore $\frac{v_e}{c} = 0,00009935$.

-> $\frac{v}{c} \cdot \frac{180^\circ}{\pi} = 20,5^{\prime\prime}$.

$k_A= 20,5^{\prime\prime}$ is dubbed constant of aberration for the annual aberration.^[5]

Stellar aberration due to Earth's rotation

An astronomer at the latitude $\varphi$ rotates in 24 hours around the axis of Earth. His speed of rotation is therefore $v_r= \frac{\cos \varphi \cdot 40000\,km}{1\,d} = \cos \varphi \cdot 463\,m/s$. Hence $\frac{v_r}{c} = \cos \varphi \cdot 1,54\cdot10^{-6}$. Form this so called diurnal aberration one gets an additional contribution of (at max.) $\cos \varphi \cdot 0,32^{\prime\prime}$.

Stellar aberration due the orbit of our solar system around the center of the Milky Way

The rest frame of the center of mass of our solar system isn't a perfect inertial frame of reference since our solar system orbits around the center of the Milky Way. An estimation for the time of period of circulation is 230 million years (estimations vary between 225 and 250 million years).^[6][7] As the estimation for the distance between our solar system and the center of the Milky Way is about 28000 Ly, the assumed orbital speed of our solar system is $\scriptstyle v = \frac{2\pi \cdot 280000 \cdot 9,461\cdot10^{15}\,m}{230\cdot 10^{6} \cdot 365,25 \cdot 24 \cdot 3600} \approx\, 230\,km/s$. This would cause an aberration ellipse with a major semiaxis of 2,6' (arcminutes). Therefore, in one year the aberration angle could change (at max.) $\scriptstyle 2,6^\prime \cdot \frac{2 \pi \cdot 1\,a}{230000000\,a}$ = 4,3 µas (microarcseconds). This very small value isn't detectable now, perhaps it's possible with the planned mission of the Gaia spacecraft.^[8]

See also

• Aberration of light
• Parallax in astronomy

Notes

1. ↑ R. Emden's (rhetorical) question "How will the direction of the light rays of stars incident upon Earth be influenced by the motion of Earth and the motion (not the location) of the star at the time of emission?" is answered by him with "Not at all."
2. ↑ I.e. the stars can rotate very fast like eclipsing binary stars as long as the change in position is much smaller than the star-earth distance.
3. ↑ It's easy to (technically) compensate the deviation of the direction of the axes due to Earth's rotation and earth's orbit but it almost impossible to (technically) compensate a velocity in the order of several 1000 m/s. This can only be achieved with help of mathematics.
4. ↑ Of course every inertial frame of reference is physically equal, but the rest system of the barycenter has the mathematical advantage, that it's the mean frame of reference and since the time of Galileo Galilei it's the preferred frame of reference
5. ↑ The lines of x-axis and x'-axis lie on top of each other
6. ↑ An even better quasi-inertial frame of reference would be a frame of reference in which the center of mass of the Milky Way would rest or be in uniform motion.
7. ↑ Here the solution with cot 0 = -π/2 has to be used, see Plot of arccot (x)
8. ↑ For the star this period lies in the past.

References

1. ↑ Herschel, John Frederick William (1844). "Schreiben des Herrn Baronets Herschel an den Herausgeber". Astronomische Nachrichten 22 (520): 249–254. Bibcode:1844AN.....22..249H. doi:10.1002/asna.18450221702.
2. ↑ Emden, R. (1926). "Aberration und Relativitätstheorie". Die Naturwissenschaften 14 (16): 329–335. Bibcode:1926NW.....14..329E. doi:10.1007/BF01506966. (14. Jahrgang, Heft 16)
3. ↑ See for example equation (4) in Klaus Kassner (2002). "Why the Bradley aberration cannot be used to measure absolute speeds. A comment.". arXiv:astro-ph/0203056 [astro-ph].
4. ↑ Spezielle Relativitätstheorie für Studienanfänger, Jürgen Freund, vdf-Verlag (Zürich), 2nd. Ed. (2005)
5. ↑ see e.g. The Constant of Aberration and the Solar Parallax
6. ↑ see hypertextbook.com
7. ↑ NASA - StarChild Question of the Month for February 2000
8. ↑ This is a hint from H.-H. Voigt: Abriss der Astronomie, 6. Auflage. Herausgegeben von H.-J. Röser und W. Tscharnuter. © 2012 Wiley-VCH Verlag GmbH & Co. KGaA., page 35 of extract