Astrophysics of Galaxies (astro-ph.GA)

Abstract: Modified Newtonian Dynamics (MOND) is one of the most popular alternative theories of dark matter to explain the missing mass problem in galaxies. Although it remains controversial regarding MOND as a fundamental theory, MOND phenomenology has been shown to widely apply in different galaxies, which gives challenges to the standard $\Lambda$ cold dark matter model. In this article, we derive analytically the galactic rotation curve gradient in the MOND framework and present a rigorous analysis to examine the MOND phenomenology in our Galaxy. By assuming a benchmark baryonic disk density profile and two popular families of MOND interpolating functions, we show for the first time that the recent discovery of the declining Galactic rotation curve in the outer region ($R \approx 17-23$ kpc) can almost rule out the MOND phenomenology at more than $5\sigma$. This strongly supports some of the previous studies claiming that MOND is neither a fundamental theory nor a universal description of galactic properties.

Abstract: Context. Several bright emission line regions have been observed in the central 100 parsecs of the active galaxy NGC 1068. Aims. We aim to determine the properties and ionization mechanism of three regions of NGC 1068: the nucleus (B) and two clouds located at 0.3" and 0.7" north of it (C and D). Methods. We combined SPHERE (0.95 - 1.65 um) and SINFONI (1.5 - 2.45 um) spectra for the three regions B, C, and D. We compared these spectra to several CLOUDY photoionization models and to the MAPPINGS III Library of Fast Radiative Shock Models. Results. The emission line spectra of the three regions are almost identical to each other and contribute to most of the emission line flux in the nuclear region. The emitting media contain multiple phases, the most luminous of which have temperatures ranging from 104.8 K to 106 K. Central photoionization models can reproduce some features of the spectra, but the fast radiative shock model provides the best fit to the data. Conclusions. The similarity between the three regions indicates that they belong to the same class of objects. Based on our comparisons, we conclude that they are shock regions located where the jet of the active galactic nucleus impacts massive molecular clouds.

Abstract: We study gas and metal outflow from massive galaxies in protocluster regions at $z=3-9$ by using the results of the FOREVER22 simulation project. Our simulations contain massive haloes with $M_{\rm h} \gtrsim 10^{13}~\rm M_{\odot}$, showing high star formation rates of $> 100~\rm M_{\odot}~yr^{-1}$ and hosting supermassive black holes with $M_{\rm BH} \gtrsim 10^{8}~\rm M_{\odot}$. We show that the mass loading factor ($\eta_{\rm M}$) sensitively depends on the halo mass and it is $\eta_{\rm M} = 1.2~(9.2)$ for $M_{\rm h} = 10^{13}~(10^{11})~\rm M_{\odot}$. Once the halo mass exceeds $\sim 10^{12.5}~\rm M_{\odot}$, the outflow velocity of the gas rapidly decreases near a virial radius, and the gas returns to a galactic centre finally as a fountain flow. Also, the metal inflow and outflow rates sensitively depend on the halo mass and redshift. At $z=3$, the inflow rate becomes larger than the outflow one if $M_{\rm h} \gtrsim 10^{13.0}~\rm M_{\odot}$. Thus, we suggest that massive haloes cannot be efficient metal enrichment sources beyond virial radii that will be probed in future observations, e.g., studies of metal absorption lines with the Prime Focus Spectrograph on the Subaru telescope.

Abstract: Growing evidence in support of a connection between Active Galactic Nuclei (AGN) activity and the Ram-Pressure Stripping (RPS) phenomenon has been found both observationally and theoretically in the past decades. In this work, we further explore the impact of RPS on the AGN activity by estimating the gas-phase metallicity of nuclear regions and the mass-metallicity relation of galaxies at $z \leq$ 0.07 and with stellar masses $\log {\rm M}_* / {\rm M}_\odot \geq 9.0 $, either experiencing RPS or not. To measure oxygen abundances, we exploit Integral Field Spectroscopy data from the GASP and MaNGA surveys, photoionization models generated with the code CLOUDY and the code Nebulabayes to compare models and observations. In particular, we build CLOUDY models to reproduce line ratios induced by photoionization from stars, AGN, or a contribution of both. We find that the distributions of metallicity and [O III]$\lambda$5007 luminosity of galaxies undergoing RPS are similar to the ones of undisturbed galaxies. Independently of the RPS, we do not find a correlation between stellar mass and AGN metallicity in the mass range $\log {\rm M}_* / {\rm M}_\odot \geq 10.4$, while for the star-forming galaxies we observe the well-known mass-metallicity relation (MZR) between $ 9.0 \leq \log \ {\rm M}_* /{\rm M}_\odot \leq 10.8$ with a scatter mainly driven by the star-formation rate (SFR) and a plateau around $\log {\rm M}_* / {\rm M}_\odot \sim 10.5$. The gas-phase metallicity in the nuclei of AGN hosts is enhanced with respect to those of SF galaxies by a factor of $\sim$ 0.05 dex regardless of the RPS.

Abstract: The spatial metallicity distribution of star clusters in the Small Magellanic Cloud (SMC) has recently been found to correlate as a V-shaped function with the semi-major axis of an elliptical framework proposed to assume a projected galaxy flattening. We report results on the impact that the use of such a framework can produce on our understanding of the SMC formation and its chemical enrichment. We show that clusters with similar semi-major axes are placed at a very different distances from the SMC centre. The recently claimed bimodal metallicity distribution of clusters projected on the innermost SMC regions and the V-shaped metallicity gradient fade away when actual distances are used. Although a large dispersion prevails, clusters older than $\sim$ 1 Gyr exhibit a shallow metallicity gradient, caused by slightly different spatial distributions of clusters younger and older than $\sim$ 4 Gyr; the former being more centrally concentrated and having a mean metallicity ([Fe/H]) $\sim$ 0.15 dex more metal-rich than that of older clusters. This metallicity gradient does not show any dependence with the position angle, except for clusters placed beyond 11 kpc, which are located in the eastern side of the galaxy.

Abstract: We present the median surface brightness profiles of diffuse Ly$\alpha$ haloes (LAHs) around star-forming galaxies by stacking 155 spectroscopically confirmed Ly$\alpha$ emitters (LAEs) at 3<z<4 in the MUSE Extremely Deep Field (MXDF), with median Ly$\alpha$ luminosity $\mathrm{L_{Ly\alpha} \approx 10^{41.1} erg\,s^{-1}}$. After correcting for a systematic surface brightness offset we identified in the datacube, we detect extended Ly$\alpha$ emission out to a distance of 270 kpc. The median Ly$\alpha$ surface brightness profile shows a power-law decrease in the inner 20 kpc, and a possible flattening trend at larger distance. This shape is similar for LAEs with different Ly$\alpha$ luminosities, but the normalisation of the surface brightness profile increases with luminosity. At distances larger than 50 kpc, we observe strong overlap of adjacent LAHs, and the Ly$\alpha$ surface brightness is dominated by the LAHs of nearby LAEs. We find no clear evidence of redshift evolution of the observed Ly$\alpha$ profiles when comparing with samples at 4<z<5 and 5<z<6. Our results are consistent with a scenario in which the inner 20 kpc of the LAH is powered by star formation in the central galaxy, while the LAH beyond a radius of 50 kpc is dominated by photons from surrounding galaxies.

Abstract: We report multi-frequency VLBI studies of the sub-parsec scale structure of the two-sided jet in the nearby radio galaxy NGC 4261. Our analyses include new observations using the Source Frequency Phase Referencing technique with the Very Long Baseline Array at 44 and 88 GHz, as well as archival data at 15 and 43 GHz. Our results show an extended double-sided structure at 43/44 GHz and provide a clear image of the nuclear region at 88 GHz, showing a core size of $\sim$0.09 mas and a brightness temperature of $\sim1.3\times10^{9}$ K. Proper motions are measured for the first time in the two-sided jet, with apparent speeds ranging from $0.31\pm0.14\,c$ to $0.59\pm0.40\,c$ in the approaching jet and $0.32\pm0.14\,c$ in the receding jet. The jet-to-counter-jet brightness ratio allows us to constrain the viewing angle to between $\sim54^{\circ}$ and $84^{\circ}$ and the intrinsic speed to between $\sim0.30\,c$ and $0.55\,c$. We confirm the parabolic shape of the upstream jet on both sides of the central engine, with a power-law index of $0.56\pm0.07$. Notably, the jet collimation is found to be already completed at sub-parsec scales, with a transition location of about 0.61 pc, which is significantly smaller than the Bondi radius of 99.2 pc. This behavior can be interpreted as the initial confinement of the jet by external pressure from either the geometrically thick, optically thin advection-dominated accretion flows (ADAF) or the disk wind launched from it. Alternatively, the shape transition may also be explained by the internal flow transition from a magnetically dominated to a particle-dominated regime.

Abstract: The conditions required for massive star formation are debated, particularly whether massive stars must form in conjunction with massive clusters. Some authors have advanced the view that stars of any mass (below the total cluster mass) can form in clusters of any mass with some probability (random sampling). Others pointed out that the scatter in the determinations of the most massive star mass for a given cluster mass was consistent with the measurement error, such that the mass of the most massive star was determined by the total cluster mass (optimal sampling). Here we investigate the relation between cluster mass (M\textsubscript{ecl}) and the maximum stellar mass (M\textsubscript{max}) using a suite of SPH simulations. Varying cloud mass and turbulence random seed results in a range of cluster masses which we compare with their respective maximum star masses. We find that more massive clusters will have, on average, higher mass stars with this trend being steeper at lower cluster masses ($M\textsubscript{max} \propto M\textsubscript{ecl}^{0.31}$ for $M\textsubscript{ecl}<500M\,_{\odot}$) and flattening at higher cluster masses ($M\textsubscript{max} \propto M\textsubscript{ecl}^{0.11}$ for $M\textsubscript{ecl}>500M\,_{\odot}$). This rules out purely stochastic star formation in our simulations. Significant scatter in the maximum masses with identical initial conditions also rules out the possibility that the relation is purely deterministic (that is that a given cluster mass will result in a specific maximum stellar mass). In conclusion our simulations disagree with both random and optimal sampling of the initial mass function.