text:
# 1 Introduction

## 1 Introduction

Our universe seems, according to the present-day evidence, to be spatially flat and to possess a nonvanishing cosmological constant . These features, while not yet rock-solid experimentally, are hardly what would have been anticipated by the founding fathers of cosmology. The cosmological constant in particular is, for cosmologists and general relativists, the Great Mistake. And for elementary particle physicists it is the Great Embarrassment. It is fair to say that each school would just as soon see it go away. But in this paper we assume that it will not do so, and that the present evidence will prevail.

The cosmological constant is a peculiar quantity. By definition it has something to do with cosmology. But it also has something to do with the local structure of elementary particle physics, where it represents the stress-energy density $`\mu `$ of the vacuum.

$$_{cc}=d^4x\sqrt{g}\mu ^4=\frac{1}{8\pi G}d^4x\sqrt{g}\mathrm{\Lambda }.$$
(1)
Instead of the parameter $`\mu `$, the cosmologist will use $`\mathrm{\Lambda }`$, as defined above, but expressed in terms of the properties of the spacetime of our distant future, a future dominated by dark energy and exponential expansion:

$$ds^2=d\widehat{t}^2e^{2H_{\mathrm{}}\widehat{t}}(d\widehat{r}^2+\widehat{r}^2d\theta ^2+\widehat{r}^2\mathrm{sin}^2\theta d\varphi ^2).$$
(2)
This matter-free spacetime is equivalent to static DeSitter space, characterized by a horizon radius $`R_{\mathrm{}}`$,

$$ds^2=\left(1\frac{r^2}{R_{\mathrm{}}^2}\right)dt^2\left(1\frac{r^2}{R_{\mathrm{}}^2}\right)^1dr^2r^2(d\theta ^2+\mathrm{sin}^2\theta d\varphi ^2)$$
(3)
via the coordinate transformations

$`r`$ $`=`$ $`\widehat{r}e^{H_{\mathrm{}}\widehat{t}}`$
$`t`$ $`=`$ $`\widehat{t}{\displaystyle \frac{1}{2H_{\mathrm{}}}}\mathrm{}n(1H_{\mathrm{}}^2r^2).`$ (4)
This horizon radius is determined in terms of the asymptotic Hubble constant

$`H_{\mathrm{}}`$ $`=`$ $`\underset{t\mathrm{}}{lim}H(t)`$
$`R_{\mathrm{}}^1`$ $`=`$ $`H_{\mathrm{}}=\sqrt{\mathrm{\Omega }_\mathrm{\Lambda }}H_0`$ (5)
which in turn is determined by gravitational dynamics of the vacuum energy characterized by $`\mu `$.

$$\frac{1}{R_{\mathrm{}}^2}=H_{\mathrm{}}^2=\frac{8\pi G}{3}\mu ^4=\frac{\mathrm{\Lambda }}{3}.$$
(6)
We see from the above equations that, despite the formally infinite extent of the spatially flat Friedmann-Robertson-Walker expansion universe, Eq. (2), the presence of a nonvanishing cosmological constant provides a way of ascribing an intrinsic, observer-independent, size parameter to our universe.Given $`\mathrm{\Omega }=1`$, the only other objective choice of size parameter would seem to be to utilize one of the several landmark times characterizing the history of our universe, e.g. the time of electroweak or strong phase transition, of matter-radiation equality, or of decoupling. Our choice is that of matter-dark energy equality, and appears to us to be the most fundamental.

In this paper we shall be considering an ensemble of universes similar to our own but with different intrinsic sizes $`R_{\mathrm{}}`$.Note that $`R_{\mathrm{}}`$ is not the scale-size $`R(t)`$ characterizing the expansion of our own Robertson-Walker universe. We are not assuming that the fundamental constants are time-dependent. Each universe in the ensemble undergoes its own Big Bang, and is characterized by distinct values of standard-model parameters. The ambivalence between the elementary particle view of the cosmological constant as vacuum energy/pressure and the cosmological view of it as a size parameter for the universe is sharpened by looking at it from this viewpoint. In particular, conventional wisdom would say that all of the basic parameters of the standard model, such as $`\mathrm{\Lambda }_{QCD}`$ or the electroweak vacuum condensate $`v`$, are to the best of our knowledge independent of each other. This means that they are also independent of the vacuum energy/pressure characterized by the scale $`\mu `$, since the cosmological term is just another term in the standard-model Lagrangian density. This in turn implies that $`\mathrm{\Lambda }_{QCD}`$ and $`v`$ are also independent of the size $`R_{\mathrm{}}`$ of the universe. By definition this is not the case for the cosmological term (cf. Fig. 1). For universes smaller than our own, the vacuum energy density grows. And for universes smaller than about 10 km, the vacuum energy density exceeds 1 GeV/fm3, the energy scale of the QCD vacuum.

Does this matter at all? It is at least arguable that it does. To really understand the vacuum state is one of the most important goals of fundamental theory, and it is often presumed that one must go to unified theories, such as string theory, to attain true enlightenment. But this would imply that, for the vacuum state, the cosmological degrees of freedom talk to the elementary particle degrees of freedom such as quark/gluon or Higgs in an essential way. This point of view is reinforced by the fact that the dark-energy term in the action is formally renormalized by the quantum corrections contributed by all the other terms in the action.

If there is an interconnection between dark energy and QCD vacuum fluctuations, we might suspect that interesting things occur when the cosmological vacuum energy scale and the QCD vacuum energy scale become comparable. Perhaps there is a discontinuous change, such as occurs for QED at the electroweak scale. Or perhaps the QCD scale does not wait for such a catastrophe to occur, but changes continuously as the size parameter of the universe changes, in a way which is similar to the way the cosmological constant itself changes. It is this latter option we entertain in this paper, an option we have in fact already suggested . What we assume is, first, that all dimensionful parameters $`X`$ of the standard model may vary with $`R_{\mathrm{}}`$, but that to leading approximation they are straight lines in a log-log plot, i.e. they satisfy a simple renormalization-group equation

$$R_{\mathrm{}}\frac{X}{R_{\mathrm{}}}=\frac{1}{2}\mu \frac{X}{\mu }=p_XX+\mathrm{}$$
(7)
We shall discuss the unspecified corrections and other details a little more in Section 2.

By itself the above assumption includes the conventional-wisdom option, illustrated in Fig. 1, and contains little news. However we in addition assume that all fundamental dimensional parameters $`X`$ flow toward a fixed point occurring for universes of approximately Planck/GUT size, and that this is the only such fixed point. Note that only two parameters are needed to describe the gross dependence of $`X`$ on $`R_{\mathrm{}}`$, and that they are fixed by the value of $`X`$ at the fixed point and the value of $`X`$ observed by us. This is illustrated in Fig. 2.

The first part of the fixed-point assumption should not be any surprise, since we do expect new physics to occur at the Planck/GUT scale. The second part of the assumption is really what provides the motive power for the remainder of this note, and it is not much more at this stage than an application of Occam’s razor. In detail there are probably exceptions to the rule, but perhaps insight can be gained even in the absence of being able to apprehend the exceptions. But leaving aside the possible, even probable complications, what we have at this point is a description of the ensemble of universes we are considering, characterized by the value of $`R_{\mathrm{}}`$, in terms of modified standard-model parameters. And we again emphasize that there are no extra parameters which have been introduced. Therefore we can hope to explore in principle the properties of such universes, using well-defined extrapolations of the laws of chemistry, atomic physics, nuclear physics, etc., with no extra arbitrary assumptions. In particular we can explore the “bandwidth” of features possessed by our own universe. That is, we may try to determine the minimum and maximum sizes $`R_{\mathrm{}}`$ for which nuclear matter exists, or for which hydrogen-burning stars exist, or for which elements as heavy as carbon exist and are produced. These and other examples will be discussed in Section 3, where we shall estimate that if the radius of a universe in our ensemble of universes is within a factor $`\sqrt{2}`$ of our own, the conditions of life as we know them appear to be satisfied.

Up to this point we may regard the ensemble of universes under consideration as an abstract set, a la Gibbs, and the study of their properties as an abstract intellectual exercise, perhaps of value in the long run in understanding either microphysics beyond the standard model, or the macrophysics of our visible universe. However, there is clear motivation to go further, and to presume that such an ensemble actually exists. In particular we may assume that our universe is one member of a multiverse, with the remaining members causally disconnected from us, as discussed extensively by Rees  and others .

If such a multiverse ensemble really exists, then a primary quantity of interest is the number distribution $`n(R_{\mathrm{}})`$ of universes of a given size, defined as

$$R_{\mathrm{}}\frac{dN}{dR_{\mathrm{}}}n(R_{\mathrm{}}).$$
(8)
As mentioned above, we will roughly determine in Section 3 the bandwidth $`\mathrm{\Delta }R_{\mathrm{}}/R_{\mathrm{}}`$ within which the conditions for life as we understand it exist; it is of order 1. Then if

$$\mathrm{\Delta }N=n(R_{\mathrm{}})\frac{\mathrm{\Delta }R_{\mathrm{}}}{R_{\mathrm{}}}1$$
(9)
we may argue that it is not improbable that life should exist in the multiverse. This is just the condition

$$n(R_{\mathrm{}})1$$
(10)
which appears not to be a heavy constraint. The above line of argument, and concomitant set of problems, parallels the lines of argument used to understand our place in our own universe. Why do we live on Earth rather than Mercury or Pluto? The former is too hot, the latter too cold. Is our existence improbable, in the sense that the parameters characterizing Planet Earth are very finely tuned? The simplest answer is that if life as we know it exists elsewhere in the universe, i.e. there is a sufficiently large population of planets to allow the replication of conditions found on Earth, no fine tuning is required. The jury is still out with respect to what that answer is . But we may argue that the question is, at least in principle, a scientific question. And indeed the hypothesis of a multiverse softens the above constraint to only require that the multiverse contain planets with conditions suitable for supporting life as we know it.

There is an even more specific—and speculative—scenario which can be entertained and which is discussed in Section 5. It is a reductionist version of evolutionary cosmology as envisaged by Smolin , utilizing a speculative model of black hole interiors dubbed gravastars by Mazur and Mottola . In this scenario, the interiors of mature black holes are nonsingular and described by the aforementioned static DeSitter metric which characterizes the future of our own universe. This strongly suggests a cosmology of nested black holes. The interiors of the black holes in our universe comprise daughter universes, within which there are granddaughter black-hole universes, etc. Going in the opposite direction, we may surmise that our universe consists of the interior of a black hole existing in a mother universe, which in turn is embedded within a grandmother universe, etc. An important issue in this picture of cosmology is the determination of the various species of daughter universes (supermassive galactic-center black holes, stellar-collapse black holes, $`\mathrm{}`$) and their size distribution relative to the size of the parent. Other important parameters are the fertilities of mothers, i.e. the number of daughter universes created per mother, as a function of $`R_{\mathrm{}}`$ and species. We will try to estimate these parameters from data and astrophysical theory, and then try to estimate, for example, the number of sister universes there are in the multiverse, and thereby to re-examine questions posed above, such as estimating the number of planets, galaxies, and/or universes within the multiverse which might support life as we know it.

By now we have clearly entered a highly speculative level. Indeed we have organized this note such that at the beginning of each new section readers making it to that point can become dismissive and bail out. But in the hope that there is at least one person left reading this paragraph, we continue on.

Why all this speculation? From the point of view of this writer, it is motivated by the gravastar scenario, and the related ideas of emergent gravity and emergent standard model, as advocated by Volovik  and others . The vacuum is visualized as similar to a quantum liquid such as helium at low temperatures. In the “gravastar” scenario, the blackhole universes are droplets of the quantum liquid, with order parameters which depend on the size of the droplet. The cosmological constant is small because in the ground state of the liquid droplet the pressure (which is measured by the cosmological constant!) vanishes, up to surface corrections. In the picture advocated in this note it is not only the cosmological-constant term in the standard-model Lagrangian density which is a size-dependent order parameter, but all the others as well. Indeed for an infinite universe, characterized by an infinite value of the DeSitter horizon radius $`R_{\mathrm{}}`$, the entire standard-model Lagrangian trivializes to a free field theory . All standard-model interactions are therefore viewed as dependent upon the existence of a boundary to our universe. Evidently in the opposite Planck/GUT limit everything becomes strongly coupled.

The crux of this set of ideas lies in the development of a microscopic theory along these lines. And the construction of such a theory may be aided by having a rough picture of the most likely cosmological context for these ideas. It is this which is our primary motivation. But there are a host of obstacles. Some of these are taken up in Section 6, which is devoted to lessons learned and to conclusions, such as they are.

## 2 Standard Model Parameters

The fundamental premise of this paper was already stated above Eq. (7), and we expect this hypothesis to be most accurate for the dimensionful parameters most closely associated with vacuum energy. These are the cosmological-constant scale $`\mu `$, $`\mathrm{\Lambda }_{QCD}`$, the electroweak condensate value $`v`$, and probably a large mass scale associated with neutrino mass, in particular the masses $`M`$ of the heavy gauge-singlet Majorana particles associated with the see-saw mechanism of neutrino mass generation. These masses appear to be in the range $`10^{13}10^{15}`$ GeV, near the GUT scale. For all these quantities, we assume that Eq. (7) holds to good accuracy.

The fact that $`\mathrm{\Lambda }_{QCD}`$ varies with $`R=R_{\mathrm{}}`$ (hereafter we drop the subscript) leads to an important consequence, namely that the strong coupling constant $`\alpha _s(q^2)`$ must also vary with $`R`$. Since

$$\frac{1}{\alpha _s(q^2)}b_s\mathrm{}n\frac{q^2}{\mathrm{\Lambda }_{QCD}^2}b_s=\frac{332n_f}{12\pi }$$
(11)
and

$$\frac{M_{pl}^2}{\mathrm{\Lambda }_{QCD}^2}(M_{pl}^2R^2)^{p_s}p_s\frac{2}{3}$$
(12)
it follows that

$$\frac{1}{\alpha _s(q^2,R^2)}=b_sp_s\mathrm{}nM_{pl}^2R^2b_s\mathrm{}n\frac{M_{pl}^2}{q^2}.$$
(13)
What we have is a new renormalization-group equation for $`\alpha _s`$

$$R^2\frac{}{R^2}\left(\frac{1}{\alpha _s}\right)=b_sp_s+O(\alpha _s)$$
(14)
which we may compare with the usual expression for the running of $`\alpha _s`$ with $`q^2`$,

$$q^2\frac{}{q^2}\left(\frac{1}{\alpha _s}\right)=b_s+O(\alpha _s).$$
(15)
It is important that this behavior holds for $`\alpha _s`$ evaluated at the GUT scale $`M`$. As long as coupling constant unification makes any sense at all, we may infer that the electroweak and electromagnetic couplings must also possess the same behavior. Since

$$\alpha _s^1(M^2,R^2)=b_s\mathrm{}n\frac{M^2}{\mathrm{\Lambda }_{QCD}^2}\alpha _1^1(M^2,R^2)\alpha _2^1(M^2,R^2)$$
(16)
it follows that the electroweak couplings at the weak scale are

$$\alpha _i^1(v^2,R^2)=\alpha _i^1(M^2,R^2)+b_i\mathrm{}n\frac{M^2}{v^2}=b_s\mathrm{}n\frac{M^2}{\mathrm{\Lambda }_{QCD}^2}+b_i\mathrm{}n\frac{M^2}{v^2}.$$
(17)
But because $`M/v`$ and $`M/\mathrm{\Lambda }_{QCD}`$ scale as powers of $`M_{pl}R`$, it follows that

$$\alpha _i^1(v^2,R^2)=(\mathrm{const})\mathrm{}nM_{pl}Ri=1,2$$
(18)
and

$$\frac{\alpha _1(v^2,R^2)}{\alpha _2(v^2,R^2)}\mathrm{tan}^2\theta _W=\mathrm{const}.$$
(19)
Therefore the weak mixing angle $`\theta _W`$ is to good approximation independent of $`R`$, as is the ratio $`\alpha /\alpha _W\alpha (R^2)/\alpha _2(v^2,R^2)`$. The usual diagram of gauge coupling running and unification is shifted in scale as $`R`$ is varied. Because the weak mixing angle does not depend upon $`R`$, the entire figure becomes self-similar (Fig. 3). To leading order, 1/$`\alpha `$ and 1/$`\alpha _{\mathrm{weak}}`$ depend linearly on $`\mathrm{}nR`$ and vanish at the Planck/GUT radius (Fig. 4).

Finally, we shall again invoke Occam and assume that the large Higgs Yukawa couplings $`\lambda `$ and $`h_{\mathrm{top}}^2`$ are no exception to the rule, and that they also obey the same rule, namely that the inverse couplings vary linearly with $`\mathrm{}nR`$ and vanish for $`R`$ at the Planck/GUT scale. The usual renormalization group equations connecting the Higgs self-coupling $`\lambda `$ to the top-quark Higgs coupling $`h_{\mathrm{top}}`$ (with important QCD corrections) remain unchanged. The new equations are again

$$R^2\frac{}{R^2}\left(\frac{1}{h_t^2}\right)=\mathrm{constant}R^2\frac{}{R^2}\left(\frac{1}{\lambda }\right)=\mathrm{constant}$$
(20)
and may be used to determine how $`\lambda `$ and $`h_{\mathrm{top}}`$, evaluated at either the GUT scale or the infrared scale, vary as $`R`$ is varied.

We now return to consideration of other standard-model parameters with dimension of mass, starting with the masses of top quark, Higgs bosons, and electroweak gauge bosons, all of which have a mass formula of the form

$$m=gv$$
(21)
where $`g`$ stands for a generic dimensionless coupling constant. This implies a renormalization-group equation of the form

$$\frac{R}{m}\frac{m}{R}=\frac{R}{v}\frac{v}{R}+\frac{R}{g}\frac{g}{R}=p_v+(\mathrm{constant})g^2$$
(22)
which clearly possesses an order $`g^2`$ “radiative correction” to the leading behavior. What is clearly happening is that the ratio $`m/v`$ is stable, and does not run as a power of $`R`$ but only as a power of $`g`$, i.e of $`\mathrm{}nR`$.

In the case of these particles the corrections are not very important, because their masses are so close to the value of the electroweak vev $`v`$. A more dramatic example is given by the electron mass, which in a sense lies at the opposite extreme. We do not know whether to regard the ratio $`m_e/v`$ as a function of dimensionless coupling constants, i.e. dependent only on $`\mathrm{}nR`$, or as a ratio of fundamental scales, i.e. dependent on a power of $`R`$. In the former case we have

$$\frac{v}{m_e}(\mathrm{}nMR)^n$$
(23)
where $`M`$ is the GUT/Planck mass scale. In the latter case the flow is

$$\frac{v}{m_e}(MR)^{p_e}p_e0.1.$$
(24)
For numerical estimation of the former case we shall choose $`n=4`$ for the electron. Our motivation is simply to assign one power of some $`g^2`$, of typical order of magnitude, i.e. $`(\mathrm{}nMR)^1`$, per mass hierarchy level. Thus for bottom, strange, down, electron, we take $`n=1,2,3,4`$ respectively.

Returning to the question of $`R`$-dependence of electron mass, we can evaluate each case. The results are shown in Fig. 6, and we may regard the shaded region as a one-parameter region of uncertainty. We also plot in Fig. 6 the $`R`$-dependence of the ratio of electron mass to proton mass, which is a crucial parameter for chemistry and condensed matter physics.

The $`R`$-dependence of other small masses present in the standard model should be similarly regarded, especially the up and down quark masses which drive chiral symmetry breaking of the strong interactions and are responsible for the pion mass. For the strong interactions there is to good approximation only the scale $`\mathrm{\Lambda }_{QCD}`$, which by itself determines all masses other than that of the pions. In particular the proton mass is proportional to $`\mathrm{\Lambda }_{QCD}`$, as well as the masses of all mesons and baryons other than the pions (and kaons). But the scale set by the pion mass, whose square varies linearly with the light quark masses and with $`\mathrm{\Lambda }_{QCD}`$, does matter. It is the $`R`$ dependence of the ratio of pion to proton mass which will be the crucial parameter for nuclear physics. Its square is plotted in Fig. 8. In Fig. 8, we also plot $`(m_K/m_\pi )^2`$, which is proportional to the ratio of strange-quark to down-quark mass.

The remaining parameters of the standard model are $`\theta `$, the $`CP`$-violating parameter within QCD, and the CKM mixing parameters, which are closely related to the small quark and lepton masses. In addition there are neutrino masses and mixings. The underlying physics still awaits better understanding, and we have little to add here. These parameters do not appear to be of great importance for what follows in the remainder of this note.

## 3 Properties of Matter in the Ensemble of Universes

One of our main goals is to investigate the properties of universes assumed to be almost the same as ours, i.e. with radii $`R`$ almost the same as ours, and with cosmological initial conditions similar to ours. However we shall also consider more extreme cases, namely universes with radii within 30 orders of magnitude of our own. One reason for considering such large bandwidth is that there will be some properties of these universes which are robust, and do not vary all that much over all those powers of 10. For example the weak and electromagnetic fine-structure constants are in this context robust, having values in this range of $`R`$ which are within a factor two of what we observe. On the other hand it is well known, especially amongst the “anthropic” community, that other properties of our universe are very finely tuned and will only exist over a quite small bandwidth. We shall pay special attention to such “anthropic” constraints, as discussed for example in the book by Barrow and Tipler , and will be interested in the bandwidth in $`R`$ for which they are satisfied.

We shall begin by considering how the properties of elementary particles, of nuclear matter, and of ordinary matter vary with $`R`$. We then investigate how the structure of astrophysical objects of interest, such as planets, stars, etc. vary as $`R`$ is varied.

### 3.1 Elementary Particle Properties

Even for the smallest universe that we shall consider, with radius 100 microns, there is good separation (a factor 10–20) between the electroweak scale and the strong interaction scale. Heavy quarks, electroweak gauge bosons, Higgs particles, etc. are still unstable and will not grossly influence the phenomenology of ordinary matter. A marginal case is that of the strange quark. An estimate of its effect is given by the ratio of kaon to pion mass exhibited in Fig. 8. We see that even for the extreme cases the ratios always stay comfortably above unity, suggesting that we do not err badly in neglecting strange-quark contributions to ordinary matter. The strange hadron masses appear to stay high enough to allow semileptonic weak decays at the very least to proceed.

### 3.2 Nuclear and Atomic Matter

Crucial to the properties of nuclear and atomic matter are the values of the fine-structure constant (here constrained to a reasonable range of values), the ratio of electron to proton mass, the ratio of pion to proton mass, and the neutron-proton mass difference. As long as the electron-proton mass ratio stays small, atomic physics and chemistry will remain recognizable. We see from Fig. 6 that this is in fact the case. Likewise, in Fig. 8 we see that the pion mass stays well below the proton mass over all the range to be considered.

We conclude that over the 60 orders of magnitude we shall consider, it would appear that chemistry and condensed matter physics will be at least existent and reasonably recognizable. As for nuclear matter, it should exist in recognizable form as the pure chiral limit is approached, because the long range force due to pion exchange is not crucial. It might bind the nucleons a little more (or less), but probably not enough to change the phase structure. However as the pion mass increases, there is more potential for trouble. According to Fig. 8, this appears to occur only in the largest or smallest universes that we shall consider. All this will be discussed in more detail in Section 3.4.

On the other hand, as the pion mass decreases, the neutron-proton mass difference varies in a nontrivial way. It is composed of two pieces. The dominant one is due to the mass difference of the up and down (current-) quarks, and the other is electromagnetic . Schematically we may write

$`\mathrm{\Delta }m=(m_nm_p)`$ $`=`$ $`a_1(m_um_d)+b(\alpha m_p)`$
$`=`$ $`a_2\left({\displaystyle \frac{m_um_d}{m_u+m_d}}\right){\displaystyle \frac{m_\pi ^2}{m_p}}+b(\alpha m_p)`$
$`=`$ $`a_3{\displaystyle \frac{m_\pi ^2}{m_p}}+b(\alpha m_p)`$
$``$ $`0.1m_p\left[\left({\displaystyle \frac{m_\pi }{m_p}}\right)^2\alpha \right]`$
where we assume that $`a_i`$, $`b`$, and the ratio of the difference of up and down quark masses to their sum are to good approximation scale-independent. In the last line we have used the accepted values of the two contributions  in approximate form to provide a useful mnemonic. Note that the electromagnetic and quark contributions to $`\mathrm{\Delta }m`$ are of opposite sign. As the chiral limit is approached, the neutron becomes stable and the proton unstable. The latter case is clearly a serious matter for atomic physics and chemistry, which might even cease to exist. However, from Fig. 9 we see that the only cases where this becomes a problem are for universes whose radii are ten to fifteen orders of magnitude larger or smaller than the radius of our universe. The cosmologies for those cases will evidently be nontrivially different from our own, and we will briefly return to this issue later.

### 3.3 Stable Cosmological Objects

Other than black holes, all the large stable cosmological objects exist as a consequence of the Pauli principle. Fermion degeneracy pressure in one form or another provides the repulsion that prevents such objects to gravitationally collapse. This mechanism is so robust that we can expect it to operate over the whole 60 orders of magnitude of radii $`R`$ which we consider. Three obvious classes to consider are planets, white dwarfs, and neutron stars. In these three cases the degeneracy pressure is provided by nonrelativistic electrons, relativistic electrons, and neutrons respectively. We begin by briefly reviewing these cases.

The density of a planet-like object is fixed by the interatomic force, and the spacing of atomic nuclei is of order $`(\alpha m_e)^1`$. This gives for the baryon number of a planet of radius $`r`$ the value

$$BA(\alpha m_er)^3$$
(26)
where $`A`$ is the mean atomic number of the nucleus. (If heavy elements are not produced in the universe of interest, then we take $`A=1`$, and limit our attention to Jupiter-like planets). The chemical binding energy per nucleus is of order the Rydberg, and from this we can determine the total chemical energy and equate it with the gravitational energy in order to determine the characteristic size $`r`$ of the planet:

$$U_{\mathrm{chem}}(\alpha ^2m_e)\frac{B}{A}\frac{B^2}{r}\left(\frac{m_p}{M_{pl}}\right)^2U_{\mathrm{grav}}.$$
(27)
Upon eliminating $`r`$, this leads to

$$B\frac{\alpha ^{3/2}}{A^2}\left(\frac{M_{pl}}{m_p}\right)^3$$
(28)
In a similar way, we may consider white dwarfs, where relativistic electron degeneracy pressure balances the gravitational energy. In that case the baryon number is given by

$$Br^3p_F^3$$
(29)
where $`p_F`$ is the Fermi momentum of the electron plasma. The energy-balance equation is

$$U_{\mathrm{degen}}p_FB\frac{B^2}{r}\left(\frac{m_p}{M_{pl}}\right)^2U_{\mathrm{grav}}$$
(30)
which simplifies to

$$B\left(\frac{M_{pl}}{m_p}\right)^3.$$
(31)
Finally, we may consider the case of the neutron star. It is similar to the white dwarf case. One simply replaces the electron Fermi momentum $`p_F`$ with $`\mathrm{\Lambda }_{QCD}`$ which characterizes the neutron Fermi momentum, and arrives at the same result.

We see that in all three cases the baryon number, hence the mass, of the object scales as the inverse third power of the proton mass, and therefore scales as the appropriate power of $`R`$. The result is shown in Fig. 10. We therefore can anticipate the existence and can understand the properties of these massive objects, throughout the 60 orders of magnitude of $`R`$ we consider. Stars, however, are another matter. The question of whether these large objects ignite and burn, and for how long, depends on details. Before addressing stellar structure we consider some of the finer points having to do with the nuclear force.

### 3.4 The Nuclear Force

The simplest system in nuclear physics is the dinucleon. It is a delicate case, especially in the context of astrophysics, because the nonexistence of bound diprotons and/or dineutrons is needed to keep stars burning , and the existence of a deuterium bound state is an essential ingredient for fusion reactions in stars as well as in big-bang nucleosynthesis. What is needed is the dependence of the binding energies upon $`(m_\pi /m_p)`$. There is no consensus on what the answer is. Depending upon the method, different magnitudes and even signs are obtained . What is important in our application is the value of $`(m_\pi /m_p)`$ for which the deuteron becomes unbound, as well as the value of $`(m_\pi /m_p)`$ for which the diproton and/or dineutron might become bound. We choose here an estimate which lies in the midrange of what is generally considered , and has the sign dictated by naive intuition; as one approaches the chiral limit, the binding energies increase. Our choices are as follows:

$`\mathrm{Deuteron}\mathrm{bound}\mathrm{if}m_\pi /m_p`$ $``$ $`0.16`$
$`\mathrm{Diproton}\mathrm{bound}\mathrm{if}m_\pi /m_p`$ $``$ $`0.08.`$ (32)
We emphasize that these choices are uncertain, but probably by not more than a factor 3. However, it is arguable  that the dinucleon remains unbound even in the chiral limit, in contradiction to the choice made above. But, it will turn out that in what follows we will not consider any region of parameter space where, given the parameters we have chosen, the diproton is bound, so that for us the issue is moot.

Finally, we may consider the mechanism for producing carbon in stars. This depends upon the existence of the anthropically famous triple-$`\alpha `$ reaction

$`{}_{}{}^{4}\mathrm{He}+^4\mathrm{He}`$ $``$ $`{}_{}{}^{8}\mathrm{Be}`$
$`{}_{}{}^{8}\mathrm{Be}+^4\mathrm{He}`$ $``$ $`{}_{}{}^{12}C+2\gamma `$ (33)
with the resonance in $`{}_{}{}^{12}C`$ predicted by Hoyle , together with the absence of a crucial level in $`{}_{}{}^{16}O`$. The parameter sensitivity of this process, which is of order $`\delta E`$ 100 keV in a system with binding energy scales in the 10 MeV range, is discussed by Oberhummer et al. , among others . The result is that an 0.3 percent variation in the overall strength of the nuclear force is enough to strongly modify this delicately balanced mechanism. If such a perturbation were applied to the deuteron, it would change its binding energy by about five percent. We conclude that at most the triple-$`\alpha `$ process represents a parameter sensitivity a factor 20 greater than what one obtains from considering dinucleon binding. However, the actual sensitivity of the triple alpha mechanism may be considerably less, because the dinucleon binding could be more sensitive to the long-range pionic tail of the force than the interactions between compact, closed-shell alpha-particles. It would be helpful to have a good description of the dependence of the intermediate-range, isosinglet, spin-independent attractive force upon pion mass. But at present this seems not to exist.

## 4 Cosmology

In this section we explore how big bang cosmological evolution depends on the ultimate “size” $`R`$ of the particular universe which is created. There are a variety of epochs in the history of a universe which are especially sensitive to parameter variations. Before going into more details we briefly sketch them here to set the stage:

1. We take as initial condition of the universe its state just after inflationary reheating (the assumption that inflation indeed occurs will not be too important), with the initial temperature taken to be of the Planck/GUT scale. The universe is always taken to be spatially flat. The magnitude of the primordial density fluctuations, which eventually account for the observed fluctuation spectrum in the 3 degree microwave background, is in principle a parameter to be specified. In practice we shall choose it to be equal to what it is in our universe, $`\delta _H=\delta \rho /\rho 2\times 10^5`$, independent of $`R`$.

2. The baryon asymmetry of the universe is assumed to be generated in some intrinsic way from unknown, extended-standard-model mechanisms at a very high temperature scale. The details of this mechanism are at present very uncertain. Therefore the $`R`$ dependence of this asymmetry will be treated in a way similar to how the electron mass was treated. We assume that for universes with a size of order the Planck radius, the baryon asymmetry is large, of order unity. We assume that the interpolation from the Planck size to large universes like our own may behave as a power of $`R`$, or as a power of $`\mathrm{}nR`$, each option taken to be an extreme case. The result is shown in Fig. 11. While the uncertainties become large for universes very different in size from our own, at least the dependence upon $`R`$ is monotonic.

3. As the universe cools, the phase transitions at electroweak and QCD scales proceed in a way similar to our universe. (The baryon asymmetry may be modified at the electroweak scale via “sphaleron” effects , and if so it is the modified asymmetry which is shown in Fig. 11.) Differences appear at the epoch of nucleosynthesis, at a temperature of order $`10^3m_p`$. The mechanisms are sensitive to the time at which neutrinos decouple from the plasma and baryonic chemical equilibrium is lost. The abundances of 4He, deuterium, and hydrogen (and even 2He) become sensitive to the parameters and require a detailed discussion.

4. At some very uncertain temperature scale, (cold) dark matter decouples from the plasma and evolves, eventually becoming a major component of the matter density. The physics of this is obscure. We shall assume that the cold dark matter is composed of WIMPs, by which we mean that their interactions with each other and with ordinary matter are characterized by a scale somewhere around the electroweak scale.

5. The epochs of matter-radiation equality and of decoupling of radiation from matter are also parameter sensitive. What happens during these periods provide initial conditions for the subsequent evolution of large-scale structure formation. All of this will require a detailed discussion.

### 4.1 Nucleosynthesis

As the universe cools below the QCD phase transition, quarks and antiquarks bind into mesons and baryons, and the mesons soon disappear. Neutrons and protons are kept in chemical equilibrium by electroweak scattering processes induced by neutrinos. Eventually the neutrinos decouple, the criterion for decoupling being that the expansion rate of the universe exceed the collision rate. The expansion rate for a radiation dominated universe is

$$H^2\frac{T^4}{M_{pl}^2}.$$
(34)
Equating $`H`$ to the collision rate gives

$$H\frac{T^5}{v^4}.$$
(35)
This leads to the criterion

$$\left(\frac{T}{m_p}\right)\left(\frac{v}{m_p}\right)\left(\frac{v}{M_{pl}}\right)^{1/3}.$$
(36)
This result is plotted in Fig. 12.

In our universe, neutrino decoupling occurs at a temperature of about 0.7 MeV. At that temperature the ratio of neutrons to protons has been depleted by about a factor seven, due to the Boltzmann factor containing the neutron-proton mass difference. The remaining neutrons capture into deuterium, which is then converted quickly to 4He by fusion reactions. The net result is a primordial helium abundance of 22 percent or so.

Had neutrino decoupling occurred much earlier, the neutron-proton ratio would have been unity. All the baryons would end up as deuterium, which would then convert via fusion to helium. Conversely, if decoupling were to occur much later, then the neutrons would be removed by the neutrino reactions, and there would be nothing left at low temperatures but hydrogen. We therefore expect the dependence of the abundance of primordial helium to change from very high for $`R`$ smaller than the radius of our universe to very low for $`R`$ greater.

However, for very large or small $`R`$ the situation is more complicated and in fact uncertain. As $`R`$ varies the ratio of pion to proton mass varies, and with it the binding energy of the dibaryons. For large values of $`(m_\pi /m_p)`$, the deuteron probably does not exist, and the fusion reactions are blocked. If the pion mass is very small, then the diproton may be bound, and instead of hydrogen in the final state of proton universes, there would initially be 2He. Further complicating the situation is the $`R`$ dependence of the ratio of neutron-proton mass difference to electron mass, which for small $`R`$ can fall below unity, leading to a stable neutron. In addition, if $`(m_\pi /m_p)`$ is sufficiently small, the neutron-proton mass difference changes sign. The situation is sketched out in Fig. 13, where we identify various regions in the two-dimensional parameter space of $`R`$ and $`(m_\pi /m_p)`$ for which the baryogenesis scenarios qualitatively change. There are seven distinct regions of the parameter space we consider. Regions I, IV, and V are characterized by an unstable deuteron. Regions V, VI, and VII are characterized by a stable neutron. In most of region VII the proton is unstable. In regions I and II neutrinos decouple at such a low temperature that the fraction of neutrons in the mix is less than 5 percent. Consequently those universes evolves into predominantly hydrogen. In region VII the opposite occurs, and the 4He fraction exceeds 90 percent. Only in region III is the situation qualitatively the same as for our universe. If $`R`$ lies between $`10^{24}`$ cm and $`10^{33}`$ cm, this is assured to be the case, although this conclusion rests heavily upon the assumption made in Eq. (32).

We now briefly describe the individual baryogenesis histories for the seven regions we have identified:

Region I: In this region the deuteron is unbound, and the decoupling of neutrinos occurs so late that the $`n/p`$ ratio is less than 5 percent. The result is a nearly pure hydrogen universe. However fusion reactions within stars will not proceed because of the absence of deuterium.

Region II: Again the $`n/p`$ ratio is less than 5 percent, and a predominantly hydrogen universe is formed. But now the deuteron exists, so that in principle stars can burn hydrogen into helium.

Region III: As noted above, this region resembles—and includes—our own universe.

Region IV: In this region the deuteron is unbound. Although neutrinos decouple relatively early, when the $`n/p`$ ratio is not too small, the fusion reactions are blocked. The extra neutrons decay, and we are again left with a hydrogen-dominated universe. But as in Region I, fusion reactions within stars are blocked.

Region V: This region differs from Region IV, because here the neutron is stable. The universe will be mixed hydrogenic and neutron, with fusion reactions again blocked because of the absence of deuterium.

Region VI: In this region, deuterium exists. The $`n/p`$ ratio at neutrino decoupling is not small, so that nucleosynthesis of 4He should proceed. The neutron is stable, but primordial neutrons are presumably found in the helium. Fusion reactions in stars should be able to proceed.

Region VII: In most of this region the proton is unstable and decays to the neutron with positron emission. The $`p/n`$ ratio is large enough (but less than unity!) so that primordial 4He will be produced. Because the deuteron is stable, fusion reactions in stars may proceed. Hydrogen-based chemistry will not exist, although perhaps some deuterium-based chemistry might survive.

It is noteworthy that in all seven regions the electron chemical potential does not vanish. An electron plasma will persist until decoupling occurs at a much lower temperature scale.

### 4.2 Stars

The properties of these regions are perhaps well enough defined that one could go further and map out the subsequent cosmological history in a little more detail. We shall not try to do so here. But before going on to more general cosmological questions, we will consider (within our Region III) the additional oft-cited constraints on the existence of long-lived stars and of the conditions appropriate to the production of carbon and other heavier elements . An immediate reason for doing so is anthropic; we would like to know the bandwidth in $`R`$ within which the changes in standard model parameters are small enough to preserve the conditions in our universe which are conducive to life as we know it.

Quite a long list of “anthropic” constraints exist. Upon examination of the items on that list, it should come as no surprise that the most restrictive by far is the existence of the triple-$`\alpha `$ fusion-reaction chain which allows the production of carbon and thereby the existence of heavier elements. We already mentioned in Section 3 that this constraint could be 20 times more sensitive than the constraints used above regarding the existence of bound deuterium. Examination of Fig. 13 shows that this enhanced sensitivity roughly translates into

$$\left|\mathrm{log}_{10}\frac{R}{R_0}\right|<0.2$$
(37)
or that

$$0.7<\frac{R}{R_0}<1.5$$
(38)
where $`R_0`$ is the radius of our universe. In other words, if the radius of a universe in our ensemble of universes is within roughly a factor $`\sqrt{2}`$ of ours, the standard model parameters are close enough to our own not to upset the conditions necessary for existence of life in that universe.

### 4.3 Large Scale Structure

As the universe continues to cool and expand, the era of matter dominance emerges. In the scenario we consider ($`\mathrm{\Lambda }CDM`$), it is the cold dark matter that is essential in initiating the growth of density fluctuations. The baryons carry less of the energy density, and they stay coupled to the photons for much longer, thereby being unable to fully participate in the growth of inhomogeneities until decoupling is reached.

We review briefly the standard calculations in order to see the parameter dependencies . The abundance of cold dark matter WIMP particles $`X`$ is estimated by equating their rate of production and/or annihilation to the Hubble expansion rate at the time/temperature of WIMP decoupling.

$$n_X\sigma vH\frac{T^2}{M_{pl}}.$$
(39)
Here $`n_X`$ is the number density and $`T`$ the temperature at decoupling. Normalizing the abundance to the abundance of photons, proportional to the cube of the temperature, gives

$$\left(\frac{n_X}{n_\gamma }\right)\frac{1}{M_{pl}T\sigma v}\frac{20}{m_XM_{pl}\sigma v}$$
(40)
where we use the fact that within a factor two the WIMP decoupling temperature is twenty times lower than the rest mass of the $`X`$ particles over a very wide range of parameters .

We may further relate this to the abundance of baryons, by introducing the baryon-to-photon entropy ratio $`\eta `$, already discussed above and depicted in Fig. 11:

$$\frac{\mathrm{\Omega }_X}{\mathrm{\Omega }_B}\left(\frac{m_X}{m_p}\right)\left(\frac{n_X}{n_\gamma }\right)\left(\frac{n_\gamma }{n_p}\right)\frac{20}{m_pM_{pl}\sigma v\eta }.$$
(41)
When this is evaluated for our universe, the cross section estimate is

$$\sigma v\sigma 10^{37}\mathrm{cm}^2.$$
(42)
This cross section is close to the electroweak scale. Therefore we assume, as do many others , that it scales with the inverse square of the electroweak vev $`v`$ and perhaps some power of a coupling constant,

$$\sigma v\frac{\alpha ^n}{v^2}$$
(43)
with $`n`$ taken to be two or three.

We may now look at the $`R`$-dependence of the ratio of dark matter to baryonic matter. It is plotted in Fig. 14. We see that dark matter will dominate over baryonic matter provided the radius of the universe is greater than $`10^4`$ of our own.

We may also determine the temperature $`T_{EQ}`$ when the contributions of matter and radiation are equal, and which signals the onset of matter-dominated expansion of the universe. From Eq. (40) we find

$$\frac{\rho _X}{\rho _\gamma }\frac{m_X}{T}\left(\frac{n_X}{n_\gamma }\right)\frac{20}{TM_{pl}\sigma v}\frac{20v^2}{TM_{pl}\alpha ^n}$$
(44)
from which it follows

$$T_{EQ}20\frac{v^2}{M_{pl}\alpha ^n}.$$
(45)
This is plotted in Fig. 15. Also shown there is the temperature at which radiation decouples from ordinary matter and the universe becomes transparent. The formula which controls this is

$$T_d\frac{\alpha ^2m_e}{40}.$$
(46)
We see that if $`R`$ is less than $`10^5`$ the size of our universe

$$T_d<T_{EQ}$$
(47)
while the opposite is true for larger universes. Therefore for small universes the WIMP degrees of freedom will first feel the Jeans instability, with the growth of fluctuations in baryonic matter occurring later. For the large universes the baryonic fluctuations grow together with the WIMP fluctuations. It is not immediately clear how much of a difference this might make in the creation of large-scale structure.

The evolution of the large scale structure is in general a complex topic. The most straightforward part of the subject consists in the growth of small density perturbations in the linear regime. As mentioned above, we assume that the typical scale of primordial perturbations, present at the earliest epoch we consider when temperatures were at the Planck/GUT scale, are of order

$$\left(\frac{\delta \rho }{\rho }\right)_02\times 10^5$$
(48)
as measured by the cosmic microwave background temperature fluctuations. These density perturbations remain frozen at more or less this value during the radiation-dominated epoch, but grow rapidly once the matter-dominated epoch begins, provided their wavelength is less than the horizon scale at the time $`t_{EQ}`$ of matter-radiation equality. The growth of the amplitude scales with the scale factor of the universe

$$\left(\frac{\delta \rho }{\rho }\right)_R\left(\frac{R}{R_{EQ}}\right)\left(\frac{\delta \rho }{\rho }\right)_0$$
(49)
provided the perturbation is small and one remains in the linear regime. After matter dominates radiation, the Robertson-Walker scale factor of the universe is given by

$$R(t)=R_0\left(\mathrm{sinh}\frac{3}{2}H_{\mathrm{}}t\right)^{2/3}.$$
(50)
When the sinh factor equals unity, one has equal amounts of ordinary (dark plus baryonic) matter and dark energy. We denote this point in time with a subscript $`\mathrm{\Lambda }`$

$$\mathrm{sinh}\frac{3}{2}H_{\mathrm{}}t_\mathrm{\Lambda }=1$$
(51)
and define this as “cosmological freezeout.” For later times, when dark energy is dominant, the growth of fluctuations will cease, and again be frozen in place. Our own universe is in this state of transition, with the present time $`t_0`$ given by

$$\mathrm{tanh}\frac{3}{2}H_{\mathrm{}}t_0=\sqrt{\mathrm{\Omega }_\mathrm{\Lambda }}0.84.$$
(52)
It will in general suffice to equate the present time $`t_0`$ with $`t_\mathrm{\Lambda }`$.

The total amount of growth of initial perturbations, from matter-radiation equality to late times, is therefore simply given (assuming linearity) by the redshift factor between $`t_{EQ}`$ and $`t_\mathrm{\Lambda }`$. Putting in the numbers for our universe, for short-wavelength modes, one finds

$$\left(\frac{\delta \rho }{\rho }\right)_\mathrm{\Lambda }\left(\frac{R_\mathrm{\Lambda }}{R_{EQ}}\right)\left(\frac{\delta \rho }{\rho }\right)_0=(1+z)_{EQ}\left(\frac{\delta \rho }{\rho }\right)_0\left(\frac{T_{EQ}}{T_\mathrm{\Lambda }}\right)\left(\frac{\delta \rho }{\rho }\right)_0(3\times 10^4)\left(\frac{\delta \rho }{\rho }\right)_0.$$
(53)
This marginally contradicts the linearity assumption. Therefore for fluctuations of wavelength large compared to the critical wavelength $`\lambda _{EQ}`$, the total amount of growth will remain in the linear regime. Consequently we expect that the largest structures exhibiting very high density contrast will be limited in size to roughly $`\lambda _{EQ}`$, defined as the wavelength or frequency (in comoving conformally flat coordinates) which is comparable to the horizon scale at $`t_{EQ}`$. Since

$$\lambda _0^t\frac{dt^{}}{R(t^{})}=\frac{3t}{R(t)}t^{1/3}R(t)^{1/2}$$
(54)
it follows that the physical size of this structure at present is

$$r_0=R_0\lambda _{EQ}=\left[\frac{12t_{EQ}}{H_{\mathrm{}}^2}\right]^{1/3}H_{\mathrm{}}^1\left(\frac{8t_{EQ}}{t_\mathrm{\Lambda }}\right)^{1/3}=\frac{8^{1/3}H_{\mathrm{}}^1}{(1+z)_{EQ}^{1/2}}$$
(55)
which scales as the inverse square root of the redshift factor, as shown. For our universe, this implies that the density contrast should be small on scales larger than about 1/400 of the size of the universe, characterized above by $`H_{\mathrm{}}^1`$, and large on scales smaller than that. This is consistent with what is observed.

As we mentioned above, we expect that these gross features of this structure formation are (for the $`\mathrm{\Lambda }CDM`$ scenario) determined by the cold dark matter which is the dominant component of the matter density. On the other hand, the structure on smaller scales may crucially depend upon the baryonic component of the ordinary matter, because, according to our WIMP hypothesis, the dark matter component acts as a collisionless dilute gas, while the baryonic component is more susceptible to nongravitational dissipative mechanisms.

We may now investigate how much things change as the radius $`R`$ of the universe is varied. We have assumed (quite arbitrarily) that the primordial fluctuation scale is $`2\times 10^5`$ independent of $`R`$. As described above, this fluctuation in general grows linearly with $`R`$ from the time $`t_{EQ}`$ of matter-radiation equality until the time $`t_\mathrm{\Lambda }`$ at which dark matter and dark energy (cosmological constant) contribute equally to the Hubble expansion. The former is given in Fig. 15 and Eq. (45), while the latter is given by Eq. (51).

$$T_{EQ}\frac{20v^2}{M_{pl}\alpha ^n}t_\mathrm{\Lambda }H_{\mathrm{}}^1R.$$
(56)
We must convert $`T_{EQ}`$ to $`t_{EQ}`$ using the fundamental relationship between them, valid in the radiation-dominated epoch:

$$t_{EQ}\frac{M_{pl}}{T_{EQ}^2}.$$
(57)
We thereby obtain

$$\left(\frac{\delta \rho }{\rho }\right)_\mathrm{\Lambda }\left(\frac{t_\mathrm{\Lambda }}{t_{EQ}}\right)^{2/3}\left(\frac{\delta \rho }{\rho }\right)_0\left(\frac{RT_{EQ}^2}{M_{pl}}\right)^{2/3}\left(\frac{\delta \rho }{\rho }\right)_0(1+z)_{EQ}\left(\frac{\delta \rho }{\rho _0}\right).$$
(58)
The size $`r_0`$ of the largest structures, relative to the size $`R`$ of the universe, follows from Eq. (55) and has a similar form:

$$r_0\left(\frac{t_{EQ}}{t_\mathrm{\Lambda }}\right)^{1/3}R\left(\frac{RT_{EQ}^2}{M_{pl}}\right)^{1/3}R(1+z)_{EQ}^{1/2}R.$$
(59)
The $`R`$ dependence of the red-shift factor $`(1+z)_{EQ}`$ is shown in Fig. 16. We see that for a rather large bandwidth, of order 10 to 15 powers of ten, the amplification of the primordial perturbations is within an order of magnitude of what is present in our own universe. We also recall that for universes larger than $`10^4`$ of ours, the baryonic fraction $`\mathrm{\Omega }_B/\mathrm{\Omega }_X`$ is small, but not negligibly so. Therefore we may surmise that, whatever the mechanism is that creates black holes in the centers of galaxies, it will probably still be operative in this class of universes as well. However, the basis for this conclusion is very fragile, since it rests upon our assumptions of the properties of dark matter.

## 5 Emergent Cosmology

If the ensemble of universes we have been considering actually exists, then there are anthropic consequences, as mentioned in the introduction. From the behavior found in the previous section, we may conclude that if the number of universes per octave (factor two) in radius $`R`$ is large compared to unity, then it is reasonable that we should be present in the ensemble . There is a caveat; if the mean number of planets per universe (of our size $`R`$) which are appropriate for the support of life as we know it is small compared to unity, then the number of universes per octave needed to make reasonable our existence must be correspondingly increased. The planetary situation is not well understood , so this option is not academic. But either way, it would not seem outrageous that enough universes exist to take care of the problem.

There are probably as many models of multiverses as there are practitioners foolish enough to deal with the idea. In this section we shall play with a specific model, motivated by the idea of emergent field theories, a concept born from analogies with condensed matter physics . The model will also be related to the ideas of evolutionary cosmology developed by Smolin , albeit in a more deterministic framework. As mentioned in the introduction, our reason for indulging in this fantasy is to try to obtain some guidance in the search for a satisfactory microscopic emergent theory.

The basic premise underlying the emergence approach is that the vacuum of particle physics and cosmology is analogous to a quantum liquid in equilibrium at very low temperature. Such a system has essentially zero pressure. But the measure of vacuum pressure is the cosmological constant itself, explaining not only why it should be zero, but why it is not quite zero: a droplet of vacuum of finite size will have pressure due to surface effects. This is just what happens in the DeSitter universe (cf. Eq. (4)).

Chapline et al. , and Mazur and Mottola , have recently carried this notion further, and argue that a black hole is to be considered a droplet of quantum liquid, with a nonsingular interior which is in fact static DeSitter space. The value of the cosmological constant in the interior differs from its value exterior to the horizon and serves as a kind of order parameter. In what we have described, this is generalized to all the standard model parameters, which evidently are also discontinuous across the horizon. This picture is ready-made for the cosmological setting in which we find ourselves: not only does our universe contain a large number of “daughter” black-hole fluid droplets, but our universe itself can be considered the interior of a much bigger droplet, which presumably exists, along with many other “sister” droplets, in a much larger “mother” universe. From this starting point, one easily sees that a genealogy can be defined. The properties of mother and daughter universes will depend upon how different in size they are from our own, and how differently the physics works at those size scales. It is this question that we take up in this section, building upon what was learned in the previous sections., We shall not venture very far beyond one generation in either direction; there will be more than enough uncertainty at this level.

It is easiest and most direct to first consider the daughter universes, because there are some data. There are at least two kinds of daughters—the supermassive, galactic black holes, of horizon size roughly 15 orders of magnitude smaller than the size of our universe, and the stellar-size black holes which are six to eight orders of magnitude smaller still. Despite the greater uncertainty in the underlying astrophysics, we specialize to the former because they are closest to us in size.

Rather than characterizing the black hole size by its horizon radius, it is also useful to give it in terms of the volume of comoving matter needed to form the black hole, which we assume is baryonic in origin. In our universe the mass of such black holes is in the range of $`10^6`$ to $`10^9`$ solar masses , or $`10^{53}`$ to $`10^{56}`$ proton masses, out of a total of about $`10^{79}`$ in the universe. So the fraction by volume of total comoving baryonic matter that ends up in one of these black holes is $`10^{23}`$ to $`10^{26}`$. Taking a cube root gives the fraction in linear scale of roughly $`10^8`$ to $`10^9`$. This should be compared with the fraction in linear scale of about $`10^3`$ for the largest scale structures found in our universe.

Our main purpose in spinning out these numbers is to try to infer the most likely size of our mother universe. Evidently the first rough guess would be 15 powers of ten larger than our own. But by the time one goes out those fifteen orders of magnitude, the cosmology has significantly changed, and it is possible that one must go even further. Let us review what was learned in the previous section for the cosmology of a candidate mother universe, say, of radius $`10^{45}`$ cm.

The early evolution of a mother universe of this size would be similar to our universe. Nucleosynthesis would occur in one of three possible scenarios (cf. Fig. 13), but in all cases hydrogen would predominate in the long run. The three cases are distinguished by the nonexistence of deuterium and/or the relative abundance of primordial helium. As the temperature decreased, decoupling of matter from radiation would occur during the radiation-dominated epoch (cf. Fig. 15). This means that the growth of density contrast would from the outset involve both the baryonic matter and the dark matter. However, the growth factor, which scales with the redshift at the time $`t_{EQ}`$ of matter-radiation equality, is less by about a factor ten than for our universe (cf. Fig. 16). In addition, the ratio of baryonic to dark matter is much less; instead of ten percent, the number is somewhere between 3 percent and 0.01 percent. All of these features will make baryonic structure formation more difficult. What is most important for our consideration here is whether these mother universes can give birth to daughter black holes. Baryonic matter has to aggregate in the potential wells created by the dark matter and undergo gravitational collapse. While it apparently is more difficult for this to happen, it is not clear that the number of black holes that might be created is in fact small compared to unity, when for our universe the corresponding number (for “galactic” black holes) is $`10^{10}`$ or so. We shall make a guess that the fertility curve looks something like what is depicted in Fig. 17; this would allow mother universes to be present, but make it unlikely that grandmothers exist. But we must emphasize the many huge uncertainties involved, not the least of which is the assumption that the primordial fluctuation spectrum does not depend upon $`R`$.

Relative to our universe, the fraction of matter in the mother universe which is baryonic is, as already mentioned, less than for our universe. This might affect not only the frequency of occurrence of black-hole formation, but also the size distribution. We do not try to estimate the effect, mainly out of lack of competence. But it is likely that the ratio of size of mother to daughter indeed grows with overall scale $`R`$, in the way sketched out in Fig. 18. But we emphasize that we are approaching a level of almost complete guesswork.

Despite all these uncertainties, it seems relatively safe to conclude, given our assumptions, that the model of nested black holes for the multiverse allows at most one or two generations of parents, with a number of sister universes small compared to $`10^{10}`$, the number of (galactic) black hole daughters in our universe. It seems very unreasonable to assume a large number of “ancestor” generations, unless the primordial density fluctuations were to increase is magnitude with $`R`$. However, intuitively we would if anything expect the opposite to occur.

Since the size distribution of sister universes span a few factors ten, the fraction with size close enough to our universe to in principle support life as we know it will be a few powers of ten less than the total population of sisters. This means that the total number of universes in such a nested-black-hole multiverse which could support life as we know it is bounded above by ten to a small power. It follows that the overall number of planets in the multiverse that are candidates for habitable environments is not all that different (on a logarithmic scale) from the number in our own universe.

The above inferences are rather strong, and therefore invite an additional critical look: are these conclusions avoidable? In such a soft topic as the contents of this paper, the answer is almost certainly yes. One assumption we have been making, mainly from a desire for simplicity and definiteness, is that the ensemble of universes we consider is parametrized only by the size $`R`$ and nothing else. The other cosmological parameters, such as magnitude of the baryon asymmetry and/or the magnitude of the primordial density fluctuations, may well represent independent initial conditions, unconstrained or at best loosely constrained by the value of the size parameter $`R`$. In such cases the conclusions about abundances of mother and grandmother universes are inoperative. Exploration of such alternatives seems however to be premature, and in any case beyond the scope of this paper.

## 6 Lessons and Challenges

While everything we have discussed is very speculative, it must be admitted that, given the starting hypothesis of size-dependent standard-model parameters, we have been able to look at old questions from a somewhat different perspective. This in itself can be a benefit, inasmuch as a fresh point of view is often a key to making progress. And in fact there are some of the classic Big Questions for which partial answers can be set forward:

1. Why are there such hierarchies in scale amongst standard model parameters?

This question includes the classic “hierarchy problem”, namely the smallness of the electroweak scale $`v`$ relative to the Planck/GUT scale, which stimulates the introduction of weak-scale supersymmetry by so many practitioners. It also includes the question of why the electron mass is so much smaller than the top-quark mass. And it even includes the question of why the cosmological-constant scale is so much smaller than the QCD and electroweak scales, not to mention the Planck/GUT scale.

The answer to this general question, given the multiverse hypothesis, may be that for most universes in the ensemble there is no such huge hierarchy. If the size distribution of universes is maximum for relatively small values of $`R`$, say the GUT scale or smaller, then the typical universe has no large disparity of scales. Only the large, rare, universes like ours enjoy that property as a consequence of the assumed scaling behavior of parameters (which of course must eventually be explained).

2. Why is the fine-structure constant 1/137 so small?

The answer to this famous old question is the same as above: in small universes $`\alpha `$ is not small; only in large ones like ours is it small. There are corollaries which are answered in the same way. The most immediate is the more modern version of the above question: why are the gauge coupling constants at the GUT unification scale so small? And directly related to this question is why the QCD scale $`\mathrm{\Lambda }_{QCD}`$ is so small relative to the GUT scale. All these questions are answered in the same way: because we live in a very large universe.

3. Why is our universe so large?

This is the obvious follow-up question to the previous ones. And the answer to this is weakly anthropic: our universe is large because we inhabit it. The discussions in the previous sections show it could not be otherwise, given the scaling assumptions underlying this note.

But in addition to these questions, there is the most important one, which remains without much of an answer:

4. Why should the assumed “fixed-point” scaling behavior be true?

One response is that it (Fig. 3) looks just as credible as the conventional-wisdom alternative (Fig. 1)—which in itself sheds no light on the above questions. But at best this response is highly subjective and leaves much to be desired. To give a more satisfactory reply would be to relate the scaling behavior to the microscopic theory. This has not been done. But there are some interesting guidelines which the assumed behavior suggests. One concerns the limit of the standard model for infinite $`R`$. In that limit all dimensionless coupling constants vanish, and the standard model becomes trivial . In other words, the presence of nontrivial interactions of the particles with each other depends upon the existence of a nonvanishing cosmological constant. In the emergent, “gravastar” scenario, this states that the standard-model interactions are present only because of the presence of the DeSitter horizon, in the neighborhood of which exists new, beyond-the-standard-model physics. It is as if all the standard-model forces are in some sense Casimir effects. However, the standard kind of Casimir effect, which depends upon size of the system as an inverse power, will not do the trick. There are terms in the standard-model action, such as the Higgs mass term and the cosmological-constant term itself, which do have the typical behavior. But most of them, after appropriate rescaling of fields, depend only logarithmically upon the size parameter $`R`$. To see this, write schematically the standard model Lagrangian as

$$=F^2+\overline{\psi }\overline{)D}\psi +(D\varphi )^2+g\overline{\psi }\psi \varphi +g^2\varphi ^4\mu M\varphi ^2$$
(60)
where the first three terms are gauge, fermion, and Higgs kinetic energy terms, and the last terms are Yukawa, quartic Higgs, and Higgs mass terms respectively.§§§We here conjecture, as in our previous note , that the Higgs mass is the geometric mean of the cosmological and Planck/GUT scales. The covariant derivative is

$$D=gA$$
(61)
and $`g`$ is a generic label for gauge or Higgs coupling; we take $`\lambda g^2`$ because the assumed $`R`$ dependence is then universal. Under the rescalings

$$A=g^1\text{ }A,\psi =g^1\text{ }\psi ,\varphi =g^1\text{ }\varphi ,\text{ }D=\text{ }A$$
(62)
we find

$``$ $`=`$ $`g^2\{\text{ }F^2+\text{ }\overline{\psi }\text{ }D\text{ }\psi +(\text{ }D\text{ }\varphi )^2+\text{ }\overline{\psi }\text{ }\psi \text{ }\varphi +\text{ }\varphi ^4\mu M\text{ }\varphi ^2\}`$ (63)
$`(\mathrm{log}MR)\text{ }.`$
The action is

$$S=\frac{1}{\mathrm{}}d^4x=\frac{1}{\mathrm{}(R)}d^4x\text{ }$$
(64)
with

$$\mathrm{}(R)\frac{1}{(\mathrm{}nMR)}.$$
(65)
The entire Lagrangian density gets multiplied by a factor $`\mathrm{}nM_{pl}R`$, as if the Planck constant itself is scale-dependent, vanishing in the limit of infinite $`R`$.

We have not considered in this note such a possibility, and have in fact essentially set the Planck constant, the speed of light, and the Planck mass to unity, not allowing them to vary with $`R`$. As long as the universes in the multiverse are causally disconnected from each other, this can be defended as no more than a convention in the choice of units . However, if there is a connection between the universes, such as in the nested-black-hole scenario, then it is no longer obvious that this is a safe assumption. Relaxation of such an assumption might in fact lead to additional insight. However, exploration of this possibility is beyond the scope of this work.

In addition to the questions above, for which our approach might provide some insight, there are others for which our present lack of understanding is highlighted, and which need better answers in order to sharpen the consequences of the scaling assumptions which we have made. These include:

5. What is the mechanism by which the electron and the light quarks get their mass?

This is often viewed as a minor detail in the grand scheme of standard model problems. But in the context of this note, the lack of understanding of the origin of light quark and lepton masses is translated into relatively great uncertainty in the understanding of the relationship of our universe to other universes of different size.

6. What is the nature of the dark matter?

7. What is the origin of the baryon asymmetry, and what determines its magnitude of $`\eta =3\times 10^{10}`$?

8. Why is the value of the primordial density fluctuations $`(\delta \rho /\rho )_0`$ equal to $`2\times 10^5`$?

These three cosmological questions are hardly novel ; they are evidently crucial to better understanding the properties of the ensemble of universes we consider. Perhaps the only novelty is that we omit (here) the question of the “small” cosmological constant, usually added to the above list.

It is also worth noting that, while the flatness problem (why $`\mathrm{\Omega }=1`$) and the scale-invariant spectrum of primordial fluctuations represent something of a triumph for the idea of inflation, there remains no good answer to the eighth question: the magnitude of the primordial fluctuation spectrum is simply fit to the data, and not understood at all from more fundamental considerations. And, as discussed at the end of the previous section, it is possible that these parameters should be considered as independent characterizations of members of the ensemble of universes, i.e. as initial conditions not strongly dependent upon the size parameter $`R`$.

Finally, there are the lessons, if any, which are learned from this exercise that may be applied to the hypothesis of emergence. The idea of emergence provides some motivation for the scaling behavior assumed from the beginning of this note. But it has many daunting problems associated with it:

1. Why are violations of Lorentz covariance so small?

Condensed-matter analogs of emergence suggest in general that symmetries such as Lorentz covariance are just low energy approximations. At high enough energies deviations are to be expected. But experiment severely limits such deviations. For example, noncovariant corrections to charge-renormalization, an ultraviolet-sensitive quantity, are limited  to less than one part in $`10^{31}`$. This comprises a staggeringly restrictive constraint. It would seem essential that there be a very small parameter which characterizes the violations. And the scaling behavior of parameters studied here suggests that a necessary (but far from sufficient) condition for the Lorentz-violating terms in the Lagrangian is that they scale as inverse powers of the radius $`R`$ of the universe. If this is so, then very small universes exhibit very little symmetry, while the very large ones like our own exhibit Lorentz symmetry, etc., with very small corrections.

2. What is the structure of event horizons?

In the nested-black-hole, or “gravastar” scenario, there is “new physics” at horizons. This is endemic in the condensed-matter analogues . And in our picture, standard-model parameters (including the cosmological constant) are discontinuous across horizons, indicating that at the surface of discontinuity conventional-physics descriptions of what is going on are incomplete. There also appear to be violations of the weak energy conditions of classical general relativity . One manifestation of this appears to be that there are large classes of null geodesics (in particular those which have nontrivial transverse motion) which are “bound” to the horizon. There is a nontrivial problem here of providing a consistent description.

In addition, if our universe is to be regarded as the DeSitter interior of a gravastar, then there must be in our universe preferred comoving observers, presumably not ourselves, with respect to which there is the “physical” horizon associated with our black-hole interior. It then becomes an interesting question as to where we should regard ourselves relative to these central observers: how far away are they, and in what direction? Might there be observational issues associated with this preferred center of our universe? While these are quite interesting questions, they also lie beyond the scope of this note.

3. How are gravastars formed?

If the gravastar picture is in fact viable, then there must be a time evolution of the “new physics” which is associated with the horizon. But for large black holes, it is hard to find an intrinsic, local parameter associated with the horizon, because classically it can be regarded as an artifact associated with a choice of coordinate systems. In the emergence scenario, general covariance is only a low energy approximation. This implies that the description of gravastar formation will require a “best” choice of coordinates. What should be chosen?

This is only one of the difficult issues involving gravastar formation. Another involves rotation: no “eternal rotating gravastar” generalization of the nonrotating case has been found.A default option is to invoke an “eye of the hurricane” model. Choose the Kerr metric for the exterior, and static DeSitter space for the interior. Then build an appropriate interpolating boundary layer with an exotic spacetime, which probably contains vorticity. In nested-black hole cosmologies, it is necessary that all the black-hole universes (including ours) are characterized by a value of spin as well as mass. This is not only a complication, but also an opportunity for linking standard model discrete symmetry violations, in particular $`CP`$, to the existence of a spin axis for the new physics at the DeSitter horizon—physics which presumably controls the nontrivial interaction features of the standard model.

On the more positive side, some insight on the history of gravastar formation might be gleaned by comparing the formation of a daughter black hole with the formation of our own parent universe. The characteristic time for the formation of the daughter can be easily taken to be at the very least many millions of years, a timescale much larger than the size of the gravastar. If we assume the same for our parent universe, it follows that the formation time for our universe should be considered to be much larger than the size parameter $`R`$—in other words orders of magnitude larger than $`10^{10}`$ years . This might imply that the formation time of the “new physics” on the DeSitter horizon likewise is long compared to $`10^{10}`$ years. Some kind of cosmological “bounce” scenario  might have the best chance of providing a concrete implementation of this inference.

4. What is the microscopic physics underlying the emergence scenario?

This question remains unanswered. Necessary conditions are that the gauge bosons of the standard model, as well as the graviton, should be considered collective modes of the presumed “quantum liquid” vacuum. They quite likely should be all considered Goldstone modes  associated with various kinds of spontaneous symmetry breakdown. The pattern of internal symmetries, especially in the fermion representations, must be an essential clue.

Finding the answer to this last question may well afford the best chance of turning the very speculative material in this note into something considerably more concrete.

## Acknowledgments

It is a pleasure to thank my colleagues at Stanford for many helpful discussions, especially R. Adler, P. Chen, and M. Weinstein. I also thank S. Beane and M. Savage for valuable discussions on the chiral limit of nuclear physics.