Multiverse Predictions for Habitability: Stellar and Atmospheric Habitability

Stellar activity and planetary atmospheric properties have the potential to strongly influence habitability. To date, neither have been adequately studied in the multiverse context, so there has been no assessment of how these effects impact the probabilities of observing our fundamental constants. Here, we consider the effects of solar wind, mass loss, and extreme ultra-violet (XUV) flux on planetary atmospheres, how these effects scale with fundamental constants, and how this affects the likelihood of our observations. We determine the minimum atmospheric mass that can withstand erosion, maintain liquid surface water, and buffer diurnal temperature changes. We consider two plausible sources of Earth's atmosphere, as well as the notion that only initially slowly rotating stars are habitable, and find that all are equally compatible with the multiverse. We consider whether planetary magnetic fields are necessary for habitability, and find five boundaries in parameter space where magnetic fields are precluded. We find that if an Earth-like carbon-to-oxygen ratio is required for life, atmospheric effects do not have much of an impact on multiverse calculations. If significantly different carbon-to-oxygen ratios are compatible with life, magnetic fields must not be essential for life, and planet atmosphere must not scale with stellar nitrogen abundance, or else the multiverse would be ruled out to a high degree of confidence.


Introduction
The multiverse hypothesis, which posits that other universes with different laws of physics exist, is an intriguing idea in theoretical cosmology that has so far proven challenging to test [1]. This paper is part of a broader series aiming to rectify this, by generating a plethora of predictions within the multiverse framework regarding the nature of habitability [2][3][4][5][6][7]. The core of this process is the requirement that the multiverse can only be a consistent theory of cosmology if it predicts that our presence in this particular universe is not too improbable; one way of falsifying the multiverse is to find that it predicts that the vast majority of complex (multicellular) life exists in universes with features different from our own. Our contribution to this procedure lies in the recognition that the distribution of complex life, and so observers, throughout the multiverse, depends heavily on the assumptions we make about the nature of habitability. Thus, certain habitability conditions, that are otherwise quite widely discussed, are incompatible with the multiverse. If we ultimately find that the requirements for complex life are incompatible with the multiverse, we will be able to falsify the theory, to a calculable level of statistical significance. Conversely, if we ultimately determine that all currently unknown habitability conditions turn out to be in line with multiverse expectations, we will accrue a long list of supporting evidence for the theory.
It remains to check the compatibility of each habitability condition with the multiverse framework by systematically incorporating them into our calculation of the distribution of observers throughout the multiverse, and the subsequent calculation of the probability of our observations. To this end, we have organized this endeavor into several papers on the topic, each dealing with a loosely overarching theme. The current paper explores several aspects relating to properties of planetary atmospheres, and stellar activity. The two are tightly related, and considered by many to be essential for the maintenance of planetary habitability.
The compatibility of a habitability condition H with the multiverse is determined by the probability of observing our values of the fundamental constants. This is communicated through the Bayes factor, which is defined relative to the baseline case where atmospheric effects are not important H 0 by BpHq " BpHq{BpH 0 q, where BpHq " Ppα|Hq Ppβ|Hq Ppγ|Hq Ppδ u |Hq Ppδ d |Hq (1) and Ppx|Hq " minpPpx ă x obs |Hq, Ppx ą x obs |Hqq, for the fine structure constant α, the electron to proton mass m e {m p " β, the proton to Planck mass m p {M pl " γ, the up quark to proton mass m u {m p " δ u , and the down quark to proton mass m d {m p " δ d . The probability of observing particular values of the constants is defined through the probability density function ppx|Hq 9 p prior pxq Hpxq, as described in more detail in [2]. For the baseline habitability condition H 0 , we take the most successful account of our observations we have considered, which is that complex life requires light from a relatively narrow spectral band for photosynthesis, and that the habitability of a planet is directly proportional to the amount of entropy it receives from incident starlight [2,4].
In Section 2, we discuss generalities of stellar properties, and how these vary with physical constants, deriving expressions that will be crucial for the rest of the paper. In Section 3, we discuss atmospheric loss processes, focusing in particular on extreme ultra-violet (XUV)driven energy limited escape. Determining the importance of this process as constants vary necessitates determination of a great many factors, including stellar spin-down history, initial atmospheric mass, and mass required for surface water retention, which we detail within. Section 4 is dedicated to stellar wind stripping present on planets without an intrinsic magnetic field, and the conditions for planetary magnetic fields to arise.
We find that the significance of atmospheric properties depends on which additional habitability assumptions are made. If we take that an Earth-like carbon-to-oxygen ratio is required for life, as is commonly assumed, then the atmospheric conditions we consider do not strongly affect the probabilities we compute, and so they are neither favored nor disfavored by the multiverse. However, if we adopt the stance that life does not depend on the carbon-tooxygen ratio, several atmospheric conditions do strongly affect the multiverse probabilities. Both the idea that atmospheric mass scales linearly with stellar nitrogen abundance and the idea that planetary magnetic fields are required for habitability cause the probability of our observations to significantly drop, and so both these conditions are incompatible with the multiverse hypothesis. The strategy to test the multiverse is then to check whether this prediction is correct; if life indeed does not depend on planetary carbon-to-oxygen ratio, but either of these other two conditions is found true, the multiverse will be ruled out to high significance.

How Do Stellar Properties Change in Other Universes?
Changes in stellar properties were among the first aspects to be investigated within a multiverse framework. Refs. [8][9][10] worked out how the properties of stars such as mass, lifetime, and luminosity change when constants vary. Ref. [11] discuss photosynthetic potential of starlight. Much work has been done on how different nuclear stability thresholds affect stellar fusion: Refs. [12][13][14] investigated the effects of diproton stability. Refs. [15][16][17][18] discuss the effects of alpha burning. Refs. [19][20][21] investigate deuteron stability. Ref. [22] investigated the consequences of tritium stability. Refs. [23,24] discuss non-nuclear energy production pathways. Ref. [25] discuss the sizes of white dwarfs and neutron stars. Refs. [2,26] discuss the entropy production as a key to habitability.
However, all of these previous studies have so far neglected some of the finer-grained stellar properties, which may nevertheless be just as important for determining the habitability of a planetary system. Among these are properties of stellar coronae, magnetic fields, Sunspot fraction, stellar wind, rotation, and X-ray luminosity. In part, this neglect may be due to prudence on the previous authors' parts, as many of these aspects remain imperfectly understood theoretically, making extrapolation of their behaviors to different universes fraught with potentially misplaced certainty. However, much progress has been made in the understanding of many of these aspects in recent years, and we take advantage of these recent advances to establish a first attempt at determining how these properties may differ in other universes.

Stellar Properties
Expressions for stellar mass, radius, temperature, luminosity, and lifetime in terms of fundamental constants are all already well known (see, e.g., [27]), so we merely reproduce them here: The symbol λ " M ‹ {M Ch is a dimensionless parameterization of stellar mass in terms of the Chandrasekhar mass M Ch " 122.4M 3 pl {m 2 P " 1.4M @ . In these and all following expressions, the functional dependence on constants is derived using physical arguments, and the coefficients are set to accurately reproduce the correct values for our Sun, for the observed values of the physical constants.
In addition, we will need the following expressions for the mass, density, orbital location, total incident power, and day length of an Earth-like planet, which is defined as both temperate (can maintain liquid surface water) and terrestrial (can retain heavy but not light atmospheric gases): Though there will be a certain tolerable range for each of these parameters, we specify to the Earth's values for our calculations. Additionally, note that the temperate requirement has dictated that the incident stellar power is evaluated at a temp ("AU for our values), making this quantity independent of stellar mass.

Speed of Stellar Wind
The escape velocity of a star is v esc " For the Sun, this is 618 km/s. The speed of solar wind is around 400-1000 km/s, roughly the same order of magnitude. This results from the fact that the escaping wind is nonthermal, as particles that have enough energy to make it off the Sun usually have a surplus of the same order. This is larger than the thermal sound speed, which invariably depends on height. For the photosphere, The sound speed of the corona is higher, as discussed below.

Scale Height
The scale height of a star is given by a competition between thermal and gravitational processes as This is 100-1000 km for the Sun, and sets the scale for many processes, including the granule size and typical magnetic flux tube length.

Stellar Magnetic Field
The magnetic field at the stellar surface is created by a highly complex and incompletely understood dynamo mechanism [28,29]. However, the details of the precise mechanism are unimportant for determining the overall field strength, which is set by equipartition of energy as [30] For the photosphere pressure, we use P photosphere " T 4 ‹ , as appropriate for a n " 3 polytrope, which describes stellar structure well [31]. The numerical value matches the observational quantity B surf " 2G. This yields an estimate for the total field strength at the surface, which consists of both open field lines that contribute to the star's long range magnetic field, as well as highly complex field configurations that do not. The long range field is related to the total strength by where f open is the fraction of field lines which are open. It is this factor that introduces rotational dependence to the stellar magnetic field.

Fraction of Open Field lines
The fraction of stellar magnetic field lines which are "open" (i.e., extend to infinity, rather than form a closed loop) depends both on stellar rotation and temperature. This was postulated to depend on Rossby number in [32] as f open " 0.55 expp´2.03 Roq (8) where Rossby number is the ratio of rotation period to convective turnover time, Ro " P rot {τ conv . For the convective turnover time, we use the expression from [33]: where we have neglected terms that cause shutoff for large temperatures. The turnover temperature is set by molecular absorption processes, p . This is normalized to yield a Rossby number of 1. 96 and an open field line fraction of 0.01 for the Sun. The coefficient p q, and is normalized to be 246.4 days for our Sun. Expressing this in terms of fundamental parameters depends on the distribution of stellar rotation periods, which is discussed below.

Corona
The corona is the hotter, much less dense outer layer of a star. Its properties are continuous with the star's extended stellar wind region of influence, and is the source region of most of the variable activity leading to space weather.

Density of Corona
In the formalism of [34], the density of the corona (at the transition region) is determined by the equilibration of heating and cooling processes. The heating rate is given by Q heat " ρ corona v 3 {λ c , where λ c is the granular scale, roughly set by the scale height H " c 2 s {g. The cooling rate for bremsstrahlung is Q cool " n e n p σ T v , where σ T " 8π{3α 2 {m 2 e is the Thomson cross section and " αm e is the typical energy exchange [10]. These are equal when This is equal to 10´1 6 g/cm 3 for the Sun.

Temperature of Corona
The corona is about two orders of magnitude hotter than the photosphere, which has proven puzzling to explain theoretically for many years. Consequently, various competing theories have been developed to explain the anomalously high temperature [35]. Perhaps the most popular account is that of Alfvén wave heating, which posits that energy is transferred to the corona from the stellar interior by turbulent plasma oscillations. In the following, we only consider this theory, which gives the heat flux as [36]: Here δv 2 " T{m p and v A " B{ ? ρ corona . This determines temperature through the diffusion equation S "´κ th ∇T " κ th T{H ‹ . From [28], the thermal conductivity of a stellar plasma is where log Λ "5-20 is the Coulomb logarithm, which has mild parameter dependence, but can be ignored. This can be solved for T to yield

Mass Loss Rate
According to [31], many analytic mass loss formulas have no strong theoretical justification. Whatever the underlying mechanism for solar wind, it is constrained by the continuity equation to obey This is normalized to yield 2ˆ10´1 4 M @ /yr for the Sun. With this, we may ponder whether in some universes the stellar wind is strong enough to deplete stellar material before the available nuclear energy is exhausted; in such universes, type II supernovae would not occur, with stars instead ending their lives having blown off material to the point where fusion ceases. We find t ‹ 9 M{M ‹ " .16α 5 β 1{4 {pλ 3{2 γq " 3ˆ10 6 , so that if α were a factor of 20 lower, this would indeed be the case. However, this may not preclude the distribution of heavy elements into the interstellar medium, if enough reach the wind-launch site. More work is needed to determine whether this mechanism can be at play, whether heavy elements collect in the stellar core, or whether the strong wind effectively extinguishes the star before any heavy elements are created. In any case, including this boundary in parameter space does not appreciably affect the probabilities we compute.

Alfvén Radius
The Alfvén radius is the point at which an appreciable azimuthal velocity component develops. This is set by Throughout we take the Sun's Alfvén radius to be R A " 24R @ , though it can vary by a factor of 2 throughout the solar cycle [37]. By the continuity equation, the quantity ρ v r 91{r 2 . The radial dependence of v r can be found using Parker's model of solar wind, which gives If we define the sonic radius R s " GM ‹ {p2c 2 s q, then for r " R s , this gives v r Ñ 2c s logpr{R s q 1{2 [38], though to first approximation the logarithmic dependence can be neglected.
If B is primarily dipolar, Bprq " B ‹ pR ‹ {rq 3 , and we find For more generic magnetic field profiles B " r´q, the fourth root is replaced by 1{p2q´2q.

X ray Luminosity
A star's X-ray luminosity, which is an important driver of planetary atmospheric loss, is greatly enhanced with respect to the thermal contribution by dynamo processes. As such, X-ray luminosity is found to correlate well with both magnetic activity and rotation speed for slowly rotating stars [39]. For stars with rotation periods less than a few days, however, the X-ray luminosity is found to saturate to about 10´3 of the bolometric luminosity. The origin of this is not well understood, but could be due either to the saturation of surface magnetic flux, or internal dynamo [40], representing a qualitatively different regime of energy transport. These two regimes can be encapsulated with the following expression which reproduces the linear rotation-activity relation between X-ray luminosity and magnetic flux found in [41]. Here, we have defined a convective speed in terms of the convective turnover time in Equation (9) as v conv " R ‹ {τ conv .

Rotation
Since stellar activity depends on rotation rate, and stellar rotation decreases over time, the majority of a planet's atmospheric loss may occur during the initial phase of stellar evolution. Here, we derive expressions for initial stellar rotation as well as spindown rate.

Initial Stellar Rotation
Stars are observed to have a spread of rotation periods within the span of several days, with periods that increase as they age [42]. At formation time, one may expect that stars inherit their rotation from their collapsed dust cloud, but an order of magnitude estimate reveals that the angular momentum of the dust cloud vastly exceeds stellar angular momentum [43]. Indeed, a star possessing that much angular momentum would exceed the critical breakup velocity, and would quickly jettison its material. Instead, the star radiates angular momentum through its surrounding disk until it drops below the breakup speed, and can coalesce [42]. This process results in initial stellar rotation frequencies being close to their breakup velocity, as observed in [44]:

Stellar Spindown Time
Stars lose angular momentum throughout their evolution via stellar wind. While a star's angular momentum is given by J " MR 2 ‹ Ω, to estimate angular momentum loss we must keep in mind that the stellar wind travels radially outward until the Alfvén radius, and so angular momentum loss is given by 9 J " 9 MR 2 A Ω [45]. This increased lever arm greatly enhances spindown, and also introduces extra rotation dependence, as the Alfvén radius depends on spin. A linear dependence R A 9Ω leads to a cubic evolution equation for Ω, as first discussed in [46]. Additionally, a qualitative shift in spindown behavior empirically occurs when the rotation frequency exceeds a critical value, akin to the convective turnover time given in Equation (18). This leads to the following equation governing the evolution of rotation [42]: This also sets the spindown time as For stars rotating more rapidly than the convective turnover time, spindown is set by the star's convective churn, rather than rotation. Below this, the evolution 9 Ω " Ω 3 leads to the well established Skumanich law, P rot " ? t [47]. For fast rotators, the decay is instead exponential. These are all the properties of stars we will need to model our habitability effects in the sections below.

How Do Atmospheric Properties Differ in Other Universes?
We now turn our attention to planetary atmospheres, and whether their character is substantially different in other universes. In particular, we ask what physics determines that the Earth's atmospheric mass is six orders of magnitude less than the planet's mass, how this compares to the minimum needed for several habitability considerations, and whether the expected atmospheric mass is lower than these thresholds for different values of the fundamental constants.
At first glance, it may seem strange to attempt to explain atmospheric mass fraction in terms of fundamental constants. After all, the solar system alone exhibits an enormous diversity of atmospheric mass fractions amongst its planets, from almost zero around small inner rocky bodies to nearly unity for the gas giants. Indeed, atmospheric mass seems to depend on a great number of variables: planetary mass, interior and surface chemistry, orbit, evolution, flux, and the presence or absence of life [48]. Even Venus, though remarkably similar to Earth in orbit and mass, has an atmosphere 90 times Earth's. However, closer inspection reveals hidden regularity; Venus's atmospheric nitrogen content is only 3-4 times that of Earth's, placing it at the same order of magnitude [49]. Its carbon dioxide content, which comprises the bulk of the atmosphere, is the same order of magnitude as that found dissolved in Earth's oceans and compressed into sedimentary rock [50]. Even Venus's initial water content is estimated to have been similar to Earth's [50] (though recent work indicates that even if its initial water content were similar, it may not have ever been able to condense from a steam atmosphere to form an ocean [51]). Evidently, this diversity stems from the different phases each species can undergo, rather than the primordial abundance of each element, giving hope that the overall mass fraction may be understood by processes operating in the early solar system, as well as galactic element abundances. Furthermore, if this is the case, we have hope of extrapolating these values to other universes.
In the following, we focus on nitrogen, as the only gas which is noncondensible under temperate conditions, and present in appreciable quantities. Its presence is essential for stability of liquid water on Earth's surface [52]. It was estimated in [53] that the nitrogen contained in the Earth's mantle is between 3-10 times that of Earth's atmospheric nitrogen, a ratio that is certainly affected by the presence of other species, but is likely to hold as a rough order of magnitude estimate under a range of conditions [49]. Earth's atmospheric nitrogen has remained constant to within a factor of two over the past 3 Gyr, as evidenced by analyzing raindrop imprint size [54] and the isotopic composition of quartz [55].
In the following, we consider two explanations for the magnitude of Earth's nitrogen abundance, corresponding to two different plausible sources: late accretion by chondrites, and initially, as dissolved material in Earth's original building blocks. Each of these hypotheses has different implications for the amounts of nitrogen on planets elsewhere in our universe, as well as throughout the multiverse. Additionally, it is an open question how planetary nitrogen abundance scales with initial stellar system nitrogen abundance, which has important implications for the multiverse, as we happen to be very close to a boundary beyond which nitrogen abundance is reduced by a factor of 270. At the two extremes, the dependence may be linear, if the nitrogen content of solar system bodies was not close to their carrying capacity, or independent, if the bodies were saturated. The dependence probably lies somewhere between these two extremes, but we report how adopting each assumption alters our multiverse probabilities, which serves to bracket the upper and lower limits for our calculations.
We then consider three atmospheric mass thresholds that are plausibly related to habitability. The first is the amount of atmosphere that can be stripped away by stellar flux. The second is related to the pressure necessary to maintain liquid surface water. Third, the mass needed to buffer diurnal temperature changes. Finally, we consider the possibility that only initially slowly rotating stars in our universe are capable of retaining their atmospheres, and assess the compatibility of this hypothesis with the multiverse.

Possible Sources of Atmosphere
The fact that Earth possesses an atmosphere containing volatile constituents is somewhat of a mystery, given that the conditions during Earth's formation were much hotter than their condensation temperatures. Naively, this would result in inner planets that are almost completely comprised of refractory elements, which is manifestly not the case. In the following, we consider two leading theories for the origin of Earth's nitrogen atmosphere: delivery during late accretion from outer system bodies, and as a result of initial accretion from nitrogen dissolved in Earth's original building blocks.

Initial Atmosphere Delivered during Accretion
The classic account for Earth's volatile budget is from planetesimals initially situated outside the solar system's ice line, where temperatures were below the condensation point of volatile species. It has been estimated that up to 7.5 atmospheric masses could have been delivered by carbonaceous chondrites after the main phase of planet formation was completed [56]. This account has the simplicity of explaining the origin of Earth's atmosphere and ocean by a single common source. Additionally, it can readily explain the hierarchy of why Earth's ocean is "100 times more massive than the atmosphere, as [57] demonstrated that the H 2 O/N 2 impact degassing ratio is "100 for a range of different chondrites. Finally, we would like to stress that in this scenario, final atmospheric mass will be highly stochastic, as the material delivered through late accretion is dominated by few large bodies [58]. Thus, while we compute the expected value, it should be kept in mind that this scenario yields a distribution of atmospheric mass ratios.
In [7], we derive the planetary ocean mass fraction delivered via planetesimal accretion during planet formation in terms of the amount of material delivered during late accretion. In this scenario, the atmospheric volatiles are delivered in the same manner. Therefore, we may posit the atmospheric mass fraction to simply be For details on how this expression was obtained, we refer the reader to [7].

Atmospheric Mass as a Result of Accretion by N-Rich Bodies
Here, we follow [59] by considering that Earth's nitrogen was delivered during accretion in the form of dissolved N inside rock and metal. We may then derive the total amount of resulting nitrogen as a function of body mass, with the presumption that only nitrogen in the interior of these planetesimals will be incorporated into the planet's final budget.
The initial nitrogen fraction of a planetesimal is f N " m N {m pp , where m pp is the mass of the planetesimal. The resultant nitrogen budget is obtained through the magma ocean and core asf For the fraction of nitrogen dissolved in the magma ocean, we use [60]: where p 1 and p 2 are coefficients, taken here to scale as p i 9Ry 4 , with Ry the Rydberg constant that dictates the electronic energy scale. The quantity f O 2 is oxygen fugacity, and will depend of the primordial abundances of the two elements. The partial pressure can be rewritten in terms of atmospheric nitrogen mass as This can then be used to find an equation determiningf N : where for cleanliness we have defined k 1 " ζτ{p 1 , k 2 " ζ 2 f O´3 {2 2 τ{p 2 , ζ " p1`ZD N q{p1`Zq, and τ " gM pp p1´f N q{A. This can be solved forf N to find We find that for large mass bodies, k 1 ,k 2 Ñ 8,f N Ñ f N , so that planetary nitrogen abundance matches the primordial value. In the limit k 2 Ñ 0, this expression simplifies significantly tof N Ñ f N k 1 {p1`k 1 q. This expression allows us to derive the final nitrogen abundance as a function of planetesimal mass, by noting that g{A " Gρ To determine the planetary nitrogen abundance fraction that results from original accretion, we need the typical planetesimal size. For this, we use the isolation mass M iso " 1.31 [3]. In the limit that k 1 ! 1 and neglecting the dependence on planetary mass of D N , this giveŝ Interestingly, the dependence on stellar mass of this quantity is practically indistinguishable from that of the alternate nitrogen source, Eqn. (22).

Which Atmospheric Thresholds Are Important for Habitability?
Earth's atmosphere is quite comfortably above any catastrophic thresholds, being about two orders of magnitude larger than needed to prevent total atmospheric escape, maintain liquid surface water, and buffer diurnal temperature changes. However, given the exponential dependence on constants of some of these conditions, we investigate the influence each exerts on our multiverse calculations.

How Much Atmospheric Loss Occurs in Other Universes?
In this paper, we restrict our attention to terrestrial planets, which are defined such that light gases such as hydrogen and helium, but not heavy gases such as water, oxygen and nitrogen, undergo Jeans escape. For these planets, the dominant form of atmospheric escape is driven by stellar XUV light, and is in the energy limited regime (for recent reviews, see [61,62]). The mass loss rate for this type of escape is given by equating the energy of UV light absorbed by the atmosphere with the energy of atmospheric particles ejected at the escape speed [63], Here, is an unimportant efficiency factor. This is independent of atmospheric mass, being limited by the amount of energy imparted in the upper atmosphere rather than the amount of material present. In [64] it was estimated that an XUV flux greater than 60 times Earth's value would be needed to induce a catastrophic mass loss rate of 1.8ˆ10 9 g/s, capable of eroding the entire atmosphere. For reference, M dwarfs and young K dwarfs are subjected to 100-400 times Earth's XUV flux [85].
The total atmospheric mass loss through X-ray flux may be found through Equation (18), taking rotation evolution into account: Using the evolution dictated by Equation (20) and in the limit t ‹ " t brake , this integral can be performed to find The condition t ‹ " t brake , which holds by three orders of magnitude in our universe, is not necessarily generic; we compute t brake {t ‹ " 0.0019λ 5{4 α 1{2 { f open , which can be much larger than 1 if no stellar magnetic field lines are open for certain parameters. However, including a more complete expression does not affect the calculated probabilities appreciably, while considerably complicating the formulae. In Figure 1, we display the atmospheric mass loss for temperate, terrestrial planets as a function of stellar mass, for three different values of the fine structure constant. The difference resulting from adopting the two alternate origin scenarios is also displayed, but is seen to be minimal. This defines some stellar mass below which more than the initial atmosphere is lost through XUV irradiation, which depends on fundamental constants, and can be larger than the solar mass (λ " 1{1.8) in some regions of parameter space.

How Much Atmosphere Is Needed to Maintain Liquid Surface Water?
Liquid surface water can exist only when atmospheric pressure exceeds that at the triple point, where the three low energy phases of water coexist in equilibrium. The location of the triple point can be determined by noting that the solid-liquid transition is almost independent of pressure, and occurs at temperature set by the vibrational molecular energy T freeze " The liquid-gas transition is given by the Clausius-Clapeyron equation as PpTq " P 0 e´L {T , where the latent heat of evaporation is L " αr H 2 O . The coefficient P 0 can be found by enforcing that the phase curve terminates at the observed critical point of water of 647 K and 22.1 MPa. Though an imperfect description, the van der Waals equation of state may be used to provide a theoretical expectation for the location of the critical point in terms of the molecular radius and energy , yielding T crit " 8{27 and P crit " {p18πr 3 q [65]. Normalizing to fit our observed values, this yields the pressure at the triple point to be This can then be related to minimal atmospheric mass capable of supporting liquid water through M min " 4πR 2 C P triple {g, giving This is about 0.006M atm for our values.

How Much Atmosphere Is Needed to Buffer Diurnal Temperature Changes?
Earth's atmosphere retains substantial heat, which buffets the day-night temperature difference from the otherwise extreme variations that would occur, such as the day-night temperature differences on the Moon and Mars which can reach hundreds of degrees Kelvin. This occurs because the relaxation time of Earth's atmosphere, estimated as the ratio of thermal energy over the power supplied, t relax " E therm {Q solar , is about 100 days. This gives For small enough atmospheric mass, this is less than half a day, and the atmosphere does not play a significant role in averaging out daily variations of stellar flux. This occurs at the threshold The exact mass depends strongly on water content, as evidenced by the extreme temperature differences present in Earth deserts, but we do not consider this here.

Are Only Slowly Rotating Stars Habitable?
There is evidence from noble gas isotopes [66], the Moon [67], and Venus [68] that the Sun began as an anomalously slow rotator. However, it is not currently possible to determine precisely how slow, and many studies only differentiate between stars in the lower 25 percentile. If true, this suggests a selection effect: ordinarily rotating stars may be incapable of hosting life, presumably due to high early atmospheric loss.
To determine the compatibility of this habitability hypothesis with the multiverse, we follow the fraction of slowly rotating stars f slow " minˆM atm ∆M XUV , 1˙ (36) This treats the initial rotation distribution as uniform up to the natural value Ω 0 , which is loosely consistent with observations of stellar populations [69]. To account for the observation that the Sun appears to be in the lower 25 percentile, we rescale the fraction of slow rotators (of Sun-like stars) in our universe to be 1/4.

Is Atmospheric Stability a Factor Determining Our Presence in This Universe?
We can now test the various atmospheric habitability thresholds, on the basis of their compatibility with our observations within the multiverse. To this end, we test the following four thresholds: loss due to XUV radiation, the minimal mass for stable liquid surface water, the minimal mass to buffer diurnal temperature changes, and the notion that only stars which are slowly rotating are habitable. In addition, we check both the early and late origin scenarios for our atmosphere, both an independent and linearly dependent abundance as a function of stellar nitrogen content, and either restricting to a narrow range of Earth-like carbon-to-oxygen values, or not. In Table 1, we display the various Bayes factors for each of these combinations.
We find that when restricting consideration to a narrow range of carbon-to-oxygen ratios, the Bayes factors for the various habitability criteria do not vary significantly. When considering the carbon-to-oxygen ratio to not play a factor in habitability, however, several of the habitability criteria are severely disfavored in the multiverse context. The disfavored criteria all have to do with the assumption that planetary nitrogen content scales linearly with stellar system nitrogen abundance, and does not depend on the atmospheric source or threshold mass. This is a consequence of our universe being situated very close to a precipitous threshold where nitrogen-14 is unstable [6], which affects the probabilities if the carbon-to-oxygen ratio is unimportant but does not if restricted to the subspace where the carbon-to-oxygen ratio is close to our observed value. We note that in [6] we found additional reasons to favor a restricted range of carbon-to-oxygen ratio based on the observed Hoyle energy value and organic to rock ratio in our universe. Apart from this insight, no strong preference can be given to the different atmospheric origin scenarios, threshold masses, or expectation on whether only slow rotators are habitable. Our conclusion is that it is certainly consistent that atmospheric mass may play a large role in the habitability of our universe, but it does not appear to be a driving factor in determining our particular observations. Table 1. Bayes factors for various atmospheric habitability criteria relative to the baseline case where atmosphere mass is unimportant for habitability. Small values indicate that a set of assumptions is disfavored to a corresponding degree in the multiverse framework. The cases considered are that atmosphere must be large enough to withstand XUV loss, be above the triple point of water, can buffer diurnal temperature changes, and that only slowly rotating stars are habitable. The late delivery vs. initial columns consider both potential origins of the atmosphere, and the N dep columns consider that planetary nitrogen abundance scales linearly with stellar system abundance.

Are Planetary Magnetic Fields Generic?
A planet's magnetic field is purported to be essential for habitability, as it shields against charged particles, preventing stellar wind stripping (see, e.g., [70]). However, it must be pointed out that magnetic fields also provide several avenues for ion escape [71], which may well represent the dominant form of atmospheric loss on Earth today [72]. Indeed, Venus has managed to retain its atmosphere without an intrinsic (as opposed to induced by the Sun's) magnetic field, despite being closer to the Sun than Earth.
Given the uncertain importance of planetary magnetic fields for habitability, we ask whether their properties change significantly in other universes, and thus whether demanding their presence influences the probabilities of any of our observables. We focus on five relevant aspects required for a magnetic field to be both present and protective: (i) The core's magnetic Reynolds number is large enough to support a dynamo. (ii) The magnetosphere must extend beyond the atmosphere, as otherwise it will have little effect on loss properties. (iii) The star's temperate zone must be outside its Alfvén zone, as otherwise the planetary and stellar magnetic field lines connect, forming a direct line of transport which dumps stellar wind onto the planet's poles, rather than act as a shield. (iv) The development of a magnetic field requires a metal core, placing limits on the oxygen content of the planet. (v) The magnetic field is generated through a dynamo, and so requires the core to remain at least partly liquid for an appreciable duration. If planetary magnetic fields are essential for habitability, all of these conditions must be met.

When Is the Magnetic Reynolds Number Large Enough to Induce a Dynamo?
Both theory and simulations of the Earth's core indicate that a dynamo will only exist when advection of the magnetic field dominates over diffusion [73]. This can be summarized as a condition on the magnetic Reynolds number R a " v core L{η ą 10-100 (Earth's magnetic Reynolds number is about 10 3 ) [74]. This condition can be used to place constraints on the fundamental constants, using the length scale L " R core , and magnetic diffusivity η " 1{p4πσ electric q with σ electric the electrical conductivity, which is related to thermal conductivity through the Wiedemann-Franz law [75] In [4], we found an expression for the thermal diffusivity in terms of fundamental constants p q, which is related to thermal conductivity through κ heat "κ heat {pc p ρ rock q. The core convective speed can be obtained from mixing length theory, v core " pLq{pρ rock H T qq 1{3 [76]. Using our expression for heat flux q " 0.58α 11{2 m 5 e {M pl from [4] and the generic expression for scale height H T " c 2 s {g, we find These scalings are not significantly altered if we instead use the magnetostrophic estimate for the core convection velocity, also from [76].

Is the Magnetosphere Always Larger Than the Atmosphere?
In order for a planetary magnetic field to be an effective shield against stellar wind, it must extend beyond the atmosphere. The size of the magnetosphere can be estimated as the point at which the magnetic pressure is equal to the stellar wind pressure, yielding for a dipole field [77] the standoff distance: To evaluate this, we use the expressions for density and speed of solar wind from Section 2. It remains to estimate B 0 , the strength of the magnetic field at the planet's surface.
There are an inordinate number of proposals for how planetary magnetic field strength depends on planetary characteristics, as reviewed in [78]. We use the Elasser number rule, which posits that the Lorentz force and Coriolis force are roughly equal, and results in Here Ω is the planet's angular rotation speed. For terrestrial planets, the atmospheric scale height is much smaller than planetary radius, and so it suffices to compare the magnetosphere size to the latter. Using our expressions above, and defining Y " pR core {R planet q 3 , we find this to be This ratio evaluates to 10 for our values of the constants and Earth's core radius. The dependence on fundamental constants is rather weak, and so it takes a drastic change to alter the conclusion that the magnetosphere extends beyond the atmosphere.

When Is the Temperate Zone Outside the Alfvén Zone?
If a planet is orbiting inside its host star's Alfvén zone, its intrinsic magnetic field lines will connect to the star's, which will result in a highly increased level of bombardment by charged particles. This is expected to be the case for the inner planets in the Trappist-1 system, for instance [79] from simulations. Given our expressions for both the temperate zone and Alfvén radius from Section 2, it is straightforward to derive their ratio: For fixed physical constants, this defines a smallest stellar mass for which this condition holds. In our universe this is about .1 M @ , in accordance with the expectation that Proxima Centauri b, which orbits a 0.12 M @ star at 0.05 AU, is outside the Alfvén zone for the most likely values inferred for its orbital parameters [80]. Our treatment ignores the nonsphericity and nonstationarity of the Alfvén zone and potential planetary eccentricity, which may cause the orbit to periodically dip into the Alfvén zone throughout the year.
Note that the dependence on starspot fraction is of crucial importance in this expression, as otherwise this threshold stellar mass would be smaller than the smallest stellar mass. As such, this condition is loosely coincident with the onset of a full stellar convection zone.

When Does a Core Stratify Geochemically?
In [81], it was pointed out that if a planet's mantle oxygen content is too high, the iron will all be in the form of iron oxide (FeO), and will not differentiate to form a core. They find that the quantity R 1 = (Mg + 2Si + Fe)/O must exceed 1 in order for a core to develop, Earth's value of this ratio being 1.23. Interestingly, about 4% of the Earth's oxygen is left over after binding with magnesium and silicon, so that only 86% of Earth's iron makes it into the core. This raises the additional possibility that if a planet's oxygen is depleted before its magnesium and silicon are consumed, no iron will be left in the mantle or crust. This could have an additional adverse effect on habitability, which would restrict the allowable oxygen content required for habitability to a narrow range, but we leave exploration of this for future work. It was argued in [82] that planets with core mass fraction below "0.24 would have much higher rates of volatile subduction, due to more extensive volcanism, thicker crust, and stabilized amphibole group. This places a potential lower limit to the allowable core mass for habitability.
Though the core development condition depends on the ratio R 1 above, this depends on the abundances of both the alpha elements and iron, which are set by two different supernova processes, and so will scale differently with fundamental constants. In [6], we found the dependence of the alpha element abundances (C, O, Mg, and Si) on the Hoyle resonance energy E R " .626pm u`md q`p0.58α´0.0042qm p , with m u , m d the masses of the up and down quarks, as found in [18]. We also found an expression for the metal to rock ratio, from which we may determine the quantity For Earth, this value is R 2 " 0.163. This assumes a linear relationship between stellar and planetary metal to rock ratios, which indeed is found [83].
In Figure 2, we plot the oxygenation ratio for various values of the metal ratio R 2 . It can be seen that with R 2 held fixed at the observed value, the oxygenation ratio is less than 1 for ∆E R ą 3.6 keV. While we are rather close to a potential anthropic boundary with metal fraction held fixed, allowing it to vary relaxes this closeness. In fact, there is a silicon and magnesium rich region of parameter space for larger values of ∆E R which also satisfy the R 1 ą 1 requirement. Above a metal fraction of 0.62, these two branches merge, and planets will always contain enough iron to form a core. As discussed in [6], such metal rich planets may be unsuitable for life for reasons other than the possession of a magnetic field, but we found that placing an upper bound on the metal content does not appreciably affect the probabilities we compute, and we do not concern ourselves with such a boundary here.

What Sets the Core Solidification Timescale?
The presence of a planetary dynamo requires a liquid convective core, which cannot be sustained indefinitely. As heat leaks from the planet, an initially liquid core will cool and solidify. If the solidification timescale is too rapid, any magnetic field will cease before life can take hold on a planet, and so one important consideration is the longevity of a liquid core.
First, we must establish that terrestrial planets possess enough heat for their cores to initially be liquid. This follows almost from our definition of a terrestrial planet, which demands that the gravitational binding energy is of the same order of magnitude as molecular binding energies, so that chemical reactions may take place on the planet's surface. Given the increased temperature and pressure of the planetary interior, the melting point will naturally be exceeded in the core.
The solidification timescale can be simply estimated as t solid " E core {Q heat , where E core is the energy required to be leached from the core for solidification to take place, and Q heat is the total core power. A proper estimate of E core would take into account the difference between the gravitational binding energy and the energy that would result in solidification; thankfully, however, these two energies are similar in magnitude, another consequence of restricting our attention to terrestrial planets. So, we may approximate the total energy in the core as E core " GM 2 core {R core . By the same token, Q heat has components due to formation and crystallization, which are roughly equal. In [4], we found that Q heat " GM planet ρ rock κ heat based on dimensional analysis. There, we also consider radiogenic heat and time dependence in more detail, which we neglect here. This may indeed be important; as discussed in [84], too much radioactive heating can prevent core convection. However, we do not consider this in detail here.
With this, the core solidification timescale is very simple: If a long-lived liquid outer core is necessary for habitability, this timescale must be larger than some timescale typical for the development of complex life, which we take here to be proportional to the stellar lifetime (see [2] for an exploration of different choices on this matter). We normalize this time to the expectation that the outer core will remain liquid for another 700 Myr from [85].
An alternative view is that a solid inner core is actually necessary for the sustenance of a magnetic field, in spite of geologic evidence to the contrary (see [86] for zircon evidence of a magnetic field at 3-4 Ga). The inner core may have developed as late as 565 Mya, based on magnetic evidence from Ediacaran rocks that record an anomalously low field strength, taken to signal a rearrangement in field configuration indicative of a recently established solid inner core [87]. This apparent incompatibility is reconciled if another mechanism generated the magnetic field before core solidification, as for example a long lived liquid mantle ocean [88]). In this case, the above timescale would need to be comparable to the evolutionary timescale, rather than simply longer than it.

Is a Planetary Magnetic Field Necessary for Habitability?
To treat intrinsic planetary magnetic fields as essential for habitability, we include the product of all five factors into the habitability condition as If we incorporate this into our calculation, we find that the Bayes factor relative to the base case where magnetic fields are not taken to be important is B " 1.52. We also probe the relative importance of each of these subconditions in Table 2 by first only incorporating each condition in isolation, and then incorporating the four others without each condition, into the calculation. Of the five factors considered, the magnetic Reynolds number, magnetosphere radius, and Alfvén zone conditions do not perceptibly alter the probabilities. The core existence condition slightly decreases the probabilities, while the core timescale condition slightly increases them. So, the notion that a magnetic field is necessary for habitability is compatible with the multiverse, and although it is even slightly preferred to the base case, the difference is not statistically meaningful enough to draw the conclusion that the converse hypothesis is disfavored.
We also remark that the base case here took the carbon-to-oxygen ratio to be important for habitability. If instead we drop this assumption, we find the Bayes factor is B " 0.050. The driving factor in making this so low is the core solidification timescale, as can be seen in Table  2. Therefore, we find that the assumption that planetary magnetic fields are important is only compatible with the multiverse if carbon-to-oxygen ratio is also important. This echoes the results of Section 3 and [6], where we found that restricting the carbon-to-oxygen ratio was important for compatibility with the multiverse on other accounts. This also suggests a test of the multiverse hypothesis, for if we find that complex life occurs only on magnetized planets but independently of carbon-to-oxygen ratio, our presence in this universe would be quite unlikely.

Discussion
Though few agree on exactly what conditions are required for habitability, it surely depends on the confluence of a great many factors. Likewise, our notions of habitability strongly affect the expectation for the distribution of life, both throughout our universe, and in others. Because of this, very fine-grained effects have the potential to radically alter our estimations of the probability of our existence in this particular universe, and our observations in general. This places us in a scenario where the importance of all discussed habitability factors must be tested before we may make any statements about multiverse probabilities with a relatively high degree of certainty. The stellar activity and atmospheric aspect of this program was undertaken in this paper. Uncertainties abound: the physics dictating the corona, stellar wind and flares, the relative importance of different atmospheric erosion rates, the ultimate source of Earth's atmosphere, the importance of planetary magnetic fields, and the distribution of all these quantities across different stellar systems are only now coming to light. While we have tried to hedge our ignorance in as many aspects as possible by contemplating competing accounts of these effects, we have necessarily restricted our attention in certain cases, and completely neglected other potentially important effects. Thus, while our work cannot claim to be a definitive exploration of stellar activity and atmospheric effects in other universes, it does represent an important first step.
Perhaps the biggest takeaway of our findings is that, if one believes that a relatively narrow carbon-to-oxygen ratio is required for complex life (as may be argued by the vastly different tectonic regimes that occur outside the interval (0.5,1)), any atmospheric habitability condition we considered had no significant bearing on multiverse probabilities. In this light, all that can be said is that atmospheric presence and stability does not appear to be a major determining factor for why we are in this universe. This is plausible, since the Earth's atmosphere is about two orders of magnitude larger than any threshold we are aware of, but many effects we consider depend exponentially on fundamental constants, so this conclusion is by no means automatic.
On the other hand, if we entertain the possibility that a carbon-to-oxygen ratio relatively close to ours is not required for habitability, altogether different conclusions are drawn. We are forced to conclude, under this assumption, that planetary magnetic fields cannot be important for life, because it renders many of the otherwise less likely regions of parameter space infertile, making us outliers. Additionally, when treating the carbon-to-oxygen ratio as unimportant, we find planetary atmospheric nitrogen must not scale with stellar system nitrogen abundance, or our presence in this universe is unlikely, independent of uncertainties about atmospheric source and lower atmospheric mass threshold.
Both of these findings also suggest potential methods for testing the multiverse hypothesis, if the true habitability conditions turn out to be incompatible with these expectations. So, if we find that an Earth-like C/O is not needed for complex life and either that atmosphere mass scales with stellar nitrogen or that planetary magnetic fields are required for life, the predictions the multiverse framework has made will be found to be incorrect. While these tests may be rather far off, the salient point is that they are possible in principle. Various biosignatures have already been proposed to help determine the distribution of life throughout the universe, for several places inside and out of our solar system, including searching for relic biomarker compounds on Mars [89], abundance ratios of organic compounds on icy moons such as Enceladus [90,91], chemical disequilibrium and even microbial absorption in Venus's atmosphere [92], and atmospheric gases such as oxygen around exoplanets [93]. In fact, it is conceivable that the next few generations of experiments will be able to measure biosignatures on exoplanet populations large enough to distinguish trends with respect to system parameters such as composition [94], that the presence of planetary magnetic fields can be measured through auroral emissions [95], and that the relation between planetary atmospheric size and stellar composition can be determined [96].