1. Let Us Now Look at That Famous “Simple Solution” Template for Massive Gravity
In  pages 212-213, we have that there is a Minkowski simple model for massive gravity, leading to
Whereas if we us , and  and really look at , we can add in the following as to the scale factor used in Equation (1) which uses at the surface of a presumed non singular start to the expansion of the universe
Leading to  the inflaton.
And what we will use later the “inflaton potential” we write as
Equation (1) and Equation (2) lead us to the interesting restraint of
Our preliminary consideration is to have the time, in Equation (5) commensurate with Planck time, whereas a radius of the “bubble” of initial space-time commensurate with Planck length, using the dimensions given in .
We will see where this becomes important and how it factors into discussion of the strength of a GW signal and possibly polarization.
2. A Review of Nonsingular Cosmologies and the Friedman Equation, and How We Utilize Equation (5)
A given in our work is that within the nucleation of space-time, that time, even as given in Equation (5) simply does not exist, but that we will be able to use the results of Freeze  as far as a nonsingular Friedman equation to come up with
Here, we have that is a space-time density function, whereas is related to the tension of a space-time bubble presumably of the order of a Planck radius. And we are also using what is given in  as far as a Dark Energy model, which we write, for energy density, as given by, if is the Ricci scalar  and we use also the notation of . We have also the discussion given in  and set as spacetime curvature, and so then we have
We cite references   and  as to the work done with regards to quantum gravity needed to formulate if we can have space-time with a non singular start to expansion of the Universe. We do this because the quantum gravity approximation used, as given in   and  meant that time as we think of it, at the onset of the bubble of nucleated space-time does not exist in the usual causal structure we identify with, in later expansion of the Universe.
In terms of the bubble of spacetime before inflation, we submit that time does not really exist and that then we will be considering a rewrite of the above as having, effectively . And the term . And if we apply Equation (6), we have within the bubble of spacetime.
Then we have the following, namely
Then we get
This value for at just about the surface of a bubble of initial space-time nucleation, would be if were small but not zero
We would then be tending toward, if we used
If we have, then say small value
Then we would have, say for small is small, then we have, that we could even see the following, namely
In either case, we will be doing our calculations to determine what this has to say as to the frequency of a signal from this event, as well as the strength of GW, and then also the possible polarization states. This would have to be contrasted with Equation (5), as to what that says about bounding values for the input into Equation (15) below.
Let us now refer to a simple but specific model given by Dr. Corda, as to the way one can implement an amplitude for GW  which the author views as commensurate with regards to amplitude for this initial system of GW generation at the start of the expansion of the Universe.
From  we have that there are minimal extensions of General Relativity for the existence of GW generation which may be thoroughly explored which will enable us to have GW early universe production. Before we get to that we will briefly review what is known about the GW generation of the classical kind which is a precursor to understanding the problems with GW amplitude calculations in the beginning of Space-time evolution.
3. Brief Review of the Basics of GW Amplitude Calculations in the Traditional GW Theories
To do this, we look at  which is Maggiore’s first volume masterpiece of GW data analysis and collation of data. In doing so, on page 211, we have that if a GW were signified to run in the Z direction, then in lowest order, assuming that we are using, if time is retarded time, then the two main polarizations are describable as
Here, we have that, usually, .
In doing so, we can write up what we have here with respect to the Einstein Energy tensor as
Here by  we can write
And the Ricci Scalar as
Whereas we can then proceed in making sense of what Equation (16), Equation (17) and Equation (18) are saying in terms of Massive Gravity, and relic GW generation.
The action will start in the interpretation of Ricci tensor, time domain, and R as a Ricci scalar term.
4. Interpreting R00 Ricci Tensor, Time Domain, and R
 gives us a different view of the Equation (16) to Equation (18) terms which we can summarize, below. Afterwards we will go to that procedure brought up by Corda in  in terms of relic conditions. Here, the Ricci scalar is given directly as
Whereas the term, given in Equation (17) gets a new paintjob which we can have as by Kolb and Turner  for a Roberson-Walker cosmology, if we stick initially with Equation (17) would have the following equality
Whereas we would re write Equation (19) as the following
In our case we would consider the Equation (17) to be dominant in our problem, and therefore choose Equation (16) to read as
i.e. of all things a negative energy density, which is weird, at the surface of the bubble of space-time, whereas if we use Equation (21) directly we would have, instead
Whereas then we would have the following to contend with, as far as
Whereas we could write, the other magnitude as
These would have to be configured to be at the surface of the space-time bubble with the Time derivative, in this case, according to Maggiore to be retarded time, and also using Equation (2) for the scale factor. Whereas we could put in different candidates as to the g00 term.
5. Comparing These Values of with the Density Values of Equation (10) for Density, as Put in Would Lead to Equation (24) and Equation (25) Re Written in Terms of the Bubble of Space-Time with a Base of Equation (13) Would Remove the Potential of Negative Values
Here is the problem if that is done. I.e. the answer would be dependent upon a reading, experimentally upon the tension of the bubble of Space-time presumed. Difficult but not impossible. But necessitating a huge amount of work. would need experimental confirmation, whereas we would be using retarded time in the time differentiation.
This though lacks the third term from an additional polarization term which was used by Corda in  which will be referred to later. Before doing that, look at the additional development we could do, with Equation (10) for a GW density term. We can though add in the cross polarization term before going to the Corda treatment in . It is his Equation (48). Which is our Equation (28)
6. Linking to the Corda Approach in  for Relic GW Conditions, i.e. Adding Equations Together
The Corda approach, in Reference  is in its Equation (48) a representation we will refer to as setting, if we used the cosmological parameter referenced in Equation.
Here we make the following identification, namely
The missing term is in . We will make some comments as to this additional term.
7. Adding in and Coming up with a Primordial GW Amplitude Generation
To do this we will consider having a look at the Equation (45) of reference  which we will write as
We claim that this would be commensurate with the Hubble expansion H0 defined by setting it to today’s space time crucial density
Whereas earlier we would have
Our approximation in our paper is to use the following substitution, namely make the following substitution
Then to use an integrated version of this in Equation (30), Equation (31), and Equation (32) and then come up with a linkage to analyzing
Re write this as follows set
Then the relic GW density goes as
This would lead to a third relic condition GW contribution we would call
8. What about Tensor-Mode, Vector-Mode and Scalar-Mode Gravitons as Spin Zero and Spin 1 “Gravitons”?
We will review our work further but it will be in strict fidelity with what Wen Hao and Fangyu Li did, which is in ARXIV and was just recently published in European Journal C, i.e. see this .
Reference  is unusually complete. What the author is arguing for is that since the details of the scalar and vector mode interactions are highly detector specific, that adding in these 2 additional polarization states, akin to the five polarization states of massive gravitons will await fine tuning the details of the material in this document, as to the first three polarizations, as well as getting more details as to the instrumentation of the Li Baker detector. This will be followed up upon, pending utilizing the details as to the following document  .
My guess is that what is done in  is for higher than four dimensional, brane physics models, whereas our starting point is for 3 + 1 dimensional models, but that the differences are due strictly to higher dimensions employed.
Not that this is alluded to, as far as future work, we will go to the issue of relic GW frequencies.
9. Order of Magnitude Estimates for Frequency, Say of Relic GW
We use the most direct route toward doing it, and say if we used
For sufficiently large value we hope we have a negative value
If instead we have a small value
These expressions should be compared to Equation (5) for which we can write, say 
We think that this means that at the boundary of a space time bubble, that this would force us to have an enormous value for . We would also, if the initial energy were just at the boundary of the bubble of space time, have the odd situation for which we would have the following, namely at the surface of the bubble, so if one used Planck units with we could have the following weird situation
To put it mildly, we would need a lot more experimental data sets! All this leads to considerable interplay between the equations given in Equation (38) to Equation (42). We will discuss this in our conclusions. Keep in mind that the interplay with all of this, we will try to right after the Bubble of initially nucleated Space-time to have a linkage to simpler Quantum mechanics involved as seen in . In addition we wish to in doing so to have fidelity with the insights of Christian Corda, as seen in  as to mixed Scalar-Tensor modes and models of Gravity.
In our write up of Equation (38) to Equation (42) we are assuming for all general purposes that the scale factor would have a minimum value. Keep in mind that the surface of the bubble would be, in fidelity with reference  be obeying having the scale factor at a minimum value and that in line with inflationary e folds of 60  , or 1.14 times 1026. If this is kept in mind, and say that we have a frequency range of about 1037 Hertz, as a result of the above, we would then have say 1011 Hertz for Earth signals taken for GW detection.
10. What Is the Strength of a Signal for Our Model?
GW signals have a simple strength of GW moniker   The simplest idea is to look at the behavior of GW from a massive black hole, and in doing so, the strain is simply M the mass of a black hole, r = distance of source to Earth from Black hole, with v as given in Equation (43) being the speed of masses, just short of the speed of light. Then
If the original universe were acting say as a black hole just before the onset of the black hole exploding via the violence of inflation, if we wish to have an observed GW signal, and have say v/c about 1, in terms of the early universe, with r being the distance of the expansion in 13.6 billion years, we could have a configuration for which if we were near the initial start point of the Universe, that Equation (43) would be of the order of 1, whereas, the e fold expansion of 1026 from an initial start point would then have say a signal strength of the order of 10−25 to 10−26, and if so a strain value of 10−26 would be commensurate with GW signal detectors being designed as of Earth orbit where we are right now. Verifying though the initial configuration of Equation (43) as approaching unity, though, would require serious investigative work and would entail perhaps spending attention to the fractal geometry given in  .
The entire reference in  is essentially a referencing of the need for fractal dimensionality, and where it fits in, as far as the geometry of space-time. On page 94 of  the Hausdorff-Besicovitch dimension, which is fractal is done, in gory detail. Why this is important? Note that in  we made reference to , which if d (dimension) is fractal, means that the temperature for the universe, so assumed, has many surprises if we wish to connect Equation (44) with a scaling of i.e. there is then likely not a simply linear relationship between temperature, and frequency, hence we have to be more careful.
Now, make use of  again and this relationship below though will stand the test of time, but it is NOT dependent upon frequency! Here is a count of gravitons which may commence from our earlier generation of early universe GW. From 
In other words, due to fractal geometry, while Equation (22) and Equation (23) may hold, as stated in , we may have some observational surprises as far as frequency spectrum data sets as far as GW from relic initial GW generating events in the early universe.
In other wise, while the scaling of 10−26 downward in terms of general frequency from a nonsingular universe starting point is likely correct, the details of the frequency spectrum in the modern era from Primordial beginnings may be a bit different from what we think.
11. Is There Another Way to Form a Quantum Wavefunction of the Universe Rather than Wheeler De Witt?
In  there is reference to a solution for which we still have
Whereas we have, and
Whereas  sets as an energy density, for a whole slew of grab bag assorted “topics”, and
Here, is a so called Order parameter, and after a small scale factor solution we would get a linear combination of Bessel and Hankel functions for the wavefunction of the universe, which for small scale factor allowed to go to zero we have, if and are constants, and the scale factor can get arbitrarily small. The Equation (49) below is specified by Reference .
Keep in mind that when we do this we are using the time dependence given by Equation (37), Equation (38) and Equation (39) whereas we will be specifying a frequency dependence given by Equation (40).
So, how could we write in a minimum time step whereas we also have a quantum function?
I will make a suggestion as such, whereas, it may be not accepted. But this is my candidate. From Gasiorowicz , for a wave function just outside the minimum scale factor, I would try using Equation (39) and Equation (49) to come up with a wavefunction at the surface of the space-time bubble which would have much the same information as given in Equation (27) above. To do this though we will refer to a result by Shestakova , and , and then re interpret it via some physics from Stephen Gasiorowicz  and compare that directly with Equation (46), and Equation (47).
12. Simple QM Wavefunction at the Surface of the Space-Time Non Singular Start Point Matched against Equation (27)
We begin with  and  by T. P. Shestakova where we have on its page 9, if we use order parameter p = 1, the following S.E. and also a solution, which is dependent upon a scale factor with a dimensional factor put in as a fudge factor
If so then, we will be looking at, using Arfken, 4th edition, page 667  so that the solution to Equation (50), given by,  and  has the form
Note that the CRC handbook,  gives us an expansion so we can do the following
Here if is very small whereas we take we would have
Leading to a comparison between two forms of the Wheeler De Witt equation give the same information, i.e. questioning if the following are equivalent
If so, and the answer is yes then we will then have to address what a minimum bubble of space-time ascertains as to future developments of space-time, and for this we will be looking at  for inspiration.
13. How to Link Equation (54) to a Schrodinger Type Equation for Planck Time Physics
The relationship can be ascertained as to stating that we can use right at the surface of the space-time bubble
Keep in mind that this is dependent upon a linkage to Equation (34), which has no explicit time dependence whereas Equation (55) has a lot of time dependence.
14. And What about Polarization States, Especially If We Have Massive Gravitons?
See . Our supposition is that as given in  that the massive gravity signature we would seek would be commensurate with “pure longitudinal and transverse breathing polarizations in the massive Horndeski theory and f(R) gravity”. The details of such are given in their  document as given just above their Equation (34) whereas we would have us calculate the electric component R (tjtk) of the Riemann tensor, as given in their document on page 10 of their arXIV document.
The details of their Equation (34) as in Reference  would in the end have to be made fiducially relevant to the experimental platform as selected by an experimental gravity research team, and in itself would require a massive summary of known and sought after we have already outlined a procedure tying in what we did with  as given by Dr. Corda. We intend to find Graviscalar and Gravivecor contributions, in line with work presented which is in fidelity with quantum gravity states as summarized already in this document.
15. Summing Up: And a Path Forward
1) I am expecting a strain of Gravitational wave strength of h ~ 10−25 to 10−26.
2) A frequency range for GW detected in Earth orbit of about 1010 Hertz, which could be commensurate with initial GW of 1036 to 1037 Hertz for reasons which are in this manuscript. There would be due to e fold values of 60 a dramatic drop in frequencies.
3) The polarization states to watch, indicating if possible, would be for massive gravitons.
4) Should there be a linkage between an initial worm hole style start to the formation of the observable universe, as a compliment to the big bank/inflation metrics, the author strongly recommends a review of the article given by Bruno J. Barros and Francisco S. N. Lobo as to relevant Wormhole geometry , and the question of negative energy states in a wormhole threat. If as may be the case in nonsingular starts to the expansion of the universe, negative energy states make an appearance, this may be a clue as to what is called f(R) geometries which in turn may have relevance to some of the polarization states in .
The wormhole analogy may be crucial as to looking at the linkage between super massive black holes and recycling of the universe, which may be investigated later.
5) Also the author recommends a doubling down on the physics of this diagram in the last  article.
See Figure 1 below which is in . We argue that the initial physics, may help us understand more of the self replication phenomena so outlined.
Figure 1. Which is from reference .
We also should understand if the physics of a linkage between Equation (33) and Equation (34) play a role as well.
6) Finally, and not least is to see if our feed into Equation (33) to Equation (34) and linkages to Figure 1, have in their genesis, the methodologies discussed in , i.e. as brought up in Phase transitions and the renormalization group as brought up by Jean-Zinn-Justin, which has a very interesting discussion of the Nancy Kerrigan problem, i.e. as of page 133 of  which states that DE did NOT interfere with, or put constraints upon structure formation. We submit that Figure 1 is a case of the DE density, initially being replaced as we did earlier, in the manuscript with cosmological energy density, using the equations we copied from Gong, et al. . This requires further investigation, and may if resolved have a tie into Figure 1 physics, of .
7) Doing all the above, would be hard, but would be the cornerstone of introducing the concept of using the construction of a non-Archimedean field to physics , which would be a way also to take into account what may be fractal geometry incorporated into space-time physics. In mathematics, a non-Archimedean ordered field is an ordered field that does not satisfy the Archimedean property. Examples are the Levi-Civita field, the hyperreal numbers, the surreal numbers, and the field of rational functions with real coefficients with a suitable order.
This would be a way to understand more of Equation (22), if the dimensional coefficient, d, is rep-resenting fractal geometry.
16. And Now for Some Very Specific Goals for the Future
We have referenced several different polarization methods for massive gravitons. Our idea is to form a bridge between our extension between  as a minimal extension of relativity and the methods brought up by Dr. Hao Wen, Dr. Fangyu Li, and others recently as published in European Physics Journal C, given by  with an eye toward coming up with specific counter parts to the work of Gong, for massive gravitons . The work in  is outstanding but in a sense very general, as we will try to look at what polarization states survive the onset of quantum gravity conditions at our nonsingular start to the universe.
In addition, the author at the 4th Zeldovich conference, November 11, at 14:00 CET time made reference to the enduring mystery of error bars in the CMBR measurement as discussed by Abhay Akshenkar  of how his loop quantum gravity may form a solution to this problem. My entire nonsingular start to the evolution of the universe was initiated as a way to attack this problem addressed by Abhay and I will endeavor to do more on that in a future publication.
Keep in mind that there are 61 or so e folds from the start to the finish of inflation, allowing the universe to expand to the size of a large grapefruit, and that the total amount of expansion in 13.4 billion years of evolution is about 141 e folds.
Finally Karen Freeze in Zeldovich 4  spoke on the physics of massive Dark Stars which she said may contribute to supermassive black hole formation in the center of elliptical galaxies. This insight will be part of an extension of our work into gravitational astronomy and will be referenced as well in a future publication.
This work is supported in part by National Nature Science Foundation of China grant No. 11375279.
Personal big thank you to Dr. Christian Corda for his outstanding questions which are highlighted in this manuscript.
 Tolley, A.J. (2015) Cosmological Applications of Massive Gravity. In: Papantonopoulos, E., Ed., Modifications of Einstein’s Theory of Gravity at Large Distances. Lecture Notes in Physics, Springer, Cham, vol. 892, 203-224.
 Amico, G.D., de Rham, C., Dubosvsky, S., Gabadadze, G., Pirtskhalava, D. and Tolley, A.J. (2011) Massive Cosmologies. Physical Review D, 84, 124046.
 Beckwith, A. (2020) How to Obtain a Mass of a Graviton, and Does This Methodology Lead to Voids? Journal of High Energy Physics, Gravitation and Cosmology, 6, 416-439.
 Padmanabhan, T. (2006) An Invitation to Astrophysics. World Press Scientific, World Scientific Series in Astronomy and Astrophysics: Volume 8. World Scientific Publishing Co. Pte Ltd., Singapore.
 Freeze, K., Brown, M. and Kinney, W. (2012) The Phantom Bounce: A New Proposal for an Oscillating Cosmology. In: Mersini-Houghton, L. and Vaas, R. Eds., The Arrows of Time. Fundamental Theories of Physics, Springer, Berlin, Heidelberg, vol. 172, 149-156.
 Corda, C. (2008) Primordial Production of Massive Relic Gravitational Waves from a Weak Modification of General Relativity. Astroparticle Physics, 30, 209-215.
 Li, F.-Y., Wen, H., Fang, Z.-Y., Li, D. and Zhang, T.-J. (2018) Electromagnetic Counterparts of High-Frequency Gravitational Waves Having Additional Polarization States: Distinguishing and Probing Tensor-Mode, Vector-Mode and Scalar-Mode Gravitons.
 Li, F.Y., Baker Jr., R.M.L., Fang, Z.Y., Stepheson, G.V. and Chen, Z.Y. (2008) Perturbative Photon Fluxes Generated by High-Frequency Gravitational Waves and Their Physical Effects. The European Physical Journal C, 56, 407-423.
 Li, F.Y., et al. (2009) Signal Photon Flux and Background Noise in a Coupling Electromagnetic Detecting System for High Frequency Gravitational Waves. Physical Review D, 86, 064013.
 Corda, C. (2009) Interferometric Detection of Gravitational Waves: The Definitive Test for General Relativity. International Journal of Modern Physics D, 18, 2275-2282.
 Vieira, H.S., Bezerra, V.B., Muniz, C.R., Cunha, M.S. and Christiansen, H.R. (2020) Class of Solutions of the Wheeler-DeWitt Equation with Ordering Parameter. Physics Letters B, 809, 135712.
 Shestakova, T.P. (2019) On the Meaning of the Wave Function of the Universe. International Journal of Modern Physics D, 28, 1941009.
 Zinn-Justin, J. (2003) Phase Transitions and Renormalization Group: From theory to Numbers. In: Duplantier, B. and Rivaseau, V., Eds., Vacuum Energy, Renormalization: Poincare Seminar 2002, Springer, Basel AG, Berlin, 213-239.
 Beckwith, A. (2020) Using “Enhanced Quantization” to Bound the Cosmological Constant, and Quantum Number n for Production of 100 Relic Mini Black Holes and Concluding with a Proof of a Minimum Time Step. Zeldovich 4. September 11 Section on Dark Matter and Dark Energy.
 Freese, K., Rindler-Daller, T., Spolyar, D. and Valluri, M. (2016) Dark Stars: A Review. Reports on Progress in Physics, 79, 066902.