abstract:  A mechanism is identified whereby dark matter (DM) in protostellar halos dramatically alters the current theoretical framework for the formation of the first stars. Heat from neutralino DM annihilation is shown to overwhelm any cooling mechanism, consequently impeding the star formation process and possibly leading to a new stellar phase. A "dark star" may result: a giant ($\mathrel{\rlap{\lower 4pt\hbox{\hskip 1pt$\sim$}}}>1$ AU) hydrogen-helium star powered by DM annihilation instead of nuclear fusion. Observational consequences are discussed.
author: and
title: Dark matter and the first stars: a new phase of stellar evolution

The first stars in the universe mark the end of the cosmic dark ages, reionize the universe, and provide the enriched gas required for later stellar generations. They may also be important as precursors to black holes that coalesce and power bright early quasars. The first stars are thought to form inside haloes of dark matter of mass $10^5
- 10^6 M_\odot$ at redshifts $z=10-50$. These haloes arose from the merging of smaller structures, as overdense regions in the universe assemble heirarchically into ever larger haloes. The haloes consist of 85% dark matter and 15% baryons in the form of pristine hydrogen and helium (from Big Bang nucleosynthesis). The baryonic matter cools and collapses via molecular hydrogen cooling into a single small protostar at the center of the halo (for reviews see e.g. ).

In this paper we consider the effect of the dark matter particles on the formation process of these first stars. We focus on the most compelling dark matter candidate, the lightest supersymmetric particle. Supersymmetry (SUSY), at this point a beautiful theoretical construct, has the capability of addressing many unanswered questions in particle theory as well as providing the underpinnings of a more fundamental theory such as string theory. If SUSY is right, then for every known particle in the universe, there is an as yet undiscovered partner. The lightest of these, the Lightest Supersymmetric Particle or LSP, would provide the dark matter in the universe. The search for SUSY is one of the motivations for building the Large Hadron Collider at CERN, and one may hope that it will be discovered as early as 2008.

The LSP is the favorite dark matter candidate of many physicists. This is true not only because of the beautiful properties of SUSY, but also because the LSP automatically has the right properties to provide 24% of the energy density of the universe. In particular, the neutralino, the SUSY partner of the W, Z, and Higgs bosons, has the required weak interaction cross section and $\sim$ GeV - TeV mass to give the correct amount of dark matter in the universe today. The SUSY particles are in thermal equilibrium in the early universe, and annihilate among themselves to produce the relic density today. It is this same annihilation process that is the basis of the work we consider here. The SUSY particles, also known as WIMPs (Weakly Interacting Massive Particles), annihilate with one another wherever their density is high enough. Such high densities are achieved in the early universe, in galactic haloes today , in the Sun and Earth , and, as we will show, also in the first stars. As our canonical values, we will use the standard value $\langle \sigma v \rangle = 3 \times 10^{-26} {\rm
cm^3/sec}$ for the annihilation cross section and $m_\chi = 100$ GeV for the SUSY particle mass; but will also consider a broader range of WIMP masses (1 GeV–10 TeV) and cross-sections [^1]. The effects we find apply equally well to other WIMP candidates, such as sterile neutrinos or Kaluza-Klein particles.

*DM density profile.* In order to study the effects of DM on the star formation process, we need to know its density profile inside the baryonic core. While simulations have obtained remarkably good density profiles for the collapsing protostellar gas, they have unfortunately (as yet) been unable to do so for the DM. Thus we use adiabatic contraction to obtain estimates of the DM profile. As our initial conditions, we take an overdense region of $10^5-10^6 M_\odot$ with an NFW profile for both DM and gas, where the gas contribution is 15% of that of the DM. (For comparison, we also use a Burkert profile , which has a DM core before contraction.) Then, as the gas collapses, we allow the DM to respond to the changing baryonic gravitational potential (gas density profiles taken from simulations of ). The final DM density profiles are computed with adiabatic contraction ($M(r)r$ = constant) and are shown in Fig. 1 for concentration parameter $c=10$ at a redshift $z=19$ and halo mass $M=10^6 M_\odot$. It is important to point out that our results do not change much for other choices of these parameters; e.g. even for $c=1$, the dark matter density only changes by a factor of 4. After contraction, the DM density at the outer edge of the baryonic core is roughly $\rho_\chi \simeq 5 {\rm GeV/cm}^{-3}
(n/{\rm cm}^{3})^{0.81}$ and scales as $\rho_\chi \propto r^{-1.9}$ outside the core.

Our adiabatically contracted NFW profiles match the DM profile obtained numerically in  (see their Fig.2). They present two DM profiles (their earliest and latest profiles), for $n \sim 10^3 {\rm cm}^{-3}$ and $n \sim 10^{13} {\rm cm}^{-3}$ as far inward as $5 \times 10^{-3}$pc and $0.1$pc. The slope of these two curves is the same as ours. $\rho_\chi \propto r^{-1.9}$. If one uses this slope and extrapolates inward to smaller radii and to higher densities, then one obtains the same DM densities as with our adiabatic contraction approach. We are encouraged by this agreement. It is interesting to note that adiabatically contracted objects of any mass, even outside of the context of Pop III stars, might lead to enhanced DM annihilation leading to photons observable by GLAST or to a modified reionization history.

As a caveat, we note that the approach of adiabatic contraction must be used with caution. It formally requires the orbital particle time to be short compared to the collapse time, though in practice the method works well beyond this limit (see e.g. ). We are also concerned about the requirement of spherical symmetry, when in fact there are filaments, clumps, and mergers, so that dynamical friction or violent relaxation may take place. The gas may form bars that could sweep out DM, though the rotation timescale for bars may be too long for them to be important. We encourage simulators to improve the DM resolution in first stars to confirm our results. The closest previous work is that of Merritt , who used initial profiles $\rho \propto r^{-\gamma}$ with $\gamma = [0,2]$ around a central black hole and found final profiles $\rho \propto
r^{-\gamma'}$ with $\gamma'=[2.25,2.5]$ (i.e. even steeper than we what we found). We caution the reader that our heating estimates below rely upon DM densities obtained with these assumptions.

*DM Heating.* WIMP annihilation produces energy at a rate per unit volume $Q_{\rm
ann} = \langle \sigma v \rangle \rho_\chi^2/m_\chi \linebreak \simeq  1.2
\times  10^{-29} {\rm erg/cm^3/s} \,\,\, (\langle \sigma v \rangle / (3
\times 10^{-26} {\rm cm^3/s}))  \linebreak(n/{\rm cm^{-3}})^{1.62} (m_\chi/(100
{\rm GeV}))^{-1}$. In the early stages of Pop III star formation, when the gas density is low, most of this energy is radiated away. However, as the gas collapses and its density increases, a substantial fraction $f_Q$ of the annihilation energy is deposited into the gas, heating it up at a rate $f_Q Q_{\rm ann}$ per unit volume.

Previous work considered the effects of DM annihilation on earlier low density phases of Pop III star formation ($n \mathrel{\rlap{\lower 4pt\hbox{\hskip 1pt$\sim$}}}<10^4 {\rm cm}^{-3}$). They rightly concluded that 100 GeV neutralinos cannot heat the protostar at these low densities because their annihilation products simply escape out of the object without depositing much energy inside. They consequently focused on lighter particles, such as 1–10 keV sterile neutrinos and 1–100 MeV light dark matter.

However, we find that, for WIMP mass $m_\chi = 100$GeV (1 GeV), a crucial transition takes place when the gas density reaches $n>
10^{13} {\rm cm}^{-3}$ ($n>10^9 {\rm cm}^{-3}$). Above this density, most of the annihilation energy remains inside the core and heats it up to the point where further collapse of the core becomes difficult. At this point the DM density is 2% (10%) of the gas density in the core, the size of the core is $\sim 17$ A.U. ($\sim 960$ A.U.), its mass is $\sim 0.6 M_\odot$ ($\sim 11 M_\odot$), $f_Q \sim 2/3$, and the energy produced by DM heating is $\sim 140 L_\odot (\sim 835
L_\odot)$.

Here we estimate the fraction $f_Q$ of DM annihilation energy that remains inside the gas core (more detailed work is in progress). $f_Q$ scales linearly with the gas density and depends on the relative number of the various annihilation products and their energy spectrum. The latter are heavily dependent on the WIMP model. From our experience with neutralino DM, we assume the following typical values: about 30% of the energy goes into neutrinos, 30% into photons, and 30% into stable charged particles like electrons and positrons. Unstable particles like neutral pions, charged pions, and muons decay into photons, neutrinos, and electrons before exiting the cloud. The energy spectrum of photons and electrons depends to some extent on the exact annihilation channels. For our purpose, we consider typical spectra produced in Pythia simulations of 500 GeV neutralino annihilation . Other spectral shapes will change the precise values of our results but not the overall effect.

Neutrinos escape from the cloud without depositing an appreciable amount of energy. Electrons below a critical energy $E_c \approx 280$ MeV in hydrogen lose energy predominantly by ionization. Higher energy electrons do so by emission of bremsstrahlung photons. As these bremsstrahlung photons pass near gas nuclei, they create electrons and positrons, which in turn may generate other bremsstrahlung photons. Thus starts a sequence of photon, electron, and positron conversions: an electromagnetic (EM) cascade. While photons above $\approx 100$ MeV also initiate an EM cascade, lower energy photons transfer most of their energy to electrons in the gas (Compton scattering).

We approximate the energy loss of electrons via ionization with $4.41{\rm MeV}/E$ (g/cm$^2$)$^{-1}$. For EM cascades, we assume a gamma distribution for the mean longitudinal profile. Thus the fraction of energy lost in traversing a thickness $X$ of gas equals $\gamma(a,bX/X_0)/\Gamma(a)$, where $\gamma(x,y)$ is the incomplete gamma function, $X_0=63$ g/cm$^2$ is the radiation length in hydrogen, $a=1+b[\ln(E/E_c)-0.5]$, and $b=0.5$ . We estimate the core thickness as $X=1.2 m_p n r_0$. Here $m_p$ is the proton mass, $r_0$ is the core radius, and the factor of 1.2 is appropriate for a uniform sphere. We model the fraction of energy loss of photons by converting each photon to an electron of the same energy after one photon attenuation length. The latter is computed from formulas in , interpolated to produce the hydrogen curve in , figure 27.16.

Annihilation of DM outside the core can also contribute to heating inside the core. The annihilation products are not stopped until they reach the core because the baryon density outside falls as $r^{-2.3}$. We integrate the effects of DM particles annihilating in the external region at $r>r_0$, whose annihilation products have a probability of hitting the core given by $\pi r_0^2 / 4 \pi r^2$. Taking $\rho_\chi \sim 1/r^2$ outside the core, we find a 25% enhancement in core heating due to the external region. Henceforth our results are obtained without this enhancement factor.

*Results.* To compare with DM heating, we include all relevant cooling mechanisms. The dominant mechanism is H$_2$ cooling (we use the rates in ), though our code also includes other effects such as H line cooling (important only at high temperatures) and Compton cooling (important only if the gas becomes very ionized) . We use the opacities from ; e.g. at $n \sim 10^{13} {\rm
cm}^{-3}$, we take a $\sim 8\%$ cooling efficiency. Setting the heating rate equal to the cooling rate gives the critical temperature $T_c(n)$ at a given density $n$ below which heating dominates.

In figure 2 we compare $T_c(n)$ to typical evolution tracks in the temperature-density phase plane. The blue (solid) and green (dotted) lines show the temperature evolution of the protostellar gas in the simulations (without DM) of and respectively. The red (dashed and dot-dashed) lines show the critical temperature: below these lines, DM heating dominates over all cooling mechanisms.

Figure 2 illustrates results for a range of WIMP masses from 1 GeV–10 TeV for a canonical $3 \times 10^{-26} {\rm cm^3}/{\rm sec}$ annihilation cross-section. Since the heating rate scales as $\langle \sigma v \rangle / m_\chi$, these same curves equivalently apply to a variety of cross-sections for a given WIMP mass.

The important result is that the blue/green (evolutionary) and red (critical temperature) lines always cross, regardless of WIMP mass or H$_2$ fraction: this is a robust result. As soon as the core density reaches this crossing point, the DM heating dominates inside the core and changes its evolution. Notice that at $m_\chi=$`<!-- -->`{=html}1 GeV, the crossing point for small ${H_2}$ fraction is at low densities, around $n \sim 10^5 {\rm cm}^{-3}$, in agreement with . If the ${H_2}$ fraction is increased, cooling dominates for a longer time, as expected, but not forever.

Our results were obtained for two possible values for the H$_2$ fraction: the value given by the simulations without DM, and the case of 100% molecular hydrogen. The fully molecular case produces the most effective cooling of the core and hence may underestimate the importance of DM heating. DM annihilation gives rise to additional electrons which could in principle produce more H$_2$ and go in the direction of cooling the star; however, for 100 GeV mass WIMPs, the DM only becomes important at such high gas densities that the core has already reached a 100% H$_2$ fraction and this effect is irrelevant. For lighter WIMPs, the transition to DM heating may occur at a lower H$_2$ fraction and a different baryon density.

The crossing point is stable against temperature changes at constant density and composition: if the temperature increases away from the critical value, cooling begins to dominate and drives the temperature back down to the critical value; if the temperature decreases, heating dominates and lowers the temperature. The reason for this behavior is that the H$_2$ cooling rate increases with temperature while the DM heating is independent of temperature.

As soon as the DM annihilation products are contained inside the protostellar core, the heating dominates over the cooling. Hence, for 100 GeV (1 GeV) neutralino DM, once the gas density reaches a critical value of $\sim 10^{13} {\rm cm}^{-3}$ ($10^9 {\rm cm}^{-3}$), the heating rate from DM annihilation exceeds the rate of hydrogen cooling. The protostellar core is *prevented from cooling and collapsing* further.

*Discussion.* Our main conclusion is that the standard picture of Pop III star formation is drastically modified by neutralino dark matter annihilation inside the protostellar object. The DM annihilation provides a heat source that exceeds any cooling mechanism and thereby hinders the further collapse of the protostar.

We propose that a new type of object is created, a "dark star" supported by DM annihilation rather than fusion. One question is how long such a phase of stellar evolution lasts. If such an object were stable for a long time period, it would even be possible for these dark stars to still exist today. Dark stars could last as long as the DM annihilation timescale, $\tau_e = m_\chi/(\rho_\chi \langle \sigma
v \rangle)$ $\simeq$ $0.6~ {\rm Gyr}$ $(n/10^{13} {\rm
cm}^{-3})^{-0.8}$ $(m_\chi/100 {\rm GeV})$ $(\sigma v/3 \times
10^{-26}{\rm cm}^3 {\rm s}^{-1})^{-1}$. For our canonical case, we find $\tau_e \sim 600$ Myr ($\sim 15$ Myr) for $n = 10^{13} {\rm
cm}^{-3}$ ($n = 10^{15} {\rm cm}^{-3}$). By comparison, the entire timescale for collapse (without taking into account DM annihilation) is $\sim 1$ Myr at $z=50$ or 100 Myr at $z=15$; even for the more recent episodes, the dynamical time at the high densities considered here is very short ($<10^3$ yr). However, after this DM annihilates away, it is possible that the DM hole in the small central core can fill in again, depending on the DM orbits at this stage. DM further out can also continue to heat the core. On the other hand, as baryons continue to accrete onto the protostar, it is possible that the annihilation shuts down sooner. The lifetime of the dark star phase is crucial to addressing the question of the effects it has on the universe.

Much further work is required to address the evolution of these objects; we here briefly speculate as to possible outcomes. One possibility is that the dark star phase still persists today, so that these objects never reach the main sequence. A second possibility would be a shorter dark star phase, during which the gas core is in a state of quasi-hydrostatic equilibrium. For our 100 GeV case, both the surface gravity and the hydrostatic force per unit mass from the temperature and pressure gradient are of the same order of magnitude, $\sim 10^{-3}$ cm/s$^2$. Outer material would continue to accrete onto the quasi-hydrostatic core , probably accompanied by the formation of an accretion shock. Alternatively, once the core's central temperature reaches $\sim 10^6$ K, deuterium burning and a pp chain may start. The star would then emit radiation and contract in a Kelvin-Helmholtz timescale, until the central density and pressure is high enough to start the CNO cycle. The star would finally reach the zero-age main sequence. In this scenario, Pop III star formation would be delayed rather than blocked. A third possibility would be that the core's contraction slows down as a consequence of DM heating, and yet the core continues to contract further. Then DM heating would continue to dominate over available cooling mechanisms at ever higher baryon densities.

The effects of a dark star phase of stellar evolution could be very interesting. The reionization of the IGM could be quite different, as would be the production of the heavy elements required to form all future generations of stars. We clearly need to recompute stellar structure with this new heat source, to see how different dark stars would look relative to the ordinary fusion-driven stars. It is possible, e.g., that dark stars are luminous but at lower temperatures, so that UV radiation does reionize the IGM yet the scenario is quite different. It is possible, e.g., that reionization is delayed. Perhaps the discrepancy between measurements of $\sigma_8$ by WMAP3 and in Lyman-$\alpha$ could be resolved in this way.

DM heating may also alter the mass of Pop III stars. Even without DM heating, the mass is still uncertain. Refs.  have explored different possibilities for the accretion process (disk vs. spherical) of baryonic matter onto the central protostar as well as sensitivity to cosmology and found great uncertainty as to the final stellar mass. DM heating could also affect the result. For example, in the case of spherical accretion, the heating might produce radiation at the Eddington luminosity whose pressure prevents further accretion. Thus, the DM heating might lead to lower mass stars than in the standard picture. Alternatively the initial protostellar object may be larger, and dark stars might accrete enough material to form large black holes en route to building the (as yet unexplained) $10^9 M_\odot$ black holes observed at $z \sim 6$.

What are other observational consequences of a "dark star"? If these objects are luminous but differ from ordinary stars (e.g. shine at lower temperatures), then James Webb Space Telescope could in principle find them (at z$\sim$`<!-- -->`{=html}10) and differentiate their spectra from those expected in the standard first star formation scenarios. In addition, one might hope to detect the DM annihilation products such as neutrinos and $\gamma$-rays. However, the angular resolution of current and planned detectors (AMANDA, ICECUBE, GLAST, HESS, VERITAS, MAGIC) is not good enough to identify an individual dark star source at $z> 10$, so that the $\nu$ and $\gamma$s would add to the extragalactic backgrounds and could provide limits at best. Remnants today of the dark star phase would be more testable, such as stellar remnants still in enhanced DM distributions. Today's remnant $10^6$ DM halos that were once the site of star formation may be modified due to the adiabatic contraction that took place earlier, so that enhanced DM annihilation might still occur today and these objects could be identified in upcoming experiments as individual $\nu$ or $\gamma$ sources.

**Acknowledgments.** This paper would not have come into existence without the help of Pierre Salati; we thank him for our conversations at the inception of this project. We are also extremely grateful to Chris McKee for his encouragement and for a series of discussions during the progress of this work. We thank A. Aguirre, P. Madau, F. Palla, J. Primack, R. Schneider, S. Stahler, G. Starkman, and N. Yoshida for discussions. K.F. acknowledges support from the DOE and MCTP via the Univ. of Michigan; from the Miller Inst. at UC Berkeley; and thanks the Physics Dept. at UCSC. D.S. and K.F. thank the Galileo Galilei Inst. in Florence, Italy, for support. P.G. acknowledges NSF grant PHY-0456825. D.S. acknowledges NSF grant AST-0507117 and GAANN.

[^1]: The details of the interactions and masses of the neutralinos depend on a large number of model parameters. In the minimal supergravity model, experimental and observational bounds restrict the neutralino mass $m_\chi$ to 50 GeV–2 TeV, while the annihilation cross section $\sigma v$ lies within an order of magnitude of $\langle \sigma v
    \rangle = 3 \times 10^{-26} {\rm cm^3/sec}$ (except at the low end of the mass range where it could be several orders of magnitude smaller). Nonthermal particles can have annihilation cross-sections that are many orders of magnitude larger (e.g. ) and would have even more drastic effects.


abstract: Motivation: In recent years, image-based biological assays have steadily become high-throughput, sparking a need for fast automated methods to extract biologically-meaningful information from hundreds of thousands of images. Taking inspiration from the success of ImageNet, we curate CytoImageNet, a large-scale dataset of openly-sourced and weakly-labeled microscopy images (890K images, 894 classes). Pretraining on CytoImageNet yields features that are competitive to ImageNet features on downstream microscopy classification tasks. We show evidence that CytoImageNet features capture information not available in ImageNet-trained features. The dataset is made available at <https://www.kaggle.com/stanleyhua/cytoimagenet>.
author: Stanley Bryan Z. Hua
Computer Science
University of Toronto
[`email@example.com`](mailto:firstname.lastname@example.com)
Alex X. Lu
Computer Science
University of Toronto
(Current Affiliation: Microsoft Research)
`firstname.lastname@example.com`
Alan M. Moses
Cell and Systems Biology
Computer Science
Center of Genome Evolution and Function
University of Toronto
`firstname.lastname@example.com`
title: CytoImageNet: A large-scale pretraining dataset for bioimage transfer learning

# Introduction

Automated microscopy has given biologists a tool for studying various cell types under varying experimental conditions with minimal assistance. In recent years, automated image analysis has proven useful for drug discovery and functional genomics research; identifying the function of novel genes and cellular response to various drugs \[1\]. With the capacity to perform large-scale screens rapidly, tens of thousands of images are being collected daily, demanding computational methods to extract biology from these images \[2\].

Methods for automated image analysis remains an active area of research. One strategy for automating the analysis of microscopy images is to use low-level features from classic computer vision strategies; extracting thousands of general features related to area, shape, texture and intensity of single-cells \[3, 4\]. However, these features usually require extensive engineering, parameter tuning, and effective single-cell segmentation, all of which must be customized to the images of interest. Furthermore, these features may not capture relevant biological signal and instead capture experimental noise \[5\]. To reduce labor and improve the sensitivity of extracted features, some studies rely on deep learning to automatically learn features. These strategies range from supervised training with labels, autoencoders to more recently, self-supervised learning \[5, 6, 7\]. However, these models require users to train on their own data, which is computationally expensive, time-consuming, and requires specialized hardware like GPUs.

As an alternative to training such models, some studies instead rely on transfer learning. Most commonly, studies re-use features from neural networks trained to classify images from ImageNet, a large-scale natural image classification dataset \[8, 9, 10\]. By training models to perform well on large datasets, they can learn useful features that may be transferred to other datasets.

ImageNet is a large-scale dataset of more than 14 million diverse and annotated natural images \[11\]. A defining factor in the success of ImageNet is their diversity of image labels and image sources. Models that perform the best on ImageNet-1K become the standard for transfer learning \[12\]. These models learn a diverse set of features that are able to capture complex and generalizable patterns inherent in their dataset. Surprisingly, ImageNet features have even shown good transferability on microscopy classification tasks \[13, 14\].

However, a recent study suggests that "ImageNet features are less general than previously suggested," and that pretraining on datasets with a closer domain match to the downstream task benefits transfer learning \[12\]. In fact, Caicedo et al. showed that pretraining convolutional networks (CNN) on weakly labeled microscopy images yields features that outperform ImageNet features on downstream classification tasks \[5\]. We hypothesize that features from pretraining on a diverse large-scale microscopy dataset will yield more domain-relevant features than ImageNet.

We present CytoImageNet, a large-scale image dataset of weakly labeled microscopy images (890K images, 894 classes). Incorporating image data and labels from 40 distinct microscopy datasets, CytoImageNet mimics the diverse and complex nature of ImageNet. We provide a framework for curating similar large-scale pretraining datasets using openly available microscopy datasets.

We believe that pretraining on CytoImageNet may yield biologically-meaningful image representations useful in any downstream biological task. Since concatenation of CytoImageNet-trained features and ImageNet-trained features performs significantly better than either one alone, we believe that CytoImageNet pretraining provides features that are useful, and different from ImageNet features.

# Methods

CytoImageNet is a dataset of 890,737 microscopy images sourced from 40 openly available datasets. The selected datasets originated from 5 databases: (1) Recursion, (2) Image Data Resource (IDR) \[15\], (3) Broad Bioimage Benchmark Collection (BBBC) \[16\], (4) Kaggle and (5) Cell Image Library (CIL). Electron microscopy and histopathology images were excluded. Code for the methods below is available at <https://github.com/stan-hua/CytoImageNet>. The dataset is available at <https://www.kaggle.com/stanleyhua/cytoimagenet>.

| **Database**                        | **Number of Labels Contributed** |
|:------------------------------------|:--------------------------------:|
| Recursion                           |               651                |
| Image Data Resource                 |               450                |
| Broad Bioimage Benchmark Collection |               202                |
| Kaggle                              |                27                |
| Cell Image Library                  |                1                 |

Number of labels contributed by each image database

## Weak label assignment & stratified downsampling

As a training target for CytoImageNet, we sought to associate each microscopy image with a class label. As expert-assigned annotations are not readily available or standardized across publicly available image datasets, we instead chose to assign images to classes based upon metadata associated with each of the image datasets we used. We considered 7 kinds of metadata: (organism, cell_type, cell_visible, phenotype, compound, gene, and sirna). Each column was searched for unique values and their counts to find potential weak labels. To create 894 labels, unique values with at least 287 images served as potential labels. Beginning with the least counts, each potential label was iterated through, searching for images that can be assigned to the label. Images that were already assigned labels are excluded, and their unique image index is stored in a hash table pointing to their label. To improve the diversity of images selected for a label, stratified sampling was done by sampling from images that were grouped by values in the remaining potential label columns.

## Image standardization & upsampling

To preprocess images, we converted all images to PNG. RGB images are converted to grayscale by averaging over channels. We normalized each channel between \[0, 255\] using the 0.1 and 99.9th percentile pixel intensity as the lower and upper ranges respectively \[6\]. Finally, for fluorescent images, we merged the channels.

We took 1-4 crops from each image of varying spatial scales where possible. Unless the image was smaller than 70x70 pixels, we cropped images into four quadrants. To introduce variations in spatial scale, we further cropped these quadrants by sampling crops of either full, 0.5, 0.25, or 0.125 of the quadrant (excluding any crop scales that would produce images smaller than 70x70 pixels). Finally, we filtered any crops that were all black, all white, or where the 75th percentile pixel intensity is equal to 0. Crop pixel intensities were re-normalized as we previously did for full images.

## Training

CytoImageNet is split into a training and validation set with 10% used for validation. This yields roughly 900 training samples for each label. Images are fed in batches of 64 with random 0 to 360 degrees rotations. We train convolutional networks (EfficientNetB0) to classify one of the 894 labels, by minimizing the categorical cross-entropy loss of predictions to ground truth labels. Randomly initialized models were trained for 24 epochs (2 weeks) on an NVIDIA Tesla K40C. The loss was optimized via the Adam optimizer with learning rate of 0.001.

## Evaluation

We compare features extracted from a CytoImageNet-pretrained model, an ImageNet-pretrained model, and a randomly-initialized model on three distinct downstream tasks: 1) BBBC021 mechanism-of-action classification, 2) Cells-Out-Of-Sample (COOS-7) protein localization classification, and 3) a CYCLoPs dataset for protein localization classification \[17, 18, 19\]. Additionally, we would like to know if CytoImageNet and ImageNet pretraining lead to different learnt features. If so, fusing the extracted features from both may increase performance on downstream tasks. To test this, we concatenate CytoImageNet and ImageNet extracted features for each image in the transfer tasks.

BBBC021 is a dataset of fully imaged human cells. Cells are treated with one of 113 small molecules at 8 concentrations, and fluorescent images are captured staining for nucleus, actin and microtubules. The phenotypic profiling problem is presented, where the goal is to extract features containing meaningful information about the cellular phenotype exhibited. Each of 103 unique compound-concentration treatment is labeled with a mechanism-of-action (MOA). The MOA is predicted for each unique treatment (averaging features over all treatment examples) by matching the MOA of the closest point excluding points of the same compound. We use a 1-nearest neighbors using cosine distance as the metric and report the not-same-compound (NSC) accuracy.

On the other hand, the WT2 dataset from the CYCLoPs database consists of 27,058 single-cell images of yeast cells. The task is to classify the subcellular localization of a fluoresced protein, given two channels staining for the protein of interest and the cytosol. An 11-nearest neighbors is used to perform leave-one-out classification, and the accuracy is reported.

Lastly, COOS-7 contains 132,209 single-cell images of mouse cells. Images are spread over 1 training set and 4 testing sets, where each single-cell image contains a protein and nucleus fluorescent channels. COOS-7 was curated, such that each succeeding testing set has a greater degree of covariate shift from the training set. Similar to the CYCLoPs dataset, the COOS-7 dataset presents a protein subcellular localization classification task. We use a 11-nearest neighbors fit on the training set to predict on the testing sets and report the accuracy.

We chose these datasets (1) to assess the generalizability of CytoImageNet-pretrained models on distinct biological image datasets and problems, and (2) to compare against ImageNet features that have already been shown to perform relatively well on these datasets \[13, 7, 19\].

Features are extracted from the penultimate layer of EfficientNetB0 models. ImageNet-pretrained EfficietNetB0 weights are loaded in using Keras, but the weights themselves belong to the source paper \[20\]. Predictions are generated using k-Nearest Neighbors (kNN), and accuracy is reported with a 95% confidence interval. BBBC021 kNN is implemented with scikit-learn, while kNN is implemented using the faiss library for COOS-7 and CYCLoPs datasets.

We compare the performance of features extracted following two optional preprocessing methods: normalization and channel merging. Normalization refers to channel normalization with the lower bound and upper bound of the 0.1 and 99.9th percentile pixel intensity respectively \[6\]. Channel merging refers to the merging of fluorescent channels into a single grayscale image before extracting features. If not, features are extracted for each fluorescent channel then concatenated.

|  | Random | ImageNet | CytoImageNet | Fusion |
|:--:|:--:|:--:|:--:|:--:|
| BBC021 | 27.18 $\pm$ 8.59% | **83.5** $\pm$ 7.17% | **83.5** $\pm$ 7.17% | **86.41** $\pm$ 6.62% |
| CYCLoPs | 53.06 $\pm$ 0.59% | **68.47** $\pm$ 0.55% | 65.19 $\pm$ 0.57% | **77.97** $\pm$ 0.49% |
| COOS7 Test Set 1 | 65.68 $\pm$ 0.91% | **88.87** $\pm$ 0.61% | 88.58 $\pm$ 0.61% | **94.80** $\pm$ 0.43% |
| COOS7 Test Set 2 | 67.81 $\pm$ 0.7% | 88.93 $\pm$ 0.47% | **89.37** $\pm$ 0.46% | **95.07** $\pm$ 0.33% |
| COOS7 Test Set 3 | 48.77 $\pm$ 0.54% | **75.91** $\pm$ 0.46% | 65.97 $\pm$ 0.51% | **78.97** $\pm$ 0.44% |
| COOS7 Test Set 4 | 51.88 $\pm$ 0.56% | **82.19** $\pm$ 0.43% | 78.52 $\pm$ 0.46% | **87.47** $\pm$ 0.37% |

kNN classification accuracy (with 95% confidence intervals) on downstream tasks, comparing EfficientNetB0 penultimate layer features. 'Random' refers to the model with randomly-initialized weights. 'Fusion' refers to the concatenation of ImageNet features and CytoImageNet features. Concatenation of randomly-initialized features with ImageNet features gave no increase in performance (data not shown).

# Discussion

CytoImageNet pretraining shows comparable performance to ImageNet pretraining on all downstream tasks. Normalizing pixel intensities for each channel and extracting features for each channel led to the best performance for all features, implying that automated methods for feature extraction still require some form of image preprocessing.

Given that CytoImageNet has a closer domain match to the downstream microscopy tasks, it is surprising that CytoImageNet features alone do not outperform ImageNet features. A recent study on ImageNet pretraining reported a strong correlation between ImageNet validation accuracy and transfer accuracy \[12\]. Our trained EfficientNetB0 model only achieved 13.42% accuracy on the training set and more importantly, 11.32% on the CytoImageNet validation set, where 402 of the 894 labels had no images correctly predicted. Of which, 375 of the 402 were compound labels. Meanwhile, EfficientNetB0 was shown to attain a top-1 accuracy of 77.1% on ImageNet-1K \[20\]. Given our low validation accuracy, we speculate that hyperparameter optimization will lead to further increases in benchmark classification accuracy. We believe this is an important area for future research.

Encouragingly, fusing extracted CytoImageNet and ImageNet features (via concatenation) performs the best on downstream tasks by a significant margin. This implies that CytoImageNet and ImageNet pretraining result in the learning of different albeit meaningful image representations. We believe this is due to CytoImageNet-trained models learning some feature of microscopy images that are not available in natural images. These results suggest that future applications of bioimage transfer learning may benefit the most from the fusion of CytoImageNet and ImageNet features.

In summary, we present CytoImageNet; a pretraining dataset of 890,737 microscopy images and 894 classes. We contribute a simple framework for combining openly available microscopy datasets, in the hopes that this sparks interest in the vast amount of rich microscopy data available online. Furthermore, we hope that this may inspire others to experiment on their own; both to create new large-scale datasets and to develop image classification architectures and training procedures that generalize to other image domains.

# Acknowledgements

Stanley Z. Hua was funded by the University of Toronto Cell & Systems Biology undergraduate research award. Alan M. Moses holds a Tier II Canada Research Chair. The authors would like to acknowledge Recursion, Image Data Resource, Broad Bioimage Benchmark Collection and Kaggle for allowing open access to their microscopy datasets, which were used to construct the CytoImageNet dataset. We would like to thank the authors responsible for uploading their data on these databases. We would also like to thank Juan Caicedo and Anne Carpenter for their valuable input regarding the Broad Institute datasets. Finally, the authors would like to acknowledge the Canadian Foundation for Innovation that sponsored the NVIDIA GPU with which the models were trained on.

# References

\[1\] Caicedo JC, Singh S, Carpenter AE. Applications in image-based profiling of perturbations. Current Opinion in Biotechnology. 2016;39:134–142. doi:10.1016/j.copbio.2016.04.003

\[2\] Wollman R, Stuurman N. High throughput microscopy: From raw images to discoveries. Journal of Cell Science. 2007;120(21):3715–3722. doi:10.1242/jcs.013623

\[3\] Danuser G. Computer Vision in cell biology. Cell. 2011;147(5):973–978. doi:10.1016/j.cell.2011.11.001

\[4\] Carpenter AE, Jones TR, Lamprecht MR, Clarke C, Kang I, Friman O, Guertin DA, Chang J, Lindquist RA, Moffat J, et al. CellProfiler: image analysis software for identifying and quantifying cell phenotypes. Genome Biology. 2006;7(10). doi:10.1186/gb-2006-7-10-r100

\[5\] Caicedo JC, McQuin C, Goodman A, Singh S, Carpenter AE. Weakly supervised learning of single-cell feature embeddings. 2018 IEEE/CVF Conference on Computer Vision and Pattern Recognition. 2018. doi:10.1109/cvpr.2018.00970

\[6\] Kraus OZ, Grys BT, Ba J, Chong Y, Frey BJ, Boone C, Andrews BJ. Automated Analysis of high‐content microscopy data with Deep Learning. Molecular Systems Biology. 2017;13(4):924. doi:10.15252/msb.20177551

\[7\] Lu AX, Kraus OZ, Cooper S, Moses AM. Learning unsupervised feature representations for single cell microscopy images with paired cell inpainting. PLOS Computational Biology. 2019;15(9). doi:10.1371/journal.pcbi.1007348

\[8\] Huh M, Agrawal P, Efros AA. What makes ImageNet good for transfer learning? arXiv. 2016;1608.08614.

\[9\] Razavian AS, Azizpour H, Sullivan J, and Carlsson S. CNN features off-the-shelf: an astounding baseline for recognition. In IEEE Conference on Computer Vision and Pattern Recognition Workshops (CVPRW),pages 512–519. IEEE, 2014.

\[10\] Donahue J, Jia Y, Vinyals O, Hoffman J, Zhang N, Tzeng E, and Darrell T. Decaf: A deep convolutional activation feature for generic visual recognition. In International Conference on Machine Learning, pages 647–655, 2014.

\[11\] Deng J, Dong W, Socher R, Li L-J, Kai Li, Li Fei-Fei. ImageNet: A large-scale hierarchical image database. 2009 IEEE Conference on Computer Vision and Pattern Recognition. 2009. doi:10.1109/cvpr.2009.5206848

\[12\] Kornblith S, Shlens J, Le QV. Do Better Imagenet Models Transfer Better? 2019 IEEE/CVF Conference on Computer Vision and Pattern Recognition (CVPR). 2019. doi:10.1109/cvpr.2019.00277

\[13\] Pawlowski N, Caicedo JC, Singh S, Carpenter AE, Storkey A. Automating morphological profiling with generic deep convolutional networks. bioRxiv. 2016. doi:10.1101/085118

\[14\] Kensert A, Harrison PJ, Spjuth O. Transfer learning with deep convolutional neural networks for classifying cellular morphological changes. SLAS DISCOVERY: Advancing the Science of Drug Discovery. 2019;24(4):466–475. doi:10.1177/2472555218818756

\[15\] Williams E, Moore J, Li SW, Rustici G, Tarkowska A, Chessel A, Leo S, Antal B, Ferguson RK, Sarkans U, et al. Image Data Resource: A Bioimage Data Integration and publication platform. Nature Methods. 2017;14(8):775–781. doi:10.1038/nmeth.4326

\[16\] Ljosa V, Sokolnicki KL, Carpenter AE. Annotated high-throughput microscopy image sets for validation. Nature Methods. 2012;9(7):637–637. doi:10.1038/nmeth.2083

\[17\] Caie PD, Walls RE, Ingleston-Orme A, Daya S, Houslay T, Eagle R, Roberts ME, Carragher NO. High-content phenotypic profiling of drug response signatures across distinct cancer cells. Molecular Cancer Therapeutics. 2010;9(6):1913–1926. doi:10.1158/1535-7163.mct-09-1148

\[18\] Koh JL, Chong YT, Friesen H, Moses A, Boone C, Andrews BJ, Moffat J. Cyclops: A comprehensive database constructed from automated analysis of protein abundance and subcellular localization patterns in saccharomyces cerevisiae. G3 Genes\|Genomes\|Genetics. 2015;5(6):1223–1232. doi:10.1534/g3.115.017830

\[19\] Lu AX, Lu AX, Schormann W, Andrews DW, Moses AM. The Cells Out of Sample (COOS) dataset and benchmarks for measuring out-of-sample generalization of image classifiers. ArXiv. 2019;abs/1906.07282.

\[20\] Tan M, Le Q. EfficientNet: Rethinking Model Scaling for Convolutional Neural Networks. In: 36th International Conference on Machine Learning (ICML 2019) Long Beach, California, USA, 9-15 June 2019. Red Hook, NY: Curran Associates, Inc; 2019.

\[21\] Zhang H, Cao X, Tang M, Zhong G, Si Y, Li H, Zhu F, Liao Q, Li L, Zhao J, et al. A subcellular map of the human kinome. eLife. 2021;10. doi:10.7554/elife.64943

# Supplementary Material

```latex
\resizebox{\textwidth}{!}{%
\begin{tabular}{@{}llll@{}}
\toprule
\textbf{Database} & \textbf{Shorthand} & \textbf{Dataset Name}                                & \textbf{Link}                                              \\ \midrule
BBBC              & bbbc012            & C. elegans infection marker                          & https://bbbc.broadinstitute.org/BBBC012                    \\ \midrule
BBBC              & bbbc010            & C. elegans live/dead assay                           & https://bbbc.broadinstitute.org/BBBC010                    \\ \midrule
BBBC              & bbbc011            & C. elegans metabolism assay                          & https://bbbc.broadinstitute.org/BBBC011                    \\ \midrule
BBBC              & bbbc048            & Cell Cycle Jurkat Cells                              & https://bbbc.broadinstitute.org/BBBC048                    \\ \midrule
BBBC              & bbbc017            & Human HT29 colon-cancer cells shRNAi screen          & https://bbbc.broadinstitute.org/BBBC017                    \\ \midrule
BBBC   & bbbc026          & Human Hepatocyte and Murine Fibroblast cells - Co-culture experiment & https://bbbc.broadinstitute.org/BBBC026                                \\ \midrule
BBBC              & bbbc006            & Human U2OS cells (out of focus)                      & https://bbbc.broadinstitute.org/BBBC006                    \\ \midrule
BBBC              & bbbc022            & Human U2OS cells - compound cell-painting experiment & https://bbbc.broadinstitute.org/BBBC022                    \\ \midrule
BBBC              & bbbc013            & Human U2OS cells cytoplasm-nucleus translocation     & https://bbbc.broadinstitute.org/BBBC013                    \\ \midrule
BBBC              & bbbc014            & Human U2OS cells cytoplasm-nucleus translocation (2) & https://bbbc.broadinstitute.org/BBBC014                    \\ \midrule
BBBC              & bbbc015            & Human U2OS cells transfluor                          & https://bbbc.broadinstitute.org/BBBC015                    \\ \midrule
BBBC              & bbbc025            & Human U2OS cells - RNAi Cell Painting experiment     & https://bbbc.broadinstitute.org/BBBC025                    \\ \midrule
BBBC              & bbbc045            & Human White Blood Cells                              & https://bbbc.broadinstitute.org/BBBC045                    \\ \midrule
BBBC              & bbbc051            & Human kidney cortex cells                            & https://bbbc.broadinstitute.org/BBBC051                    \\ \midrule
BBBC              & bbbc038            & Kaggle 2018 Data Science Bowl                        & https://bbbc.broadinstitute.org/BBBC038                    \\ \midrule
BBBC              & bbbc020            & Murine bone-marrow derived macrophages               & https://bbbc.broadinstitute.org/BBBC020                    \\ \midrule
BBBC              & bbbc041            & P. vivax (malaria) infected human blood smears       & https://bbbc.broadinstitute.org/BBBC041                    \\ \midrule
BBBC              & bbbc005            & Synthetic cells                                      & https://bbbc.broadinstitute.org/BBBC005                    \\ \midrule
CIL               & cil\_kinome        & Kinome Atlas                                         & http://cellimagelibrary.org/pages/kinome\_atlas            \\ \midrule
IDR               & idr0081            & Adenovirus                                           & https://idr.openmicroscopy.org/webclient/?show=screen-2406 \\ \midrule
IDR               & idr0017            & Chemical-Genetic Interaction Map                     & https://idr.openmicroscopy.org/webclient/?show=screen-1151 \\ \midrule
IDR               & idr0016            & Compound Profiling                                   & https://idr.openmicroscopy.org/webclient/?show=screen-1251 \\ \midrule
IDR               & idr0009            & Early Secretory Pathway                              & https://idr.openmicroscopy.org/webclient/?show=screen-803  \\ \midrule
IDR               & idr0041            & Mitotic Atlas                                        & https://idr.openmicroscopy.org/webclient/?show=project-404 \\ \midrule
IDR               & idr0021            & Pericentriolar Material                              & https://idr.openmicroscopy.org/webclient/?show=project-51  \\ \midrule
IDR               & idr0080            & Perturbation                                         & https://idr.openmicroscopy.org/webclient/?show=screen-2701 \\ \midrule
IDR               & idr0088            & Phenomic Profiling                                   & https://idr.openmicroscopy.org/webclient/?show=screen-2651 \\ \midrule
IDR               & idr0003            & Plasticity                                           & https://idr.openmicroscopy.org/webclient/?show=screen-51   \\ \midrule
IDR               & idr0072            & Subcellular Localization                             & https://idr.openmicroscopy.org/webclient/?show=screen-2952 \\ \midrule
IDR               & idr0037            & Variation in Human iPSC lines                        & https://idr.openmicroscopy.org/webclient/?show=screen-2051 \\ \midrule
IDR               & idr0067            & Yeast Meiosis                                        & https://idr.openmicroscopy.org/webclient/?show=project-904 \\ \midrule
IDR    & idr0093          & siRNA screen for cell size and RNA production                        & https://idr.openmicroscopy.org/webclient/?show=screen-2751             \\ \midrule
Kaggle & kag\_cell\_cycle & Cell Cycle Experiments                                               & https://www.kaggle.com/paultimothymooney/cell-cycle-experiments        \\ \midrule
Kaggle & kag\_hpa\_single & Human Protein Atlas - Single Cell Classification                     & https://www.kaggle.com/c/hpa-single-cell-image-classification/data     \\ \midrule
Kaggle & kag\_hpa         & Human Protein Atlas Image Classification                             & https://www.kaggle.com/c/human-protein-atlas-image-classification/data \\ \midrule
Kaggle            & kag\_leukemia      & Leukemia Classification                              & https://www.kaggle.com/andrewmvd/leukemia-classification   \\ \midrule
Recursion         & rec\_rxrx1         & RxRx1                                                & https://www.rxrx.ai/rxrx1\#the-data                        \\ \midrule
Recursion         & rec\_rxrx19a       & RxRx19a                                              & https://www.rxrx.ai/rxrx19a                                \\ \midrule
Recursion         & rec\_rxrx19b       & RxRx19b                                              & https://www.rxrx.ai/rxrx19b                                \\ \midrule
Recursion         & rec\_rxrx2         & RxRx2                                                & https://www.rxrx.ai/rxrx2                                  \\ \bottomrule
\end{tabular}%
}
```


abstract: Supermassive black holes appear to be generic components of galactic nuclei. Following their formation in the early universe, black holes should often find themselves in bound pairs as a consequence of galaxy mergers. The greatest uncertainty in estimating the coalescence time of black hole binaries is the degree to which a binary wanders about the center of the galactic potential. A simple model for binary decay is presented which qualitatively reproduces the evolution observed in the $N$-body simulations. The model predicts binary coalescence times that are never less than several billion years. Mass ejection by a decaying black hole binary should substantially lower the density of its host nucleus. The weak density cusps of bright ellipticals may be explained in this way if it is assumed that these galaxies formed from nearly gas-free mergers.
author: David Merritt
title: Black Holes and Galaxy Evolution

# Introduction

The case for supermassive black holes (BHs) in galactic nuclei is now very strong, and in a few galaxies almost irrefutable. The evidence is most compelling in galaxies where the kinematics of stars or gas can be traced into sub-parsec scales; examples are the Milky Way (Genzel et al. 1997; Ghez et al. 1998), where stellar proper motions can be measured at distances of $\mathrel{\mathpalette\fun <}0.02$ pc, and NGC 4258 (Miyoshi et al. 1995), where the rotation curve of H$_2$O maser sources obeys Kepler's law into $\sim 0.13$ pc. For a larger sample of galaxies, kinematical data extending into $\sim 10$ or $\sim 100$ pc provide suggestive though not yet irrefutable evidence of supermassive compact objects (Kormendy & Richstone 1995).

Although the masses of the detected BHs comprise on average only $\sim 0.3\%$ of the mass of their host spheroids (Ho 1998), there is a growing body of work suggesting that the dynamical influence of a supermassive BH can extend far beyond the nucleus. Furthermore the formation and growth of BHs may be intimately connected with the evolution of galaxies on larger scales. For instance, mergers between galaxies containing nuclear BHs would produce supermassive BH binaries which would eventually coalesce via the emission of gravitational radiation (Begelman, Blandford & Rees 1980). The formation and decay of these binaries may be relevant to a wide range of phenomena, from the wiggling of radio jets (Kaastra & Roos 1992) to the destruction of stellar density cusps (Ebisuzaki et al. 1991).

This review discusses the time scale for coalescence of supermassive BH binaries (§2) and the dynamical effect of the BHs on stellar nuclei (§3) and the large-scale structure of galaxies (§4).

# Supermassive Black Hole Binaries

Supermassive BHs appear to be ubiquitous components of galactic nuclei. Since galaxies often merge, one would expect to form BH binaries at a rate that is roughly equal to the galaxy merger rate. The formation and coalescence of a BH binary is believed to take place via three, fairly distinct stages (Begelman, Blandford & Rees 1980):

1\. Two parent galaxies interact; the BHs – surrounded by their dense star clusters – sink to the center of the common potential well via dynamical friction, forming a binary.

2\. The BH binary shrinks by ejecting stars from the nucleus via three-body interactions.

3\. When the binary separation has fallen to the value at which gravitational radiation becomes efficient, the BHs coalesce.

In the absence of gas, the rates of the first two processes are in principle possible to calculate via $N$-body simulations. However the problem is a difficult one to simulate due to the wide range of length and time scales; none of the published simulations have succeeded in following the decay to the point where gravitational radiation would become important. Here we derive a simple analytical model for the decay of BH binaries which reproduces the results of the $N$-body simulations, then use it to estimate the coalescence time.

Consider a binary with total mass $M_{12}=m_1+m_2$, $m_1\geq m_2$ and semimajor axis $a(t)$. We assume that the background galaxy is initially a singular isothermal sphere, with density and mass profiles $$\rho_0(r) = {\sigma^2\over 2\pi G r^2}, \ \ \ \ M_0(r) = {2\sigma^2r\over
G};$$ $\sigma$ is the one-dimensional stellar velocity dispersion. The binary becomes "hard" when $a$ falls below $\sim a_h=Gm_2/4\sigma^2$ (Quinlan 1996); subsequent evolution is driven by the capture and ejection of stars that interact with the binary. The hardening rate $$H = {\sigma\over G\rho}{d\over dt}\left({1\over a}\right)$$ depends only weakly on $m_1/m_2$ and on $a$ for $a\mathrel{\mathpalette\fun <}a_h$; in the limit $a\ll a_h$, $H\approx 16$ (Hills 1992; Mikkola & Valtonen 1992; Quinlan 1996). Mass ejection occurs at a rate $$J = {1\over M_{12}}{dM_{ej}\over d\ln(1/a)}$$ where $J\approx 1$ is again nearly independent of ($m_1/m_2$, $a$) for $a\ll a_h$ (Quinlan 1996). Thus $$M_{ej} \approx J M_{12} \ln(a_h/a).$$

If the binary remained fixed at the center of a spherical potential, it would quickly eject all stars on orbits with pericenters less than $\sim a$. The local density would fall to zero and the binary would cease to harden. We assume instead that mass removal produces a constant-density core of stars out to some radius $r_c(t)$, where $$M_{ej} = M_0(r_c) - {4\pi\over 3}\rho_0(r_c)r_c^3 = {4\over
3}{\sigma^2 r_c\over G}.$$ In order to achieve this, the binary must interact with stars at radii $r\mathrel{\mathpalette\fun <}r_c\gg a$. This may seem an optimistic assumption, but it turns out that something very similar happens in the $N$-body simulations, as discussed below.

The binary separation rapidly falls below $r_c$ and subsequent evolution occurs against an evolving background density $\rho(t) = \sigma^2/2\pi G r_c^2(t)$. Combining equations (2), (4) and (5), we find for the hardening rate $${d(a_h/a)\over d(t/t_0)} = {1\over\ln^2(a_h/a)},\ \ \ \ \
t_0={9\pi J^2\over H}{\left(M_{12}\over 2m_2\right)}
{\left(GM_{12}\over\sigma^3\right)}.$$ The logarithm represents the decrease in the hardening rate due to the declining stellar density. Integrating with respect to time, $${t-t_h\over t_0} = {a_h\over a}\left[\ln^2\left({a_h\over
a}\right) -
\ln\left({a_h\over a}\right)^2 + 2\left(1-{a\over
a_h}\right)\right]$$ where $a(t_h)=a_h$. The functions $t(a)$ and $a(t)$ may be approximated as $${t-t_h\over t_0} \approx 0.6\left({a_h\over a}\right)\ln^2
\left({a_h\over a}\right),\ \ \ \ \ \
{a_h\over a} \approx
{4(t-t_h)/t_0\over\ln^2\left[(t-t_h)/t_0\right]}$$ for $t-t_h\mathrel{\mathpalette\fun >}$ a few $\times\ t_0$.

We may compare the predictions of this simple model to the results of $N$-body simulations of BH binary decay. Quinlan & Hernquist (1997) (QH) followed the evolution of a BH binary in a galaxy with an initial $\rho\propto r^{-2}$ density cusp. Expressed in their units, the scaling parameters in eq. (7) are $a_h=1/200$, $t_0\approx 1/10$ and $t_h\approx 8$. Fig. 1 is a plot of eq. (7) scaled to these units. The predicted time dependence $a(t)$ is qualitatively very similar to that found by QH (their Fig. 2). However, those authors found that the rate of binary decay depended systematically on the number $N$ of particles used to represent the galaxy: larger values for $N$ gave slower evolution. For $N=6250$, the smallest number considered by them, Fig. 1 overestimates the evolution rate by a factor of only $\sim 2$, remarkably close given the approximations made. For $N=2.5\times 10^4$ and $10^5$, the adjustment factors are $\sim 3$ and $\sim 5$ respectively.

QH attributed the $N$-dependence of the hardening rate to wandering of the binary about the center of the galaxy potential. In a core of density $\rho$, an rms velocity $v_b$ of the binary's center of mass will produce a wandering with amplitude $r_w \sim v_b/(G\rho)^{1/2}$. If the source of the binary's motion were elastic encounters with stars of mass $m_*$, equipartition arguments would give $v_b/\sigma\approx (m_*/M_{12})^{1/2}\approx 10^{-4}$. However the binary converts some of its binding energy into bulk motion when it ejects stars, with momenta $\sim
m_*(a_h/a)^{1/2}\sigma$. This "super-elastic scattering" allows the binary to sample the stellar density over a much wider range of radii than if it remained stationary; QH in fact verified that fixing the binary to the center of the potential caused the hardening to cease after the binary had ejected most of the stars within $r=a_h$.

If the background density within the wandering radius $r_w$ were to fall, the gravitational restoring force acting on the binary would also drop and the wandering amplitude would increase, roughly in proportion to $\rho^{-1/2}$. This feedback mechanism should tend to inhibit the formation of a central "hole" in the stellar density (Mikkola & Valtonen 1992). The crude model derived above mimics this process by forcing the background density to remain constant for $r< r_c$; the model could be improved by incorporating the wandering radius explicitly.

In the $N$-body simulations, the perturbations that the binary experiences from the "stars" become less noisy as $N$ is increased, implying a reduced wandering radius and a lower hardening rate. Makino (1997) found that the wandering amplitude varied roughly as $N^{-1/2}$ for $N<3\times 10^5$. However QH found that the hardening rate did not change appreciably when $N$ was increased from $10^5$ to $2\times 10^5$. They argued that $r_w$ might reach a limiting amplitude for large $N$ due to the feedback mechanism discussed above. This hypothesis is important to check; however doing so will require simulations with very large particle numbers.

The binary will rapidly coalesce when the time scale for emission of gravitational radiation $$t_{gr} = {5\over 256}{c^5a^4\over G^3 m_1m_2M_{12}}$$ (Peters 1964, for a circular binary) equals the hardening time $|a/\dot a|$. Using eq. (6), this occurs when $a=a_{gr}$ where $${\left(a_{gr}/a_h\right)^5\over \ln^2\left(a_{gr}/
a_h\right)} = {9\pi\times 16^5 J^2\over 20H}
\left({m_1\over m_2}\right) \left({M_{12}\over
2m_2}\right)^2\left({\sigma\over c}\right)^5$$ or $$a_{gr}/a_h \approx A |\ln A |^{0.4},\ \ \ \ \ \
 A = 9.85\left({m_1\over m_2}\right)^{0.2} \left({M_{12}\over
2m_2}\right)^{0.4} \left({\sigma\over c}\right).$$ For $\sigma/c = 300/(3\times 10^5) \approx 0.001$, $A\approx 0.01$ and $$a_{gr}/a_h\approx 0.018,$$ i.e. the binary must shrink by a factor of $\sim 50$ beyond the hardening radius for coalescence to occur. The simulations of QH followed the decay only over a factor of $\sim 12$ beyond $a_h$.

We can use the analytic expressions derived here to extend the $N$-body simulations until a time corresponding to coalescence; the correction factor of 5 mentioned above is applied to $(a/a_h)$ in order to reproduce the hardening rate observed in the largest-$N$ simulations. Fig. 2a shows the time from formation of a hard binary until coalescence, $\Delta t_1$, as a function of $\sigma$ and $m_2$ for a binary with $m_2/m_1=1/2$. An approximate expression is $$\Delta t_1 \approx 8t_0 A^{-1}|\ln A|^{8/5};$$ if we further approximate the logarithm as a constant, $|\ln A|\approx |\ln
0.01| \approx 4.61$, then $$\Delta t_1 \approx 1.4\times 10^{10} {\rm yr} \left({m_2\over m_1}\right)^{0.2}
\left({M_{12}\over 2m_2}\right)^{0.6} \left({M_{12}\over
10^9M_{\odot}}\right) \left({\sigma\over 200 {\rm \ km\ s}^{-1}}\right)^{-4}.$$ While eq. (14) is an accurate extrapolation of the evolution seen in the $N$-body experiments, it should probably be interpreted as a lower limit to the coalescence time expected in a real spherical galaxy due to the likely overestimate of the wandering radius in the simulations. However the dependence of $\Delta t_1$ on $m_2$ and $\sigma$ is probably more robust. Since for real galaxies $L\sim \sigma^4$ and $M_{12}\sim L$, eq. (14) predicts a coalescence time that is almost independent of $\sigma$ and $M_{12}$. The solid line in Fig. 2a is the "Faber-Jackson law" for BHs, $${M_{\bullet}\over 10^8M_{\odot}} \approx 3.1
\left({\sigma\over 200{\rm \ km\ s^{-1}}}\right)^4;$$ this relation is consistent with three galaxies whose BH masses are well-determined: M84 (Bower et al. 1998), M87 (Macchetto et al. 1997) and NGC 4258 (Miyoshi et al. 1995). Fig. 2a suggests a typical coalescence time of several billion years.

Very low mass binaries would evolve more rapidly (eq. 14), but another factor would delay their coalescence: the necessity of the two BHs finding their way to the potential center following a galaxy merger. Governato et al. (1994) carried out merger simulations of galaxies containing BHs; they concluded that inspiral times could be very long if the infalling galaxy had a sufficiently low mean density to be tidally disrupted. A more likely scenario is that the infalling BH retains some of the mass of its host galaxy. We can estimate the inspiral time by assuming that the smaller galaxy is tidally limited by the larger galaxy as its orbit decays. Following Merritt (1984), the tidal radius is located at the point where the effective potential (gravitational plus centrifugal) has a saddle point. Approximating both galaxies as singular isothermal spheres and assuming that the smaller galaxy (of mass $m_g$) is on a circular orbit of radius $r$, its tidal radius and mass become $$r_g^3 \approx  {Gm_gr^2\over 4\sigma^2}, \ \ \ \
m_g \approx {\sigma_g^3r\over 2G\sigma}$$ with $\sigma$ and $\sigma_g$ the velocity dispersions of the larger and smaller galaxies respectively. Chandrasekhar's formula then gives for the orbital decay rate and infall time $${dr\over dt} = -0.30{Gm_g\over\sigma r} \ln\Lambda \approx
-0.151{\sigma^3_g\over\sigma^2}\ln\Lambda, \ \ \ \
\Delta t_2 \approx 0.30 {r_e\sigma^2\over\sigma^3_g}$$ (Tremaine 1990), where $\ln\Lambda\approx \ln(r/0.5r_g)\approx \ln(4\sigma/\sigma_g)\approx
2$. Taking $r(0) = r_e\approx 2.6\ {\rm kpc} (\sigma/200
\ {\rm km\ s}^{-1})^3$ (Valluri & Merritt 1998), the effective radius of the larger galaxy, the inspiral time becomes $$\Delta t_2 \approx 9.6\times 10^7 {\rm yr} \left({\sigma\over 200
\ {\rm km s}^{-1}}\right)^5\left({m_2\over
10^8M_{\odot}}\right)^{-3/4}$$ where eq. (15) has been used to relate $\sigma_g$ to $m_2$. This relation is plotted in Fig. 2b and the total time $t_{tot} = \Delta t_1 + \Delta t_2$ in Fig. 2c. Remarkably, Fig. 2c suggests that real BHs have roughly the mass that would be required to minimize their total coalescence time. Nevertheless, this time appears to be very long.

Other mechanisms might accelerate the coalescence. If the potential of the stellar nucleus is non-axisymmetric, orbits will not conserve angular momentum and stars with a much wider range of energies may interact with the binary. However it is not clear whether self-consistent triaxiality can be maintained in a nucleus where the potential is dominated by a BH (Sridhar & Touma 1999; Merritt & Valluri 1999). Gas if present would also exert a drag force on the binary. "Loss-cone refilling" by two-body interactions between stars (Frank & Rees 1976) is almost certainly unimportant due to the long relaxation times in nuclei (Valtonen 1996).

If coalescence times are comparable to the mean time between galaxy mergers, a third BH would often be introduced into a nucleus containing an uncoalesced pair. The most likely outcome is ejection of two BHs in opposite directions with the third remaining at the center (Valtonen 1996). This process could substantially reduce the mean ratio of BH mass to galaxy mass over time.

# Formation and Destruction of Nuclei

Combining eqs. (4) and (12), the total mass ejected by the BH binary between $t_h$ and $t_{gr}$ is $\sim 4 M_{12}$. The structure of the pre-existing nucleus should therefore be disturbed out to a radius where the enclosed stellar mass was a few times $M_{12}$. A reasonable guess for the initial density profile is $\rho\approx
r^{-\gamma}, \gamma\approx 2$, as assumed in the model above; adiabatic growth of a BH in a nucleus with a shallower profile generically leads to a power-law cusp with index near 2 (Quinlan et al. 1995; Merritt & Quinlan 1998), and this is also the slope observed in the faintest ellipticals (Kormendy et al. 1995) which are the least likely to have experienced the sort of core disruption discussed above.

Simply removing all of the stars with energies below some minimum $E_{min}$ produces a density profile near the BH of $$\rho(r) = 4\pi\int_{E_{min}}^0
f(E)\sqrt{2\left[E-\Phi(r)\right]}dE \approx
C\times f(E_{min})\sqrt{{GM_{12}\over r}} \propto r^{-1/2}$$ assuming an isotropic velocity distribution $f(E)$. Nakano & Makino (1999a,b) found this to be a good description of what happens in $N$-body simulations of a single BH spiralling into a pre-existing core. In the case of a decaying BH binary, stars on eccentric orbits are most likely to interact with the binary and be removed, leaving behind a nucleus in which the stellar motions are strongly biased toward circular. The resulting density profile is difficult to calculate. The $N$-body simulations discussed above (Makino 1997; QH) yield central profiles that are crudely describable as power laws, with indices $0\mathrel{\mathpalette\fun <}\gamma\mathrel{\mathpalette\fun <}1$; however these results are probably dependent on the degree of BH wandering and hence on $N$.

Which galaxies in the current universe would be expected to contain these low-density cores? Kauffmann & Haehnelt (1999) present a semi-analytic model for the formation of galaxies in a cold-dark-matter universe; they assume that supermassive BHs form from gas that is driven into the centers of galaxies during mergers. The predicted ratio of gas mass to stellar mass during the last major merger is a steep function of galaxy luminosity in their model; this ratio drops from $\sim 3$ for $M_v=-18$ to $\sim 0.3$ for $M_v=-21$, albeit with considerable scatter (their Fig. 14). Gas-rich mergers would be expected to form dense nuclei while gas-free mergers should produce low-density cores after the BHs coalesce. Figure 3 shows that the variation of $\gamma$ with galaxy luminosity is roughly consistent with Kauffmann & Haehnelt's model. Faint ellipticals, $M_v\mathrel{\mathpalette\fun >}-19$, have $\gamma\approx 2$, while for brighter galaxies $\gamma$ decreases with increasing luminosity. This trend is often described as a dichotomy although the data of Fig. 3 suggest a continuous distribution, as expected in a model like Kauffmann & Haehnelt's.

# Large-Scale Evolution

The gravitational influence of a supermassive BH can extend far beyond the nucleus in a galaxy that is not axisymmetric (Gerhard & Binney 1985). The mechanism is dynamical chaos induced in the stellar orbits by close passages to the BH (Merritt & Valluri 1996). Although relatively little $N$-body work has been done on this problem, one study (Merritt & Quinlan 1998) found that nuclear point masses can cause initially triaxial galaxies to evolve to globally axisymmetric shapes in little more than a crossing time when the ratio of BH mass to galaxy mass exceeds $\sim 2\%$. This mass ratio is consistent with the value that induces a transition to global stochasticity and rapid mixing in the phase space of box-like orbits (Valluri & Merritt 1998).

Mergers of stellar disks produce generically triaxial objects (Barnes 1996); adding a dissipative component to the simulations produces end states that are much more nearly axisymmetric (Barnes & Hernquist 1996). The evolution in shape occurs rapidly once a few percent of the total mass has accumulated in the center (Barnes 1999). The stars in these simulations respond to the "gas" only insofar as the latter affects the potential; thus the change in shape is a purely stellar-dynamical phenomenon, and we would expect to see similar evolution even in gas-free mergers if the BHs are sufficiently massive.

The typical ratio of BH mass to galaxy mass is $\sim
0.003$ (Ho 1998), somewhat less than the value required to induce a rapid transition to axisymmetry. For $M_{\bullet}/M_{gal}=0.003$, Merritt & Quinlan (1998) found a time scale of $\sim 10^2$ crossing times for the evolution to axisymmetry. This exceeds a Hubble time for galaxies with luminosities above $M_v\approx -19$ (Valluri & Merritt 1998); hence we might expect bright ellipticals to often retain their merger-induced triaxial shapes. In fact there is evidence for a systematic change in the shape of ellipticals at about this magnitude (Tremblay & Merritt 1996); bright ellipticals as a class show evidence for triaxiality, while the axis-ratio distribution of faint ellipticals is consistent with axisymmetry.

The same orbital evolution that destroys triaxiality also tends to produce a smoother phase-space density. One consequence is that "boxiness" in the isodensity contours – a natural consequence of a non-smooth phase space distribution (Binney & Petrou 1985) – should be destroyed along with triaxiality. This too has been observed in simulations of gaseous mergers; the effect is sometimes attributed to "dissipation" (e.g. Bekki & Shioya 1997) but it is again almost certainly a purely stellar dynamical effect. A prediction is that triaxiality should correlate with boxiness, since a central mass concentration or BH should tend to destroy both following a merger. Evidence for such a correlation has been noted by Kormendy & Bender (1996).

I thank M. Haehnelt for useful discussions.


abstract: One of the tests of nucleosynthesis theory is the $^3$He abundance in the Galaxy. $^3$He$^+$ is observed through its 3.46 cm hyperfine level in HII regions and the $^3$He/H ratio compares well with theory. Since $^3$He  can be created or destroyed in nuclear reactions, one would expect that its abundance shows a trend with the amount of such reactions, so with distance to the Center of the Galaxy and with metallicity. Such trends are lacking in observations. This is explained by assuming that the HII clouds are recently formed out of the primordial micro brown dwarfs of earth mass predicted by gravitational hydrodynamics. If indeed existing, they would preserve their primordial $^3$He/H ratio and spread this when evaporating into HII clouds, independent of the location in the Galaxy.
  .
  In the development of the argument, it is also explained that wide binaries do not rule out the MACHO dark matter predicted by gravitational hydrodynamics, but are rather immersed as visible partners in Jeans clusters of dark micro brown dwarfs.
author: Theo M. Nieuwenhuizen$^{1,2}$[^1]
$^{1}$ Center for Cosmology and Particle Physics, New York University, 4 Washington Place, New York, NY 10003, USA
$^{2}$Institute for Theoretical Physics, University of Amsterdam, Science Park 904, P.O. Box 94485, 1090 GL Amsterdam, The Netherlands
date: Accepted ... Received ...; in original form ...
title: Explanation of the Helium-3 problem

# Introduction

Nucleosynthesis in the WMAP era is said to be a parameter-free theory, see e. g. . With cosmology described by the Friedman-Lemaitre-Robertson-Walker metric and the number of neutrino species known, accurate predictions can be made, ready for comparison with observations. We will discuss one of these tests, the one for $^3$He, that is in agreement with theory, but has a puzzling aspect.

# The helium-3 problem

$^3$He is produced in nucleosynthesis and the observed abundance \[$^3$He/H\]$\approx1.0\cdot 10^{-5}$ agrees with theory. But not all is well, so let us explain the way to observe it. The $^3$He nucleus consists of two protons and a neutron, so it has spin $\frac{3}{2}$ or ${\frac{1}{2}}$. The singly ionized $^3$He$^+$ contains one electron, thus the singlet spin-${\frac{1}{2}}$ state brings an analogy with the neutral H-atom. Indeed, there is an analog of the 21-cm spin-flip transition, which occurs roughly at 21/$Z^2$=21/4 cm, more precisely at 3.46 cm. The emission by this transition line provides an observational signature of singly ionized $^3$He located in regions of ionized H gas (so-called HII regions) and in planetary nebulae that also contain ionized H and $^3$He.

Hii  regions have been created rather recently in the cosmic history (), and are thus expected to reflect the chemical history of the surroundings. $^3$He can be created or destroyed in nuclear reactions. One would thus expect its local concentration in the Galaxy to depend on the amount of reactions that have taken place, a measure for which is the local metallicity. Since near the Center more reactions took place, there is a clear gradient in the metallicity, so one also expects a significantly higher $^3$He concentration near the Center than in the outskirts. But though observations show fluctuating concentrations, there is no statistical trend with distance or with metallicity. Specifically, low values of \[$^3$He/H\]$\approx1.0\cdot 10^{-5}$ are observed at 4, 9, 10, 12 and 16 kpc from the Center, while up to two-three times higher values are observed in between ().

The important point for the present paper is that the HII regions show no evidence for stellar $^3$He enrichment during the last 4.5 Gyr; the low values thus being found at any distance from the Center is called *the* $^3$He *problem* (Bania et al. 2007). This situation creates a paradox, since it calls for a delicate balance between creation and destruction *independent of the amount of reactions* (the metallicity). To be compatible with this result, Galactic chemical evolution models () require that $\sim90\%$ of solar analog stars are non-producers of $^3$He. But the two best observed planetary nebulae do have higher concentration (), and consequently should then both be members of the 10% class of stars that do produce $^3$He – which seems unlikely. For a review of the current status of $^3$He evolution, see .

# Micro brown dwarfs from Gravitational hydrodynamics

Gravitational hydrodynamics is an approach that stresses the role of turbulence in cases of small viscosity (more precisely, large Reynolds number) and the possibility of viscous structure formation of turbulent flows on scales where the Reynolds number is low. It is well understood that at the transition from plasma to neutral gas (decoupling or recombination), the newly formed gas breaks up at the Jeans scale, forming Jeans clumps of some 600,000 $M_\odot$. It was put forward by that the Jeans clumps themselves fragment at the viscous scale into objects of earth mass, turning the Jeans gas clumps into Jeans clusters of some $2\cdot10^{11}$ micro brown dwarfs ($\mu$BDs) of earth mass. The Galaxy is predicted to contain about two million Jeans clusters in its halo that make up the full Galactic dark matter (, ).

These objects are belong to the class called MACHOs (Massive Astrophysical Compact Halo Objects). They have been detected in quasar micro lensing () and the observed signature has approximately equal positive and negative events (), uniquely indicating microlensing by a population at unit optical depth. However, in direct searches in front of the Magellanic clouds they have not been detected. Inspection of the literature in this field reveals that this low mass scale was covered only in one paper (). This reports observations carried out in the early and mid nineties using a telescope of 40 cm in diameter. Later MACHO searches did not cover this mass range, so the controversy between the Schild detection and the Renault et al. non-detection was never resolved. It is planned to redo the MACHO search in front of the Magellanic clouds using a much larger telescope (Schild 2011). They will be searched in what is normally called cirrus dust clouds, but what could just be an agglomeration of dark Jeans clusters. One further indication for this bold assumption is that the temperature of "cirrus dust" is about 15 K (), i. .e., near the H triple point of 13.8 K. While cold dust theories have no explanation why the dust should minimally have this temperature, the required release of latent heat of these compact hydrogen clouds would then keep them long at this temperature, before finally freezing.

Analysis of wide binaries in the Galaxy halo has ruled out any MACHO dark matter component more heavy than 43 $M_\odot$ (Yoo et al. 2004). That this conclusion is somewhat too strong ( ), will be of little help in our situation. For a typical wide separation of the components, the major axis is of the order $a=$ 18,000 AU $=0.09$ pc. In our case the relevant MACHO  objects for this application are the individual Jeans clusters of 600,000 $M_\odot$ and radius of 1.4 pc, so they would completely disrupt the binaries – if they were on their own. However, we now point out that these halo wide binaries - being much smaller than the Jeans clusters - must lie inside Jeans clusters themselves, as a result of which they must be much more stable than commonly expected. So we predict that these wide binaries have a center of mass motion of ca 200 km/s, a mutual speed of ca 20 km/s and are embedded in a baryonic dark matter matrix of ca 600,000 $M_\odot$. These features can be tested. Clearly, stabilized inside Jeans clusters, wide binaries do not rule out the MACHO dark matter.

One would expect that the fragmentation of the Jeans cloud has been seen in numerics. When Truelove et al. (1997) noticed instabilities in their simulations, they dismissed them as unphysical and built in Jeans filters to remove them. Recently instabilities at the $10^{-3}-10^{-4}M_\odot$ level were observed in the simulation of the first stars in turbulent gas clouds (Clark et al. 2011): the gas first fragments and then the fragments aggregate to form the stars. The authors carefully explain how the gas can cool, namely grace to the formation of molecular hydrogen (H$_2$) with help of a small fraction of free electrons still present, so that heat can be radiated away through H$_2$ levels. This allows the gas to cool locally after which it can fragment. Though the authors did not carefully sort out what would be the minimal fragmentation scale, their fragmentation at the $10^{-3}-10^{-4}M_\odot$ level is a already strong support for the gravitational hydrodynamics picture and indeed points towards its predicted earth mass MACHOs.

Nuclear synthesis attributes enough matter for the $\mu$BDs. About 4.5% of the critical density of the Universe is baryonic. At best 0.5 % is luminous, and some 2–3% is observed in X-ray gas. The *missing baryon problem* refers to the fact that most of the baryons are unaccounted for, about 60% is missing at cosmic scales; an inventory of baryons in and around the Milky Way reveals at best 25% of the expected baryons. Though they are believed to be located in unobserved relatively cool X-ray clouds, this need not be the full or the only explanation. Indeed, the missing dark baryons may be locked up in $\mu$BDs. Radiation from what is commonly called "cold cirrus clouds" or "cold cirrus dust" is observed to have a temperature of 15 K (Veneziani et al. 2010), or, more generally, temperatures between 40 K and 15 K (). Dust models have no explanation for the minimum of 15 K. But the radiation may actually arise from the thermal atmospheres of $\mu$BDs (). The fact that H has a critical point at 33 K and a triple point at 13.8 K coincides with many observations of "cold dust" temperatures by the Herschel telescope between 15 and 40 K, with the lower values condensing near 15 K, see Figures 2 and 3 of .

The picture of baryonic dark matter locked up in $\mu$BDs grouped in – mostly – dark Jeans clusters (JCs) explains another, not often mentioned problem: *the missing Jeans clusters problem*. Indeed, it is agreed that after the decoupling all gas fragments in Jeans clumps, but where are they now? Gravitational hydrodynamics asserts that they just constitute the halo dark matter of galaxies.

From the point of view of galactic structures there is a lot of support for the picture of Jeans clouds consisting of micro brown dwarfs. Galactic rotation curves flatten when the ca $2\cdot 10^6$ JCs of the Galaxy have an isothermal distribution $\rho(r)\approx v^2/2\pi G r^2$, where $v$ is the velocity dispersion, about 200 km/s for Jeans clusters in the Galaxy and 20 km/s for $\mu$BDs inside a Jeans cluster. The Tully-Fisher and and Faber-Jackson relations follow if one assumes that star formation arises when JCs heat each others mBDs  by tidal forces when they come within a certain radius of each other, which brings the proportionality *Luminosity* $\sim \#\, {\rm stars}\, \sim \int{\rm d}^3r\,\rho^2(r)\sim v^4/R_\ast$ for some characteristic radius $R_\ast$  ().

In galaxy merging the observed young globular clusters may not or not only appear due to tidal disruption, and its unknown star formation process, but also due to tidal heating of the mBDs, which expand and can form stars millions of years after the merging process has taken place (). Thus the galaxy merging process can transform, along the merging path, dark Jeans clusters in situ in the observed bright young globular clusters.

Mysterious radio events were reported by Ofek et al. (2010). They are frequent (1000/square degree/year), radio loud ($>$ 1Jy), and have neither a precursor nor a follow up, and have no detectable counterparts in the infrared, visible or X-ray spectrum. Within isothermal modeling, they have been connected to merging of $\mu$BDs inside Jeans clusters, and the event duration of more than half an hour to several days allowed to estimate their radius as 3 $\times$ (the duration in days) $\times$ (the solar radius) ().

The theory of star formation that has been inconclusive for long, strongly benefits from the principle that star (and planet) formation arises from aggregation of $\mu$BDs (Gibson and Schild 2011). The *iron planet core problem* acknowledges that iron planet cores, like in the earth, are difficult to explain, since iron is mostly observed in oxides. But $\mu$BDs with their size larger than the Sun and weighing only as much as the Earth, can collect the intergalactic iron dust as "vacuum cleaners", while the H atmosphere is all too eager to dissolve the oxides for making water, after which the iron can sink to the center. Iron cores are then explained from glueing small iron cores in the aggregation of the $\mu$BDs to form planets ().

A relation between globular clusters and black holes was discussed () and will be expanded elsewhere, in particular in connection with role of Jeans clusters in solving the so-called last parsec problem in black hole merging processes.

Nowadays galaxies are observed at redshifts basically up to $z=10$. This being considered before the reionization era, they should not be visible but immersed in H clouds, so that many unobserved weak galaxies are invoked, which should ionize the gas and create a "pencil of visibility" in our direction. The galaxy UDFy-38135539 was established to have redshift $z=8.55$, so we see light that it emitted 600 million years after the big bang (). The authors say in the abstract that a significant contribution from other, probably fainter galaxies nearby, is needed to explain its visibility. This is a deus ex machina. $\mu$BDs, on the other hand, offer a more obvious explanation: the hydrogen is locked up in condensed objects, so most of the space is empty and transparent ().

One observes that the amount of dust and the rate of star formation had a peak at redshift $z\sim2$ and are smaller in recent times. This higher amount of dust content in massive galaxies at higher redshift is difficult to explain in standard dust evolution models (). But it finds an easy explanation if one accepts that the "dust" radiation is produced by the atmospheres of $\mu$BDs. In the course of time, they are more and more used up because they coagulate to form heavier objects and stars. Consequently, less fuel for star formation is ultimately available and less $\mu$BD outer surface exists to emit "dust" radiation.

# The answer to the Helium-3 problem

Within gravitational hydrodynamics it is natural to identify the HII clouds with partly or fully evaporated Jeans clusters, in which a large fraction of the $\mu$BDs, if not all, evaporated by heating in a strong star forming region. In such a case the $\mu$BDs of a Jeans cluster can expand into a gas cloud. Their would have kept the primordial $^3$He content up to then, and it would not change when they evaporated into HII clouds, so the primordial value would appear independent of their location in the Galaxy and independent of the local metallicity, which explains the $^3$He problem.

The most massive and largest HII region in the Local Group is the region 30 Doradus in the Large Magellanic Cloud. It has been carefully studied recently (Lopez et al 2008). In addition to a large ionized gas mass $\sim 8 \cdot 10^5M_\odot$ (Kennicutt 1984), the 30 Doradus nebula also has $\sim 10^6M_\odot$ of CO (Johansson et al. 1998). This HII mass is approximately the mass of a single Jeans cluster, supporting the just mentioned gravitational hydrodynamics picture, while the CO content exhibits a large metallicity.

In conclusion, in gravitational hydrodynamics many problems, like the 15 K cold dust temperature, the visibility of early galaxies and the $^3$He problem, find a simple explanation.

## Acknowledgment

It is a pleasure to thank David Hogg for pointing us at the wide binary problem.

[^1]: E-mail: email@example.com