Young Moving Groups (YMGs) are close (<100pc), coherent collections of young (<100Myr) stars that appear to have formed in the same star-forming molecular cloud. As such we would expect their individual initial mass functions (IMFs) to be similar to other star-forming regions, and by extension the Galactic field. Their close proximity to the Sun and their young ages means that YMGs are promising locations to search for young forming exoplanets. However, due to their low numbers of stars, stochastic sampling of the IMF means their stellar populations could vary significantly. We determine the range of planet-hosting stars (spectral types A, G and M) possible from sampling the IMF multiple times, and find that some YMGs appear deficient in M-dwarfs. We then use these data to show that the expected probability of detecting terrestrial magma ocean planets is highly dependent on the exact numbers of stars produced through stochastic sampling of the IMF.

Studies investigating the relationship between AGN power and the star formation rates (SFRs) of their host galaxies often rely on averaging techniques -- such as stacking -- to incorporate information from non-detections. However, averages, and especially means, can be strongly affected by outliers and can therefore give a misleading indication of the "typical" case. Recently, a number of studies have taken a step further by binning their sample in terms of AGN power (approximated by the 2-10keV luminosity of the AGN), and investigating how the SFR distribution differs between these bins. These bin thresholds are often weakly motivated, and binning implicitly assumes that sources within the same bin have similar (or even identical) properties. In this paper, we investigate whether the distribution of host SFRs -- relative to the locus of the star-forming main sequence (i.e., $R_{\rm MS}$) -- changes continuously as a function of AGN power. We achieve this by using a hierarchical Bayesian model that completely removes the need to bin in AGN power. In doing so, we find strong evidence that the $R_{\rm MS}$ distribution changes with 2-10keV X-ray luminosity. The results suggest that higher X_ray luminosity AGNs have a tighter physical connection to the star-forming process than lower X-ray luminosity AGNs, at least within the $0.8 < z < 1.2$ redshift range considered here.

We present a new modeling of the X-ray luminosity function (XLF) of Active Galactic Nuclei (AGN) out to z$\sim$3, dissecting the contribution of main-sequence (MS) and starburst (SB) galaxies. For each galaxy population, we convolved the observed galaxy stellar mass (M$_{\star}$) function with a grid of M$_{\star}$-independent Eddington ratio ($\lambda_{\rm EDD}$) distributions, normalised via empirical black hole accretion rate (BHAR) to star formation rate (SFR) relations. Our simple approach yields an excellent agreement with the observed XLF since z$\sim$3. We find that the redshift evolution of the observed XLF can only be reproduced through an intrinsic flattening of the $\lambda_{\rm EDD}$ distribution, and with a positive shift of the break $\lambda^{*}$, consistent with an anti-hierarchical behavior. The AGN accretion history is predominantly made by massive (10$^{10}<$M$_{\star}<$10$^{11}$ M$_{\odot}$) MS galaxies, while SB-driven BH accretion, possibly associated with galaxy mergers, becomes dominant only in bright quasars, at $\log$(L$_{\rm X}$/erg s$^{-1}$)$>$44.36 + 1.28$\cdot$(1+z). We infer that the probability of finding highly-accreting ($\lambda_{\rm EDD}>$ 10%) AGN significantly increases with redshift, from 0.4% (3.0%) at z=0.5 to 6.5% (15.3%) at z=3 for MS (SB) galaxies, implying a longer AGN duty cycle in the early Universe. Our results strongly favor a M$_{\star}$-dependent ratio between BHAR and SFR, as BHAR/SFR $\propto$ M$_{\star}^{0.73[+0.22,-0.29]}$, supporting a non-linear BH buildup relative to the host. Finally, this framework opens potential questions on super-Eddington BH accretion and different $\lambda_{\rm EDD}$ prescriptions for understanding the cosmic BH mass assembly.

Supermassive black holes (SMBHs) are tightly correlated with their hosts but the origin of such connection remains elusive. To explore the cosmic build-up of this scaling relation, we present an empirically-motivated model that tracks galaxy and SMBH growth down to z=0. Starting from a random mass seed distribution at z=10, we assume that each galaxy evolves on the star-forming "main sequence" (MS) and each BH follows the recently-derived stellar mass (M$_{\star}$) dependent ratio between BH accretion rate and star formation rate, going as BHAR/SFR$\propto$M$_{\star}^{0.73[+0.22,-0.29]}$. Our simple recipe naturally describes the BH-galaxy build-up in two stages. At first, the SMBH lags behind the host that evolves along the MS. Later, as the galaxy grows in M$_{\star}$, our M$_{\star}$-dependent BHAR/SFR induces a super-linear BH growth, as M$_{\rm BH}$$\propto$M$_{\star}^{1.7}$. According to this formalism, smaller BH seeds increase their relative mass faster and earlier than bigger BH seeds, at fixed M$_{\star}$, thus setting along a gradually tighter M$_{\rm BH}$-M$_{\star}$ locus towards higher M$_{\star}$. Assuming reasonable values of the radiative efficiency $\epsilon \sim$0.1, our empirical trend agrees with both high-redshift model predictions and intrinsic M$_{\rm BH}$-M$_{\star}$ relations of local BHs. We speculate that the observed non-linear BH-galaxy build-up is reflected in a twofold behavior with dark matter halo mass (M$_{\rm DM}$), displaying a clear turnover at M$_{\rm DM}\sim$2$\times$10$^{12}$M$_{\odot}$. While Supernovae-driven feedback suppresses BH growth in smaller halos (BHAR/SFR$\propto$M$_{\rm DM}^{1.6}$), above the M$_{\rm DM}$ threshold cold gas inflows possibly fuel both BH accretion and star formation in a similar fashion (BHAR/SFR$\propto$M$_{\rm DM}^{0.3}$).

We use $N$-body simulations to probe the early phases of the dynamical evolution of star-forming regions and focus on mass and velocity distributions of unbound stars. In this parameter space study, we vary the initial virial ratio and degree of spatial and kinematic substructure and analyse the fraction of stars that become unbound in two different mass classes (above and below 8 M$_{\odot}$). We find that the fraction of unbound stars differs depending on the initial conditions. After 10 Myr, in initially highly subvirial, substructured simulations, the high-mass and lower-mass unbound fractions are similar at $\sim$23 per cent. In initially virialised, substructured simulations, we find only $\sim$16 per cent of all high-mass stars are unbound, whereas $\sim$37 per cent of all lower-mass stars are. The velocity distributions of unbound stars only show differences for extremely different initial conditions. The distributions are dominated by large numbers of lower-mass stars becoming unbound just above the escape velocity of $\sim$3 km s$^{-1}$ with unbound high-mass stars moving faster on average than lower-mass unbound stars. We see no high-mass runaway stars (velocity > 30 km s$^{-1}$) from any of our initial conditions and only an occasional lower-mass runaway star from initially subvirial/substructured simulations. In our simulations, we find a small number of lower-mass walkaway stars (with velocity 5-30 km s$^{-1}$) from all of our initial conditions. These walkaway stars should be observable around many nearby star-forming regions with Gaia.

We infer and compare the specific X-ray luminosity distributions for a sample of massive (i.e. $\log_{10} (M*/M\odot) > 10.5$) galaxies split according to their far-infrared-derived star-forming properties (i.e., starburst and non-starburst) and redshift. We model each distribution as a power-law with an upper and lower turnover, and adopt a maximum likelihood method to include information from non-detections in the form of upper limits. When we use our inferred distributions to calculate the ratios of high to low sLx AGN (corresponding to above and below $0.1\lambda_{\text{Edd}}$, respectively) we find that starbursts have significantly higher proportions of high sLx AGN compared to their non-starburst counterparts. These findings help explain the increase in average X-ray luminosity in bins of increasing SFR reported by previous studies.

We explore the distribution of RMS=SFR/SFR_MS (where SFR_MS is the star formation rate of "Main Sequence" star-forming galaxies) for AGN hosts at z=1. We split our sample into two bins of X-ray luminosity divided at Lx=2x10^43erg s-1 to investigate whether the RMS distribution changes as a function of AGN power. Our main results suggest that, when the RMS distribution of AGN hosts is modelled as a log-normal distribution (i.e. the same shape as that of MS galaxies), galaxies hosting more powerful X-ray AGNs (i.e. Lx>2x10^43erg s-1) display a narrower RMS distribution that is shifted to higher values compared to their lower Lx counterparts. In addition, we find that more powerful X-ray AGNs have SFRs that are more consistent with that of MS galaxies compared to lower Lx AGNs. Despite this, the mean SFRs (as opposed to RMS) measured from these distributions are consistent with the previously observed flat relationship between SFR and Lx. Our results suggest that the typical star-forming properties of AGN hosts change with Lx , and that more powerful AGNs typically reside in more MS-like star-forming galaxies compared to lower Lx AGNs.

The lack of a strong correlation between AGN X-ray luminosity ($L_X$; a proxy for AGN power) and the star formation rate (SFR) of their host galaxies has recently been attributed to stochastic AGN variability. Studies using population synthesis models have incorporated this by assuming a broad, universal (i.e. does not depend on the host galaxy properties) probability distribution for AGN specific X-ray luminosities (i.e. the ratio of $L_X$ to host stellar mass; a common proxy for Eddington ratio). However, recent studies have demonstrated that this universal Eddington ratio distribution fails to reproduce the observed X-ray luminosity functions beyond z$\sim$1.2. Furthermore, empirical studies have recently shown that the Eddington ratio distribution may instead depend upon host galaxy properties, such as SFR and/or stellar mass. To investigate this further we develop a population synthesis model in which the Eddington ratio distribution is different for star-forming and quiescent host galaxies. We show that, although this model is able to reproduce the observed X-ray luminosity functions out to z$\sim$2, it fails to simultaneously reproduce the observed flat relationship between SFR and X-ray luminosity. We can solve this, however, by incorporating a mass dependency in the AGN Eddington ratio distribution for star-forming host galaxies. Overall, our models indicate that a relative suppression of low Eddington ratios ($\lambda_{\rm Edd}\lesssim$0.1) in lower mass galaxies (M<$10^{10-11}$Msun) is required to reproduce both the observed X-ray luminosity functions and the observed flat SFR/X-ray relationship.