Search

Advances in modern technologies enable the characterisation of exoplanetary atmospheres, most efficiently exploiting the transmission spectroscopy technique. We performed visible (VIS) and near infrared (nIR) high-resolution spectroscopic observations of one transit of HD 149026b, a close-in orbit sub Saturn exoplanet. We first analysed the radial velocity data, refining the value of the projected spin-orbit obliquity. Then we performed transmission spectroscopy, looking for absorption signals from the planetary atmosphere. We find no evidence for H$\alpha$, NaI D2 - D1, MgI and LiI in the VIS and metastable helium triplet HeI(2$^3$S) in the nIR using a line-by-line approach. The non-detection of HeI is also supported by theoretical simulations. With the use of the cross-correlation technique, we do not detect TiI, VI, CrI, FeI and VO in the visible, and CH$_4$, CO$_2$, H$_2$O, HCN, NH$_3$, VO in the nIR. Our non-detection of TiI in the planetary atmosphere is in contrast with a previous detection. We performed injection-retrieval tests, finding that our dataset is sensitive to our TiI model. The non-detection supports the TiI cold-trap theory, which is valid for planets with $T_{\rm eq} <$ 2200 K like HD 149026b. Even if we do not attribute it directly to the planet, we find a possibly significant TiI signal highly redshifted ($\simeq$+20 km s$^{-1}$) with respect to the planetary restframe. Redshifted signals are also found in the FeI and CrI maps. While we can exclude an eccentric orbit to cause it, we investigated the possibility of material accretion falling onto the star, possibly supported by the presence of strong LiI in the stellar spectrum, without finding conclusive results. The analysis of multiple transits datasets could shed more light on this target.

One of the primary sources of gravitational waves (GWs) anticipated to be detected by the Laser Interferometer Space Antenna (LISA) are Galactic double white dwarf binaries (DWDs). However, most of these binaries will be unresolved, and their GWs will overlap incoherently, creating a stochastic noise known as the Galactic foreground. Similarly, the population of unresolved systems in the Milky Way's (MW) satellites is expected to contribute to a stochastic gravitational wave background (SGWB). Due to their anisotropy and the annual motion of the LISA constellation, both the Galactic foreground and the satellite SGWB fall into the category of cyclostationary processes. Leveraging this property, we develop a purely frequency-based method to study LISA's capability to detect the MW foreground and SGWBs from the most promising MW satellites. We analyze both mock data generated by an astrophysically motivated SGWB spectrum, and realistic ones from a DWD population generated via binary population synthesis. We are able to recover or put constrains on the candidate foregrounds, reconstructing -- in the presence of noise uncertainties -- their sky distribution and spectrum. Our findings highlight the significance of the interplay between the astrophysical spectrum and LISA's sensitivity to detect the satellites' SGWB. Considering an astrophysically motivated prior on the satellite positions improves their detectability, which becomes otherwise challenging in the presence of the Galactic foreground. Furthermore, we explore the potential to observe a hypothetical satellite located behind the Galactic disk. Our results suggest that a Large Magellanic Cloud-like satellite could indeed be observable by LISA., Comment: 19 pages, 13 figures, 2 tables

A simple model for the X-ray weakness of JWST-selected broad-line AGNs is proposed under the assumption that the majority of these sources are fed at super-Eddington accretion rates. In these conditions, the hot inner corona above the geometrically thin disk that is responsible for the emission of X-rays in "normal" AGNs will be embedded instead in a funnel-like reflection geometry. The coronal plasma will Compton upscatter optical/UV photons from the underlying thick disk as well as the surrounding funnel walls, and the high soft-photon energy density will cool down the plasma to temperatures in the range 30-40 keV. The resulting X-ray spectra are predicted to be extremely soft, with power-law photon indices Gamma=2.8-4.0, making high-z super-Eddington AGNs largely undetectable by Chandra., Comment: 7 pages, 4 figures, matches version accepted for publication in The Astrophysical Journal Letters

Deep surveys with the James Webb Space Telescope (JWST) have revealed an emergent population of moderate-luminosity, broad-line active galactic nuclei (AGNs) at 4< z< 14 powered by accretion onto early massive black holes. The high number densities reported, together with the large Lyman-continuum (LyC) production efficiency and leakiness into the intergalactic medium (IGM) that are typical of UV-selected AGNs, lead us to reassess a scenario where AGNs are the sole drivers of the cosmic hydrogen/helium reionization process. Our approach is based on the assumptions, grounded in recent observations, that: (a) the fraction of broad-line AGNs among galaxies is around 10-15%; (b) the mean escape fraction of hydrogen LyC radiation is high, >80%, in AGN hosts and is negligible otherwise; and (c) internal absorption at 4 ryd or a steep ionizing EUV spectrum delay full reionization of HeII until z~2.8-3.0, in agreement with observations of the HeII Lyman-alpha forest. In our fiducial models: 1) hydrogen reionization is 99% completed by redshift z~5.3-5.5, and reaches its midpoint at z~6.5-6.7; (2) the integrated Thomson scattering optical depth to reionization is ~0.05, consistent with constraints from cosmic microwave background (CMB) anisotropy data; and (3) the abundant AGN population detected by JWST does not violate constraints on the unresolved X-ray background., Comment: 13 pages, 5 figures, matches version accepted for publication in The Astrophysical Journal

We investigate the relevance of kinematically identified bulges, discs and their role relative to galaxy quenching. We utilize an analysis of the SDSS-MaNGA survey conducted with the GPU-based code BANG which simultaneously models galaxy photometry and kinematics to decompose galaxies into their structural components. Below M~1011 Msun, galaxies exhibit a wide range of dynamical properties, determined by the relative prominence of a dispersion-supported inner region and a rotationally-supported disc. Our analysis reveals a natural separation between these classes, with only a minor fraction of stellar mass retained by structures exhibiting intermediate dynamical support. When examining galaxies in terms of their star formation activity, an apparent decrease in rotational support is observed as they move below the star-forming main sequence. This behaviour is evident with luminosity-weighted tracers of kinematics, while it almost vanishes with mass-weighted tracers. Luminosity-weighted quantities not only capture differences in kinematics but also in the stellar population, potentially leading to biased interpretations of galaxy dynamical properties and quenching. Our findings suggest that quenching does not imply almost any structural transformation in galaxies below M~10^11 Msun. Processes as disc fading more likely account for observed differences in mass-weighted and luminosity-weighted galaxy properties; when the galactic disc ceases star formation, its mass-to-light ratio grows without any significant morphological transformation. The picture is remarkably different above M~10^11 Msun. Regardless of the tracer used, a substantial increase in galaxy dispersion support is observed along with a significant structural change. A different quenching mechanism, most likely associated with mergers, dominates. Notably, this mechanism is confined to a very limited range of high masses., Comment: Accepted for publication in A&A. 16 pages, 10 figures

Motivated by the recent surge in interest concerning white dwarf (WD) planets, this work presents the first numerical exploration of WD-driven atmospheric escape, whereby the high-energy radiation from a hot/young WD can trigger the outflow of the hydrogen-helium envelope for close-in planets. As a pilot investigation, we focus on two specific cases: a gas giant and a sub-Neptune-sized planet, both orbiting a rapidly cooling WD with mass $M_\ast$ = 0.6 \msun\ and separation $a$ = 0.02 AU. In both cases, the ensuing mass outflow rates exceed $10^{14}$ g sec$^{-1}$ for WD temperatures greater than $T_{\rm WD} \simeq$ 50,000 K. At $T_{\rm WD} \simeq$ 18,000 K [/22,000 K], the sub-Neptune [/gas giant] mass outflow rate approaches $10^{12}$ g sec$^{-1}$, i.e., comparable to the strongest outflows expected from close-in planets around late main-sequence stars. Whereas the gas giant remains virtually unaffected from an evolutionary standpoint, atmospheric escape may have sizable effects for the sub-Neptune, depending on its dynamical history, e.g., assuming that the hydrogen-helium envelope makes up 1 [/4] per cent of the planet mass, the entire envelope would be evaporated away so long as the planet reaches 0.02 AU within the first 230 [/130] Myr of the WD formation. We discuss how these results can be generalized to eccentric orbits with effective semi-major axis $a'=a/(1-e^2)^{1/4}$, which receive the same orbit-averaged irradiation. Extended to a much broader parameter space, this approach can be exploited to model the expected demographics of WD planets as a function of their initial mass, composition and migration history, as well as their potential for habitability., Comment: ApJ Letters, in press

Detecting stochastic gravitational wave backgrounds (SGWBs) with The Laser Interferometer Space Antenna (LISA) is among the mission science objectives. Disentangling SGWBs of astrophysical and cosmological origin is a challenging task, further complicated by the noise level uncertainties. In this study, we introduce a Bayesian methodology to infer upon SGWBs, taking inspiration from Gaussian stochastic processes. We investigate the suitability of the approach for signal of unknown spectral shape. We do by discretely exploring the model hyperparameters, a first step towards a more efficient transdimensional exploration. We apply the proposed method to a representative astrophysical scenario: the inference on the astrophysical foreground of Extreme Mass Ratio Inspirals, recently estimated in~\cite{Pozzoli2023}. We find the algorithm to be capable of recovering the injected signal even with large priors, while simultaneously providing estimate of the noise level., Comment: 16 pages, 10 figures, 0 tables

We present a novel implementation for active galactic nucleus (AGN) feedback through ultra-fast winds in the code gizmo. Our feedback recipe accounts for the angular dependence of radiative feedback upon black hole spin. We self-consistently evolve in time i) the gas accretion process from resolved scales to an unresolved AGN disc, ii) the evolution of the spin of the massive black hole (MBH), iii) the injection of AGN-driven winds into the resolved scales, and iv) the spin-induced anisotropy of the overall feedback process. We test our implementation by following the propagation of the wind-driven outflow into a homogeneous medium, and we compare the results against simple analytical models. Then, we consider an isolated galaxy setup and there we study the impact of the AGN feedback on the evolution of the MBH and the of the host galaxy. We find that: i) AGN feedback limits the gas inflow that powers the MBH, with a consequent weak impact on the host galaxy characterized by a star formation (SF) suppression of about a factor of two in the nuclear region; ii) the impact of AGN feedback on the host galaxy and on MBH growth is primarily determined by the AGN luminosity, rather than by its angular pattern set by the MBH spin; iii) the imprint of the angular pattern of the AGN radiation emission manifest in a more clear way at high accretion rates. At such high rates the more isotropic angular patterns, proper to higher spin values, sweep away gas in the nuclear region more easily, hence causing a slower MBH mass and spin growths and a higher quenching of the SF. We argue that the influence of spin-dependent anisotropy of AGN feedback on MBH and galaxy evolution is likely to be relevant in those scenarios characterized by high and prolonged MBH accretion episodes and by high AGN wind-galaxy coupling. Such conditions are more frequently met in galaxy mergers and/or high redshift galaxies.

Absorption of stellar X-ray and Extreme Ultraviolet radiation in the upper atmosphere of close-in exoplanets can give rise to hydrodynamic outflows, which may lead to the gradual shedding of their primordial, light element envelopes. Excess absorption by neutral helium atoms in the metastable state has recently emerged as a viable diagnostic of atmospheric escape. Here we present a public module to the 1D photo-ionization hydrodynamic code ATES, designed to calculate the HeI triplet transmission probability for a broad range of planetary parameters. By relaxing the isothermal outflow assumption, the code enables a self-consistent assessment of the HeI triplet absorption depth along with the atmospheric mass loss rate and the outflow temperature profile, which strongly affects the recombination rate of HeII into HeI triplet. We investigate how the transit signal can be expected to depend upon known system parameters, including host spectral type, orbital distance, as well as planet gravity. At variance with previous studies, which identified K-type stars as favorable hosts, we conclude that late M-dwarfs with Neptune-sized planets orbiting at ~0.05-0.1 AU can be expected to yield the strongest transit signal well in excess of 30% for near-cosmological He/H abundances. More generally, we show that the physics which regulates the population and depletion of the metastable state, combined with geometrical effects, can yield somewhat counter-intuitive results, such as a non-monotonic dependence of the transit depth on orbital distance. These are compounded by a strong degeneracy between the stellar EUV flux intensity and the atmospheric He/H abundance, both of which are highly uncertain. Compared against spectroscopy data our modelling suggests that either a large fraction of the targets have helium depleted envelopes, or, that the input stellar EUV spectra are systematically overestimated., Comment: A&A accepted. The ATES code and the Transmission Probability Module (TPM) can be found at https://github.com/AndreaCaldiroli/ATES-Code

In this work, we use general relativistic magnetohydrodynamics simulations to explore the effect of spin orientation on the dynamics of gas in the vicinity of merging black holes. We present a suite of eight simulations of unequal-mass, spinning black hole binaries embedded in magnetized clouds of matter. Each binary evolution covers approximately 15 orbits before the coalescence. The geometry of the accretion flows in the vicinity of the black holes is significantly altered by the orientation of the individual spins with respect to the orbital angular momentum, with the primary black hole dominating the mass accretion rate $\dot{M}$. We observe quasiperiodic modulations of $\dot{M}$ in most of the configurations, whose amplitude is dependent on the orientation of the black hole spins. We find the presence of a relation between the average amplitude of $\dot{M}$ and the spin precession parameter $\chi_{\mathrm{p}}$ showing that spin misalignment systematically leads to stronger modulation, whereas configurations with spins aligned to the orbital angular momentum damp out the quasiperiodicity. This finding suggests a possible signature imprinted in the accretion luminosity of precessing binaries approaching merger and has possible consequences on future multimessenger observations of massive binary black hole systems., Comment: 12 pages, 8 figures. Submitted to PRD

We present the results obtained from follow-up observations of the MUSE Ultra Deep Field (MUDF) at X-ray energies with XMM-Newton. The MUDF is centred on a unique field with two bright, physically associated quasars at $z\simeq3.23$, separated by $\sim$500 kpc in projection. Both quasars are embedded within extended Ly$\alpha$ nebulae ($\gtrsim 100~\rm kpc$ at a surface brightness flux level of $\approx 6\times 10^{-19} \rm erg~s^{-1}~cm^{-2}~arcsec^{-2}$), whose elongated morphology is suggestive of an extended filament connecting the quasar haloes. The new X-ray observations presented here allow us to characterise the physical properties (e.g. X-ray slope, luminosities, gas column densities) in the innermost region of the MUDF quasars. We find that both quasars are X-ray underluminous compared to objects at similar ultraviolet luminosities. Based on our X-ray spectral analysis, absorbing columns of $N_H(z)\gtrsim$ 10$^{23}$ cm$^{-2}$ appear unlikely, therefore such a weakness is possibly intrinsic. When also including literature data, we do not observe any detectable trend between the area of the nebulae and nuclear luminosities at both the rest-frame 2 keV and 2500 $\rm \mathring{A}$. The area is also not correlated with the X-ray photon index nor with the integrated band flux in the hard band (2$-$10 keV). We also do not find any trend between the extended Ly$\alpha$ emission of the nebulae and the nuclear X-ray luminosity. Finally, the properties of the MUDF quasars' nebulae are consistent with the observed relation between the Ly$\alpha$ integrated luminosity of the nebulae and their area. Our results suggest that the quasar ionization power is not a strong driver of the morphology and size of the nebulae., Comment: 15 pages, 9 figures, reference added, published in MNRAS

We investigate the impact of field-to-field variation, deriving from cosmic variance, in measured Lyman-$\alpha$ emitter (LAE) luminosity functions (LFs) and this variation's impact on inferences of the neutral fraction of the intergalactic medium (IGM) during reionization. We post-process a z=7 IGM simulation to populate the dark matter halos with LAEs. These LAEs have realistic UV magnitudes, Ly$\alpha$ fluxes, and Ly$\alpha$ line profiles. We calculate the attenuation of Ly$\alpha$ emission in universes with varying IGM neutral fraction, $\bar{\rm{x}}_{\rm{HI}}$. In a $\bar{\rm{x}}_{\rm{HI}}=0.3$ simulation, we perform 100 realizations of a mock 2 square degree survey with a redshift window $\Delta z = 0.5$ and flux limit $\rm{f}_{Ly\alpha}>1\times10^{-17}\:\rm{ergs}\:\: \rm{s}^{-1} \: \rm{cm}^{-2}$; such a survey is typical in depth and volume of the largest LAE surveys conducted today. For each realization, we compute the LAE LF and use it to recover the input $\bar{\rm{x}}_{\rm{HI}}$. Comparing the inferred values of $\bar{\rm{x}}_{\rm{HI}}$ across the ensemble of the surveys, we find that cosmic variance, deriving from large-scale structure and variation in the neutral gas along the sightline, imposes a floor in the uncertainty of $\Delta \bar{\rm{x}}_{\rm{HI}} \sim 0.2$ when $\bar{\rm{x}}_{\rm{HI}}$ $=0.3$. We explore mitigation strategies to decrease this uncertainty, such as increasing the volume, decreasing the flux limit, or probing the volume with many independent fields. Increasing the area and/or depth of the survey does not mitigate the uncertainty, but composing a survey with many independent fields is effective. This finding highlights the best strategy for LAE surveys aiming at constraining $\bar{\rm{x}}_{\rm{HI}}$ of the universe during reionization., Comment: 17 pages, 13 figures

From a purely photometric perspective galaxies are generally decomposed into a bulge+disc system, with bulges being dispersion-dominated and discs rotationally-supported. However, recent observations have demonstrated that such a framework oversimplifies complexity, especially if one considers galaxy kinematics. To address this issue we introduced with the GPU-based code \textsc{bang} a novel approach that employs analytical potential-density pairs as galactic components, allowing for a computationally fast, still reliable fit of the morphological and kinematic properties of galaxies. Here we apply \textsc{bang} to the SDSS-MaNGA survey, estimating key parameters such as mass, radial extensions, and dynamics, for both bulges and discs of +10,000 objects. We test our methodology against a smaller subsample of galaxies independently analysed with an orbit-based algorithm, finding agreement in the recovered total stellar mass. We also manage to reproduce well-established scaling relations, demonstrating how proper dynamical modelling can result in tighter correlations and provide corrections to standard approaches. Finally, we propose a more general way of decomposing galaxies into "hot" and "cold" components, showing a correlation with orbit-based approaches and visually determined morphological type. Unexpected tails in the "hot-to-total" mass-ratio distribution are present for galaxies of all morphologies, possibly due to visual morphology misclassifications., Comment: 13 pages, 10 figures, accepted for pubblication on MNRAS

The dozens of rocky exoplanets discovered in the Circumstellar Habitable Zone (CHZ) currently represent the most suitable places to host life as we know it outside the Solar System. However, the presumed presence of liquid water on the CHZ planets does not guarantee suitable environments for the emergence of life. According to experimental studies, the building blocks of life are most likely produced photochemically in presence of a minimum ultraviolet (UV) flux. On the other hand, high UV flux can be life-threatening, leading to atmospheric erosion and damaging biomolecules essential to life. These arguments raise questions about the actual habitability of CHZ planets around stars other than Solar-type ones, with different UV to bolometric luminosity ratios. By combining the "principle of mediocricy" and recent experimental studies, we define UV boundary conditions (UV-habitable Zone, UHZ) within which life can possibly emerge and evolve. We investigate whether exoplanets discovered in CHZs do indeed experience such conditions. By analysing Swift-UV/Optical Telescope data, we measure the near ultraviolet (NUV) luminosities of 17 stars harbouring 23 planets in their CHZ. We derive an empirical relation between NUV luminosity and stellar effective temperature. We find that eighteen of the CHZ exoplanets actually orbit outside the UHZ, i.e., the NUV luminosity of their M-dwarf hosts is decisively too low to trigger abiogenesis - through cyanosulfidic chemistry - on them. Only stars with effective temperature >3900 K illuminate their CHZ planets with enough NUV radiation to trigger abiogenesis. Alternatively, colder stars would require a high-energy flaring activity., Comment: 8 pages, 3 figures

Primordial black holes (PBH), supposedly formed in the very early Universe, have been proposed as a possible viable dark matter candidate. In this work we characterize the expected gravitational wave (GW) losses from a population of PBHs orbiting Sgr A$^{\star}$, the super-massive black hole at the Galactic center (GC), and assess the signal detectability by the planned space-borne interferometer LISA and by the proposed next generation space-borne interferometer $\mu$Ares. Assuming that PBHs indeed form the entire diffuse mass allowed to reside within the orbit of the S2 star, we compute an upper limit to the expected GW signal both from resolved and non-resolved sources, under the further assumptions of monochromatic mass function and thermally distributed eccentricities. By comparing with our previous work where PBHs on circular orbits were assumed, we show for 1 M$_{\odot}$ PBHs how the GW signal from high harmonics over a 10 year data stream increases by a factor of six the chances of LISA detectability, from the $\approx 10\%$ of the circular case, to $\approx 60\%$, whereas multiple sources can be identified in $20\%$ of our mock populations. The background signal, made by summing up all non resolved sources, should be certainly detectable thanks to the PBHs with higher eccentricity evolving under two body relaxation. In the case of $\mu$Ares, because of its improved sensitivity in the $\mu$Hz band, one third of the entire population of PBHs orbiting Sgr A$^{\star}$ would be resolved. The background noise from the remaining non resolved sources should be detectable as well. Finally we present the results for different PBH masses., Comment: 11 pages, 9 figures

The evolution of the luminosity function (LF) of Active Galactic Nuclei (AGNs) at $z \gtrsim 5$ represents a key constraint to understand their contribution to the ionizing photon budget necessary to trigger the last phase transition in the Universe, i.e. the epoch of Reionization. Recent searches for bright high-z AGNs suggest that the space densities of this population at $z>4$ has to be revised upwards, and sparks new questions about their evolutionary paths. Gas accretion is the key physical mechanism to understand both the distribution of luminous sources and the growth of central Super-Massive Black Holes (SMBHs). In this work, we model the high-z AGN-LF assuming that high-z luminous AGN shine at their Eddington limit: we derive the expected evolution as a function of the ``duty-cycle'' ($f_{\rm dc}$), i.e. the fraction of life-time that a given SMBH spends accreting at the Eddington rate. Our results show that intermediate values ($f_{\rm dc} \simeq 0.1$) predict the best agreement with the ionizing background and photoionization rate, but do not provide enough ionizing photons to account for the observed evolution of the hydrogen neutral fraction. Smaller values ($f_{\rm dc} \lesssim 0.05$) are required for AGNs to be the dominant population responsible for Hydrogen reionization in the Early Universe. We then show that this low-$f_{\rm dc}$ evolution can be reconciled with the current constraints on Helium reionization, although it implies a relatively large number of inactive SMBHs at $z\gtrsim5$, in tension with SMBH growth models based on heavy seeding., Comment: 11 pages, 5 figures, MNRAS accepted

The putative ubiquity of massive black holes (MBH) at the center of galaxies, and the hierarchical progress of structure formation along the cosmic history, together necessarily imply the existence of a large population of cosmic MBH binaries. Such systems are understood to be the loudest sources of gravitational waves at mHz frequencies, the regime that will be probed by the next Laser Interferometer Space Antenna (LISA). It has been proposed that the rate at which MBHs pair and then bind to form binaries is critically dependent upon the feedback exerted by the MBHs on the surrounding gaseous environment. Using the publicly available code GIZMO, we perform a suite of simulations aimed at studying the dynamics of a MBH pair embedded in a gaseous disk on 100 pc scale. By means of dedicated modules, we follow the dynamics of MBHs in the presence of different spin-dependent radiative feedback models, and compare the results to a benchmark case with no feedback at all. Our main finding is that feedback causes the secondary MBH to shrink its orbit at a reduced pace, when compared to models where feedback is absent. Moreover, such slower inspiral occurs on eccentric orbits, as feedback has the net effect of hampering the circularization process. Though idealized in many aspects, our study highlights and quantifies the importance of including spin-dependent feedback recipes in hydrodynamic simulations of MBH pairs, and ultimately in assessing the cosmological coalescence rate of such systems in view of their detection through gravitational waves.

Dissecting the underlying structure of galaxies is of main importance in the framework of galaxy formation and evolution theories. While a classical bulge+disc decomposition of disc galaxies is usually taken as granted, this is only rarely solidly founded upon the full exploitation of the richness of data arising from spectroscopic studies with integral field units. In this work we describe a fully Bayesian estimation method of the global structure of disc galaxies which makes use of the wealth of photometric, kinematic, and mass-to-light ratio data, and that can be seen as a first step towards a machine-learning approach, certainly needed when dealing with larger samples of galaxies. Ours is a novel, hybrid line of action in tackling the problem of galactic parameter estimation, neither purely photometric nor orbit-based. Being rooted on a nested sampler, our code, which is available publicly as an online repository, allows for a statistical assessment of the need for multiple components in the dissecting process. As a first case-study the GPU-optimized code is applied to the S0 galaxy NGC-7683, finding that in this galaxy a pseudo-bulge, possibly the remnant of a bar-like structure, does exist in the center of the system. These results are then tested against the publicly available, orbit-based code DYNAMITE, finding substantial agreement., Comment: 14 pages, 7 figures, accepted for publication in MNRAS

A recent analysis of data from the ESA Gaia mission demonstrated that the kinematics of stars in the Milky Way can be modelled without invoking the presence of dark matter whatsoever. Indeed, the higher-than-Keplerian velocities observed in outer stars can be ascribed to the properties of the general relativistic (GR) metric assumed to describe the Galaxy. Here we generalize the concept, and derive the most general exact GR model for a stationary axi-symmetric dynamically cold dust structure. We explicitly show how deviations from the commonly adopted Newtonian dynamics are indeed manifest even at low velocities and low densities. We provide for the first time a detailed description of the frequency shift experienced by photons travelling from any emission site within an external disc galaxy to the observer, relating the outcome of the shift measurement to the gravitational properties of the galaxy. Finally, we devise a novel, groundbreaking observational test potentially able to fully characterize the GR metric under the minimal set of assumptions mentioned above. The proposed experiment exploits the effects non-diagonal GR terms have on the frequency shift of photons, ultimately providing a test to evaluate whether dark matter is actually required by disc galaxy kinematics., Comment: 7 pages, 0 figures

The Laser Interferometer Space Antenna (LISA) will be a transformative experiment for gravitational wave astronomy, and, as such, it will offer unique opportunities to address many key astrophysical questions in a completely novel way. The synergy with ground-based and space-born instruments in the electromagnetic domain, by enabling multi-messenger observations, will add further to the discovery potential of LISA. The next decade is crucial to prepare the astrophysical community for LISA's first observations. This review outlines the extensive landscape of astrophysical theory, numerical simulations, and astronomical observations that are instrumental for modeling and interpreting the upcoming LISA datastream. To this aim, the current knowledge in three main source classes for LISA is reviewed; ultracompact stellar-mass binaries, massive black hole binaries, and extreme or intermediate mass ratio inspirals. The relevant astrophysical processes and the established modeling techniques are summarized. Likewise, open issues and gaps in our understanding of these sources are highlighted, along with an indication of how LISA could help making progress in the different areas. New research avenues that LISA itself, or its joint exploitation with upcoming studies in the electromagnetic domain, will enable, are also illustrated. Improvements in modeling and analysis approaches, such as the combination of numerical simulations and modern data science techniques, are discussed. This review is intended to be a starting point for using LISA as a new discovery tool for understanding our Universe.

In this work we characterize the expected gravitational wave signal detectable by the planned space-borne interferometer LISA and the proposed next generation space-borne interferometer $\mu$Ares arising from a population of primordial black holes orbiting Sgr A*, the super-massive black hole at the Galactic center. Assuming that such objects indeed form the entire diffuse mass allowed by the observed orbit of S2 in the Galactic center, under the simplified assumption of circular orbits and monochromatic mass function, we assess the expected signal in gravitational waves, either from resolved and non-resolved sources. We estimate a small but non negligible chance of $\simeq$ 10% of detecting one single 1 M$_{\odot}$ primordial black hole with LISA in a 10-year-long data stream, while the background signal due to unresolved sources would essentially elude any reasonable chance of detection. On the contrary, $\mu$Ares, with a $\simeq$ 3 orders-of-magnitude better sensitivity at $\simeq$ 10$^{-5}$ Hz, would be able to resolve $\simeq$ 140 solar mass primordial black holes in the same amount of time, while the unresolved background should be observable with an integrated signal-to-noise ratio $\gtrsim$ 100. Allowing the typical PBH mass to be in the range 0.01-10 M$_{\odot}$ would increase LISA chance of detection to $\simeq$ 40% towards the lower limit of the mass spectrum. In the case of $\mu$Ares, instead, we find a "sweet spot" just about 1 M$_{\odot}$, a mass for which the number of resolvable events is indeed maximized., Comment: 10 pages, 11 figures; Referee reviewed

We present general relativistic magneto-hydrodynamical simulations of equal-mass spinning black hole binary mergers embedded in a magnetized gas cloud. We focus on the effect of the spin orientation relative to the orbital angular momentum on the flow dynamics, mass accretion rate and Poynting luminosity. We find that, across the inspiral, the gas accreting onto the individual black holes concentrates into disk-like overdensities, whose angular momenta are oriented towards the spin axes and which persist until merger. We identify quasi-periodic modulations occurring in the mass accretion rate at the level of 1-20%, evolving in parallel with the gravitational wave chirp. The similarity between the accretion rate time-series and the gravitational strain is a consequence of the interplay between strong, dynamical gravitational fields and magnetic fields in the vicinity of the inspiralling black holes. This result suggests that quasi-periodicity in the pre-merger accretion rate of massive binaries is not exclusive of environments in which the black holes are embedded in a circumbinary accretion disk, and could provide an additional useful signature of electromagnetic emission concurrent to low-frequency gravitational wave detection., Comment: 9 pages, 4 figures. Accepted for publication in The Astrophysical Journal Letters

Making use of the publicly available 1D photoionization hydrodynamics code ATES we set out to investigate the combined effects of planetary gravitational potential energy ($\phi_p\equiv GM_p/R_p$) and stellar X-ray and Extreme Ultraviolet (XUV) irradiation ($F_{\rm XUV}$) on the evaporation efficiency ($\eta$) of moderately-to-highly irradiated gaseous planets, from sub-Neptunes through hot Jupiters. We show that the (known) existence of a threshold potential above which energy-limited escape (i.e., $\eta\simeq 1$) is unattainable can be inferred analytically. For $\log \phi_p\gtrsim \log \phi_p^{\rm thr}\approx [12.9-13.2]$ (in cgs units), most of the energy absorption occurs where the average kinetic energy acquired by the ions through photo-electron collisions is insufficient for escape. This causes the evaporation efficiency to plummet with increasing $\phi_p$,. Whether or not planets with $\phi_p\lesssim \phi_p^{\rm thr}$ exhibit energy-limited outflows is regulated primarily by the stellar irradiation level. Specifically, for low-gravity planets, above $F_{\rm XUV}\simeq 10^{4-5}$ erg cm$^{-2}$s$^{-1}$ Ly$\alpha$ losses overtake adiabatic and advective cooling and the evaporation efficiency of low-gravity planets drops below the energy-limited approximation, albeit remaining largely independent of $\phi_p$Further, we show that whereas $\eta$ increases as $F_{\rm XUV}$ increases for planets above $\phi^{\rm thr}_p$, the opposite is true for low-gravity planets. This behavior can be understood by examining the relative fractional contributions of advective and radiative losses as a function of atmospheric temperature. This novel framework enables a reliable, physically motivated prediction of the expected evaporation efficiency for a given planetary system; an analytical approximation of the best-fitting $\eta$ is given in the appendix., Comment: 12 pages, 8 figures. Accepted for publication in A&A. In the new version, section 4 and 5 have been completely revised after the referee revision. Comments are welcome

Intense X-ray and ultraviolet stellar irradiation can heat and inflate the atmospheres of closely orbiting exoplanets, driving mass outflows that may be significant enough to evaporate a sizable fraction of the planet atmosphere over the system lifetime. The recent surge in the number of known exoplanets, together with the imminent deployment of new ground and space-based facilities for exoplanet discovery and characterization, requires a prompt and efficient assessment of the most promising targets for intensive spectroscopic follow-ups. To this purpose, we developed ATES (ATmospheric EScape); a new hydrodynamics code that is specifically designed to compute the temperature, density, velocity and ionization fraction profiles of highly irradiated planetary atmospheres, along with the current, steady-state mass loss rate. ATES solves the one-dimensional Euler, mass and energy conservation equations in radial coordinates through a finite-volume scheme. The hydrodynamics module is paired with a photoionization equilibrium solver that includes cooling via bremsstrahlung, recombination and collisional excitation/ionization for the case of a primordial atmosphere entirely composed of atomic hydrogen and helium, whilst also accounting for advection of the different ion species. Compared against the results of 14 moderately-to-highly irradiated planets simulated with The PLUTO-CLOUDY Interface (TPCI), ATES yields remarkably good agreement at a significantly smaller fraction of the computational time. A convergence study shows that ATES recovers stable, steady-state hydrodynamic solutions for systems with $\log(-\phi_p) \lesssim 12.9 + 0.17\log F_{\rm XUV}$. Incidentally, atmospheres of systems above this threshold are generally thought to be undergoing Jeans escape. The code, which also features a user-friendly graphic interface, is available publicly as an online repository., Comment: 17 pages, 11 figures. Revised version. Accepted by A&A. Minor changes only; this new version includes a convergence study. Comments are welcome

We present a definitive assessment of the role of Inverse Compton scattering of Cosmic Microwave Background photons (IC/CMB) in the context of radio galaxies. Owing to the steep increase of the CMB radiation energy density, IC/CMB is supposed to become progressively more important with respect to radio synchrotron cooling as the redshift increases. For typical energies at play, this process will up-scatter the CMB photons into the X-ray band, and is thus expected to yield a redshift-dependent, concurrent X-ray brightening and radio dimming of the jet-powered structures. Here we show how a conclusive proof of this effect hinges on high-resolution imaging data in which the extended lobes can be distinguished from the compact hot spots where synchrotron-self-Compton dominates the X-ray emission regardless of redshift. We analyze Chandra and Very Large Array data of 11 radio galaxies between 1.33 radio galaxies below <1 mJy, IC/CMB alone can not be responsible for a deficit in high-z, radio-loud AGN if--as we argue--such AGN typically have bright, compact hot spots., Comment: 16 pages, 5 figures, 5 tables. Accepted for publication in MNRAS

Intermediate mass black holes (IMBHs) in the mass range $10^2-10^5\,\mathrm{M_{\odot}}$ bridge the gap between stellar black holes (BHs) and supermassive BHs. Here, we investigate the possibility that IMBHs form in young star clusters via runaway collisions and BH mergers. We analyze $10^4$ simulations of dense young star clusters, featuring up-to-date stellar wind models and prescriptions for core collapse and (pulsational) pair instability. In our simulations, only 9 IMBHs out of 218 form via binary BH mergers, with a mass $\sim{}100-140$ M$_\odot$. This channel is strongly suppressed by the low escape velocity of our star clusters. In contrast, IMBHs with masses up to $\sim{}438$ M$_{\odot}$ efficiently form via runaway stellar collisions, especially at low metallicity. Up to $\sim{}0.2$~% of all the simulated BHs are IMBHs, depending on progenitor's metallicity. The runaway formation channel is strongly suppressed in metal-rich ($Z=0.02$) star clusters, because of stellar winds. IMBHs are extremely efficient in pairing with other BHs: $\sim{}70$% of them are members of a binary BH at the end of the simulations. However, we do not find any IMBH-BH merger. More massive star clusters are more efficient in forming IMBHs: $\sim{}8$% ($\sim{}1$%) of the simulated clusters with initial mass $10^4-3\times{}10^4$ M$_\odot$ ($10^3-5\times{}10^3$ M$_\odot$) host at least one IMBH., Comment: 12 pages, 16 figures, 2 tables. Comments welcome

We perform the first suite of fully general relativistic magnetohydrodynamic simulations of spinning massive black hole binary mergers. We consider binary black holes with spins of different magnitudes aligned to the orbital angular momentum, which are immersed in a hot, magnetized gas cloud. We investigate the effect of the spin and degree of magnetization (defined through the fluid parameter $\beta^{-1}\equiv p_{\mathrm{mag}}/p_{\mathrm{fluid}}$) on the properties of the accretion flow. We find that magnetized accretion flows are characterized by more turbulent dynamics, as the magnetic field lines are twisted and compressed during the late inspiral. Post-merger, the polar regions around the spin axis of the remnant Kerr black hole are magnetically dominated, and the magnetic field strength is increased by a factor $\sim$10$^2$ (independently from the initial value of $\beta^{-1}$). The magnetized gas in the equatorial plane acquires higher angular momentum, and settles in a thin circular structure around the black hole. We find that mass accretion rates of magnetized configurations are generally smaller than in the unmagnetized cases by up to a factor $\sim$3. Black hole spins have also a suppressing effect on the accretion rate, as large as $\sim$48\%. As a potential driver for electromagnetic emission we follow the evolution of the Poynting luminosity, which increases after merger up to a factor $\sim2$ with increasing spin, regardless of the initial level of magnetization of the fluid. Our results stress the importance of taking into account both spins and magnetic fields when studying accretion processes onto merging massive black holes., Comment: 15 pages, 12 figures

The time frame in which hydrogen reionization occurred is highly uncertain, but can be constrained by observations of Lyman-alpha (Ly$\alpha$) emission from distant sources. Neutral hydrogen in the intergalactic medium (IGM) attenuates Ly$\alpha$~photons emitted by galaxies. As reionization progressed the IGM opacity decreased, increasing Ly$\alpha$~visibility. The galaxy Ly$\alpha$~luminosity function (LF) is thus a useful tool to constrain the timeline of reionization. In this work, we model the Ly$\alpha$~LF as a function of redshift, $z=5-10$, and average IGM neutral hydrogen fraction, $\overline{x}_\textsc{hi}$. We combine the Ly$\alpha$~luminosity probability distribution obtained from inhomogeneous reionization simulations with a model for the UV LF to model the Ly$\alpha$~LF. As the neutral fraction increases, the average number density of Ly$\alpha$~emitting galaxies decreases, and are less luminous, though for $\overline{x}_\textsc{hi} \lesssim 0.4$ there is only a small decrease of the Ly$\alpha$~LF. We use our model to infer the IGM neutral fraction at $z=6.6, 7.0, 7.3$ from observed Ly$\alpha$~LFs. We conclude that there is a significant increase in the neutral fraction with increasing redshift: $\overline{x}_\textsc{hi}(z=6.6)=0.08^{+ 0.08}_{- 0.05}, \, \overline{x}_\textsc{hi}(z=7.0)=0.28 \pm 0.05$ and $\overline{x}_\textsc{hi}(z=7.3)=0.83^{+ 0.06}_{- 0.07}$. We predict trends in the Ly$\alpha$~luminosity density and Schechter parameters as a function of redshift and the neutral fraction. We find that the Ly$\alpha$~luminosity density decreases as the universe becomes more neutral. Furthermore, as the neutral fraction increases, the faint-end slope of the Ly$\alpha$~LF steepens, and the characteristic Ly$\alpha$~luminosity shifts to lower values, concluding that the evolving shape of the Ly$\alpha$~LF -- not just its integral -- is an important tool to study reionization., Comment: 20 pages, 10 figures, submitted to ApJ

Third Generation ground based Gravitational Wave Interferometers, like the Einstein Telescope (ET), Cosmic Explorer (CE), and the Laser Interferometer Space Antenna (LISA) will detectcoalescing binary black holes over a wide mass spectrum and across all cosmic epochs. We track the cosmological growth of the earliest light and heavy seeds that swiftly transit into the supermassive domain using a semi analytical model for the formation of quasars at $z=6.4$, 2 and $0.2$, in which we follow black hole coalescences driven by triple interactions. We find that light seed binaries of several $10^2$ M$_\odot$ are accessible to ET with a signal-to-noise ratio ($S/N$) of $10-20$ at $620$). Mergers involving heavy seeds ($\sim 10^5 M_\odot - 10^6 M_\odot$) would be within reach up to $z=20$ in the LISA frequency domain. The lower-z model predicts $11.25(18.7)$ ET(LISA) events per year, overall., Comment: 16 pages, 8 figures accepted for publication in MNRAS

We report on deep Chandra X-ray Telescope imaging observations of 4C 63.20, one of the few known radio galaxies at z>3.5. The X-ray counterpart is resolved into a core plus two off-nuclear sources that (combined) account for close to 30% of the total X-ray flux. Their morphology and orientation are consistent with a diffuse, lobe-like nature, albeit compact hotspots cannot be ruled out. The broadband spectral energy distribution of 4C 63.20 can be reproduced with a jet model where the majority of the radio flux can be ascribed to synchrotron emission from the hotspots, whereas the (non-nuclear) X-ray emission is produced via Inverse Compton (IC) off of Cosmic Microwave Background (CMB) photons within the extended lobes. This scenario is broadly consistent with the expectation from highly magnetized lobes in a hotter CMB, and supports the view that IC/CMB may quench less extreme radio lobes at high redshifts., Comment: 10 pages, 8 figures, accepted to MNRAS

We present a new constraint on the Lyman Continuum (LyC) escape fraction at z~1.3. We obtain deep, high sensitivity far-UV imaging with the Advanced Camera for Surveys (ACS) Solar Blind Channel (SBC) on the Hubble Space Telescope (HST), targeting 11 star-forming galaxies at 1.2 190 \AA, low stellar mass (M* < 10^10 M_sun) and U-band magnitude of U<24.2. These criteria identify young, low metallicity star bursting populations similar to the primordial star-forming galaxies believed to have reionized the universe. We do not detect any LyC signal (with S/N >3) in the individual galaxies or in the stack in the far-UV images. We place $3\sigma$ limits on the relative escape fraction of individual galaxies to be f_{esc,rel}<[0.10-0.22] and a stacked $3\sigma$ limit of f_{esc,rel}<0.07. Comparing to the confirmed LyC emitters from the literature, the galaxies in our sample span similar ranges of various galaxy properties including stellar mass, dust attenuation, and star formation rate (SFR). In particular, we compare the distribution of H$\alpha$ and [OIII] EWs of confirmed LyC emitters and non-detections including the galaxies in this study. Finally, we discuss if a dichotomy seen in the distribution of H$\alpha$ EWs can perhaps distinguish the LyC emitters from the non-detections., Comment: 16 pages, 10 Figures, 2 Tables

Young star clusters are the most common birth-place of massive stars and are dynamically active environments. Here, we study the formation of black holes (BHs) and binary black holes (BBHs) in young star clusters, by means of 6000 N-body simulations coupled with binary population synthesis. We probe three different stellar metallicities (Z=0.02, 0.002 and 0.0002) and two initial density regimes (density at the half-mass radius $\rho_{\rm h}\ge{}3.4\times10^4$ and $\ge{1.5\times10^2}$ M$_\odot$ pc$^{-3}$ in dense and loose star clusters, respectively). Metal-poor clusters tend to form more massive BHs than metal-rich ones. We find $\sim{}6$, $\sim{}2$, and $<1$% of BHs with mass $m_{\rm BH}>60$ M$_\odot$ at Z=0.0002, 0.002 and 0.02, respectively. In metal-poor clusters, we form intermediate-mass BHs with mass up to $\sim{}320$ M$_\odot$. BBH mergers born via dynamical exchanges (exchanged BBHs) can be more massive than BBH mergers formed from binary evolution: the former (latter) reach total mass up to $\sim{}140$ M$_\odot$ ($\sim{}80$ M$_\odot$). The most massive BBH merger in our simulations has primary mass $\sim{}88$ M$_\odot$, inside the pair-instability mass gap, and a mass ratio of $\sim{}0.5$. Only BBHs born in young star clusters from metal-poor progenitors can match the masses of GW170729, the most massive event in O1 and O2, and those of GW190412, the first unequal-mass merger. We estimate a local BBH merger rate density $\sim{}110$ and $\sim{}55$ Gpc$^{-3}$ yr$^{-1}$, if we assume that all stars form in loose and dense star clusters, respectively., Comment: 14 pages, 10 figures, 5 tables, accepted for publication in MNRAS

We study the dynamical evolution of eccentric massive black hole binaries (MBHBs) interacting with unbound stars by means of an extensive set of three body scattering experiments. Compared to previous studies, we extend the investigation down to a MBHB mass ratio of $q=m_2/m_1=10^{-4}$, where $m_1$ and $m_2$ are the masses of the primary and secondary hole respectively. Contrary to a simple extrapolation from higher mass ratios, we find that for $q\lesssim 10^{-3}$ the eccentricity growth rate becomes negative, i.e., the binary {\it circularises} as it shrinks. This behaviour is due to the subset of interacting stars captured in metastable counter-rotating orbits; those stars tend to extract angular momentum from the binary, promoting eccentricity growth for $q>10^{-3}$, but tend to inject angular momentum into the binary driving it towards circularisation for $q<10^{-3}$. The physical origin of this behaviour requires a detailed study of the orbits of this subset of stars and is currently under investigation. Our findings might have important consequences for intermediate MBHs (IMBHs) inspiralling onto MBHs (e.g. a putative $10^3\rm M_{\odot}$ black hole inspiralling onto SgrA$^*$)., Comment: 6 pages, 5 figures; accepted for publication in MNRAS

Pair instability (PI) and pulsational PI prevent the formation of black holes (BHs) with mass $\gtrsim{}60$ M$_\odot$ from single star evolution. Here, we investigate the possibility that BHs with mass in the PI gap form via stellar mergers and multiple stellar mergers, facilitated by dynamical encounters in young star clusters. We analyze $10^4$ simulations, run with the direct N-body code nbody6++gpu coupled with the population synthesis code MOBSE. We find that up to $\sim{}6$~% of all simulated BHs have mass in the PI gap, depending on progenitor's metallicity. This formation channel is strongly suppressed in metal-rich (Z = 0.02) star clusters, because of stellar winds. BHs with mass in the PI gap are initially single BHs but can efficiently acquire companions through dynamical exchanges. We find that $\sim{}$21%, 10% and 0.5% of all binary BHs have at least one component in the PI mass gap at metallicity Z = 0.0002, 0.002 and 0.02, respectively. Based on the evolution of the cosmic star formation rate and metallicity, and under the assumption that all stars form in young star clusters, we predict that $\sim{}5$~% of all binary BH mergers detectable by advanced LIGO and Virgo at their design sensitivity have at least one component in the PI mass gap., Comment: 7 pages, 3 figures, 2 tables, MNRAS, in press

Searching for distinctive signatures, which characterize different formation channels of binary black holes (BBHs), is a crucial step towards the interpretation of current and future gravitational wave detections. Here, we investigate the demography of merging BBHs in young star clusters (SCs), which are the nursery of massive stars. We performed $4\times{} 10^3$ N-body simulations of SCs with metallicity $Z=0.002$, initial binary fraction $0.4$ and fractal initial conditions, to mimic the clumpiness of star forming regions. Our simulations include a novel population-synthesis approach based on the code MOBSE. We find that SC dynamics does not affect the merger rate significantly, but leaves a strong fingerprint on the properties of merging BBHs. More than 50 % of merging BBHs in young SCs form by dynamical exchanges in the first few Myr. Dynamically formed merging BBHs are significantly heavier than merging BBHs in isolated binaries: merging BBHs with total mass up to $\sim{}120$ M$_\odot$ form in young SCs, while the maximum total mass of merging BBHs in isolated binaries with the same metallicity is only $\sim{}70$ M$_\odot$. Merging BBHs born via dynamical exchanges tend to have smaller mass ratios than BBHs in isolated binaries. Furthermore, SC dynamics speeds up the merger: the delay time between star formation and coalescence is significantly shorter in young SCs. In our simulations, massive systems such as GW170729 form only via dynamical exchanges. Finally $\sim{}2$ % of merging BBHs in young SCs have mass in the pair-instability mass gap ($\sim{}60-120$ M$_\odot$). This represents a unique fingerprint of merging BBHs in SCs., Comment: 15 pages, 11 figures, 5 tables

Coalescing massive black hole binaries (MBHBs) of $10^{4-7} \rm M_{\odot}$, forming in the aftermath of galaxy mergers, are primary targets of the space mission LISA, the {\it Laser Interferometer Space Antenna}. An assessment of LISA detection prospects requires an estimate of the abundance and properties of MBHBs that form and evolve during the assembly of cosmic structures. To this aim, we employ a semi-analytic model to follow the co-evolution of MBHBs within their host galaxies. We identify three major evolutionary channels driving the binaries to coalescence: two standard paths along which the binary evolution is driven by interactions with the stellar and/or gaseous environment, and a novel channel where MBHB coalescence occurs during the interaction with a third black hole. For each channel, we follow the orbital evolution of MBHBs with physically motivated models that include a self-consistent treatment of the orbital eccentricity. We find that LISA will detect between $\approx 25$ and $\approx 75$ events per year depending on the seed model. We show that triple-induced coalescences can range from a few detected events up to $\sim 30\%$ of the total detected mergers. Moreover, even if the standard gas/stars-driven evolutionary channels should fail and MBHBs were to stall, triple interactions would still occur as a result of the hierarchical nature of galaxy formation, resulting in about $\approx 10$ to $\approx 20$ LISA detections per year. Remarkably, triple interactions among the black holes can produce coalescing binaries with large eccentricities ($\gtrsim 0.9$) upon entrance into the LISA band. This eccentricity will remain significant ($\sim 0.1$) also at merger, requiring suitable templates for parameter estimation., Comment: 17 pages, 10 figures, Accepted MNRAS

Bars are a key factor in the long-term evolution of spiral galaxies, in their unique role in redistributing angular momentum and transporting gas and stars on large scales. The Eris-suite simulations are cosmological zoom-in, N-body, smoothed-particle hydrodynamic simulations built to follow the formation and evolution of a Milky Way-sized galaxy across the build-up of the large scale structure. Here we analyse and describe the outcome of two particular simulations taken from the Eris suite - ErisBH and Eris2k - which mainly differ in the prescriptions employed for gas cooling, star formation, and feedback from supernovae and black holes. Our study shows that the enhanced effective feedback in Eris2k, due to the collective effect of the different micro-physics implementations, results in a galaxy which is less massive than its ErisBH counterpart till z~1. However, when the stellar content is large enough so that global dynamical instabilities can be triggered, the galaxy in Eris2k develops a stronger and more extended bar with respect to ErisBH. We demonstrate that he structural properties and time evolution of the two bars are very different. Our results highlight the importance of accurate sub-grid prescriptions in cosmological zoom-in simulations of the process of galaxy formation and evolution, and the possible use of a statistical sample of barred galaxies to assess the strength of the stellar feedback., Comment: 16 pages, 14 Figure; accepted for publication in MNRAS; added references and figures, more details, results unchanged

In a recent paper (Fumagalli et al. 2017) we reported on the detection of a diffuse H$\alpha$ glow in the outskirts of the nearby, edge-on disc galaxy UGC 7321 observed with the Multi Unit Spectroscopic Explorer (MUSE) at the ESO Very Large Telescope. By interpreting the H$\alpha$ emission as fluorescence arising from hydrogen ionised by an external (i.e., extragalactic) radiation field, we estimated the UV background (UVB) intensity in terms of HI ionisation rate (per ion) at $z\simeq 0$ to be in the range $\Gamma_{\rm HI}\sim 6-8\times 10^{-14}$ s$^{-1}$. In the present work, by performing radiative transfer calculations over a large set of models of the gaseous disc of UGC 7321, we refine our estimate and through an MCMC analysis derive a value for the photoionisaton rate of $\Gamma_{\rm HI}=7.27^{+2.93}_{-2.90}\times 10^{-14}$ s$^{-1}$. In particular, our analysis demonstrates that this value is robust against large variations in the galaxy model and that the uncertainties are mainly driven by the errors associated with the observed H$\alpha$ surface brightness. Our measurement is consistent with several recent determinations of the same quantity by a completely independent technique (i.e., flux decrement analysis of the Ly$\alpha$-forest), and support the notion that the low redshift UVB is largely dominated by active galactic nuclei (AGNs), possibly with no need of further contribution from star forming galaxies., Comment: 6 pages, 2 figures, submitted to MNRAS

The grand challenges of contemporary fundamental physics---dark matter, dark energy, vacuum energy, inflation and early universe cosmology, singularities and the hierarchy problem---all involve gravity as a key component. And of all gravitational phenomena, black holes stand out in their elegant simplicity, while harbouring some of the most remarkable predictions of General Relativity: event horizons, singularities and ergoregions. The hitherto invisible landscape of the gravitational Universe is being unveiled before our eyes: the historical direct detection of gravitational waves by the LIGO-Virgo collaboration marks the dawn of a new era of scientific exploration. Gravitational-wave astronomy will allow us to test models of black hole formation, growth and evolution, as well as models of gravitational-wave generation and propagation. It will provide evidence for event horizons and ergoregions, test the theory of General Relativity itself, and may reveal the existence of new fundamental fields. The synthesis of these results has the potential to radically reshape our understanding of the cosmos and of the laws of Nature. The purpose of this work is to present a concise, yet comprehensive overview of the state of the art in the relevant fields of research, summarize important open problems, and lay out a roadmap for future progress., Comment: White Paper for the COST action "Gravitational Waves, Black Holes, and Fundamental Physics", 272 pages, 12 figures; v4: updated references and author list. Overall improvements and corrections. To appear in Classical and Quantum Gravity

It has been proposed that close interactions with satellite galaxies can significantly perturb the morphology of the main galaxy. However, the dynamics of an already formed bar following the interaction with the external environment has not been studied in detail in a fully cosmological context. In this work, analysing the cosmological zoom-in simulation Eris2k, we study the effects that a very unequal-mass flyby crossing the stellar disc has on the stability of the pre-existing bar. We characterize the evolution of the bar strength and length showing that the perturbation exerted by the flyby shuffles the orbits of stars for less than one Gyr. After this time, the bar shows a remarkable resilience, reforming with properties comparable to those it had before the interaction. Our work shows that close unequal-mass encounters, the most frequent interactions occurring during the evolution of cosmic structures, have (i) an overall minor impact on the global evolution of the bar in the long term, still (ii) the effect is destructive and (iii) a very weak interaction is sufficient to dismantle a strong bar leading to its "apparent death". As a consequence, due to the non-negligible duration of the bar-less period, a fraction of observed spiral galaxies classified as non-barred could be prone to bar formation., Comment: 6 pages, 3 Figure; Accepted for publication in MNRAS

Recent observations suggest that hydrogen reionization ends late ($z \simeq 6$) and proceeds quickly. We present here a new model of the meta-galactic UV/X-ray background (UVB) that is consistent with this. It adopts the most recent determinations of the ionizing emissivity due to stars and AGN, as well as of the HI absorber column density distribution. Another major improvement is a new treatment of the intergalactic medium (IGM) opacity for ionizing photons that is able to consistently capture the transition from a neutral to an ionized IGM. Previous synthesis models of the UVB, when used in simulations, yield reionization and thermal histories that are inconsistent with the assumed ionizing emissivities. With our new treatment, this discrepancy is fully resolved. In our fiducial model, galaxies leaking $\lesssim 18\%$ of their Lyman continuum emission drive HI reionization, while AGN drive HeII reionization (completing at $z \simeq 6.2$ and $2.8$ respectively). Due to the limited time available for cooling between HI and HeII reionization, higher IGM temperatures are predicted for late reionization scenarios. In our fiducial model, the predicted temperatures agree well with observational constraints at $z \lesssim 4$, while being slightly high compared to (somewhat uncertain) data above that. Models with a larger contribution of AGN are instead disfavoured by the temperature data, as well as by measurements of the HI and HeII Lyman-$\alpha$ forest opacities. We also present "equivalent-equilibrium" ionization/heating rates that mimic our fiducial UVB model for use in simulation codes that assume ionization equilibrium., Comment: 23 pages, 11 figures, submitted to MNRAS, see ancillary files for machine-readable rate tables

Inspiraling massive black-hole binaries (MBHBs) forming in the aftermath of galaxy mergers are expected to be the loudest gravitational-wave (GW) sources relevant for pulsar-timing arrays (PTAs) at nHz frequencies. The incoherent overlap of signals from a cosmic population of MBHBs gives rise to a stochastic GW background (GWB) with characteristic strain around $h_c\sim10^{-15}$ at a reference frequency of 1 yr$^{-1}$, although uncertainties around this value are large. Current PTAs are piercing into the GW amplitude range predicted by MBHB-population models, but no detection has been reported so far. To assess the future success prospects of PTA experiments, it is therefore important to estimate the minimum GWB level consistent with our current understanding of the formation and evolution of galaxies and massive black holes (MBHs). To this purpose, we couple a semianalytic model of galaxy evolution and an extensive study of the statistical outcome of triple MBH interactions. We show that even in the most pessimistic scenario where all MBHBs stall before entering the GW-dominated regime, triple interactions resulting from subsequent galaxy mergers inevitably drive a considerable fraction of the MBHB population to coalescence. At frequencies relevant for PTA, the resulting GWB is only a factor of 2-to-3 suppressed compared to a fiducial model where binaries are allowed to merge over Gyr timescales. Coupled with current estimates of the expected GWB amplitude range, our findings suggest that the minimum GWB from cosmic MBHBs is unlikely to be lower than $h_c\sim 10^{-16}$ (at $f = 1$ yr$^{-1}$), well within the expected sensitivity of projected PTAs based on future observations with FAST, MeerKAT and SKA., Comment: 16 pages, 9 figures, accepted for publication in MNRAS

Massive black hole binaries (MBHBs) are expected to form at the centre of merging galaxies during the hierarchical assembly of the cosmic structure, and are expected to be the loudest sources of gravitational waves (GWs) in the low frequency domain. However, because of the dearth of energy exchanges with background stars and gas, many of these MBHBs may stall at separations too large for GW emission to drive them to coalescence in less than a Hubble time. Triple MBH systems are then bound to form after a further galaxy merger, triggering a complex and rich dynamics that can eventually lead to MBH coalescence. Here we report on the results of a large set of numerical simulations, where MBH triplets are set in spherical stellar potentials and MBH dynamics is followed through 2.5 post-Newtonian order in the equations of motion. From our full suite of simulated systems we find that a fraction $\simeq 20-30$% of the MBH binaries that would otherwise stall are led to coalesce within a Hubble time. The corresponding coalescence timescale peaks around 300 Myr, while the eccentricity close to the plunge, albeit small, is non-negligible ($\lesssim 0.1$). We construct and discuss marginalised probability distributions of the main parameters involved and, in a companion paper of the series, we will use the results presented here to forecast the contribution of MBH triplets to the GW signal in the nHz regime probed by Pulsar Timing Array experiments., Comment: 19 pages, 14 figures, accepted for publication in MNRAS

We highlight some subtleties that affect naive implementations of quadrupolar and octupolar gravitational waveforms from numerically-integrated trajectories of three-body systems. Some of those subtleties arise from the requirement that the source be contained in its "coordinate near zone" when applying the standard PN formulae for gravitational-wave emission, and from the need to use the non-linear Einstein equations to correctly derive the quadrupole emission formula. We show that some of these subtleties were occasionally overlooked in the literature, with consequences for published results. We also provide prescriptions that lead to correct and robust predictions for the waveforms computed from numerically-integrated orbits., Comment: 20 pages, 6 figures

Recent work suggests that strong emission line, star-forming galaxies may be significant Lyman Continuum leakers. We combine archival HST broadband ultraviolet and optical imaging (F275W and F606W, respectively) with emission line catalogs derived from WFC3 IR G141 grism spectroscopy to search for escaping Lyman Continuum (LyC) emission from homogeneously selected $z\simeq$2.5 SFGs. We detect no escaping Lyman Continuum from SFGs selected on [OII] nebular emission (N=208) and, within a narrow redshift range, on [OIII]/[OII]. We measure 1$\sigma$ upper limits to the LyC escape fraction relative to the non-ionizing UV continuum from [OII] emitters, $f_{esc}<$5.6%, and strong [OIII]/[OII]$>$5 ELGs, $f_{esc}<$14.0%. Our observations are not deep enough to detect $f_{esc}\lesssim$10% typical of the low redshift Lyman continuum emitters. However, we find that this population represents a small fraction of the star-forming galaxy population at $z\simeq$2. Thus, unless the number of extreme emission line galaxies grows substantially to z$>$6, such galaxies may be insufficient for reionization. Deeper survey data in the rest-frame ionizing UV will be necessary to determine whether strong line ratios could be useful for pre-selecting LyC leakers at high redshift., Comment: accepted to ApJ Letters

The emergence of a large-scale stellar bar is one of the most striking features in disc galaxies. By means of state-of-the-art cosmological zoom-in simulations, we study the formation and evolution of bars in Milky Way-like galaxies in a fully cosmological context, including the physics of gas dissipation, star formation, and supernova feedback. Our goal is to characterise the actual trigger of the non-axisymmetric perturbation that leads to the strong bar observable in the simulations at z=0, discriminating between an internal/secular versus an external/tidal origin. To this aim, we run a suite of cosmological zoom-in simulations altering the original history of galaxy-satellite interactions at a time when the main galaxy, though already bar-unstable, does not feature any non-axisymmetric structure yet. We find that the main effect of a late minor merger and of a close fly-by is to delay the time of bar formation and those two dynamical events are not directly responsible for the development of the bar and do not alter significantly its global properties (e.g. its final extension). We conclude that, once the disc has grown to a mass large enough to sustain global non-axisymmetric modes, then bar formation is inevitable., Comment: 15 pages, 6 Figure; Accepted for publication in MNRAS - Full-resolution images are available at http://www.dfm.uninsubria.it/fh/FHpages/Research.html

In an effort to understand the cause of the apparent depletion in the number density of radio-loud AGNs at $z>3$, this work investigates the viability of the so-called Cosmic Microwave Background (CMB) quenching mechanism of intrinsically jetted, high-z AGNs, whereby Inverse Compton scattering of CMB photons off electrons within the extended lobes results in a substantial dimming of the lobe synchrotron emission at GHz frequencies, while simultaneously boosting their diffuse X-ray signal. We focus on five $z>3.5$ radio galaxies that have sufficiently deep Chandra exposure (> 50 ks) to warrant a meaningful investigation of any extended X-ray emission. For those objects with evidence for statistically significant extended X-ray lobes (4C 41.17 and 4C 03.24), we combine the Chandra measurements with literature data at lower frequencies to assemble the systems' Spectral Energy Distributions (SEDs), and utilize state-of-the-art SED modelling (Ghisellini et al. 2015) -- including emission from the disk, torus, jet, hotspots, and lobes -- to infer their physical parameters. For both radio galaxies, the magnetic energy density in the hotspots is found to exceed the energy density in CMB photons, wheres the opposite is true for the lobes. This implies that any extended synchrotron emission likely originates from the hotspots themselves, rather than the lobes. Conversely, Inverse Compton scattering of CMB photons dominates the extended X-ray emission from the lobes, which are effectively "radio-quenched". As a result, CMB quenching is effective in these systems in spite of the fact that the observed X-ray to radio luminosity ratio does not bear the signature $(1+z)^4$ dependence of the CMB energy density., Comment: MNRAS in press; 17 pages, 8 figures, 4 tables

We report the detection of extended Halpha emission from the tip of the HI disk of the nearby edge-on galaxy UGC 7321, observed with the Multi Unit Spectroscopic Explorer (MUSE) instrument at the Very Large Telescope. The Halpha surface brightness fades rapidly where the HI column density drops below N(HI) = 10^19 cm^-2 , consistent with fluorescence arising at the ionisation front from gas that is photoionized by the extragalactic ultraviolet background (UVB). The surface brightness measured at this location is (1.2 +/- 0.5)x10^-19 erg/s/cm^2/arcsec^2, where the error is mostly systematic and results from the proximity of the signal to the edge of the MUSE field of view, and from the presence of a sky line next to the redshifted Halpha wavelength. By combining the Halpha and the HI 21 cm maps with a radiative transfer calculation of an exponential disk illuminated by the UVB, we derive a value for the HI photoionization rate of Gamma ~ (6-8)x10^-14 1/s . This value is consistent with transmission statistics of the Lyalpha forest and with recent models of a UVB which is dominated by quasars., Comment: 15 pages, 11 figures. Published in MNRAS