Mayeul Arminjon

Charge conservation in a gravitational field in the scalar ether theory

A modification of the Maxwell equations due to the presence of a gravitational field was formerly proposed for a scalar theory with a preferred reference frame. With this modification, the electric charge is not conserved. The aim of the present work was to numerically assess the amount of charge production or destruction. We propose an asymptotic scheme for the electromagnetic field in a weak and slowly varying gravitational field. This scheme is valid independently of the theory and the "gravitationally-modified" Maxwell equations. Then we apply this scheme to plane waves and to a group of Hertzian dipoles in the scalar ether theory. The predicted amounts of charge production/destruction discard the formerly proposed gravitationallymodified Maxwell equations. The theoretical reason for that is the assumption that the total energy tensor is the sum of the energy tensor of the medium producing the electromagnetic (e.m.) field and the e.m. energy tensor. This means that an additional, "interaction" tensor has to be present. With this assumption, the standard Maxwell equations in a curved spacetime, which predict charge conservation, are compatible with the investigated theory. We find that the interaction energy might contribute to the dark matter.

Introduction and summary

Since the standard Maxwell equations apply in special relativity i.e. in a flat Minkowski spacetime, any theory of gravitation with a curved spacetime has to modify them, while ensuring that the modification reduces to the standard form in the particular case of flat spacetime. In a foregoing work [START_REF] Arminjon | Continuum dynamics and the electromagnetic field in the scalar ether theory of gravitation[END_REF], such a modification of the Maxwell equations in a gravitational field has been proposed for a scalar theory with a preferred reference frame, or (in short) with an "ether". (That theory is compatible with special relativity, however, and has passed a number of tests; see the Introduction of Ref. [START_REF] Arminjon | Continuum dynamics and the electromagnetic field in the scalar ether theory of gravitation[END_REF], where the motivations for the theory are also displayed.) It turns out that this modification leads to a violation of charge conservation if the gravitational field is time-dependent in the preferred frame [START_REF] Arminjon | Continuum dynamics and the electromagnetic field in the scalar ether theory of gravitation[END_REF]. However, the amount of that violation, which is the really important thing, could not be assessed at that stage.

The purpose of the present work, therefore, was to provide such assessment on a solid basis. To begin with, Section 2 summarizes the equations for the electromagnetic (e.m.) field in the theory, as given in Ref. [START_REF] Arminjon | Continuum dynamics and the electromagnetic field in the scalar ether theory of gravitation[END_REF]. Then Section 3 establishes the exact equations for charge production/destruction which result from those equations. In Section 4, we build an asymptotic approximation scheme for the electromagnetic field in a weak and slowly-varying gravitational field. This scheme has a general value, its main equations being valid for quite any theory of gravitation with a curved spacetime and the corresponding modification of the Maxwell equations. By using this scheme, in Sect. [START_REF] Arminjon | Space isotropy and weak equivalence principle in a scalar theory of gravity[END_REF] we obtain an explicit expression for the charge production rate that results from the investigated modification of the Maxwell equations in a gravitational field. We assess the time variation of the Newtonian potential and gravity acceleration, which enter that expression. We apply this to a plane wave in Sect. 6, and to a group of Hertzian dipoles in Sect. 7, thus getting figures. To do that, we follow a domain close to the source of the e.m. field in its motion through the "ether" of the theory and we integrate the charge production rate in small subdomains, using Lorentz transformations forth and back between the moving frame and the ether. It turns out that, for a reasonable velocity through the ether, the amounts thus assessed seem much too high to be compatible with the experimental facts.

In Sects. 8 and 9, we find the theoretical reason for this failure of the formerly proposed modification of the Maxwell equations in the investigated theory of gravitation: it is the assumption that the sum of two energy tensors, that of the charged medium and that of the e.m. field, is the total energy tensor which obeys the dynamical equation in a gravitational field. This as-sumption indeed leads to theoretical conclusions which cannot be generally true. Thus the only solution to this problem is to abandon this assumption, which means to introduce an additional, "interaction" tensor (Sect. 9). Then, the Maxwell equations in a gravitational field are not determined any more by the equation for continuum dynamics satisfied by the charged medium subjected to the Lorentz force. In particular, one may assume the standard gravitationally-modified Maxwell equations used in general relativity, which, we show in Appendix A, have a simple interpretation in the framework of the investigated theory. Section 9 ends by showing that, very generally, the gravitationally-modified Maxwell equations remain compatible with the geometrical optics of the theory, i.e., with the dynamics of a photon subjected to a gravitational and a non-gravitational force, which is detailed in Appendix B. Section 10 discusses the meaning of our findings. In particular, we suggest that the "interaction" tensor, which we find to be necessarily present according to the present theory, should at least contribute to the "missing mass" which has been invoked to explain motion at a galactic scale.

2 Main equations for the e.m. field in the scalar theory

In Ref. [START_REF] Arminjon | Continuum dynamics and the electromagnetic field in the scalar ether theory of gravitation[END_REF], the equations have been written in the Gauss system of units, for simplicity. However, due to the fact that in the Gauss units the Coulomb law is written as

F = q q r r 3 , (1)
in these units the charge dimension is

Q = (M L 3 T -2 ) 1 2
, where M , L and T represent mass, length, and time, respectively. This would introduce an undesirable coupling when writing weak-field asymptotic expansions, for in fact charge varies independently of mass, length, and time. Therefore, while presenting now the main equations [START_REF] Arminjon | Continuum dynamics and the electromagnetic field in the scalar ether theory of gravitation[END_REF], we rewrite them in the so-called SI units, i.e., in the "meter, kilogram, second, Ampère" (MKSA) system. The electromagnetic field is defined by the antisymmetric space-time tensor F : F µν = -F νµ , that obeys the standard first group of the Maxwell equations:

F λµ ,ν + F µν,λ + F νλ,µ = F λµ ;ν + F µν;λ + F νλ;µ = 0. ( 2
)
The expression of the Lorentz force is1

F i = q c F i 0 β + F i j v j c = q F i µ dx µ dt x . (3)
Here v j ≡ dx j /dt x (j = 1, 2, 3) is the velocity of the particle, measured with the local time t x , with

dt x dt = β(t, x), (4)
where t = x 0 /c is the coordinate time in a coordinate system (x µ ) adapted to (or bound with) the preferred reference fluid E assumed by the theory, and such that the synchronization condition γ 0j = 0 is true [START_REF] Arminjon | Continuum dynamics and the electromagnetic field in the scalar ether theory of gravitation[END_REF]; x ≡ (x i ); and (in any such coordinates)

β ≡ √ γ 00 . (5)
Equation ( 3) may be rewritten in space-vector form as

F = q (E + v ∧ B) , (a ∧ b) i ≡ e i jk a j b k , (6)
where the electric and magnetic vector fields are the spatial vector fields with components

E i ≡ c F i 0 β , B k ≡ - 1 2 e ijk F ij . (7)
In Eqs. ( 6) and ( 7), e ijk is the usual antisymmetric spatial tensor, its indices being raised or lowered using the spatial metric g in the preferred frame E; in spatial coordinate systems whose natural basis is direct, we have

e ijk = √ g ε ijk , e ijk = 1 √ g ε ijk , (8)
with ε ijk the signature of the permutation (i j k) and

g ≡ det(g ij ). ( 9
)
For a continuous charged medium, the electric charge density is defined as

ρ el ≡ δq/δV, ( 10
)
where δV is the volume element measured with the physical volume measure:

δV = √ g dx 1 dx 2 dx 3 . (11)
The 4-current is

J µ ≡ ρ el dx µ /dt x . (12)
The Lorentz force density is written, in accordance with (3), as

f i ≡ δF i δV = F i µ J µ . ( 13
)
Dynamics of a test particle is defined by an extension to curved spacetime of the special-relativistic form of Newton's second law [START_REF] Arminjon | On the extension of Newton's second law to theories of gravitation in curved space-time[END_REF]. This applies to the non-interacting particles that constitute a dust. The following dynamical equation for a continuous medium with velocity field v, subjected to a nongravitational external force density field f , has thus been derived for a dust from Newton's second law, and has been assumed to stay valid for a general continuum [START_REF] Arminjon | Continuum dynamics and the electromagnetic field in the scalar ether theory of gravitation[END_REF]:

T 0ν medium ;ν = b 0 (T medium ) + f .v cβ , T iν medium ;ν = b i (T medium ) + f i , ( 14
)
where T medium is the energy-momentum tensor of the continuous medium and b 0 (T ) ≡

1 2 γ 00 g ij,0 T ij , b i (T ) ≡ 1 2 g ij g jk,0 T 0k . ( 15
)
For a charged continuum, f i is given by Eq. ( 13), and the energy-momentum tensor T charged medium has to be substituted for T medium . Then Eq. ( 14) for this medium can be rewritten as:

T µν charged medium ;ν = b µ (T charged medium ) + F µ ν J ν . ( 16
)
It seems natural, almost obvious, to assume (i) that the total energy-momentum is the sum

T = T charged medium + T field , (17)
where T field is the energy-momentum tensor of the electromagnetic field [START_REF] Landau | The classical theory of fields, 3rd English edition[END_REF][START_REF] Fock | The theory of space, time and gravitation[END_REF]:

T µν field ≡ -F µ λ F νλ + 1 4 γ µν F λρ F λρ /µ 0 , (18) 5
or

T field ≡ (T µ field ν ) = F 2 - 1 4 (tr F 2 )I /µ 0 , F ≡ (F µ ν ); (19)
and (ii) that the total tensor T obeys the general equation ( 14) for continuum dynamics, without any non-gravitational external force:

T µν ;ν = b µ (T ). ( 20
)
These two assumptions have indeed been made in Ref. [START_REF] Arminjon | Continuum dynamics and the electromagnetic field in the scalar ether theory of gravitation[END_REF]. Combining Eqs. ( 16) and (20) using Eqs. ( 17) and ( 18), one derives

F µ λ F λν ;ν = µ 0 b µ (T field ) -F µ λ J λ , (21)
where b µ (T field ) is given by Eqs. ( 15) and (18). Under the two assumptions right above, this is the second group of the gravitationally-modified Maxwell equations in the investigated theory -at least for the generic case where the field tensor F is invertible (det

F ≡ det (F µ ν ) = 0). Indeed, if the matrix (F µ ν ) is invertible, (21) is equivalent to F µν ;ν = µ 0 (G µ ν b ν (T field ) -J µ ) , (G µ ν ) ≡ (F µ ν ) -1 . (22)
(Note that G, like F , is an antisymmetric tensor, G νµ = -G µν .)

Charge balance: exact equations

Due to the antisymmetry of F , we have F µν ;ν;µ = 0, and we get from Eq. (22

): ρ ≡ (J µ ) ;µ = (G µ ν b ν (T field )) ;µ . (23)
Thus, J µ ;µ = 0 -which, as is well known, means exact charge conservation -is not true in general, according to Eq. (22). To relate ρ with the charge production or destruction in some "substantial" domain Ω of the charged continuum, we note that, in coordinates adapted to (or bound with) this continuum C, we have

√ g C = √ -γ U 0 , (24)
with g C ≡ det(g ij C ), g C being the spatial metric in the reference fluid C, γ ≡ det(γ µν ), and U µ ≡ dx µ /ds the four-velocity field of the continuum.

(We used the well-known general relation [START_REF] Landau | The classical theory of fields, 3rd English edition[END_REF] γ = -γ 00 g.)

Therefore, noting ρ el ≡ δq/δV C the proper charge density, we have

d dt Ω δq = Ω ∂(ρ el √ g C ) ∂t d 3 x = c Ω (ρ el U µ ) ;µ √ -γ d 3 x. ( 26
)
The two equalities are true (only) in coordinates adapted to C. However, the rightmost integral is invariant under any change of the spatial coordinates, i.e., x i = ψ i (x 0 , (x j )). Also the time coordinate x 0 , with t ≡ x 0 /c, is arbitrary. Thus we have in any coordinates x µ :

d dt Ω δq = Ω ρ √ -γ d 3 x. ( 27
)
Of course the domain Ω as well as its boundary depend on t in general spatial coordinates x i .

Weak field approximation

We assume that the system of interest, S, e.g. the solar system or even only the Earth, has a weak and slowly varying gravitational field. To take benefit of this, we use a post-Newtonian (PN) approximation scheme that is in accordance with the general principles of asymptotic analysis. See Ref. [START_REF] Arminjon | Space isotropy and weak equivalence principle in a scalar theory of gravity[END_REF], Sect. V, and references therein. To do this, we conceptually associate with S a family (S λ ) of gravitating systems, depending smoothly on a parameter λ, so that the fields can be assumed to have asymptotic expansions as λ → 0. In fact we need only low-order Taylor expansions at λ = 0, which necessarily exist if the fields indeed depend smoothly on λ. Our assumption of a weak field then amounts to say that S corresponds to a small value of λ, say λ 0 1, so S = S λ 0 . Thus we can use the asymptotic expansions to approximate the values of the fields in that system by neglecting the remainder term. Each system is made of perfect fluids. The orders in λ of the corresponding fields: pressure, density, velocity, and the scalar gravitational field V ≡ -c 2 Logβ with β ≡ √ γ 00 , are the same as in a remarkable Newtonian similarity transformation [START_REF] Futamase | Newtonian and post-Newtonian approximations are asymptotic to general relativity[END_REF][START_REF] Arminjon | Asymptotic expansions for relativistic celestial mechanics[END_REF]. It follows that by adopting [M] λ = λ[M] and [T] λ = [T]/ √ λ as the new units for the system S λ (where [M] and [T] are the starting units of mass and time), all these fields become order λ 0 , and the small parameter λ is proportional to 1/c 2 ; indeed λ = (c 0 /c λ ) 2 , where c 0 is the velocity of light in the starting time unit [T], and c or c λ is with the time unit [T] λ . It thus becomes easy to derive asymptotic expansions. In particular, we have [START_REF] Arminjon | Space isotropy and weak equivalence principle in a scalar theory of gravity[END_REF]:

β ≡ √ γ 00 = 1 -U/c 2 + O(c -4 ). ( 28
)
where U is the zero-order term in the expansion of V . Moreover, the expansions are written in these λ-dependent units, in particular with the time T (the preferred time of the theory [START_REF] Arminjon | Continuum dynamics and the electromagnetic field in the scalar ether theory of gravitation[END_REF][START_REF] Arminjon | Space isotropy and weak equivalence principle in a scalar theory of gravity[END_REF]) being counted with the unit [T] λ for the system S λ , thus (28) should be written more precisely:

β (λ) (T, x) = 1 -U (T, x)/c 2 λ + O(c -4 λ ). ( 29
)
Note that U (T, x) is a coefficient in this asymptotic expansion and thus does not depend on λ. One may differentiate the PN expansions like (29) with respect to the corresponding time variable, here T ≡ x 0 /c (as well as with respect to the spatial coordinates x i ) [START_REF] Arminjon | Space isotropy and weak equivalence principle in a scalar theory of gravity[END_REF]. This automatically accounts for the slow time variation. By doing this in the assumed wave equation for V , we obtain that U in Eq. ( 28) or (29) obeys the Poisson equation with, on its r.h.s., the zero-order term in the expansion of the active mass density, hence U is the Newtonian potential [START_REF] Arminjon | Space isotropy and weak equivalence principle in a scalar theory of gravity[END_REF]. We obtain also that for the spatial metric, assumed in the theory to have the form

g = β -2 g 0 (30)
(with g 0 an invariable Euclidean metric), we have:

∂g ij ∂T = 2c -2 ∂ T U δ ij + O(c -4 ) (31)
(adopting Cartesian coordinates for the Euclidean metric g 0 until the end of Sect. 7).

Regarding now the electromagnetic field: we assume that the field tensor F , as well as the current 4-vector J, also depend smoothly on λ, hence they too admit low-order Taylor expansions at λ = 0 -but a priori we cannot say anything about the order of the main terms, hence we simply set (in the λ-dependent time and mass units, see above):

F = c n 0 F + c -2 1 F + O(c -4 ) (32)
and

J = c m 0 J + c -2 1 J + O(c -4 ) , (33)
for some integers n and m (positive, negative, or zero). Moreover, in contrast to the case with the gravitational field, we do not assume that F is slowly varying (nor weak, as witnessed by the arbitrariness of n). This means that the expansions (32)-( 33) are not PN ones but post-Minkowskian (PM) expansions. We developed an asymptotic PM approximation scheme [START_REF] Arminjon | Gravitational effects on light rays and binary pulsar energy loss in a scalar theory of gravity[END_REF] using the same method as the one we developed for the PN approximation [START_REF] Arminjon | Space isotropy and weak equivalence principle in a scalar theory of gravity[END_REF][START_REF] Arminjon | Asymptotic expansions for relativistic celestial mechanics[END_REF], in particular the time T is counted with the unit [T] λ for the system S λ as is the case with PN expansions, but the time variable now is x 0 = cT . 2One may differentiate PM expansions with respect to that new time variable. Therefore, in the modified Maxwell equation ( 21), the term F λν ;ν is of order c n as is also the term F µ λ . On the r.h.s., we note that µ 0 has dimension M LQ -2 , hence with the mass unit [M] λ = λ[M] we have µ 0 = µ 00 λ -1 , with µ 00 the value of µ 0 in the starting units. That is,

µ 0 = µ 00 c 2 (34)
if we take c 0 = 1 for the simplicity of writing. Thus the r.h.s. of ( 21) is of order c n+m+2 while the l.h.s. is of order c 2n , hence we must have

m + n + 2 = 2n, or m = n -2. (35)
From ( 15), (31), and (32), we get without difficulty

µ 0 b µ (T field ) = c 2n-3 ∂ T U × 0 T jj (µ = 0) or - 0 T 0i (µ = i) + O(c 2n-5 ), (36)
where, in accordance with (19), the matrix of the mixed components

0 T ≡ ( 0 T λ ν ) is 0 T = 0 F 2 - 1 4 (tr 0 F 2 )I. ( 37
)
By entering the expansions (32)-(33) of F and J into (21) using this, we get the lowest-order term in the weak-field expansion of (21) as

0 F µ λ 0 F λν ,ν = -µ 00 0 F µ λ 0 J λ . (38)
(We used also the identity

F λν ;ν = ( √ -γF λν ,ν )/ √ -γ [3]
and the fact that, from ( 25) and (30), we have

√ -γ = 1 + O(c -2 ).) (39)
Equation ( 38) is an exact equation for the expansion coefficient fields 0 F and 0 J, as is usual with an asymptotic expansion [START_REF] Arminjon | Asymptotic expansions for relativistic celestial mechanics[END_REF]. Hence, when the field matrix

0 F ≡ ( 0 F λ ν ) is invertible, that field is an exact solution of the flat-spacetime Maxwell equation: 0 F λν ,ν = -µ 00 0 J λ . ( 40
)
On the other hand, using (32), (36), and (39) in ( 23) gives us

ρ = c n-5 µ -1 00 0 G µ0 0 T jj - 0 G µi 0 T 0i ∂ T U ,µ + O(c n-7 ), (41)
where

0 G ≡ ( 0 G µ ν ) ≡ 0 F -1 . Due to (32)-(33), 0 F , 0 G, 0
T and 0 J do not have the physical dimensions of the corresponding fields F , G, etc.. Therefore, in accordance with (32)-(33), let us define

F 1 ≡ c n 0 F and J 1 ≡ c m 0 J. (42)
In view of (34) and (40), F 1 and J 1 are solutions of the flat-spacetime Maxwell equation with the correct dimensions in the SI units:

(F 1 ) λν ,ν = -µ 0 (J 1 ) λ . (43)
Obviously, the matrix F 1 is invertible when 0 F is, and its inverse is 34) and (37), the e.m. T-tensor T 1 associated with F 1 by ( 19) is

G 1 ≡ c -n 0 G. Also, from (
T 1 = µ -1 00 c 2n-2 0 T . ( 44
)
With this, we can rewrite (41) as

ρ = c -3 G 1 µ0 T 1 jj -G 1 µi T 1 0i ∂ T U ,µ + O(c -5 ). ( 45
)
Note that this is independent of the integers n and m in Eqs. (32)-(33).

Also note that, substituting (32) into the first group (2), one finds that 0 F and hence F 1 are also exact solutions of (2). To use the result (45) so as to assess the charge production that could be expected as a consequence of the modified second group (21), we have conversely to assume the following: Let F 1 and J 1 be solutions of the complete flat-spacetime Maxwell equations, i.e., the first group (2) and the second group (43). (Note that the first group (2) is the same independently of the presence or absence of a gravitational field.) In a given weak gravitational field, to any such pair (F 1 , J 1 ) there corresponds a unique solution (F , J) of the first group (2) and the modified second group (21), such that (F 1 , J 1 ) be the first terms in the PM expansion of (F , J), i.e., we have (42). This is expectable from perturbative arguments.

Explicit expressions

To use (45) we need explicit expressions, in terms of the components E i and B i of the electric and magnetic fields, for the field matrix F ≡ (F λ ν ), the corresponding T-tensor T , the inverse matrix G ≡ (G λ ν ) ≡ F -1 , and the 4-vector e with 3 e µ ≡ G µ0 T jj -G µi T 0i .

(We omit the index 1 used in (43)-(45).) We have

F ≡ (F µ ν ) =     0 E 1 /c E 2 /c E 3 /c E 1 /c 0 B 3 -B 2 E 2 /c -B 3 0 B 1 E 3 /c B 2 -B 1 0     . (47)
When detF = 0, the inverse matrix G = F -1 has a known expression for a general 4 × 4 matrix:

G = 1 detF 1 6 (tr F ) 3 -3tr F tr F 2 + 2tr F 3 I - 1 2 (tr F ) 2 -tr F 2 + F 2 tr F -F 3 . ( 48
)
3 Like in Sect. 4, the equations in this section are valid in spatial coordinates that are adapted to the frame E and, more specifically, are Cartesian for the Euclidean spatial metric g 0 ; and with the time coordinate being x 0 = cT , the preferred time coordinate of the theory. However, Eqs. (47)-( 52) are valid more generally in any Cartesian coordinates for the flat spacetime metric γ 0 that is built [START_REF] Arminjon | Space isotropy and weak equivalence principle in a scalar theory of gravity[END_REF] with g 0 and T . We use E i ≡ E i and B i ≡ B i , which means lowering the index i for E and B with the spatial, Euclidean part of metric γ 0 in such Cartesian coordinates.

Here, F has the following block structure:

F = 0 l l T a , ( 49
)
where l is a 1×3 matrix (row) and a is a 3×3 antisymmetric matrix. Clearly, tr F = 0. After computing F 3 in terms of this block structure, one finds that F 3 has all diagonal terms zero, so tr F 3 = 0, too. Using this in (47) and (48), with Matlab Symbolic Toolbox we obtain for G ≡ (G µν ) the following matrix:

    0 - B 1 c B 1 E 1 +B 2 E 2 +B 3 E 3 - B 2 c B 1 E 1 +B 2 E 2 +B 3 E 3 - B 3 c B 1 E 1 +B 2 E 2 +B 3 E 3 B 1 c B 1 E 1 +B 2 E 2 +B 3 E 3 0 E 3 B 1 E 1 +B 2 E 2 +B 3 E 3 - E 2 B 1 E 1 +B 2 E 2 +B 3 E 3 B 2 c B 1 E 1 +B 2 E 2 +B 3 E 3 - E 3 B 1 E 1 +B 2 E 2 +B 3 E 3 0 E 1 B 1 E 1 +B 2 E 2 +B 3 E 3 B 3 c B 1 E 1 +B 2 E 2 +B 3 E 3 E 2 B 1 E 1 +B 2 E 2 +B 3 E 3 - E 1 B 1 E 1 +B 2 E 2 +B 3 E 3 0     .
(50) (We checked that the product F .G is exactly the identity matrix I.) For the combinations (45) of the T-tensor (19), we get from (47) well-known results:

T jj = B 1 2 c 2 + B 2 2 c 2 + B 3 2 c 2 + E 1 2 + E 2 2 + E 3 2 2 c 2 µ 0 ≡ W field , (51)
T 0 i = -B 2 E 3 -B 3 E 2 c µ 0 B 1 E 3 -B 3 E 1 c µ 0 -B 1 E 2 -B 2 E 1 c µ 0 . ( 52
)
Putting this in (46) using (50) gives us

e µ =      0 B 1 3 c 2 +B 1 B 2 2 c 2 +B 1 B 3 2 c 2 +B 1 E 1 2 -B 1 E 2 2 -B 1 E 3 2 +2 B 2 E 1 E 2 +2 B 3 E 1 E 3 2 c µ 0 (B 1 E 1 +B 2 E 2 +B 3 E 3 ) B 1 2 B 2 c 2 +2 B 1 E 1 E 2 +B 2 3 c 2 +B 2 B 3 2 c 2 -B 2 E 1 2 +B 2 E 2 2 -B 2 E 3 2 +2 B 3 E 2 E 3 2 c µ 0 (B 1 E 1 +B 2 E 2 +B 3 E 3 ) B 1 2 B 3 c 2 +2 B 1 E 1 E 3 +B 2 2 B 3 c 2 +2 B 2 E 2 E 3 +B 3 3 c 2 -B 3 E 1 2 -B 3 E 2 2 +B 3 E 3 2 2 c µ 0 (B 1 E 1 +B 2 E 2 +B 3 E 3 )      . (53) Note that e 0 = 0 comes simply from the fact that B.(E ∧ B) = 0. Thus (45) rewrites as ρ = c -3 e i ∂ T U ,i + O(c -5 ). ( 54
)
To assess ∂ T U and ∂ T (∇U ), we note that these time derivatives must be evaluated in the preferred reference frame E assumed by the theory, and that the system of interest producing the e.m. field is expected to have a motion through that "ether", with a velocity field v whose modulus may be in the range 10-1000 km/s. In that situation, the main contribution to ∂ T U is that which is due to the mere translation motion of the relevant astronomical body -say the Earth -through E. Indeed, from the well-known integral giving U in terms of the Newtonian (i.e., zero-order PN) mass-energy density ρ, one gets exactly dU/dT ≡ ∂ T U + v.∇U = 0 (55)

if the system producing U has a rigid motion. In particular, it is exact for the self potential of a body whose motion is rigid. This is true for the Earth to a very good approximation, moreover the potential due to the Sun is nearly constant also on Earth; the most important departure from dU/dT = 0 should come from the Moon. For a rigidly rotating spherical body, we have moreover

v.∇U = V.∇U, (56)
with V ≡ ȧ, a(T ) being the center of the body. Thus to a good approximation (except for an unexpectedly small velocity V ), largely sufficient to get an order-of-magnitude estimate, we may consider the Earth as an isolated, rigidly moving, spherically symmetric body. For such a body we have:

∂ T U = -V.∇U = GM (r) r 2 V.e r ,
) 58
Let us compute ∂ T (∇U ), with ∇U = -GM (r) r 2 e r . To do that we may assume that a(T ) = VT . We then get

∂ T ∇U = -G V 2 T -V.x 4π ρ(r) r - 3M (r) r 4 e r - M (r)V r 3 . (59)
On the r.h.s., the two terms inside the large parentheses have the same order of magnitude as the last term inside the braces. Moreover those two terms cancel one another if the rotating spherical body is homogeneous. In that case we thus have:

∂ T ∇U = GM (r) r 3 V. ( 60
)
We shall use this approximation to get an order-of-magnitude estimate. On the Earth's surface, (57) gives 6 Case of a plane wave

∂ T U gV r 10V 10
Let us consider a general plane e.m. wave, whose propagation direction may be assumed parallel to the Cartesian basis vector i ≡ ∂ 1 :

E 1 = 0, E = E(x 1 ), B 1 = 0, B = B(x 1 ), cB = i ∧ E. ( 63
)
This is a solution of the flat Maxwell equations in vacuo. The electric and magnetic fields are orthogonal for such a wave. For the field matrix F ≡ (F µ ν ), most generally, the condition det F = 0 is equivalent to E.B ≡ g(E, B) = 0 [START_REF] Arminjon | Continuum dynamics and the electromagnetic field in the scalar ether theory of gravitation[END_REF]. It follows that F is not invertible for a plane wave. But the e.m. field ( Ẽ, B)

≡ (E + E , B + B ), (64)
with (E , B ) any constant e.m. field, is still a vacuum solution of the flat Maxwell equations depending only on x 1 -in fact, it is still a "plane wave" to the same extent as (63). For that field, generically, F is invertible. In that case, e i [Eq. ( 53)] is well defined and has e i ,i = 0, because (imposing E 1 = B 1 = 0) e 1 = 0 and e i = e i (x 1 ). Neglecting the term c -3 e i (∂ T U ) ,i in view of (62) and the extreme smallness of the c -3 factor, we thus get from (54) for the field (64):

ρ = 0 (Plane wave, c -3 e i (∂ T U ) ,i neglected). ( 65
)
Since this is independent of the constant e.m. field (E , B ), it remains true at the limit of a "pure" plane e.m. wave (63), for which cB = i∧E. However, the value of the neglected term does depend on the constant e.m. field (E , B ), and can be large if that field is very strong.

Case of a group of Hertzian dipoles

Hertz's famous oscillating dipole is the electric charge distribution

ρ el = T d,b,ω ≡ -e -iωt d.∇δ b (66)
(more exactly, ρ el is the real part of the r.h.s.). Here: b is the dipole's position; d is the dipole vector, and δ b is Dirac's measure at b. The associated 3-current is (the real part of the r.h.s. in):

j = -iωd e -iωt δ b , (67)
and the corresponding conserved 4-current is J ≡ (cρ el , j). The following electric and magnetic fields (or the associated field matrix (47)) provide an exact solution of the flat Maxwell equations in the distributional sense, for that current J:

4 E = α k 2 r (d -(r.d)r) cos ϕ + [3(r.d)r -d] cos ϕ r 3 + k sin ϕ r 2 , (68)
B = α k 2 (r ∧ d) cos ϕ r - sin ϕ kr 2 , k = ω c , ϕ ≡ kr -ωt. (69)
Here, α ≡ 1 4π 0 = 9 × 10 6 , α ≡ c We thus consider a group of dipoles that all are at rest in a common frame E V moving uniformly at V with respect to the preferred reference frame E. To calculate ρ in the vicinity of those moving dipoles, we consider successively each cube C in a regular mesh of small cubes at rest in E V . Each cube C is defined by |x i -a i | ≤ h/2, where x i (i = 1, 2, 3) are Cartesian coordinates for 4 This solution can be easily found in the literature, e.g. Jackson [START_REF] Jackson | Classical electrodynamics[END_REF], McDonald [START_REF] Mcdonald | Radiation in the near zone of a Hertzian dipole[END_REF], but whether or not it is an exact solution does not appear clearly. We can prove that it is an exact solution; here is a summary. In the Lorenz gauge, the vector potential A obeys the wave equation: A = µ 0 j ([9], Eq. (6.16)). For the case of a stationary source, the integral solution A simplifies to Eq. (9.3) in Ref. [START_REF] Jackson | Classical electrodynamics[END_REF]. The potential (9.3) rewrites exactly as (9.11) [START_REF] Jackson | Classical electrodynamics[END_REF]. When applied to the "Dirac current" (67), we may easily show that (9.11) rewrites also exactly as (9.13). In turn, (9.13) leads to the solution (68)-(69) [START_REF] Jackson | Classical electrodynamics[END_REF]. That (9.3) applies to the singular current (67), can be checked from the general formula A = µ 0 E S, with E the elementary solution of the wave operator that has support in the half-space x 0 ≥ 0 of R 4 -that is the distribution given by E, h = R 3 h(|x| , x) dx/(4π |x|), h ∈ D(R 4 ) [START_REF] Guichardet | Calcul intégral -maîtrise de mathématiques C2[END_REF]. That general formula for A is here applied with S the current (67): as a distribution acting on functions in D(R 4 ), this current must be more precisely defined as S ≡ -iωd U ⊗ δ b , with U the distribution on R associated with the function t → e -iωt .

the Euclidean metric g 0 that is the spatial part of the flat Minkowski metric γ 0 (Note 3) in the inertial frame E V . (We mean inertial in the Minkowski space (V, γ 0 ), with V the spacetime manifold. The uniformity of V is meant in the same sense.) Neglecting the O(c -5 ) remainder in (54), we have:

C E (T ) ρ(T, X) d 3 X = c -3 C E (T ) e i ∂ T U ,i d 3 X = c -3 ∂C E (T ) g 0 (e, n)∂ T U dS, ( 70
)
where C E (T ) is the cube C, as it appears at the time T in the frame E, the latter being endowed with Cartesian coordinates X i (i = 1, 2, 3) for the Euclidean metric g 0 (that is the spatial part of metric γ 0 in the inertial frame E); we denote X ≡ (X i ); n is the external normal (for the metric g 0 ) to the boundary ∂C E (T ), and dS is the Euclidean surface element. The Cartesian coordinates (X i ) = (X, Y, Z) and (x i ) = (x, y, z) that we consider are such that the Ox axis coincides with OX and is parallel to V, so we are in the conditions of a special Lorentz transformation L V from E to E V :

L V : t = γ V (T -V X/c 2 ), x = γ V (X -V T ), y = Y, z = Z, (71)
with t the inertial time in the frame E V and γ V the Lorentz factor.

Let x ≡ (x, y, z) be a given point of the domain C that is at rest in the frame E V . At the time T of the preferred inertial frame, that point corresponds in the frame E to a spatial position X ≡ (X, Y, Z) such that

(t, x) = L V (T, X). (72)
The unknowns are X and the value t of the time of E V . (Here Y = y and Z = z from (71), so the unknowns are t and X.) Equation (71) leads to

t = T γ V - V x c 2 . ( 73
)
Knowing t and x we compute the sum of the fields (68)-(69) from the different dipoles, which is thus got in the moving frame E V , say E , B . We then transform the field to the frame E by the inverse Lorentz transformation:

E 1 = E 1 , E 2 = γ V (E 2 + V B 3 ), E 3 = γ V (E 3 -V B 2 ), ( 74
) B 1 = B 1 , B 2 = γ V (B 2 -V E 3 /c 2 ), B 3 = γ V (B 3 + V E 2 /c 2 ), (75)
and hence can compute the components e i of e by Eq. (53). We can thus calculate the surface integral on the r.h.s. of (70). From (70), the value of the field at the center a of the cube C is then approximately5 ρ(T, a) c -3

∂C E (T ) e i n i ∂ T U dS/(h 3 γ V ), (76)
since h 3 γ V is the Euclidean volume of the rectangle parallelepiped C E (T ), deduced from the small cube C by Lorentz contraction.

For three dipoles with d = 100 nC.m, ν = 100 MHz (λ = 3 m), situated at λ from one another, we get fields with moduli E a few 10 5 V/m, B 15 T (in the moving frame). With V = 10 km/s, ρ(T, x) (counted in electrons per period per cubic meter) then has peaks at ≈ ±2 × 10 8 e/m 3 /period. (The peaks are very sharp and their values depend somewhat on the discretization. We also integrated in time and the values keep very high.) This seems untenable, even though the sign of the predicted charge production alternates in space. Therefore, the version proposed in Ref. [START_REF] Arminjon | Continuum dynamics and the electromagnetic field in the scalar ether theory of gravitation[END_REF] of the gravitationally-modified second group of Maxwell equations, Eq. ( 21) here, looks like being discarded.

The reason for the problem

So it seems that Eq. ( 21) is not the right Maxwell second group of the theory. Why does this happen? As explained in Sect. 2, Eq. ( 21) has been deduced [START_REF] Arminjon | Continuum dynamics and the electromagnetic field in the scalar ether theory of gravitation[END_REF], under two assumptions, from the general dynamical equation for a continuum subjected to a non-gravitational external force density field f i , Eq. ( 14) -when the latter is applied to the case that the continuum is a charged one and the external force density f i is the Lorentz force density [START_REF] Stephani | General relativity -an introduction to the theory of the gravitational field[END_REF]. Those two assumptions are:

(i) The total energy-momentum tensor is the sum T = T charged medium +T field .

(ii) The total tensor T obeys the general equation for continuum dynamics, without any non-gravitational force, Eq. (20).

As it has been noted in Ref. [START_REF] Arminjon | Continuum dynamics and the electromagnetic field in the scalar ether theory of gravitation[END_REF], Eq. ( 14) applied to the charged continuum, plus Assumptions (i) and (ii), lead straightforwardly, in view of the linearity of the dependence b µ = b µ (T ), to

T 0ν field ;ν = b 0 (T field ) - f .v cβ , T iν field ;ν = b i (T field ) -f i . ( 77
)
It is precisely this equation which, combined with Eq. ( 13) giving the Lorentz force density and with the expression (18) of the energy-momentum tensor, leads to the second group (21) [START_REF] Arminjon | Continuum dynamics and the electromagnetic field in the scalar ether theory of gravitation[END_REF]. Now we observe that Eq. ( 77) 2 has exactly the form of ( 14) 2 as applied not to the charged medium but to the electromagnetic "field continuum" itself, if and only if the density field of the non-gravitational external force on the field continuum is given by

f i field = -f i ≡ -f i charged medium . (78)
In addition, we observe that Eq. (77) 1 has exactly the form of ( 14) 1 as applied to the "field continuum", if and only if the velocity field of that continuum is well-defined and verifies the following relation:

f field .v field = -f .v ≡ -f charged medium .v charged medium , (79)
that is, with Eq. ( 78), if and only if

f charged medium . (v field -v charged medium ) = 0. ( 80
)
The force density f field should represent the reaction of the charged continuum to the Lorentz force exerted by the field continuum. Thus (78) means the opposition of action and reaction. Poincaré has shown that this does not apply to "matter" (the charged medium, e.g. a Hertz oscillator), but this was in the sense that matter emitting e.m. radiation does not conserve its momentum unless one counts also the momentum of the emitted e.m. field, i.e., the total momentum is in fact conserved [START_REF] Poincaré | La théorie de Lorentz et le principe de réaction[END_REF]. In the present case, the actio-reactio opposition suggested by Eq. (77) 2 concerns the charged medium on one hand and the e.m. field on the other hand, thus it means that there is no net force on the combined medium: "charged medium plus e.m. field" (except for gravitation). This is consistent with the conservation of the total momentum (in the absence of gravitation). 6 As to the velocity of the field continuum: there is one situation for which it can be naturally deduced from its energy-momentum tensor, namely the case that the energy-momentum tensor is a tensor product, i.e. has the bilinear form

T µν = V µ V ν . ( 81
)
Equation ( 81) is true for a dust made of ordinary matter (that is composed of non-interacting particles with non-zero rest mass, and that behaves as a perfect fluid with zero pressure): for such an "ordinary dust" we have

T µν = ρ * U µ U ν , ( 82
)
where ρ * ≥ 0 is the proper rest-mass density, and where

U µ ≡ dx µ /ds (83)
is the four-velocity field, so that (81) applies with V µ ≡ √ ρ * U µ . For ordinary dust we deduce from (83) that the coordinate 3-velocity is

u i ≡ dx i dt = dx i ds ds dt = c U i U 0 , (84)
hence from (82)

u i = c T 0i T 00 . ( 85
)
The latter relation is not generally true for a continuous medium, e.g. it does not apply to a perfect fluid if the pressure is not zero; simultaneously, in that case, Eq. ( 81) is not true either. But it is natural from the case of a dust to assume that (85) applies to a continuous medium when Eq. ( 81) is true. If one applies (85) to the energy-momentum tensor T field of an e.m. field for which Eq. ( 81) is true, then one finds from tr T field = 0 that the velocity v field defined with physical clocks, with components v i field = 1 β u i , has modulus c as determined with the physical space metric g ([1], Eq. ( 70)):

v 2 field ≡ g(v field , v field ) = c 2 . ( 86
)
In turn, this is consistent with the fact that when (81) is true we have a "dust of photons", which behaves so also dynamically [START_REF] Arminjon | Continuum dynamics and the electromagnetic field in the scalar ether theory of gravitation[END_REF] (see Sect. 9 below). Moreover, it can be proved that, in order that an e.m. field verify Eq. ( 81), it is necessary and sufficient that this be a "null field", i.e., that both invariants be zero [START_REF] Arminjon | Continuum dynamics and the electromagnetic field in the scalar ether theory of gravitation[END_REF][START_REF] Stephani | General relativity -an introduction to the theory of the gravitational field[END_REF]:

E.B = 0, E 2 = c 2 B 2 . ( 87
)
In a coordinate system adapted to E and such that, at the event X considered, we have β(X) = 1 and g ij (X) = δ ij , the tensor F has the form (47). It follows that T 0i is given by ( 52) [in which now the indices of the components E i and B i have been lowered with the "physical" spatial metric g], and we get from this by (85):

v field = c(E ∧ B)/ |E ∧ B| ≡ ck. ( 88
)
Thus, v field = c for an e.m. field whose T-tensor has the form (81). Surprisingly, for such a "null field", both v field = c and Eq. ( 80) are true, although v charged medium < c and even v charged medium c in usual conditions. To see this, remember that f ≡ f charged medium is the Lorentz force density, given by

f = ρ el (E + v ∧ B) (89)
(now we note again v ≡ v charged medium for shortness). We get from this and from (88):

f .v = ρ el E.v, f .v field = cρ el (v ∧ B) .k. (90)
By ( 87) and (88), we can take (at the event X) a spatial basis (i, j, k) which is orthonormal for g and such that

E = cBi, B = Bj. ( 91
) Setting v = v 1 i + v 2 j + v 3 k, we thus get: E.v = cBv 1 , v ∧ B = B(v 1 k -v 3 i), c (v ∧ B) .k = cBv 1 . (92)
That is, from (90), we have Eq. (80). Clearly, the foregoing proof depends in an essential way on the assumption that the e.m. field is a null field, so if that is not the case Eq. (80) has no reason to apply. In fact for a general e.m. field it is not clear at all how one should define its velocity field v field .

In summary: the equation of continuum dynamics [START_REF] Truesdell | Mechanical basis of diffusion[END_REF] applied to the charged continuum, together with Assumptions (i) and (ii) above, imply Eq. (77) -from which the discarded second group (21) follows. Equation (77) 2 has the form of the equation of continuum dynamics (14) 2 as applied to the "field continuum" having the energy-momentum tensor (18), with the force density on the field continuum being the opposite of the Lorentz force exerted on the charged continuum by the field continuum, Eq. (78). Moreover, Eq. ( 77) 1 also has the form of the equation of continuum dynamics (14) 1 as applied to the "field continuum", if the velocity of the field continuum is well defined and its projection on the direction of the Lorentz force is the same as that of the velocity of the charged continuum, Eq. (80). While Eq. ( 80) turns out to apply in the case of a "null field", it has no reason to be true for a general e.m. field, for which one does not even know how to define the velocity field v field . So Eq. (77) 1 has no reason to be true for a general e.m. field.

The latter is one conceptual reason, admittedly not very strong, why Eq. ( 21) is not the right Maxwell second group of the theory -in addition to the hard fact that it leads to charge production/destruction at high rates.

Recall that, at least for a dust, the equation of continuum dynamics ( 14) is derived from Newton's second law, hence it should apply to the charged medium. Thus to avoid Eq. (77) and the discarded second group (21), either Assumption (i) or Assumption (ii) has to be abandoned.

The solution of the problem

We can't leave Assumption (ii), because the concept of a "total" energymomentum tensor obeying Eq. ( 20) is necessary to the theory of gravitation. So it is Assumption (i) that has to be abandoned. This means that there must exist an additional energy-momentum-stress tensor, let us call it the "interaction tensor", such that the total tensor obeying Eq. ( 20) is given by

T = T charged medium + T field + T inter . (93)
We note that, in general, at an event X for which T charged medium = 0, we have also T field = 0, so that we are in the presence of a mixture (in the precise sense used in the theory of diffusion): the two constituents of that mixture are the charged medium and the e.m. field. It is then standard that indeed the effective energy-momentum-stress tensor of the mixture as a whole is not the sum of the energy-momentum-stress tensors of its constituents [START_REF] Truesdell | Mechanical basis of diffusion[END_REF][START_REF] Müller | A thermodynamic theory of mixtures of fluids[END_REF]. (It thus does not mean that there is an additional medium present beyond the charged medium and the e.m. field.) Now, given the necessary Assumption (ii) and the equation of continuum dynamics ( 14), Eq. ( 21) is equivalent to the opposite, i.e., to Assumption (i). This is another conceptual reason why Eq. ( 21) is not the right Maxwell second group of the theory.

With Eq. ( 93) replacing Eq. ( 17), it is clear that the equation for continuum dynamics [START_REF] Truesdell | Mechanical basis of diffusion[END_REF] applied to the charged medium, together with Assumption (ii), do not imply Eq. (77) any more; hence they do not determine the modified Maxwell second group any more. Therefore, the usual "gravitationally-modified" second group valid in GR and in the other "metric theories of gravitation":

F µν ;ν = -µ 0 J µ , (94)
is not precluded any more, as it was before in the investigated theory. We show in Appendix A that Eq. ( 94), as well as the first group (2), can be written in terms of the spatial metric and the local time in the synchronized preferred reference frame E, and then take nearly the usual form of the flatspacetime Maxwell equations for 3-vectors. Thus the standard 2nd group (94) is well compatible with the present theory. As is known, it leads to exact charge conservation. At the present stage, other forms of the modified Maxwell second group can not be excluded either, provided they would be found to lead to low-enough charge production in usual situations. We will now show how the compatibility with geometrical optics, which was proved in Sect. 6 of Ref. [START_REF] Arminjon | Continuum dynamics and the electromagnetic field in the scalar ether theory of gravitation[END_REF] with the discarded second group (21), holds true in a more general situation, with emphasis on the case of the standard second group (94).

The main modification to be made to the argument there, is that one needs to assume that an external force density f field is indeed acting on the field continuum, as the "reaction" of the charged medium to the Lorentz force exerted on it by the field continuum. For the link with geometrical optics i.e. with Newton's second law applied to individual photons, we have only to consider the case of a "null field", i.e. the case that the T-tensor (19) of the e.m. field has the form (81). In that case, as shown in Ref. [START_REF] Arminjon | Continuum dynamics and the electromagnetic field in the scalar ether theory of gravitation[END_REF], the spatial part (14) 2 of the continuum dynamical equation is equivalent to Newton's second law for a "substantial" volume element of the field continuum, in the form ( [START_REF] Arminjon | Continuum dynamics and the electromagnetic field in the scalar ether theory of gravitation[END_REF], Eq. ( 25) with f = f for a dust):

7 δF + δE c 2 g = 1 β D Dt δE c 2 v , (96)
where

δF ≡ f δV, δE ≡ T 0 0 δV, v ≡ 1 β dx dt ≡ 1 β u. ( 97
)
This equivalence is valid if the T-tensor has the form (81), independently of the nature of the continuum. 8 Therefore, as noted in Ref. [START_REF] Arminjon | Continuum dynamics and the electromagnetic field in the scalar ether theory of gravitation[END_REF], it is valid also for a null e.m. field, and this indeed behaves as a "dust of photons". In Ref.

[1], it was also proved that the time part ( 14) 1 of the continuum dynamical equation is equivalent to the energy equation. The latter is transposed to a continuum from the following form valid for a test particle [START_REF] Arminjon | Continuum dynamics and the electromagnetic field in the scalar ether theory of gravitation[END_REF]:

d (E β) dt = E ∂ β ∂ t + β 2 F.v, F.v ≡ g(F, v) ≡ g ij F i v j . (98)
While the proof of that equivalence in Ref. [START_REF] Arminjon | Continuum dynamics and the electromagnetic field in the scalar ether theory of gravitation[END_REF] was limited to the case without external force density, i.e. the case that f = 0 in (96), it is straightforward to extend it to show that (14) 1 is equivalent to

d (δE β) dt = δE ∂β ∂t + β 2 δF.v, (99)
in which δE and δF are defined by Eq. (97). 9

7 Here the spatial vector gravity acceleration g is given by [START_REF] Arminjon | Continuum dynamics and the electromagnetic field in the scalar ether theory of gravitation[END_REF] g ≡ -c 2 grad g β β

with (grad g β) i ≡ g ij β ,j where (g ij ) is the inverse matrix of matrix (g ij ), and D Dt is the relevant time derivative ensuring that Leibniz' rule for the derivation of a scalar product is verified [START_REF] Arminjon | Continuum dynamics and the electromagnetic field in the scalar ether theory of gravitation[END_REF][START_REF] Arminjon | On the extension of Newton's second law to theories of gravitation in curved space-time[END_REF]. In terms of the T-tensor, the coordinate velocity u is given by (85). 8 This was noted in Ref. [START_REF] Arminjon | Continuum dynamics and the electromagnetic field in the scalar ether theory of gravitation[END_REF], Sect. 6, though in the case that the external force density field is f = 0. However, the proof given in Ref. [START_REF] Arminjon | Continuum dynamics and the electromagnetic field in the scalar ether theory of gravitation[END_REF], Sect. 3.2, applies with the definitions ρ ≡ T 00 /( β2 c 2 ) and u i ≡ cT i0 /T 00 to a medium for which the T-tensor has the form (81), independently of its nature, provided there is no force internal to this continuum, f = f . 9 By the way, it is not difficult to check that Eqs. ( 96)-(99) are covariant under any change (110), thus including any change x 0 = φ(x 0 ), even though for (99) this has not been noted in Ref. [START_REF] Arminjon | Continuum dynamics and the electromagnetic field in the scalar ether theory of gravitation[END_REF].

Thus, when the T-tensor of the e.m. field has the form (81), the dynamics of a volume element that is followed in its motion is just the same as that of an individual photon subjected to the gravitation and to an external force (detailed in Appendix B), and it may equivalently be defined by the dynamical equation ( 14), as applied to the field continuum:

T 0ν field ;ν = b 0 (T field ) + f field .v field cβ , T iν field ;ν = b i (T field ) + f i field . ( 100
)
Note that in the latter equation we have always, in view of (19):

µ 0 T µν field ;ν = -F µ λ F νλ ;ν . ( 101
)
(This is Eq. (61) of Ref. [START_REF] Arminjon | Continuum dynamics and the electromagnetic field in the scalar ether theory of gravitation[END_REF], rewritten in the MKSA system.) Therefore, with the standard second group (94), we get the well-known equation

T µν field ;ν = -F µ λ J λ . (102)
In vacuo (J µ = 0), Eqs. ( 100) and (102) give us:

0 = T 0ν field ;ν = b 0 (T field ) + f field .v field cβ , 0 = T iν field ;ν = b i (T field ) + f i field .
(103) This shows that there must indeed be an external force acting on the photon dust (the "reaction" to the Lorentz force), in addition to the gravitation.

Discussion

The main conclusions of this work are as follows:

(i) The formerly proposed modification of Maxwell's second group in a gravitational field [START_REF] Arminjon | Continuum dynamics and the electromagnetic field in the scalar ether theory of gravitation[END_REF] in the investigated theory predicts unrealistically high rates of production/destruction of electric charge. Therefore, that gravitationally-modified Maxwell 2nd group has to be discarded.

(ii) The theoretical reason for that is the former assumption [START_REF] Arminjon | Continuum dynamics and the electromagnetic field in the scalar ether theory of gravitation[END_REF] according to which the total energy tensor which obeys the dynamical equation in a gravitational field (20), is the sum of the energy tensor of the charged medium and that of the e.m. field, Eq. ( 17). This assumption is not consistent with the fact that these two media form a mixture and, in addition, has a consequence which has no reason to be verified in general.

(iii) Therefore, one must assume an additional, "interaction" energymomentum tensor, such that Eq. ( 17) is replaced by Eq. ( 93). With this, Maxwell's second group in a gravitational field is less constrained, in particular the standard version (94) valid in GR becomes well compatible with the investigated theory.

(iv) Also, one must assume that, at least for a null e.m. field (which indeed can be considered as a continuous medium with a well-defined velocity field, and to which one may definitely apply Newton's second law of the present theory), there is a reaction force exerted on the e.m. field by the charged medium. Equation (93) means that the presence of usual matter producing an e.m. field necessarily gives rise (according to the present theory) to the presence of an additional kind of energy, with energy tensor T inter . The latter does not generally vanish outside the charged medium that emits the e.m. field. If the standard 2nd group (94) is assumed, we get immediately from ( 16), (20), and (102):

T µν inter ;ν = b µ (T field ) + b µ (T inter ). ( 104
)
Without the b µ (T field ) term, this equation would be identical with the dynamical equation in a gravitational field (20), with the energy tensor T inter in the place of the total tensor T . The time component of Eq. (104) rewrites as:10 T 00 inter ,0 + T 0j inter ,j = (Log β) ,0 T 00 inter + b 0 (T field ),

or

T 00 inter ,0 + T 0j inter ,j = (Log β) ,0 T 00 inter -β -4 W field . (106)
Each of the two source terms on the r.h.s. is proportional to the variation of the gravitational field in the preferred reference frame: (Log β) ,0 -c -3 ∂ T U in a weak field. However, the first term (with T 00 inter ) expresses the usual energy conservation in the investigated theory: the conservation of the total energy has just the same form as (106), without the term involving W field and with the total tensor T in the place of T inter . That usual energy conservation is a balance between matter energy and gravitational energy [START_REF] Arminjon | Space isotropy and weak equivalence principle in a scalar theory of gravity[END_REF]. Thus, one may consider that the source of the interaction energy really is the second term that is proportional to the e.m. energy W field , Eq. ( 51). The e.m. energy should have grossly an ellipsoidal distribution around the center of a galaxy, which should be much less flat than the mass distribution of the luminous objects, since the e.m. energy decreases as 1/r 2 from each of these objects. Hence we may a priori expect that the same applies to the interaction energy T 00 inter . As is well known, to explain the motions at a galactic scale one is led to assume a distribution of unseen, "dark" matter, which should fill a "halo" around the center of a galaxy or a cluster of galaxies. Therefore, it seems natural to conjecture that the interaction energy, which is necessarily present according to the present theory, is a contribution to dark matter. This is worth a further investigation.

A Meaning of the standard gravitationallymodified Maxwell equations

First, let us observe that, with the definitions (7) for the electric and magnetic fields, Maxwell's first group (2) can be rewritten almost exactly in the usual form of the Maxwell-Gauss and Maxwell-Faraday equations, namely:

div g B = 0, rot g E = - ∂B ∂t x , (107)
the mere difference being thus in the use of the local time (4) (i.e., ∂ ∂tx = 1 β(t,x) ∂ ∂t ) of the synchronized reference frame E (i.e., γ 0i = 0) and in the fact that the operators are defined with the help of the spatial metric g in the frame E:

div g B ≡ B i |i = 1 √ g √ gB i ,i , (rot g E) i ≡ e ij k E k |j , (108)
with |j the covariant derivative associated with g. Indeed, as one easily checks from the definition (7), Eq. (107) coincides with (2) in coordinates x µ such that, at the event X considered, we have

g ij (X) = δ ij , g ij,k (X) = 0, β(X) = 1 (109)
(see Eq. (24.14) of Fock [START_REF] Fock | The theory of space, time and gravitation[END_REF]). Starting from one coordinate system that is adapted to the preferred frame and that verifies the synchronization condition γ 0i = 0, one can get to another one that in addition verifies (109), by a change

x i = ψ i (x 1 , x 2 , x 3 ), x 0 = φ(x 0 ). ( 110
)
Since each of the two equations in (107) is invariant under such a change, our statement is proved. In nearly the same way, from the relation valid in any coordinates for the antisymmetric tensor F µν :

F µν ;ν = 1 √ -γ √ -γF µν ,ν , (111)
and using Eq. ( 25) for γ, we find that the µ = 0 component of Eq. ( 94) [the standardly-modified second group] rewrites as the Maxwell-Poisson equation in terms of metric g:

div g E = µ 0 c 2 ρ el ≡ ρ el 0 , (112)
whereas the spatial components of Eq. ( 94) rewrite as a space vector equation involving an additional term as compared with the flat-spacetime Maxwell-Ampère equation:

rot g B - 1 c 2 ∂E ∂t x - 1 c 2 g ∧ B = µ 0 j, (113)
where j i ≡ J i and the spatial vector g is given by Eq. (95). To rewrite (94) as ( 112) and (113), we use coordinates that, in addition to (109), are such that (∂β/∂x 0 )(X) = 0. (114)

The full set (109) and (114) can be fulfilled by a change (110); cf. Ref.

[3], around Eq. (85.18). Equations ( 112) and (113) are invariant under a change (110). Note that the derivation applies in any synchronized reference frame, but the gravity acceleration vector (95) makes little sense in a general situation unless one assumes the preferred-frame dynamics of the investigated theory.

B Dynamics of a photon under gravitational and non-gravitational forces

Our extension of Newton's second law has exactly the same form for a mass particle and for a photon, i.e. [START_REF] Arminjon | Continuum dynamics and the electromagnetic field in the scalar ether theory of gravitation[END_REF][START_REF] Arminjon | On the extension of Newton's second law to theories of gravitation in curved space-time[END_REF][START_REF] Arminjon | Energy and equations of motion in a tentative theory of gravity with a privileged reference frame[END_REF]]

F + E c 2 g = 1 c 2 D(Ev) Dt x ≡ 1 c 2 β D(Ev) Dt . ( 115
)
For a photon, we define E = hν, h being Planck's constant and ν the frequency as measured with the local time: ν ≡ dn/dt x ≡ (1/β)dn/dt with n the number of periods. The energy equation derived from (115) is also the same for a mass particle and for a photon, i.e., Eq. ( 98). In the case without an external force F, the proof has been given in full for a mass particle in Ref. [START_REF] Arminjon | Energy and equations of motion in a tentative theory of gravity with a privileged reference frame[END_REF], and has been outlined also in Ref. [START_REF] Arminjon | Energy and equations of motion in a tentative theory of gravity with a privileged reference frame[END_REF] for a photon. Here we give the proof with F for a photon, for completeness. Equation ( 115) is equivalent to:

Eg + c 2 F = E Dv Dt x + dE dt x v, (116)
whence by taking the scalar product g with v, using Leibniz' rule verified [START_REF] Arminjon | On the extension of Newton's second law to theories of gravitation in curved space-time[END_REF][START_REF] Arminjon | Energy and equations of motion in a tentative theory of gravity with a privileged reference frame[END_REF] by the D/Dt derivative:

Eg.v = E d dt x v 2 2 + dE dt x v 2 -c 2 F.v. ( 117
)
Since v 2 = c 2 for a photon and since v ≡ (1/β)u, this rewrites as

E β g.u = c 2 dE dt x -F.v . (118)
From (95), we have

g.u = -c 2 g(grad g β, u) β = - c 2 β β ,i dx i dt = - c 2 β dβ dt - ∂β ∂t . (119)
Hence, (118) is

-E dβ dt - ∂β ∂t = β dE dt -β 2 F.v, (120)
whence follows Eq. (98).

Thus, for a photon, to deduce the energy equation (98) from Newton's second law (115), we used the property v 2 = c 2 . Now we show that the energy equation allows us to rewrite Newton's second law in a form which ensures that v 2 = c 2 is indeed maintained at all times. Equation ( 98) is equivalent to

β dE dt + E ∂β ∂t + β ,i u i = E ∂β ∂t + β 2 F.v. (121)
Removing the term present on both sides and using again (119), this gives us

dE dt = β F.v + E c 2 g.v . ( 122
)
With this, Newton's second law (116) rewrites as

E Dv Dt x + F.v + E c 2 g.v v = Eg + c 2 F, (123)
or Dv Dt x = g - g.v c 2 + F.v E v + c 2 E F. (124)
It follows from the latter equation that

d dt x v 2 2 = v. Dv Dt x = (c 2 -v 2 ) g.v c 2 + F.v E . (125)
Note that all of this is true for a mass particle as well as for a photon. Thus, if we have v 2 = c 2 at the initial time, this condition is maintained at all times. In general relativity, the condition v 2 = c 2 (i.e., ds 2 = 0) for a photon is not dynamically implied by the (geodesic) law of motion, instead it is assumed from the outset as one considers a "null geodesic".

0 u 2

 02 r ≡ |x -a(T )| , e r ≡ (x -a(T ))/r, (57) where M (r) ≡ 4π r ρ(u) du. (

  5 (MKSA) for V = 10 km/s (10 4 MKSA), (61) and (60) gives ∂ T ∇U gV/R, |∂ T ∇U | 1.5 × 10 -2 (MKSA) for V = 10 km/s. (62)

4π 2 .

 2 39 × 10 7 (MKSA). Thus, E.B = 0 since E is in the plane containing d and r ≡ (x-b)/ |x -b|, while B ∝ r∧d. However, we may consider a group of dipoles with different b's and d's, and by adding the corresponding fields (68)-(69) we get exact solutions of the flat Maxwell equations for which, generically, we have E.B = 0.

Greek indices vary from 0 to 3, Latin ones from 1 to 3 (spatial indices). Semi-colon means covariant derivative using the connection associated with the "physical" (curved) space-time metric, the latter being denoted by γ. Indices are raised and lowered with the help of this metric, unless explicitly mentioned otherwise.

Let S be a self-gravitating system of perfect fluids producing a weak gravitational field. The PM approximation associates with S a family (S λ ) of systems for which the velocity field u (λ) is order zero in λ, not order √ λ, hence it does not fit with the Newtonian limit. Nevertheless, it can be applied to a given such physical system S even if the fields really are slowly time-varying: the latter assumption just is not used. See the discussion in Ref.[START_REF] Arminjon | Gravitational effects on light rays and binary pulsar energy loss in a scalar theory of gravity[END_REF], §4.1. However, in the case of an e.m. wave, we expect that it really is not varying slowly, of course.

We approximate the integral of w.n (with w a spatial vector field, here w j = e j ∂ T U ) over an i face of the rectangle parallelepiped C E (T ) as I ±hh (w i (S 1 ) + ... + w i (S 4 ))/4, where h and h are the rectangular face's sides and S k (k = 1, ..., 4) are the vertices of that face. This algorithm has been tested by applying it to simple vector fields, e.g. w(x, y, z) = (-xy, xy, z 2 ). As expected, we found that ∂C w i n i dS/h 3 , thus calculated, approximates closely the value of div w = w j ,j at the centre a of a small cube C with side h.

See §III.C in Ref.[START_REF] Arminjon | Space isotropy and weak equivalence principle in a scalar theory of gravity[END_REF] for a discussion of the total momentum in the presence of gravitation in the present theory.

One sees that by using Eqs. (23) and (25) in Ref.[START_REF] Arminjon | Space isotropy and weak equivalence principle in a scalar theory of gravity[END_REF] plus Eqs. (15) and (30) above. The second step uses also Eq. (51).