1. First Let Us View the Penrose Suggestion as to Black Holes from a Prior Universe
In order to see this, consider a suggestion as to black holes, being the template for a start to the present universe given by , and also  which has the Penrose suggestion of an imprint of a prior Universe black hole having an effect upon the CMBR spectrum. The CMBR spectrum is a real datum, but the worth of getting this information would be in terms of having what was said in  as to the “ghost” of prior universe black hole radiation. This idea is presented in our Figure 1, which is copied from page 2 of reference .
Figure 1. A conformal diagram representing the effect of a highly energetic event occurring at the space-time point H. In CCC, H is taken to be a Hawking point, where virtually the entire Hawking radiation of a previous-aeon supermassive black hole is concentrated at H by the conformal compression of the hole’s radiating future.
i.e. the image of a supermassive black hole is then superimposed upon the CMBR and then the question is, does this have any other effects than a signature on the CMBR?
See a summary of this in  which is a way of saying that the CMBR can have traces, via CCC recycling of Black holes.
The CMBR can be affected by Gravitational radiation, as seen in , with more details in . We argue in doing all this that we have in actuality two choices, the first one, on if we have a singular start to the universe, as given by the Penrose singularity theorem as given in , and the next one if we have a nonsingular modification of the Friedman equation as given in .
Having a non singular scale factor will allow for the 2nd case which we argue is in tandem with the treatment of the ghosts of the supermassive black holes is a way to initiate putting the ghosts of prior supermassive black holes into the present CMBR, which has to be reconciled with the generalized facts given by Zeldovich, in , and congruent with  .
2. The Basic Facts of the Penrose Singularity Theorem and Why We Do Not Assume a Point Scale Factor, as an Initial Starting Point
 has this foundational structure which is relevant to out inquiry
8.2. Quantum effects a very important line of research arises from the tension between the singularity theorems and the (yet unfound) theory of quantum gravity. It is widely accepted that the existence of classical singularities signal a breakdown of the classical theory at extreme conditions, which is precisely when gravitational quantum effects will become relevant. Thus, there is a need to clarify if the singularity theorems, or part of them, can survive when entering into a quantum regime, or if they then simply vanish altogether. For a general discussion, see . A first step towards the analysis of singularity theorems in this respect is the weakening of the “energy conditions”—also relevant in the classical Singularity theorems 38 regime, that is to say, finding an appropriate version of the curvature condition in the theorems. Early results in this direction include the theorems based on averaged energy conditions  as discussed in subsection 5.1, which were used to deal with the quantum violations of the energy conditions in , later improved in . A larger discussion can be found in  (and references therein) and has been recently newly considered in , where an analysis of Raychadhuri-like equations is performed proving that it is viable to have energy-momentum tensors which fail to satisfy even averaged energy conditions as long as an appropriate version of them with an exponential damping factor are in place. This leads to a proof of a version of the Penrose singularity theorem allowing for global violations of the energy conditions.
End of quote
What we will be analyzing is a definite quantum violation of the energy condition as given by .
Having said that,  has some of the basic facts to consider with a statement of the Euclidian Universe which is very pertinent to our situation.
Unlike the black hole pair creation, one couldn’t say that the de Sitter universe was created out of field energy in a pre-existing space. Instead, it would quite literally be created out of nothing: not just out of the vacuum but out of absolutely nothing at all because there is nothing outside the universe. In the Euclidean regime, the de Sitter universe is just a closed space like the surface of the Earth but with two more dimensions. If the cosmological constant is small compared to the Planck value, the curvature of the Euclidean four sphere should be small. This will mean that the saddle point approximation to the path integral should be good, and that the calculation of the wave function of the universe won’t be affected by our ignorance of what happens in very high curvatures.
End of quote
Our point of divergence from this statement is A. We do NOT assume that the Universe is made out of “nothing” because we assume via a version of CCC cosmology  that there has been regular inputs into this present universe; B. We do assume that we may have a wavefunction of the Universe which does not take into accord high curvature regimes; C. We also though in making the wavefunction of the Universe assume what was given in  so as to avoid having a scale factor involving a singular point, which would make our suppositions unworkable.
3. Stating the Wave Function Treatment of Acceleration from the Initial Phase Transition
In a word we are assuming that there is a marked phase transition between early time pre expansion physics, and then the expansion physics. To do so we are looking at . i.e. what we are doing is akin to the following description which was pitched to be relevant to the electroweak regime. With certain caveats, we hold that it is similar to what we would have at the very start of expansion of the universe. i.e. change the description from the electroweak to the initial nonsingular starting point of the universe in line with  and we are set to go for our cosmological evolution. See .
Phase transitions are a generic, but not universal, feature of gauge field theories, like the Standard Model, which are based on elementary particle mass generation by spontaneous symmetry-breaking  . When there is a phase transition in a gauge theory it is (except for special parameter choices) first-order, which means that just below the critical temperature, the universe transitions from a metastable quasi-equilibrium state into a stable equilibrium state, through a process of bubble nucleation, growth, and merger    . Such a first-order phase transition in the early universe naturally leads to the production of gravitational waves  . If it took place around the electroweak scale, by which we mean temperatures in the range 100 - 1000 GeV, the gravitational wave signal could lie in the frequency range of the upcoming space-based gravitational wave detector LISA (Laser Interferometer Space Antenna) . The approval of the mission, and the detection of gravitational waves , has generated enormous interest in phase transitions in the early universe.
End of quote
The evolution of energy section, in , about pp. 34-39 will be obeyed once we commence once we get through the initial phase transitions. In the mean time, once we get to the initial phase transitions,  we have that  describes well the initial gravitational wave generation well.
We are concerned though about if there is not just CMBR support for the Penrose conjecture    and  but also GW signatures from Pre Big Bang black holes disturbing the COMBR leading to gravitational signatures. which outline this addendum to classical GR.
 Hee-Boon Low and Chi Xiong, outline on page 60, Equation (13) an additional term, Vinens equation for a vortex tangle to be added to the usual cosmological evolution equations, which may account for what was brought up in , namely a reconciliation of the Higgs hypothesis given in pp35 to 44, but that would be a bridge for the creation of matter-energy after a first order phase transition. What we will be concerned with next will be the bridge from CCP inputs into the generation of the present cosmos, to our present universe. From between a CCC cosmology construction, with an additional modification put in, to the physics of the question we put in initially which will be of with regards to signatures of Black Holes, from a prior universe to our present cosmos. This will be the reason to use   .
See    as to what we are trying to ascertain exists in terms of GW signatures for black holes in the CMBR   whereas referencing the superfluid universe, as a candidate, after the phase transition  The phase transition details, plus preserving the inputs from  will be the conclusion of this manuscript.
4. Emergent Tunneling Wave Space-Time Equation Used to Compare the Start of Expansion
We use the construction from  as to, if the initial “potential” is very large, how to isolate the form of the wavefunction, especially if , even if a is the initial value, i.e. very, very small, even if , and then by page 269 of  go to the following formulation. Namely that we look at
Our entries into all the above will be the subject of the next several sessions of our document and we will endeavor to explain how this all fits into the early universe modeling we will be working with as far as foundational research into Quantum gravity.
In terms of a contribution from the modified Einstein Equations, we will be considering the time components of the Einstein equation with the stress and strain component in quantum expectation language, i.e. as in Colins, Martin and Squires, written as, due to emergent space-time from pre Planckian to Planckian regimes. i.e. look at the time related transitions, and this is what we get. i.e. the basic equations are related as follow 
Here we will have the following approximations, used. Namely look at
which are a result of the following equations used when we are looking at a tunneling Wave Space-time equation. Then, using this procedure above is using Padmabhan’s inputs into the inflaton potential and also getting initial space-time data into the inflaton itself, as given in .
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 .
Furthermore, we have from  pp. 212-213, and .
And then, a minimum time step we define via a minimum time step of
Note that if the time as defined by Equation (8) is on the order of Planck time, i.e. 10−44 seconds, we have then that .
We then close with a statement, if
If t he wave function is taken via Equation (1) put in, we have that the effective wave function for the “Quantum” part in the RHS of Equation (9) is then
Whereas we will be looking at the following approximation for the right hand side of Equation (10), namely
The above Equation (11) is set equal to zero.
Hence also note we are going to do the following for breaking down the parts of what would otherwise be an ENORMOUS integrand expression to evaluate. Namely the midpoint rule for numerical evaluation of the integral, using . We will justify this procedure due to the “radius” of the initial “space-bubble” is so close to a singularity, and we are talking about emergent spacetime. Hence, this idea as of  is used whereas I would NEVER otherwise even think of it. The centering of the initial problem, about its time component should also be seen in the context that the idea as will be used in further developments is that time, within the space-time bubble itself does NOT exist, and that there is instead a construction of time on the SURFACE of the initial near singularity.
The author is well aware of the Penrose singularity theorem . At the end of this document, there will be a discussion of where I think it (the singularity theorem) holds, and what is missed, which is part of the multiverse discussion which is explicitly referred to and put in our text, in order to achieve having a constant value of cosmological cycle per cycle without referencing the insanity of the Anthropic principle which has been used by Tippler and others to insert a religious dimension into Cosmology, which I view as missing the point entirely as to the evolution of the Universe itself . This is different from an argument given by Tegmark, as of  which gives a surprisingly well motivated argument as to spatial dimensions. I think the argument is well worth investigating. But the Strong Antropic principle is quite another matter entirely, to be avoided.
From there we will use simple approximations to ascertain as to how to analyze Equation (11) and Equation (12) and proceed with the rest of this article mainly using the midpoint rule for numerical approximation.
To do so, we look at the following approximation, namely if
In addition the fact I will refer to a very strong initial domain wall, of a boundary between the early universe, and also state the following.
Pick the following in order to insure the initial wave function as given by Equation (1) remains valid within the context of a bubble of space-time at the start of inflation, namely have:
This restriction is in tandem with having t to the 61st power at the start of expansion of the Planckian space-time regime as reflecting the hyper rapid onset of inflation itself.
5. Utilizing Equation (2) and Equation (3), as to Nonsingular Friedman Equations to Obtain One Value for the Cosmological Constant
To do this, we will bring up several considerations about the Cosmological Constant, within our framework of a non singular start to the expansion of our universe.
A given in our work is that within the nucleation of space-time, that time, even as given in Equation (2) 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
The easiest thing to do would be to make the following identification, namely writing a quadratic equation in terms of density, which we would write as the following form, namely for an unspecified for now the following quadratic equation in terms of a general density function which we would write as
Then we could write the following, namely
We will from here isolate the Cosmological constant, which we write as
The simplification if we have, say with nonzero come up with
We will, instead, in comparing this value of the cosmological constant go to a method of comparing two first integrals, and from there compare that value with Equation (21) and Equation (22).
6. Why Is a Nonsingular Start to Cosmology Important
In Figure 2 we copy what was done by Ashtekar, in Zelsovich 4 as to what was part of anisotropic fits to the E and B polarization, as given.
Questioning by me, the attendee of the Zeldovich 4 conference (I had my own presentation) revealed that the correction Abhay inserted to the left of the “temperatire” spectrum was not decisively impacted in his LQG treatment of the E and B polarizations. I view the item as far from closed, and I also told Zeldovich 4 that a non singular start to the expansion of the universe makes this problem far more tractable. FTR.
We review  and in and in addition look at  I also proceed to Klauder Enhanced quantization and the Cosmological constant which is in the next section.
7. Klauder Enhanced Quantization and the Cosmological Constant
Our methodology is to compare what we are doing with the work done by the author initially without enhanced quantization and to argue that the results are qualitatively a close match, with the Klauder.
We are going to go to page 78 by Klauder   of what he calls on page 78 a restricted Quantum action principle which he writes as: where we write a 1-1 equivalence as in    .
Our assumption is that is a constant, hence we assume then the following approximation, from  .
Our innovation is to then equate and to assume small time step values. Then .
Figure 2. Now looking at what was discussed by Abhay Ashtekar in Zeldovich 4. On September 7, 2020 .
Let us compare for a moment Equation (25) with what is derived in Equation (22). Depending upon what we put in for the density function , what jumps out is that we have, then almost the same positive and negative contributions which will in fact lead to a nonperturbative solution to the Cosmological constant which we claim should, with further research converge to the same values. We will attend to this in the future. We will comment on comparing Equation (22) and Equation (25) in the end of the document with recommendations. See also   and .
Now we have a candidate balance of contributions to Λ at the start of the universe. Assuming DE is Λ how do we account for a possible reduction of the role of DM? From the early universe?
According to Freese, Rindler-DallerSpolyarand Valluri massive BH of a mass 108 solar masses formed early on, in the center of spiral galaxies, forming from DM stars as presented, here . The point is this, that the following pie chart probably had more DM than today, with the Dark Stars perhaps shrinking the amount of DM and boosting the total budget of DR in the cosmos.
We do not exactly know what the DE is, but here is a clue as to its role in the present universe. And DM may have been considerably more prevalent than today, if the SM black holes from the center of galaxies are from giant DM stars then if Equation (22) and Equation (25) are prevalent DE, as we suggest,. Then the giant DM stars as postulated in ;  may be our preferred venue with the formation of elliptical galaxies in setting up conditions for reacceleration of the Universe. This is something to consider and is a research direction which the author will pursue in the future.
Figure 3. Attributed to commonly available descriptions    .
8. Resetting Dark Energy as the Cosmological Constant and Reacceleration of the Universe 1 Billion Years Ago
Here we are using inputs from     , with  in particular] in particular by Peebles and Ratra giving the case that if the mathematics of  and  are taken care of, that indeed the Equation (22) and Equation (25) may be tending toward reacceleration of the universe, in line with a discussion given by Thanu Padmabhan in  which is simple and elegant. For additional work, we need to refer to  and the idea of a “jerk” calculation   where re-acceleration of the universe is in tandem with the Jerk becoming negative valued, signifying re acceleration.
The author in  specifically did calculations as to the Jerk (re-acceleration) lead to a speed up of cosmological expansion, attributing this development to DE. What the author sees as a necessary extension of his research is to reconcile  - ; with the given phenomenology of  where reacceleration of the universe is stated to be a fact which should be investigated in terms of the cosmological constant, that is, if Equation (22) and Equation (25) are giving the same information. That has to be confirmed via rigorous calculations. Also, look at , and .
9. Dark Energy Density as Affected by Equation (22) and Equation (25)
We have also the discussion given in  and set as spacetime curvature, and so then we have
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 (26) we have within the bubble of spacetime This value for at just about the surface of a bubble.
We can then use a minimum timestep defined by using  . If so then, if the initial volume is of the cube of a Planck length, we have a time step defined via, if the value, were small, we would be probably be looking at
A limiting solution, of all things turns out to be if the curvature term is small, with, among other approximations, to be, surprise, to have the bubble tension term dropping out. But this is not a GENERAL solution.
10. Does a Physics Equivalence between Equation (28) and Equation (29) Tie in with the Physics Stated in Equation (8)?
This really has to be ascertained. Equation (8) result is a very clean derivation involving a so called “trivial” solution to the massive gravity problem, which leads to yet another question to be confirmed, if possibly in experimental data sets. Does the following hold? In terms of massive Gravity? It is from  , Novello.
The above formula (30) has to be confirmed and in some sense proved, as being not in conflict with Equation (22) and Equation (25) of this document.
11. Does the Following Arguments Give Equivalent Information, to What We Obtain Comparing Equation (22) and Equation (25)?
See this final take away as to what the cosmological constant is equivalent to, i.e. does the following make sense? If we multiply Equation (31) say by the cube of a Planck volume, and then can we work with the following derivation?
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. The main thing, can this be confirmed as far as initial GW?
The above is with regards to Relic GW generated at the start of inflation.
12. Summing Up and a Path Forward
1) I am expecting a strain of Gravitational wave strength of h ~ 10−25 to 10−26;
2) The polarization states to watch, indicating if possible, would be for massive gravitons. The author is trying to ascertain massive graviton formation via the confirmation of the result given by Novello, in  about massive gravity, and confirmation of  in terms of a weak modification of GR, if feasible. Upon confirmation of that, we will next go to…
3) The ultimate goal would be also to confirm if the multiverse generalization of the Penrose Cyclic Conformal cosmology which the author introduced in  is verifiable and does away with the anthropic principle. i.e. we have a statistical averaging procedure for keeping the Planck constant, namely the same cycle per cycle;
4) As mentioned above, in the beginning of this document also asks about the initial Penrose singularity. For reasons we will state in a future publication, we hold that if the multiverse generalization of the Cyclic Conformal cosmology of Penrose is confirmed, it may aid in falsifying the Penrose singularity. In particular, we hope we can confirm, one way or another   so far as to obtain conditions for quantum gravity which may, in part reflect the conditions given in  as to primordial GW states and their polarization commensurate with being able to answer if gravity is classical and /or quantum in nature;
5) Finally, we refer to Figure 2 as far as reproduction of structure in the universe as given by , and we ask the foundational question of questions as follows;
6) In  and .
We have wormhole geometries. Is the initial configuration of the Universe in terms of energy conditions, i.e. the negative energy condition at the throat of a worm hole, in any case even partly relevant to the Pre Planckian space-time situation so outlined in this paper?
Finally, as we implied in referencing  we find that in order to do away with the Anthropic principle, the following references in terms of Ergodic mixing of the Partition function of the Universe was utilized, as far as a multiverse.   . This is by all accounts a classical physics mixing procedure, as is well understood. Does this application of Ergodic mixing in any way negate the existence of initial conditions pertinent to Quantum gravity?
13. Conclusion. Do Virtual Images of Prior Supermassive Black Holes on the CMBR Have a Counterpart in Gravitational Wave Signatures and Signals?
What we can focus upon in part is the initial phases of when we look at the simpler Starobinsky model as to , as to a simpler Cosmology without a Penrose singularity as a start point, as seen in , as a compliment to what we have been referencing.
With the Singularity, at the start of Space-time, we have that we will have to revert to Figure 1, of this manuscript, and  and likely have only the CMBR signatures as a record of the “ghost” black holes, i.e. no GW values could possibly be obtained from the CMBR, whereas if we have a nonsingular start to expansion of the universe, there are mechanisms which exist which may allow for GW signatures which would exist.
Keep in mind this means that we will likely use the DE construction from , whereas the Penrose suggestion we are using can be reviewed using  - .
We will review next the variant of CCC cosmology used in our model, and conclude with final remarks as to the possibility of obtaining GW signatures of the “phantom” or pre universe black holes.
14. A Multiverse Generalization of CCC Cosmology
We are extending Penrose’s suggestion of cyclic universes , We are extending Penrose’s suggestion of cyclic universes, black hole evaporation, and the embedding structure our universe is contained within, This multiverse embeds BHs and may resolve what appears to be an impossible dichotomy. The following is largely taken from  and has serious relevance to the final part of the conclusion. That there are no fewer than N universes undergoing Penrose ‘infinite expansion’ (Penrose)  contained in a mega universe structure. Furthermore, each of the N universes has black hole evaporation, with the Hawking
Figure 4. Which is from Mukhanov .
radiation from decaying black holes. If each of the N universes is defined by a partition function, called , then there exist an information ensemble of mixed minimum information correlated as about 107 - 108 bits of information per partition function in the set , so minimum information is conserved between a set of partition functions per universe .
However, there is non-uniqueness of information put into each partition function . Furthermore Hawking radiation from the black holes is collated via a strange attractor collection in the mega universe structure to form a new big bang for each of theN universes represented by . Verification of this mega structure compression and expansion of information with a non-uniqueness of information placed in each of theN universes favors ergodic mixing treatments of initial values for each of N universes expanding from a singularity beginning. The value, will be using (Ng, 2008) . . How to tie in this energy expression, utilized in the partition function will be to look at the formation of a nontrivial gravitational measure as a new big bang for each of theN universes as by the density of states at a given energy for partition function .
Each of identified with Equation (35) above, are with the iteration for N universes  (Penrose, 2006) Then the following holds, by asserting the following claim to the universe, as a mixed state, with black holes playing a major part.
For N number of universes, with each for j = 1 to N being the partition function of each universe just before the blend into the RHS of Equation (36) above for our present universe. Also, each of the independent universes given by are constructed by the absorption of one to ten million black holes taking in energy. i.e. (Penrose)  . Furthermore, the main point is similar to what was done in  in terms of general ergodic mixing.
What is done in Claim 1 and Claim 2 is to come up as to how a multi dimensional representation of black hole physics enables continual mixing of spacetime   largely as a way to avoid the Anthropic principle, as to a preferred set of initial conditions. Claim 2 is to use what is known as CCC cosmology, which can be thought of as the following. First. Have a big bang ( initial expansion) for the universe. FOR SAY redshift z = 10, a ten billion years ago, SMBH formation starts. Matter-energy is vacuumed up by the SMBHs, then
if there is quintessence, constant otherwise (38)
is, here a constant. Then the main methodology in the Penrose proposal has been evaluating a change in the metric by a conformal mapping to
Penrose’s suggestion has been to utilize the following .
In fall into cosmic black hopes has been the main mechanism which the author asserts would be useful for the recycling apparent in Equation (40) above with the caveat that is kept constant from cycle to cycle as represented by We claim that the invariance of the Planck combined with Equation (36) above gives a good indication of a uniform mass to a graviton, per cycle, as far as heavy gravity, provided that holds’.
15. GW Signatures from SM Black Holes from a Prior Universe as a Way to Obtain a Proof of CCC Cosmology
To do this, would mean we review in dept  as mentioned earlier. Do the conditions for  show up at the surface of the space-time bubble?
The question delves in the matter of how CMBR patterns reflect GW values. Meaning that a future research project has to be initiated as to do we have conditions at the surface of a nonsingular start to the universe for generating GW waves?
From  we have the following as to generation of GWs, and CMBR.
Quote from page 23
LENSING-INDUCED B MODES Above we showed that density perturbations do not induce a curl in the polarization. However, that derivation assumed only linear perturbations, in which each Fourier mode of the density field is considered independently. B modes may still arise from the density at higher order. The largest relevant nonlinear effect in the CMB is weak gravitational lensing (or cosmic shear) of the primordial CMB temperature-polarization pattern by density perturbations between us and the CMB surface of last scatter (Zaldarriaga & Seljak 1998, Lewis & Challinor 2006). Lensing displaces the temperature and polarization from a given direction θ at the surface of last scatter to an adjacent position θ + δθ:
End of quote
Do we have additional density perturbations? On the surface of the “bubble” of initial space-time? If the conditions of the quote, as from page 23 of  are met, we likely have GW signals.
That should be doable as an experimental project on the surface of a bubble, provided higher order perturbative modes are identified.
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.
 Kamionkowski, M. and Kovetz, E.D. (2016) The Quest for B Modes from Inflationary Gravitational Waves. Annual Review of Astronomy and Astrophysics, 54, 227-269.
 Jeong, D.H. and Kamionkowski, M. (2020) Gravitational Waves, CMB Polarization, and the Hubble Tension. Physical Review Letters, 124, Article ID: 041301
 Giovannini, M. (1999) Production and Detection of Relic Gravitions in Quintesential Inflationary Models. Physical Review D, 60, Article ID: 123511.
 Low, H.B. and Xiong, C. (2018) Fantasia of a Superfluid Universe, as Part of the Special Volume Put out by World Scientific entitled—In Memory of Kerson Huang. World Scientific, Memorial Volume for Kerson Huang, 55-60.
 Huang, K., Low, H.B. and Tung, R.S. (2012) Scalar Field Cosmology II: Superfluidity, Quantum Turbulence, and Inflation. International Journal of Modern Physics A, 27, Article ID: 1250154.
 Tolley, A.J. (2015) Cosmological Applications of Massive Gravity. In: Papantonopoulos, E., Ed., Modifications of Einstein’s Theory of Gravity at large distances, Vol. 892 Springer, Cham, 203-224.
 D’Amico, G., de Rham, C., Dubosvsky, S., Gabadadze, G., Pirtskhalava, D. and Tolley, A.J. (2011) Massive Cosmologies. Physical Review D, 84, Article ID: 124046.
 Press, W.H., Teukolsky, S.A., Vetterling, W.T. and Flannery, B.P. (2007) Chapter 4. Integration of Functions. In: Numerical Recipes: The Art of Scientific Computing, 3rd Edition, Cambridge University Press, New York.
 Freeze, K., Brown, M. and Kinney, W. (2012) The Phantom Bounce, A New Proposal for an Oscillating Cosmology. In: Mercini-Houghton, L, and Vaas, R., Eds., The Arrows of Time, Vol. 172, Springer, Berlin, Heidelberg, 149-156.
 Abhay, A., Brajesh, G., Jeong, D.H. and Sreenath, V. (2020) Alleviating the Tension in the Cosmic Microwave Background Using Planck-Scale Physics. Physical Review Letters, 125, Article ID: 051302.
 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.
 Freese, K., Rindler-Daller, T., Spolyar, D. and Valluri, M. (2016) Dark Stars: A Review. Reports on Progress in Physics, 79, Article ID: 066902.
 Peebles, P.J.E. and Ratra, B. (2003) The Cosmological Constant and Dark Energy. Reviews of Modern Physics, 75, 559-606.
 Visser, M. (2004) Jerk, Snap, and the Cosmological Equation of State. Classical and Quantum Gravity, 21, 2603-2616.
 Beckwith, A. (2019) Identifying a Kaluza Klein Treatment of a Graviton Permitting a Deceleration Parameter q(z) as an Alternative to Standard DE. Journal of High Energy Physics, Gravitation and Cosmology, 5, 193-217.
 Corda, C. (2008) Primordial Production of Massive Relic Gravitational Waves from a Weak Modification of General Relativity. Astroparticle Physics, 30, 209-215.
 Ahmedov, H. and Aliev, A.N. (2011) More on New Massive Gravity: Exact Solutions. Physical Review Letters, 106, Article ID: 021301.
 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.
 Hajian, A. (2010) Are There Echoes From The Pre-Big Bang Universe? A Search for Low Variance Circles in the CMB Sky. The Astrophysical Journal, 740, 52.
 DeAbreu, A., Contreras, D. and Scott, D. (2015) Searching for Concentric Low Variance Circles in the Cosmic Microwave Background. Journal of Cosmology and Astroparticle Physics, 2015, 031.
 Gurzadyan, V.G. and Penrose, R. (2013) On CCC-Predicted Concentric Low-Variance Circles in the CMB Sky. The European Physical Journal Plus, 128, Article No. 22.
 Beckwith, A. (2014) Analyzing Black Hole Super-Radiance Emission of Particles/Energy from a Black Hole as a Gedankenexperiment to Get Bounds on the Mass of a Graviton. Advances in High Energy Physics, 2014, Article ID: 230713.