In this paper we study gravitational lensing of gravitational wave events. The main point of the present work is to introduce a semi-analytic approach so that each ingredient can be varied and tested individually. Our analytic model for the source population is motivated by a numerical study and we compare semi-analytical and numerical results. We determine the expected magnification for events seen at a given observed luminosity distance. We find that while the probability of significant magnification of the observed LIGO-Virgo events is very small, the probability distribution of the magnification has a significant tail to high $\mu$ such that e.g. the variance of the magnification is very large and even diverges in the geometric optics approximation. For the 10 binary black hole mergers observed by LIGO-Virgo in the O1+O2 observation campaigns, the probability that one of them has been magnified with magnification of 5 or bigger is $\mathcal{P}_{\text{obs}}(>5)\sim 0.01$ while the probability of magnification by 50 or bigger is $\mathcal{P}_{\text{obs}}(>50)\sim 0.005$.

We revisit the predictions for the merger rate of massive black hole binaries detectable by the Laser Interferometer Space Antenna (LISA) and their background signal for pulsar-timing arrays. We focus on the effect of the delays between the merger of galaxies and the final coalescence of black hole binaries, and on the effect of supernova feedback on the growth of black holes. By utilizing a semi-analytic galaxy formation model, not only do we account for the processes that drive the evolution of binaries at separations $\lesssim 1$ pc (gas-driven migration, stellar hardening and triple/quadruple massive black hole systems), but we also improve on previous studies by accounting for the time spent by massive black hole pairs from kpc down to a few pc separation. We also include the effect of supernova feedback, which may eject a substantial amount of gas from the nuclear region of low-mass galaxies, thus hampering the growth of black holes via accretion and suppressing their orbital migration in circumbinary disks. In spite of the inclusion of these novel physical effects, we predict that the LISA detection rate should still be in excess of $\sim 2 \mbox{yr}^{-1}$, irrespective of the model for the seeds of the black hole population at high redshifts. However, scenarios in which black holes form from $\sim100 M_\odot$ seeds are more significantly impacted by the feedback from supernovae. We also present predictions for the mass ratio distribution of the merger population, and find that binaries typically have mass ratios between $\sim 0.1$ and $1$. Predictions for the stochastic background in the band of pulsar-timing array experiments are instead rather robust, and show only a mild dependence on the model.

We show that the anisotropies of the astrophysical stochastic gravitational wave background in the mHz band have a strong dependence on the modelling of galactic and sub-galactic physics. We explore a wide range of self-consistent astrophysical models for stellar evolution and for the distribution of orbital parameters, all calibrated such that they predict the same number of resolved mergers to fit the number of detections during LIGO/Virgo O1+O2 observations runs. We show that different physical choices for the process of black hole collapse and cut-off in the black hole mass distribution give fractional differences in the angular power spectrum of anisotropies up to 50\% on all angular scales. We also point out that the astrophysical information which can be extracted from anisotropies is complementary to the isotropic background and individual mergers. These results underline the interest in the anisotropies of the stochastic gravitational wave background as a new and potentially rich field of research, at the cross-road between astrophysics and cosmology.

Since the very beginning of astronomy the location of objects on the sky has been a fundamental observational quantity that has been taken for granted. While precise two dimensional positional information is easy to obtain for observations in the electromagnetic spectrum, the positional accuracy of current and near future gravitational wave detectors is limited to between tens and hundreds of square degrees, which makes it extremely challenging to identify the host galaxies of gravitational wave events or to confidently detect any electromagnetic counterparts. Gravitational wave observations provide information on source properties and distances that is complementary to the information in any associated electromagnetic emission and that is very hard to obtain in any other way. Observing systems with multiple messengers thus has scientific potential much greater than the sum of its parts. A gravitational wave detector with higher angular resolution would significantly increase the prospects for finding the hosts of gravitational wave sources and triggering a multi-messenger follow-up campaign. An observatory with arcminute precision or better could be realised within the Voyage 2050 programme by creating a large baseline interferometer array in space and would have transformative scientific potential. Precise positional information of standard sirens would enable precision measurements of cosmological parameters and offer new insights on structure formation; a high angular resolution gravitational wave observatory would allow the detection of a stochastic background and resolution of the anisotropies within it; it would also allow the study of accretion processes around black holes; and it would have tremendous potential for tests of modified gravity and the discovery of physics beyond the Standard Model.

This article explores the properties (amplitude and shape) of the angular power spectrum of the anisotropies of the astrophysical gravitational wave background (AGWB) focusing on the signatures of the astrophysical models describing sub-galactic physics. It demonstrates that while some parameters have negligible impact others, and in particular the stellar evolution models, the metallicity and the merger time delay distribution can result in relative differences of order 40% in the angular power spectrum of anisotropies in both the LIGO/Virgo and LISA frequency bands. It is also shown that the monopole and the anisotropic components of the AGWB are complementary and sensitive to different astrophysical parameters. It follows that AGWB anisotropies are a new observable with the potential to provide new astrophysical information that can not be accessed otherwise.

The LIGO discoveries have rekindled suggestions that primordial black holes (BHs) may constitute part to all of the dark matter (DM) in the Universe. Such suggestions came from 1) the observed merger rate of the BHs, 2) their unusual masses, 3) their low/zero spins, and 4) also from the independently uncovered cosmic infrared background (CIB) fluctuations signal of high amplitude and coherence with unresolved cosmic X-ray background (CXB). Here we summarize the prospects to resolve this important issue with electromagnetic observations using the instruments and tools expected in the 2020's. These prospects appear promising to make significant, and potentially critical, advances. We demonstrate that in the next decade, new space- and ground-borne electromagnetic instruments, combined with concurrent theoretical efforts, should shed critical light on the long-considered link between primordial BHs and DM. Specifically the new data and methodologies under this program will involve: I) Probing with high precision the spatial spectrum of source-subtracted CIB with Euclid and WFIRST, and its coherence with unresolved cosmic X-ray background using eROSITA and Athena, II) Advanced searches for microlensing of Galactic stars by the intervening Galactic Halo BHs with OGLE, Gaia, LSST and WFIRST, III) Supernovae (SNe) lensing in the upcoming surveys with WFIRST, LSST and also potentially with Euclid and JWST, IV) Advanced theoretical work to understand the details of PBH accretion and evolution and their influence on cosmic microwave background (CMB) anisotropies in light of the next generation CMB experiments, V) Better new samples and theoretical understanding involving stability and properties of ultra faint dwarf galaxies, pulsar timing, and cosmological quasar lensing.

We investigate the discrepancy pointed out by Jenkins et al. in Ref. [1] between the predictions of anisotropies of the astrophysical gravitational wave (GW) background, derived using different methods in Cusin et al. [2] and in Jenkins et al. [3]. We show that this discrepancy is not due to our treatment of galaxy clustering, contrary to the claim made in Ref. [1] and we show that our modeling of clustering gives results in very good agreement with observations. Furthermore we show that the power law spectrum used in Refs. [1] and [3] to describe galaxy clustering is incorrect on large scales and leads to a different scaling of the multipoles $C_\ell$. Moreover, we also explain that the analytic derivation of the gravitational wave background correlation function in Refs. [1] and [3] is mathematically ill-defined and predicts an amplitude of the angular power spectrum which depends on the (arbitrary) choice of a non-physical cut-off.

Gravitational-wave astronomy has the potential to explore one of the deepest and most puzzling aspects of Einstein's theory: the existence of black holes. A plethora of ultracompact, horizonless objects have been proposed to arise in models inspired by quantum gravity. These objects may solve Hawking's information-loss paradox and the singularity problem associated with black holes, while mimicking almost all of their classical properties. They are, however, generically unstable on relatively short timescales. Here, we show that this "ergoregion instability" leads to a strong stochastic background of gravitational waves, at a level detectable by current and future gravitational-wave detectors. The absence of such background in the first observation run of Advanced LIGO already imposes the most stringent limits to date on black-hole alternatives, showing that certain models of "quantum-dressed" stellar black holes can be at most a small percentage of the total population. The future LISA mission will allow for similar constraints on supermassive black-hole mimickers.

Recent detections of gravitational waves from merging binary black holes opened new possibilities to study the evolution of massive stars and black hole formation. In particular, stellar evolution models may be constrained on the basis of the differences in the predicted distribution of black hole masses and redshifts. In this work we propose a framework that combines galaxy and stellar evolution models and use it to predict the detection rates of merging binary black holes for various stellar evolution models. We discuss the prospects of constraining the shape of the time delay distribution of merging binaries using just the observed distribution of chirp masses. Finally, we consider a generic model of primordial black hole formation and discuss the possibility of distinguishing it from stellar-origin black holes.

The synthesis of $r$-process elements is known to involve extremely energetic explosions. At the same time, recent observations find significant $r$-process enrichment even in extremely small ultra-faint dwarf (UFD) galaxies. This raises the question of retainment of those elements within their hosts. We estimate the retainment fraction and find that it is large $\sim 0.9$, unless the $r$-process event is very energetic ($\gtrsim 10^{52}$erg) and / or the host has lost a large fraction of its gas prior to the event. We estimate the $r$-process mass per event and rate as implied by abundances in UFDs, taking into account imperfect retainment and different models of UFD evolution. The results are consistent with previous estimates (Beniamini et al. 2016) and with the constraints from the recently detected macronova accompanying a neutron star merger (GW170817). We also estimate the distribution of abundances predicted by these models. We find that $\sim 0.07$ of UFDs should have $r$-process enrichment. The results are consistent with both the mean values and the fluctuations of [Eu/Fe] in galactic metal poor stars, supporting the possibility that UFDs are the main 'building blocks' of the galactic halo population.

We present the first predictions for the angular power spectrum of the astrophysical gravitational wave background constituted of the radiation emitted by all resolved and unresolved astrophysical sources. Its shape and amplitude depend on both the astrophysical properties on galactic scales and on cosmological properties. We show that the angular power spectrum behaves as $C_{\ell}\propto 1/{\ell}$ on large scales and that relative fluctuations of the signal are of order 30% at 100 Hz. We also present the correlations of the astrophysical gravitational wave background with weak-lensing and galaxy distribution. These numerical results pave the way to the study of a new observable at the crossroad between general relativity, astrophysics and cosmology.

The abundance of nitrogen in the interstellar medium is a powerful probe of star for- mation processes over cosmological timescales. Since nitrogen can be produced both in massive and intermediate-mass stars with metallicity-dependent yields, its evolution is challenging to model, as evidenced by the differences between theoretical predictions and observations. In this work we attempt to identify the sources of these discrepancies using a cosmic evolution model. To further complicate matters, there is considerable dispersion in the abundances from observations of DLAs at redshift 2 - 3. We study the evolution of nitrogen with a detailed chemical evolution model and find good agreement with observations, including the relative abundances of N/O and N/Si ratios. We find that the principal contribution of nitrogen comes from intermediate mass stars, with the exception of systems with the lowest N/H, where nitrogen production might possibly be dominated by massive stars. This last result could be strengthened if stellar rotation which is important at low metallicity can produce significant amounts of nitrogen. Moreover, these systems likely reside in host galaxies with stellar masses below 10**8.5 solar mass. We also study the origin of the observed dispersion in nitrogen abundances using the cosmological hydrodynamical simulations Horizon-AGN. We conclude that this dispersion can originate from two effects: differ- ence in the masses of the DLA host galaxies, and difference in the their position inside the galaxy.

Ultralight bosons can induce superradiant instabilities in spinning black holes, tapping their rotational energy to trigger the growth of a bosonic condensate. Possible observational imprints of these boson clouds include (i) direct detection of the nearly monochromatic (resolvable or stochastic) gravitational waves emitted by the condensate, and (ii) statistically significant evidence for the formation of "holes" at large spins in the spin versus mass plane (sometimes also referred to as "Regge plane") of astrophysical black holes. In this work, we focus on the prospects of LISA and LIGO detecting or constraining scalars with mass in the range $m_s\in [10^{-19},\,10^{-15}]$ eV and $m_s\in [10^{-14},\,10^{-11}]$ eV, respectively. Using astrophysical models of black-hole populations calibrated to observations and black-hole perturbation theory calculations of the gravitational emission, we find that, in optimistic scenarios, LIGO could observe a stochastic background of gravitational radiation in the range $m_s\in [2\times 10^{-13}, 10^{-12}]$ eV, and up to $10^4$ resolvable events in a $4$-year search if $m_s\sim 3\times 10^{-13}\,{\rm eV}$. LISA could observe a stochastic background for boson masses in the range $m_s\in [5\times 10^{-19}, 5\times 10^{-16}]$, and up to $\sim 10^3$ resolvable events in a $4$-year search if $m_s\sim 10^{-17}\,{\rm eV}$. LISA could further measure spins for black-hole binaries with component masses in the range $[10^3, 10^7]~M_\odot$, which is not probed by traditional spin-measurement techniques. A statistical analysis of the spin distribution of these binaries could either rule out scalar fields in the mass range $\sim [4 \times 10^{-18}, 10^{-14}]$ eV, or measure $m_s$ with ten percent accuracy if light scalars in the mass range $\sim [10^{-17}, 10^{-13}]$ eV exist.

Ultralight scalar fields around spinning black holes can trigger superradiant instabilities, forming a long-lived bosonic condensate outside the horizon. We use numerical solutions of the perturbed field equations and astrophysical models of massive and stellar-mass black hole populations to compute, for the first time, the stochastic gravitational-wave background from these sources. In optimistic scenarios the background is observable by Advanced LIGO and LISA for field masses $m_s$ in the range $\sim [2\times 10^{-13}, 10^{-12}]\,{\rm eV}$ and $\sim 5\times[ 10^{-19}, 10^{-16}]\,{\rm eV}$, respectively, and it can affect the detectability of resolvable sources. Our estimates suggest that an analysis of the stochastic background limits from LIGO O1 might already be used to marginally exclude axions with mass $\sim 10^{-12.5}{\rm eV}$. Semicoherent searches with Advanced LIGO (LISA) should detect $\sim 15~(5)$ to $200~(40)$ resolvable sources for scalar field masses $3\times 10^{-13}$ ($10^{-17}$) eV. LISA measurements of massive BH spins could either rule out bosons in the range $\sim [10^{-18}, 2\times 10^{-13}]$ eV, or measure $m_s$ with ten percent accuracy in the range $\sim[10^{-17}, 10^{-13}]$ eV.

Massive black-hole binaries, formed when galaxies merge, are among the primary sources of gravitational waves targeted by ongoing Pulsar Timing Array (PTA) experiments and the upcoming space-based LISA interferometer. However, their formation and merger rates are still highly uncertain. Recent upper limits on the stochastic gravitational-wave background obtained by PTAs are starting being in marginal tension with theoretical models for the pairing and orbital evolution of these systems. This tension can be resolved by assuming that these binaries are more eccentric or interact more strongly with the environment (gas and stars) than expected, or by accounting for possible selection biases in the construction of the theoretical models. However, another (pessimistic) possibility is that these binaries do not merge at all, but stall at large ($\sim$ pc) separations. We explore this extreme scenario by using a galaxy-formation semi-analytic model including massive black holes (isolated and in binaries), and show that future generations of PTAs will detect the stochastic gravitational-wave background from the massive black-hole binary population within $10-15$ years of observations, even in the "nightmare scenario" in which all binaries stall at the hardening radius. Moreover, we argue that this scenario is too pessimistic, because our model predicts the existence of a sub-population of binaries with small mass ratios ($q \lesssim 10^{-3}$) that should merge within a Hubble time simply as a result of gravitational-wave emission. This sub-population will be observable with large signal-to-noise ratios by future PTAs thanks to next-generation radio telescopes such as SKA or FAST, and possibly by LISA.

Modeling the stochastic gravitational wave background from various astrophysical sources is a key objective in view of upcoming observations with ground- and space-based gravitational wave observatories such as Advanced LIGO, VIRGO, eLISA and PTA. We develop a synthetic model framework that follows the evolution of single and binary compact objects in an astrophysical context. We describe the formation and merger rates of binaries, the evolution of their orbital parameters with time and the spectrum of emitted gravitational waves at different stages of binary evolution. Our approach is modular and allows us to test and constrain different ingredients of the model, including stellar evolution, black hole formation scenarios and the properties of binary systems. We use this framework in the context of a particularly well-motivated astrophysical setup to calculate the gravitational wave background from several types of sources, including inspiraling stellar-mass binary black holes that have not merged during a Hubble time. We find that this signal, albeit weak, has a characteristic shape that can help constrain the properties of binary black holes in a way complementary to observations of the background from merger events. We discuss possible applications of our framework in the context of other gravitational wave sources, such as supermassive black holes.

The recent detection of the binary black hole merger GW150914 demonstrates the existence of black holes more massive than previously observed in X-ray binaries in our Galaxy. This article explores different scenarios of black hole formation in the context of self-consistent cosmic chemical evolution models that simultaneously match observations of the cosmic star formation rate, optical depth to reionization and metallicity of the interstellar medium. This framework is used to calculate the mass distribution of merging black hole binaries and its evolution with redshift. We also study the implications of the black hole mass distribution for the stochastic gravitational wave background from mergers and from core collapse events.

Deuterium is created during Big Bang Nucleosynthesis, and, in contrast to the other light stable nuclei, can only be destroyed thereafter by fusion in stellar interiors. In this paper we study the cosmic evolution of the deuterium abundance in the interstellar medium and its dispersion using realistic galaxy evolution models. We find that models that reproduce the observed metal abundance are compatible with observations of the deuterium abundance in the local ISM and z ~ 3 absorption line systems. In particular, we reproduce the low astration factor which we attribute to a low global star formation efficiency. We calculate the dispersion in deuterium abundance arising from different structure formation histories in different parts of the Universe. Our model also predicts an extremely tight correlation between deuterium and metal abundances which could be used to measure the primordial deuterium abundance.

Recent chemical abundance measurements of damped Ly-alpha absorbers (DLAs) revealed an intrinsic scatter in their metallicity of ~0.5 dex out to z~5. In order to explore the origin of this scatter, we build a semi-analytic model which traces the chemical evolution of the interstellar matter in small regions of the Universe with different mean density, from over- to underdense regions. We show that the different histories of structure formation in these regions, namely halo abundance, mass and stellar content, is reflected in the chemical properties of the protogalaxies, and in particular of DLAs. We calculate mean metallicity-redshift relations and show that the metallicity dispersion arising from this environmental effect amounts to ~0.25 dex and is an important contributor to the observed overall intrinsic scatter.

The mass fraction of hot gas in clusters is a basic quantity whose level and dependence on the cluster mass and redshift are intimately linked to all cluster X-ray and SZ measures. Modeling the evolution of the gas fraction is clearly a necessary ingredient in the description of the hierarchical growth of clusters through mergers of subclumps and mass accretion on the one hand, and the dispersal of gas from the cluster galaxies by tidal interactions, galactic winds, and ram pressure stripping on the other hand. A reasonably complete description of this evolution can only be given by very detailed hydrodynamical simulations, which are, however, resource-intensive, and difficult to implement in the mapping of parameter space. A much more practical approach is the use of semi-analytic modeling that can be easily implemented to explore a wide range of parameters. We present first results from a simple model that describes the build up of the gas mass fraction in clusters by following the overall impact of the above processes during the merger and accretion history of each cluster in the ensemble. Acceptable ranges for model parameters are deduced through comparison with results of X-ray observations. Basic implications of our work for modeling cluster statistical properties, and the use of these properties in joint cosmological data analyses, are discussed.

Galaxy formation is at the forefront of observation and theory in cosmology. An improved understanding is essential for improving our knowledge both of the cosmological parameters, of the contents of the universe, and of our origins. In these lectures intended for graduate students, galaxy formation theory is reviewed and confronted with recent observational issues. In Lecture 1, the following topics are presented: star formation considerations, including IMF, star formation efficiency and star formation rate, the origin of the galaxy luminosity function, and feedback in dwarf galaxies. In Lecture 2, we describe formation of disks and massive spheroids, including the growth of supermassive black holes, negative feedback in spheroids, the AGN-star formation connection, star formation rates at high redshift and the baryon fraction in galaxies.

Statistical measures of galaxy clusters are sensitive to neutrino masses in the sub-eV range. We explore the possibility of using cluster number counts from the ongoing PLANCK/SZ and future cosmic-variance-limited surveys to constrain neutrino masses from CMB data alone. The precision with which the total neutrino mass can be determined from SZ number counts is limited mostly by uncertainties in the cluster mass function and intracluster gas evolution; these are explicitly accounted for in our analysis. We find that projected results from the PLANCK/SZ survey can be used to determine the total neutrino mass with a (1\sigma) uncertainty of 0.06 eV, assuming it is in the range 0.1-0.3 eV, and the survey detection limit is set at the 5\sigma significance level. Our results constitute a significant improvement on the limits expected from PLANCK/CMB lensing measurements, 0.15 eV. Based on expected results from future cosmic-variance-limited (CVL) SZ survey we predict a 1\sigma uncertainty of 0.04 eV, a level comparable to that expected when CMB lensing extraction is carried out with the same experiment. A few percent uncertainty in the mass function parameters could result in up to a factor \sim 2-3 degradation of our PLANCK and CVL forecasts. Our analysis shows that cluster number counts provide a viable complementary cosmological probe to CMB lensing constraints on the total neutrino mass.

We have recently developed an extended merger-tree model that efficiently follows hierarchical evolution of galaxy clusters and provides a quantitative description of both their dark matter and gas properties. We employed this diagnostic tool to calculate the thermal SZ power spectrum and cluster number counts, accounting explicitly for uncertainties in the relevant statistical and intrinsic cluster properties, such as the halo mass function and the gas equation of state. Results of these calculations are compared with those obtained from a direct analytic treatment and from hydrodynamical simulations. We show that under certain assumptions on the gas mass fraction our results are consistent with the latest SPT measurement. Our approach can be particularly useful in predicting cluster number counts and their dependence on cluster and cosmological parameters.

We model the formation and evolution of galaxy clusters in the framework of an extended dark matter halo merger-tree algorithm that includes baryons and incorporates basic physical considerations. Our modified treatment is employed to calculate the probability density functions of the halo concentration parameter, intracluster gas temperature, and the integrated Comptonization parameter for different cluster masses and observation redshifts. Scaling relations between cluster mass and these observables are deduced that are somewhat different than previous results. Modeling uncertainties in the predicted probability density functions are estimated. Our treatment and the insight gained from the results presented in this paper can simplify the comparison of theoretical predictions with results from ongoing and future cluster surveys.