The existence of massive stellar black hole binaries (MBHBs), with primary black hole (BH) masses $\ge 31 \, M_\odot$, was proven by the detection of the gravitational wave (GW) event GW150914 during the first LIGO/Virgo observing run (O1), and successively confirmed by seven additional GW signals discovered in the O1 and O2 data. By adopting the galaxy formation model \texttt{GAMESH} coupled with binary population synthesis (BPS) calculations, here we investigate the origin of these MBHBs by selecting simulated binaries compatible in mass and coalescence redshifts. We find that their cosmic birth rates peak in the redshift range $6.5 \leq z \leq 10$, regardless of the adopted BPS. These MBHBs are then old systems forming in low-metallicity ($Z \sim [0.01-0.1] \, Z_{\odot}$), low-stellar-mass galaxies, before the end of cosmic reionization, i.e. significantly beyond the peak of cosmic star formation. GW signals generated by coalescing MBHBs open up new possibilities to probe the nature of stellar populations in remote galaxies, at present too faint to be detected by available electromagnetic facilities.

Coalescences of binary white dwarfs represent a copious source of information for gravitational wave interferometers operating in the decihertz band. Moreover, according to the double degenerate scenario, they have been suggested to be possible progenitors of supernovae (SNe) Type Ia events. In this paper we discuss the detectability of gravitational waves emitted by the inspiral of double white dwarfs. We focus on the constraints that can be derived on the source's luminosity distance, and on other binary's parameters, such as the angular momentum orientation. We explore the possibility of coincident detections of gravitational and electromagnetic signals; the latter comes from the observation of the supernova counterpart. Confirmation of the double degenerate scenario would allow one to use distances inferred in the gravitational wave channel to consistently calibrate SNe as standard candles. We find that decihertz gravitational wave interferometers can measure the luminosity distance with relative accuracy better than $1\%$ for binaries at 100 Mpc. We show how multimessenger observations can put strong constraints on the Hubble constant, which are tighter than current bounds at low redshift, and how they can potentially shed new light on the differences with early-universe measurements.

Here we present a detailed analysis of the properties and evolution of different dwarf galaxies, candidate to host the coalescence of black hole binary systems (BHB) generating GW150914-like events. By adopting a novel theoretical framework coupling the binary population synthesis code \texttt{SeBa} with the Galaxy formation model \texttt{GAMESH}, we can investigate the detailed evolution of these objects in a well resolved cosmological volume of 4~cMpc, having a Milky Way-like (MW) galaxy forming at its center. We identify three classes of interesting candidate galaxies: MW progenitors, dwarf satellites and dwarf galaxies evolving in isolation. We find that: (i) despite differences in individual histories and specific environments the candidates reduce to only nine representative galaxies; (ii) among them, $\sim44\%$ merges into the MW halo progenitors by the redshift of the expected signal, while the remaining dwarfs are found as isolated or as satellites of the MW and their evolution is strongly shaped by both peculiar dynamical history and environmental feedback; (iii) a stringent condition for the environments where GW150914-like binaries can form comes from a combination of the accretion history of their DM halos and the radiative feedback in the high redshift universe; (iv) by comparing with the observed catalogues from DGS and ALLSMOG surveys we find two observed dwarfs respecting the properties predicted by our model. We finally note how the present analysis opens the possibility to build future strategies for host galaxy identification.

Supernovae (SNe) are considered to have a major role in dust enrichment of high redshift galaxies and, due to the short lifetimes of interstellar grains, in dust replenishment of local galaxies. Here we explore how SN dust yields depend on the mass, metallicity, and rotation rate of the progenitor stars, and on the properties of the explosion. To this aim, assuming uniform mixing inside the ejecta, we quantify the dust mass produced by a sample of SN models with progenitor masses $13~M_{\odot} \leq M \leq 120~M_{\odot}$, metallicity $\rm -3 \leq [Fe/H] \leq 0$, rotation rate $\rm v_{\rm rot} = 0$ and $300$~km/s, that explode with a fixed energy of $1.2 \times 10^{51}$~erg (FE models) or with explosion properties calibrated to reproduce the $\rm ^{56}Ni$ - $M$ relation inferred from SN observations (CE models). We find that rotation favours more efficient dust production, particularly for more massive, low metallicity stars, but that metallicity and explosion properties have the largest effects on the dust mass and its composition. In FE models, SNe with $M \leq 20 - 25 ~M_{\odot}$ are more efficient at forming dust: between 0.1 and 1 $M_\odot$ is formed in a single explosion, with a composition dominated by silicates, carbon and magnetite grains when $\rm [Fe/H] = 0$, and by carbon and magnetite grains when $\rm [Fe/H] < 0$. In CE models, the ejecta are massive and metal-rich and dust production is more efficient. The dust mass increases with $M$ and it is dominated by silicates, at all [Fe/H].

With the recent discovery of a dozen dusty star-forming galaxies and around 30 quasars at z>5 that are hyper-luminous in the infrared ($\mu$$L_{\rm IR}>10^{13}$ L$_{\odot}$, where $\mu$ is a lensing magnification factor), the possibility has opened up for SPICA, the proposed ESA M5 mid-/far-infrared mission, to extend its spectroscopic studies toward the epoch of reionization and beyond. In this paper, we examine the feasibility and scientific potential of such observations with SPICA's far-infrared spectrometer SAFARI, which will probe a spectral range (35-230 $\mu$m) that will be unexplored by ALMA and JWST. Our simulations show that SAFARI is capable of delivering good-quality spectra for hyper-luminous infrared galaxies (HyLIRGs) at z=5-10, allowing us to sample spectral features in the rest-frame mid-infrared and to investigate a host of key scientific issues, such as the relative importance of star formation versus AGN, the hardness of the radiation field, the level of chemical enrichment, and the properties of the molecular gas. From a broader perspective, SAFARI offers the potential to open up a new frontier in the study of the early Universe, providing access to uniquely powerful spectral features for probing first-generation objects, such as the key cooling lines of low-metallicity or metal-free forming galaxies (fine-structure and H2 lines) and emission features of solid compounds freshly synthesized by Population III supernovae. Ultimately, SAFARI's ability to explore the high-redshift Universe will be determined by the availability of sufficiently bright targets (whether intrinsically luminous or gravitationally lensed). With its launch expected around 2030, SPICA is ideally positioned to take full advantage of upcoming wide-field surveys such as LSST, SKA, Euclid, and WFIRST, which are likely to provide extraordinary targets for SAFARI.

IR spectroscopy in the range 12-230 micron with the SPace IR telescope for Cosmology and Astrophysics (SPICA) will reveal the physical processes that govern the formation and evolution of galaxies and black holes through cosmic time, bridging the gap between the James Webb Space Telescope (JWST) and the new generation of Extremely Large Telescopes (ELTs) at shorter wavelengths and the Atacama Large Millimeter Array (ALMA) at longer wavelengths. SPICA, with its 2.5-m telescope actively-cooled to below 8K, will obtain the first spectroscopic determination, in the mid-IR rest-frame, of both the star-formation rate and black hole accretion rate histories of galaxies, reaching lookback times of 12 Gyr, for large statistically significant samples. Densities, temperatures, radiation fields and gas-phase metallicities will be measured in dust-obscured galaxies and active galactic nuclei (AGN), sampling a large range in mass and luminosity, from faint local dwarf galaxies to luminous quasars in the distant Universe. AGN and starburst feedback and feeding mechanisms in distant galaxies will be uncovered through detailed measurements of molecular and atomic line profiles. SPICA's large-area deep spectrophotometric surveys will provide mid-IR spectra and continuum fluxes for unbiased samples of tens of thousands of galaxies, out to redshifts of z~6. Furthermore, SPICA spectroscopy will uncover the most luminous galaxies in the first few hundred million years of the Universe, through their characteristic dust and molecular hydrogen features.

Here we investigate the origin of the dust mass (Mdust) observed in the Milky Way (MW) and of dust scaling relations found in a sample of local galaxies from the DGS and KINGFISH surveys. To this aim, we model dust production from Asymptotic Giant Branch (AGB) stars and supernovae (SNe) in simulated galaxies forming along the assembly of a Milky Way-like halo in a well resolved cosmic volume of 4cMpc using the GAMESH pipeline. We explore the impact of different sets of metallicity and mass-dependent AGB and SN dust yields on the predicted Mdust. Our results show that models accounting for grain destruction by the SN reverse shock predict a total dust mass in the MW that is a factor of ~4 lower than observed, and can not reproduce the observed galaxy-scale relations between dust and stellar masses, and dust-to-gas ratios and metallicity, with a smaller discrepancy in galaxies with low metallicity (12 + log(O/H) < 7.5) and low stellar masses (Mstar < 10^7 Msun). In agreement with previous studies, we suggest that competing processes in the interstellar medium must be at play to explain the observed trends. Our result reinforces this conclusion by showing that it holds independently of the adopted AGB and SN dust yields.

We present a novel theoretical model to characterize the formation and coalescence sites of compact binaries in a cosmological context. This is based on the coupling between the binary population synthesis code SeBa with a simulation following the formation of a Milky Way-like halo in a well resolved cosmic volume of 4 cMpc, performed with the GAMESH pipeline. We have applied this technique to investigate when and where systems with properties similar to the recently observed LIGO/VIRGO events are more likely to form and where they are more likely to reside when they coalesce. We find that more than 70% of GW151226 and LVT151012-like systems form in galaxies with stellar mass M* > 10^8 Msun in the redshift range [0.06 - 3] and [0.14 - 11.3], respectively. All GW150914-like systems form in low-metallicity dwarfs with M* < 5 10^6 Msun at 2.4 < z < 4.2. Despite these initial differences, by the time they reach coalescence the observed events are most likely hosted by star forming galaxies with M* > 10^{10} Msun. Due to tidal stripping and radiative feedback, a non negligible fraction of GW150914-like candidates end-up in galaxies with properties similar to dwarf spheroidals and ultra-faint satellites.

Dust grains are classically thought to form in the winds of asymptotic giant branch (AGB) stars. However, there is increasing evidence today for dust formation in supernovae (SNe). To establish the relative importance of these two classes of stellar sources of dust, it is important to know the fraction of freshly formed dust in SN ejecta that is able to survive the passage of the reverse shock and be injected in the interstellar medium. We have developed a new code (GRASH\_Rev) which follows the newly-formed dust evolution throughout the supernova explosion until the merging of the forward shock with the circumstellar ISM. We have considered four well studied SNe in the Milky Way and Large Magellanic Cloud: SN1987A, CasA, the Crab Nebula, and N49. For all the simulated models, we find good agreement with observations and estimate that between 1 and 8$\%$ of the observed mass will survive, leading to a SN dust production rate of $(3.9 \pm 3.7) \times 10^{-4}$ M$_{\odot}$yr$^{-1}$ in the Milky Way. This value is one order of magnitude larger than the dust production rate by AGB stars but insufficient to counterbalance the dust destruction by SNe, therefore requiring dust accretion in the gas phase.

The direct discovery of gravitational waves has finally opened a new observational window on our Universe, suggesting that the population of coalescing binary black holes is larger than previously expected. These sources produce an unresolved background of gravitational waves, potentially observables by ground-based interferometers. In this paper we investigate how modified theories of gravity, modeled using the ppE formalism, affect the expected signal, and analyze the detectability of the resulting stochastic background by current and future ground-based interferometers. We find the constraints that AdLIGO would be able to set on modified theories, showing that they may significantly improve the current bounds obtained from astrophysical observations of binary pulsars.

Dust grains are classically thought to form in the winds of AGB stars. However, nowadays there is increasing evidence for dust formation in SNe. In order to establish the relative importance of these two classes of stellar sources of dust it is important to know what is the fraction of freshly formed dust in SN ejecta that is able to survive the passage of the reverse shock and be injected in the interstellar medium. With this aim, we have developed a new code, GRASH_Rev, that allows to follow the dynamics of dust grains in the shocked SN ejecta and to compute the time evolution of the mass, composition and size distribution of the grains. We consider four well studied SNe in the Milky Way and LMC: SN 1987a, Cas A, the Crab Nebula, and N49. For all the simulated models, we find good agreement with observations. Our study suggests that SN 1987A is too young for the reverse shock to have affected the dust mass. Conversely, in the other three SNe, the reverse shock has already destroyed between 10 and 40% of the initial dust mass. However, the largest dust mass destruction is predicted to occur between 10^3 and 10^5 yr after the explosions. Since the oldest SN in the sample has an estimated age of 4800 yr, current observations can only provide an upper limit to the mass of SN dust that will enrich the interstellar medium, the so-called effective dust yields. We find that only between 1 and 8% of the currently observed mass will survive. This is in good agreement with the values adopted in chemical evolution models which consider the effect of the SN reverse shock. We discuss the astrophysical implications of our results for dust enrichment in local galaxies and at high redshift.

We quantify the role of Population (Pop) III core-collapse supernovae (SNe) as the first cosmic dust polluters. Starting from a homogeneous set of stellar progenitors with masses in the range [13 - 80] Msun, we find that the mass and composition of newly formed dust depend on the mixing efficiency of the ejecta and the degree of fallback experienced during the explosion. For standard Pop III SNe, whose explosions are calibrated to reproduce the average elemental abundances of Galactic halo stars with [Fe/H] < -2.5, between 0.18 and 3.1 Msun (0.39 - 1.76 Msun) of dust can form in uniformly mixed (unmixed) ejecta, and the dominant grain species are silicates. We also investigate dust formation in the ejecta of faint Pop III SN, where the ejecta experience a strong fallback. By examining a set of models, tailored to minimize the scatter with the abundances of carbon-enhanced Galactic halo stars with [Fe/H ] < -4, we find that amorphous carbon is the only grain species that forms, with masses in the range 2.7 10^{-3} - 0.27 Msun (7.5 10^{-4} - 0.11 Msun) for uniformly mixed (unmixed) ejecta models. Finally, for all the models we estimate the amount and composition of dust that survives the passage of the reverse shock, and find that, depending on circumstellar medium densities, between 3 and 50% (10 - 80%) of dust produced by standard (faint) Pop III SNe can contribute to early dust enrichment.

The possible formation of Direct Collapse Black Holes (DCBHs) in the first metal-free atomic cooling halos at high redshifts ($z > 10$) is nowadays object of intense study and several methods to prove their existence are currently under development. The abrupt collapse of a massive ($\sim 10^4 - 10^5 \, \mathrm{M_{\odot}}$) and rotating object is a powerful source of gravitational waves emission. In this work, we employ modern waveforms and the improved knowledge on the DCBHs formation rate to estimate the gravitational signal emitted by these sources at cosmological distances. Their formation rate is very high ($\sim 10^4 \, \mathrm{yr^{-1}}$ up to $z\sim20$), but due to a short duration of the collapse event ($\sim 2-30\, \mathrm{s}$, depending on the DCBH mass) the integrated signal from these sources is characterized by a very low duty-cycle (${\cal D}\sim 10^{-3}$), i.e. a shot-noise signal. Our results show that the estimated signal lies above the foreseen sensitivity of the Ultimate-DECIGO observatory in the frequency range $(0.8-300) \, \mathrm{mHz}$, with a peak amplitude $\Omega_{gw} = 1.1 \times 10^{-54}$ at $\nu_{max} = 0.9 \, \mathrm{mHz}$ and a peak Signal-to-Noise Ratio $\mathrm{SNR}\sim 22$ at $\nu = 20 \, \mathrm{mHz}$. This amplitude is lower than the Galactic confusion noise, generated by binary systems of compact objects in the same frequency band. For this reason, advanced techniques will be required to separate this signal from background and foreground noise components. As a proof-of-concept, we conclude by proposing a simple method, based on the auto-correlation function, to recognize the presence of a ${\cal D} \ll 1$ signal buried into the continuous noise. The aim of this work is to test the existence of a large population of high-z DCBHs, by observing the gravitational waves emitted during their infancy.

We investigate the condition for the formation of low-mass second-generation stars in the early universe. It has been proposed that gas cooling by dust thermal emission can trigger fragmentation of a low-metallicity star-forming gas cloud. In order to determine the critical condition in which dust cooling induces the formation of low-mass stars, we follow the thermal evolution of a collapsing cloud by a one-zone semi-analytic collapse model. Earlier studies assume the dust amount in the local universe, where all refractory elements are depleted onto grains, and/or assume the constant dust amount during gas collapse. In this paper, we employ the models of dust formation and destruction in early supernovae to derive the realistic dust compositions and size distributions for multiple species as the initial conditions of our collapse calculations. We also follow accretion of heavy elements in the gas phase onto dust grains, i.e., grain growth, during gas contraction. We find that grain growth well alters the fragmentation property of the clouds, and that this still does not approach to the value in the local universe. The critical conditions can be written by the gas metallicity Zcr and the initial depletion efficiency fdep,0 of gas-phase metal onto grains, or dust-to-metal mass ratio, as (Zcr/10^{-5.5} Zsun) = (fdep,0/0.18)^{-0.44} with small scatters in the range of Zcr = [0.06--3.2]x10^{-5} Zsun. We also show that the initial dust composition and size distribution are important to determine Zcr.

We investigate the origin of carbon-enhanced metal-poor (CEMP) stars starting from the recently discovered $\rm [Fe/H]<-7.1$ star SMSS J031300 (Keller et al. 2014). We show that the elemental abundances observed on the surface of SMSS J031300 can be well fit by the yields of faint, metal free, supernovae. Using properly calibrated faint supernova explosion models, we study, for the first time, the formation of dust grains in such carbon-rich, iron-poor supernova ejecta. Calculations are performed assuming both unmixed and uniformly mixed ejecta and taking into account the partial destruction by the supernova reverse shock. We find that, due to the paucity of refractory elements beside carbon, amorphous carbon is the only grain species to form, with carbon condensation efficiencies that range between (0.15-0.84), resulting in dust yields in the range (0.025-2.25)M$_{\odot}$. We follow the collapse and fragmentation of a star forming cloud enriched by the products of these faint supernova explosions and we explore the role played by fine structure line cooling and dust cooling. We show that even if grain growth during the collapse has a minor effect of the dust-to-gas ratio, due to C depletion into CO molecules at an early stage of the collapse, the formation of CEMP low-mass stars, such as SMSS J031300, could be triggered by dust cooling and fragmentation. A comparison between model predictions and observations of a sample of C-normal and C-rich metal-poor stars supports the idea that a single common pathway may be responsible for the formation of the first low-mass stars.

In this paper we perform a detailed analysis of the effect of various approximations which have been used in the literature to compute the detection rates of compact binary coalescences for interferometric gravitational wave detectors. We evaluate the detection rates for the coalescence of BH-BH, NS-NS, and BH-NS binaries taking into account their specific statistical properties obtained from population synthesis models, the cosmic star formation rate history, and the effects of redshift on the emitted gravitational wave signals. The results are compared with those obtained with procedures that are based on different levels of approximations, such as adopting averaged values for the total mass and symmetric mass ratio for all the systems of a binary population, using these to compute the horizon distance for individual detectors, or estimating the coalescence rate density within this distance by its local value. We show that these approximations introduce a bias which depends not only on the approximation which is used, but also on the statistical properties of the considered binary population.

We compute the gravitational wave background (GWB) generated by a cosmological population of (BH-BH) binaries using hybrid waveforms recently produced by numerical simulations of (BH-BH) coalescence, which include the inspiral, merger and ring-down contributions. A large sample of binary systems is simulated using the population synthesis code SeBa, and we extract fundamental statistical information on (BH-BH) physical parameters (primary and secondary BH masses, orbital separations and eccentricities, formation and merger timescales). We then derive the binary birth and merger rates using the theoretical cosmic star formation history obtained from a numerical study which reproduces the available observational data at redshifts $z < 8$. We evaluate the contributions of the inspiral, merger and ring-down signals to the GWB, and discuss how these depend on the parameters which critically affect the number of coalescing (BH-BH) systems. We find that Advanced LIGO/Virgo have a chance to detect the GWB signal from the inspiral phase with a $(S/N)=10$ only for the most optimistic model, which predicts the highest local merger rate of 0.85 Mpc$^{-3}$ Myr$^{-1}$. Third generation detectors, such as ET, could reveal the GWB from the inspiral phase predicted by any of the considered models. In addition, ET could sample the merger phase of the evolution at least for models which predict local merger rates between $[0.053 - 0.85]$ Mpc$^{-3}$ Myr$^{-1}$, which are more than a factor 2 lower the the upper limit inferred from the analysis of the LIGO S5 run\cite{Abadieetal2011}.

This paper has been withdrawn by the authors because it had major revisions.

Two classes of high energy sources in our galaxy are believed to host magnetars, neutron stars whose emission results from the dissipation of their magnetic field. The extremely high magnetic field of magnetars distorts their shape, and causes the emission of a conspicuous gravitational waves signal if rotation is fast and takes place around a different axis than the symmetry axis of the magnetic distortion. Based on a numerical model of the cosmic star formation history, we derive the cosmological background of gravitational waves produced by magnetars, when they are very young and fast spinning. We adopt different models for the configuration and strength of the internal magnetic field (which determines the distortion) as well as different values of the external dipole field strength (which governs the spin evolution of magnetars over a wide range of parameters). We find that the expected gravitational wave background differs considerably from one model to another. The strongest signals are generated for magnetars with very intense toroidal internal fields ($\sim 10^{16}$ G range) and external dipole fields of $\sim 10^{14}$, as envisaged in models aimed at explaining the properties of the Dec 2004 giant flare from SGR 1806-20. Such signals should be easily detectable with third generation ground based interferometers such as the Einstein Telescope.

Astrophysical sources emit gravitational waves in a large variety of processes occurred since the beginning of star and galaxy formation. These waves permeate our high redshift Universe, and form a background which is the result of the superposition of different components, each associated to a specific astrophysical process. Each component has different spectral properties and features that it is important to investigate in view of a possible, future detection. In this contribution, we will review recent theoretical predictions for backgrounds produced by extragalactic sources and discuss their detectability with current and future gravitational wave observatories.

Using the results of a numerical simulation which follows the evolution, metal enrichment and energy deposition of both Population III and Population II stars, we predict the redshift dependence of the formation rate of black hole remnants of Population III stars with masses 100- 500Msun and of neutron stars(black holes) remnants of Population II stars with masses 8-20Msun (20-40Msun). We describe the gravitational wave spectrum produced by Population III and Population II sources adopting the most appropriate signals available in the literature and we compute the stochastic backgrounds resulting from the cumulative emission of these sources throughout the history of the Universe. With the aim of assessing whether these backgrounds might act as foregrounds for signals generated in the Inflationary epoch, we compare their amplitudes with the sensitivity of currently planned and future ground/space-based interferometers.

We propose a new method to study the quasi-normal modes of rotating relativistic stars. Oscillations are treated as perturbations in the frequency domain of the stationary, axisymmetric background describing a rotating star. The perturbed quantities are expanded in circular harmonics, and the resulting 2D-equations they satisfy are integrated using spectral methods in the (r,theta)-plane. The asymptotic conditions at infinity, needed to find the mode frequencies, are implemented by generalizing the standing wave boundary condition commonly used in the non rotating case. As a test, the method is applied to find the quasi-normal mode frequencies of a slowly rotating star.

We discuss the possibility that the detection of gravitational waves emitted by compact stars may allow to constrain the MIT bag model of quark matter equation of state. Our results show that the combined knowledge of the frequency of the emitted gravitational wave and of the mass, or the radiation radius, of the source allows one to discriminate between strange stars and neutron stars and set stringent bounds on the bag constants.

We construct models of rotating stars using the perturbative approach introduced by J. Hartle in 1967, and a set of equations of state proposed to model hadronic interactions in the inner core of neutron stars. We integrate the equations of stellar structure to third order in the angular velocity and show, comparing our results to those obtained with fully non linear codes, to what extent third order corrections are needed to accurately reproduce the moment of inertia of a star which rotates at rates comparable to that of the fastest isolated pulsars.