On the dynamics of the electron: Introduction, §§1, 9

Henri Poincaré Translated by Scott A. Walter from Rendiconti del Circolo Matematico di Palermo 21, 1906, 129–176.
Published in J. Renn and M. Schemmel (eds.), The Genesis of General Relativity Vol. 3: Theories of Gravitation in the Twilight of Classical Physics; Part I (Boston Studies in the Philosophy of Science 250), Springer, 2007, 253–271

Translator’s note: The original notation is faithfully reproduced, including the use of ‘‘$d$’’ for both ordinary and partial differentiation. The translator’s footnote calls use arabic numbers; Poincaré’s single footnote call is marked by an asterisk. For alternative translations of Poincaré’s memoir see C. W. Kilmister (Special Theory of Relativity, Oxford: Pergamon, 1970, 145–185), and H. M. Schwartz (American Journal of Physics 39:1287–1294; 40:862–872, 1282–1287).

Introduction.

It seems at first that the aberration of light and related optical and electrical phenomena will provide us with a means of determining the absolute motion of the Earth, or rather its motion with respect to the ether, as opposed to its motion with respect to other celestial bodies. Fresnel pursued this idea, but soon recognized that the Earth’s motion does not alter the laws of refraction and reflection. Analogous experiments, like that of the water-filled telescope, and all those considering terms no higher than first order relative to the aberration, yielded only negative results; the explanation was soon discovered. But Michelson, who conceived an experiment sensitive to terms depending on the square of the aberration, failed in turn.

It appears that this impossibility to detect the absolute motion of the Earth by experiment may be a general law of nature; we are naturally inclined to admit this law, which we will call the Postulate of Relativity and admit without restriction. Whether or not this postulate, which up to now agrees with experiment, may later be corroborated or disproved by experiments of greater precision, it is interesting in any case to ascertain its consequences.

An explanation was proposed by Lorentz and FitzGerald, who introduced the hypothesis of a contraction of all bodies in the direction of the Earth’s motion and proportional to the square of the aberration. This contraction, which we will call the Lorentzian contraction, would explain Michelson’s experiment and all others performed up to now. The hypothesis would become insufficient, however, if we were to admit the postulate of relativity in full generality.

Lorentz then sought to extend his hypothesis and to modify it in order to obtain perfect agreement with this postulate. This is what he succeeded in doing in his article entitled Electromagnetic phenomena in a system moving with any velocity smaller that that of light (Proceedings of the Amsterdam Academy, May 27, 1904).

The importance of the question persuaded me to take it up in turn; the results I obtained agree with those of Mr. Lorentz on all the significant points. I was led merely to modify and extend them only in a few details; further on we will see the points of divergence, which are of secondary importance.

Lorentz’s idea may be summed up like this: if we are able to impress a translation upon an entire system without modifying any observable phenomena, it is because the equations of an electromagnetic medium are unaltered by certain transformations, which we will call Lorentz transformations. Two systems, one of which is at rest, the other in translation, become thereby exact images of each other.

Langevin**Langevin was anticipated by Mr. Bucherer of Bonn, who earlier advanced the same idea. (See: Bucherer, Mathematische Einführung in die Elektronentheorie, August, 1904. Teubner, Leipzig).) sought to modify Lorentz’s idea; for both authors, the moving electron takes the form of a flattened ellipsoid. For Lorentz, two axes of the ellipsoid remain constant, while for Langevin, ellipsoid volume remains constant. The two scientists also showed that these two hypotheses are corroborated by Kaufmann’s experiments to the same extent as the original hypothesis of Abraham (rigid-sphere electron).

The advantage of Langevin’s theory is that it requires only electromagnetic forces, and bonds; it is, however, incompatible with the postulate of relativity. This is what Lorentz showed, and this is what I found in turn using a different method, which calls on principles of group theory.

We must return therefore to Lorentz’s theory, but if we want to do this and avoid intolerable contradictions, we must posit the existence of a special force that explains both the contraction, and the constancy of two of the axes. I sought to determine this force, and found that it may be assimilated to a constant external pressure on the deformable and compressible electron, whose work is proportional to the electron’s change in volume.

If the inertia of matter is exclusively of electromagnetic origin, as generally admitted in the wake of Kaufmann’s experiment, and all forces are of electromagnetic origin (apart from this constant pressure that I just mentioned), the postulate of relativity may be established with perfect rigor. This is what I show by a very simple calculation based on the principle of least action.

But that is not all. In the article cited above, Lorentz judged it necessary to extend his hypothesis in such a way that the postulate remains valid in case there are forces of non-electromagnetic origin. According to Lorentz, all forces are affected by the Lorentz transformation (and consequently by a translation) in the same way as electromagnetic forces.

It was important to examine this hypothesis closely, and in particular to ascertain the modifications we would have to apply to the laws of gravitation.

We find first of all that it requires us to assume that gravitational propagation is not instantaneous, but occurs with the speed of light. One might think that this is reason enough to reject the hypothesis, since Laplace demonstrated that this cannot be the case. In reality, however, the effect of this propagation is compensated in large part by a different cause, in such a way that no contradiction arises between the proposed law and astronomical observations.

Is it possible to find a law satisfying Lorentz’s condition, and reducing to Newton’s law whenever the speeds of celestial bodies are small enough to allow us to neglect their squares (as well as the product of acceleration and distance) with respect to the square of the speed of light?

To this question we must respond in the affirmative, as we will see later.

Modified in this way, is the law compatible with astronomical observations?

It seems so on first sight, but the question will be settled only after an extended discussion.

Suppose, then, that this discussion is settled in favor of the new hypothesis, what should we conclude? If propagation of attraction occurs with the speed of light, it could not be a fortuitous accident. Rather, it must be because it is a function of the ether, and then we would have to try to penetrate the nature of this function, and to relate it to other fluid functions.

We cannot be content with a simple juxtaposition of formulas that agree with each other by good fortune alone; these formulas must, in a manner of speaking, interpenetrate. The mind will be satisfied only when it believes it has perceived the reason for this agreement, and the belief is strong enough to entertain the illusion that it could have been predicted.

But the question may be viewed from a different perspective, better shown via an analogy. Let us imagine a pre-Copernican astronomer who reflects on Ptolemy’s system; he will notice that for all the planets, one of two circles – epicycle or deferent – is traversed in the same time. This fact cannot be due to chance, and consequently between all the planets there is a mysterious link we can only guess at.

Copernicus, however, destroys this apparent link by a simple change in the coordinate axes that were considered fixed. Each planet now describes a single circle, and orbital periods become independent (until Kepler reestablishes the link that was believed to have been destroyed).

It is possible that something analogous is taking place here. If we were to admit the postulate of relativity, we would find the same number in the law of gravitation and the laws of electromagnetism—the speed of light—and we would find it again in all other forces of any origin whatsoever. This state of affairs may be explained in one of two ways: either everything in the universe would be of electromagnetic origin, or this aspect—shared, as it were, by all physical phenomena—would be a mere epiphenomenon, something due to our methods of measurement. How do we go about measuring? The first response will be: we transport objects considered to be invariable solids, one on top of the other. But that is no longer true in the current theory if we admit the Lorentzian contraction. In this theory, two lengths are equal, by definition, if they are traversed by light in equal times.

Perhaps if we were to abandon this definition Lorentz’s theory would be as fully overthrown as was Ptolemy’s system by Copernicus’s intervention. Should that happen some day, it would not prove that Lorentz’s efforts were in vain, because regardless of what one may think, Ptolemy was useful to Copernicus.

I, too, have not hesitated to publish these few partial results, even if at this very moment the discovery of magneto-cathode rays seems to threaten the entire theory.

§ 1. — Lorentz Transformation.

Lorentz adopted a certain system of units in order to do away with $4\pi$ factors in formulas. I will do the same, and in addition, select units of length and time in such a way that the speed of light equals 1. Under these conditions, and denoting electric displacement $f$, $g$, $h$, magnetic intensity $\alpha$, $\beta$, $\gamma$, vector potential $F$, $G$, $H$, scalar potential $\psi$, charge density $\rho$, electron velocity $\xi$, $\eta$, $\zeta$, and current $u$, $v$, $w$, the fundamental formulas become:

 \left.\begin{aligned} \displaystyle u&\displaystyle=\frac{df}{dt}+\rho\xi=% \frac{d\gamma}{dy}-\frac{d\beta}{dz},\quad\alpha=\frac{dH}{dy}-\frac{dG}{dz},% \quad f=-\frac{dF}{dt}-\frac{d\psi}{dx},\\ \displaystyle\frac{d\alpha}{dt}&\displaystyle=\frac{dg}{dz}-\frac{dh}{dy},% \quad\frac{d\rho}{dt}+\sum\frac{d\rho\xi}{dx}=0,\quad\sum\frac{df}{dx}=\rho,% \quad\frac{d\psi}{dt}+\sum\frac{dF}{dx}=0,\\ \displaystyle\Box&\displaystyle=\Delta-\frac{d^{2}}{dt^{2}}=\sum\frac{d^{2}}{% dx^{2}}-\frac{d^{2}}{dt^{2}},\qquad\Box\psi=-\rho,\qquad\Box F=-\rho\xi.\end{% aligned}\right\} (1)

An elementary particle of matter of volume $dxdydz$ is acted upon by a mechanical force, the components of which are derived from the formula:

 $X=\rho f+\rho(\eta\gamma-\zeta\beta).$ (2)

These equations admit a remarkable transformation discovered by Lorentz, which owes its interest to the fact that it explains why no experiment can inform us of the absolute motion of the universe. Let us put:

 $x^{\prime}=kl(x+\varepsilon t),\qquad t^{\prime}=kl(t+\varepsilon x),\qquad y^% {\prime}=\ell y,\qquad z^{\prime}=\ell z,$ (3)

where $\ell$ and $\varepsilon$ are two arbitrary constants, such that

 $k=\frac{1}{\sqrt{1-\varepsilon^{2}}}.$

Now if we put:

 $\Box^{\prime}=\sum\frac{d^{2}}{dx^{\prime}{}^{2}}-\frac{d^{2}}{dt^{\prime}{}^{% 2}},$

we will have:

 $\Box^{\prime}=\Box\ell^{-2}.$

Let a sphere be carried along with the electron in uniform translation, and let the equation of this mobile sphere be:

 $(x-\xi t)^{2}+(y-\eta t)^{2}+(z-\zeta t)^{2}=r^{2},$

and the volume of the sphere be $\frac{4}{3}\pi r^{3}$.11The original reads: ‘‘$\frac{4}{3}\pi r^{2}$’’.

The transformation will change the sphere into an ellipsoid, the equation of which is easy to find. We thus deduce easily from (3):

 $x=\frac{k}{\ell}(x^{\prime}-\varepsilon t^{\prime}),\quad t=\frac{k}{\ell}(t^{% \prime}-\varepsilon x^{\prime}),\quad y=\frac{y^{\prime}}{\ell},\quad z=\frac{% z^{\prime}}{\ell}.{}$ (3’)

The equation of the ellipsoid then becomes:

 $k^{2}(x^{\prime}-\varepsilon t^{\prime}-\xi t^{\prime}+\varepsilon\xi x^{% \prime})^{2}+(y^{\prime}-\eta kt^{\prime}+\eta k\varepsilon x^{\prime})^{2}+(z% ^{\prime}-\zeta kt^{\prime}+\zeta k\varepsilon x^{\prime})^{2}=\ell^{2}r^{2}.$

This ellipsoid is in uniform motion; for $t^{\prime}=0$, it reduces to

 $k^{2}x^{\prime}{}^{2}(1+\xi\varepsilon)^{2}+(y^{\prime}+\eta k\varepsilon x^{% \prime})^{2}+(z^{\prime}+\zeta k\varepsilon x^{\prime})^{2}=\ell^{2}r^{2},$

and has a volume:

 $\frac{4}{3}\pi r^{3}\frac{\ell^{3}}{k(1+\xi\varepsilon)}.$

If we want electron charge to be unaltered by the transformation, and if we designate the new charge density $\rho^{\prime}$, we will find:

 $\rho^{\prime}=\frac{k}{\ell^{3}}(\rho+\varepsilon\rho\xi).$ (4)

What will be the new velocity components $\xi^{\prime}$, $\eta^{\prime}$ and $\zeta^{\prime}$? We should have:

 $\displaystyle\xi^{\prime}=\frac{dx^{\prime}}{dt^{\prime}}=\frac{d(x+% \varepsilon t)}{d(t+\varepsilon x)}=\frac{\xi+\varepsilon}{1+\varepsilon\xi},$ $\displaystyle\eta^{\prime}=\frac{dy^{\prime}}{dt^{\prime}}=\frac{dy}{kd(t+% \varepsilon x)}=\frac{\eta}{k(1+\varepsilon\xi)},\qquad\zeta^{\prime}=\frac{% \zeta}{k(1+\varepsilon\xi)},$

whence:

 $\rho^{\prime}\xi^{\prime}=\frac{k}{\ell^{3}}(\rho\xi+\varepsilon\rho),\qquad% \rho^{\prime}\eta^{\prime}=\frac{1}{\ell^{3}}\rho\eta,\qquad\rho^{\prime}\zeta% ^{\prime}=\frac{1}{\ell^{3}}\rho\zeta.{}$ (4’)

Here is where I must point out for the first time a difference with Lorentz. In my notation, Lorentz put (l.c., page 813, formulas 7 and 8):

 $\rho^{\prime}=\frac{1}{k\ell^{3}}\rho,\qquad\xi^{\prime}=k^{2}(\xi+\varepsilon% ),\qquad\eta^{\prime}=k\eta,\qquad\zeta^{\prime}=k\zeta.$

In this way we recover the formulas:

 $\rho^{\prime}\xi^{\prime}=\frac{k}{\ell^{3}}(\rho\xi+\varepsilon\rho),\qquad% \rho^{\prime}\eta^{\prime}=\frac{1}{\ell^{3}}\rho\eta,\qquad\rho^{\prime}\zeta% ^{\prime}=\frac{1}{\ell^{3}}\rho\zeta;$

although the value of $\rho^{\prime}$ differs.

It is important to notice that the formulas (4) and ($4^{\prime}$) satisfy the condition of continuity

 $\frac{d\rho^{\prime}}{dt^{\prime}}+\sum\frac{d\rho^{\prime}\xi^{\prime}}{dx^{% \prime}}=0.$

To see this, let $\lambda$ be an undetermined coefficient and $D$ the Jacobian of

 $t+\lambda\rho,\qquad x+\lambda\rho\xi,\qquad y+\lambda\rho\eta,\qquad z+% \lambda\rho\zeta$ (5)

with respect to $t$, $x$, $y$ and $z$. It follows that:

 $D=D_{0}+D_{1}\lambda+D_{2}\lambda^{2}+D_{3}\lambda^{3}+D_{4}\lambda^{4},$

with $D_{0}=1$, $D_{1}=\displaystyle\frac{d\rho}{dt}+\sum\frac{d\rho\xi}{dx}=0$.

Let $\lambda^{\prime}=\ell^{4}\rho^{\prime}$;22The original reads: ‘‘$\lambda^{\prime}=\ell^{2}\rho^{\prime}$’’. then the 4 functions

 $t^{\prime}+\lambda^{\prime}\rho^{\prime},\qquad x^{\prime}+\lambda^{\prime}% \rho^{\prime}\xi^{\prime},\qquad y^{\prime}+\lambda^{\prime}\rho^{\prime}\eta^% {\prime},\qquad z^{\prime}+\lambda^{\prime}\rho^{\prime}\zeta^{\prime}{}$ (5’)

are related to the functions (5) by the same linear relationships as the old variables to the new ones. Therefore, if we denote $D^{\prime}$ the Jacobian of the functions ($5^{\prime}$) with respect to the new variables, it follows that:

 $D^{\prime}=D,\qquad D^{\prime}=D_{0}^{\prime}+D_{1}^{\prime}\lambda^{\prime}+% \cdots+D_{4}^{\prime}\lambda^{\prime}{}^{4},$

and thereby:33The original reads: ‘‘$D_{1}^{\prime}=\ell^{-2}D_{1}$’’.

 $D_{0}^{\prime}=D_{0}=1,\qquad D_{1}^{\prime}=\ell^{-4}D_{1}=0=\frac{d\rho^{% \prime}}{dt^{\prime}}+\sum\frac{d\rho^{\prime}\xi^{\prime}}{dx^{\prime}}.% \qquad\text{Q.E.D.}$

Under Lorentz’s hypothesis, this condition would not be met since $\rho^{\prime}$ has a different value.

We will define the new vector and scalar potentials in such a way as to satisfy the conditions

 $\Box^{\prime}\psi^{\prime}=-\rho^{\prime},\qquad\Box^{\prime}F^{\prime}=-\rho^% {\prime}\xi^{\prime}.$ (6)

From this we deduce:

 $\psi^{\prime}=\frac{k}{\ell}(\psi+\varepsilon F),\quad F^{\prime}=\frac{k}{% \ell}(F+\varepsilon\psi),\quad G^{\prime}=\frac{1}{\ell}G,\quad H^{\prime}=% \frac{1}{\ell}H.$ (7)

These formulas differ noticeably from those of Lorentz, although the divergence stems ultimately from the definitions employed.

New electric and magnetic fields are now chosen in order to satisfy the equations

 $f^{\prime}=-\frac{dF^{\prime}}{dt^{\prime}}-\frac{d\psi^{\prime}}{dx^{\prime}}% ,\qquad\alpha^{\prime}=\frac{dH^{\prime}}{dy^{\prime}}-\frac{dG^{\prime}}{dz^{% \prime}}.$ (8)

It is easy to see that:

 $\frac{d}{dt^{\prime}}=\frac{k}{\ell}\left(\frac{d}{dt}-\varepsilon\frac{d}{dx}% \right),\quad\frac{d}{dx^{\prime}}=\frac{k}{\ell}\left(\frac{d}{dx}-% \varepsilon\frac{d}{dt}\right),\quad\frac{d}{dy^{\prime}}=\frac{1}{\ell}\frac{% d}{dy},\quad\frac{d}{dz^{\prime}}=\frac{1}{\ell}\frac{d}{dz}$

and we deduce thereby:

 \left.\begin{aligned} \displaystyle f^{\prime}&\displaystyle=\frac{1}{\ell^{2}% }f,&\hskip 10.0pt\qquad g^{\prime}&\displaystyle=\frac{k}{\ell^{2}}(g+% \varepsilon\gamma),&\displaystyle\qquad h^{\prime}&\hskip 10.0pt=\frac{k}{\ell% ^{2}}(h-\varepsilon\beta),\\ \displaystyle\alpha^{\prime}&\displaystyle=\frac{1}{\ell^{2}}\alpha,&\hskip 10% .0pt\qquad\beta^{\prime}&\displaystyle=\frac{k}{\ell^{2}}(\beta-\varepsilon h)% ,&\displaystyle\qquad\gamma^{\prime}&\hskip 10.0pt=\frac{k}{\ell^{2}}(\gamma+% \varepsilon g).\end{aligned}\quad\right\} (9)

These formulas are identical to those of Lorentz.

Our transformation does not alter (1). In fact, the condition of continuity, as well as (6) and (8) were already featured in (1) (neglecting the primes).

Combining (6) with the condition of continuity, we obtain:

 $\frac{d\psi^{\prime}}{dt^{\prime}}+\sum\frac{dF^{\prime}}{dx^{\prime}}=0.$ (10)

It remains for us to establish:

 $\frac{df^{\prime}}{dt^{\prime}}+\rho^{\prime}\xi^{\prime}=\frac{d\gamma^{% \prime}}{dy^{\prime}}-\frac{d\beta^{\prime}}{dz^{\prime}},\qquad\frac{d\alpha^% {\prime}}{dt^{\prime}}=\frac{dg^{\prime}}{dz^{\prime}}-\frac{dh^{\prime}}{dy^{% \prime}},\qquad\sum\frac{df^{\prime}}{dx^{\prime}}=\rho^{\prime}$

and it is easy to see that these are necessary consequences of (6), (8) and (10).

We must now compare forces before and after the transformation.

Let $X$, $Y$, $Z$ be the force prior to the transformation, and $X^{\prime}$, $Y^{\prime}$, $Z^{\prime}$ the force after the transformation, both forces being per unit volume. In order for $X^{\prime}$ to satisfy the same equations as before the transformation, we must have:

 $\displaystyle X^{\prime}$ $\displaystyle=\rho^{\prime}f^{\prime}+\rho^{\prime}(\eta^{\prime}\gamma^{% \prime}-\zeta^{\prime}\beta^{\prime}),$ $\displaystyle Y^{\prime}$ $\displaystyle=\rho^{\prime}g^{\prime}+\rho^{\prime}(\zeta^{\prime}\alpha^{% \prime}-\xi^{\prime}\gamma^{\prime}),$ $\displaystyle Z^{\prime}$ $\displaystyle=\rho^{\prime}h^{\prime}+\rho^{\prime}(\xi^{\prime}\beta^{\prime}% -\eta^{\prime}\alpha^{\prime}),$

or, replacing all quantities by their values (4), ($4^{\prime}$) and (9), and in light of (2):

 \left.\begin{aligned} \displaystyle X^{\prime}&\displaystyle=\frac{k}{\ell^{5}% }(X+\varepsilon\sum X\xi),\\ \displaystyle Y^{\prime}&\displaystyle=\frac{1}{\ell^{5}}Y,\\ \displaystyle Z^{\prime}&\displaystyle=\frac{1}{\ell^{5}}Z.\end{aligned}\quad\right\} (11)

Instead of representing the components of force per unit volume by $X_{1}$, $Y_{1}$, $Z_{1}$, we now let these terms represent the force per unit electron charge, and we let $X_{1}^{\prime}$, $Y_{1}^{\prime}$, $Z_{1}^{\prime}$ represent the latter force after transformation. It follows that:

 $X_{1}=f+\eta\gamma-\zeta\beta,\quad X_{1}^{\prime}=f^{\prime}+\eta^{\prime}% \gamma^{\prime}-\zeta^{\prime}\beta^{\prime},\quad X=\rho X_{1},\quad X^{% \prime}=\rho X_{1}^{\prime},$

and we obtain the equations:

 \left.\begin{aligned} \displaystyle X_{1}^{\prime}&\displaystyle=\frac{k}{\ell% ^{5}}\frac{\rho}{\rho^{\prime}}(X_{1}+\varepsilon\sum X_{1}\xi),\\ \displaystyle Y_{1}^{\prime}&\displaystyle=\frac{1}{\ell^{5}}\frac{\rho}{\rho^% {\prime}}Y_{1},\\ \displaystyle Z_{1}^{\prime}&\displaystyle=\frac{1}{\ell^{5}}\frac{\rho}{\rho^% {\prime}}Z_{1}.\end{aligned}\quad\right\}{} (11’)

Lorentz found (page 813, equation (10) with different notation):

 \left.\begin{aligned} \displaystyle X_{1}&\displaystyle=\ell^{2}X_{1}^{\prime}% -\ell^{2}\varepsilon(\eta^{\prime}g^{\prime}+\zeta^{\prime}h^{\prime}),\\ \displaystyle Y_{1}&\displaystyle=\frac{\ell^{2}}{k}Y_{1}^{\prime}+\frac{\ell^% {2}\varepsilon}{k}\xi^{\prime}g^{\prime},\\ \displaystyle Z_{1}&\displaystyle=\frac{\ell^{2}}{k}Z_{1}^{\prime}+\frac{\ell^% {2}\varepsilon}{k}\xi^{\prime}h^{\prime}.\end{aligned}\quad\right\}{} (11’’)

Before going any further, it is important to locate the source of this significant divergence. It obviously springs from the fact that the formulas for $\xi^{\prime}$, $\eta^{\prime}$ and $\zeta^{\prime}$ are not the same, while the formulas for the electric and magnetic fields are the same.

If electron inertia is exclusively of electromagnetic origin, and if electrons are subject only to forces of electromagnetic origin, then the conditions of equilibrium require that:

 $X=Y=Z=0$

inside the electrons.

According to (11), these relationships are equivalent to

 $X^{\prime}=Y^{\prime}=Z^{\prime}=0.$

The electron’s equilibrium conditions are therefore unaltered by the transformation.

Unfortunately, such a simple hypothesis is inadmissible. In fact, if we assume $\xi=\eta=\zeta=0$, the condition $X=Y=Z=0$ leads necessarily to $f=g=h=0$, and consequently, to $\displaystyle\sum\frac{df}{dx}=0$, i.e., $\rho=0$. Similar results obtain for the most general case. We must then admit that in addition to electromagnetic forces there are either non-electromagnetic forces or bonds. Therefore, we need to identify the conditions that these forces or these bonds must satisfy for electron equilibrium to be undisturbed by the transformation. This will be the object of an upcoming section.

§ 9. — Hypotheses Concerning Gravitation.

In this way Lorentz’s theory would fully explain the impossibility of detecting absolute motion, if all forces were of electromagnetic origin.

But there exist other forces to which an electromagnetic origin cannot be attributed, such as gravitation, for example. It may in fact happen, that two systems of bodies produce equivalent electromagnetic fields, i.e., exert the same action on electrified bodies and on currents, and at the same time, these two systems do not exert the same gravitational action on Newtonian masses. The gravitational field is therefore distinct from the electromagnetic field. Lorentz was obliged thereby to extend his hypothesis with the assumption that forces of any origin whatsoever, and gravitation in particular, are affected by a translation (or, if one prefers, by the Lorentz transformation) in the same manner as electromagnetic forces.

It is now appropriate to enter into the details of this hypothesis, and to examine it more closely. If we want the Newtonian force to be affected by the Lorentz transformation in this fashion, we can no longer suppose that it depends only on the relative position of the attracting and attracted bodies at the instant considered. The force should also depend on the velocities of the two bodies. And that is not all: it will be natural to suppose that the force acting on the attracted body at the instant $t$ depends on the position and velocity of this body at this same instant $t$, but it will also depend on the position and velocity of the attracting body, not at the instant $t$, but at an earlier instant, as if gravitation had taken a certain time to propagate.

Let us now consider the position of the attracted body at the instant $t_{0}$, and let $x_{0}$, $y_{0}$, $z_{0}$ be its coordinates, and $\xi$, $\eta$, $\zeta$ its velocity components at this instant; let us consider also the attracting body at the corresponding instant $t_{0}+t$, and let its coordinates be $x_{0}+x$, $y_{0}+y$, $z_{0}+z$, and its velocity components be $\xi_{1}$, $\eta_{1}$, $\zeta_{1}$ at this instant.

First we should have a relationship

 $\varphi\,(t,\ x,\ y,\ z,\ \xi,\ \eta,\ \zeta,\ \xi_{1},\ \eta_{1},\ \zeta_{1})=0$ (1)

in order to define the time $t$. This relationship will define the law of propagation of gravitational action (I do not constrain myself by any means to a propagation velocity equal in all directions).

Now let $X_{1}$, $Y_{1}$, $Z_{1}$ be the three components of the action exerted on the attracted body at the instant $t_{0}$;44The original reads: ‘‘à l’instant $t$’’. we want to express $X_{1}$, $Y_{1}$, $Z_{1}$ as functions of

 $t,\ x,\ y,\ z,\ \xi,\ \eta,\ \zeta,\ \xi_{1},\ \eta_{1},\ \zeta_{1}.$ (2)

What conditions must be satisfied?

1° The condition (1) should not be altered by transformations of the Lorentz group.

2° The components $X_{1}$, $Y_{1}$, $Z_{1}$ should be affected by transformations of the Lorentz group in the same manner as the electromagnetic forces designated by the same letters, i.e., in accordance with ($11^{\prime}$) of section 1.

3° When the two bodies are at rest, the ordinary law of attraction will be recovered.

It is important to note that in the latter case, the relationship (1) vanishes, because if the two bodies are at rest the time $t$ plays no role.

Posed in this fashion the problem is obviously indeterminate. We will therefore seek to satisfy to the utmost other, complementary conditions.

4° Since astronomical observations do not seem to show a sensible deviation from Newton’s law, we will choose the solution that differs the least with this law for small velocities of the two bodies.

5° We will make an effort to arrange matters in such a way that $t$ is always negative. Although we can imagine that the effect of gravitation requires a certain time in order to propagate, it would be difficult to understand how this effect could depend on the position not yet attained by the attracting body.

There is one case where the indeterminacy of the problem vanishes; it is the one where the two bodies are in mutual relative rest, i.e., where

 $\xi=\xi_{1},\qquad\eta=\eta_{1},\qquad\zeta=\zeta_{1};$

this is then the case we will examine first, by supposing that these velocities are constant, such that the two bodies are engaged in a common uniform rectilinear translation.

We may suppose that the $x-$axis is parallel to this translation, such that $\eta=\zeta=0$, and we will let $\varepsilon=-\xi$.

If we apply the Lorentz transformation under these conditions, after the transformation the two bodies will be at rest, and it follows that:

 $\xi^{\prime}=\eta^{\prime}=\zeta^{\prime}=0.$

The components $X_{1}$, $Y_{1}$, $Z_{1}$ should then agree with Newton’s law and we will have, apart from a constant factor:

 $X_{1}^{\prime}=-\frac{x}{r^{\prime}{}^{3}},\quad Y_{1}^{\prime}=-\frac{y}{r^{% \prime}{}^{3}},\quad Z_{1}^{\prime}=-\frac{z}{r^{\prime}{}^{3}},\quad r^{% \prime}{}^{2}=x^{\prime}{}^{2}+y^{\prime}{}^{2}+z^{\prime}{}^{2}.$ (3)

But according to section 1 we have:

 \begin{gathered}\displaystyle\begin{aligned} \displaystyle x^{\prime}&% \displaystyle=k(x+\varepsilon t),\quad y^{\prime}=y,\quad z^{\prime}=z,\quad t% ^{\prime}=k(t+\varepsilon x),\\ \displaystyle\frac{\rho^{\prime}}{\rho}&\displaystyle=k(1+\xi\varepsilon)=k(1-% \varepsilon^{2})=\frac{1}{k},\quad\sum X_{1}\xi=-X_{1}\varepsilon,\end{aligned% }\\ \displaystyle\begin{aligned} \displaystyle X_{1}^{\prime}&\displaystyle=k\frac% {\rho}{\rho^{\prime}}(X_{1}+\varepsilon\sum X_{1}\xi)=k^{2}X_{1}(1-\varepsilon% ^{2})=X_{1},\\ \displaystyle Y_{1}^{\prime}&\displaystyle=k\frac{\rho}{\rho^{\prime}}Y_{1}=kY% _{1}\\ \displaystyle Z_{1}^{\prime}&\displaystyle=kZ_{1}.\end{aligned}\end{gathered}

 $x+\varepsilon t=x-\xi t,\qquad r^{\prime}{}^{2}=k^{2}(x-\xi t)^{2}+y^{2}+z^{2}$

and

 $X_{1}=\frac{-k(x-\xi t)}{r^{\prime}{}^{3}},\quad Y_{1}=\frac{-y}{kr^{\prime}{}% ^{3}},\quad Z_{1}=\frac{-z}{kr^{\prime}{}^{3}};$ (4)

which may be written:

 $X_{1}=\frac{dV}{dx},\quad Y_{1}=\frac{dV}{dy},\quad Z_{1}=\frac{dV}{dz};\quad V% =\frac{1}{kr^{\prime}}.{}$ (4’)

It seems at first that the indeterminacy remains, since we made no hypotheses concerning the value of $t$, i.e., the transmission speed; and that besides, $x$ is a function of $t$. It is easy to see, however, that the terms appearing in our formulas, $x-\xi t$, $y$, $z$, do not depend on $t$.

We see that if the two bodies translate together, the force acting on the attracted body is perpendicular to an ellipsoid, at the center of which lies the attracting body.

To advance further, we need to look for the invariants of the Lorentz group.

We know that the substitutions of this group (assuming $\ell=1$) are linear substitutions that leave unaltered the quadratic form

 $x^{2}+y^{2}+z^{2}-t^{2}.$

Let us also put:

 $\displaystyle\xi$ $\displaystyle=\frac{\delta x}{\delta t},$ $\displaystyle\qquad\eta$ $\displaystyle=\frac{\delta y}{\delta t},$ $\displaystyle\qquad\zeta$ $\displaystyle=\frac{\delta z}{\delta t},$ $\displaystyle\xi_{1}$ $\displaystyle=\frac{\delta_{1}x}{\delta_{1}t},$ $\displaystyle\quad\eta_{1}$ $\displaystyle=\frac{\delta_{1}y}{\delta_{1}t},$ $\displaystyle\quad\zeta_{1}$ $\displaystyle=\frac{\delta_{1}z}{\delta_{1}t};$

we see that the Lorentz transformation will make $\delta x$, $\delta y$, $\delta z$, $\delta t$, and $\delta_{1}x$, $\delta_{1}y$, $\delta_{1}z$, $\delta_{1}t$ undergo the same linear substitutions as $x$, $y$, $z$, $t$.

Let us regard

 $\begin{matrix}x,&y,&z,&t\sqrt{-1},\\ \delta x,&\delta y,&\delta z,&\delta t\sqrt{-1},\\ \delta_{1}x,&\delta_{1}y,&\delta_{1}z,&\delta_{1}t\sqrt{-1},\end{matrix}$

as the coordinates of 3 points $P$, $P^{\prime}$, $P^{\prime\prime}$ in space of 4 dimensions. We see that the Lorentz transformation is merely a rotation in this space about the origin, assumed fixed. Consequently, we will have no distinct invariants apart from the 6 distances between the 3 points $P$, $P^{\prime}$, $P^{\prime\prime}$, considered separately and with the origin, or, if one prefers, apart from the two expressions

 $x^{2}+y^{2}+z^{2}-t^{2},\qquad x\delta x+y\delta y+z\delta z-t\delta t,$

or the 4 expressions of like form deduced from an arbitrary permutation of the 3 points $P$, $P^{\prime}$, $P^{\prime\prime}$.

But what we seek are invariants that are functions of the 10 variables (2). Therefore, among the combinations of our 6 invariants we must find those depending only on these 10 variables, i.e., those that are 0th degree homogeneous with respect both to $\delta x$, $\delta y$, $\delta z$, $\delta t$, and to $\delta_{1}x$, $\delta_{1}y$, $\delta_{1}z$, $\delta_{1}t$. We will then be left with 4 distinct invariants:

 $\sum x^{2}-t^{2},\quad\frac{t-\sum x\xi}{\sqrt{1-\sum\xi^{2}}},\quad\frac{t-% \sum x\xi_{1}}{\sqrt{1-\sum\xi_{1}^{2}}},\quad\frac{1-\sum\xi\xi_{1}}{\sqrt{% \left(1-\sum\xi^{2}\right)\left(1-\sum\xi_{1}^{2}\right)}}.$ (5)

Next let us see how the force components are transformed; we recall the equations (11) of section 1, that refer not to the force $X_{1}$, $Y_{1}$, $Z_{1}$ considered at present, but to the force per unit volume: $X$, $Y$, $Z$.

We designate moreover

 $T=\sum X\xi;$

we will see that (11) can be written ($\ell=1$):

 \left.\begin{aligned} \displaystyle X^{\prime}&\displaystyle=k(X+\varepsilon T% ),&\hskip 10.0pt\qquad T^{\prime}&\displaystyle=k(T+\varepsilon X),\\ \displaystyle Y^{\prime}&\displaystyle=Y,&\hskip 10.0pt\qquad Z^{\prime}&% \displaystyle=Z;\end{aligned}\quad\right\} (6)

in such a way that $X$, $Y$, $Z$, $T$ undergo the same transformation as $x$, $y$, $z$, $t$. Consquently, the group invariants will be

 $\sum X^{2}-T^{2},\quad\sum Xx-Tt,\quad\sum X\delta x-T\delta t,\quad\sum X% \delta_{1}x-T\delta_{1}t.$

However, it is not $X$, $Y$, $Z$ that we need, but $X_{1}$, $Y_{1}$, $Z_{1}$, with

 $T_{1}=\sum X_{1}\xi.$

We see that

 $\frac{X_{1}}{X}=\frac{Y_{1}}{Y}=\frac{Z_{1}}{Z}=\frac{T_{1}}{T}=\frac{1}{\rho}.$

Therefore, the Lorentz transformation will act in the same manner on $X_{1}$, $Y_{1}$, $Z_{1}$, $T_{1}$, as on $X$, $Y$, $Z$, $T$, except that these expressions will be multiplied moreover by

 $\frac{\rho}{\rho^{\prime}}=\frac{1}{k(1+\xi\varepsilon)}=\frac{\delta t}{% \delta t^{\prime}}.$

Likewise, the Lorentz transformation will act in the same way on $\xi$, $\eta$, $\zeta$, $1$ as on $\delta x$, $\delta y$, $\delta z$, $\delta t$, except that these expressions will be multiplied moreover by the same factor:

 $\frac{\delta t}{\delta t^{\prime}}=\frac{1}{k(1+\xi\varepsilon)}.$

Next we consider $X$, $Y$, $Z$, $T\sqrt{-1}$ as the coordinates of a fourth point $Q$; the invariants will then be functions of the mutual distances of the five points

 $0,\quad P,\quad P^{\prime},\quad P^{\prime\prime},\quad Q$

and among these functions we must retain only those that are 0th degree homogeneous with respect, on one hand, to

 $X,\quad Y,\quad Z,\quad T,\quad\delta x,\quad\delta y,\quad\delta z,\quad\delta t$

(variables that can be replaced further by $X_{1}$, $Y_{1}$, $Z_{1}$, $T_{1}$, $\xi$, $\eta$, $\zeta$, 1), and on the other hand, with respect to55The original reads ‘‘$\delta_{1}x$, $\delta_{1}y$, $\delta_{1}z$, 1.’’

 $\delta_{1}x,\qquad\delta_{1}y,\qquad\delta_{1}z,\qquad\delta_{1}t$

(variables that can be replaced further by $\xi_{1}$, $\eta_{1}$, $\zeta_{1}$, 1).

In this way we find, beyond the four invariants (5), four distinct new invariants:

 $\frac{\sum X_{1}^{2}-T_{1}^{2}}{1-\sum\xi^{2}},\quad\frac{\sum X_{1}x-T_{1}t}{% \sqrt{1-\sum\xi^{2}}},\quad\frac{\sum X_{1}\xi_{1}-T_{1}}{\sqrt{1-\sum\xi^{2}}% \sqrt{1-\sum\xi_{1}^{2}}},\quad\frac{\sum X_{1}\xi-T_{1}}{1-\sum\xi^{2}}.$ (7)

The latter invariant is always null according to the definition of $T_{1}$.

These terms being settled, what conditions must be satisfied?

1° The first term of (1), defining the velocity of propagation, has to be a function of the 4 invariants (5).

A wealth of hypotheses can obviously be entertained, of which we will examine only two:

A) We can have

 $\sum x^{2}-t^{2}=r^{2}-t^{2}=0,$

from whence $t=\pm r$, and, since $t$ has to be negative, $t=-r$. This means that the velocity of propagation is equal to that of light. It seems at first that this hypothesis ought to be rejected outright. Laplace showed in effect that the propagation is either instantaneous or much faster than that of light. However, Laplace examined the hypothesis of finite propagation velocity ceteris non mutatis; here, on the contrary, this hypothesis is conjoined with many others, and it may be that between them a more or less perfect compensation takes place. The application of the Lorentz transformation has already provided us with numerous examples of this.

B) We can have

 $\frac{t-\sum x\xi_{1}}{\sqrt{1-\sum\xi_{1}^{2}}}=0,\qquad t=\sum x\xi_{1}.$

The propagation velocity is therefore much faster than that of light, but in certain cases $t$ could be positive, which, as we mentioned, seems hardly admissible.66The original reads ‘‘$t$ pourrait être négatif.’’ We will therefore stick with hypothesis (A).

2° The four invariants (7) ought to be functions of the invariants (5).

3° When the two bodies are at absolute rest, $X_{1}$, $Y_{1}$, $Z_{1}$ ought to have the values given by Newton’s law, and when they are at relative rest, the values given by (4).

For the case of absolute rest, the first two invariants (7) ought to reduce to

 $\sum X_{1}^{2},\qquad\sum X_{1}x,$

or, by Newton’s law, to

 $\frac{1}{r^{4}},\qquad-\frac{1}{r};$

in addition, according to hypothesis (A), the $2^{\text{d}}$ and $3^{\text{rd}}$ invariants in (5) become:

 $\frac{-r-\sum x\xi}{\sqrt{1-\sum\xi^{2}}},\quad\frac{-r-\sum x\xi_{1}}{\sqrt{1% -\sum\xi_{1}^{2}}},$

that is, for absolute rest,

 $-r,\qquad-r.$

We may therefore admit, for example, that the first two invariants in (7) reduce to77The original has (4) instead of (7).

 $\frac{\left(1-\sum\xi_{1}^{2}\right)^{2}}{\left(r+\sum x\xi_{1}\right)^{4}},% \quad-\frac{\sqrt{1-\sum\xi_{1}^{2}}}{r+\sum x\xi_{1}},$

although other combinations are possible.

A choice must be made among these combinations, and furthermore, we need a $3^{\text{rd}}$ equation in order to define $X_{1}$, $Y_{1}$, $Z_{1}$. In making such a choice, we should try to come as close as possible to Newton’s law. Let us see what happens when we neglect the squares of the velocities $\xi$, $\eta$, etc. (still letting $t=-r$). The 4 invariants (5) then become:

 $0,\quad-r-\sum x\xi,\quad-r-\sum x\xi_{1},\quad 1$

and the 4 invariants (7) become:

 $\sum X_{1}^{2},\quad\sum X_{1}(x+\xi r),\quad\sum X_{1}(\xi_{1}-\xi),\quad 0.$

Before we can make a comparison with Newton’s law, another transformation is required. In the case under consideration, $x_{0}+x$, $y_{0}+y$, $z_{0}+z$, represent the coordinates of the attracting body at the instant $t_{0}+t$, and $r=\sqrt{\sum x^{2}}$. With Newton’s law we have to consider the coordinates of the attracting body $x_{0}+x_{1}$, $y_{0}+y_{1}$, $z_{0}+z_{1}$ at the instant $t_{0}$, and the distance $r_{1}=\sqrt{\sum x^{2}}$.

We may neglect the square of the time $t$ required for propagation, and proceed, consequently, as if the motion were uniform; we then have:

 $x=x_{1}+\xi_{1}t,\quad y=y_{1}+\eta_{1}t,\quad z=z_{1}+\zeta_{1}t,\quad r(r-r_% {1})=\sum x\xi_{1}t;$

or, since $t=-r$,

 $x=x_{1}-\xi_{1}r,\quad y=y_{1}-\eta_{1}r,\quad z=z_{1}-\zeta_{1}r,\quad r=r_{1% }-\sum x\xi_{1};$

such that our 4 invariants (5) become:

 $0,\quad-r_{1}+\sum x(\xi_{1}-\xi),\quad-r_{1},\quad 1$

and our 4 invariants (7) become:

 $\sum X_{1}^{2},\quad\sum X_{1}[x_{1}+(\xi-\xi_{1})r_{1}],\quad\sum X_{1}(\xi_{% 1}-\xi),\quad 0.$

In the second of these expressions I wrote $r_{1}$ instead of $r$, because $r$ is multiplied by $\xi-\xi_{1}$, and because I neglect the square of $\xi$.

For these 4 invariants (7), Newton’s law would yield

 $\frac{1}{r_{1}^{4}},\qquad-\frac{1}{r_{1}}-\frac{\sum x_{1}(\xi-\xi_{1})}{r_{1% }^{2}},\qquad\frac{\sum x_{1}(\xi-\xi_{1})}{r_{1}^{3}},\qquad 0.$

Therefore, if we designate the $2^{\text{nd}}$ and $3^{\text{rd}}$ of the invariants (5) as $A$ and $B$, and the first 3 invariants of (7) as $M$, $N$, $P$, we will satisfy Newton’s law to first-order terms in the square of velocity by setting:

 $M=\frac{1}{B^{4}},\quad N=\frac{+A}{B^{2}},\quad P=\frac{A-B}{B^{3}}.$ (8)

This solution is not unique. Let $C$ be the $4^{\text{th}}$ invariant in (5); $C-1$ is of the order of the square of $\xi$, and it is the same with $(A-B)^{2}$.

The solution (8) appears at first to be the simplest, nevertheless, it may not be adopted. In fact, since $M$, $N$, $P$ are functions of $X_{1}$, $Y_{1}$, $Z_{1}$, and $T_{1}=\sum X_{1}\xi$, the values of $X_{1}$, $Y_{1}$, $Z_{1}$ can be drawn from these three equations (8), but in certain cases these values would become imaginary.

To avoid this difficulty we will proceed in a different manner. Let us put:

 $k_{0}=\frac{1}{\sqrt{1-\sum\xi^{2}}},\quad k_{1}=\frac{1}{\sqrt{1-\sum\xi_{1}^% {2}}},$

which is justified by analogy with the notation

 $k=\frac{1}{\sqrt{1-\sum\xi^{2}}}$

featured in the Lorentz substitution.

In this case, and in light of the condition $-r=t$, the invariants (5) become:

 $0,\quad A=-k_{0}(r+\sum x\xi),\quad B=-k_{1}(r+\sum x\xi_{1}),\quad C=k_{0}k_{% 1}(1-\sum\xi\xi_{1}).$

Moreover, we notice that the following systems of quantities:

 $\begin{matrix}x,&y,&z,&-r=t,\\ k_{0}X_{1},&k_{0}Y_{1},&k_{0}Z_{1},&k_{0}T_{1},\\ k_{0}\xi,&k_{0}\eta,&k_{0}\zeta,&k_{0},\\ k_{1}\xi_{1},&k_{1}\eta_{1},&k_{1}\zeta_{1},&k_{1}\end{matrix}$

undergo the same linear substitutions when the transformations of the Lorentz group are applied to them. We are led thereby to put:

 \left.\begin{aligned} \displaystyle X_{1}&\displaystyle=x\frac{\alpha}{k_{0}}+% \xi\beta+\xi_{1}\frac{k_{1}}{k_{0}}\gamma,\\ \displaystyle Y_{1}&\displaystyle=y\frac{\alpha}{k_{0}}+\eta\beta+\eta_{1}% \frac{k_{1}}{k_{0}}\gamma,\\ \displaystyle Z_{1}&\displaystyle=z\frac{\alpha}{k_{0}}+\zeta\beta+\zeta_{1}% \frac{k_{1}}{k_{0}}\gamma,\\ \displaystyle T_{1}&\displaystyle=-r\frac{\alpha}{k_{0}}+\beta+\frac{k_{1}}{k_% {0}}\gamma.\end{aligned}\quad\right\} (9)

It is clear that if $\alpha$, $\beta$, $\gamma$ are invariants, $X_{1}$, $Y_{1}$, $Z_{1}$, $T_{1}$ will satisfy the fundamental condition, i.e., the Lorentz transformations will make them undergo an appropriate linear substitution.

However, for equations (9) to be compatible we must have

 $X_{1}\xi-T_{1}=0,$

which becomes, replacing $X_{1}$, $T_{1}$, $Z_{1}$, $T_{1}$ with their values in (9) and multiplying by $k_{0}^{2}$:

 $-A\alpha-\beta-C\gamma=0.$ (10)

What we would like is that the values of $X_{1}$, $Y_{1}$, $Z_{1}$ remain in line with Newton’s law when we neglect (as above) the squares of velocities $\xi$, etc. with respect to the square of the velocity of light, and the products of acceleration and distance.

We could select

 $\beta=0,\qquad\gamma=-\frac{A\alpha}{C}.$

To the adopted order of approximation, we obtain

 $k_{0}=k_{1}=1,\ C=1,\ A=-r_{1}+\sum x(\xi_{1}-\xi),\ B=-r_{1},\ x=x_{1}+\xi_{1% }t=x_{1}-\xi_{1}r.$

The $1^{\text{st}}$ equation in (9) then becomes

 $X_{1}=\alpha(x-A\xi_{1}).$

But if the square of $\xi$ is neglected, $A\xi_{1}$ can be replaced by $-r_{1}\xi_{1}$, or by $-r\xi_{1}$, which yields:

 $X_{1}=\alpha(x+\xi_{1}r)=\alpha x_{1}.$

Newton’s law would yield

 $X_{1}=-\frac{x_{1}}{r_{1}^{3}}.$

Consequently, we must select a value for the invariant $\alpha$ which reduces to $\displaystyle-\frac{1}{r_{1}^{3}}$ in the adopted order of approximation, that is, $\displaystyle\frac{1}{B^{3}}$. Equations (9) will become:

 \left.\begin{aligned} \displaystyle X_{1}&\displaystyle=\frac{x}{k_{0}B^{3}}-% \xi_{1}\frac{k_{1}}{k_{0}}\frac{A}{B^{3}C},\\ \displaystyle Y_{1}&\displaystyle=\frac{y}{k_{0}B^{3}}-\eta_{1}\frac{k_{1}}{k_% {0}}\frac{A}{B^{3}C},\\ \displaystyle Z_{1}&\displaystyle=\frac{z}{k_{0}B^{3}}-\zeta_{1}\frac{k_{1}}{k% _{0}}\frac{A}{B^{3}C},\\ \displaystyle T_{1}&\displaystyle=-\frac{r}{k_{0}B^{3}}-\frac{k_{1}}{k_{0}}% \frac{A}{B^{3}C}.\end{aligned}\quad\right\} (11)

We notice first that the corrected attraction is composed of two components: one parallel to the vector joining the positions of the two bodies, the other parallel to the velocity of the attracting body.

Remember that when we speak of the position or velocity of the attracting body, this refers to its position or velocity at the instant the gravitational wave takes off; for the attracted body, on the contrary, this refers to the position or velocity at the instant the gravitational wave arrives, assuming that this wave propagates with the velocity of light.

I believe it would be premature to seek to push the discussion of these formulas further; I will therefore confine myself to a few remarks.

1° The solutions (11) are not unique; we may, in fact, replace the the global factor $\displaystyle\frac{1}{B^{3}}$ by

 $\frac{1}{B^{3}}+(C-1)f_{1}\,(A,\ B,\ C)+(A-B)^{2}f_{2}\,(A,\ B,\ C),$

where $f_{1}$ and $f_{2}$ are arbitrary functions of $A$, $B$, $C$. Alternatively, we may forgo setting $\beta$ to zero, but add any complementary terms to $\alpha$, $\beta$, $\gamma$ that satisfy condition (10) and are of second order with respect to the $\xi$ for $\alpha$, and of first order for $\beta$ and $\gamma$.

2° The first equation in (11) may be written:

 $X_{1}=\frac{k_{1}}{B^{3}C}\left[x\left(1-\sum\xi\xi_{1}\right)+\xi_{1}\left(r+% \sum x\xi\right)\right]{}$ (11’)

and the quantity in brackets itself may be written:

 $(x+r\xi_{1})+\eta(\xi_{1}y-x\eta_{1})+\zeta(\xi_{1}z-x\zeta_{1}),$ (12)

such that the total force may be separated into three components corresponding to the three parentheses of expression (12); the first component is vaguely analogous to the mechanical force due to the electric field, the two others to the mechanical force due to the magnetic field; to extend the analogy I may, in light of the first remark, replace $\displaystyle\frac{1}{B^{3}}$ in (11) by $\displaystyle\frac{C}{B^{3}}$, in such a way that $X_{1}$, $Y_{1}$, $Z_{1}$ are linear functions of the attracted body’s velocity $\xi$, $\eta$, $\zeta$, since $C$ has vanished from the denominator of ($11^{\prime}$).

Next we put:

 \left.\begin{aligned} \displaystyle k_{1}(x+r\xi_{1})&\displaystyle=\lambda,&% \hskip 10.0pt\quad k_{1}(y+r\eta_{1})&\displaystyle=\mu,&\displaystyle\quad k_% {1}(z+r\zeta_{1})&\hskip 10.0pt=\nu,\\ \displaystyle k_{1}(\eta_{1}z-\zeta_{1}y)&\displaystyle=\lambda^{\prime},&% \hskip 10.0pt\quad k_{1}(\zeta_{1}x-\xi_{1}z)&\displaystyle=\mu^{\prime},&% \displaystyle\quad k_{1}(\xi_{1}y-x\eta_{1})&\hskip 10.0pt=\nu^{\prime};\\ \end{aligned}\right\} (13)

and since $C$ has vanished from the denominator of ($11^{\prime}$), it will follow that:

 \left.\begin{aligned} \displaystyle X_{1}&\displaystyle=\frac{\lambda}{B^{3}}-% \frac{\eta\nu^{\prime}-\zeta\mu^{\prime}}{B^{3}},\\ \displaystyle Y_{1}&\displaystyle=\frac{\mu}{B^{3}}-\frac{\zeta\lambda^{\prime% }-\xi\nu^{\prime}}{B^{3}},\\ \displaystyle Z_{1}&\displaystyle=\frac{\nu}{B^{3}}-\frac{\xi\mu^{\prime}-\eta% \lambda^{\prime}}{B^{3}};\end{aligned}\quad\right\} (14)

and we will have moreover:

 $B^{2}=\sum\lambda^{2}-\sum\lambda^{\prime}{}^{2}.$ (15)

Now $\lambda$, $\mu$, $\nu$, or $\displaystyle\frac{\lambda}{B^{3}}$, $\displaystyle\frac{\mu}{B^{3}}$, $\displaystyle\frac{\nu}{B^{3}}$, is an electric field of sorts, while $\lambda^{\prime}$, $\mu^{\prime}$, $\nu^{\prime}$, or rather $\displaystyle\frac{\lambda^{\prime}}{B^{3}}$, $\displaystyle\frac{\mu^{\prime}}{B^{3}}$, $\displaystyle\frac{\nu^{\prime}}{B^{3}}$ is a magnetic field of sorts.

3° The postulate of relativity would compel us to adopt solution (11), or solution (14), or any solution at all among those derived on the basis of the first remark. However, the first question to ask is whether or not these solutions are compatible with astronomical observations. The deviation from Newton’s law is of the order of $\xi^{2}$, i.e., 10000 times smaller than if it were of the order of $\xi$, i.e., if the propagation were to take place with the velocity of light, ceteris non mutatis; consequently, it is legitimate to hope that it will not be too large. To settle this question, however, would require an extended discussion.

Paris, July, 1905.

H. Poincaré