Arxiv:2103.05449V1 [Astro-Ph.EP] 9 Mar 2021

Draft version September 3, 2021 Typeset using LATEX twocolumn style in AASTeX63

Assessing the Formation of Solid Hydrogen Objects in Starless Molecular Cloud Cores

W. Garrett Levine 1 and Gregory Laughlin 1

1Yale University 52 Hillhouse, New Haven, CT 06511, USA

ABSTRACT The properties of the first-discovered interstellar object (ISO), 1I/2017 (‘Oumuamua), differ from both Solar System asteroids and comets, casting doubt on a protoplanetary disk origin. In this study, we investigate the possibility that it formed with a substantial H2 ice component in the starless core of a giant molecular cloud. While interstellar solid hydrogen has yet to be detected, this constituent would explain a number of the ISO’s properties. We consider the relevant processes required to build decameter-sized, solid hydrogen bodies and assess the plausibility of growth in various size regimes. Via an energy balance argument, we find the most severe barrier to formation is the extremely low temperature required for the favorability of molecular hydrogen ice. However, if deposition occurs, we find that the turbulence within starless cores is conducive for growth into kilometer-sized bodies on sufficiently short timescales. Then, we analyze mass loss in the interstellar medium and determine the necessary size for a hydrogen object to survive a journey to the Solar System as a function of ISO age. Finally, we discuss the implications if the H2 explanation is correct, and we assess the future prospects of ISO science. If hydrogen ice ISOs do exist, our hypothesized formation pathway would require a small population of porous, 100 µm dust in a starless core region that has cooled to 2.8 K via adiabatic expansion of the surrounding gas and excellent shielding from electromagnetic radiation and cosmic rays. 1. INTRODUCTION ample, Bialy & Loeb(2018) propose an unusually-thin The observed properties of the interstellar object minor axis, whereas others investigate ultra-porous ob- −5 3 (ISO) 1I/2017 U1 (‘Oumuamua) do not fit neatly within jects with bulk density ρ ∼ 10 g/cm . Specifically, the known astronomical taxonomy. In the Solar Sys- Moro-Mart´ın(2019) considers coagulation of icy grains tem, minor bodies are usually classified as asteroids or in a protoplanetary disk, and Luu et al.(2020) examines comets, with the latter often experiencing force from diffusion-limited aggregation in an exo-comet’s tail. solar-driven sublimation of volatiles. Optical and IR Radiation pressure explanations generally assume im- images of ‘Oumuamua found no coma (Meech et al. plicitly that ‘Oumuamua’s origin was an exoplane- 2017; Jewitt et al. 2017; Trilling et al. 2018), yet analysis tary system: a “disk product” hypothesis. However, of its trajectory revealed non-gravitational acceleration this birthplace has difficulties with explaining either (Micheli et al. 2018). Additionally, dramatic light-curve ‘Oumuamua’s kinematics or the inferred number den- variations implied that its shape was approximately that sity as representative of a Galactic reservoir of ISOs. of a 6:6:1 oblate ellipsoid (Mashchenko 2019). Together, First, an inbound velocity near the Local Standard of Rest (LSR) is consistent with an object age of less than

arXiv:2103.05449v1 [astro-ph.EP] 9 Mar 2021 these attributes suggest that ‘Oumuamua is a novel type of object and distinct from both Solar System comets as 100 Myr (Mamajek 2017; Hallatt & Wiegert 2020). If well as the extrasolar comet 2I/Borisov (Jewitt & Luu disk products form consistently throughout Galactic his- 2019). tory, ‘Oumuamua would have coincidentally been se- To reconcile ‘Oumuamua’s trajectory with a non- lected from the population’s youngest 1%. Further- volatile composition, some studies have suggested ra- more, Pan-STARRS’ detection of any ISO within 0.2 diation pressure as the non-gravitational force. For ex- AU of Earth implies number densities such that the average star must eject at least 20 M⊕ of material in disk- product hypotheses (Do et al. 2018). This large produc- Corresponding author: W. Garrett Levine tivity is in marginal tension with planetary formation [email protected] models. 2 Levine and Laughlin

As an alternative to radiation pressure, Seligman & the H2 solid-vapor equilibrium and producing objects Laughlin(2020) propose that outgassing of molecular with a variety of characteristic sizes. In this study, we hydrogen acquired in a giant molecular cloud (GMC) focus on the feasibility of producing ‘Oumuamua-like ob- could account for ‘Oumuamua’s non-gravitational ac- jects in starless cores. celeration and non-detection of common volatiles. Addi- Observationally, Sandford et al.(1993) attributed a tionally, it explains the inbound kinematics and extreme detected 2.415 µm feature in the direction of ρ-Oph to aspect ratio as consequences of erosion during the finite solid H2 embedded in interstellar ices. However, the lifetime of solid H2 in the interstellar medium (ISM) result is unconfirmed and has recently been noted as (Domokos et al. 2017). Because GMCs are a larger mass an upper bound on concentration (Boogert et al. 2015). source than exoplanetary systems, number density con- Our study here, however, considers bulk object forma- straints are easily satisfied. Finally, a “cloud product” tion as opposed to processed ices. would corroborate ‘Oumuamua’s strong dynamical cor- We divide hypothetical cloud product formation into relation with the surrounding Columba and Carina stel- distinct size regimes, then consider the processes affect- lar associations (Gaidos 2018; Hallatt & Wiegert 2020), ing growth in each one. In Section2, we consider the the latter of which may be as young as 25 Myr as deter- starless core energy balance and the prospects for solid mined by recent spectroscopic evidence (Schneider et al. hydrogen deposition on dust. Then, Section3 constrains 2019). Furthermore, the cloud product hypothesis im- the rates of H2 and dust accretion onto an isolated ob- plies an origin for the ISO that was independent of a ject. Next, Section4 analyzes systems with multiple particular host star. aggregating “hailstones,” which may form planetesimal- Because H2 is the only proposed composition that can sized, snowball-like ISOs. Section5 describes erosion in explain all of ‘Oumuamua’s properties as generic to a the ISM, which determines the object’s metallicity upon population of objects, it is imperative that potential encountering the Solar System. Finally, Sections6 and formation pathways be closely-examined. This study 7 address the ramifications of ‘Oumuamua as a cloud evaluates the plausibility of H2 ice ISOs originating product. Because of the stringent requirements for de- in GMCs, a possibility raised soon after ‘Oumuamua’s position, confirmation of H2 in future ISOs would probe detection by Füglistaler& Pfenniger(2018). Specifi- previously-undetected interstellar conditions. cally, we focus on regions that are unproductive in star- formation since these will host the coldest temperatures. 2. SOLID HYDROGEN DEPOSITION These comparatively-quiescent GMCs, in contrast to the solar birthplace, host numerous starless cores (Vázquez- 2.1. Thermodynamic Requirements for Solid Hydrogen Semadeni et al. 2005). Figure1 shows the phase diagram for pure H 2. If Solid hydrogen has previously been examined as a ‘Oumuamua-like objects contain molecular hydrogen, possible source of dark matter (Pfenniger & Combes the outgassing material could only be acquired where 1994; White 1996). However, current ΛCDM cosmol- the solid phase is favored, motivating our scrutiny of ogy requires a substantial non-baryonic component to the coldest, densest interstellar environments: starless the Universe’s matter density in addition to any remain- cores in GMCs. Although the presence of He may mod- ing unidentified hydrogen (Bergström 2000). Assuming ify the equilibrium lines, we assume the effect is neg- ISO number density consistent with Do et al.(2018), ligible because of helium’s inert nature. In support of a population of ‘Oumuamua-like objects would not be this premise, Füglistaler& Pfenniger(2016) consider de- dynamically important for the Galaxy. Thus, the chan- position in two-component fluids and find qualitatively nel of H2 production that we investigate would have no similar phenomena compared to a pure system. immediate ramifications for the topic of dark matter. Nonetheless, the approximately exponential temper- Theoretically, Füglistaler& Pfenniger(2018) analyze ature dependence of solid-vapor equilibrium (Walker gravitationally-unstable fluids near solid-vapor equilib- 2013) means that slight deviations from Figure1 may rium, arguing that starless cores may spawn objects critically affect the feasibility of deposition. While pre- ranging from small condensates via deposition to multi- vious studies of H2-He mixtures have often considered kilometer ISOs via coagulation. They find scaling laws conditions analogous to giant planet interiors (Militzer at high pressure using molecular dynamics simulation, 2005), laboratory data for ISM-like temperature and then extrapolate to assess GMCs. Similarly, Walker & pressure will be necessary if hydrogen cloud products Wardle(2019) consider self-gravitating interstellar gas are found in the future. clouds, showing that heat transport via convection may be conducive to lowering ambient temperatures beneath 2.2. Extreme Conditions in Starless Cores Solid Hydrogen Objects from Starless Cores 3

Table 1. Summary of assumed conditions in an extremely cold, dark, dense starless core.

Symbol Quantity Value Source

Tsc Gas Temperature 2.8 K Kong et al.(2016) 6 -3 nH2 Number Density of H2 3.0 × 10 cm Ivlev et al.(2015) −2 φd Dust/Gas Mass Ratio 1.99 × 10 Draine(2011)

fHe Helium Number Fraction 0.2 Draine(2011)

Av Visual Extinction 200 mag Draine(2011) 3 Rsc Starless Core Radius 5.0 × 10 AU Weidenschilling & Ruzmaikina(1994) 4 cs Sound Speed 1.0 × 10 cm/s Weidenschilling & Ruzmaikina(1994) 6 τsc Starless Core Lifetime 3 × 10 yr Kong et al.(2016) 7 τcl GMC Lifetime 3 × 10 yr Kong et al.(2016)

ments other than H and He are in the solid phase. While much ISM hydrogen is atomic, catalysis on dust grains leads to the molecular form dominating in cores (Hol- lenbach & Salpeter 1971). Importantly, our assumed temperature favors the solid phase of molecular hydrogen, as indicated by the phase diagram in Figure1. Sustained Tsc this close to that of the cosmic microwave background (CMB) has yet to be seen, but excitation temperatures around 3.5 K in M31 and 4.0 K in our Galaxy have been reported (Loinard & Allen 1998; Kong et al. 2016). GMC observations continue to find density structure at the smallest resolvable scales; thus, the coldest, densest regions have yet to be constrained (Vázquez-Semadeniet al. 2005). Because observations have yet to detect temperatures below the Figure 1. Phase diagram for H2 (Air Liquide 1976; Ander- son 1989), with the solid circular mark corresponding to the solid-vapor equilibrium line in starless cores, we must fiducial starless core from Table1. instead examine the most frigid conditions that these regions could theoretically support. GMC morphology is often governed by the turbulence that is ubiquitous in the ISM (Vázquez-Semadeniet al. 2.3. Dust Grain Energy Balance in Starless Cores 2005), so density and pressure within a given cloud can Even if thermodynamics supports the formation and vary by orders-of-magnitude. Typical 30 Myr GMC life- persistence of solid H2, there must be surfaces on which times are too short for conditions to reach thermody- deposition begins. Dust is the most viable candidate, namic equilibrium (Shu et al. 1987), and the gas is often but its temperature Tgr will not necessarily equilibrate described as “clumpy” (Draine 2011). Because starless to that of the surrounding gas. In starless cores, grains cores are often colder and better-shielded from electro- experience the following energy sources (Draine 2011): magnetic radiation than the surrounding cloud, they are the optimal locations to investigate solid hydrogen de- 1. Any power from stars in the Galaxy – the inter- position (Füglistaler& Pfenniger 2018). stellar radiation field (ISRF) – at wavelengths for Table1 details our fiducial assumed conditions, and which the surrounding GMC is optically-thin. the corresponding point is marked in Figure1. Nom- inally, these environments would be unstable to the 2. Blackbody CMB radiation with Tcmb = 2.725 K. canonical Jeans’ gravitational collapse, but we assume 3. Infrared (IR) photons emitted by warm dust in the that magnetic field pressure supports the core for a outer cloud layers. few 106 yr before re-expansion (Vázquez-Semadeni et al. 2005; Kong et al. 2016). We equate the dust/gas mass 4. Inelastic gas-grain collisions working to equilibrate ratio to the cloud metallicity, which implies that all ele- Tgr and Tsc. 4 Levine and Laughlin

5. Bombardment by relativistic cosmic rays (CRs). While this ISRF prohibits solid H2 deposition, most incident radiation is eliminated in the outer cloud lay- 6. Heat of decay from short-lived radionuclides ers. Hoang & Loeb(2020) consider attenuation of all (SLRs) delivered from supernovae. incident radiation by a single scaling factor, but we use wavelength-dependent extinction A from Mathis et al. Broadly, we classify these sources into incident radia- λ (1983). Specifically, the quantities and are less tion (items 1-3) and other mechanisms (items 4-6). We a e than unity for MRN-sized grains and given by Draine assume there are no field stars providing anisotropic (2011) for silicate grains as flux, consistent with examining clouds with low star- formation rates. Furthermore, because the only heat h λ i−2h a i ≈ 1.4 × 10−3 . (2) sources for the gas are the ISRF, CRs, and collisions 100 µm 0.1 µm with dust, T drives the overall starless core energy bal- gr Thus, the grain temperature from Equation1 is also ance at our assumed number density (Goldsmith 2001). size-dependent. With Equation2, we can compute Therefore, constraining the dust temperature is imper- the average efficiency over the Planck function to find ative to assessing the formation of solid H . 2 (Draine 2011) In general, these energy fluxes can vary with position in the Galaxy. Nonetheless, a young inferred ISO age h a ih T i2 ≈ 1.3 × 10−6 , (3) and a low velocity versus the LSR indicate that ‘Oumua- pl 0.1 µm 1 K mua has likely experienced interstellar conditions similar to those surrounding the Solar System today. Thus, which is e for ideal particles. With these efficiencies, the magnitudes of the incident power sources are better we plot the resulting Tgr for varying grain sizes on the constrained than they would be for an old object. solid lines in Figure2. MRN dust is too warm, so we display a & 1000 µm where deposition may be feasible. 2.4. Photon Radiation Bath for Ideal Grains In addition, we must consider the CMB with Tcmb = 2.725 K. Attenuation is negligible in starless cores at We first consider ideal grains characterized by peak wavelength λ ∼ 1000 µm. We plot the modified the Mathis-Rumpl-Nordsieck (MRN) distribution, with cmb T with the combined flux on the dashed lines of Figures dn/da ∝ a−3.5 describing dust size a = 0.005 − 0.25 µm gr 2. As expected, the CMB has a pronounced effect on the (Mathis et al. 1977). Most of the mass resides in dust energy balance only as grain sizes approach λ . large grains, while nearly all of the surface area is cmb Finally, we consider another source of photons: warm in small grains (Draine 2011). Furthermore, we con- dust from the outer GMC layers. Grains located at lower sider spherical, non-porous particles with bulk density A are responsible for optical darkness in starless cores, ρ = 3 g/cm3. While both silicate and graphite dust v d and they re-emit energy primarily as IR. The spectrum should reside in starless cores, the former will maintain and strength of this radiation depend on cloud geometry a slightly lower temperature and be more permissive of and other considerations. For example, Mathis et al. deposition (Mathis et al. 1983). With φ from Table d (1983) superimposes two blackbody fields from graphite 1, we find number density n = 1 × 10−6 cm−3 for the d and silicate grains. Similarly, we use fields with T = largest grains of mass m = 2 × 10−13 g. 1 d 16 K and T = 7 K and dilution factors W = 1 × 10−3 At steady-state, the fluxes sum to zero: P dE/dt = 0. 2 1 and W = 6 × 10−3, with our deviations from Mathis For an isolated grain in a radiation bath, we equate the 2 et al.(1983) attributed to galactocentric distance. absorbed (dE /dt) and emitted (dE /dt) power as a e The effect of warm dust flux added to the ISRF and dE −dE CMB on the energy balance of starless core grains is a = f A = σT 4 A = e , (1) dt r g a gr s e dt shown in the dotted lines of Figure2. Deposition cannot occur on ideal grains in our fiducial starless core’s where f is the incident radiation flux, A and A r g s radiation bath. are the cross-section and surface areas of the particle, a and e are the absorption and emission efficiencies, 2.5. Other Processes’ Effects on Ideal Grains and σ is the Stefan-Boltzmann constant. We assume the ISRF at galactocentric distance 10 kpc from Mathis Next, we consider inelastic gas-grain collisions reduc- −2 2 ing the difference between Tgr and Tsc (Burke & Hollen- et al.(1983), with fr = fISRF = 2.7 × 10 erg/cm /s over all wavelengths. Without considering other power bach 1983). The effect is quantified by Draine(2011) sources such as the CMB, a perfectly-efficient blackbody as

(a = e = 1) in only the Galactic ISRF would have tem- 1/2 dEgas X h8kTsc i perature T = 3.3 K. = n A α 2k∆T, (4) gr dt i πm g i i i Solid Hydrogen Objects from Starless Cores 5 where the sum is over all i gaseous species with masses mi, and ∆T = Tsc − Tgr. Here, αi quantifies the degree of inelasticity, the accommodation coefficient of Burke & Hollenbach(1983). For sticking fraction of the ith species si and the rest of the collisions elastic, we find αi = si; otherwise, energy conservation requires that αi > si. We consider the gas to contain only H2 and He since other molecules freeze around Av & 25 (Mathis et al. 1983; Goldsmith 2001). With αi = α = 0.5 (Draine 2011), the H2/He ratio from Table1, and spherical grains, we simplify Equation4 to

" # dE h n i gas H2 2 p 3 = 1.16√ a α 32πk Tsc ∆T. (5) dt mH2

For MRN-sized dust, inelastic collisions do not effec- tively regulate Tgr for ∆T < 10 K. Furthermore, the a2 dependence of both Equations1 and5 in the ideal blackbody limit implies that gas-grain interactions are negligible for all particle sizes. A violent means of kinetic energy transfer comes from cosmic ray (CR) bombardment, which heats the ambient gas and strikes the grains. White(1996) considers heating of ISM dust and assumes fCR = 1.6 × 10−3 erg/cm2/s. Because the GMC attenuates some low-energy CRs, this estimate is an upper bound. Since fCR is smaller than the power delivered by the photon bath, it can be neglected in the energy balance. Finally, we consider short-lived radionuclides (SLRs) synthesized in supernovae. For example, approximately 8 × 10−6 of aluminum is Al-26, with half-life 7 × 105 yr and heat of decay 4 MeV (Gounelle & Mei- bom 2008). Given cosmic Al abundance of 3 × 10−6 Figure 2. Temperature of spherical ideal grains in a photon (Draine 2011), we find that the typical grain receives bath for various sources of incident radiation. flux 1.5×10−5 erg/cm2/s. Because clouds with low star- formation rates are not enhanced in radioactive mate- Non-ideal particles might populate starless cores rial, we can safely ignore SLRs. through either in-situ coagulation or delivery from an Both kinetic and radiogenic power are inconsequen- exogenous source. In particular, subsonic turbulence imparts relative velocities that usually do not shatter tial, so the photon bath drives Tgr in the ideal regime. Thus, the energy balance imposes a severe constraint for colliding dust aggregates, so dense cores naturally host solid hydrogen deposition on MRN-like dust. grains larger than the MRN upper bound appropriate for other ISM dust (Ormel et al. 2009). For example, 2.6. Formation of Non-Ideal Grains numerical models from Ormel et al.(2011) and Wei- Thus far, our underlying assumption has been ideal, denschilling & Ruzmaikina(1994) find that clouds may 6 sub-micron grains. While the MRN distribution fits in- produce 100 µm conglomerates within a few 10 yr. terstellar extinction well, it does not preclude the ex- Other studies have analyzed enhanced coagulation istence of larger, porous, non-spherical dust. Instead from electrostatic effects (Spitzer Jr 1941; Ivlev et al. of “reddening” the ISRF, sub-millimeter particles are 2015). Because the rates of positive charging via the in the regime of geometric optics, extinguish all vi- photoelectric effect and negative charging via the ambi- 2 sual wavelengths equally, and appear spectroscopically ent cold plasma both scale as a , the mean grain charge “grey,” as understood via Equation2(Draine 2011). is independent of particle size and near zero for a range 6 Levine and Laughlin of plausible conditions. However, the random-walk na- be lower than in Figure2. In addition, a porous, fracture of these charging pathways means that individual tal configuration may increase the area radiating power grains are charged; this condition is optimal for coagu- while holding the cross-section intercepting flux close to lation (Ivlev et al. 2015). Oppositely-charged dust has constant. increased collisional cross-sections from monopole effects Furthermore, conglomerates acquire impurities be- at large separation due to net charge (Spitzer Jr 1941) cause dust does not necessarily coagulate with grains of and dipole effects at shorter distances due to induced similar composition. With all elements except H and He charge distributions on conductive grains (Draine 2011). in the solid phase, the resulting imperfections may per- Another potential source of non-ideal particles could mit more vibrational modes for the grain than an ideal- be delivery from the planetary nebulae of AGB stars. ized sphere. Furthermore, some of the impurity particles Specifically, Sahai et al.(2017) detected continuum may efficiently emit IR radiation (Werner & Salpeter millimeter emission in the Boomerang Nebula, which 1969). It must be noted, however, that grains which they attribute to sub-centimeter grains with total mass support H2 ice due to inhomogeneities would have diffi- −4 5 × 10 M . Expulsion into the ISM may fracture par- culty retaining these irregularities at macroscopic size. ticles, although conglomerate survival is unconstrained. Additionally, it’s reasonable to assume that icy star- Coincidentally, Sahai et al.(2017) also detect 1 K excita- less core grains may have non-zero albedos and reflect tion temperature in the adiabatically-expanding nebula, more incident radiation than their soot-like ISM coun- although the rapid outflow is unlikely to support depo- terparts. If particles hold sufficient solid hydrogen, the sition. object albedo will be determined by that of H2 ice. Characterization of a diffuse reservoir of coagulated Each of the aforementioned imperfections lowers Tgr. interstellar particles is challenging, but Spitzer has pro- The largest coagulated grains (a & 100 µm) may be close vided indirect evidence via GMC observations. A sys- to the deposition point even as spherical, non-porous tem of aggregates has smaller cross-sectional area versus blackbodies, so it’s possible that some of the largest ag- free grains, so the dust extinction per hydrogen nucleon gregates can support the phase transition. While depo- decreases with coagulation (Kim & Martin 1996; Ormel sition can be theoretically ruled-out in the ideal regime, et al. 2011). Additionally, limited direct measurements the observational fit of a solid hydrogen component for of presolar grains also have indicated that coagulated ag- ‘Oumuamua warrants further investigation into the con- gregates must exist. For example, Ulysses and Stardust figuration of sub-millimeter, non-ideal dust in starless both recovered porous, micron-sized aggregates (Land- cores. graf et al. 2000; Westphal et al. 2014). While the assumption of idealized MRN grains pro- 2.8. Spacing of Deposition Sites within Starless Cores vides a good model for diffuse interstellar space, starless If hydrogen deposition does occur, the phase transi- cores may host conglomerates with more complicated tion could only begin on non-ideal grain surfaces. Al- geometries and compositions. Thus, we must assess the though the existence of sub-millimeter particles is un- ability of these aggregates to support hydrogen deposi- constrained, we can place lower bounds on their typical tion. separation. We assume that a & 100 µm conglomerates permit deposition through a combination of the afore- 2.7. Energy Balance of Non-Ideal Grains mentioned non-ideal effects. While these aggregates do Non-ideal dust conglomerates, in addition to their not reach the minimum temperatures from Figure2, this large size, will be porous, anisotropic, and inhomoge- size is a reasonable upper bound on coagulation (Wei- neous. Because the material strength of individual small denschilling & Ruzmaikina 1994). grains resists compaction in sub-millimeter aggregates, To begin, we express the number of dust grains with coagulated grains are characterized by fractal dimension size ad needed to build a conglomerate with size ac and fractal dimension Df as (Moro-Mart´ın 2019) Df ≈ 2 (Kataoka et al. 2013). While the precise modifi- cations to the energy balance are grain-specific, we can D h ac i f estimate the magnitude of the population-level devia- N ≈ . (6) ad tions required for H2 production. For example, non-spherical geometry modifies the ra- If fraction χ of the original MRN dust participates in coagulation and the resulting conglomerates are tio As/Ag, thus altering Tgr. Greenberg & Shah(1971) consider cylindrical grains, showing that temperatures uniformly-spaced, then the average separation lc is can decrease by 10% versus the spherical limit. If this 1/3 D /3 1/3 h(N/χ)i h ac i f h 1 i effect persists for irregularly-shaped dust, then Tgr will lc ≈ ≈ . (7) nd ad χnd Solid Hydrogen Objects from Starless Cores 7

Table 2. Summary of deposition site spacing lc(ac) √ for noteworthy values of coagulation fraction χ, as- 2 2 −1 λH2 = ( 2π(dH nH2 + dHenHe)) , (8) suming conglomerates form from the largest MRN con- 2

stituents. WR1994 references Weidenschilling & Ruz- where dH2 = 290 pm and dHe = 260 pm are the effec- maikina(1994). tive sizes of H2 and He. In the fiducial core from Ta- 7 ble1, λH2 = 8 × 10 cm. For objects with a λH2 , Model χ lc(100 µm) lc(1000 µm) growth is not diffusion-limited. Hydrogen deposition, Extended MRN 1.0 5.4 × 103 cm 2.5 × 104 cm instead of forming dendritic “branches” akin to terres- WR1994 High χ 10−2 2.5 × 104 cm 1.2 × 105 cm trial snowflakes (Libbrecht 2017), would create spherical WR1994 Low χ 10−8 2.5 × 106 cm 1.2 × 107 cm “hailstones.”. If solid hydrogen freezes on dust grains, any initial anisotropies, porous regions, and other deformities will be covered by ice. Therefore, in order to continue deposition, either the incident radiation bath must be lower First, we consider an “extended MRN model” with than we estimated in Section2 or the albedo of H 2 χ = 1, where nearly all of the dust mass resides in ice must be sufficiently high. Otherwise, the object’s aggregates. Because most of the MRN mass is in temperature would rise above 2.8 K as its geometry ap- 0.25 µm grains, we consider conglomerate growth from proaches spherical symmetry. only these constituents. Using φd from Table1, we find Applying Equation8 to dust-dust encounters gives l = 5400 cm for a = 100 µm deposition sites. 14 c λd & 10 cm for the largest MRN grains. We consider In starless cores, the observation of the 9.7 µm ab- only pairwise collisions with dust of equal or greater sorption feature attributed to MRN-sized silicate grains size, so conglomerates will have even longer λ. Thus, implies that not all of the dust participates in coagula- turbulent gas dynamics dictate the particles’ motion. tion: χ 6= 1 (Draine 2003). Generally, older clouds have Sub-micron dust may align with ambient magnetic fields higher χ and larger conglomerates than younger ones as observed in ISM polarization signatures, but bigger (Ormel et al. 2009). Numerical models from Weiden- grains are unaffected (Draine 2011). schilling & Ruzmaikina(1994) find a range of feasible χ, ranging from 10−2 to 10−8 for the most evolved and 3.2. Effects of Starless Core Turbulence least evolved starless cores. Reynolds Numbers in starless cores are among the Characteristic separations are shown in Table2 for a highest of any astrophysical flow, so turbulence is dif- −1/3 range of deposition site sizes and χ. The lc ∝ χ ficult to resolve numerically without draconian approx- dependence means that separation of aggregates varies imations (Draine 2011). Additionally, effects such as slowly for changing coagulation assumptions. Therefore, ambipolar diffusion mean that the requirements for the better-constrained φd and Df lead to reasonable es- ideal magnetohydrodynamics are unfulfilled (Vázquez- timates of distance between deposition sites. Semadeni et al. 2005). Nonetheless observations show Alternatively, non-MRN particles may be delivered that starless core turbulence empirically fits the sub- from AGB star outflows and other exogenous phenom- sonic Kolmogorov spectrum (Larson 1981), so we adopt ena. Spacing between these potential deposition sites this characterization for our analysis. depends primarily on the starless core’s proximity to the Conveniently, the Kolmogorov framework determines source. While nearness to a luminous source would re- the statistics of turbulent feature strength, size, and sult in temperatures too high for deposition, AGB stars frequency (Weidenschilling & Ruzmaikina 1994). The 5 evolve into White Dwarfs within a few 10 yr and then largest starless core eddies are of characteristic size would not affect ambient temperature. llg ∼ 0.2 Rsc and have flow velocities vlg ∼ βcs with 3. ACCRETION ONTO CLOUD PRODUCTS β . 1. Thus, they have turnover times tlg ∼ llg/(βcs) dictating their survival. Kinetic energy is E ∼ v2 Provided that deposition occurs, we consider the ac- k lg and cascades to smaller features until dissipative viscous cretion of H ice and the remaining free dust onto an iso- 2 forces become efficient (Weidenschilling & Ruzmaikina lated cloud product. We discuss growth to centimeter- 1994). These minimum-size eddies are characterized by sized bodies whose dynamics are independent of the turbulence. 3 3 1/4 3.1. Mean Free Path and Local Vapor Pressure lsm ∼ (η /ρg) , 1/4 The mean free path λ between gas-gas collisions for a vsm ∼ (η/ρg) , 1/2 given H2 molecule is calculated by tsm ∼ (η/ρg) , (9) 8 Levine and Laughlin

Table 3. Turbulent eddy conditions in starless cores. Here, growth-modifying phenomena, an object of radius a ex- leddy, veddy, and teddy are the length scale, flow velocity, and periences an average collision rate of hydrogen molecules turnover time of the turbulent features, respectively, and with thermal velocity vth as across is the minimum size for a hydrogen hailstone to be unaffected by Epstein drag in the eddy. Γ = nH2 Agvth . (13) Eddy Size leddy veddy teddy across With sticking fraction s, we find the volume growth rate Largest 1.5 × 1016 cm 104 cm/s 5 × 104 yr 5 cm dV Smallest 1.8 × 109 cm 30 cm/s 3 yr 3 µm = Γ(m /ρ )s , (14) dt H2 ice

where mH2 is the mass of one molecule of H2. Then, letting dV/dt = As(da/dt), we ﬁnd the radial growth rate (Seligman & Laughlin 2020). With appro- where ρg is the gas density, η = 0.1 µPa·s is the dynamic priate numerical values, we obtain the timescale tdep for 3 viscosity of cold H2 (Souers 1986), and ∼ vlg/llg is the accretion to create spherical hailstones that cross all ed- energy dissipation rate (Weidenschilling & Ruzmaikina dies 1994). Numerical values for ﬁducial starless core eddies 2a ρ are given in Table3. The length scale lsm must be larger lg ice tdep = (15) n sp2kT m than λH2 , so the Kolmogorov spectrum provides a min- H2 sc H2 imum eddy size that is slightly above this theoretical limit. h a ih nH i−1 t ≈ 3.5 × 105 yr lg 2 Now, we consider the dynamics of a particle with mass dep 5 cm 3 × 106 cm−3 mp, radius ap, and density ρp, traveling at velocity vp h s i−1h T i−1/2 sc . (16) through an eddy. Its momentum changes via Epstein 0.5 2.8 K drag, |F | = (4πρ a2 v v /3), on response timescale D g p th p After this time, Epstein drag no longer dominates the (Weidenschilling & Ruzmaikina 1994) hailstone dynamics.

mpvp apρp In Equation 16, we have assumed constant colli- tD ∼ ∼ . (10) FD vthρg sion rate and sticking fraction over the entire time pe- riod. However, deposition releases latent heat Ldep ∼ The particle crosses an eddy with turnover time teddy 0.8 kJ/mol (Roder 1973), so the hailstone must radiate if tD & teddy (Weidenschilling & Ruzmaikina 1994). this energy to maintain its frigid temperature and add Micron-sized dust is too small, and conglomerates are more hydrogen ice. Thus, the heat released by the phase too porous to cross even the smallest eddies. However, transition itself appends another source term to the en- if deposition proceeds and ρp approaches the density ρice ergy balance in Equation1 for cloud products that are of solid H2, then the critical object sizes asm and alg to growing via accretion. cross the smallest and largest eddies are For hailstones where a λcmb, the eﬃciencies of absorption a and emission e approach unity. If hydro- s 2 ρgvthη gen ice objects can decouple from the entire turbulent asm ∼ 2 , (11) ρice spectrum on the lifetime of a starless core, they must dissipate the latent heat such that their growth is lim-

llgρgvth ited by the flux of incident gas particles. Because of alg ∼ . (12) 4 ρicecsβ the T dependence in the Stefan-Boltzmann blackbody emission law, the accretion rate is quite sensitive to the Assuming aggregates a & 100 µm serve as deposition ice temperature. sites, any conglomerates mantled by H2 cross small ed- This thermal constraint reiterates the requirement for dies unimpeded, yet hailstones are somewhat affected the starless core to have excellent shielding from elec- by the turbulence until they reach alg ∼ 5 cm. Objects tromagnetic radiation and possibly pockets of adiabatic larger than alg move independently of all features in the expansion in order to accommodate the latent heat of Kolmogorov spectrum. deposition. As with the dust conglomerate energy balance, the presence of impurity particles may aid the hail- 3.3. Constraints on Deposition Rates stone in efficiently radiating infrared photons. An upper bound on the deposition rate can be found Thus far, we have neglected the dissipation of grav- from the flux of gaseous H2 molecules. Barring other itational potential energy. While unimportant at the Solid Hydrogen Objects from Starless Cores 9 centimeter scale, it provides an upper bound on cloud product size. We equate the latent heat of sublimation for a cloud product to its gravitational binding energy and simplify, finding

3 L R2 = sub . (17) 4π Gρice From Equation 17, we see that gravitational potential energy only restricts the growth of cloud products larger than 4 × 108 cm. This size is much larger than the hailstones that could form via accretion on the characteristic lifetimes of starless cores. 3.4. Gravitational Focusing of Dust Hydrogen deposition dominates for small hailstones, but gravitational focusing of slow-moving dust may meaningfully contribute to mass accretion for large objects. Using conservation of energy and momentum, the impact parameter b of dust with mass md attracted to a body of radius a ad and density ρice is

πGm ρ a4 b2 = a2 + d ice . (18) 3kTsc Here, we have adopted constant vth at Tsc and assumed low cloud product metallicity via ρ = ρice. Figure3 quantifies b and mass flux versus hailstone size. We continue to assume sH2 = αH2 = 0.5, but expect dust sd = 1 from low collision velocities and electrostatic forces. Next, we incorporate gravitational focusing into the radial growth model. We substitute Equation 18 into Equation 13 for particles of mass m, number density n, and sticking fraction s colliding with a hailstone of H2 ice to find

√ √ Figure 3. Panel (a) shows the impact parameter for a given ! 3 ! da 2nsp 2m πGns 2 cloud product radius and mass, assuming ρ = ρice. We plot b = kTscm + √ a . (19) dt 2ρice 6 kTsc for both hydrogen gas and MRN dust. Panel (b) shows hailstone mass ﬂux given the impact parameter from Equation Thus far, we’ve assumed accretion of only one compo- 18 and the ﬁducial starless core conditions in Table1. nent, such as either dust or molecular hydrogen. Setting the constant terms in parentheses as k1 and k2, we get the subscripts correspond to the aforementioned con-

stituents. Because k1,d k1,H2 and k2,d k2,H2 , we da 2 = k1 + k2a . (20) set k1 ≈ k1,H2 and k2 ≈ k2,d. We note that k2, which dt represents only dust flux, does not depend on ρ . Physically, the k term refers to the simple case of accre- ice 1 Figure4 confirms the approximation’s validity by tion onto the geometric cross-section, and the k term 2 comparing Equation 21 to a numerical solution for Equa- adds the impact of gravitational focusing. Solving with tion 18. The latter method considers a single, isolated initial condition a(0) = 0 since the cloud products begin cloud product and computes the mass flux using the small, the hailstone size grows as fourth-order Runge-Kutta method. While the analytic r solution cannot track hailstone metallicity, the numeri- k1 p a(t) = tan( k k t) . (21) cal scheme records this information. For both methods, k 1 2 2 we assume that the accretion rate is limited by the flux For two-component accretion of H and dust, we 2 of colliding gas particles. Equation 21 eventually leads let k1 = k1,H2 + k1,d and k2 = k2,H2 + k2,d, where 10 Levine and Laughlin

rest over which two objects, each with mass M, can gravitationally aggregate in time τ as √ 2GMτ 2/3 r = 2 . (22) 0 π Upon contact, resurfacing merges the bodies. How- ever, constant mass is not necessarily true in starless cores because accretion and hailstone-hailstone gravitation may be coincident. To investigate the assumption, we extend the Runge-Kutta model from Section3 to encompass two cloud products simultaneously accreting mass according to Equation 19 and gravitationally at- tracting. Objects are initialized at alg, so we may ignore Figure 4. Hailstone mass (left-hand axis) versus time for turbulence. analytic and Runge-Kutta solutions for an isolated body as- Figure5(a) shows the coagulation time versus initial suming that growth is limited by collisions with ambient gas separation for eddy-crossing hailstones and the output of particles. Cloud product metallicity is shown on the right- Equation 22 for M = m = 45 g, the mass of hailstones hand axis. Hailstones are initialized at 10−9 g in the numer- lg ical scheme, but the analytic solution begins with zero mass. with radius alg. For fiducial separations of 100 µm grains However, the hailstone masses found from the two methods in Table2, the constant mass assumption is fine, but quickly converge for larger cloud product size. objects with larger r0 accrete significant material and coagulate in less time than predicted by Equation 22. to runaway growth, but objects do not attain the rele- Similarly, panel (b) shows the mass mc of the individual vant sizes on starless core lifetimes. bodies at coagulation versus the initial separation r0. Gravitational focusing leads to negligible dust ac- Previously, we determined that cloud products must cretion at the centimeter-scale, although it is impor- become centimeter-sized before moving independently of tant if cloud products become planetesimal-sized bodies. the turbulence. Now, we’ve shown that the coagulation Nonetheless, linear radial growth for collision-limited ac- time for two such decoupled objects in starless cores is cretion of small hailstones means that growth to eddy- short enough to treat accretion and pairwise aggregation crossing objects is nearly independent of initial conglom- as distinct processes, consistent with Füglistaler& Pfen- erate size. However, it also implies that accretion alone niger(2018). Thus, Equation 22 holds for our fiducial cannot explain a decameter-scale cloud product on star- separations. less core lifetimes. 4.2. Behavior of Multi-Body Systems 4. GRAVITATIONAL ATTRACTION AND Now, we consider a system of N interacting equal mass AGGREGATION OF CLOUD PRODUCTS bodies. Rearranging Equation 22, we find the character- Here, we examine hailstones that are unimpeded by istic attraction timescale τ for one eddy-crossing object drag. Namely, pairwise gravitational forces monopolize and its nearest neighbor as the dynamics and may permit aggregation to kilometer- 3/2 sized “snowballs” (Füglistaler& Pfenniger 2018). We π r0 τ = p , (23) find different characteristic timescales for accretion of 4 Gmlg H2 and gravitational aggregation, permitting us to treat For relevant initial separations, mass aggregation by ac- these as distinct, separate, and non-overlapping pro- cretion is negligible on timescale τ. cesses. In a system of hailstones with average initial sepa- −3 ration r0 = lc, the number density is n0 ∼ lc be- 4.1. Coagulation Time for Two Objects fore coagulation occurs. After one timescale τ, the At small size, pairwise impacts between hailstones are masses of the remaining bodies are approximately dou- determined by an expression with the same form as bled (m1 → 2 mlg), and the number density is roughly Equation 13, leading to negligible growth in this stochas- halved (n1 → 0.5 n0) if the volume in which they tic manner. For larger objects, however, gravitation in- are contained is unchanged. Extending these relation- creases the collision rate. To fix ideas, we adopt the as- ships to find the characteristic timescale τi between the sumption of constant mass hailstones. Conservation of (i−1)th and ith aggregation, we substitute into Equation energy and momentum determine the distance r0 from 23 and simplify to find Solid Hydrogen Objects from Starless Cores 11

Table 4. Estimated snowball radii resulting from the aggregation of all solid H2 hailstones in eddies of size lsm and 10 lsm.

Hailstone lc Eddy Size Snowball Size af 5.0 × 104 cm 3.0 × 109 cm 3.11 km 5.0 × 105 cm 3.0 × 109 cm 0.31 km 5.0 × 106 cm 3.0 × 109 cm 0.03 km 5.0 × 104 cm 3.0 × 1010 cm 31.1 km 5.0 × 105 cm 3.0 × 1010 cm 3.11 km 5.0 × 106 cm 3.0 × 1010 cm 0.31 km

Next, we solve for the ﬁnal radius af of objects formed after Nf mass doubling events as

Nf 1/3 3mlg2 af = . (25) 4πρice

These aggregated H2 objects, “snowballs,” should have final sizes set by the turbulent spectrum in which the original hailstones are contained. Because small eddies are the most plentiful and the weakest, we expect that features of a few lsm will determine the most common snowball sizes. Ambient gas in these eddies flows with approximate velocity 1 m/s, so the frontal collisional heating from Epstein drag described in Hoang & Loeb(2020) is unimportant. Table4 shows the resulting object radii from coagulat- ing hailstones with size alg in eddies of size leddy, assuming no additional accretion. We find that kilometer-scale cloud products could aggregate within a few 104 yr. If the hailstones are attracted to a common center of mass, Figure 5. Panel (a) shows the coagulation time τ versus then the time to coagulate all of the objects in an eddy initial separation r , assuming the pairwise gravitation of 0 would be shorter than this estimate. After aggregation, cloud products with mass mlg ∼ 45 g. Points corresponding 11 to fiducial deposition site spacing values χ are marked for the planetesimal-like objects will collect a few 10 g of reference. Panel (b) shows the mass of objects at coagulation dust via gravitational focusing on GMC lifetimes. While for the same initial separation. the newly-aggregated cloud products would still be neg- ligibly metallic, dust will impact the sublimation rate for hydrogen located near these impurities. π Finally, we must consider the attraction of bodies to τ = (n Gm )−1/2 = τ , (24) i 4 0 lg the center of mass in the starless core. For an object of mass M located distance d from the core’s cen- which is independent of i. Thus, once pairwise gravi- ter with mass M interior to d, the Hill Sphere is tation dominates, we expect roughly constant time be- in r ≈ d(M/3M )1/3. We expect snowballs to form in tween each mass-doubling event. Equation 24 takes the H in √ dark regions with small M and thus, r l . Fur- form of Jeans’ collapse: τ ∝ 1/ Gρ (Jeans 1902). Es- in H c J thermore, the ratio of r to the average snowball sepa- sentially, this result can be understood as the free-fall H ration is unchanged for successive mass-doubling events, of a region of hailstones that are unaffected by turbu- 4 so coagulation is always unaffected by cloud gravity. lence. For lc = r0 = 2.5 × 10 cm in Table2, we find characteristic τ ≈ 60 yr. We also note that this result is 4.3. Revisiting the Turbulent Spectrum independent of the ambient starless core temperature. 12 Levine and Laughlin

So far, we have characterized turbulence via the (Füglistaler& Pfenniger 2018) which could scatter snow- Kolmogorov framework, despite knowing that the gas balls, the drag heating induced on the tenuous H2 ice is compressible and partially ionized. Instead, one makes this event less favored than more gentle disper- could reasonably choose a sub-Alfvénic regime or sal. Moreover, a scattering event may impart a relative numerically-simulate the non-ideal magnetohydrody- velocity to the snowball versus the environment. As- namics. Nonetheless, only the following general prop- suming the GMC travels with the LSR, an ejected body erties of subsonic 3-D flow are important to building could have different kinematics compared to ‘Oumua- hypothetical cloud products: mua (Mamajek 2017; Hallatt & Wiegert 2020). Due to the barriers imposed by ejections, snowballs 1. Kinetic energy should cascade from eddies with must be released by the dissipation of the ambient gas. large sizes, high velocities, and long turnover times Re-expansion of cores can occur, especially in unproduc- to weaker, but more plentiful features. tive star-forming regions like the Columba and Carina associations (Vázquez-Semadeniet al. 2005). Denser 2. The smallest eddies should be a few λH in size. 2 cores are more apt to collapse, but because their struc- 3. Turbulence should be sufficiently mild so that ture is governed by the turbulence from which they orig- it neither shatters dust conglomerates nor subli- inate, a diversity of behavior is expected (Draine 2011). mates H2 ice via heating from Epstein drag. 5.2. Thermal Sublimation in the Interstellar Medium Since cores are the most quiescent GMC regions and have among the lowest ion fractions in the ISM (Draine Once the starless core dissipates via radiation pressure 2011), these attributes should apply for any appropriate (Blitz & Shu 1980), net sublimation begins. Specifically, choice of turbulence. While the Kolmogorov cascade ap- Hoang & Loeb(2020) quantify erosion for a pure H 2 proximately fits observations, an alternate characteriza- snowball with temperature Tice via the crystal lattice 12 −1 tion would still find similar plausibility of cloud product oscillation frequency (ν0 ∼ 10 s ) as formation. 1/3 da mH2 −Eb = −ν0 exp , (27) 5. PROCESSING IN INTERSTELLAR SPACE dt ρice kTice Now, we examine the feasibility of liberating snowballs where Eb is the binding energy and related to Lsub. from the starless core and discuss processing in interstel- Evaporative cooling would consistently maintain ice temperature around 3 K (Hoang & Loeb 2020), so Equa- lar space. We find a lower bound on the size of H2 cloud products that could survive to the Solar System. tion 27 gives the minimum radius of a hydrogen snowball that would survive an ISM journey time of ttravel as 5.1. Release of Snowballs into Interstellar Space ttravel Once formed, snowballs must leave the starless core a ∼ 3 km , (28) min 30 Myr to become true interstellar objects. Cloud products in a collapsing, star-forming structure will be destroyed, in approximate agreement with Hoang & Loeb(2020). so only failed cores can produce long-lasting H2 ice. This characteristic timescale independently corrobo- Furthermore, we can eliminate the possibility of a vio- rates the observation of an extreme aspect ratio for lent ejection because of the energies involved. Assuming ‘Oumuamua despite its young kinematic age, implying snowballs form at the core center, we can find the energy that objects of its class are transient with respect to the needed for ejection per unit mass Eej/M as Galactic lifetime. Additionally, Hoang & Loeb(2020) discuss enhanced collisional heating on the frontal snow- Z ∞ Eej Min(r)dr ball cross-section due to motion relative to the ambi- = G 2 , (26) M 0 r ent gas. However, this assumes a high velocity versus where Min(r) is the core mass interior to radius the LSR, inconsistent with observations of ‘Oumuamua’s r. Na¨ıvely assuming a homogeneous sphere, we find kinematics. Eej/M ∼ 9.8 kJ/g. For realistic density distributions, energies remain on the same order of magnitude. Be- 5.3. Hydrogen Loss for High Metallicity Remnants cause (Eej/M) > Lsub, any ejection event would need to Hydrogen loss increases snowball metallicity over time provide multiple times the sublimation energy without because the non-volatile dust grains remain attached to destroying the body. the ISO. If ‘Oumuamua hosted H2 ice, the surface cover- While gravitational instabilities in starless cores age fraction must have only been around 6% in the inner can assume a variety of behaviors and morphologies Solar System (Seligman & Laughlin 2020). In contrast, Solid Hydrogen Objects from Starless Cores 13 a pure H2 object’s non-gravitational acceleration would individually, which is insufficient to explain ‘Oumua- have exceeded observations by a factor of 16. mua’s non-gravitational acceleration. While the small- Because the binding energy of H2 is higher to im- est MRN grains could store sufficient hydrogen, there purities than to itself, we expect negligible sublima- is no known mechanism to create decameter-sized ISOs tion from silicate surfaces thinly-covered in H2 (Walker from mantled dust grains. Thus, if ‘Oumuamua were a 2013). Thus, Equation 27 is invalid for non-negligible cloud product, it must have originated with more than metallicity, even at a temperature above 3 K. Most dust a thin H2 coating. collection occurs when the cloud product is close to its maximum size, so its outer layers will be more metallic. 6. DISCUSSION Thus, we could expect shielding of inner H2 and slower 6.1. Future Detection of Molecular Hydrogen mass loss, but it’s unlikely to be an order-of-magnitude There has yet to be an indisputable discovery of in- effect. terstellar solid H2, and most searches have targeted em- Regardless of cloud product metallicity, however, cos- bedded hydrogen within common ISM ices (Sandford mic rays destroy solid H2 and chemically process non- et al. 1993). While Lin et al.(2011) notes that detec- volatile constituents. Thus, for objects dominated by + tion of the H6 radical would confirm the presence of refractory substances, bombardment by energetic par- UV-bombarded H2, the precise spectroscopic signature ticles sets the timescale for hydrogen ice loss. Selig- has not been experimentally-determined. Nonetheless, man & Laughlin(2020) consider an ellipsoidal object, the future confirmation of processed H2 ice would not finding that cosmic ray flux with strength similar to necessarily imply the existence of macroscopic ISOs. the value within the Solar neighborhood can destroy To confirm that ‘Oumuamua-like objects have an H2 7 1 an ‘Oumuamua-sized snowball within a few 10 yr . component, in-situ missions would provide the most Thus, cloud products with aggregated silicates will erode definitive answers (Snodgrass & Jones 2019). Addition- rapidly at first from sublimation, then more slowly as ally, the forthcoming Vera Rubin Observatory should cosmic rays become the primary loss process. identify ISOs before their closest approach to Earth and Unless they acquire substantial impurities, kilometer- spur follow-up studies of both their intrinsic properties scale snowballs would vanish within 10 Myr. Thus, and their dynamics (Cook et al. 2016). Narrower density any H2 cloud products passing through the Solar Sys- constraints and analysis of a longer trajectory than was tem must have been originally larger and formed early available for ‘Oumuamua could differentiate between ra- enough in GMC lifetime to attract sufficient dust. Con- diation pressure, solid hydrogen, and other hypotheses. sequently, observed interstellar interlopers would be bi- ased towards initially planetesimal-sized bodies. 6.2. Solid Hydrogen and Cosmological Evolution

5.4. Preservation of a Thin Hydrogen Coating Because CMB temperature scales with cosmological redshift z as Tcmb = (1 + z) 2.725 K, solid hydrogen for- Thin hydrogen coatings could survive almost indef- mation is hindered at z & 0.1 (Hoang & Loeb 2020). initely on cloud products (Walker 2013). Nonetheless, Figure6 plots the number density “ice line” versus red- 11 we must reconcile the preserved H2 mass with the 10 g shift, showing that solid H2 could form only within the of outgassing material required to explain ‘Oumuamua’s past 750 Myr2. If solid hydrogen ISOs occur in the non-gravitational acceleration (Seligman & Laughlin present epoch, we would expect increasing production ˚ 2020). For example, 10A binding sites (Hoang & Loeb in the Universe’s far future (White 1996). 2020) covering a smooth 115 m×115 m×20 m ellipsoid in a monolayer of molecular hydrogen could only preserve 6.3. GMC Constraints from Solid Hydrogen ISOs a few grams of H2. Thus, a cloud product ‘Oumuamua If ‘Oumuamua-like objects contain dynamically- could not have been this weathered. important H2, it would illuminate previously- Given the strict thermodynamic constraints on form- unobserved conditions in the coldest starless cores. ing hydrogen ice cloud products, we evaluate the forma- Furthermore, we could constrain the eﬃciency of the tion of snowballs directly from H2-coated dust. To as- formation process. Do et al.(2018) estimates the semble an ISO with mass 1012 g directly from MRN par- −3 ISO number density as niso ∼ 0.2 AU with 3.5 24 ticles, approximately 10 grains with radius 0.25 µm are years of Pan-STARRS data, so we update this to 7 required. Only 10 g of H2 could be stored on the grains

2 Lookback times were found via the Cosmological Calculator at 1 Simulations of ‘Oumuamua’s energy balance can be accessed at https://home.fnal.gov/∼gnedin/cc/ for the default ﬂat Universe https://github.com/DSeligman/Oumuamua Hydrogen. (Hubble’s constant 68.14 km/s/Mpc, matter density 0.3036). 14 Levine and Laughlin

Table 5. Summary of relevant volumes and timescales to form 15 an ISO with 10 g of H2 mass via an unknown catalytic process.

Dust Size Formation Time Box Side l Box Side l/lsm 0.005 µm 3 × 105 yr 5.0 × 109 cm 2.8 0.25 µm 3 × 105 yr 9.5 × 109 cm 5.3 0.005 µm 3 × 106 yr 2.3 × 109 cm 1.3 0.25 µm 3 × 106 yr 4.5 × 109 cm 2.5

Assuming the MRN distribution of dust sizes (dn/da ∝ a−3.5), nearly all of the surface area will be contained in the smallest grains. Therefore, we consider only these grains as catalysis sites in our order-of-magnitude com- putation. Thus, we can determine the GMC volume Figure 6. Saturation number density of H at T , plotted 2 cmb with characteristic side length l required to produce the for given redshift and approximate time before present. amount of H2 required for ‘Oumuamua’s survival in the ISM. −3 niso ∼ 0.1 AU to account for additional time without 2 another ‘Oumuamua-like object. We consider the mass hnH πa vH mH stnd i1/3 l = 2 d 2 2 (32) of snowballs miso at their maximum, before processing M in the ISM. Due to ‘Oumuamua’s low velocity relative to Table5 shows a few example box sizes required for the the LSR, we assume that the Earth is passing through production of 1015 g ISOs. Coincidentally, these results a group of cloud products that have barely moved from are close to the H2 mean free path that creates natu- that fiducial Galactic orbit since their formation. There- ral length scales. Despite this correspondence, there is fore, the efficiency e is no known chemical process that would both freeze and n m release H2. Furthermore, free hydrogen crystals would e ≈ iso iso , (29) ρcl experience the ambient photon bath discussed Section 2 and rapidly sublimate. Thus, catalytic processes, if where ρcl is the GMC mass density. Thus, we find they occur at all, probably cannot produce the multi- kilometer objects required for survival to the Solar Sys- h ih i3h 3 −3 i −5 niso riso 10 cm tem. e ∼ 2.5 × 10 −3 .(30) 0.1 AU 2 km nH2 6.5. Implications for the Population of ISOs Thus, ISO formation in a snowball-spawning cloud would be approximately 400 times less efficient than star We wish not to speculate on the broader population formation (Blitz & Shu 1980). If ‘Oumuamua-like ob- of objects from just ‘Oumuamua, but rather hypothe- jects have solid H2, population-level statistics will con- size generic attributes based on the characteristics of siderably refine our estimate. our outlined formation mechanism. Because the approximate lifetime of the solid H2 component is 50 Myr, cloud 6.4. Catalytic Processes as an Alternative Mechanism products of older age would be left solely with refractory In Section3, we considered formation where freezing substances. The remnant population’s number density H2 attaches mostly onto other H2 molecules. However, depends primarily on the exact lookback time at which another hypothetical ISO assembly mechanism may in- the formation of their solid hydrogen predecessors be- volve a catalytic process where deposition only occurs came feasible. on dust grains, perhaps at slightly higher temperature Additionally, ISOs with active H2 sublimation would than Tgr = 3 K. If solid hydrogen particles are formed be concentrated near GMCs with failed cores. Thus, and released, then the cross-section Ag in Equation 13 the spatial distribution of young H2 snowballs would is constant. In time t, one grain of radius a can catalyze be anisotropic, whereas the eroded remnants would deposition for mass form a homogeneous Galactic background. Both new and old cloud products should pass through the So- 2 Mcat = nH2 πa vthmH2 st . (31) lar System, but only the former should accelerate non- Solid Hydrogen Objects from Starless Cores 15 gravitationally. Within the next decade, a larger sam- mum radius for survival given a travel time. The H2 ple of objects will further constrain the population-level fraction at the Solar encounter is determined by the statistics of ISOs. interstellar journey time and the mass of silicates attracted via gravitational focusing. Young objects typ- 7. SUMMARY AND CONCLUSIONS ically have inbound velocities near the Local Standard

Beginning from the observation that H2 outgassing of Rest, as they have not experienced strong scattering would explain many of ‘Oumuamua’s exceptional at- events (Hallatt & Wiegert 2020). tributes (Seligman & Laughlin 2020), we analyze the fea- Currently, solid hydrogen is the only hypothesis that sibility of a generic ISO formation mechanism in GMCs. explains all of ‘Oumuamua’s observed properties – non- Because of the severe thermodynamic criteria for de- gravitational acceleration, lack of cometary volatiles, position, only failed cores in unproductive star-forming high aspect ratio, inbound kinematics, and number den- regions could possibly spawn these bodies. sity – as generic to a class of naturally-forming objects. Table6 summarizes the size regimes and physical pro- Thus, despite its exotic nature, the H2 cloud product cesses we have considered in the hypothetical production remains as compelling as any other interpretation. The of generic ‘Oumuamua-like ISOs. Specifically, we exam- composition of ‘Oumuamua-like ISOs will be testable ine starless cores for the possibility that they may spawn with forthcoming observational surveys, and the answer solid H2 cloud products. We find that the strongest bar- may illuminate a previously-unknown avenue of small rier is the temperature requirement for H2 deposition. body formation. Because molecular clouds are not opaque in the infrared, both the CMB and radiation from warm dust can pene- We thank Daniel Pfenniger for a detailed peer-review, trate to all relevant visual extinctions. Nonetheless, it is which greatly improved the manuscript. Additionally, possible that deposition could occur on non-ideal grain we thank Bruce Draine and Dave Hollenbach for the conglomerates. engaging conversations on ISM physics and astrophys- If deposition occurs, hailstones that cross all turbulent ical dust. We also thank Joel Ong for the comment eddies in the Kolmogorov spectrum could form within on the figure caption corresponding to the accretion 350, 000 yr given efficient dissipation of the latent heat. model. Finally, we are grateful to Darryl Seligman, At this macroscopic size, pairwise attraction between Mark Walker, Fred Adams, and Christopher Lindsay cloud products becomes the catalyst for mass aggrega- for providing comments from closely-reading our draft tion. In the ISM, thermal sublimation determines the paper. This research used the numpy (Harris et al. 2020) hydrogen component’s lifetime and introduces a mini- and matplotlib (Hunter 2007) packages in python.

REFERENCES

Air Liquide. 1976, Gas Encyclopedia (Amsterdam: Elsevier) Füglistaler,A., & Pfenniger, D. 2016, A&A, 591, A100 Anderson, H. L. 1989, A Physicist’s Desk Reference, 2nd —. 2018, A&A, 613, A64 edn. (American Institute of Physics) Gaidos, E. 2018, MNRAS, 477, 5692 Bergström,L. 2000, RPPh, 63, 793 Goldsmith, P. F. 2001, ApJ, 557, 736 Bialy, S., & Loeb, A. 2018, ApJL, 868, L1 Gounelle, M., & Meibom, A. 2008, ApJ, 680, 781 Blitz, L., & Shu, F. H. 1980, ApJ, 238, 148 Greenberg, J., & Shah, G. 1971, A&A, 12, 250 Boogert, A. A., Gerakines, P. A., & Whittet, D. C. 2015, Hallatt, T., & Wiegert, P. 2020, AJ, 159, 147 ARA&A, 53, 541 Harris, C. R., Millman, K. J., van der Walt, S. J., et al. Burke, J., & Hollenbach, D. 1983, ApJ, 265, 223 2020, Nautr, 585, 357 Cook, N. V., Ragozzine, D., Granvik, M., & Stephens, Hoang, T., & Loeb, A. 2020, ApJL, 899, L23 D. C. 2016, ApJ, 825, 51 Hollenbach, D., & Salpeter, E. 1971, ApJ, 163, 155 Do, A., Tucker, M. A., & Tonry, J. 2018, ApJL, 855, L10 Hunter, J. D. 2007, Computing in science & engineering, 9, Domokos, G., Sipos, A. A.,´ Szabó,G. M., & Várkonyi, P. L. 90 2017, RNAAS, 1, 50 Ivlev, A., Padovani, M., Galli, D., & Caselli, P. 2015, ApJ, Draine, B. T. 2003, ARA&A, 41, 241 812, 135 Draine, B. T. 2011, Physics of the Interstellar and Jeans, J. H. 1902, Philosophical Transactions of the Royal Intergalactic Medium (Princeton, NJ: Princeton Society of London. Series A, Containing Papers of a University Press) Mathematical or Physical Character, 199, 1 16 Levine and Laughlin

Table 6. Summary of size regimes considered in the formation of solid H2 cloud products. For each characteristic object, we state the physical processes, relevant formation timescales, and potential barriers to formation.

Size Regime Object Type Relevant Processes Timescale Possible Barriers

0.005 µm < a < 0.25 µm MRN-Sized ISM GMC complex and starless 3 Myr None. Dust Grains core formation.

0.25 µm < a < 100 µm Large Dust Grain coagulation and 1 Myr Starless cores must support Conglomerates possible exogenous delivery. formation of large grains.

100 µm < a < 5 cm Solid Hydrogen Deposition limited by 350 kyr Tsc must be below 2.8 K. Hailstones hydrogen gas ﬂux. Latent heat must dissipate.

5 cm < a < 5 km Aggregating Gravitational coagulation of 10 kyr Initial separation must be Snowballs eddy-crossing hailstones. small for gravitation.

100 m < a < 5 km Interstellar Thermal sublimation and 50 Myr ISO must aggregate enough Snowballs cosmic ray bombardment. dust to preserve H2.

a ∼ 100 m ‘Oumuamua-Like Non-gravitational 5 yr None. Interlopers acceleration, outgassing.

Jewitt, D., & Luu, J. 2019, ApJL, 886, L29 Meech, K. J., Weryk, R., Micheli, M., et al. 2017, Natur, Jewitt, D., Luu, J., Rajagopal, J., et al. 2017, ApJL, 850, 552, 378 L36 Micheli, M., Farnocchia, D., Meech, K. J., et al. 2018, Kataoka, A., Tanaka, H., Okuzumi, S., & Wada, K. 2013, Natur, 559, 223 A&A, 557, L4 Militzer, B. 2005, Journal of low temperature physics, 139, Kim, S.-H., & Martin, P. 1996, ApJ, 462, 296 739 Kong, S., Tan, J. C., Caselli, P., et al. 2016, ApJ, 821, 94 Moro-Mart´ın,A. 2019, ApJL, 872, L32 Landgraf, M., Baggaley, W., Grün,E., Krüger,H., & Ormel, C., Min, M., Tielens, A., Dominik, C., & Paszun, D. Linkert, G. 2000, JGRA, 105, 10343 2011, A&A, 532, A43 Larson, R. B. 1981, MNRAS, 194, 809 Ormel, C., Paszun, D., Dominik, C., & Tielens, A. 2009, Libbrecht, K. G. 2017, AnRMS, 47, 271 A&A, 502, 845 Lin, C. Y., Gilbert, A. T., & Walker, M. A. 2011, ApJ-, Pfenniger, D., & Combes, F. 1994, A&A, 285, 94 736, 91 Roder, H. M. 1973, Survey of the Properties of the Loinard, L., & Allen, R. J. 1998, ApJ, 499, 227 Hydrogen Isotopes Below Their Critical Temperatures, Luu, J. X., Flekkøy, E. G., & Toussaint, R. 2020, ApJL, Vol. 641 (US Government Printing Office) 900, L22 Sahai, R., Vlemmings, W., & Nyman, L.-A.˚ 2017, ApJ, 841, Mamajek, E. 2017, RNAAS, 1, 21 110 Mashchenko, S. 2019, MNRAS, 489, 3003 Sandford, S. A., Allamandola, L. J., & Geballe, T. R. 1993, Mathis, J., Mezger, P., & Panagia, N. 1983, A&A, 128, 212 Sci, 262, 400 Mathis, J. S., Rumpl, W., & Nordsieck, K. H. 1977, ApJ, Schneider, A. C., Shkolnik, E. L., Allers, K. N., et al. 2019, 217, 425 AJ, 157, 234 Solid Hydrogen Objects from Starless Cores 17

Seligman, D., & Laughlin, G. 2020, ApJL, 896, L8 V´azquez-Semadeni,E., Kim, J., Shadmehri, M., & Shu, F. H., Adams, F. C., & Lizano, S. 1987, ARA&A, 25, Ballesteros-Paredes, J. 2005, ApJ, 618, 344 23 Walker, M. A. 2013, MNRAS, 434, 2814 Walker, M. A., & Wardle, M. J. 2019, ApJ, 881, 69 Snodgrass, C., & Jones, G. H. 2019, NatCo, 10, 1 Weidenschilling, S., & Ruzmaikina, T. 1994, ApJ, 430, 713 Souers, P. C. 1986, Hydrogen Properties for Fusion Energy Werner, M. W., & Salpeter, E. 1969, Monthly Notices of (Berkeley, CA: University of California Press) the Royal Astronomical Society, 145, 249 Spitzer Jr, L. 1941, ApJ, 93, 369 Westphal, A. J., Stroud, R. M., Bechtel, H. A., et al. 2014, Trilling, D. E., Mommert, M., Hora, J. L., et al. 2018, AJ, Sci, 345, 786 White, R. S. 1996, Ap&SS, 240, 75 156, 261