Dynamics and Merger Rate of Primordial Black Holes in a Cluster

The PBH clusters can be sources of gravitational waves, and the merger rate depends on the spatial distribution of PBHs in the cluster which changes over time. It is well known that gravitational collisional systems experience the core collapse that leads to significant increase of the central density and shrinking of the core. After core collapse, the cluster expands almost self-similarly (i.e., density profile extends in size without changing its shape). These dynamic processes affect the merger rate of PBHs. In this paper, the dynamics of the PBH cluster is considered using the Fokker–Planck equation. We calculate the merger rate of PBHs on cosmic time scales and show that its time dependence has a unique signature. Namely, it grows by about an order of magnitude at the moment of core collapse which depends on the characteristics of the cluster, and then decreases according to the dependence R∝t−1.48. It was obtained for monochromatic and power-law PBH mass distributions with some fixed parameters. Obtained results can be used to test the model of the PBH clusters via observation of gravitational waves at high redshift.

Shedding Light On An Undiscovered Aspect That Perfectly Mimics Cosmic Acceleration

The comparison of redshift-distance relationship for high and low-redshift supernovae revealed the surprising transition of the Universe’s expansion from deceleration to acceleration. As compared to local supernovae, remote supernovae appear 10% to 25% dimmer as they are further away than expected. The expansion rate obtained for local supernovae is higher with low redshifts as compared to the expansion rate obtained for remote supernovae with high redshifts. Since observed redshifts in an expanding Universe provide an estimate of recession velocities, therefore, it is very disturbing to find that low recession velocities (just 1% of speed of light) indicate a faster rate of expansion (acceleration), whereas high recession velocities (60% of speed of light) indicate a slower rate of expansion (deceleration). In this paper, I unravel an undiscovered aspect that perfectly mimics cosmic acceleration. Rather than “cosmic deceleration that preceded the current epoch of cosmic acceleration”, I show in this paper, that “consecutive expansion epochs of the Universe that preceded the current epoch of cosmic expansion” were responsible for placing remote supernovae further away than expected. As a consequence of consecutive expansion, expansion began for remote structures in preceding expansion epochs before it did for local structures in the current (or more recent) expansion epoch; remote supernovae, quasars, and gamma-ray bursts are therefore not only further away than expected, but they also happen to yield a slower rate of expansion, thereby suggesting their deceleration even with “superluminal expansion”. As a result of consecutive expansion, preceding expansion epochs appear to be decelerating as compared to the expansion epoch that succeeds them. The analysis is based on the redshift-distance relationship plotted for 580 type Ia supernovae from the Supernova Cosmology Project, 7 additional high-redshift type Ia supernovae discovered through the Advanced Camera for Surveys on the Hubble Space Telescope from the Great Observatories Origins Deep Survey Treasury program, and 1 additional very high-redshift type Ia supernova discovered with Wide Field and Planetary Camera 2 on the Hubble Space Telescope. The results obtained by the High-Z Supernova Search Team through observations of type Ia supernovae have also been analysed. Studies incorporating quasars and gamma-ray bursts to determine how the expansion of the Universe has changed over time have been taken into consideration as well. The results obtained in this paper have been confirmed by plotting velocity-distance relationship, expansion rate vs. time relationship, expansion factor vs. time relationship, scale factor vs. time relationship, scale factor vs. distance relationship, distance-redshift relationship, and distance modulus vs. redshift relationship, moreover, the deceleration parameter (q0) is also found to be negative (q0 < 0).

The Black Hole Mass Function Across Cosmic Times. I. Stellar Black Holes and Light Seed Distribution

This is the first paper in a series aimed at modeling the black hole (BH) mass function, from the stellar to the intermediate to the (super)massive regime. In the present work, we focus on stellar BHs and provide an ab initio computation of their mass function across cosmic times; we mainly consider the standard, and likely dominant, production channel of stellar-mass BHs constituted by isolated single/binary star evolution. Specifically, we exploit the state-of-the-art stellar and binary evolutionary code SEVN, and couple its outputs with redshift-dependent galaxy statistics and empirical scaling relations involving galaxy metallicity, star formation rate and stellar mass. The resulting relic mass function dN / dVd log m • as a function of the BH mass m • features a rather flat shape up to m • ≈ 50 M ⊙ and then a log-normal decline for larger masses, while its overall normalization at a given mass increases with decreasing redshift. We highlight the contribution to the local mass function from isolated stars evolving into BHs and from binary stellar systems ending up in single or binary BHs. We also include the distortion on the mass function induced by binary BH mergers, finding that it has a minor effect at the high-mass end. We estimate a local stellar BH relic mass density of ρ • ≈ 5 × 107 M ⊙ Mpc−3, which exceeds by more than two orders of magnitude that in supermassive BHs; this translates into an energy density parameter Ω• ≈ 4 × 10−4, implying that the total mass in stellar BHs amounts to ≲1% of the local baryonic matter. We show how our mass function for merging BH binaries compares with the recent estimates from gravitational wave observations by LIGO/Virgo, and discuss the possible implications for dynamical formation of BH binaries in dense environments like star clusters. We address the impact of adopting different binary stellar evolution codes (SEVN and COSMIC) on the mass function, and find the main differences to occur at the high-mass end, in connection with the numerical treatment of stellar binary evolution effects. We highlight that our results can provide a firm theoretical basis for a physically motivated light seed distribution at high redshift, to be implemented in semi-analytic and numerical models of BH formation and evolution. Finally, we stress that the present work can constitute a starting point to investigate the origin of heavy seeds and the growth of (super)massive BHs in high-redshift star-forming galaxies, that we will pursue in forthcoming papers.

Driving Galactic Outflows with Magnetic Fields at Low and High Redshift

Although galactic outflows play a key role in our understanding of the evolution of galaxies, the exact mechanism by which galactic outflows are driven is still far from being understood and, therefore, our understanding of associated feedback mechanisms that control the evolution of galaxies is still plagued by many enigmas. In this work, we present a simple toy model that can provide insight on how non-axisymmetric instabilities in galaxies (bars, spiral arms, warps) can lead to local exponential magnetic field growth by radial flows beyond the equipartition value by at least two orders of magnitude on a timescale of a few 100 Myr. Our predictions show that the process can lead to galactic outflows in barred spiral galaxies with a mass-loading factor η ≈ 0.1, in agreement with our numerical simulations. Moreover, our outflow mechanism could contribute to an understanding of the large fraction of barred spiral galaxies that show signs of galactic outflows in the chang-es survey. Extending our model shows the importance of such processes in high-redshift galaxies by assuming equipartition between magnetic energy and turbulent energy. Simple estimates for the star formation rate in our model together with cross correlated masses from the star-forming main sequence at redshifts z ∼ 2 allow us to estimate the outflow rate and mass-loading factors by non-axisymmetric instabilities and a subsequent radial inflow dynamo, giving mass-loading factors of η ≈ 0.1 for galaxies in the range of M ⋆ = 109–1012 M ⊙, in good agreement with recent results of sinfoni and kmos 3D.

COLDz: Probing Cosmic Star Formation With Radio Free–Free Emission

Radio free–free emission is considered to be one of the most reliable tracers of star formation in galaxies. However, as it constitutes the faintest part of the radio spectrum—being roughly an order of magnitude less luminous than radio synchrotron emission at the GHz frequencies typically targeted in radio surveys—the usage of free–free emission as a star formation rate tracer has mostly remained limited to the local universe. Here, we perform a multifrequency radio stacking analysis using deep Karl G. Jansky Very Large Array observations at 1.4, 3, 5, 10, and 34 GHz in the COSMOS and GOODS-North fields to probe free–free emission in typical galaxies at the peak of cosmic star formation. We find that z ∼ 0.5–3 star-forming galaxies exhibit radio emission at rest-frame frequencies of ∼65–90 GHz that is ∼1.5–2 times fainter than would be expected from a simple combination of free–free and synchrotron emission, as in the prototypical starburst galaxy M82. We interpret this as a deficit in high-frequency synchrotron emission, while the level of free–free emission is as expected from M82. We additionally provide the first constraints on the cosmic star formation history using free–free emission at 0.5 ≲ z ≲ 3, which are in good agreement with more established tracers at high redshift. In the future, deep multifrequency radio surveys will be crucial in order to accurately determine the shape of the radio spectrum of faint star-forming galaxies, and to further establish radio free–free emission as a tracer of high-redshift star formation.

Cosmology at high redshift — a probe of fundamental physics

An observational program focused on the high redshift (2<z<6) Universe has the opportunity to dramatically improve over upcoming LSS and CMB surveys on measurements of both the standard cosmological model and its extensions. Using a Fisher matrix formalism that builds upon recent advances in Lagrangian perturbation theory, we forecast constraints for future spectroscopic and 21-cm surveys on the standard cosmological model, curvature, neutrino mass, relativistic species, primordial features, primordial non-Gaussianity, dynamical dark energy, and gravitational slip. We compare these constraints with those achievable by current or near-future surveys such as DESI and Euclid, all under the same forecasting formalism, and compare our formalism with traditional linear methods. Our Python code FishLSS — used to calculate the Fisher information of the full shape power spectrum, CMB lensing, the cross-correlation of CMB lensing with galaxies, and combinations thereof — is publicly available.

Mapping the Morphology and Kinematics of a Lyα-selected Nebula at z = 3.15 with MUSE

Recent wide-field integral-field spectroscopy has revealed the detailed properties of high-redshift Lyα nebulae, most often targeted due to the presence of an active galactic nucleus (AGN). Here, we use VLT/MUSE to resolve the morphology and kinematics of a nebula initially identified due to strong Lyα emission at z ∼ 3.2 (LABn06). Our observations reveal a two-lobed Lyα nebula, at least ∼173 pkpc in diameter, with a light-weighted centroid near a mid-infrared source (within ≈17.2 pkpc) that appears to host an obscured AGN. The Lyα emission near the AGN is also coincident in velocity with the kinematic center of the nebula, suggesting that the nebula is both morphologically and kinematically centered on the AGN. Compared to AGN-selected Lyα nebulae, the surface-brightness profile of this nebula follows a typical exponential profile at large radii (>25 pkpc), although at small radii, the profile shows an unusual dip at the location of the AGN. The kinematics and asymmetry are similar to, and the C iv and He ii upper limits are consistent with, other AGN-powered Lyα nebulae. Double-peaked and asymmetric line profiles suggest that Lyα resonant scattering may be important in this nebula. These results support the picture of the AGN being responsible for powering a Lyα nebula that is oriented roughly in the plane of the sky. Further observations will explore whether the central surface-brightness depression is indicative of either an unusual gas or dust distribution or variation in the ionizing output of the AGN over time.

Finding High-redshift Galaxies with JWST

One of the primary goals for the upcoming James Webb Space Telescope is to observe the first galaxies. Predictions for planned and proposed surveys have typically focused on average galaxy counts, assuming a random distribution of galaxies across the observed field. The first and most-massive galaxies, however, are expected to be tightly clustered, an effect known as cosmic variance. We show that cosmic variance is likely to be the dominant contribution to uncertainty for high-redshift mass and luminosity functions, and that median high-redshift and high-mass galaxy counts for planned observations lie significantly below average counts. Several different strategies are considered for improving our understanding of the first galaxies, including adding depth, area, and independent pointings. Adding independent pointings is shown to be the most efficient both for discovering the single highest-redshift galaxy and also for constraining mass and luminosity functions.

Multiphase Outflows in High-redshift Quasar Host Galaxies

We present Atacama Large Millimeter/submillimeter Array (ALMA) observations of six radio-loud quasar host galaxies at z = 1.4–2.3. We combine the kiloparsec-scale resolution ALMA observations with high spatial resolution adaptive optics integral field spectrograph data of the ionized gas. We detect molecular gas emission in five quasar host galaxies and resolve the molecular interstellar medium using the CO (3–2) or CO (4–3) rotational transitions. Clumpy molecular outflows are detected in four quasar host galaxies and a merger system 21 kpc away from one quasar. Between the ionized and cold molecular gas phases, the majority of the outflowing mass is in a molecular phase, while for three out of four detected multiphase gas outflows, the majority of the kinetic luminosity and momentum flux is in the ionized phase. Combining the energetics of the multiphase outflows, we find that their driving mechanism is consistent with energy-conserving shocks produced by the impact of the quasar jets with the gas in the galaxy. By assessing the molecular gas mass to the dynamics of the outflows, we estimate a molecular gas depletion timescale of a few megayears. The gas outflow rates exceed the star formation rates, suggesting that quasar feedback is a major mechanism of gas depletion at the present time. The coupling efficiency between the kinetic luminosity of the outflows and the bolometric luminosity of the quasar of 0.1%–1% is consistent with theoretical predictions. Studying multiphase gas outflows at high redshift is important for quantifying the impact of negative feedback in shaping the evolution of massive galaxies.

Thermal Instabilities and Shattering in the High-redshift WHIM: Convergence Criteria and Implications for Low-metallicity Strong H i Absorbers

Using a novel suite of cosmological simulations zooming in on a megaparsec-scale intergalactic sheet (pancake) at z ∼ (3–5), we conduct an in-depth study of the thermal properties and H i content of the warm-hot intergalactic medium (WHIM) at those redshifts. The simulations span nearly three orders of magnitude in gas cell mass, ∼(7.7 × 106–1.5 × 104)M ⊙, one of the highest-resolution simulations of such a large patch of the intergalactic medium (IGM) to date. At z ∼ 5, a strong accretion shock develops around the pancake. Gas in the postshock region proceeds to cool rapidly, triggering thermal instabilities and generating a multiphase medium. We find the mass, morphology, and distribution of H i in the WHIM to all be unconverged, even at our highest resolution. Interestingly, the lack of convergence is more severe for the less-dense, metal-poor intrapancake medium (IPM) in between filaments and far outside galaxies. With increased resolution, the IPM develops a shattered structure with most of the H i in kiloparsec-scale clouds. From our lowest-to-highest resolution, the covering fraction of metal-poor (Z < 10−3 Z ⊙) Lyman-limit systems (N H I > 1017.2cm−2) in the z ∼ 4 IPM increases from ∼(3–15)%, while that of metal-poor damped Lyα absorbers (N H I > 1020cm−2) increases from ∼(0.2–0.6)%, with no sign of convergence. We find that a necessary condition for the formation of a multiphase shattered structure is resolving the cooling length, l cool = c s t cool, at T ∼ 105 K. If this is unresolved, gas “piles up” at T ≲ 105 K and further cooling becomes very inefficient. We conclude that state-of-the-art cosmological simulations are still unable to resolve the multiphase structure of the WHIM, with potentially far-reaching implications.