Appraisal of a New Gravitational Constant
51A, S.P. 57, Accesso a M., 03017 Morolo ItalyAbstract
The need of extending the theory of relativity has led M.Tailherer to the hypothesis of a new fundamental equation and constant, embodying in a unique wave equation for the graviton the link between gradient of curvature and deformation of metric. As direct continuation of a preceding work, here a new assessment of the constant S in the Vortex Theory of gravitation is given in a more direct approach than 1st approximation yielding S =(2.5±1.2)E-19 m-1. Issues are concerned fitting by Maple four binary systems data, also allowing to assign a meaningful inertial mass to the graviton (5.5±2.6)E-61 Kg confirming known heuristic bounding. In Appendix an easy way of getting the vortex’s gradient formula is shown along with the whole action of the model and the description of the tide effect on a test mass with respect to a x polarized gravitational wave in the case of an asymmetric source.
Keywords: Alternative theories to General Relativity (TGR), Gravitational waves, Vorticity Tensor, Binary Star Systems, massive gravity. PACS. 04.30
Received May 05, 2015; Revised May 08, 2015; Accepted May 21, 2015
Copyright © 2015 Science and Education Publishing. All Rights Reserved.Cite this article:
- Sandro Antonelli. Appraisal of a New Gravitational Constant. International Journal of Physics. Vol. 3, No. 4, 2015, pp 139-149. https://pubs.sciepub.com/ijp/3/4/1
- Antonelli, Sandro. "Appraisal of a New Gravitational Constant." International Journal of Physics 3.4 (2015): 139-149.
- Antonelli, S. (2015). Appraisal of a New Gravitational Constant. International Journal of Physics, 3(4), 139-149.
- Antonelli, Sandro. "Appraisal of a New Gravitational Constant." International Journal of Physics 3, no. 4 (2015): 139-149.
Import into BibTeX | Import into EndNote | Import into RefMan | Import into RefWorks |
At a glance: Figures
1. Introduction
It is well-known as TGR is related to some compelling questions leading towards ambiguities [3]: does not it seem a contradiction indeed that a body subject to mere gravitational fields travel a geodesic of space-time losing in the meanwhile gravitational energy? Furthermore, the expression of gravitational damping in the known formula of third derivative of quadrupole moment (cf.eq.(7)), although proven valid by now also for strong fields by the post-Newtonian formalism with so great an approximation [18], on the other hand breaks down the general covariance in favour of linear one, just because the PN (v/c)2 series expansion cannot avoid the lost of the full symmetry of the exact Einstein’s Equations. This is demonstrated by the fact that in TGR any scalar quantity has a physical content only if has an invariant character with respect to general transformations of coordinates, but not so gravitational energy fluxed away, seen that it is always possible to choose a suitable frame for which two gravitating bodies are at rest [2].
Hence I think that from such dilemmas one cannot get out without taking the general relativity as only an aspect of an extended theory as proposed recently by the intriguing tensor Vortex gravity model [6], any advance being subordinate to the determination of the constant S, a very cumbersome goal related to the evaluation of the following delayed integro-differential equation, nothing but the solution{1} at great distance of the double curl in curved Space-Time of the main equation recalled shortly in Appendix B (see also [4, 5, 6]):
(1) |
with and where
is the energy momentum tensor (EM) of two stars masses (in our case the system PSR 1913+16 {2}) and the curvature Riemann tensor, the slash standing for contravariant derivative.
We can well distinguish in eq. (1) two terms whose the first, bound recursively to LHS , the deformation of metric with respect to proper time, proved at a great extent pre-eminent over the second one, (that is the double curl of the skew-symmetric tensor we called C Riemann tensor ) by tens of order of magnitude both for the former (see [1]) as for the latter solution by successive approximations, even though it has been evidenced that the contribution of the first term dropped in a clearly way in this second approximation. The reader is strongly recommended to go through the two papers [1, 6].
2. Methods
Reviewing the antecedents in the previous calculus, as reported in [1], we chose the Minkowski metric as start metric in the successive approximations procedure, i.e. a constant cartesian metric But already to the second approximation in which the new metric is taken out integrating numerically with respect to the proper time
(2) |
the calculus proved incredibly CPU-consuming: the chief term in Eq.(1) yields elements of the deformation tensor as large as 125 Gbytes each, unmanageable for any computing software. The choice has been thereby to simplify as much as possible to getting the 2nd order solution renouncing to covariant derivatives in the new metric concerning the term including the double curl (the term in Eq.(1)). Worth of noting the first term in Eq.(1), encompassing the curvature tensor is therefore really a crucial one and together the most burdensome needing hundreds of GBytes. Thus some kind of ansatz is required that should be chosen with insight.
Because of the way which is has been got the C tensor (see footnote 1) we have attempted to value the curvature tensor as to first approximation. Once evaluated the new through integration of according to Eq.(2), as done in the paper [1], we held at this stage the new Riemann tensor not to depend explicitly on the factors k nor on S because of the Christoffel symbols that it is composed of, which involves terms of the form However this choice is only of little interest in expressing the mean energy loss through the semi-axis and the simple 1/S2 dependence, so an efficacious value, but it makes absolutely unreliable the physical estimate of the constant S as already seen in [1] of order of 1014 m-1, while a just simple model to the second order has given the more acceptable result of S=0.3261 m-1 which has allowed already to attribute the graviton the mass 10-43Kg. Therefore, in effect we should expect some hinging on the ratio k/S. Taking as working hypothesis for the new Riemann tensor the more appealing makeshift expanding
(3) |
improves the estimate of S to 3.4E-14 m-1 for the three binary at issue, a value not far from the outcome of the complete direct calculus, still approximated, that we are going to carry out. Even though the Riemann tensor is denoted by the apex (2) the solution is always considered at 1st order because all the tensors are raised/lowered by the flat metric. We get for this model (see the whole power at https://1drv.ms/1NOnDU6 ; there “asse” stands for major semi-axis, “th” for colatitudes, “f” for orbital phase)
= -7.891E24 J/s against -7.959E24 J/s for Einstein’s and = -2.423E-12 to within 0.7% with the experimental one against [13,sect. 4.1.2]- as seen later too:
(4) |
In the following we are referring to the model of Eq.(3) as 1st approximation Vortex Model. As checked with CILEA supercomputer section the direct calculus of energy loss via equ.(1) is very beyond standard RAM performance (see footnote 3) without taking in account the lack of manageability of the output files to average in turn. The strategy that has been adopted here therefore was to guess that the passage of integration of Eq.(2) results in a staircase-like behaviour for each of the components of that can be well interpolated by the straight line of the kind being the orbital phase.
This fetch is put on to obtainwhile fora 3 steps Simpson algorithm has been adopted. Moreover, taking the direction of propagation of the wave from coordinate system of the orbit plane, as can be verified heuristically on data also from [1] concerning B1913+16, it is a good approximation to put d , especially to reduce the burden of calculation (especially introduced by the inverse ).This may be justified by seeing that the distribution function may be roughly shaped as two triangles each centered about a unit interval on . Then, since it is noted the differential energy loss not to depend significantly on the colatitude has been fixed to(edge-on line of sight). We have over-marked E because the power should be indeed averaged over time, so a factor to multiply is needed to turn back the average on the orbital anomaly, namely
(5) |
with and
The term dividing in Eq.(5) is omitted however in that the integration over the anomaly cancels it because the asymptotic metric to manage the strength tensor of Vortex theory (see later) is the usual Minkowski’s which does not depend on .Thus in the following plots for it is meant the integrand of Eq.(5) that we could style as “effective” differential power. Besides, the wrinkle has been preferred to focus the reckoning in correspondence with each one of the two main resonance peaks wherein to fit the two values of S. Actually, if one proceeded the whole calculus with the value of S sought for the 1st peak, the second one on the right would be found unequal of several orders, just because of the very rough but necessary simplification in getting the new metric in Eq.(2) through the numeric integral (as said, only one step per orbital interval in the azimuthal variable{3}) of deformation tensor carried out in such an approximation, inasmuch as the new metric becomes less accurate as the interval of integration extends over the whole angular orbital journey, so occasioning a greater growth in systematic errors. For this purpose the interval has been broken in the two and in each of which has been chosen to calculate by Maple the effective power function with respect to the anomaly in steps of 0.1-0,2 rad in the orbital plane. In short, this trick saves from a awkward integral over Worth to note is the fact that this does not imply loss of generality in that is seen that in Vortex Theory there is not any uncoupling of the from the mode in the orbital plane line of sight as in TGR (see probl.3.2 in [13]) although the mode here is differently defined. But let us just recall now the main steps to get the expression for dE/dt. We have applied firstly this method to the classical B1913+16 binary system while comparing the results with the other more recent B1534+12 [17], J0737-3039 [16] and B2127+11C[19] discovered in 1991 , 2003 and 1989. In carrying out the strength tensor (see Appendix B and [6]) for the energy loss calculus{4} with, the azimuthal variable has to be taken as a function of delayed time and so also of distance from the observer, i.e. which raises of a great deal the burden of calculus. This has led to a function of S f(S) to be inserted in the differential relation of third Kepler's law outputting a polynomial equation of which to find the root:
(6) |
Fortunately, the function f(S) , being monotonic, has been able to be tackled through the dichotomic method in figuring out the root (other methods as secants or Newton'one would introduce propagation errors).
3. Results
The reckoning has yielded the following double S constant series in the two intervals:
Taking the geometric average gives for S = 2.47E-19 m-1. Taking also the semi-dispersion for error estimate about each binary system, the relative error may be expressed asgetting for =±1.16E-19 m-1.
With these value we are able now to draw the energy loss with respect to the orbital phase.
We report the instantaneous emitted power for a given direction at polar angle and azimuth corrected of the time averaging factor in the direct calculus versus first approximation (Figure 2) for the Vortex Model for the B1913+16 system. As a collation, we report also the Einstein's one as known proportional to the products of third derivatives of the quadrupole moment [9]:
(7) |
It is worth mentioning that a single point in the direct calculus takes about three hours of processing. We note in this direct calculus a moving to the right of the secondary negative peak but we have had to cope with a larger reckoning approximation, although more realistic, with respect to an ad hoc model such as that of Eq.(3). In the direct calculus the peaks are obtained correspondently to the values of Table 1.
Our strict criterion in calculus of precision, as recalled in [1], is given by the effective symmetry of the deformation tensor (Figure3(a)) and the skew-symmetry of C Riemann tensor (Figure 3(b)) along the orbit, i.e. to the semi-axis extent concerning the two tensors that appear in the evaluation of the integral in Eq.(1), as included in the constraint on the variables of the Dirac delta of the EM tensor of the two point masses of the binary system. The accuracy is shown bettered with respect to our model Eq.(3).
We then report the B1913+16 orbital period decay advancing at periastron time in Figure 4. Expressing it{5} [13] as at the modified julian epoch tn (MJD), we may collate directly the Vortex Model with TGR tendency for in the matching of the experimental fitted curve. We remember The TGR prediction is to within the celebrated 0.2 percent [14].
By inserting the mean energy as from the virial theorem in the 3th Kepler’s law and differentiating, the variation of the orbital period comes out as , ; By replacing the values, one gets = -2.4042(2)E-12 for TGR versus -2.417E-12 for the Vortex Model against = -2.4056(51)E-12 [14]. We so have 0.13±0.21% agreement for TGR against 0.5% for Vortex. However better can be found for the other two systems.
For what concerns the other star two systems B1534 and the double pulsar J0737 we too report as before the differential energy loss in the Vortex direct calculus and 1st approximation paralleled with the respective Einstein’s prevision. The differential functions prepared to be averaged by the usual factor along the orbit appear to disagree of at least two order of magnitude as to B1534 and J0737 as well as for the other two systems as seen below. However what physically matters is the integral over the solid angle to get < dE/dt > and this points out all the more we are dealing with two obviously different models. Vortex Theory is seen to be more reliable in predicting period shift as for these two systems. Data values other than available are extrapolated assuming as before the rotator model of the accumulated phase in the relative motion .
About the relativistic system J0737 we calculated = -1.2523E-12 for Vortex vs = -1.2592E-12 for Einstein’s theory against the intrinsic experimental value = -1.252(12)E-12 [16]. Hence the most precise test ever obtained for Vortex’s Model.
We may again observe in Figure6(b) the overcoming of Vortex theory to TGR. The outcome = -0.173E-12 seems indeed to resolve the small discrepancy puzzle between TGR and observed value of period variation, proper of this binary system. In comparison with =-0.174(11)E-12 we have in fact =-0.192E-12 [17].
At last we have reviewed the B2127+C11 dwarf-NS system very similar to B1913+16 for which = -3.98E-12 suits to = (-3.95±0.13)E-12 as well as = -3.98E-12. However uncertainties prevent us to usefully collect directly the models. We report the angular distribution of the radiated energy against Einstein’s one:
4. Discussion and Conclusions
In spite of all glitches of numerical quadratures, we deem the new value of the constant S =(2.5±1.2)E-19 m-1 more meaningful than our model at first approximation (S|1 =3E-14 m-1). As has been said in the reference article [6], it can be related to the mean free path of graviton through a Compton-like wavelenght = S -1 =h/mc meant to attribute a value to the "inertial mass" of graviton, in this case amounting to (5.5±2.6)E-61 Kg, not far from former bounds [7, 8, 11, 12]. Applying this Vortex Model to keplerian problems at galactic distances as for galaxy rotational curves or solar system issues is a very entangled challenge to be faced with intensive amount of skill analysis. Further verifications and bearing out for consistency of this result other than bettering of statistics should also be provided from other binary NS-NS or dwarf/NS systems such as for instance J1756-2251 and J1518+4904. An indirect confirmation should come also from the energy radiated from the problem of Black Hole in-falling into another [15], excepting over dimensioned RAM required.
Appendix A: Tidal deformations in Vortex Theory
Because instructive how Vortex Model works, we want to get the strain in a given direction as expressed by the formula (17) of [1] caused by the tidal force of a gravitational wave impinging on a certain region free-falling in the space-time. We will follow the notations therein by leaning on the derivation from the quoted books of Ferrarese and Sokolnikoff [4, 10]. Let us consider the evolution with respect to the proper time of a continuum of identifiable material particles from the reference or initial configuration C0 to the actual one , which we refer to the reference frames X and Y respectively, (they are different in general, see the two terns of Figure8), whose a typical point P will have coordinates and in passing from the undeformed state to the deformed of final one of the continuum. Thus, we may take the deformation of C0 in to be a continuous and one-to-one function so that the transformation of the curvilinear coordinates of the material particle P is singled-valued:
(A.1) |
depending on whether we use the initial or final coordinates as independent variables we will have the so-called eulerian or lagrangian motion standpoint.
Therefore, once fixed the initial element of coordinate the evolution of is determined uniquely by knowing the proper time (molecular point of view). Without loss of meaning O could coincide with P0 whose coordinates for convenience we choose hereafter to be the orthogonal Cartesian ones with respect to the tern (this is equivalent to refer them to a locally geodetic frame) whence the notation .
Moreover, since we are interested in detectable effects we can also omit the temporal components and in our reasoning because decoupled from the space coordinates as discussed in the getting the expression for the deformation velocity Kσ in [1, 2], and because the Minkowski tensor is of signature (+,-,-,-) it suffices to change sign of according to Eq.(2):
(2b) |
to work with space distance positive in the next formulas. So, beware here is positive definite .
Therefore, for the line element in we have with , and unit vector of dP0 . On the other hand, by definition, so that
(A.2) |
Similarly, we can write for the line element in P :dl2 = dP · dP, where dP, referred to the orthogonal tern in P0 reads as from the relation . Hence
(A.3) |
We define the physical component of the strain in a given direction as the elongation for unit length in that direction. Taking in account (A.2) e (A.3) we have:
or through Eq.(2b)
(A.4) |
i.e. Eq.(17) of [1]. This formula can help us to infer the formula for angle deformation between coordinate axes and to interpret the off-diagonal components of the strain tensor (or tensor of tidal forces if we think of this model); in fact let consider line elements being parallel to axes before deformation, say: after the deformation, by following the reasoning leading to (A.3), we shall have for their scalar product:
(A.5) |
and because of Eq.(2b) sign changed we get:
(A.6) |
Finally, on using (A.4) that is to say:
(A.7) |
and inserting it into Eq.(A.6) via (A.5) we are able to express
(A.8) |
so that if we term the amount to complete a right angle after deformation of coordinate axes and by taking the sinus function of we obtain the (A.8) again, and therefore:
(A.9) |
namely the decreasing of angle between axes direction is related to so called the shearing strain or off-diagonal element of strain tensor. To conclude we show in Figure 9 the effect of the presence of shearing strains on an unit vector according to Eq.(A.4) disregarding angular axes deformation: it is plain how cylindrical symmetry gets broken. Of course, polarized produces in-phase elongations in the direction of coordinate axes as shown in [1] according (A.7).
Appendix B: Action of the model and a cartesian derivation of main equations
We can infer the main equations of the Vortex Model from a variational action principle in a form remembering that for the coupling of charged matter current interacting with an electromagnetic field:
whereupon the strength tensor in which play the role of the field potential as in electromagnetism and the second term in brackets the matter source term. From that the EM tensor in free space reads:
We already know the variation leads to wave equation in presence of matter as seen in the reference paper [6]. By the way, let us show that the classical deformation speed corresponds to the deformation tensor of mechanics. Defined the local basis of vectors at each point of the cronotope , differentiating with respect to , we get for definition of Christoffel symbols : . From that, it turns out that the gradient of the 4-velocity leading to the classical expression:
We know that although the vortex and the deformation tensor are both independent, they are differentially bound through the Vorticity formula of continua that we now derive in 3 dimension cartesian coordinates as due to Zhen Li (Tsinghua Univ.), addressing the reader to the more thorough approach of Ferrarese and Stazi [4]. Standing the definition of vorticity:
the Jacobi relation turns out:
Given that, then we have:
and remembering that :
hence the vorticity formula for its gradient:
which holds in V4 general coordinates as well. We make notice that in the definition of vorticity and deformation tensors [1, 6], both sides of previous equation are scaled by the speed of light c in order to introduce the constant S of dimensions of [lenght]-1, whence the ansatz: , , so we have the main equation:
Acknowledgement
I'm grateful to CINECA supercomputing section for supplying its 192Gb RAM HP-DL980 hardware.
Notes
1. The reader should notice the change of notation with respect to [6] to banish any misleading of the skew symmetric Riemann tensor obtained as (respectively in [6] Eq. 3.3; 4.12) with that of Einstein's equations Here is meant derived by contraction of the first order estimate of the Riemann tensor as stressed in [1].
2. We recall some hallmarks of the B1913+16 binary system: m1:1.4414(2) Msun; m2:1.3867(2) Msun ; e: 0.617; orbital period T : 7.72 hours; major semi-axis a: 1.950E9 m; distance: 5kpc (16300 LY); dT/dt: -2.4056(51)E-12 (see also ATNF catalogue)
3. We know from the sampling theorem that in order to rebuild a periodic function, there is need for a sampling frequency at least double of that maximum of its spectrum (Nyquist frequency). Nevertheless taking yet three step per interval entailing a course estimate of the numerical integration need about 500Gbytes, while four steps per interval would involve too much bigger files to manage for , of order of Terabytes.
4. as known from the free divergence of EM tensor. Indeed, we should strictly have written with L(g) the restriction of the metric over a surface enclosing the binary system, but we take it instead as the usual jacobian in spherical coordinates in account of the sidereal distance from it.
5. The formula is got by expressing the accumulated phase as function of time , expanding the frequency v to 2nd order and drawing out the difference between the times of passage at periastron and the n-th proper period time n/v0
References
[1] | Antonelli,S., “Outstanding Outcomes from a Recent Theory of Gravity”, International Journal of Physics, Sciepub, 2,6,2014, pp.267-276, 2014. [Online] Available: https://www.sciepub.com/portal/downloads?doi=10.12691/ijp-2-6-10&filename=ijp-2-6-10.pdf. | ||
In article | |||
[2] | Landau, L., Lifchitz, E., The Classical Theory of Fields, vol.2, sec.82, Butterworth-Heinemann, 1975. | ||
In article | |||
[3] | Loinger,A., “Gravitational waves are fictitious entities II”, https://xxx.lanl.gov/abs/astro-ph/9904207 , 1999. | ||
In article | |||
[4] | Ferrarese,G., Stazi, L., Lezioni di Meccanica Razionale , vol 2 chap.VIII, §1.10, Pitagora Eds, Bologna, 1989 . | ||
In article | |||
[5] | Ferrarese, G., Bini, G.., Introduction to relativistic continuum mechanics, Springer, Heidelberg, 2008. | ||
In article | View Article | ||
[6] | Tailherer, M., “A Critical Reading on the Theory of Gravitational Wave Propagation”, Scientific Journals, Journal of Physical & Natural Sciences 1, 1 , 2007. | ||
In article | |||
[7] | Gershtein, S.S., Logunov, A.A.,Mestvirishvili, M.A., “The upper limit on the graviton mass”, 1997. [Online]. Available: arXiv:hep-th/9711147. | ||
In article | |||
[8] | Visser, M., “Mass for the Graviton” ,1998, arXiv:gr-qc/9705051. | ||
In article | |||
[9] | Ohanian, H.C., Ruffini, R., Gravitation and Spacetime, Chap.5 §5.2, W.W. Norton & Company, 1994. | ||
In article | |||
[10] | Sokolnikoff, I.S., Tensor Analysis, Chap. 6, §105, Wiley, New York, 1951. | ||
In article | |||
[11] | Finn, L.S., Sutton, P.J., “Bounding the mass of the graviton using binary pulsar observations”, Phys. Rev. D65, 044022, 2002. [Online]. Available: arXiv:gr-qc/0109049. | ||
In article | |||
[12] | Will, C.,W., “Bounding the mass of the graviton using gravitational-wave observations of inspiralling compact binaries. Phys. Rev. D57, 2061. (1998) [Online]. Available: arXiv:gr-qc/9709011. | ||
In article | |||
[13] | Maggiore, M., Gravitational waves, vol. 1, §6.2.3, OUP Oxford, 2007. | ||
In article | |||
[14] | Weisberg, J.M.,Taylor, J.H., “Relativistic Binary Pulsar B1913+16: Thirty Years of Observations and Analysis”,2004 . [Online] . Available: arXiv:astro-ph/0407149 v1 | ||
In article | |||
[15] | Antonelli. S., “Importance of the Gauss' 2nd differential form in Relativity”, viXra e-Print archive, 2013, [Online]. Available: https://www.vixra.org/abs/1310.0190. | ||
In article | |||
[16] | Kramer, M., et al.,“Tests of general relativity from timing the double pulsar”, 2006, [Online]. Available: arXiv:0609417 [astro-ph]. | ||
In article | |||
[17] | Wolszczan, A., “A nearby 37.9-ms radio pulsar in a relativistic binary system”, Nature, 350, 688 ,1991. | ||
In article | View Article | ||
[18] | Blanchet,L., “Post-Newtonian theory and the two-body problem”, 2010. [Online]. Available: arXiv.org/gr-qc/0907.3596 | ||
In article | |||
[19] | Jacoby, B. A, et al. “Measurement of orbital decay in the double neutron star binary B2127+11C”, The Astroph. J., 644, L113-L116, 2006. [Online]. Available: arXiv:astro-ph/0605375. | ||
In article | |||