Astrophysics of Galaxies (astro-ph.GA)

Abstract: The inner kiloparsec regions surrounding sub-Eddington (luminosity less than 10$^{-3}$ in Eddington units, L$_{Edd}$) supermassive black holes (BHs) often show a characteristic network of dust filaments that terminate in a nuclear spiral in the central parsecs. Here we study the role and fate of these filaments in one of the least accreting BHs known, M31 (10$^{-7}$ L$_{Edd}$) using hydrodynamical simulations. The evolution of a streamer of gas particles moving under the barred potential of M31 is followed from kiloparsec distance to the central parsecs. After an exploratory study of initial conditions, a compelling fit to the observed dust/ionized gas morphologies and line-of-sight velocities in the inner hundreds of parsecs is produced. After several million years of streamer evolution, during which friction, thermal dissipation, and self-collisions have taken place, the gas settles into a disk tens of parsecs wide. This is fed by numerous filaments that arise from an outer circumnuclear ring and spiral toward the center. The final configuration is tightly constrained by a critical input mass in the streamer of several 10$^3$ M$_{\odot}$ (at an injection rate of 10$^{-4}$ M$_{\odot}$ yr$^{-1}$); values above or below this lead to filament fragmentation or dispersion respectively, which are not observed. The creation of a hot gas atmosphere in the region of $\sim$10$^6$ K is key to the development of a nuclear spiral during the simulation. The final inflow rate at 1pc from the center is $\sim$1.7 $\times$ 10$^{-7}$ M$_{\odot}$ yr$^{-1}$, consistent with the quiescent state of the M31 BH.

Abstract: Asymmetric emission of gravitational waves during mergers of black holes (BHs) produces a recoil kick, which can set a newly formed BH on a bound orbit around the center of its host galaxy, or even completely eject it. To study this population of recoiling BHs we extract properties of galaxies with merging BHs from Illustris TNG300 simulation and then employ both analytical and numerical techniques to model unresolved process of BH recoil. This comparative analysis between analytical and numerical models shows that, on cosmological scales, numerically modeled recoiling BHs have a higher escape probability and predict a greater number of offset active galactic nuclei (AGN). BH escaped probability $>40~ \%$ is expected in 25 $\%$ of merger remnants in numerical models, compared to 8$\%$ in analytical models. At the same time, the predicted number of offset AGN at separations $>5$ kpc changes from 58 $\%$ for numerical models to 3 $\%$ for analytical models. Since BH ejections in major merger remnants occur in non-virialized systems, static analytical models cannot provide an accurate description. Thus we argue that numerical models should be used to estimate the expected number density of escaped BHs and offset AGN.

Abstract: Among the diverse progenitor channels leading to Type Ia Supernovae (SNe Ia), there are explosions originating from white dwarfs with sub-Chandrasekhar masses. These white dwarfs undergo detonation and explosion triggered by primary detonation in the helium shell, which has been accreted from a companion star. The double-detonation model predicts a correlation between the age of the progenitor system and the near peak brightness: the younger the exploding progenitors, the brighter the SNe. In this paper, we present our recent achievements on the study of SNe Ia properties in different locations within host galactic discs and the estimation of their progenitor population ages. Observationally, we confirm the validity of the anticipated correlation between the SN photometry and the age of their progenitors.

Abstract: We present a comprehensive search and analysis of high-redshift galaxies in a suite of nine public JWST extragalactic fields taken in Cycle 1, covering a total effective search area of $\sim358{\rm arcmin^2}$. Through conservative ($8\sigma$) photometric selection, we identify 339 galaxies at $5<z<14$, with 109 having spectroscopic redshift measurements from the literature, including recent JWST NIRSpec observations. Our regression analysis reveals that the rest-frame UV size-stellar mass relation follows $R_{\rm eff}\propto M_*^{0.20\pm0.03}$, similar to that of star-forming galaxies at $z\sim3$, but scaled down in size by $\sim0.7$dex. We find a much slower rate for the average size evolution over the redshift range, $R_{\rm eff}\propto(1+z)^{-0.4\pm0.2}$, than that derived in the literature. A fraction ($\sim13\,\%$) of our sample are marginally resolved even in the NIRCam imaging ($<100$pc), located at $>1.5\,\sigma$ below the derived size-mass slope. These compact sources exhibit a high star formation surface density $\Sigma_{\rm SFR}>10\,M_\odot\,{\rm yr^{-1}\,kpc^{-2}}$, a range in which only $<0.01\,\%$ of the local star-forming galaxy sample is found. For those with available NIRSpec data, no evidence of ongoing supermassive black hole accretion is observed. A potential explanation for the observed high [OIII]-to-Hbeta ratios could be high shock velocities, likely originating within intense star-forming regions characterized by high $\Sigma_{\rm SFR}$. Lastly, we find that the rest-frame UV and optical sizes of our sample are comparable. Our results are consistent with these early galaxies building up their structures inside-out and yet to exhibit the strong color gradient seen at lower redshift.

Abstract: Detections of extended Ly$\alpha$ halos (LAHs) around Ly$\alpha$ emitters (LAEs) have lately been reported on a regular basis, but their origin is still under investigation. Simulation studies predict that the outer regions of the extended LAHs contain a major contribution from the Ly$\alpha$ emission of faint, individually undetected LAEs. To address this matter from an observational angle, we use halo occupation distribution (HOD) modeling to reproduce the clustering of a spectroscopic sample of 1265 LAEs at $3<z<5$ from the MUSE-Wide survey. We integrate the Ly$\alpha$ luminosity function (LF) to estimate the background surface brightness due to discrete faint LAEs. We then extend the HOD statistics inwards towards small separations and compute the factor by which the measured Ly$\alpha$ surface brightness (SB) is enhanced by undetected close physical neighbors. We consider various clustering scenarios for the undetected sources and compare the corresponding radial profiles. The resulting inferred Ly$\alpha$ SB of faint LAEs ranges between $(0.4-2)\times10^{20}\;\rm{erg}\;\rm{s}^{-1}~\rm{cm}^{-2}~\rm{arcsec}^{-2}$, with a very slow radial decline outwards. Our results suggest that the outer regions of observed LAHs ($R\gtrsim50$~pkpc) could indeed contain a strong component from external (but physically associated) LAEs, possibly even be dominated by them. Only for a relatively shallow faint-end slope of the Ly$\alpha$ LF would this contribution from clustered LAEs become unimportant. We also confirm that the observed emission from the inner regions ($R\le20-30$~pkpc) is too bright to be significantly affected by clustering. We compare our findings with predicted profiles from simulations and find good overall agreement. We outline possible future measurements to further constrain the impact of discrete undetected LAEs on observed extended LAHs.

Abstract: The mass assembly history (MAH) of dark matter halos plays a crucial role in shaping the formation and evolution of galaxies. MAHs are used extensively in semi-analytic and empirical models of galaxy formation, yet current analytic methods to generate them are inaccurate and unable to capture their relationship with the halo internal structure and large-scale environment. This paper introduces FLORAH, a machine-learning framework for generating assembly histories of ensembles of dark matter halos. We train FLORAH on the assembly histories from the GUREFT and VSMDPL N-body simulations and demonstrate its ability to recover key properties such as the time evolution of mass and concentration. We obtain similar results for the galaxy stellar mass versus halo mass relation and its residuals when we run the Santa Cruz semi-analytic model on FLORAH-generated assembly histories and halo formation histories extracted from an N-body simulation. We further show that FLORAH also reproduces the dependence of clustering on properties other than mass (assembly bias), which is not captured by other analytic methods. By combining multiple networks trained on a suite of simulations with different redshift ranges and mass resolutions, we are able to construct accurate main progenitor branches (MPBs) with a wide dynamic mass range from $z = 0$ up to an ultra-high redshift $z \approx 20$, currently far beyond that of a single N-body simulation. FLORAH is the first step towards a machine learning-based framework for planting full merger trees; this will enable the exploration of different galaxy formation scenarios with great computational efficiency at unprecedented accuracy.

Abstract: We present a study of the molecular gas of seven early-type galaxies with high angular resolution data obtained as part of the mm-Wave Interferometric Survey of Dark Object Masses (WISDOM) project with the Atacama Large Millimeter/submillimeter Array. Using a fixed spatial scale approach, we study the mass surface density ($\Sigma$) and velocity dispersion ($\sigma$) of the molecular gas on spatial scales ranging from $60$ to $120$pc. Given the spatial resolution of our data ($20$ - $70$pc), we characterise these properties across many thousands of individual sight lines ($\approx50,000$ at our highest physical resolution). The molecular gas along these sight lines has a large range ($\approx2$dex) of mass surface densities and velocity dispersions $\approx40\%$ higher than those of star-forming spiral galaxies. It has virial parameters $\alpha_\mathrm{vir}$ that depend weakly on the physical scale observed, likely due to beam smearing of the bulk galactic rotation, and is generally super-virial. Comparing the internal turbulent pressure ($P_\mathrm{turb}$) to the pressure required for dynamic equilibrium ($P_\mathrm{DE}$), the ratio $P_\mathrm{turb}$/$P_\mathrm{DE}$ is significantly less than unity in all galaxies, indicating that the gas is not in dynamic equilibrium and is strongly compressed, in apparent contradiction to the virial parameters. This may be due to our neglect of shear and tidal forces, and/or the combination of three-dimensional and vertical diagnostics. Both $\alpha_\mathrm{vir}$ and $P_\mathrm{turb}$ anti-correlate with the global star-formation rate of our galaxies. We therefore conclude that the molecular gas in early-type galaxies is likely unbound, and that large-scale dynamics likely plays a critical role in its regulation. This contrasts to the giant molecular clouds in the discs of late-type galaxies, that are much closer to dynamical equilibrium.

Abstract: I report the discovery of a stellar stream (Sutlej) using Gaia DR3 proper motions and XP metallicities located ~15 degrees north of the Small Magellanic Cloud (SMC). The stream is composed of two parallel linear components ("branches") approximately ~8 x 0.6 degrees in size and separated by 2.5 degrees. The stars have a mean proper motion of (pmra,pmdec)=(+0.08 mas/yr,-1.41 mas/yr) which is quite similar to the proper motion of stars on the western side of the SMC. The color magnitude diagram of the stream stars has a clear red giant branch, horizontal branch, and main sequence turnoff that is well-matched by a PARSEC isochrone of 10 Gyr, [Fe/H]=-1.8 at 32 kpc and a total stellar mass of ~33,000 Msun. The stream is spread out over an area of 9.6 square degrees and has a surface brightness of 32.5 mag/arcsec^2. The metallicity of the stream stars from Gaia XP spectra extend over -2.5 < [M/H] < -1.0 with a median of [M/H]=-1.8. The tangential velocity of the stream stars is 214 km/s compared to the values of 448 km/s for the Large Magellanic Cloud and 428 km/s for the SMC. While the radial velocity of the stream is not yet known, a comparison of the space velocities using a range of assumed radial velocities, shows that the stream is unlikely to be associated with the Magellanic Clouds. The tangential velocity vector is misaligned with the stream by ~25 degrees which might indicate an important gravitational influence from the nearby Magellanic Clouds.

Abstract: The quasar HE0230$-$2130 is lensed by two galaxies at similar redshifts into four observed images. Using modeled quasar positions from fitting the brightness of the quasar images in ground-based imaging data from the Magellan telescope, we find that lens mass models where both galaxies are each parametrized with a singular power-law (PL) profile predict five quasar images. One of the predicted images is unobserved even though it is distinctively offset from the lensing galaxies and is bright enough to be observable. This missing image gives rise to new opportunities to study the galaxies' mass distribution. To interpret the quad configuration of this system, we test different profile assumption with the aim to obtain lens mass models that predicts correctly only four observed images. We test the effect of adopting cored profiles for the lensing galaxies, of external shear, and of additional profiles to represent a dark matter clump. By comparing the Bayesian evidence of different model parametrizations, we favor the model class that consists of two singular PL profiles for the lensing galaxies and a cored isothermal sphere in the region of the previously predicted fifth images (rNIS profile). We estimate the mass of the rNIS clumps inside its Einstein radius and find that 18\% are in the range $10^6 M_{\odot} \leq M_{\rm rNIS}\leq 10^9 M_{\odot}$, which is the predicted mass range of dark matter subhalos in cold dark matter simulations, or the mass of low-mass dark matter satellite galaxies. The second most likely model class, with a relative probability of 94\%, is the model where the smaller lensing galaxy is described by a cored PL profile with external shear. Our study demonstrates that lensed quasar images are sensitive to dark matter structure in the gravitational lens.

Abstract: Understanding the spectral properties of sources is crucial for the characterization of the radio source population. In this work, we have extensively studied the ELAIS N1 field using various low-frequency radio observations. For the first time, we present the 1250\,MHz observations of the field using the upgraded Giant Meterwave Radio Telescope (uGMRT) that reach a central off-source RMS noise of $\sim 12\,\mu$Jy\,beam$^{-1}$. A source catalogue of 1086 sources is compiled at $5\sigma$ threshold ($>60\,\mu$Jy) to derive the normalized differential source counts at this frequency that is consistent with existing observations and simulations. We present the spectral indices derived in two ways: two-point spectral indices and by fitting a power-law. The latter yielded a median $\alpha = -0.57\pm 0.14$, and we identified nine ultra-steep spectrum sources using these spectral indices. Further, using a radio colour diagram, we identify the three mega-hertz peaked spectrum (MPS) sources, while three other MPS sources are identified from the visual inspection of the spectra, the properties of which are discussed. In our study of the classified sources in the ELAIS N1 field, we present the relationship between $\alpha$ and $z$. We find no evidence of an inverse correlation between these two quantities and suggest that the nature of the radio spectrum remains independent of the large-scale properties of the galaxies that vary with redshifts.