Next Article in Journal
Reflections on the Covariance of Modified Teleparallel Theories of Gravity
Next Article in Special Issue
Multiverse Predictions for Habitability: Fraction of Planets that Develop Life
Previous Article in Journal
Detecting the Hadron-Quark Phase Transition with Gravitational Waves
Previous Article in Special Issue
Time Reversal Symmetry in Cosmology and the Creation of a Universe–Antiuniverse Pair
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Multiverse Predictions for Habitability: Number of Potentially Habitable Planets

by
McCullen Sandora
1,2
1
Institute of Cosmology, Department of Physics and Astronomy, Tufts University, Medford, MA 02155, USA
2
Center for Particle Cosmology, Department of Physics and Astronomy, University of Pennsylvania, Philadelphia, PA 19104, USA
Universe 2019, 5(6), 157; https://doi.org/10.3390/universe5060157
Submission received: 14 May 2019 / Revised: 18 June 2019 / Accepted: 21 June 2019 / Published: 25 June 2019
(This article belongs to the Special Issue The Multiverse)

Abstract

:
How good is our universe at making habitable planets? The answer to this depends on which factors are important for life: Does a planet need to be Earth mass? Does it need to be inside the temperate zone? are systems with hot Jupiters habitable? Here, we adopt different stances on the importance of each of these criteria to determine their effects on the probabilities of measuring the observed values of several physical constants. We find that the presence of planets is a generic feature throughout the multiverse, and for the most part conditioning on their particular properties does not alter our conclusions much. We find conflict with multiverse expectations if planetary size is important and it is found to be uncorrelated with stellar mass, or the mass distribution is too steep. The existence of a temperate circumstellar zone places tight lower bounds on the fine structure constant and electron to proton mass ratio.

1. Introduction

This paper is a continuation of [1], which aims to use our current understanding from a variety of disciplines to estimate the number of observers in a universe N obs , and track how this depends on the most important microphysical quantities such as the fine structure constant α = e 2 / ( 4 π ) , the ratio of the electron to proton mass β = m e / m p , and the ratio of the proton mass to the Planck mass γ = m p / M p l . Determining these dependences as accurately as possible allows us to compare the measured values of these constants with the multiverse expectation that we are typical observers within the ensemble of allowable universes [2]. In doing so, there remain key uncertainties that reflect our ignorance of what precise conditions must be met in order for intelligent life to arise. Rather than treating this obstacle as a reason to delay this endeavor until we have reached a more mature understanding of all the complex processes involved, we instead view this as a golden opportunity: since the assumptions we make alter how habitability depends on parameters, sometimes drastically, several of the leading schools of thought for what is required for life are incompatible with the multiverse hypothesis. Generically, if we find that our universe is no good at a particular thing, then it should not be necessary for life because, if it were, we would most likely have been born in a universe which is better at that thing. Conversely, if our universe is preternaturally good at something, we expect it to play a role in the development of complex intelligent life, otherwise there would be no reason we would be in this universe. The requirements for habitability are in the process of being determined with much greater rigor through advances in astronomy, exoplanet research, climate modeling, and solar system exploration, so we expect that in the not too distant future our understanding of the requirements for intelligent life will be much more complete. At this stage of affairs, then, we are able to use the multiverse hypothesis to generate predictions for which of these habitability criteria will end up being true. These will either be vindicated, lending credence to the multiverse hypothesis, or not, thereby falsifying it.
In this work, we define the habitability of a universe H as the total number of observers it produces. In estimating the number of observers the universe contains, a great many factors must be taken into consideration. Thankfully, there has long been a useful way of organizing these factors: the Drake equation, which in a slightly modified form along the lines of [3] reads
H = λ min d λ p IMF ( λ ) × N × f p ( λ ) × n e ( λ ) × f bio ( λ ) × f int ( λ ) × N obs ( λ ) .
Here, N is the number of stars in the universe, f p is the fraction of stars containing planets, n e is the average number of habitable planets around planet-bearing stars, f bio is the fraction of planets that develop life, f int is the fraction of life bearing worlds that develop intelligence, and N obs is the number of intelligent observers per civilization. We have included dependence of the size of the host star λ = M / ( ( 8 π ) 3 / 2 M p l 3 / m p 2 ) , in order to more accurately reflect the fact that these quantities may depend on this. We then integrate over the stellar initial mass function given in [4], which approximates a broken power law with turnover at 0.2 M . This can be related to the probability of being in our universe by incorporating the relative occurrence rates of different universes as P p prior H . It was argued in [1] that a reasonable choice of prior is given by p prior 1 / ( β γ ) . However, this is ultimately set by physics at high energies, and so may in principle be something else.
Previously, the fact that the strength of gravity γ can be two orders of magnitude higher caused the biggest problems with finding a successful criterion, since most stars are in universes with stronger gravity. Though we had set out to focus solely on the properties of stars, it was only when we weighted the habitability of a system by the total entropy processed by its planets over its entire lifetime that we hit upon a fully satisfactory criterion. This was also reliant on the condition that starlight be in the photosynthetic range, colloquially referred to as ‘yellow’ light here: conservatively, this corresponds to the 600–750 nm range. Relaxing this range to be from 400–1100 nm will not qualitatively affect our results. The other factors we considered may be freely included at will without hindering this conclusion. For the majority of this paper, we take the entropy and yellow light conditions as our baseline minimal working model, and incorporate factors that influence the availability and properties of planets to determine how these alter our estimates for habitability. While our previous analysis was not heavily reliant on cutting edge results from the field of astronomy, our understanding of planets, from their population statistics to their formation pathways, has undergone rapid expansion in the past decade, and a state-of-the-art analysis needs to reflect that.
To this end, we begin by estimating the fraction of stars with planets in Section 2. Most notable is the recent determination of a threshold metallicity below which rocky planets are not found [5], as well as the understanding of the origin of this threshold. We also find the conditions for the lifetime of massive stars to be shorter than the star formation time and the fraction of galaxies able to retain metals to be sizable, but find that these only impose mild constraints. Using these allows us to determine the dependence of f p on the underlying physical parameters. Additionally, we incorporate the fraction of systems that host hot Jupiters into our analysis, and find conditions for this process to not wreck all planetary systems.
In Section 3, we turn to the average number of habitable planets n e . Two commonly used requirements for a planet to be habitable are that it needs to be both terrestrial and temperate, and so we optionally include both of these when estimating habitability. Recent results indicate that the distribution of rocky planets peaks at 1.3 R [6,7,8,9], which is rather close to the terrestrial radius capable of supporting an Earthlike atmosphere. We track how this quantity changes in alternate universes, and what implications this effect has on the number of habitable planets in these universes. Additionally, we track the location and width of the temperate zone and compare this to the typical inter-planet spacing that results from the dynamical evolution of stellar systems, which provides a rough estimate for the probability that a planet will end up in the temperate zone. Finally, we discuss the importance of planet migration and how this changes in other universes.
An appendix is provided to collect the relevant formulas for the dependence on the physical constants on the variety of processes and quantities that are needed, in order to avoid distraction from the main text.
The overarching message we derive from this analysis is that the presence of planets is not that important in determining our location in this universe, a direct consequence of the fact that the presence of planets is a nearly universal phenomenon throughout the multiverse. Including these effects barely alters the probabilities we derived before. We find that most effects act as thresholds, serving to limit the allowed parameter range rather than alter the probability distribution of observing any particular value. We find several new anthropic bounds, including the most stringent lower bound on the electron to proton mass ratio in the literature. We find that these results are relatively insensitive to the exact models of planet formation and occurrence rate used, and so are robust to these current uncertainties.
This is not the first work to address the question of whether planets are still present for alternative values of the physical parameters. Limitations on the strength of gravity and electromagnetism imposed by the existence of habitable planets was investigated in [10], where it was found that long lived, temperate, terrestrial planets can exist over a wide range of parameter space. This was continued in [11] with the investigation of the influence of the density of galaxies on planetary stability, again finding a broad allowable parameter region. Our current work is novel in not just examining the possibility of the existence of planets with desirable properties, but also taking care to incorporate modern theories of planet formation into determining whether planets with these properties are indeed produced.
Taken together, we analyze 12 distinct possible criteria for habitability in this paper (not counting migration or the different views in the planetary size distribution we consider). Coupled to the 40 we considered in [1], this represents a total of 480 different hypotheses to compare. The quantities used to compute these are displayed in Table 1 for convenience.

2. Fraction of Stars with Planets f p

Two of the factors in the Drake equation regard the existence of planets, which will be considered in turn in this paper. The first quantity to determine is f p , the fraction of stars that form planets. Here, we represent this as a product of factors:
f p = f gal × f 2 nd gen × f Z × ( f hj ) p hj .
In succession, we have: f gal , the fraction of stars in galaxies large enough to retain supernova ejecta, f 2 nd gen , the fraction of stars born after supernova enrichment, f Z , the fraction of stars born with high enough metallicity for planets to form, and an optional f hj , the fraction of stars that do not produce hot Jupiters. Here, the exponent p hj { 0 , 1 } is introduced as a choice of whether to include this last criterion or not. The other two are not treated as optional.

2.1. What Sets the Size of the Smallest Metal-Retaining Galaxy?

The requirement for a galaxy to be habitable is that it must retain its supernova ejecta in order to reprocess it into another round of metal-rich stars [12]. This sets a minimum galactic mass by the condition that the velocity of supernova ejecta is less than the escape velocity,
v SN 2 G M ret R ret .
This is the asymptotic speed the supernova ejecta attains, and, to find this, a bit of the ejecta dynamics must be used. The initial speed can be found from energy balance [13]: this can be written in the form
v SN 0 T SN A m p ,
where the temperature of the supernova is set by the Gamow energy, which is the amount required to overcome the repulsive nuclear barrier and force fusion, T SN α 2 m p . The mass of a typical particle in the ejecta is related to the atomic number A 50 , which cannot conceivably vary by much. We find that v SN 0 0.03 . However, as the ejecta moves through the intergalactic medium, it cools and slows until it merges completely. The asymptotic speed is that at which the temperature becomes equal to (about an order of magnitude less than) the hydrogen binding energy T H 2 α 2 m e / 32 10 4 K, below which Hydrogen becomes predominantly neutral and no more cooling takes place [14].
How far does the ejecta of a supernova spread before it completely merges with the interstellar medium, and, more importantly for our purposes, why? The observed value is around 100 pc [15], and this is after the blast has gone through several successive phases. The first is known as the blast wave phase, where the ejecta spread out at their initial velocity for roughly 100 yr, traveling a total of a few pc. After the amount of interstellar material encountered rivals the initial mass of the ejecta, which occurs at d ST ( M ej / ρ gal ) 1 / 3 , the blast enters the self-similar Sedov–Taylor phase, where the blast slows and expands considerably. According to standard theory [15], the self-similarity of the dynamics dictates that the temperature of the blast falls off as T d 3 . During this phase, the velocity will decrease with distance as v ( d ) = 2 / 5 ( d ST / d ) 5 / 2 until the snowplow phase begins, at which point the speed essentially does not decrease any further. The snowplow phase occurs after the temperature reaches the molecular cooling threshold, where it expands by a factor of a few until the density of the material falls to that of the surrounding medium, at which point it is completely merged. Since the bulk of the expansion takes place in the Sedov–Taylor phase, the size of the blast will be dictated by the dynamics that take place there, and so the ultimate speed is given by
v SN = v SN 0 T SN T H 2 5 / 6 .
Using the above relations, the asymptotic speed of supernova ejecta is found to be
v SN = 2.4 α β 5 / 6 .
Now, we can find an expression for the minimum mass by using M min ρ gal R 3 . Using the density of galaxies given in the appendix, this gives
M ret M p l 3 v SN 3 ρ gal 1 / 2 = 90.1 α 3 m e 5 / 2 M p l 3 κ 3 / 2 m p 9 / 2 ,
where the coefficient in the last expression has been chosen to reproduce the observed minimal mass of M ret = 10 9.5 M [16]. As explained in the appendix, the quantity κ determines the density of galaxies in terms of cosmological parameters.
While this critical mass is important conceptually, it is not as relevant to the retention of ejecta as the gravitational potential itself, which directly sets the escape velocity of the overdensity. The fraction of initial overdensities that exceed a potential of a given strength Φ is given by the Press–Schechter formalism as f = erfc ( Φ / ( 2 Q ) ) [17], where Q is the primordial amplitude of perturbations. Usually, this expression is immediately expressed in terms of mass and density, but this more primitive form will suffice for our purposes. With this, we can derive the fraction of matter that resides in potential wells deep enough to produce a second generation of stars as
f gal = erfc 4.1 α 2 β 5 / 3 Q .
Perhaps somewhat interestingly, this does not depend on the strength of gravity γ at all. This effectively acts as a step function, severely diminishing the habitability of universes where M ret > M typ . Note that this is a pessimistic estimate, as it ignores the potential for subsequent evolution that causes potential wells to deepen with time. Nevertheless, it is a very mild bound, and so a more thorough treatment is not called for.

2.2. Is Massive Star Lifetime Always Shorter Than Star Formation Time?

Since the formation of metal-rich systems is reliant on the evolution of the first stars through to their completion, if the lifetime of massive stars exceeds the duration of star formation, then no systems will form with any substantial metallicity. The second generation stars are not necessarily enriched enough to produce planets, but this serves as a sufficient condition for planets to be formed at all. Then, the fraction of stars with planets can be estimated as those that are born after a few massive stellar lifetimes have elapsed. It is worth considering how these two timescales compare for general values of the physical constants.
The star formation rate, averaged throughout the universe, is found to decline exponentially, as gas is depleted from the initial reservoir [18]. The most naive treatment one can perform is to relate the timescale of this depletion to the free fall time of a galaxy, t dep = 1 / ( ϵ SFR G ρ ) , where for simplicity the efficiency coefficient is taken to not vary with parameters. Then, the fraction of stars born after a time t is f ( t ) = e t / t dep .
This needs to be compared to the lifetime of massive stars, where massive here is taken to mean large enough to become a type II supernova. This threshold is eight solar masses in our universe, and is set by the inner core having a high enough temperature to undergo carbon fusion [19]. As such, this threshold is parametrically similar to the minimum stellar mass, which was shown in [20] to scale as λ α 3 / 2 β 3 / 4 , which equated the Gamow energy of fusion with the internal temperature of the star. The resultant mass is about two orders of magnitude larger than the minimal mass, stemming from the larger repulsive barrier for large nuclei. Then, using the stellar lifetime from [1], we find
t SN t dep = 0.24 κ 3 / 2 α 7 / 4 β 1 / 8 γ .
The normalization has been set to the value of 0.01. This has been a bit loose on several counts, namely the assumption that the second stars are always metallic enough to form planets, and the neglect of even higher mass stars, which have correspondingly shorter lifetimes. However, in practice, these worries are of no consequence, precisely because the scales are separated by such a large amount. We find that, for all practical values of the parameters, the star formation timescale exceeds the lifetime of massive stars, so that f 2 nd gen f ( t SN ) is always close to 1.

2.3. What Is the Metallicity Needed to Form Planets?

Terrestrial planets must be formed out of heavy elements, and though even a minuscule amount would suffice in terms of actual planetary mass, the planet formation processes within the protoplanetary accretion disk can only occur after a high enough metallicity is reached. A recent analysis [5] has done a nice job characterizing the metallicity needed (as a function of stellar mass, and distance to the host star) in order for planets to form out of the initial protoplanetary disk.
The underlying physical picture is that there are two timescales, the lifetime of the disk t disk and the dust settling time t dust . If the second exceeds the first, the disk will dissipate before sufficient dust may settle into the midplane, and will disperse without forming planets. Both of these depend on metallicity monotonically, and so only above some critical value will the conditions for planetary formation hold. We detail the scaling of each in turn.
The physics of accretion disks was first laid out in [21]. An initially uncollapsed cloud first condenses, and then through angular momentum transfer begins to form a disk. The disk remains around the star until either UV light from the star photoevaporates the gas or it escapes thermally, and the disk evaporates. This phase of evolution occurs on the viscous timescale of the disk, which is much smaller than the total disk lifetime [22].
Though the precise time at which this crossover occurs depends on the exact mechanism of photoevaporation (X-ray versus extreme ultraviolet, other stellar sources in clusters, presence of high activity tau phase) the drop-off of accretion is set by the disk’s secular evolution time [23], with crossover occurring on the order of this time. Therefore, the only relevant physics that needs to be kept track of is the free fall time. Since this is given by a cloud that has approximately virialized after Jeans collapse, it is simply set by the condition that the gravitational energy is equal to the thermal energy. Then, the accretion rate is given by M ˙ = c s 3 / G 10 6 M / yr , and the timescale is given by t disk = M / M ˙ Myr. The sound speed and the temperature are given by bremsstrahlung from molecular line cooling, as found in the appendix. Then,
t disk = 5.5 × 10 5 λ m p 1 / 4 M p l α 3 m e 9 / 4 Z Z p .
In [24], the metallicity dependence was studied, where it was found that cooling quickens with metallicity Z. The exact process was open to interpretation, with a handful of viable candidate mechanisms, but the overarching explanation is that the less shielding there is, the faster the disk will dissipate. Observationally and theoretically, it was found that the scaling with metallicity is consistent with p = 1 / 2 , and so we will adopt this for our analysis.
The dust settling timescale is dictated by the rate at which dust grains sink into the midplane from their initial positions in the protoplanetary disk. It is given in terms of the Keplerian timescale of the disk, t dust = 0.72 t Kepler / Z , as derived in [5]. There it was related to the growth rate of grains, and found to depend inversely on the metallicity, since only dust can participate in the accretion process. This is a function of orbital distance, but specifying to planets that form within what will become a temperate orbit for simplicity, we have
t dust = 3.0 × 10 5 λ 17 / 8 m p 7 / 4 M p l 1 / 4 α 15 / 2 m e 3 Z Z .
Here, we are explicitly assuming that significant migration does not occur, which would alter the metallicity needed for planet formation, given its dependence on orbit.
Equating these two timescales yields the critical metallicity
Z min = 4.2 × 10 4 λ 3 / 4 γ 1 / 2 α 3 β 1 / 2 .
The critical value is found to be Z min = 6.3 × 10 4 . As compared to the solar value Z = 0.02 , we find Z min = 0.03 Z [5].
This may be used to determine the fraction of stars that host planets by considering the amount of stars that are formed above this metallicity. This requires a model for how metallicity builds up inside a galaxy (of a given mass, as well as the distribution of galaxy masses). A full calculation of this sort would take us too far afield here so we will return to it in a later publication. However, a very reasonable approximation is to compare this threshold metallicity to the asymptotic metallicity a galaxy attains—this usually affords percent-level accuracy. In our universe, this is found to be Z = 0.011 [16,25], which is actually a factor of two below solar due to the intrinsic scatter in metallicities. This asymptotic value is set by the fraction of stellar mass that is transformed into heavy elements during stellar fusion, and is not expected to depend sensitively on the fundamental constants over the ranges considered.
If attention is restricted to solar mass stars, the requirement that the threshold metallicity be lower than the asymptotic value equates to α 3 β 1 / 2 γ 1 / 2 > 0.80 , which is most sensitive to the fine structure constant. Leaving the other parameters fixed, this gives α > 1 / 361 , which is stronger than the bound α > 1 / 685 based on galactic cooling found in [26,27]. A more forgiving bound is found using the smallest stellar mass, which is α 5 / 2 β 17 / 12 γ 2 / 3 > 0.32 . This latter boundary is displayed in Figure 1, where the distribution of observers is plotted as a function of the three variables α , β and γ .
Incorporating the above effects determining the fraction of stars that host planets, and again using the entropy and yellow light conditions, we can find that the probability of observing our particular values of each constant. These are defined as P ( x obs ) = min { P ( x > x obs ) , P ( x < x obs ) } for any observable quantity x, the others being integrated over, and the probability of measuring any value given by Equation (1). With this, we have
P ( α obs ) = 0.19 , P ( β obs ) = 0.45 , P ( γ obs ) = 0.32 .
These values are indistinguishable from those that were found in [1] without including these effects. We conclude that the presence of planets is a generic feature throughout most of the multiverse, and it does not alter where we expect to be situated in the slightest1.

2.4. Are Hot Jupiter Systems Habitable?

The constraints above were all quite mild, indicating that the presence of planets is a fairly generic feature of the multiverse. However, we can make a further refinement by incorporating one of the most famous statistical correlations in the field of exoplanets, the hot Jupiter–metallicity correlation [28]. This finds that the fraction of stars that possess hot Jupiters, that is, Jupiter sized planets on orbits extremely close to the star increase with metallicity as Z 2 . The general (though not universal— [29]) consensus is that these planets must have formed in the outer system before moving inward, to avoid the necessity of a disk that would be so massive as to be unstable [30]. In this scenario, the migration of the planet through the inner solar system would have certainly ejected any preexisting planets from their orbits (or worse), precluding them from sustaining life. However, the authors in [31] hypothesize that this entire process could happen early enough that the main stage of planet formation could occur after this migration had already taken place. If this turns out to be the case, then there may in fact be no correlation between hot Jupiters and habitability (barring other factors that may impact habitability [32]). If one wishes to include this effect, however, then the fraction of stars that host Earths and not Hot Jupiters is given by
f hj = 1 Z ¯ 2 Z max 2 .
Here, we have used the mean metallicity rather than averaging this fraction over the metallicity distribution, but this approximation is sufficient for a first analysis. The normalization Z max is the threshold above which the stellar system is almost assured to possess a hot Jupiter. It is set to reproduce the observed abundance of hot Jupiter systems of 3 % as Z max = 5.77 Z . In a multiverse setting, we would expect to inhabit a universe where Z max is safely above the average metallicity, beyond which any further increase would result in little increase in habitability. Somewhat in line with this expectation, then, is the fact that in our universe hot Jupiters exist in only a few percent of systems, and mainly in those that are highly metal-rich.
To investigate whether this is the result of some selection effect, we must know what determines this metallicity. The functional dependence is a clue: as the effect becomes more pronounced with the square of the nongaseous material present, this is indicative of an interaction process. What remains, however, is the question of whether this migration is a result of planet–planet or planet–disk interactions. In fact, both explanations have been considered in the literature: references may be found in [33] 2. The planet–planet hypothesis is supported by the fact that the eccentricities of observed hot Jupiters are correlated with metallicity as well, indicating a more chaotic, violent origin, rather than the steady, deterministic process indicative of planet–disk interaction. Additionally, [35] note a substantial misalignment between the orbits of known hot Jupiters and the spins of their host stars, which is most easily explained through a violent migration scenario. However, systems like WASP-47 [36], which possess both a hot Jupiter and smaller companions, demonstrate that more dynamically quiet migration pathways are possible, if not necessarily the norm. Here, we only expound upon the planet–planet scenario, though the others could just as readily be incorporated into our analysis.
We start by determining the value of Z max due to planet–planet interactions, as proposed in [37]. As noted in [38], the coexistence of hot Jupiters and low mass planets is impossible in this paradigm, as migration occurs after the disk has dissipated. From [33], this is set by the expected number of Jupiter mass planets initially formed in the outer system. They provide a framework for estimating this by determining the probability that a core will attain Jupiter mass as a result of planetesimal accretion as p Δ M / M crit , where M crit 10 M is the mass above which runaway gas accretion is possible and Δ M is the typical total accreted mass. For this, we use the analytic expression for the accretion rate from [39],
M ˙ ( 18 π ) 1 / 3 G M crit M 1 / 3 Σ planetesimals ρ 1 / 3 a .
This employs the strong gravitational focusing limit, and treats the typical relative velocities of planetesimals as roughly given by the Hill velocity. This can be used to determine the total mass accreted by simply multiplying by the disk lifetime given in Equation (10) (making use of the simplifications that the nonlinear oligarchic regime is not quite reached, and the initial isolation timescale is small compared to the disk lifetime). The probability that there will be at least two gas giants to trigger the instability will scale as p ( 2 ) N jup 2 p 2 , where N jup is the typical number of planets in the outer system, and we have assumed that p is small. The quantity N can be found by dividing the total mass of the planetary disk by the typical mass of a planet at the typical location of formation. Here, we use the initial seed being set by the isolation mass, have fixed the orbital radius to be given by the snow line, and have taken the disk temperature to be given by viscous accretion, all of which are discussed in the appendix. This gives
N jup M disk M crit 2.2 × 10 5 λ α 3 / 2 β 3 / 4 .
This scales linearly with stellar mass, in agreement with the observations in [40]. The maximal value of metallicity is found to be
Z max 1 N jup M crit M ˙ t disk = 2.3 × 10 8 α 13 / 6 β 3 / 2 κ λ .
This quantity is somewhat sensitive to both α and β , but not at all to γ . This also defines a stellar mass λ hj = 2.3 × 10 6 α 13 / 6 β 3 / 2 / κ : stars above this mass, equal to 11 M for our values, will always host hot Jupiters. This will be below the smallest stellar mass if 789 α 8 / 27 β < 1 , which will occur when the electron to proton mass ratio is about 10 times smaller. However, this criteria does not alter the probabilities much:
P ( α o b s ) = 0.18 , P ( β o b s ) = 0.44 , P ( γ o b s ) = 0.31 .
With this, notice that the probabilities are only changed from Equation (13) by a few percent. Thus, even when including the demand that hot Jupiter systems are uninhabitable, the fraction of systems with planets seems to be relatively insensitive to the physical constants.
We now turn to the next factor in the Drake equation, which deals with the characteristics of planets, rather than just their presence.

3. Number of Habitable Planets per Star n e

Next, we focus on the number of habitable planets per star, n e . The determination of habitability may depend on many factors, such as amount of water, eccentricity, presence of any moons, magnetic field, distance from its star, atmosphere, composition, etc. Here, we focus on two: temperature and size, and determine the fraction of stars that have planets with each of these characteristics.
As usual, it is possible that habitability is completely independent from these properties; which viewpoint one adopts depends on how habitable one expects environments without liquid surface water and thin atmospheres can be. In this work, we remain agnostic to either expectation and report the number of observers for all combinations of choices, where a planet will only be habitable if it is approximately Earthlike, and where the size and/or temperature of the planet have no effect on its habitability.
It should be noted that we are restricting our attention here to surface dwelling life on planets orbiting their star. Thus, life in subsurface oceans and/or on icy moons like Enceladus [41], or even more exotic types of life (e.g., [42,43]), are not considered. It is our plan to consider these alternative environments in future work.
The number of habitable planets can be broken down as
n e = n p × ( f terr ) p terr × ( f temp ) p temp .
Here, the average total number of planets around a star is n p . The fraction of terrestrial mass planets is denoted f terr , and f temp is the fraction of planets that reside within the temperate zone. The exponents p i { 0 , 1 } parameterize the choice of whether to include these conditions in the definition of habitability or not. These quantities all depend on stellar mass, giving preference to large stars because they make larger planets, and small stars in that their temperate zone is wider compared to the interplanetary spacing.
Estimating these quantities is somewhat muddled by the current uncertainties in planet formation theory. Not only is the distribution of planet masses contested in the literature, but the exact formation pathways, as well as the physics that dictates the results, is not completely settled. Where we come across disagreement, we separately try each proposal, in order to understand the sensitivity of our analysis to present uncertainties. While the results for the overall probabilities can vary by a factor of 2, the upshot is that our estimates are relatively robustx.

3.1. Why Does Our Universe Naturally Make Terrestrial Planets?

The size of a planet is of crucial importance because it dictates what kind of atmosphere it can retain. If it is too small, all atmospheric gases will eventually escape, whereas if it is too large, it will retain a thick hydrogen and helium envelope, leading to a runaway growth process. Terrestrial planets must have a very specific size in order that the escape velocity exceeds the thermal velocity for heavy gases such as H 2 O, CO 2 and N 2 , but not that of the lightest gas H and He. In our universe, and for temperatures within the range where liquid surface water is possible, this restricts the range of planetary radii to be within 0.7 and 1.6 that of Earth’s [44]. This is a narrow sliver compared to the eight orders of magnitude mass range of spherical, non-fusing bodies, ranging from the potato radius of 200 km to 10 Jupiter masses. In terms of fundamental parameters, this requirement gives the mass to be
M terr = 92 α 3 / 2 M p l 3 m e 3 / 4 m p 11 / 4 .
The coefficient has been set to reproduce Earth’s mass, but the allowed spread in masses is taken to be between 0.3–4 M.
There are compelling arguments that complex may only be possible on terrestrial planets, with atmospheres composed of only heavy gases [45]. Any smaller, and the planet would be Marslike, an apparently barren wasteland incapable of sustaining any appreciable liquid ar atmosphere. The other extreme would be Neptunelike, with its hellish surface temperatures, pressures and wind speeds. Of course, these arguments may be misguided, but here we explore the consequences of adopting them for the multiverse computations we perform.
It is important to note that the conditions that set the presence of atmospheres are completely separate from the physics that dictates the size of planets, which is set by the clumping of the initial circumstellar disk3. Nonetheless, the observed population of rocky planets is thought to peak at only slightly super-Earth mass, making the production of terrestrial planets the norm for stars throughout the universe. To be fair, the current exoplanet samples are biased towards large mass planets and become very incomplete below Earth mass [48], but a number of different groups have concluded that a detectable turnover is present near Earth masses: the authors of Ref. [9] find a good fit to a log-normal distribution that peaks at 1.3 R . In Ref. [7] a Rayleigh distribution with width 3 M is used. The authors of Ref. [8] advocate for a broken power law with turnover at 5 M . It was noted in Ref. [6] that the distribution appears to be flat below 2.8 R . Use of these differing proposed distributions make very little difference to our final outcome.
However, not everyone is convinced that the mass distribution exhibits a peak, and even if there is one, it is just as reasonable to assume that there are many more smaller mass planets for every planet of Earth size, as the plethora of small asteroids and comets in our system indicates. Because of the incompleteness of current exoplanet surveys for small mass planets, there is room for such disagreement at the current moment. Additionally, even if the mass peak is real, it is only observed for close in exoplanets, and so requires an extrapolation to Earthlike orbits, where different dynamics may be at play. One possibility is that the peak at super-Earth mass may be due to their enhanced migration capability [49]. We will consider each scenario in turn.

3.1.1. What Sets the Size of Planets?

Why is the turnover so nearly equal to Earth mass planets, out of the potentially eight orders of magnitude that could have been selected instead? Simulations provide a means to address this question: it was found in Ref. [50] that the mass of planets is directly proportional to the amount of initial material present in the disk, so that increasing disk mass makes larger, rather than more, planets. In this scenario, nearly all the material initially present in the disk eventually gets constituted into planets, with negligible (perhaps a factor of two, but not an order of magnitude) losses throughout the evolution of the system. Determining the final planet mass in this setup requires knowledge of the initial disk mass, as well as the fraction of material within the inner solar system. This boundary is set by the snow line, the difference in composition interior and exterior to which dictates the formation of rocky versus icy planets. In the following, we adopt the conventional view that little migration takes place during planet formation, and comment further on alternative pathways when migration is discussed.
Current observations indicate that the initial mass of heavy elements in protoplanetary disks is roughly proportional to disk mass, M disk = 0.01 M [51]. In Ref. [52], a stronger dependence was found, but the scatter of 0.5 dex is larger than the trend, and the observed trend was suggested to possibly be due to a selection effect arising from processing into undetectably large grains, so we omit this stronger scaling for now. For solar mass stars, this works out to be roughly 10 M Jupiter , or 3300 M ⨁. . This mass is distributed out to a radius of ∼100 AU, which is set by the conservation of angular momentum, and the initial size of the collapsing cloud. From the appendix, we arrive at the following expression for disk size:
r disk = 3.6 × 10 6 λ 1 / 3 κ M p l m p 2 .
This may be significantly altered if the young star is in a dense environment [53], but this scaling will suffice for our purposes.
To determine the amount of material present within the snow line, one must take the surface density profile into account. For this, we use the expression for Σ ( a ) found in the appendix, which is inversely proportional to a. Using these, the typical mass of an interior planet becomes
M inner η < a snow R disk M disk ,
where we have included the quantity η < 0.25 to account for the difference in composition interior to the snow line [54].
The location of the snow line is the point in the disk beyond which water condenses, equal to 2.7 AU in our solar system [55]. In order to determine its location, the temperature as a function of radius must be used, which will depend on the dominant source of heating. For most of the present paper, we use the temperature given by viscous accretion, and find
a snow = 30.4 λ M p l α 5 / 3 m e 5 / 4 m p 3 / 4 .
Note that, perhaps unsurprisingly, the snow line is situated outside the temperate zone for all relevant parameter values. The dependence on stellar mass was found taking M ˙ λ 2 , but is generically found to scale as a snow λ 2 / 3 2 [56]. This is now enough to determine the parameter dependence of M inner , which will ultimately be used to derive the expected number of planets per star as a function of these parameters.
Though we tend to favor accretion dominated disks throughout this work, irradiation from the central star can actually play a significant role as well [22]. If this is the dominant mode of heating, then the snow line will instead be given by the expression a snow = 297 λ 43 / 30 α 4 m e 4 / 3 m p 1 / 3 M p l 2 / 3 . Though this is functionally quite similar to the accretion dominated case from above, in Table 2, we investigate the effect of assuming this form instead. Actually, both are almost equally relevant in determining the position of the snow line, which helps to greatly complicate the disk structure (as well as enhance the variability between different star systems [57]). Additionally, irradiation from other neighboring stars may be important as well, especially in clusters [58], but we do not consider this contribution in this work.
Determining the average planet mass is somewhat involved, given the many distinct stages of growth that occur as microscopic dust grains agglomerate to the size of planets. For reviews of this multi-stage process, see [59,60,61]. In brief, planetesimals form characteristic masses as the final outcome of pebble (or dust) accretion. This arrangement is unstable, and ultimately leads to a phase of giant impacts, wherein planetesimals collide together to form full planets.
An estimate for the maximum mass a planet can attain after a phase of growth through chaotic giant impacts was found in [62]. There, they assumed no migration, small eccentricity, and determined the width of the ‘feeding zone’ to be Δ a = 2 v esc / Ω by noting that within this region planet–planet interactions result in collisions rather than velocity exchange. This yields
M planet 4 π Σ a 5 / 2 ρ 1 / 6 M 1 / 2 3 / 2 ,
where ρ is the average density of the planet and Σ is the disk density. With this, we can use the appendix to reinstate parameter dependence into the expressions for the average number of planets (normalized to 3 for the solar system) as well as the typical planet mass:
n p = 0.0061 α 4 / 3 β 13 / 16 κ 1 / 2 λ 5 / 6 ,
M planet M terr = 1.6 × 10 6 κ 3 / 2 λ 5 / 2 α 9 / 2 β 45 / 16 .
The latter has quite a steep dependence on stellar size. This is expected from the simulations [50], with a dependence closer to linear, but is not particularly observed in exoplanet catalogs due to large scatter [63]. In our calculations, we explore the effect of ignoring this dependence altogether, displayed in Table 2. A full treatment would take the dependence on semimajor axis into account, rather than simply evaluating at the snow line: in fact, this would be somewhat unnecessary, as the scatter observed in exoplanet surveys, simulations, and indeed within the solar system masks any dependence that may exist.
For the fraction of planets that are terrestrial, we use the log-normal distribution of [9] since this is what is expected of the core accretion process, though we check that the actual distribution used does not alter the outcome much. Under these assumptions, we can find the probability that a planet will be terrestrial as
f terr = 1 2 erf log 4 M terr M planet + 1 2 σ M 2 2 σ M 1 2 erf log . 3 M terr M planet + 1 2 σ M 2 2 σ M .
This function peaks at M planet M terr , and approaches are when M planet is very different from M terr . Being a two-parameter distribution, this requires not just the mean, but also the variance. As it is not currently known what sets this quantity, here we explore two options: the first uses maximization of entropy production to set σ M = 1 / 6 , which is fairly widely observed in natural processes [64]. This estimate should occur for large systems, but for small systems one would expect the variance to be set by shot noise instead, σ M iso / M inner , making use of the isolation mass defined in the next section. The dependence on the various parameters is displayed in Figure 2.
The probabilities of observing the observed values of our constants are computed for the various choices we made in Table 2. These can be compared with Equation (13) that only considered the fraction of stars with planets. The most significant change for our most favored prescription is the probability of observing our electron to proton mass ratio, which is decreased by less than a factor of 2. Such insensitivity hardly constitutes any evidence for whether life should only appear on terrestrial planets. More interestingly, if the dependence on stellar mass is neglected, our strength of gravity becomes quite uncommon to observe. This gives us strong reason to suppose that, if life requires terrestrial planets, we will begin to see a correlation between the two soon, and, if we don’t, then life requiring Earth mass planets is incompatible with the multiverse at the 2.7 σ level.
Aside from this, though, the choices we made to come up with these estimates do not affect the outcome very much at all. On the one hand, this is disappointing, as the stronger the dependence these probabilities have on the assumptions of planet formation, the stronger our predictions can be about which to expect to be dominant. However, this is also heartening: because the current uncertainties about planet formation do not affect the outcome all that much, we are able to trust the broad conclusions we have reached a bit better.

3.1.2. Is Life Possible on Planetesimals?

The mass discussed above really refers to the maximal planet mass of a system, which form as a result of the secondary stage of collisions after the isolated planetesimals form. However, this agglomeration will likely not completely deplete the system of its primordial planetesimals, and so there are also expected to be numerous smaller planets accompanying each large one, as is the case in and around the solar system’s asteroid belt. If these smaller bodies are considered as potential abodes for life as well, the distribution continues past Earth masses, rather than having a peak there.
As a planet is condensing out of the protoplanetary disk, it eventually reaches what is termed as the pebble isolation mass, which from the appendix is given by
M iso = 2.0 × 10 8 κ 3 / 2 λ 2 M p l 3 α 5 / 2 m e 15 / 8 m p 1 / 8 .
The isolation mass is a function of the semi-major axis, but if we evaluate it at the edge of the inner system given by Equation (23), we find, using the value for the disk surface density from the Appendix A,
M iso M terr = 2.2 × 10 6 κ 3 / 2 λ 2 α 4 β 21 / 8 .
Here, the density of the galaxy comes into play in setting the outer edge of the disk. The dependence on stellar mass in this expression is quite close to that found in [55], M iso λ 7 / 4 .
It remains to specify the distribution of planetary masses in order to find the fraction that are terrestrial in this picture. It is generically expected to take a power law form that continues to the small mass cutoff, so that N ( M ) = ( M iso / M ) q . However, different authors prefer different values for the slope: Ref. [65] find q = 0.31 ± 0.2. Ref. [66] find a nearly scale invariant distribution for the radius, which translates into q = 0.30 ± 0.03 if we use M r 3 . Refs. [39,67] find q = 0.6 for simulations of planetesimals, and Ref. [54] extrapolate from the known asteroid population to find q = 1 . Ref. [68] favors a value of q = 2 from population synthesis methods. Here, we report with various values of q to investigate its influence on the probabilities. The fraction of terrestrial planets is then
f terr = min 1 , M iso 0.3 M terr q min 1 , M iso 4 M terr q .
This is an increasing function of stellar mass until M iso > 0.3 M terr , reflecting the expectation [69] that earthlike planets should be rare among low mass stars. With this prescription, stars above a certain mass will not produce any earthlike planets because the isolation mass will exceed the largest terrestrial planet size. This defines the largest viable stellar mass λ iso = 0.0013 κ 3 / 4 α 2 β 21 / 16 , which corresponds to 6.3 M in our universe. This will only exceed the minimal stellar mass if α 1 / 2 β 33 / 16 κ 3 / 4 > 169 . This condition is most sensitive to the electron to proton mass ratio and, holding the other constants fixed, will be violated if it drops to about 11 % of its observed value.
For q < 1 , the average mass for this distribution is formally infinite, which presents a problem for using our expression for the expected number of planets in a system as given by n p = M inner / M p . However, the total mass in the inner disk introduces a large mass cutoff, for which we have
n p = q 1 q M inner M iso M inner M iso 1 q 1 M inner M iso 1 q .
This interpolates between ( q 1 ) / q M inner / M iso for q > 1 and ( q 1 ) / q for q < 0 . The ratio of masses is
M inner M iso = 0.0017 α 5 / 6 β 5 / 8 κ 1 / 2 λ 1 / 3 ,
thus that this equates to 30 with our constants. With this viewpoint, the distribution of observers is plotted in Figure 3.
The overall probabilities are calculated for different representative values of the power law in Table 3. It can be seen that, for increasing q, the probability for α increases, while the other two decrease. In particular, for q = 2 , the probability of observing our value of the electron to proton mass ratio is disfavored by 2 σ .

3.2. Why Is the Interplanet Spacing Equal to the Width of the Temperate Zone?

Perhaps the most commonly employed habitability criteria is that a planet must be positioned a suitable distance away from its host star to maintain liquid water on its surface. This assumption is so pervasive that this region is usually referred to as the circumstellar habitable zone. In the spirit of remaining agnostic toward the conditions required for life, we will adhere to the recently proposed renaming as the ‘temperate zone’ [70]. If this is indeed essential for life, then the expected number of habitable planets orbiting a star will depend both on the interplanetary spacing, as well as the width and location of the temperate zone. A rather clement feature of our universe is that the width of the temperate zone is comparable to the interplanetary spacing (for sunlike stars). Because of this, it is relatively common that one of the planets in any stellar system is situated inside the temperate zone, no matter its particular arrangement. This could be contrasted to the hypothetical case where the temperate zone were much narrower than the interplanetary spacing, in which case the odds of a planet being situated inside it would be quite low. However, this coincidence of distance scales is not automatic: both these quantities are dependent on the underlying physical parameters, and so in universes with different parameter values the expected number of potentially habitable planets per star will be altered. We go through these length scales in turn, and then fold them into our estimate for the overall habitability of the universe.
The boundaries of the temperate zone depend on the characteristics of the planet in question, such as the atmospheric mass and composition [71], its orbital period [72], etc. However, these details only alter the location of the temperate zone to subleading order, and do not affect the scaling with fundamental parameters we are interested in. If the planet is assumed to be a simple blackbody, then the temperature is set solely by the amount of incident flux. In this case, the location of the temperate zone will be
a temp = 1 2 T 2 E H 2 O 2 R = 7.6 λ 7 / 4 m p 1 / 2 M p l 1 / 2 α 5 m e 2 .
Albedo and greenhouse effects will change the coefficient, but not the overall scaling. Determining the width of the temperate zone entails finding the temperatures at which runaway climate processes occur, but we will now argue that both the inner and outer edge are dictated by (broadly speaking) the same underlying physical process of phase change, and so the width of the temperate zone will scale in the same way as its location.
The inner edge of the temperate zone is set by the runaway greenhouse effect: this occurs when a temperature threshold is crossed that allows an appreciable amount of water vapor to be sustained in the atmosphere. Since this serves to trap infrared light from escaping to space, this will increase the temperature further, in turn driving a further increase in atmospheric water vapor. Once the atmosphere is comprised primarily of water, it will be photodissociated and/or escape to space, leaving the Earth in a dry and Venus-like state [73]. Since there is always some level of outward flux, the ocean will escape into space eventually if a long enough time has elapsed. Therefore, the exact threshold for this process is defined as when the timescale for this process is of the order of one billion years, which in turn will depend on the mass of the planet’s ocean. However, key to our discussion is that the change in the atmospheric water fraction occurs very abruptly when the temperature crosses the latent heat of vaporization, going from equilibrium values of 10 5 1 in the span of 250–420 K. Because of this, the exact mass of the planet, ocean or atmosphere will only play a subleading role in determining this threshold, which instead is dictated solely by molecular processes. Since this transition is set by what is essentially the intermolecular binding energy, this scales identically with the condition for liquid water. For Earth, this occurs at a temperature of 330 K, which corresponds to 0.95 AU [74].
Similarly, the outer edge of the temperate zone is set by the runaway icehouse process. The temperature of a planet is set not only by the incident flux, but also by the carbon dioxide content content of the atmosphere and albedo. Carbonate weathering is an important regulatory feedback mechanism that serves to stabilize the temperature of a planet to remain within the temperate range by adjusting atmospheric carbon dioxide content to compensate for a change in stellar flux [75] (it does this because liquid water is necessary for efficient weathering to occur). However, this only works to an extent, since beyond a critical concentration atmospheric CO 2 increases the albedo of a planet, leading to a runaway icehouse effect [73]. Because the amount of atmospheric carbon dioxide is also set by reactions due to intermolecular forces, it scales the same as above. Therefore, both the inner and outer edges of the temperate zone are set by molecular binding energies, and so the width will always be of the same order as the mean4.
The exact delineations are subject to the uncertainties of the atmospheric model used, but, for definiteness, we take Δ R HZ = 0.73 AU, from [71].
The next task is to determine how far apart planets typically reside from each other. It was hypothesized in [77] that stellar systems are dynamically packed, in that they are filled to capacity, and the insertion of any additional planet would render the system unstable. Though this may not strictly hold all the time [78], this spacing is roughly observed in our solar system [79], in Kepler data on multiplanetary systems [80], and found to occur naturally in simulations in [81]. The essential idea is that planetary scattering either ejects or collides excess planets from an initially overpacked system until the remainder are far enough apart to be stable on the timescale of the system. The stability condition is that planets are further than a certain multiple of Hill radii away, nominally around 10, and it was found in [80] that the distribution of separations was a shifted Rayleigh distribution with usual separation 21.7 ± 9.5 R Hill .
As mentioned in [62], the typical multiple of mutual Hill radii is not universal, but can be shown to depend on the width of a planet’s ‘feeding zone’ to be given by ( ρ a 3 / M ) 1 / 4 , where ρ is the density of matter and a is the semimajor axis. Though above we were interested in the typical size of a planet during this process, here we may use this condition to find the typical spacing as
a spacing 16 π Σ a 11 / 2 ρ 1 / 2 M 3 / 2 1 / 2 = 1430 κ 1 / 2 λ 169 / 48 γ 1 / 8 α 21 / 2 β 15 / 4
This was evaluated at the center of the temperate zone to give the ratio of these two length scales in terms of fundamental constants. Since this ratio loosely sets the fraction of planets within the temperate zone, we arrive at
f temp Δ a temp a spacing = 0.0053 α 11 / 2 β 7 / 4 κ 1 / 2 λ 85 / 48 γ 5 / 8 .
Set this to the value of 1.6 in our solar system. Notice that, for α or β much smaller or γ much larger, the fraction of temperate planets will be diminished. Additionally, this quantity is larger for smaller mass stars, reflecting the comparatively broad temperate zones there.
This distribution is illustrated in Figure 4. The corresponding probabilities for observing our values of the constants are
P ( α o b s ) = 0.24 , P ( β o b s ) = 0.37 , P ( γ o b s ) = 0.15 .
Though the probabilities shift around by a factor of a few, the inclusion of this criterion does not actually affect the results very much. This is a generic conclusion, even when this is included with other criteria. Whether or not this condition is essential for life, the multiverse hypothesis is relatively insensitive to it.
Let us also take this opportunity to determine what values of parameters would render disks smaller than the habitable orbit, thereby precluding temperate planets from ever forming. Based off the expressions above and in the appendix, we have
a temp R disk = 4.0 × 10 5 κ λ 17 / 12 γ 1 / 2 α 5 β 2 .
Aside from the obvious condition that, if the mean free path of the galaxy were 100 times smaller, the protoplanetary disk would be shrunk by the same amount, this also places restrictions on the fundamental parameters. Though these introduce lower bounds for α and β , this is the region with very few temperate observers anyway, and so this does not change the habitability estimate appreciably, especially since these scales would have to shift by two orders of magnitude before the threshold is reached. A more stringent boundary would require that disks be larger than the orbits of the gas giants. If not, they would be precluded from forming in the first place, and, if these were essential for life on Earth, universes like this would be altogether barren.

3.3. Planet Migration

Traditional models of planet formation assume that little orbital migration takes place during planet formation. However, the large number of gas and ice giants found situated extremely close to their host stars [82], the presence of orbital resonances found in some exoplanet systems [83], and icy composition of planets within the snow line [84] all point to the presence of migration. Migration may strongly affect the habitability of planetary systems, as if it proceeds for too long, all inner planets will drift toward the inner edge of the disk at around 0.014 AU [85], and giant planets migrating across the temperate zone would destabilize the orbits of any existing planets. The vast diversity of systems found, as well as analytic and numeric simulations of protoplanetary disks [86], point to an exquisitely complex array of migration scenarios, which will be selectively operational depending on the characteristics of the initial system such as disk mass, viscosity, and the size and locations of the planets. Nevertheless, it is possible to perform rough estimates for when migration will be present in a given system. In this section, we do not attempt to characterize all the complex features of planet migration in universes with alternate physical constants, but rather wish to provide a useful diagnostic for when migration will be important. We will find conditions that the physical constants must satisfy so that all planets do not migrate into their host stars at an early stage of evolution, our observed values being intermediate such that this scenario only afflicts some percentage of planetary systems with unlucky characteristics. It should be noted that some amount of migration appears to have occurred even in the outer solar system [87]. However, we evidently ended up with at least one habitable location. In our particular instance, the particular positions of the giant planets ultimately halted migration, preventing the obliteration of the inner solar system [88].
Migration occurs when a planet’s influence on the surrounding disk produces a net torque on the planet itself, usually driving it inward. As such, migration halts after the time t disk when the disk has cleared out. This may be compared to the migration timescale
t mig a a ˙ L Γ Ω a 2 M Γ ,
where M is the planet’s mass, L its angular momentum and Γ the torque it experiences. A heuristic condition for when migration will be significant was found in [89] as t mig 10 t disk . There are various contributions to the torque, and which is dominant organizes migration into several different types, which depend on the circumstances of the case at hand. These are classified into type I, which arises when a trailing overdensity of dust behind a planet (and preceding in front) exerts a torque, and type II, in which the planet is capable of opening up a gap in the disk, resulting in a torque imbalance from the absence of material (for a recent review see [86]). The latter is more relevant for larger planets capable of significantly altering the disk structure, and the former more relevant for smaller planets, which are not. They will each be considered in turn, resulting in conditions on the fundamental constants that must be satisfied in order for a system with given characteristics to retain its initially habitable planets.
For type I migration, the timescale is derived in [90]:
t I = h 2 M 2 Σ Ω a 4 M ,
where h is the disk height, set by the sound speed. Then, when using the expressions from the Appendix A and Equation (10), and specifying to terrestrial planets situated in the temperate zone, we find
t I t disk λ 43 / 48 β 9 / 16 γ 3 / 8 κ α 3 / 2 .
If this quantity becomes too large, then Earthlike planets in all systems will migrate into their stars, and the universe would be uninhabitable in the traditional sense. Note the dependence λ 0.9 , indicating that this type of migration is more important for low mass stars.
This estimate would also be relevant for the production mechanism for Trappist-1 type planets proposed in [91], whereby earthlike planets form behind the snow line before migrating inwards. Such an unconventional pathway is needed to explain this system [92], which has multiple Earth sized planets orbiting the temperate zone of a 0.08 solar mass star with planets near mean motion resonances [93] and possessing icy compositions [84]. This will not be undertaken here, however.
Even if we restrict our attention to parameter values where Earthlike planets do not undergo substantial migration, we may still run into trouble if Jupiter-like planets routinely barrel through their planetary systems. As this is governed by the physics of type II migration, the conditions that must be satisfied for (the majority of) these planets to stay put are somewhat different. Here, we restrict our attention to the case where a full gap is opened in the disk, and where the mass of the planet is substantially smaller than the mass contained in the disk. In this case, the timescale of this migration process is given by [86]:
t II = 2 3 a 2 ν ,
where ν is the viscosity. Strictly speaking, these quantities are supposed to be evaluated not at the position of the planet, but rather the place of maximum angular momentum deposition, which can be approximated as a + 2.5 R Hill . For our purposes, however, this correction is negligible for the scaling arguments. Then, the ratio of this timescale to the disk lifetime is
t II t disk λ 15 / 16 γ 5 / 8 α disk α 4 β 13 / 16 .
Here, α disk is the standard parameterization of disk viscosity, discussed further in the appendix. Since this is inversely proportional to disk mass, type II migration is more important for smaller disks, opposite to the type I case. Additionally, since only the scaling with β is flipped from before, the conjunction of these two migration scenarios is capable of putting an upper bound on α and a lower bound on γ , if the other quantity is fixed. The thresholds for both types of migration are displayed in Figure 5. However, the absolute normalization of each of these timescales is uncertain, so we do not derive how the probabilities are altered due to these effects.

4. Discussion: Comparing 480 Hypotheses

Having spent the last two sections detailing the physics behind a multitude of processes that may influence the habitability of a system, we now synthesize these into an estimate of the probability of observing our measured parameter values, for each combination of individual habitability hypotheses. To summarize our results so far, we have firstly included several conditions necessary for planet formation, namely that the majority of galaxies should be larger than the minimal retentive mass that massive stars should have shorter lifetimes than the star formation timescale, and that the minimum metallicity needed to form planets should be smaller than the asymptotic value. Of these three, the third was most constraining. These are also presumably not optional, as opposed to the rest of the criteria that were considered: these were the absence of hot Jupiters, the production of terrestrial planets both through giant impact and isolation formation pathways, and the fraction of planets that end up in the temperate zone. Whether these are necessary for life are still intensely debated, and so each provides a prime opportunity to determine its compatibility with the multiverse hypothesis. Because they are all independent criteria, taken in conjunction they lead to a total of 2 × 3 × 2 = 12 possibilities. This does not include the 16 different choices we made in terms of planet formation, as well as the potentially continuous parameter signifying the slope of the power law for smaller planets. We display the probabilities of observing our measured values in Table 4 for the criteria mentioned in this paper. Note that, though the spread in probabilities is around 2–3 for each, all choices are well within an acceptable range to explain our observations within the multiverse context.
When combined with the 40 additional combinations from [1] (including combinations of the tidal locking criterion, which posits that only planets that are not tidally locked are habitable, the biological timescale criterion, positing that only stars that last several billion years are habitable, the convective criterion, which states that only stars which are not purely convective are habitable, the photosynthesis and yellow criteria, which state that photosynthesis is necessary for complex life, defined with an optimistic and pessimistic wavelength range, respectively, and the entropy criterion, where habitability is proportional to the total amount of entropy processed by a system), there are a total of 480 separate habitability criteria that may be considered. Of these, only 43 % of them give rise to probabilities of observing all three constants we consider of greater than 1 % . The full suite of criteria is displayed in Table 5 at the end of this manuscript. For brevity, we omit the convective criteria of [1] because it only ever marginally changes the numerical values, and in all instances its inclusion does not affect the viability of the combination of other hypotheses one way or the other. Of the 190 habitability criteria which give probabilities of over 10 % , all make use of the entropy condition. A further 16 which do not include the entropy condition have probabilities greater than 1 % —all of them benefit from an interplay between the yellow, tidal locking and biological timescale criteria, which place both upper and lower bounds on the types of allowed stars. The rest of the habitability criteria can safely be regarded as incompatible with the multiverse hypothesis.
The inclusion of multiple hypotheses leads to nonlinear effects, as the interplay between the distribution of purported observers and the anthropic boundaries alter the overall probabilities in sometimes surprising ways. That being said, none of the criteria have that drastic of an effect on the probabilities, especially when including the entropy condition.
Some criteria, namely the terrestrial and temperate conditions, introduce lower bounds to some combination of α and β . In fact, lower bounds on these quantities are somewhat hard to come by in the anthropics literature, though it has always been clear that they should exist, as a world with massless electrons or no electromagnetism would certainly be very different from our own. The bounds we find are stronger than those that exist in the literature.
We also introduce a new measure of our universe’s fitness, which we term the luxuriance. This is defined as the expected number of observers in our universe divided by the average number of observers per universe, restricting to universes that do have observers:
L = P ( α obs , β obs , γ obs ) d α θ ( P ( α , β , γ ) ) d α P ( α , β , γ ) .
Here, the integration is over all three constants, and θ ( x ) is the Heaviside step function: θ ( 0 ) = 0 , θ ( x ) = 1 for x > 0 . The rationale for including this is that, for most habitability criteria, the vast majority of universes will be sterile, obfuscating comparisons between different criteria. Restricting to universes that only contain life gives a better feeling for how good our universe is at satisfying the chosen criteria. If our universe is better than typical at making life, then this quantity will be greater than 1. While this is not actually the guiding principle for evaluating whether our observations are consistent with the multiverse, it is a somewhat interesting quantity to consider. It gives some indication of how strongly observers may cluster within the multiverse- and the strong dependence of the properties we discuss on physical constants leads us to expect that they will, so that the majority of observers do find themselves in overly productive universes. The luxuriance ranges by two orders of magnitude for the different possibilities we consider, but the maximum is L = 78.4 for the condition that includes the yellow, tidal locking, biological timescale, hot Jupiter, isolation planet production mechanism, temperate, and entropy criteria.
In addition, recall that certain choices for the physical processes involved in determining the structure of planets greatly affected some of the probabilities: the absence of a dependence on the mass of a planet with stellar mass, and too large a slope for the isolation mass, are both disfavored from the multiverse perspective. These scenarios are not otherwise excluded, but it would count as strong evidence against the multiverse if either of these were verified to be the case.
The one conclusion that should be drawn from this study is that planets are a rather generic feature, and not especially atypical for our particular values of the fundamental constants. Metal buildup is generic, massive stars burn out quickly, and disks tend to clump faster than they dissipate, so the inclusion of these criteria barely influenced the numerical values of the probabilities at all, apart from adding some mild anthropic boundaries. Furthermore, even specifying the characteristics of the planets that are formed, or considering different scenarios of their formation, did not alter the probabilities by more than a factor of few for the most part. This points to a reassuring robustness of our predictions that they are not as highly sensitive to the vagaries of incompletely known planet formation processes as one may have feared.
While it is simple to imagine a universe where planets are almost never made (and indeed before the current plethora of exoplanets was discovered many wondered if our universe was of this character), the parameter values needed to realize this possibility are quite extreme. Additionally, there are several features of planets in our universe that are tantalizingly coincidental, seeming to beckon for an anthropic explanation: the tendency to produce Earth mass planets, the similarity of the interplanetary spacing around sunlike stars and the width of the habitable zone, and the small but nonzero fraction of stars with hot Jupiters. Nevertheless, this reasoning was not borne out when incorporating these criteria into our analysis. The multiverse hypothesis is consistent with the expectation that life may only arise on temperate, terrestrial planets in systems without hot Jupiters, but it is essentially as compatible with the complete converse. The existence of suitable planetary environments thus does not seem to be the most important factor in determining which of the potential universes we find ourselves situated in. The other factors of the Drake equation which we will explore in subsequent works [94,95] will uncover many additional predictions for the requirements of life.

5. Conclusions

We have demonstrated that there are plenty of habitability conditions that are completely incompatible with the multiverse: what this illustrates is that, if any of the ones we have uncovered so far are shown to be the correct condition for the emergence of intelligent life, then we will be able to conclude to a very high degree of confidence (up to 5.2 σ ) that the multiverse must be wrong. It should be stressed that there is a great deal more that these conditions omit: nothing at all is said about how habitability is affected by things like planetary eccentricity, elemental composition, water abundance, or a host of other potentially paramount aspects of a planetary system [96,97]. The de facto stance on all omissions is that they have no bearing on habitability, and it will only be through future work, including all possibly relevant aspects that a fully coherent list of predictions may be assembled. Further still, placing a priority on the relative availability of each type of universe based on reasonably generic arguments, the precise probabilities, and the conclusions that follow will tremendously benefit from a way of being able to derive this prior with absolute surety. Since only a single one of the myriad habitability criteria is ultimately true, and since we will eventually be able to determine which one that is once we have a large enough sample of life-bearing planets, this demonstrates that the multiverse is capable of generating strong experimentally testable predictions that are capable of being verified or falsified on a reasonable timescale, the hallmark of a sensible scientific theory.

Funding

This research received no external funding.

Acknowledgments

I would like to thank Cullen Blake, Diana Dragomir, Scott Kenyon, Jabran Zahid, and Li Zeng for useful discussions.

Conflicts of Interest

The author declares no conflict of interest.

Appendix A. Planetary Parameters

In this appendix, we collect results on how various quantities relevant to our estimates in this work depend on orbital parameters of the stellar system, as well as fundamental quantities. To begin, we display the typical molecular binding energy:
T mol = α m mol r mol 3 = 0.037 α 2 m e 3 / 2 m p 1 / 2 .
This also defines the temperature required for liquid water.
Planets: The mass of a terrestrial planet, based on the criteria that carbon dioxide but not helium is gravitationally bound to the surface at these temperatures, is
M terr = 9.2 202.2 α 3 / 2 m e 3 / 4 M p l 3 m p 11 / 4 .
Next, we display the range of habitable orbits from the host star, using estimates from [1,98]:
a temp = 7.2 12.3 λ 7 / 4 m p 1 / 2 M p l 1 / 2 α 5 m e 2 ,
which is based off of the temperature of a perfect blackbody located at that position,
T ( a ) = T R 2 a .
The speed of a circularly orbiting planet at a given location is
v = G M a = Ω a
and Ω is the angular frequency. This also defines the orbital period as t Kepler = 2 π / Ω .
The region of influence of a planet of mass M and orbiting a star at semimajor axis a is known as the Hill sphere, and has the radius
R Hill = M 3 M 1 / 3 a .
Galaxies: the average galactic density is set by the galaxy at time of virialization, with a subsequent era of contraction due to equilibration [11,99]:
ρ gal = 5.6 × 10 7 κ 3 m p 4
corresponding to 10 23 g/cm 3 . Here, κ = Q ( η ω ) 4 / 3 = 10 16 is a composite of the primordial amplitude of perturbations Q = 1.8 × 10 5 , the baryon to photon ratio η = 6 × 10 10 , and the total matter to baryon ratio ω = 6.4 . This defines the typical freefall timescale as
t ff = 1 G ρ gal = 6.7 × 10 4 M p l κ 3 / 2 m p 2 ,
which equates to 3 × 10 7 yr. In addition, the typical temperature of the interstellar gas is relevant, which is set by the threshold for H 2 cooling, at roughly 10 4 K:
T H 2 = 0.032 α 2 m e .
This also sets the velocity dispersion in the galaxy as v 5 T H 2 / m p 20 km/s.
Disks: We take the disk mass to be set by
M disk = 0.01 M .
The size of the disk is set by the conservation of angular momentum of the initial collapsing cloud. If it has a typical angular momentum of L 0.01 M v r , with v given by the dispersion velocity of molecular clouds, we find
r disk M ρ gal 1 / 3 = 3.6 × 10 6 λ 1 / 3 M p l κ m p 2
set to be 100 AU for sunlike stars.
For infalling dust, the accretion rate M ˙ is just given in terms of the typical free fall timescale as [100]
M ˙ = c s 3 G N = 7.2 × 10 4 λ 2 α 3 m e 9 / 4 M p l 2 m p 9 / 4 .
In evaluating this quantity, the molecular energy has been used, reflecting that star formation exclusively forms within molecular clouds, where increased levels of cooling facilitate collapse. A disk whose dominant form of heating is given by accretion will have a temperature set by [55]
T accr 4 = 3 8 π G M M ˙ a 3 .
The speed of sound in the disk is given by c s T ( a ) / m p . This sets the height of the disk as the typical sound dispersion during one oscillation period. Using the accretion temperature,
h c s Ω = . 34 α 3 / 8 m e 9 / 32 a 9 / 8 λ 1 / 8 m p 1 / 32 M p l 1 / 8 .
This is normalized to 0.08 AU for earthlike orbits around sunlike stars. The viscosity of the disk is usually parameterized as
ν = α disk c s 2 Ω .
Typical values for the coefficient of proportionality are α disk 10 3 10 2 . In equilibrium, the surface density of the disk can be found to be given by
Σ = M ˙ 3 π ν .
The mass loss is usually taken to be independent of position [60], so that the surface density profile is dictated by the temperature. Early models took Σ a 3 / 2 , but this is now considered unlikely for equilibrium disks. The now standard dependence is Σ 1 / a [51], which is both observed in simulations and understood from a theoretical perspective [54]. We will take our profile to be
Σ ( a ) = M disk 2 r disk a θ ( r disk a ) .
The snow line is the distance beyond which the disk is cool enough for water to condense into a solid phase, which is found by setting the temperature equal to the vibrational energy of water molecules,
a snow = 30.4 λ M p l α 5 / 3 m e 5 / 4 m p 3 / 4 .
The size a planetesimal attains is given by the isolation mass, M iso 10 × 2 π Σ a R Hill ,60]. The rationale behind this quantity is that, once a planet is above this mass, it will have already depleted all the material within its Hill radius, and so it will no longer have a supply to continue its growth. Note that the Hill radius is itself a function of the planet mass, so that solving this equation for the mass yields the expression:
M iso ( 2 π Σ a 2 ) 3 / 2 M 1 / 2 = 1.2 × 10 6 κ 3 / 2 λ 1 / 2 m p M p l 3 / 2 a 3 / 2 .
The coefficient is set by noting that the isolation mass is about Mars sized for the solar system.

References

  1. Sandora, M. Multiverse Predictions for Habitability I: The Number of Stars and Their Properties. arXiv 2019, arXiv:1901.04614. [Google Scholar] [CrossRef]
  2. Vilenkin, A. Predictions from Quantum Cosmology. Phys. Rev. Lett. 1995, 74, 846–849. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Frank, A.; Sullivan, W.T., III. A New Empirical Constraint on the Prevalence of Technological Species in the Universe. Astrobiology 2016, 16, 359–362. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Maschberger, T. On the function describing the stellar initial mass function. Mon. Not. R. Astron. Soc. 2012, 429, 1725–1733. [Google Scholar] [CrossRef] [Green Version]
  5. Johnson, J.L.; Li, H. The first planets: The critical metallicity for planet formation. Astrophys. J. 2012, 751, 81. [Google Scholar] [CrossRef]
  6. Petigura, E.A.; Marcy, G.W.; Howard, A.W. A plateau in the planet population below twice the size of Earth. Astrophys. J. 2013, 770, 69. [Google Scholar] [CrossRef]
  7. Owen, J.E.; Wu, Y. The evaporation valley in the Kepler planets. Astrophys. J. 2017, 847, 29. [Google Scholar] [CrossRef]
  8. Ginzburg, S.; Schlichting, H.E.; Sari, R. Core-powered mass loss sculpts the radius distribution of small exoplanets. arXiv 2017, arXiv:1708.01621. [Google Scholar]
  9. Zeng, L.; Jacobsen, S.B.; Sasselov, D.D.; Vanderburg, A. Survival function analysis of planet size distribution with Gaia Data Release 2 updates. Mon. Not. R. Astron. Soc. 2018, 479, 5567–5576. [Google Scholar] [CrossRef]
  10. Adams, F.C. Constraints on Alternate Universes: Stars and habitable planets with different fundamental constants. J. Cosmol. Astropart. Phys. 2016, 2016, 042. [Google Scholar] [CrossRef]
  11. Adams, F.C.; Coppess, K.R.; Bloch, A.M. Planets in other universes: Habitability constraints on density fluctuations and galactic structure. J. Cosmol. Astropart. Phys. 2015, 2015, 030. [Google Scholar] [CrossRef]
  12. Weinberg, S. Anthropic bound on the cosmological constant. Phys. Rev. Lett. 1987, 59, 2607–2610. [Google Scholar] [CrossRef] [PubMed]
  13. Thielemann, F.K.; Nomoto, K.; Hashimoto, M.A. Core-collapse supernovae and their ejecta. Astrophys. J. 1996, 460, 408. [Google Scholar] [CrossRef]
  14. Rees, M.J.; Ostriker, J. Cooling, dynamics and fragmentation of massive gas clouds: Clues to the masses and radii of galaxies and clusters. Mon. Not. R. Astron. Soc. 1977, 179, 541–559. [Google Scholar] [CrossRef]
  15. Padmanabhan, T. Theoretical Astrophysics: Volume 2, Stars and Stellar Systems; Cambridge University Press: Cambridge, UK, 2001. [Google Scholar]
  16. Tremonti, C.A.; Heckman, T.M.; Kauffmann, G.; Brinchmann, J.; Charlot, S.; White, S.D.; Seibert, M.; Peng, E.W.; Schlegel, D.J.; Uomoto, A.; et al. The origin of the mass-metallicity relation: Insights from 53,000 star-forming galaxies in the sloan digital sky survey. Astrophys. J. 2004, 613, 898. [Google Scholar] [CrossRef]
  17. Press, W.H.; Schechter, P. Formation of galaxies and clusters of galaxies by self-similar gravitational condensation. Astrophys. J. 1974, 187, 425–438. [Google Scholar] [CrossRef]
  18. Dayal, P.; Ward, M.; Cockell, C. The habitability of the Universe through 13 billion years of cosmic time. arXiv 2016, arXiv:1606.09224. [Google Scholar]
  19. Woosley, S.E.; Heger, A.; Weaver, T.A. The evolution and explosion of massive stars. Rev. Mod. Phys. 2002, 74, 1015–1071. [Google Scholar] [CrossRef]
  20. Burrows, A.S.; Ostriker, J.P. Astronomical reach of fundamental physics. Proc. Natl. Acad. Sci. USA 2014, 111, 2409–2416. [Google Scholar] [CrossRef] [Green Version]
  21. Shakura, N.I.; Sunyaev, R.A. Black holes in binary systems. Observational appearance. Astron. Astrophys. 1973, 24, 337–355. [Google Scholar]
  22. Alexander, R.; Pascucci, I.; Andrews, S.; Armitage, P.; Cieza, L. The dispersal of protoplanetary disks. arXiv 2013, arXiv:1311.1819. [Google Scholar]
  23. Apai, D.; Lauretta, D.S. Protoplanetary Dust: Astrophysical and Cosmochemical Perspectives; Cambridge University Press: Cambridge, UK, 2010; Volume 12. [Google Scholar]
  24. Ercolano, B.; Clarke, C. Metallicity, planet formation and disc lifetimes. Mon. Not. R. Astron. Soc. 2010, 402, 2735–2743. [Google Scholar] [CrossRef] [Green Version]
  25. Zahid, H.J.; Dima, G.I.; Kudritzki, R.P.; Kewley, L.J.; Geller, M.J.; Hwang, H.S.; Silverman, J.D.; Kashino, D. The universal relation of galactic chemical evolution: The origin of the mass-metallicity relation. Astrophys. J. 2014, 791, 130. [Google Scholar] [CrossRef]
  26. Schellekens, A.N. Life at the Interface of Particle Physics and String Theory. Rev. Mod. Phys. 2013, 85, 1491–1540. [Google Scholar] [CrossRef]
  27. Tegmark, M.; Rees, M.J. Why Is the Cosmic Microwave Background Fluctuation Level 10-5? Astrophys. J. 1998, 499, 526–532. [Google Scholar] [CrossRef]
  28. Fischer, D.A.; Valenti, J. The Planet-Metallicity Correlation. Astrophys. J. 2005, 622, 1102–1117. [Google Scholar] [CrossRef]
  29. Batygin, K.; Bodenheimer, P.H.; Laughlin, G.P. In situ formation and dynamical evolution of hot Jupiter systems. Astrophys. J. 2016, 829, 114. [Google Scholar] [CrossRef]
  30. Dawson, R.I.; Johnson, J.A. Origins of Hot Jupiters. arXiv 2018, arXiv:1801.06117. [Google Scholar] [CrossRef]
  31. Raymond, S.N.; Mandell, A.M.; Sigurdsson, S. Exotic Earths: Forming Habitable Worlds with Giant Planet Migration. Science 2006, 313, 1413–1416. [Google Scholar] [CrossRef] [Green Version]
  32. Smallwood, J.L.; Martin, R.G.; Lepp, S.; Livio, M. Asteroid impacts on terrestrial planets: The effects of super-Earths and the role of the ν 6 resonance. Mon. Not. R. Astron. Soc. 2017, 473, 295–305. [Google Scholar] [CrossRef]
  33. Buchhave, L.A.; Bitsch, B.; Johansen, A.; Latham, D.W.; Bizzarro, M.; Bieryla, A.; Kipping, D.M. Jupiter Analogues Orbit Stars with an Average Metallicity Close to that of the Sun. arXiv 2018, arXiv:1802.06794. [Google Scholar]
  34. Ndugu, N.; Bitsch, B.; Jurua, E. Planet population synthesis driven by pebble accretion in cluster environments. Mon. Not. R. Astron. Soc. 2017, 474, 886–897. [Google Scholar] [CrossRef]
  35. Fabrycky, D.; Tremaine, S. Shrinking binary and planetary orbits by Kozai cycles with tidal friction. Astrophys. J. 2007, 669, 1298. [Google Scholar] [CrossRef]
  36. Becker, J.C.; Vanderburg, A.; Adams, F.C.; Rappaport, S.A.; Schwengeler, H.M. WASP-47: A hot Jupiter system with two additional planets discovered by K2. Astrophys. Lett. 2015, 812, L18. [Google Scholar] [CrossRef]
  37. Chatterjee, S.; Ford, E.B.; Matsumura, S.; Rasio, F.A. Dynamical outcomes of planet-planet scattering. Astrophys. J. 2008, 686, 580. [Google Scholar] [CrossRef]
  38. Spalding, C.; Batygin, K. A Secular Resonant Origin for the Loneliness of Hot Jupiters. Astron. J. 2017, 154, 93. [Google Scholar] [CrossRef] [Green Version]
  39. Johansen, A.; Lambrechts, M. Forming Planets via Pebble Accretion. Annu. Rev. Earth Planet. Sci. 2017, 45, 359–387. [Google Scholar] [CrossRef]
  40. Johnson, J.A.; Aller, K.M.; Howard, A.W.; Crepp, J.R. Giant planet occurrence in the stellar mass-metallicity plane. Publ. Astron. Soc. Pac. 2010, 122, 905. [Google Scholar] [CrossRef]
  41. Taubner, R.S.; Pappenreiter, P.; Zwicker, J.; Smrzka, D.; Pruckner, C.; Kolar, P.; Bernacchi, S.; Seifert, A.H.; Krajete, A.; Bach, W.; et al. Biological methane production under putative Enceladus-like conditions. Nat. Commun. 2018, 9, 748. [Google Scholar] [CrossRef] [Green Version]
  42. Bains, W. Many chemistries could be used to build living systems. Astrobiology 2004, 4, 137–167. [Google Scholar] [CrossRef]
  43. Schulze-Makuch, D.; Irwin, L.N. The prospect of alien life in exotic forms on other worlds. Naturwissenschaften 2006, 93, 155–172. [Google Scholar] [CrossRef] [PubMed]
  44. Rogers, L.A. Most 1.6 Earth-radius planets are not rocky. Astrophys. J. 2015, 801, 41. [Google Scholar] [CrossRef]
  45. Ward, P.D.; Brownlee, D. Rare Earth: Why Complex Life Is Uncommon in the Universe; Copernicus Books: New York, NY, USA, 2003. [Google Scholar]
  46. Owen, J.E.; Lai, D. Photoevaporation and high-eccentricity migration created the sub-Jovian desert. Mon. Not. R. Astron. Soc. 2018, 479, 5012–5021. [Google Scholar] [CrossRef] [Green Version]
  47. Fulton, B.J.; Petigura, E.A.; Howard, A.W.; Isaacson, H.; Marcy, G.W.; Cargile, P.A.; Hebb, L.; Weiss, L.M.; Johnson, J.A.; Morton, T.D.; et al. The California-Kepler survey. III. A gap in the radius distribution of small planets. Astron. J. 2017, 154, 109. [Google Scholar] [CrossRef]
  48. Fressin, F.; Torres, G.; Charbonneau, D.; Bryson, S.T.; Christiansen, J.; Dressing, C.D.; Jenkins, J.M.; Walkowicz, L.M.; Batalha, N.M. The false positive rate of Kepler and the occurrence of planets. Astrophys. J. 2013, 766, 81. [Google Scholar] [CrossRef]
  49. Raymond, S.N.; Boulet, T.; Izidoro, A.; Esteves, L.; Bitsch, B. Migration-driven diversity of super-Earth compositions. Mon. Not. R. Astron. Soc. Lett. 2018, 479, L81–L85. [Google Scholar] [CrossRef]
  50. Kokubo, E.; Kominami, J.; Ida, S. Formation of terrestrial planets from protoplanets. I. Statistics of basic dynamical properties. Astrophys. J. 2006, 642, 1131. [Google Scholar] [CrossRef]
  51. Williams, J.P.; Cieza, L.A. Protoplanetary disks and their evolution. Annu. Rev. Astron. Astrophys. 2011, 49, 67–117. [Google Scholar] [CrossRef]
  52. Pascucci, I.; Testi, L.; Herczeg, G.; Long, F.; Manara, C.; Hendler, N.; Mulders, G.; Krijt, S.; Ciesla, F.; Henning, T.; et al. A steeper than linear disk mass–stellar mass scaling relation. Astrophys. J. 2016, 831, 125. [Google Scholar] [CrossRef]
  53. Bate, M.R. On the diversity and statistical properties of protostellar discs. Mon. Not. R. Astron. Soc. 2018, 475, 5618–5658. [Google Scholar] [CrossRef] [Green Version]
  54. Morbidelli, A.; Lambrechts, M.; Jacobson, S.; Bitsch, B. The great dichotomy of the Solar System: Small terrestrial embryos and massive giant planet cores. Icarus 2015, 258, 418–429. [Google Scholar] [CrossRef] [Green Version]
  55. Kennedy, G.M.; Kenyon, S.J. Planet formation around stars of various masses: The snow line and the frequency of giant planets. Astrophys. J. 2008, 673, 502. [Google Scholar] [CrossRef]
  56. Kennedy, G.M.; Kenyon, S.J. Planet formation around stars of various masses: Hot Super-Earths. Astrophys. J. 2008, 682, 1264. [Google Scholar] [CrossRef]
  57. Ida, S.; Guillot, T.; Morbidelli, A. The radial dependence of pebble accretion rates: A source of diversity in planetary systems-I. Analytical formulation. Astron. Astrophys. 2016, 591, A72. [Google Scholar] [CrossRef]
  58. Adams, F.C.; Hollenbach, D.; Laughlin, G.; Gorti, U. Photoevaporation of circumstellar disks due to external far-ultraviolet radiation in stellar aggregates. Astrophys. J. 2004, 611, 360. [Google Scholar] [CrossRef]
  59. Morbidelli, A.; Lunine, J.I.; O’Brien, D.P.; Raymond, S.N.; Walsh, K.J. Building terrestrial planets. Annu. Rev. Earth Planet. Sci. 2012, 40, 251–275. [Google Scholar] [CrossRef]
  60. Youdin, A.N.; Kenyon, S.J. From disks to planets. In Planets, Stars and Stellar Systems; Springer: New York, NY, USA, 2013; pp. 1–62. [Google Scholar]
  61. Izidoro, A.; Raymond, S.N. Formation of Terrestrial Planets. In Handbook of Exoplanets; Springer: New York, NY, USA, 2018; pp. 1–59. [Google Scholar] [Green Version]
  62. Schlichting, H.E. Formation of close in super-Earths and mini-Neptunes: Required disk masses and their implications. Astrophys. J. Lett. 2014, 795, L15. [Google Scholar] [CrossRef]
  63. Sinukoff, E.; Fulton, B.; Scuderi, L.; Gaidos, E. Below One Earth: The Detection, Formation, and Properties of Subterrestrial Worlds. Space Sci. Rev. 2013, 180, 71–99. [Google Scholar] [CrossRef] [Green Version]
  64. Wu, Z.N.; Li, J.; Bai, C.Y. Scaling relations of lognormal type growth process with an extremal principle of entropy. Entropy 2017, 19, 56. [Google Scholar] [CrossRef]
  65. Cumming, A.; Butler, R.P.; Marcy, G.W.; Vogt, S.S.; Wright, J.T.; Fischer, D.A. The Keck planet search: Detectability and the minimum mass and orbital period distribution of extrasolar planets. Publ. Astron. Soc. Pac. 2008, 120, 531. [Google Scholar] [CrossRef]
  66. Zeng, L.; Jacobsen, S.B.; Sasselov, D.D.; Vanderburg, A. Survival Function Analysis of Planet Orbit Distribution and Occurrence Rate Estimate. arXiv 2018, arXiv:1801.03994. [Google Scholar]
  67. Simon, J.B.; Armitage, P.J.; Li, R.; Youdin, A.N. The mass and size distribution of planetesimals formed by the streaming instability. I. The role of self-gravity. Astrophys. J. 2016, 822, 55. [Google Scholar] [CrossRef]
  68. Mordasini, C. Planetary population synthesis. In Handbook of Exoplanets; Springer: New York, NY, USA, 2018; pp. 1–50. [Google Scholar]
  69. Raymond, S.N.; Scalo, J.; Meadows, V.S. A decreased probability of habitable planet formation around low-mass stars. Astrophys. J. 2007, 669, 606. [Google Scholar] [CrossRef]
  70. Tasker, E.; Tan, J.; Heng, K.; Kane, S.; Spiegel, D.; Brasser, R.; Casey, A.; Desch, S.; Dorn, C.; Hernlund, J.; et al. The language of exoplanet ranking metrics needs to change. Nat. Astron. 2017, 1, 0042. [Google Scholar] [CrossRef] [Green Version]
  71. Kopparapu, R.K.; Ramirez, R.; Kasting, J.F.; Eymet, V.; Robinson, T.D.; Mahadevan, S.; Terrien, R.C.; Domagal-Goldman, S.; Meadows, V.; Deshpande, R. Habitable zones around main-sequence stars: New estimates. Astrophys. J. 2013, 765, 131. [Google Scholar] [CrossRef]
  72. Yang, J.; Boué, G.; Fabrycky, D.C.; Abbot, D.S. Strong dependence of the inner edge of the habitable zone on planetary rotation rate. Astrophys. J. Lett. 2014, 787, L2. [Google Scholar] [CrossRef]
  73. Kasting, J.F.; Whitmire, D.P.; Reynolds, R.T. Habitable zones around main sequence stars. Icarus 1993, 101, 108–128. [Google Scholar] [CrossRef] [PubMed]
  74. Leconte, J.; Forget, F.; Charnay, B.; Wordsworth, R.; Pottier, A. Increased insolation threshold for runaway greenhouse processes on Earth-like planets. Nature 2013, 504, 268. [Google Scholar] [CrossRef]
  75. Walker, J.C.; Hays, P.; Kasting, J.F. A negative feedback mechanism for the long-term stabilization of Earth’s surface temperature. J. Geophys. Res. Ocean. 1981, 86, 9776–9782. [Google Scholar] [CrossRef]
  76. Rushby, A.J.; Johnson, M.; Mills, B.J.; Watson, A.J.; Claire, M.W. Long-Term Planetary Habitability and the Carbonate-Silicate Cycle. Astrobiology 2018, 18, 469–480. [Google Scholar] [CrossRef]
  77. Barnes, R.; Quinn, T. The (in) stability of planetary systems. Astrophys. J. 2004, 611, 494. [Google Scholar] [CrossRef]
  78. Dawson, R.I. Tightly Packed Planetary Systems. In Handbook of Exoplanets; Springer: New York, NY, USA, 2017; pp. 1–18. [Google Scholar]
  79. Holman, M.J.; Wisdom, J. Dynamical stability in the outer solar system and the delivery of short period comets. Astron. J. 1993, 105, 1987–1999. [Google Scholar] [CrossRef]
  80. Fang, J.; Margot, J.L. Are planetary systems filled to capacity? A study based on Kepler results. Astrophys. J. 2013, 767, 115. [Google Scholar] [CrossRef]
  81. Raymond, S.N.; Barnes, R.; Veras, D.; Armitage, P.J.; Gorelick, N.; Greenberg, R. Planet-planet scattering leads to tightly packed planetary systems. Astrophys. J. Lett. 2009, 696, L98. [Google Scholar] [CrossRef]
  82. Borucki, W.J.; Koch, D.; Basri, G.; Batalha, N.; Brown, T.; Caldwell, D.; Caldwell, J.; Christensen-Dalsgaard, J.; Cochran, W.D.; DeVore, E.; et al. Kepler planet-detection mission: Introduction and first results. Science 2010, 327, 977–980. [Google Scholar] [CrossRef] [PubMed]
  83. Snellgrove, M.; Papaloizou, J.; Nelson, R. On disc driven inward migration of resonantly coupled planets with application to the system around GJ876. Astron. Astrophys. 2001, 374, 1092–1099. [Google Scholar] [CrossRef]
  84. Unterborn, C.T.; Desch, S.J.; Hinkel, N.R.; Lorenzo, A. Inward migration of the TRAPPIST-1 planets as inferred from their water-rich compositions. Nat. Astron. 2018, 2, 297–302. [Google Scholar] [CrossRef] [Green Version]
  85. Trilling, D.E.; Lunine, J.I.; Benz, W. Orbital migration and the frequency of giant planet formation. Astron. Astrophys. 2002, 394, 241–251. [Google Scholar] [CrossRef] [Green Version]
  86. Baruteau, C.; Masset, F. Recent developments in planet migration theory. In Tides in Astronomy and Astrophysics; Springer: New York, NY, USA, 2013; pp. 201–253. [Google Scholar]
  87. Morbidelli, A.; Crida, A. The dynamics of Jupiter and Saturn in the gaseous protoplanetary disk. Icarus 2007, 191, 158–171. [Google Scholar] [CrossRef] [Green Version]
  88. Masset, F.; Snellgrove, M. Reversing type II migration: Resonance trapping of a lighter giant protoplanet. Mon. Not. R. Astron. Soc. 2001, 320, L55–L59. [Google Scholar] [CrossRef]
  89. Ida, S.; Lin, D.N. Toward a deterministic model of planetary formation. I. A desert in the mass and semimajor axis distributions of extrasolar planets. Astrophys. J. 2004, 604, 388. [Google Scholar] [CrossRef]
  90. Tanaka, H.; Takeuchi, T.; Ward, W.R. Three-dimensional interaction between a planet and an isothermal gaseous disk. I. Corotation and Lindblad torques and planet migration. Astrophys. J. 2002, 565, 1257. [Google Scholar] [CrossRef]
  91. Ormel, C.W.; Liu, B.; Schoonenberg, D. Formation of TRAPPIST-1 and other compact systems. Astron. Astrophys. 2017, 604, A1. [Google Scholar] [CrossRef] [Green Version]
  92. Gillon, M.; Triaud, A.H.; Demory, B.O.; Jehin, E.; Agol, E.; Deck, K.M.; Lederer, S.M.; De Wit, J.; Burdanov, A.; Ingalls, J.G.; et al. Seven temperate terrestrial planets around the nearby ultracool dwarf star TRAPPIST-1. Nature 2017, 542, 456. [Google Scholar] [CrossRef] [PubMed]
  93. Tamayo, D.; Rein, H.; Petrovich, C.; Murray, N. Convergent Migration Renders TRAPPIST-1 Long-lived. Astrophys. J. Lett. 2017, 840, L19. [Google Scholar] [CrossRef] [Green Version]
  94. Sandora, M. Multiverse Predictions for Habitability III: Fraction of Planets That Develop Life. arXiv 2019, arXiv:1903.06283. [Google Scholar]
  95. Sandora, M. Multiverse Predictions for Habitability IV: Fraction of Life that Develops Intelligence. arXiv 2019, arXiv:1904.11796. [Google Scholar]
  96. Schulze-Makuch, D.; Méndez, A.; Fairén, A.G.; Von Paris, P.; Turse, C.; Boyer, G.; Davila, A.F.; António, M.R.d.S.; Catling, D.; Irwin, L.N. A two-tiered approach to assessing the habitability of exoplanets. Astrobiology 2011, 11, 1041–1052. [Google Scholar] [CrossRef]
  97. Cockell, C.S.; Bush, T.; Bryce, C.; Direito, S.; Fox-Powell, M.; Harrison, J.; Lammer, H.; Landenmark, H.; Martin-Torres, J.; Nicholson, N.; et al. Habitability: A review. Astrobiology 2016, 16, 89–117. [Google Scholar] [CrossRef]
  98. Press, W.H.; Lightman, A.P. Dependence of macrophysical phenomena on the values of the fundamental constants. Philos. Trans. R. Soc. Lond. Ser. A 1983, 310, 323–334. [Google Scholar] [CrossRef]
  99. Tegmark, M.; Aguirre, A.; Rees, M.J.; Wilczek, F. Dimensionless constants, cosmology, and other dark matters. Phys. Rev. D 2006, 73, 023505. [Google Scholar] [CrossRef] [Green Version]
  100. Shu, F.H. Self-similar collapse of isothermal spheres and star formation. Astrophys. J. 1977, 214, 488–497. [Google Scholar] [CrossRef]
1
The code to compute all probabilities discussed in the text is made available at https://github.com/mccsandora/Multiverse-Habitability-Handler.
2
A third explanation was additionally given in [34] that the disk dispersal timescale increases with metallicity, allowing a longer period of accretion onto seed cores.
3
This is not entirely true: there is known to be some feedback between the irradiation of initial atmospheres from the host star that favors both small and large atmospheres [46]. While this leads to the interesting bimodal distribution of observed radii [47], this is driven by atmospheric size, and is certainly not enough to affect the terrestrial core.
4
This neglects other potential thresholds, such as the inner boundary set by the photosynthetic threshold of carbon dioxide abundance, which is remarkably close to the inner boundary discussed here [76].
Figure 1. The distribution of observers when taking metallicity effects into account, in the α γ , α β , and β γ subplanes. What is plotted is the logarithm of the probability of measuring any value of the coupling constants, with red being more probable than blue, and the black dot corresponding to our observed values. Included thresholds are an upper bound on α for hydrogen stability, as well as a lower bound for galactic cooling, and an upper bound on β for proton–proton fusion. The teal and brown curves correspond to the metallicity and supernova lifetime thresholds, respectively, and the supernova retention threshold is not relevant in this range.
Figure 1. The distribution of observers when taking metallicity effects into account, in the α γ , α β , and β γ subplanes. What is plotted is the logarithm of the probability of measuring any value of the coupling constants, with red being more probable than blue, and the black dot corresponding to our observed values. Included thresholds are an upper bound on α for hydrogen stability, as well as a lower bound for galactic cooling, and an upper bound on β for proton–proton fusion. The teal and brown curves correspond to the metallicity and supernova lifetime thresholds, respectively, and the supernova retention threshold is not relevant in this range.
Universe 05 00157 g001
Figure 2. The distribution of observers if life can only arise on a terrestrial planet in the α γ , α β , and β γ subplanes, size being dictated by a giant impact phase. Here, we have used a log-normal distribution with σ = 1 / 6 , excluded λ dependence on the average planet mass, and assumed accretion dominated disks.
Figure 2. The distribution of observers if life can only arise on a terrestrial planet in the α γ , α β , and β γ subplanes, size being dictated by a giant impact phase. Here, we have used a log-normal distribution with σ = 1 / 6 , excluded λ dependence on the average planet mass, and assumed accretion dominated disks.
Universe 05 00157 g002
Figure 3. The distribution of terrestrial planetesimals that result from isolated accretion in the α γ , α β , and β γ subplanes, without the subsequent phase of giant impacts. The slope here is q = 1 / 3 . The teal line represents the region where the largest star capable of hosting terrestrial planets is smaller than the smallest possible star.
Figure 3. The distribution of terrestrial planetesimals that result from isolated accretion in the α γ , α β , and β γ subplanes, without the subsequent phase of giant impacts. The slope here is q = 1 / 3 . The teal line represents the region where the largest star capable of hosting terrestrial planets is smaller than the smallest possible star.
Universe 05 00157 g003
Figure 4. Distribution of observers if life may only exist within the temperate zone, in the α γ , α β , and β γ subplanes. The teal line is the boundary for which planetary disks are smaller than the temperate zone.
Figure 4. Distribution of observers if life may only exist within the temperate zone, in the α γ , α β , and β γ subplanes. The teal line is the boundary for which planetary disks are smaller than the temperate zone.
Universe 05 00157 g004
Figure 5. Display of when runaway migration takes hold, in the α γ , α β , and β γ subplanes. The teal line is for type I migration, the orange line is for type II, and the purple line is when the mass of the fastest migrating body is equal to the terrestrial mass. All curves have been normalized so the timescales in our universe are five times larger than the critical value.
Figure 5. Display of when runaway migration takes hold, in the α γ , α β , and β γ subplanes. The teal line is for type I migration, the orange line is for type II, and the purple line is when the mass of the fastest migrating body is equal to the terrestrial mass. All curves have been normalized so the timescales in our universe are five times larger than the critical value.
Universe 05 00157 g005
Table 1. The quantities computed in this work. Here, Q is the amplitude of perturbations, κ parameterizes the density of galaxies, λ parameterizes stellar size, GI stands for giant impact formation mechanism, and iso stands for isolation production mechanism.
Table 1. The quantities computed in this work. Here, Q is the amplitude of perturbations, κ parameterizes the density of galaxies, λ parameterizes stellar size, GI stands for giant impact formation mechanism, and iso stands for isolation production mechanism.
QuantityDescriptionExpression
f gal fraction of stars in galaxies that retain supernova ejecta erfc ( 4.1 Q 1 α 2 β 5 / 3 )
f 2 nd gen fraction of stars born after supernova enrichment exp ( 0.24 κ 3 / 2 α 7 / 4 β 1 / 8 γ 1 )
f Z fraction of stars with high enough Z for planets θ ( 1 0.038 λ 3 / 4 α 3 β 1 / 2 γ 1 / 2 )
f hj fraction of stars without hot Jupiters 1 2.3 × 10 8 κ 2 λ 2 α 13 / 8 β 3
n p average number of planets around a starGI: Equation (25),   iso: Equation (31)
f terr fraction of terrestrial planetsGI: Equation (27),   iso: Equation (30)
f temp fraction of temperate planets 0.0053 κ 1 / 2 λ 1.77 α 11 / 2 β 7 / 4 γ 5 / 8
Table 2. Display of the insensitivity of the probabilities of our observables to the choices made regarding planet formation. In addition to the options displayed, they may have been combined, but this would belabor the point.
Table 2. Display of the insensitivity of the probabilities of our observables to the choices made regarding planet formation. In addition to the options displayed, they may have been combined, but this would belabor the point.
Choices P ( α obs ) P ( β obs ) P ( γ obs )
standard0.2290.2600.409
Rayleigh distribution0.2290.2510.411
irradiation0.2550.3200.426
with no λ dependence0.3800.2540.007
shot noise0.2320.2600.306
Table 3. Display of the probabilities of our observables for different values of the power law slope.
Table 3. Display of the probabilities of our observables for different values of the power law slope.
Exponent P ( α obs ) P ( β obs ) P ( γ obs )
q = 1 / 3 0.2340.4190.436
q = 2 / 3 0.3020.2600.424
q = 1 0.3620.1570.333
q = 2 0.4690.0440.228
Table 4. Probabilities of different combinations of the habitability hypotheses discussed in the text. Here, yellow stands for the photosynthesis criterion, S the entropy criterion (both explained below Equation (1)), temp for the temperate zone, GI and iso for the terrestrial planet criterion with the giant impact and isolation production mechanisms, and HJ the hot Jupiter condition.
Table 4. Probabilities of different combinations of the habitability hypotheses discussed in the text. Here, yellow stands for the photosynthesis criterion, S the entropy criterion (both explained below Equation (1)), temp for the temperate zone, GI and iso for the terrestrial planet criterion with the giant impact and isolation production mechanisms, and HJ the hot Jupiter condition.
Criteria P ( α obs ) P ( β obs ) P ( γ obs ) L
yellow S0.1890.4370.3182.57
yellow temp S0.2410.3630.1511.54
yellow GI S0.2290.260.4096.72
yellow GI temp S0.3770.1960.426.65
yellow iso S0.2340.4190.4365.02
yellow iso temp S0.3130.480.2453.54
yellow HJ S0.1810.4280.3082.59
yellow HJ temp S0.2370.3590.1471.54
yellow HJ GI S0.230.2610.416.72
yellow HJ GI temp S0.3780.1970.4196.64
yellow HJ iso S0.2310.4240.4315.05
yellow HJ iso temp S0.3110.4830.2413.54
Table 5. Probabilities of various hypotheses, including those from [1] (continued on following pages). In addition to the abbreviations from Table 4, the shorthand is: photo: photosynthesis (optimistic), yellow: photosynthesis (conservative), TL: tidal locking, bio: biological timescale, S: entropy, temp: temperate zone, GI and iso: terrestrial planet with giant impact and isolation production mechanisms, resp. and HJ: hot Jupiter.
Table 5. Probabilities of various hypotheses, including those from [1] (continued on following pages). In addition to the abbreviations from Table 4, the shorthand is: photo: photosynthesis (optimistic), yellow: photosynthesis (conservative), TL: tidal locking, bio: biological timescale, S: entropy, temp: temperate zone, GI and iso: terrestrial planet with giant impact and isolation production mechanisms, resp. and HJ: hot Jupiter.
Criteria P ( α obs ) P ( β obs ) P ( γ obs ) L
number of stars0.3810.3558.06  × 10 7 0.000199
temp0.1750.05248.87  × 10 6 0.000923
GI0.4240.2811.16  × 10 6 0.00021
GI temp0.2270.2799.7  × 10 6 0.00269
iso0.4220.4938.69  × 10 7 0.000239
iso temp0.1980.1128.9  × 10 6 0.00158
HJ0.3660.3367.14  × 10 7 0.000199
HJ temp0.1750.05238.87  × 10 6 0.00092
HJ GI0.4250.2821.15  × 10 6 0.00021
HJ GI temp0.2270.2799.7  × 10 6 0.0027
HJ iso0.4190.4998.53  × 10 7 0.00024
HJ iso temp0.1980.1128.9  × 10 6 0.00158
bio0.2840.1062.32  × 10 5 0.00379
bio temp0.1450.01523.95  × 10 5 0.00282
bio GI0.2110.01610.0004050.0562
bio GI temp0.2820.01620.001530.324
bio iso0.2440.2165.28  × 10 5 0.00961
bio iso temp0.1670.03345.42  × 10 5 0.00661
bio HJ0.2840.1062.03  × 10 5 0.00378
bio HJ temp0.1450.01523.95  × 10 5 0.00281
bio HJ GI0.2110.01620.0004050.0565
bio HJ GI temp0.2820.01630.001530.325
bio HJ iso0.2440.2155.13  × 10 5 0.00958
bio HJ iso temp0.1670.03345.41  × 10 5 0.00661
TL0.330.4554.57  × 10 7 3.39  × 10 5
TL temp0.4550.2571.93  × 10 7 4.61  × 10 5
TL GI0.4210.2836.42 × 10 7 3.09 × 10 5
TL GI temp0.2370.298.16 × 10 7 0.000146
TL iso0.3980.4334.6 × 10 7 5.6 × 10 5
TL iso temp0.4540.3563.63 × 10 7 0.000107
TL HJ0.3060.4333.11 × 10 7 3.15 × 10 5
TL HJ temp0.4540.2561.8 × 10 7 4.52 × 10 5
TL HJ GI0.4220.2846.38 × 10 7 3.08 × 10 5
TL HJ GI temp0.2360.2918.12 × 10 7 0.000146
TL HJ iso0.3930.4394.39 × 10 7 5.53 × 10 5
TL HJ iso temp0.4550.3553.55 × 10 7 0.000106
TL bio0.0110.3657.19 × 10 5 0.00304
TL bio temp0.01540.2764.16 × 10 5 0.00788
TL bio GI0.05620.01840.0002770.0109
TL bio GI temp0.180.0230.0004290.0631
TL bio iso0.02370.4850.0001070.0087
TL bio iso temp0.04030.4079.85 × 10 5 0.0238
TL bio HJ0.009680.3624.84 × 10 5 0.00301
TL bio HJ temp0.01490.2753.89 × 10 5 0.00779
TL bio HJ GI0.05610.01850.0002760.0109
TL bio HJ GI temp0.180.02310.0004280.063
TL bio HJ iso0.02280.4880.0001010.00861
TL bio HJ iso temp0.03940.4069.61 × 10 5 0.0235
photo0.4390.1838.43 × 10 7 0.000127
photo S0.2410.3820.3818.95
photo temp0.4030.1213.11 × 10 7 7.34 × 10 5
photo temp S0.3380.2920.2077.07
photo GI0.0160.2968.09 × 10 6 0.00169
photo GI S0.3350.1420.3334.8
photo GI temp0.009420.2768.96 × 10 6 0.00513
photo GI temp S0.4970.1130.42445.7
photo iso0.4560.2671.75 × 10 6 0.000346
photo iso S0.3050.4480.49413.7
photo iso temp0.3290.186.74 × 10 7 0.000207
photo iso temp S0.4030.4580.32112.7
photo HJ0.4390.1838.22 × 10 7 0.000126
photo HJ S0.2370.3770.3769.0
photo HJ temp0.4030.1213.08 × 10 7 7.31 × 10 5
photo HJ temp S0.3370.2910.2057.08
photo HJ GI0.0160.2978.04 × 10 6 0.00168
photo HJ GI S0.3360.1420.33135.0
photo HJ GI temp0.009420.2768.92 × 10 6 0.00513
photo HJ GI temp S0.4980.1130.42545.9
photo HJ iso0.4560.2671.72 × 10 6 0.000342
photo HJ iso S0.3020.4520.49713.8
photo HJ iso temp0.3290.186.68 × 10 7 0.000206
photo HJ iso temp S0.4020.4550.31912.8
photo bio0.06310.1031.75 × 10 5 0.00197
photo bio S0.180.4090.437.89
photo bio temp0.1230.08495.35 × 10 6 0.000987
photo bio temp S0.2220.3350.2626.98
photo bio GI0.3530.01620.0009240.166
photo bio GI S0.3160.1310.30830.9
photo bio GI temp0.1540.01410.0006880.34
photo bio GI temp S0.4840.10.40840.3
photo bio iso0.1040.1885.14 × 10 5 0.00757
photo bio iso S0.2620.420.44811.6
photo bio iso temp0.1850.1521.51 × 10 5 0.00365
photo bio iso temp S0.3220.4890.38211.8
photo bio HJ0.06310.1031.71 × 10 5 0.00195
photo bio HJ S0.1750.4040.4267.94
photo bio HJ temp0.1230.08485.31 × 10 6 0.000983
photo bio HJ temp S0.220.3330.266.99
photo bio HJ GI0.3530.01620.0009210.167
photo bio HJ GI S0.3170.1320.30931.0
photo bio HJ GI temp0.1540.01410.0006870.34
photo bio HJ GI temp S0.4850.1010.40940.5
photo bio HJ iso0.1040.1885.06 × 10 5 0.0075
photo bio HJ iso S0.2590.4240.45111.7
photo bio HJ iso temp0.1850.1511.5 × 10 5 0.00363
photo bio HJ iso temp S0.320.4920.37911.8
photo TL0.4780.2327.36 × 10 7 0.000102
photo TL S0.3760.4230.38210.8
photo TL temp0.4120.1972.61 × 10 7 6.01 × 10 5
photo TL temp S0.4530.2880.28310.8
photo TL GI0.0160.3124.34 × 10 6 0.000271
photo TL GI S0.3120.4310.44713.9
photo TL GI temp0.00960.3061.93 × 10 6 0.000332
photo TL GI temp S0.4850.3170.34516.7
photo TL iso0.410.3081.46 × 10 6 0.000264
photo TL iso S0.3840.460.43515.9
photo TL iso temp0.3230.265.48 × 10 7 0.000154
photo TL iso temp S0.4470.3310.36515.5
photo TL HJ0.4780.2326.94 × 10 7 9.94 × 10 5
photo TL HJ S0.360.3940.35511.2
photo TL HJ temp0.4120.1972.51 × 10 7 5.92 × 10 5
photo TL HJ temp S0.4580.2760.27210.9
photo TL HJ GI0.0160.3124.31 × 10 6 0.000269
photo TL HJ GI S0.3120.430.44613.9
photo TL HJ GI temp0.009590.3061.92 × 10 6 0.00033
photo TL HJ GI temp S0.4840.3160.34416.7
photo TL HJ iso0.410.3071.42 × 10 6 0.000258
photo TL HJ iso S0.3820.4520.4316.0
photo TL HJ iso temp0.3230.265.37 × 10 7 0.000151
photo TL HJ iso temp S0.4480.3250.36215.5
photo TL bio0.03540.290.0001490.017
photo TL bio S0.2520.4950.42314.9
photo TL bio temp0.04060.2910.0001160.0252
photo TL bio temp S0.330.3760.41320.9
photo TL bio GI0.3170.03290.001210.0717
photo TL bio GI S0.1040.4130.38517.9
photo TL bio GI temp0.3710.03570.000860.14
photo TL bio GI temp S0.1960.2830.39727.3
photo TL bio iso0.0870.430.0003250.0495
photo TL bio iso S0.2490.50.39320.4
photo TL bio iso temp0.09820.4320.000260.07
photo TL bio iso temp S0.3530.3830.36625.1
photo TL bio HJ0.03470.2890.0001410.0167
photo TL bio HJ S0.2180.4560.44916.0
photo TL bio HJ temp0.03990.290.0001120.025
photo TL bio HJ temp S0.3110.3540.42421.6
photo TL bio HJ GI0.3160.0330.001210.0715
photo TL bio HJ GI S0.1030.4110.38517.9
photo TL bio HJ GI temp0.370.03590.0008590.14
photo TL bio HJ GI temp S0.1960.2820.39827.2
photo TL bio HJ iso0.08540.4290.0003180.0488
photo TL bio HJ iso S0.2430.490.39720.7
photo TL bio HJ iso temp0.09660.430.0002550.0693
photo TL bio HJ iso temp S0.3480.3750.36725.3
yellow0.4860.1628.78 × 10 7 9.64 × 10 5
yellow temp0.4920.1612.31 × 10 7 2.61 × 10 5
yellow GI0.005550.2926.88 × 10 5 0.0127
yellow GI temp0.003290.2722.81 × 10 5 0.00588
yellow iso0.3910.2192.03 × 10 6 0.000336
yellow iso temp0.3870.2185.35 × 10 7 9.14 × 10 5
yellow HJ0.4860.1628.64 × 10 7 9.55 × 10 5
yellow HJ temp0.4920.1612.27 × 10 7 2.58 × 10 5
yellow HJ GI0.005540.2926.83 × 10 5 0.0126
yellow HJ GI temp0.003290.2722.79 × 10 5 0.00584
yellow HJ iso0.3910.2191.99 × 10 6 0.000333
yellow HJ iso temp0.3870.2185.26 × 10 7 9.06 × 10 5
yellow bio0.03510.1021.72 × 10 5 0.00178
yellow bio S0.1230.470.382.88
yellow bio temp0.03240.09125.21 × 10 6 0.000555
yellow bio temp S0.1260.4060.2152.05
yellow bio GI0.05030.01330.0007360.13
yellow bio GI S0.1770.250.3676.9
yellow bio GI temp0.03160.008180.0003150.0632
yellow bio GI temp S0.2990.1770.4747.21
yellow bio iso0.0810.1875.23 × 10 5 0.00821
yellow bio iso S0.1740.3860.4985.37
yellow bio iso temp0.07410.1681.62 × 10 5 0.00262
yellow bio iso temp S0.2060.4240.3244.38
yellow bio HJ0.03510.1021.69 × 10 5 0.00176
yellow bio HJ S0.1120.460.372.91
yellow bio HJ temp0.03240.09125.14 × 10 6 0.00055
yellow bio HJ temp S0.120.40.2092.05
yellow bio HJ GI0.05030.01340.0007320.13
yellow bio HJ GI S0.1770.2510.3686.9
yellow bio HJ GI temp0.03160.008190.0003130.0629
yellow bio HJ GI temp S0.2990.1770.4737.21
yellow bio HJ iso0.08090.1875.13 × 10 5 0.00814
yellow bio HJ iso S0.1690.3920.4935.4
yellow bio HJ iso temp0.0740.1681.59 × 10 5 0.0026
yellow bio HJ iso temp S0.2020.4280.3194.39
yellow TL0.03030.03081.63 × 10 6 0.000763
yellow TL S0.4570.3770.43132.6
yellow TL temp0.02290.02514.55 × 10 7 0.000253
yellow TL temp S0.3280.2810.24226.5
yellow TL GI0.003450.384.84 × 10 5 0.0256
yellow TL GI S0.350.2990.49643.7
yellow TL GI temp0.001960.3722.04 × 10 5 0.0129
yellow TL GI temp S0.4810.3110.32644.9
yellow TL iso0.02490.04123.78 × 10 6 0.00193
yellow TL iso S0.4880.4240.46944.7
yellow TL iso temp0.01860.03351.06 × 10 6 0.000636
yellow TL iso temp S0.3530.3530.29439.2
yellow TL HJ0.03030.03081.58 × 10 6 0.000756
yellow TL HJ S0.4740.3510.40733.8
yellow TL HJ temp0.02290.02514.43 × 10 7 0.000251
yellow TL HJ temp S0.3340.2660.22726.9
yellow TL HJ GI0.003450.384.82 × 10 5 0.0254
yellow TL HJ GI S0.350.2980.49543.6
yellow TL HJ GI temp0.001960.3722.03 × 10 5 0.0128
yellow TL HJ GI temp S0.4810.310.32544.7
yellow TL HJ iso0.02490.04123.69 × 10 6 0.00191
yellow TL HJ iso S0.4930.4160.46345.1
yellow TL HJ iso temp0.01860.03351.04 × 10 6 0.00063
yellow TL HJ iso temp S0.3550.3460.28939.4
yellow TL bio0.3240.3350.01145.54
yellow TL bio S0.3730.4960.32351.3
yellow TL bio temp0.3320.3350.007864.56
yellow TL bio temp S0.3990.4740.42563.1
yellow TL bio GI0.4720.4460.01266.64
yellow TL bio GI S0.1260.2780.32359.4
yellow TL bio GI temp0.4520.4440.00935.88
yellow TL bio GI temp S0.1830.2830.43377.6
yellow TL bio iso0.4350.4460.01588.4
yellow TL bio iso S0.3490.4920.32664.5
yellow TL bio iso temp0.4440.4450.01086.82
yellow TL bio iso temp S0.4040.4990.42377.2
yellow TL bio HJ0.3220.3330.01115.5
yellow TL bio HJ S0.3320.4610.34354.5
yellow TL bio HJ temp0.3310.3330.007694.53
yellow TL bio HJ temp S0.3690.4470.44466.0
yellow TL bio HJ GI0.4720.4460.01266.61
yellow TL bio HJ GI S0.1260.2770.32359.3
yellow TL bio HJ GI temp0.4530.4450.009285.85
yellow TL bio HJ GI temp S0.1830.2810.43477.4
yellow TL bio HJ iso0.4330.4440.01558.36
yellow TL bio HJ iso S0.3390.4990.3365.5
yellow TL bio HJ iso temp0.4420.4430.01066.8
yellow TL bio HJ iso temp S0.3950.4920.42878.4

Share and Cite

MDPI and ACS Style

Sandora, M. Multiverse Predictions for Habitability: Number of Potentially Habitable Planets. Universe 2019, 5, 157. https://doi.org/10.3390/universe5060157

AMA Style

Sandora M. Multiverse Predictions for Habitability: Number of Potentially Habitable Planets. Universe. 2019; 5(6):157. https://doi.org/10.3390/universe5060157

Chicago/Turabian Style

Sandora, McCullen. 2019. "Multiverse Predictions for Habitability: Number of Potentially Habitable Planets" Universe 5, no. 6: 157. https://doi.org/10.3390/universe5060157

APA Style

Sandora, M. (2019). Multiverse Predictions for Habitability: Number of Potentially Habitable Planets. Universe, 5(6), 157. https://doi.org/10.3390/universe5060157

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop