1We use the tensor calculus including the sum convention and restrict ourselves to Cartesian coordinates for simplicity; means the partial derivative .
In a previous paper  we have proposed in view of the dark matter problem a modification of the Newtonian gravity theory in such a way that the assumption of the existence of dark matter is not necessary. In detail we have given a theoretical explanation of MOND  following the structure of the electrodynamics assuming also in the case of Newtonian gravity a difference between intensive and extensive field quantities. The intensive field quantity described by obeys the homogeneous field equation1
where is the gravitational potential; simultaneously determines the gravitational force on a massive point particle (mass m) at the position
according to ( )
using the weak equivalence principle. Thus is the field strength. In consequence of (1) and (2) energy conservation is guaranteed for . On the other hand there exists the extensive field quantity determined by the mass density of the matter distribution according to the inhomogeneous field equation
(G Newtonian gravitational constant). Accordingly is the field excitation, because it is determined by the excitation Equation (3). Between both quantities the “material” equation is mediated:
with the “material” quantity , which may depend on the field strength itself in consequence e.g. of induced vacuum polarisations. In our previous paper we have assumed
where is a critical field strength, under which the value of increases drastically, so that the field strength increases compared with the Newtonian case ( ), and the assumption of dark matter will be superfluous.
In the present paper we show that the ansatz (5) can be enlarged in a very simple way in such a direction that the field strength can also not be larger than a critical field strength also in consequence of vacuum polarisations so that singularities (e. g. black holes) may be avoided. Doing this we enlarge relation (5) to
This procedure corresponds exactly to the idea of Born and Infeld in the electrodynamics  avoiding there electromagnetic field strengths larger than a critical value. Similar impressive influences of vacuum polarisations and fluctuations on the classical physics are e.g. the Casimir effect  and the Scharnhorst-Barton effect  .
1Usually Lenz’s rule is connected with time-variable electromagnetic fields. But it has a deeper meaning for all physical systems. If induced quantities, induced by time-variations or variable field strengths in space, would reinforce the original cause, the system would be unstable. Thus Lenz’s rule is rather a general stability law.
According to our proposal (6) we assume that the gravitational field strength will be weakened by increasing field-strength and enforced by decreasing field-strength in consequence e.g. of vacuum polarisations induced by the gravitational field-strength itself following Lenz’s rule. Herewith Lenz’s fundamental rule is implemented in the whole range of field strengths.2 The Newtonian gravity is valid in the very large intermediate range where is valid and where also the Newtonian gravitational constant G is determined.3 The value of amounts to 2 × 10−8 cm/sec2 in view of the flat rotation curves of the spiral galaxies (see discussions in  and footnote 5), whereas the value of may be of the order of 1030 cm/sec2, where the vacuum will be unstable e.g. with respect to spontaneous particle-antiparticle generations.4 Of course an exact theoretical derivation of the ansatz (6) does not exist until now because a quantum theory of gravity is missing. In a quantum theory of gravity vacuum polarisation effects would be included, so that a relation like (6), comparable with the running coupling constant in the non-Abelian QCD, would be expectable in the classical limit. In the electromagnetic case Heisenberg and Euler could show that the Born-Infeld ansatz follows from the quantum electrodynamics  . In the same sense we consider also our proposal as a phenomenological extension of gravity taking into account expectable quantum effects on the classical level.
2. The Integration Procedure
For solving the equation of motion (2) for a point-like test particle the knowledge of the field strength is necessary. For this we find from the field Equation (1) and Equation (3) with the use of (4) for the potential the differential equation:
Herein is given by (see (6))
3General Relativity would be also valid only in this intermediate range.
4By this acceleration an electron reaches approximately a distance of a Compton wave length during a Compton time necessary for spontaneous electron-positron pair creation.
Now we restrict ourselves for simplicity to the centrally symmetric case of a mass-sphere of mass M with constant density and radius R and consider at first the area outside the sphere, i.e. . Than the field Equation (7) reads
with the first integral:
(A integration constant). With
according to (8) we obtain from (10) the following quadratic equation for :
with the solution:
For we obtain
so that the flat rotation curves of the spiral galaxies are guaranteed with the constant orbital velocity .5 This is in accordance with our previous paper and confirms also the Tully-Fisher law  . The exact value of the constant A can be only determined by the connection of the solution (13) with the inner one for , see Equation (28).
Thus we consider now the case . With the abbreviation
we obtain from (7) the differential equation
which goes over in the centrally symmetric case into:
We solve this inhomogeneous differential equation by the method of the variation of the constants. The solution of the homogeneous equation belonging to (17) is already known and given by (10) as
5Herewith the value can be estimated, c.f. (28) and  .
where is now the integration constant, which is to vary for solving (17). This gives the differential equation for
with the solution
(B new integration constant), so that Equation (18) goes over into:
where is finally given by (11). Herewith we obtain once more a quadratic equation for with the solution:
Setting we obtain the solution (13), where B plays the role of A. Accordingly B has the meaning of
where is the mass value of an additional central point mass (see Equation (28)). For one gets for the case
Obviously this is the solution for an additional point mass at , where the field-strength possesses no singularity but reaches its maximum possible value for independently of the mass value . Therefore also the existence of black holes should be avoidable.
Avoiding in the following this central point mass we have to set . Herewith it follows from (22):
Now we obtain for after inserting of C according to (15) ( )
which means that the orbital velocity v has the value
This is greater than the Newtonian one by the factor . Therefore
the high star velocities near the center of the galaxy should be discussed newly avoiding the assumption of a black hole in the galaxy center.
Finally we have to determine the meaning of the integration constant A in the solution (13). It follows from the condition, that the solution (25) for goes steadily over into the solution (13) for at . On this way one finds immediately
Now we can give the complete solution for the massive sphere of radius R and mass M. It is valid for (see (25)):
and for (see (13)):
The square root of (29) and (30) is given by (32) and (33). In the case of the solution (29) results in the finite value (compare (24))
so that black holes may be avoided, although the field strength is very high. But for confirming this consequence finally a general relativistic investigation of the situation may be necessary.
3. Alternative Solution of the Centrally Symmetric Case
Considering the Equation (4) and Equation (6) one can solve these with respect to in the centrally symmetric case, i.e. with respect to . After that one can insert for the Newtonian expression according to (3). In the centrally symmetric case the condition (1) is fulfilled as one can test easily. In this way one finds:
For the mass-sphere with constant mass density and radius R one has to insert
Squaring of (32) results into (29) and (30).
We mention however explicitly that this procedure is in general only applicable in the centrally symmetric case because of the condition (1) and in the case, where the mass density distribution is not determined by the gravitational field strength itself. For the latter case, where is determined e.g. by the Euler equation
( pressure of the substratum) as e.g. in stars and galaxies or polytropic gas spheres the procedure of chapt. 2 is to be used.
Considering induced vacuum polarisations we could show in the framework of Newtonian gravity, that the assumption of dark matter and the existence of gravitational field singularities can be avoided. However our proposal discussed in this paper should be translated finally into a general relativistic form. This is not possible in an immediate way, because in the general theory of relativity the distinction between extensive and intensive field quantities is impossible. Therefore the modification of the theory has to start from the very beginning by a modification of the Lagrangian in the form of a -theory  as it is done also in the Born-Infeld electrodynamics. This will be our further aim.
 Milgrom, M. (1994) Dynamics with a Nonstandard Inertia-Acceleration Relation: An Alternative to Dark Matter in Galactic Systems. Annals of Physics, 229, 384-415.
Bekenstein, J.D. (2006) The Modified Newtonian Dynamics—MOND and Its Implications for New Physics. Contemporary Physics, 47, 387-483.
 Scharnhorst, K. and Barton, G. (1993) QED between Parallel Mirrors: Light Signals Faster Than c, or Amplified by the Vacuum. Journal of Physics A, 26, 2037.