Astrophysics of Galaxies (astro-ph.GA)

Abstract: We use a complete set of deep narrow-band imaging data for 384 galaxies gathered during the VESTIGE survey to derive the first Halpha luminosity function (LF) of the Virgo cluster within R200. The data allow us to cover the whole dynamic range of the Halpha LF (10^36<LHa<10^42 erg s^-1). After they are corrected for [NII] contamination and dust attenuation, the data are used to derive the SFR function in the range 10^-4<SFR<10 Mo yr^-1. These LF are compared to those derived at other frequencies or using different tracers of star formation in Virgo, in other nearby and high-z clusters, in the field, and to those predicted by the IllustrisTNG cosmological hydrodynamical simulations. The Halpha LF of the Virgo cluster is fairly flat (a=-1.07) in the range 10^38.5<LHa<10^40.5 erg s^-1, and it abruptly decreases at lower luminosities. When compared to those derived for other nearby clusters and for the field, the slope and the characteristic luminosity of the Schechter function change as a function of the dynamical mass of the system, of the temperature of the X-rays gas, and of the dynamical pressure exerted on the interstellar medium of galaxies moving at high velocity within the intracluster medium. All these trends can be explained in a scenario in which the activity of SF is reduced in massive clusters due to their hydrodynamical interaction with the surrounding medium, suggesting once again that ram-pressure stripping is the dominant mechanism affecting galaxy evolution in local clusters of dynamical mass M200>10^14 Mo. The comparison with the IllustrisTNG cosmological hydrodynamical simulations shows a more pronounced decrease at the faint end of the distribution. If Virgo is representative of typical nearby clusters of similar mass, this difference suggests that the stripping process in simulated galaxies in these environments is more efficient than observed.

Abstract: What is the nature of a star forming clump? Observations reveal these to be chaotic environments being modified and influenced by many physical processes. However, numerical simulations often define these initial star forming clumps to be idealised objects. In this paper, we define and analyse 109 star forming clumps extracted from our previous low-mass star cluster simulations. To define a clump, we identify all the gas in a simulation that ever becomes bound to or accreted onto a star, then follow the gas backwards in time until it decreases to a critical density. This gas, and its neighbouring gas, is defined as our star forming clump. Our clumps span a mass range of $0.15 \lesssim M/$M$_\odot \lesssim 10.2$, while the density range within each clump spans 2--4 orders of magnitude. The gas density distribution is not smooth, indicating that it is highly structured. The clumps are turbulent, with no coherent rotation. Independent of the initial magnetic field strength of the parent cloud, all clumps yield a similar range of field strengths. The clump magnetic field is ordered, but not reflective of the initial field geometry of the parent cloud. In general, most clump properties have a slight trend with clump mass but are independent of (or only very weakly dependent on) the properties of the parent cloud. We conclude that stars are born from a wide variety of environments and there is not a single universal star forming clump.

Abstract: Stars evolving around a supermassive black hole see their orbital orientations diffuse efficiently, a process called "vector resonant relaxation". In particular, stars within the same disc, i.e. neighbors in orientations, will slowly diffuse away from one another through this stochastic process. We use jointly (i) detailed kinetic predictions for the efficiency of this dilution and (ii) the recent observation of a stellar disc around SgrA*, the supermassive black hole at the centre of the Milky-Way, to constrain SgrA*'s unobserved stellar cluster. Notably, we investigate quantitatively the impact of a population of intermediate mass black holes on the survivability of the stellar disc.

Abstract: Carbon ($^3$P) atom is a reactive species that, according to laboratory experiments and theoretical calculations, condensates with interstellar ice components. This fact is of uttermost importance for the chemistry in the interstellar medium (ISM) because the condensation reaction is barrierless and the subsequent species formed are still reactive given their open-shell character. Carbon condensation on CO-rich ices forms the \ch{C=C=O} ($^3$$\Sigma$$^-$) species, which can be easily hydrogenated twice to form ketene (H$_2$CCO). Ketene is very reactive in terrestrial conditions, usually found as an intermediate hard to be isolated in chemical synthesis laboratories. These characteristics suggest that ketene can be a good candidate to form interstellar complex organic molecules (iCOMs) via a two-step process, i.e., its activation followed by a radical-radical coupling. In this work, reactions between ketene and atomic H, and the OH and NH$_2$ radicals on a CO-rich ice model have been explored by means of quantum chemical calculations complemented by kinetic calculations to evaluate if they are favourable in the ISM. Results indicate that H addition to ketene (helped by tunneling) to form the acetyl radical (CH$_3$CO) is the most preferred path, as the reactions with OH and NH$_2$ possess activation energies ($\geq$ 9kJ/mol) hard to surmount in the ISM conditions, unless external processes provide energy to the system. Thus, acetaldehyde (CH$_3$CHO) and, probably, ethanol (CH$_3$CH$_2$OH) formation via further hydrogenations are the possible unique operating synthetic routes. Moreover, from the computed relatively large binding energies of OH and NH$_2$ on CO ice, slow diffusion is expected, hampering possible radical-radical couplings with CH$_3$CO. The astrophysical implications of these findings are discussed considering the incoming James Webb Space Telescope observations.

Abstract: We analyze a sample of 3300 galaxies between redshifts z~3.5 and z~8.5 selected from JWST images in the Hubble Ultra Deep Field (HUDF) and UKIDSS Ultra Deep Survey field, including objects with stellar masses as low as ~ 10^8 Msun up to z~8. The depth and wavelength coverage of the JWST data allow us, for the first time, to derive robust stellar masses for such high-z, low stellar-mass galaxies on an individual basis. We compute the galaxy stellar mass function (GSMF), after complementing our sample with ancillary data from CANDELS to constrain the GMSF at high stellar masses (M > M*). Our results show a steepening of the low stellar-mass end slope (a) with redshift, with a = -1.61 (+/-0.05) at z~4 and a = -1.98 (+/-0.14) at z~7. We also observe an evolution of the normalization phi* from z~7 to z~4, with phi*(z~4)/phi*(z~7)= 130 (+210/-50). Our study incorporates a novel method for the estimation of the Eddington bias that takes into account its possible dependence both on stellar mass and redshift, while allowing for skewness in the error distribution. We finally compute the resulting cosmic stellar mass density and find a flatter evolution with redshift than previous studies.

Abstract: I review methods and techniques to build mass models of disk galaxies from gas dynamics. I focus on two key steps: (1) the derivation of rotation curves using 3D emission-line datacubes from HI, CO, and/or H-alpha observations, and (2) the calculation of the gravitational field from near-infrared images and emission-line maps, tracing the stellar and gas mass distributions, respectively. Mass models of nearby galaxies led to the establishment of the radial acceleration relation (RAR): the observed centripetal acceleration from rotation curves closely correlates with that predicted from the baryonic distribution at each galaxy radius, even when dark matter supposedly dominates the gravitational field. I conclude by discussing the (uncertain) location of Local Group dwarf spheroidal galaxies on the RAR defined by more massive disk galaxies.

Abstract: We report on the discovery of cool gas inflows towards three star-forming galaxies at $\left<z\right>\sim$ 2.30. Analysis of Keck Low-Resolution Imaging Spectrometer spectroscopy reveals redshifted low-ionisation interstellar (LIS) metal absorption lines with centroid velocities of 60 - 130 km $\rm{s}^{-1}$. These inflows represent some of the most robust detections of inflowing gas into isolated, star-forming galaxies at high redshift. Our analysis suggests that the inflows are due to recycling metal-enriched gas from previous ejections. Comparisons between the galaxies with inflows and a larger parent sample of 131 objects indicate that galaxies with detected inflows may have higher specific star-formation rates (sSFR) and star-formation-rate surface densities. However, when additional galaxies without robustly detected inflows based on centroid velocity but whose LIS absorption line profiles indicate large red-wing velocities are considered, galaxies with inflows do not show unique properties relative to those lacking inflows. Additionally, we calculate the covering fraction of cool inflowing gas as a function of red-wing inflow velocity, finding an enhancement in high sSFR binned galaxies, likely due to an increase in the amount of recycling gas. Together, these results suggest that the low detection rate of galaxies with cool inflows is primarily related to the viewing angle rather than the physical properties of the galaxies.

Abstract: The gravitational wave (GW) antenna LISA will detect the signal from coalescing massive black hole binaries (MBHBs) of $\rm 10^4\,{-}\,10^7\, M_{\odot}$, providing clues on their formation and growth along cosmic history. Some of these events will be localized with a precision of several to less than a deg$^2$, enabling the possible identification of their host galaxy. This work explores the properties of the host galaxies of LISA MBHBs below $z\,{\lesssim}\,3$. We generate a simulated lightcone by using the semi-analytical model $\mathrm{\texttt{L-Galaxies}}$ applied on the merger trees of the high-resolution N-body cosmological simulation $\mathrm{\texttt{Millennium-II}}$. The model shows that LISA MBHBs are expected to be found in optically dim ($r\,{>}\,20$), star-forming ($\rm sSFR\,{>}\,10^{-10}\, \rm yr^{-1}$), gas-rich ($f_{\rm gas}\,{>}\,0.6$) and disc-dominated ($\rm B/T\,{<}\,0.7$) \textit{low-mass galaxies} of stellar masses $10^8\,{-}\,10^9 M_{\odot}$. However, these properties are indistinguishable from those of galaxies harboring single massive black holes with comparable mass, making difficult the selection of LISA hosts among the whole population of low-mass galaxies. Motivated by this, we explore the possibility of using merger signatures to select LISA hosts. We find that 40-80% of the galaxies housing LISA MBHBs display merger features related to the interaction which brought the secondary MBH to the galaxy. Despite this, around 60% of dwarf galaxies placed in the surroundings of the LISA hosts will show such kind of features as well, challenging the unequivocal detection of LISA hosts through the search of merger signatures. Consequently, the detection of an electromagnetic transient associated with the MBHB merger will be vital to pinpoint the star-forming dwarf galaxy where these binary systems evolve and coalesce.

Abstract: We have analyzed the thermal and turbulent properties of the Low-Latitude Intermediate-Velocity Arch 1 (LLIV1). This was accomplished using archival H I emission and absorption data from two 21,cm line surveys: GHIGLS at $9.^\prime$4 resolution and DHIGLS at $1^\prime$ resolution. The spectral decomposition code $\tt{ROHSA}$ was used to model the column density of different thermal phases and also to analyze an absorption measurement against the radio source 4C~+66.09. From the latter we found spin temperature $T_{\mathrm{s}} \sim 75$K, cold gas mass fraction $f\sim0.5$, and turbulent sonic Mach number $M_t\sim3.4$. Similar to the absorption line modeling against 4C~+66.09, our best emission line decomposition model has no unstable gas across the whole field of view, suggesting that the thermal condensation and phase transition are not on-going but rather have reached an equilibrium state. The cold phase of LLIV1 appears as a collection of elongated filaments that forms a closed structure within the field decomposed. These substructures follow the orientation of the overall large scale cloud, along the diagonal of the GHIGLS field from north-west to south-east (in Galactic coordinates). The angular power spectrum of the cold phase is slightly shallower than that of the warm phase, quantifying that the cold phases have relatively more structure on small scales. Our spatially resolved map of the cold gas mass fraction in LLIV1 from DHIGLS reveals significant variations spanning the possible range of $f$, with mean and standard deviation 0.33 and 0.19, respectively.

Abstract: We present the construction of stationary boson-fermion spherically symmetric configurations governed by Newtonian gravity. Bosons are described in the Gross-Pitaevskii regime and fermions are assumed to obey Euler equations for an inviscid fluid with polytropic equation of state. The two components are coupled through the gravitational potential. The families of solutions are parametrized by the central value of the wave function describing the bosons and the central denisty of the fluid. We explore the stability of the solutions using numerical evolutions that solve the time dependent Schr\"odinger-Euler-Poisson system, using the truncation error of the numerical methods as the perturbation. We find that all configurations are stable as long as the polytropic equation of state (EoS) is enforced during the evolution. When the configurations are evolved using the ideal gas EoS they all are unstable that decay into a sort of twin solutions that approach a nearly stationary configuration. We expect these solutions and their evolution serve to test numerical codes that are currently being used in the study of Fuzzy Dark Matter plus baryons.