Search

The fraction of low-mass galaxies hosting an intermediate-mass black hole (IMBH, with masses $M_{\rm BH} \approx 10^2-10^5$ M$_\odot$), is sensitive to how black hole seeds formed in the early Universe but is observationally still unconstrained. In this paper, we assemble a sample of dwarf galaxies within 10 Mpc hosting bright nuclear star clusters (NSCs) that could host IMBHs. For a subset of them, we use their observed surface brightness from {\it Hubble Space Telescope} (\hst) images, an assumed synthetic spectrum of their stellar population, Jeans Anisotropic Model (JAM) of the stellar dynamics, and the {\tt HSIM} simulator software to create mock observations with the High Angular Resolution Monolithic Optical and Near-infrared Integral (HARMONI) field spectrograph for the Extremely Large Telescope (ELT). We analyze the simulated data cube like real data, using JAM to infer the IMBH mass and its error in a Bayesian framework. Our simulations show that the ELT/HARMONI instrument can clearly detect the existence of IMBH demographics in NSCs down to a mass of about 0.5\% of the NSC., Comment: 33 pages, 19 figures, 9 tables, submitted to MNRAS

NGC~613 is a nearby barred spiral galaxy with a nuclear ring. Exploiting high spatial resolution ($\approx20$ pc) Atacama Large Millimeter/sub-millimeter Array $^{12}$CO(1-0) observations, we study the giant molecular clouds (GMCs) in the nuclear ring and its vicinity, identifying $158$ spatially- and spectrally-resolved GMCs. The GMC sizes ($R_{\mathrm{c}}$) are comparable to those of the clouds in the Milky Way (MW) disc, but their gas masses, observed linewidths ($\sigma_{\mathrm{obs,los}}$) and gas mass surface densities are larger. The GMC size -- linewidth relation ($\sigma_{\mathrm{obs,los}}\propto R_{\mathrm{c}}^{0.77}$) is steeper than that of the clouds of the MW disc and centre, and the GMCs are on average only marginally gravitationally bound (with a mean virial parameter $\langle\alpha_{\mathrm{obs,vir}}\rangle\approx1.7$). We discuss the possible origins of the steep size -- linewidth relation and enhanced observed linewidths of the clouds and suggest that a combination of mechanisms such as stellar feedback, gas accretion and cloud-cloud collisions, as well as the gas inflows driven by the large-scale bar, may play a role., Comment: 15 pages, 8 figures, accepted for publication in MNRAS. arXiv admin note: text overlap with arXiv:2304.10471

We measure the age-velocity relationship from the lag between ionized gas and stellar tangential speeds in ~500 nearby disk galaxies from MaNGA in SDSS-IV. Selected galaxies are kinematically axisymmetric. Velocity lags are asymmetric drift, seen in the Milky Way's (MW) solar neighborhood and other Local Group galaxies; their amplitude correlates with stellar population age. The trend is qualitatively consistent in rate (d(sigma)/dt) with a simple power-law model where sigma is proportional to t^b that explains the dynamical phase-space stratification in the solar neighborhood. The model is generalized based on disk dynamical times to other radii and other galaxies. We find in-plane radial stratification parameters sigma_(0,r} (dispersion of the youngest populations) in the range of 10-40 km/s and 0.2

The mass ($M_\mathrm{BH}$) of a supermassive black hole (SMBH) can be measured using spatially-resolved kinematics of the region where the SMBH dominates gravitationally. The most reliable measurements are those that resolve the smallest physical scales around the SMBHs. We consider here three metrics to compare the physical scales probed by kinematic tracers dominated by rotation: the radius of the innermost detected kinematic tracer $R_\mathrm{min}$ normalised by respectively the SMBH's Schwarzschild radius ($R_\mathrm{Schw}\equiv 2GM_\mathrm{BH}/c^2$, where $G$ is the gravitational constant and $c$ the speed of light), sphere-of-influence (SOI) radius ($R_\mathrm{SOI}\equiv GM_\mathrm{BH}/\sigma_\mathrm{e}^2$, where $\sigma_\mathrm{e}$ is the stellar velocity dispersion within the galaxy's effective radius) and equality radius [the radius $R_\mathrm{eq}$ at which the SMBH mass equals the enclosed stellar mass, $M_\mathrm{BH}=M_*(R_\mathrm{eq})$, where $M_*(R)$ is the stellar mass enclosed within the radius $R$]. All metrics lead to analogous simple relations between $R_\mathrm{min}$ and the highest circular velocity probed $V_\mathrm{c}$. Adopting these metrics to compare the SMBH mass measurements using molecular gas kinematics to those using megamaser kinematics, we demonstrate that the best molecular gas measurements resolve material that is physically closer to the SMBHs in terms of $R_\mathrm{Schw}$ but is slightly farther in terms of $R_\mathrm{SOI}$ and $R_\mathrm{eq}$. However, molecular gas observations of nearby galaxies using the most extended configurations of the Atacama Large Millimeter/sub-millimeter Array can resolve the SOI comparably well and thus enable SMBH mass measurements as precise as the best megamaser measurements., Comment: 13 pages, 9 figures. Accepted by MNRAS

Supermassive black hole (SMBH) masses can be measured by observing the impacts of the SMBHs on dynamical tracers around them. We present high angular resolution ($0.19$ arcsec or $\approx24$ pc) Atacama Large Millimeter/submillimeter Array observations of the $^{12}$CO(3-2) line emission of the early-type galaxy NGC 4751, which reveal a highly-inclined regularly-rotating molecular gas disc with clear central Keplerian motions. Using a Hubble Space Telescope image to constrain the stellar mass distribution, we forward model the molecular gas kinematics and data cube in a Bayesian framework using the Kinematic Molecular Simulation code. Assuming a constant mass-to-light ratio ($M/L$), we infer a SMBH mass $M_\text{BH}=3.43^{+0.45}_{-0.44}\times10^9$ $\text{M}_\odot$ and a F160W filter stellar $M/L$ $M/L_\text{F160W}=(2.68\pm0.11)$ $\text{M}_\odot/\text{L}_{\odot,\text{F160W}}$ (all quoted uncertainties are at $3\sigma$ confidence). Assuming a linearly spatially-varying $M/L$, we infer $M_\text{BH}=2.79_{-0.57}^{+0.75}\times10^9$ $\text{M}_\odot$ and $\left(M/L_\text{F160W}\right)/\left(\text{M}_\odot/\text{L}_{\odot,\text{F160W}}\right)=3.07^{+0.27}_{-0.35}-0.09^{+0.08}_{-0.06}\,\left(R/\text{arcsec}\right)$, where $R$ is the galactocentric radius. We also present alternative SMBH mass estimates using the Jeans Anisotropic Modelling (JAM) method and SINFONI stellar kinematics. Assuming a cylindrically-aligned velocity ellipsoid (JAM$_\text{cyl}$) we infer $M_\text{BH}=(2.52\pm 0.36)\times10^9$ $\text{M}_\odot$, while assuming a spherically-aligned velocity ellipsoid (JAM$_\text{sph}$) we infer $M_\text{BH}=(3.24\pm0.87)\times10^9$ $\text{M}_\odot$. Our derived masses are all consistent with one another, but they are larger than (and inconsistent with) one previous stellar dynamical measurement using Schwarzschil's method and the same SINFONI kinematics., Comment: 15 pages, 13 figures, submitted to MNRAS. arXiv admin note: text overlap with arXiv:2401.16376

We have re-observed $\rm\sim$40 low-inclination, star-forming galaxies from the MaNGA survey ($\upsigma\sim65$~\kms) at $\sim$6.5 times higher spectral resolution ($\upsigma\sim10$~\kms) using the HexPak integral field unit on the WIYN 3.5m telescope. The aim of these observations is to calibrate MaNGA's instrumental resolution and to characterize turbulence in the warm interstellar medium and ionized galactic outflows. Here we report the results for the H$\rm\upalpha$ region observations as they pertain to the calibration of MaNGA's spectral resolution. Remarkably, we find that the previously-reported MaNGA line-spread-function (LSF) Gaussian width is systematically underestimated by only 1\%. The LSF increase modestly reduces the characteristic dispersion of HII regions-dominated spectra sampled at 1-2 kpc spatial scales from 23 to 20 km s$^{-1}$ in our sample, or a 25\% decrease in the random-motion kinetic energy. This commensurately lowers the dispersion zeropoint in the relation between line-width and star-formation rate surface-density in galaxies sampled on the same spatial scale. This modest zero-point shift does not appear to alter the power-law slope in the relation between line-width and star-formation rate surface-density. We also show that adopting a scheme whereby corrected line-widths are computed as the square root of the median of the difference in the squared measured line width and the squared LSF Gaussian avoids biases and allows for lower SNR data to be used reliably.

The INSPIRE project has built the largest sample of ultra-compact massive galaxies (UCMGs) at 0.12, and the other stellar population parameters.We use the Full-Index-Fitting method to fit the INSPIRE spectra to single stellar population (SSP) models. This allows us to measure, for the first time, the low-mass end slope of the IMF, as well as stellar metallicity [M/H], [Mg/Fe], [Ti/Fe] and [Na/Fe] ratios, and study correlations between them and the DoR. Similarly to normal-sized galaxies, UCMGs with larger stellar masses have overall higher metallicities. We found a correlation between the low-mass end of the IMF slope and the DoR, that, however, breaks down for systems with a more extended SFH. An even stronger dependency is found between the IMF and the fraction of mass formed at high-z. At equal velocity dispersion and metallicity, galaxies with a higher DoR have a dwarf-richer IMF than that of low-DoR counterparts. This might indicate that the cosmic epoch and formation mechanisms influence the fragmentation of the star formation cloud and hence might be the explanation for IMF variations detected in massive ETGs., Comment: 16 pages, 11 figures, submitted to MNRAS

We use high-resolution data from the millimetre-Wave Interferometric Survey of Dark Object Masses (WISDOM) project to investigate the connection between circumnuclear gas reservoirs and nuclear activity in a sample of nearby galaxies. Our sample spans a wide range of nuclear activity types including radio galaxies, Seyfert galaxies, low-luminosity active galactic nuclei (AGN) and inactive galaxies. We use measurements of nuclear millimetre continuum emission along with other archival tracers of AGN accretion/activity to investigate previous claims that at, circumnuclear scales (<100 pc), these should correlate with the mass of the cold molecular gas. We find that the molecular gas mass does not correlate with any tracer of nuclear activity. This suggests the level of nuclear activity cannot solely be regulated by the amount of cold gas around the supermassive black hole (SMBH). This indicates that AGN fuelling, that drives gas from the large scale galaxy to the nuclear regions, is not a ubiquitous process and may vary between AGN type, with timescale variations likely to be very important. By studying the structure of the central molecular gas reservoirs, we find our galaxies have a range of nuclear molecular gas concentrations. This could indicate that some of our galaxies may have had their circumnuclear regions impacted by AGN feedback, even though they currently have low nuclear activity. On the other hand, the nuclear molecular gas concentrations in our galaxies could instead be set by secular processes., Comment: 15 pages plus 3 in the appendix, 8 figures plus 1 in the appendix, 3 tables plus 4 in the appendix

The co-evolution of galaxies and supermassive black holes (SMBHs) underpins our understanding of galaxy evolution, but different methods to measure SMBH masses have only infrequently been cross-checked. We attempt to identify targets to cross-check two of the most accurate methods, megamaser and cold molecular gas dynamics. Three promising galaxies are selected from all those with existing megamaser SMBH mass measurements. We present Atacama Large Millimeter/sub-millimeter Array (ALMA) 12CO(2-1) and 230-GHz continuum observations with angular resolutions of about 0.5". Every galaxy has an extended rotating molecular gas disc and 230-GHz continuum source(s), but all also have irregularities and/or non-axisymmetric features: NGC1194 is highly inclined and has disturbed and lopsided central 12CO(2-1) emission; NGC3393 has a nuclear disc with fairly regular but patchy 12CO(2-1) emission with little gas near the kinematic major axis, faint emission in the very centre and two brighter structures reminiscent of a nuclear ring and/or spiral; NGC5765B has a strong bar and very bright 12CO(2-1) emission concentrated along two bisymmetric offset dust lanes and two bisymmetric nuclear spiral arms. 12CO(2-1) and 12CO(3-2) observations with the James Clerk Maxwell Telescope are compared with the ALMA observations. Because of the disturbed gas kinematics and the impractically long integration times required for higher angular resolution observations, none of the three galaxies is suitable for a future SMBH mass measurement. Nonetheless, increasing the number of molecular gas observations of megamaser galaxies is valuable, and the ubiquitous disturbances suggest a link between large-scale gas properties and the existence of megamasers., Comment: 17 pages, 5 figures, accepted by MNRAS

We try to understand the trends in the mass density slopes as a function of galaxy properties. We use the results from the best Jeans Anisotropic Modelling (JAM) of the integral-field stellar kinematics for near 6000 galaxies from the MaNGA DynPop project, with stellar masses of $10^{9-12}\ {\rm M_{\odot}}$, including both early-type and late-type galaxies. We use the mass-weighted density slopes for the stellar $\overline{\gamma}_*$, dark $\overline{\gamma}_{\rm DM}$, and total $\overline{\gamma}_{\rm T}$ mass from the MaNGA DynPop project. The $\overline{\gamma}_{\rm T}$ approaches a constant value of 2.2 for high $\sigma_{\rm e}$ galaxies, and flattens for lg$(\sigma_{\rm e}/{\rm km\ s^{-1}})\lesssim2.3$, reaching 1.5 for lg$(\sigma_{\rm e}/{\rm km\ s^{-1}})\approx1.8$. The total and stellar slopes track each other tightly, with $\overline{\gamma}_{\rm T}\approx\overline{\gamma}_*-0.174$ over the full $\sigma_{\rm e}$ range. This confirms the dominance of stellar matter within $R_{\rm e}$. We also show that there is no perfect conspiracy between baryonic and dark matter, as $\overline{\gamma}_*$ and $\overline{\gamma}_{\rm DM}$ do not vary inversely within the $\sigma_{\rm e}$ range. We find that the central galaxies from TNG50 and TNG100 simulations do not reproduce the observed galaxy mass distribution, which we attribute to the overestimated dark matter fraction, possibly due to a constant IMF and excessive adiabatic contraction effects in the simulations. Finally, we present the stacked dark matter density profiles and show that they are slightly steeper than the pure dark matter simulation prediction of $\overline{\gamma}_{\rm DM}\approx1$, suggesting moderate adiabatic contraction in the central region of galaxies. Our work demonstrate the power of stellar dynamics modelling for probing the interaction between stellar and dark matter and testing galaxy formation theories., Comment: 17 pages, 14 figures; accepted by MNRAS

Observations indicate that the central gas discs are smoother in early-type galaxies than their late-type counterparts, while recent simulations predict that the dynamical suppression of star formation in spheroid-dominated galaxies is preceded by the suppression of fragmentation of their interstellar media. The mass surface density power spectrum is a powerful tool to constrain the degree of structure within a gas reservoir. Specifically here, we focus on the power spectrum slope and aim to constrain whether the shear induced by a dominant spheroidal potential can induce sufficient turbulence to suppress fragmentation, resulting in the smooth central gas discs observed. We compute surface density power spectra for the nuclear gas reservoirs of fourteen simulated isolated galaxies and twelve galaxies observed as part of the mm-Wave Interferometric Survey of Dark Object Masses (WISDOM) project. Both simulated and observed galaxies range from disc-dominated galaxies to spheroids, with central stellar mass surface densities, a measure of bulge dominance, varying by more than an order of magnitude. For the simulations, the power spectra steepen with increasing central stellar mass surface density, thereby clearly linking the suppression of fragmentation to the shear-driven turbulence induced by the spheroid. The WISDOM observations show a different (but potentially consistent) picture: while there is no correlation between the power spectrum slopes and the central stellar mass surface densities, the slopes scatter around a value of 2.6. This is similar to the behaviour of the slopes of the simulated galaxies with high central stellar mass surface densities, and could indicate that high shear eventually drives incompressible turbulence., Comment: 14 pages, 7 figures + appendices (7 pages, 9 figures); accepted by MNRAS (October 09th, 2023)

Using the final MaNGA sample of 10K galaxies, we investigate the dark matter fraction $f_{\rm DM}$ within one half-light radius $R_{\rm e}$ for about 6K galaxies with good kinematics spanning a wide range of morphologies and stellar velocity dispersion. We employ two techniques to estimate $f_{\rm DM}$: (i) Jeans Anisotropic Modelling (JAM), which performs dark matter decomposition based on stellar kinematics and (ii) comparing the total dynamical mass-to-light ratios $(M/L)_{\rm JAM}$ and $(M_{\ast}/L)_{\rm SPS}$ from Stellar Population Synthesis (SPS). We find that both methods consistently show a significant trend of increasing $f_{\rm DM}$ with decreasing $\sigma_{\rm e}$ and low $f_{\rm DM}$ at larger $\sigma_{\rm e}$. For 235 early-type galaxies with the best models, we explore the variation of stellar initial mass function (IMF) by comparing the stellar mass-to-light ratios from JAM and SPS. We confirm that the stellar mass excess factor $\alpha_{\rm IMF}$ increases with $\sigma_{\rm e}$, consistent with previous studies that reported a transition from Chabrier-like to Salpeter IMF among galaxies. We show that the $\alpha_{\rm IMF}$ trend cannot be driven by $M_{\ast}/L$ or IMF gradients as it persists when allowing for radial gradients in our model. We find no evidence for the total $M/L$ increasing toward the centre. We detect weak positive correlations between $\alpha_{\rm IMF}$ and age, but no correlations with metallicity. We stack galaxy spectra according to their $\alpha_{\rm IMF}$ to search for differences in IMF-sensitive spectral features (e.g. the $\rm Na_{\rm I}$ doublet). We only find marginal evidence for such differences, which casts doubt on the validity of one or both methods to measure the IMF., Comment: 20 pages, 12 figures, 2 tables; Accepted for publication in MNRAS on 22 April 2024

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., Comment: 32 pages (16 of Appendices), 39 Figures (27 in Appendices). Accepted for publication in MNRAS

We report the discovery of the ``mm fundamental plane of black-hole accretion'', which is a tight correlation between the nuclear 1 mm luminosity ($L_{\rm \nu, mm}$), the intrinsic $2$ -- $10$~keV X-ray luminosity ($L_{\rm X,2-10}$) and the supermassive black hole (SMBH) mass ($M_{\rm BH}$) with an intrinsic scatter ($\sigma_{\rm int}$) of $0.40$ dex. The plane is found for a sample of 48 nearby galaxies, most of which are low-luminosity active galactic nuclei (LLAGN). Combining these sources with a sample of high-luminosity (quasar-like) nearby AGN, we find that the plane still holds. We also find that $M_{\rm BH}$ correlates with $L_{\rm \nu, mm}$ at a highly significant level, although such correlation is less tight than the mm fundamental plane ($\sigma_{\rm int}=0.51$ dex). Crucially, we show that spectral energy distribution (SED) models for both advection-dominated accretion flows (ADAFs) and compact jets can explain the existence of these relations, which are not reproduced by the standard torus-thin accretion disc models usually associated to quasar-like AGN. The ADAF models reproduces the observed relations somewhat better than those for compact jets, although neither provides a perfect prediction. Our findings thus suggest that radiatively-inefficient accretion processes such as those in ADAFs or compact (and thus possibly young) jets may play a key role in both low- and high-luminosity AGN. This mm fundamental plane also offers a new, rapid method to (indirectly) estimate SMBH masses., Comment: 9 pages, 5 figures

The second moment of the stellar velocity within the effective radius, denoted by $\sigma_{\rm e}^2$, is a crucial quantity in galaxy studies as it provides insight into galaxy properties and their mass distributions. However, large spectroscopic surveys typically do not measure $\sigma_{\rm e}$ directly, instead providing $\sigma_{\rm aper}$, the second moment of the stellar velocity within a fixed fiber aperture. In this paper, we derive an empirical aperture correction formula, given by $\sigma_{\rm aper}/\sigma_{\rm e}=(R_{\rm aper}/R_{\rm e})^{\alpha}$, using spatially resolved stellar kinematics extracted from approximately 10,000 Sloan Digital Sky Survey-Mapping Nearby Galaxies at Apache Point Observatory (SDSS-MaNGA) integral field unit observations. Our analysis reveals a strong dependence of $\alpha$ on the $r$-band absolute magnitude $M_{\rm r}$, $g-i$ color, and Sersic index $n_{\rm Ser}$, where $\alpha$ values are lower for brighter, redder galaxies with higher Sersic indices. Our results demonstrate that the aperture correction derived from previous literature on early-type galaxies cannot be applied to predict the aperture corrections for galaxies with intermediate Sersic indices. We provide a lookup table of $\alpha$ values for different galaxy types, with parameters in the ranges of $-18>M_{\rm r}>-24$, $0.4

We analyse the integral-field spectroscopy data for the $\approx10,000$ galaxies in final data release of the MaNGA survey. We identify 188 galaxies for which the emission lines cannot be described by single Gaussian components. These galaxies can be classified into (1) 38 galaxies with broad $H\alpha$ and [OIII] $\lambda$5007 lines, (2) 101 galaxies with broad $H\alpha$ lines but no broad [OIII] $\lambda$5007 lines, and (3) 49 galaxies with double-peaked narrow emission lines. Most of the broad line galaxies are classified as Active Galactic Nuclei (AGN) from their line ratios. The catalogue helps us further understand the AGN-galaxy coevolution through the stellar population of broad-line region host galaxies and the relation between broad lines' properties and the host galaxies' dynamical properties. The stellar population properties (including mass, age and metallicity) of broad-line host galaxies suggest there is no significant difference between narrow-line Seyfert-2 galaxies and Type-1 AGN with broad $H\alpha$ lines. We use the broad-$H\alpha$ line width and luminosity to estimate masses of black hole in these galaxies, and test the $M_{BH}-\sigma_{e}$ relation in Type-1 AGN host galaxies. Furthermore we find three dual AGN candidates supported by radio images from the VLA FIRST survey. This sample may be useful for further studies on AGN activities and feedback processes., Comment: 21 pages, 17 figures, LaTeX. Accepted by MNRAS

We present the measurement of total and stellar/dark matter decomposed mass density profile around a sample of galaxy groups and clusters with dynamical masses derived from integral-field stellar kinematics from the MaNGA survey in Paper I and weak lensing derived from the DECaLS imaging survey. Combining the two data sets enables accurate measurement of the radial density distribution from several kpc to Mpc scales. Intriguingly, we find that the excess surface density derived from stellar kinematics in the inner region cannot be explained by simply adding an NFW dark matter halo extrapolated from lensing measurement at a larger scale to a stellar mass component derived from the NASA-Sloan Atlas (NSA) catalogue. We find that a good fit to both data sets requires a stellar mass normalization about 3 times higher than that derived from the NSA catalogue, which would require an unrealistically too-heavy initial mass function for stellar mass estimation. If we keep the stellar mass normalization to that of the NSA catalogue but allow a varying inner dark matter density profile, we obtain an asymptotic slope of $\gamma_{\rm gnfw}$= $1.82_{-0.25}^{+0.15}$ and $\gamma_{\rm gnfw}$= $1.48_{-0.41}^{+0.20}$ for the group bin and the cluster bin respectively, significantly steeper than the NFW case. We also compare the total mass inner density slopes with those from TNG300 and find that the values from the simulation are lower than the observation by about $2\sigma$ level., Comment: 12 pages, 9 figures, accepted by MNRAS

We present dynamical scaling relations, combined with the stellar population properties, for a subsample of about 6000 nearby galaxies with the most reliable dynamical models extracted from the full MaNGA sample of 10K galaxies. We show that the inclination-corrected mass plane (MP) for both early-type galaxies (ETGs) and late-type galaxies (LTGs), which links dynamical mass, projected half-light radius $R_{\rm e}$, and the second stellar velocity moment $\sigma_{\rm e}$ within $R_{\rm e}$, satisfies the virial theorem and is even tighter than the uncorrected one. We find a clear parabolic relation between $\lg(M/L)_{\rm e}$, the total mass-to-light ratio within a sphere of radius $R_{\rm e}$, and $\lg\sigma_{\rm e}$, with the $M/L$ increasing with $\sigma_{\rm e}$ and for older stellar populations. However, the relation for ETGs is linear and the one for the youngest galaxies is constant. We confirm and improve the relation between average logarithmic total density slopes $\overline{\gamma_{\rm T}}$ and $\sigma_{\rm e}$: $\overline{\gamma_{\rm T}}$ become steeper with increasing $\sigma_{\rm e}$ until $\lg(\sigma_{\rm e}/{\rm km\,s^{-1}})\approx 2.2$ and then remain constant around $\overline{\gamma_{\rm T}}\approx -2.2$. The $\overline{\gamma_{\rm T}}-\sigma_{\rm e}$ variation is larger for LTGs than ETGs. At fixed $\sigma_{\rm e}$ the total density profiles steepen with galaxy age and for ETGs. We find generally low dark matter fractions, median $f_{\rm DM}(

This is the second paper of the MaNGA DynPop series, which analyzes the global stellar population, radial gradients, and non-parametric star-formation history of $\sim 10$K galaxies from the MaNGA Survey final data release 17 (DR17) and relates them with dynamical properties of galaxies. We confirm the correlation between the stellar population properties and the stellar velocity dispersion $\sigma_{\rm e}$, but also find that younger galaxies are more metal-poor at fixed $\sigma_{\rm e}$. Stellar age, metallicity, and mass-to-light ratio $M_{\ast}/L$ all decrease with increasing galaxy rotation, while their radial gradients become more negative (i.e. lower value at the outskirts). The exception is the slow rotators, which also appear to have significantly negative metallicity gradients, confirming the mass-metallicity gradient correlation. Massive disk galaxies in the green valley, on the $(\sigma_{\rm e},\rm Age)$ plane, show the most negative age and metallicity gradients, consistent with their old central bulges surrounded by young star-forming disks and metal-poor gas accretion. Galaxies with high $\sigma_{\rm e}$, steep total mass-density slope, low dark matter fraction, high $M_{\ast}/L$, and high metallicity have the highest star-formation rate at earlier times, and are currently quenched. We also discover a population of low-mass star-forming galaxies with low rotation but physically distinct from the massive slow rotators. A catalogue of these stellar population properties is provided publicly., Comment: 26 pages, 20 figures, 1 table. Accepted for publication in MNRAS. The catalogue is available at https://manga-dynpop.github.io/

This is the first paper in our series on the combined analysis of the Dynamics and stellar Population (DynPop) for the MaNGA survey in the final SDSS Data Release 17 (DR17). Here we present a catalogue of dynamically-determined quantities for over 10000 nearby galaxies based on integral-field stellar kinematics from the MaNGA survey. The dynamical properties are extracted using the axisymmetric Jeans Anisotropic Modelling (JAM) method, which was previously shown to be the most accurate for this kind of study. We assess systematic uncertainties using eight dynamical models with different assumptions. We use two orientations of the velocity ellipsoid: either cylindrically-aligned JAM$_{\rm cyl}$ or spherically-aligned JAM$_{\rm sph}$. We also make four assumptions for the models' dark vs. luminous matter distributions: (1) mass-follows-light, (2) free NFW dark halo, (3) cosmologically-constrained NFW halo, (4) generalized NFW dark halo, i.e. with free inner slope. In this catalogue, we provide the quantities related to the mass distributions (e.g. the density slopes and enclosed mass within a sphere of a given radius for total mass, stellar mass, and dark matter mass components). We also provide the complete models which can be used to compute the full luminous and mass distribution of each galaxy. Additionally, we visually assess the qualities of the models to help with model selections. We estimate the observed scatter in the measured quantities which decreases as expected with improvements in quality. For the best data quality, we find a remarkable consistency of measured quantities between different models, highlighting the robustness of the results., Comment: 30 pages, 15 figures, 4 tables, accepted by MNRAS, the catalogue is available here https://manga-dynpop.github.io/

We present high spatial resolution ($\approx24$ pc) Atacama Large Millimeter/sub-millimeter Array $^{12}$CO(2-1) observations of the central region of the nearby barred spiral galaxy NGC 5806. NGC 5806 has a highly structured molecular gas distribution with a clear nucleus, a nuclear ring and offset dust lanes. We identify $170$ spatially- and spectrally-resolved giant molecular clouds (GMCs). These clouds have comparable sizes ($R_{\mathrm{c}}$) and larger gas masses, observed linewidths ($\sigma_{\mathrm{obs,los}}$) and gas mass surface densities than those of clouds in the Milky Way disc. The size -- linewidth relation of the clouds is one of the steepest reported so far ($\sigma_{\mathrm{obs,los}}\propto R_{\mathrm{c}}^{1.20}$), the clouds are on average only marginally bound (with a mean virial parameter $\langle\alpha_{\mathrm{vir}}\rangle\approx2$), and high velocity dispersions are observed in the nuclear ring. These behaviours are likely due to bar-driven gas shocks and inflows along the offset dust lanes, and we infer an inflow velocity of $\approx120$ kms$^{-1}$ and a total molecular gas mass inflow rate of $\approx5$ M$_\odot$ yr$^{-1}$ into the nuclear ring. The observed internal velocity gradients of the clouds are consistent with internal turbulence. The number of clouds in the nuclear ring decreases with azimuthal angle downstream from the dust lanes without clear variation of cloud properties. This is likely due to the estimated short lifetime of the clouds ($\approx6$ Myr), which appears to be mainly regulated by cloud-cloud collision and/or shear processes. Overall, it thus seems that the presence of the large-scale bar and gas inflows to the centre of NGC 5806 affect cloud properties., Comment: Accepted for publication in MNRAS, 20 pages, 16 figures

We present a CO dynamical estimate of the mass of the super-massive black hole (SMBH) in three nearby early-type galaxies: NGC0612, NGC1574 and NGC4261. Our analysis is based on Atacama Large Millimeter/submillimeter Array (ALMA) Cycle 3-6 observations of the $^{12}$CO(2-1) emission line with spatial resolutions of $14-58$ pc ($0.01"-0.26"$). We detect disc-like CO distributions on scales from $\lesssim200$ pc (NGC1574 and NGC4261) to $\approx10$ kpc (NGC0612). In NGC0612 and NGC1574 the bulk of the gas is regularly rotating. The data also provide evidence for the presence of a massive dark object at the centre of NGC1574, allowing us to obtain the first measure of its mass, $M_{\rm BH}=(1.0\pm0.2)\times10^{8}$ M$_{\odot}$ (1$\sigma$ uncertainty). In NGC4261, the CO kinematics is clearly dominated by the SMBH gravitational influence, allowing us to determine an accurate black hole mass of $(1.62{\pm 0.04})\times10^{9}$ M$_{\odot}$ ($1\sigma$ uncertainty). This is fully consistent with a previous CO dynamical estimate obtained using a different modelling technique. Signs of non-circular gas motions (likely outflow) are also identified in the inner regions of NGC4261. In NGC0612, we are only able to obtain a (conservative) upper limit of $M_{\rm BH}\lesssim3.2\times10^{9}$ M$_{\odot}$. This has likely to be ascribed to the presence of a central CO hole (with a radius much larger than that of the SMBH sphere of influence), combined with the inability of obtaining a robust prediction for the CO velocity curve. The three SMBH mass estimates are overall in agreement with predictions from the $M_{\rm BH}-\sigma_{\star}$ relation., Comment: Main text: 20 pages, 14 Figures; Appendix: 7 pages, 6 Figures. Accepted for publication in MNRAS on 2023 March 28

The massive elliptical galaxy M87 has been the subject of several supermassive black hole mass measurements from stellar dynamics, gas dynamics, and recently the black hole shadow by the Event Horizon Telescope (EHT). This uniquely positions M87 as a benchmark for alternative black hole mass determination methods. Here we use stellar kinematics extracted from integral-field spectroscopy observations with Adaptive Optics (AO) using MUSE and OASIS. We exploit our high-resolution integral field spectroscopy to spectrally decompose the central AGN from the stars. We derive an accurate inner stellar-density profile and find it is flatter than previously assumed. We also use the spectrally-extracted AGN as a reference to accurately determine the observed MUSE and OASIS AO PSF. We then perform Jeans Anisotropic Modelling (JAM), with a new flexible spatially-variable anisotropy, and measure the anisotropy profile, stellar mass-to-light variations, inner dark matter fraction, and black hole mass. Our preferred black hole mass is $M_{\rm BH}=(8.7\pm1.2) \times 10^9 \ M_\odot $. However, using the inner stellar density from previous studies, we find a preferred black hole mass of $M_{\rm BH} = (5.5^{+0.5}_{-0.3}) \times 10^9 \ M_\odot $, consistent with previous work. We conduct numerous systematic tests of the kinematics and model assumptions and conclude that uncertainties in the black hole mass of M87 from previous determinations may have been underestimated and further analyses are needed., Comment: 19 pages, 16 figures, 4 tables, Submitted to MNRAS

The radio jets of an active galactic nucleus (AGN) can heat up the gas around a host galaxy and quench star formation activity. The presence of a radio jet could be related to the evolutionary path of the host galaxy and may be imprinted in the morphology and kinematics of the galaxy. In this work, we use data from the Sloan Digital Sky Survey's Mapping Nearby Galaxies at Apache Point Observatory survey and the Low Frequency Array (LOFAR) Two-Metre Sky Survey as well as the National Radio Astronomy Observatory (NRAO) the Karl G. Jansky Very Large Array (VLA) Sky Survey and the Faint Images of the Radio Sky at Twenty Centimeter survey. We combine these integral field spectroscopic data and radio data to study the link between stellar kinematics and radio AGNs. We find that the luminosity-weighted stellar angular momentum $\lambda_{Re}$ is tightly related to the range of radio luminosity and the fraction of radio AGNs F radio present in galaxies, as high-luminosity radio AGNs are only in galaxies with a small $\lambda_{Re}$, and the $F_{radio}$ at a fixed stellar mass decreases with $\lambda_{Re}$. These results indicate that galaxies with stronger random stellar motions with respect to the ordered motions might be better breeding grounds for powerful radio AGNs. This would also imply that the merger events of galaxies are important in the triggering of powerful radio jets in our sample., Comment: 10 pages,7 figures. Accepted in A&A

According to the $\Lambda$CDM cosmology, present-day galaxies with stellar masses $M_\star>10^{11} {\rm M}_\odot$ should contain a sizable fraction of dark matter within their stellar body. Models indicate that in massive early-type galaxies (ETGs) dark matter should account for $\sim60\%$ of the dynamical mass within five effective radii ($5 R_{\rm e}$). Most massive ETGs have been shaped through a two-phase process: the rapid growth of a compact core was followed by the accretion of an extended envelope through mergers. The exceedingly rare galaxies that have avoided the second phase, the so-called relic galaxies, are thought to be the frozen remains of the massive ETG population at $z\gtrsim2$. The best relic galaxy candidate discovered to date is NGC 1277, in the Perseus cluster. We used deep integral field GCMS data to revisit NGC 1277 out to an unprecedented radius of 6 kpc (corresponding to $5 R_{\rm e}$). By using Jeans anisotropic modelling we find a negligible dark matter fraction within $5 R_{\rm e}$ ($f_{\rm DM}(5 R_{\rm e})<0.05$; two-sigma confidence level), which is in tension with the expectation. Since the lack of an extended envelope would reduce dynamical friction and prevent the accretion of an envelope, we propose that NGC 1277 lost its dark matter very early or that it was dark matter deficient ab initio. We discuss our discovery in the framework of recent proposals suggesting that some relic galaxies may result from dark matter stripping as they fell in and interacted within galaxy clusters. Alternatively, NGC 1277 might have been born in a high-velocity collision of gas-rich proto-galactic fragments, where dark matter left behind a disc of dissipative baryons. We speculate that the relative velocities of $\approx2000 {\rm km/s}$ required for the latter process to happen were possible in the progenitors of the present-day rich galaxy clusters., Comment: Accepted for publication in A&A

As the earliest relics of star formation episodes of the Universe, the most massive galaxies are the key to our understanding of the stellar population, cosmic structure, and SMBH evolution. However, the details of their formation histories remain uncertain. We address these problems by planning a large survey sample of 101 ultramassive galaxies ($z\le0.3$, $|\delta+24^{\circ}|<45$\deg, $|b|>8$\deg), including 76\% ellipticals, 17\% lenticulars, and 7\% spirals brighter than $M_K\le-27$~mag (stellar mass $2\times10^{12}\lesssim M_\star\lesssim5\times10^{12}$~\Msun) with ELT/HARMONI. Our sample comprises diverse galaxy environments ranging from isolated to dense-cluster galaxies. The primary goals of the project are to (1) explore the stellar dynamics inside galaxy nuclei and weigh SMBHs, (2) constrain the black hole scaling relations at the highest mass, and (3) probe the late-time assembly of these most massive galaxies through the stellar population and kinematical gradients. We describe the survey, discuss the distinct demographics and environmental properties of the sample, and simulate their HARMONI $I_z$, $I_z+J$, and $H+K$-band observations via combining the inferred stellar-mass models from Pan-STARRS observations, an assumed synthetic spectrum of stars, and SMBHs with masses estimated based on different black hole scaling relations. Our simulations produce excellent state-of-art IFS and stellar kinematics ($\Delta V_{\rm rms}\lesssim1.5$\%) in a relatively short exposure time. We use these stellar kinematics in combination with the JAM to reconstruct the SMBH mass and its error using the MCMC simulation. Thus, these simulations and modelings can be benchmarks to evaluate the instrument models and pipelines dedicated to HARMONI to exploit the unprecedented capabilities of ELT., Comment: Accepted to MNRAS

Strong-lensing time delays enable measurement of the Hubble constant ($H_{0}$) independently of other traditional methods. The main limitation to the precision of time-delay cosmography is mass-sheet degeneracy (MSD). Some of the previous TDCOSMO analyses broke the MSD by making standard assumptions about the mass density profile of the lens galaxy, reaching 2% precision from seven lenses. However, this approach could potentially bias the $H_0$ measurement or underestimate the errors. For this work, we broke the MSD for the first time using spatially resolved kinematics of the lens galaxy in RXJ1131$-$1231 obtained from the Keck Cosmic Web Imager spectroscopy, in combination with previously published time delay and lens models derived from Hubble Space Telescope imaging. This approach allowed us to robustly estimate $H_0$, effectively implementing a maximally flexible mass model. Following a blind analysis, we estimated the angular diameter distance to the lens galaxy $D_{\rm d} = 865_{-81}^{+85}$ Mpc and the time-delay distance $D_{\Delta t} = 2180_{-271}^{+472}$ Mpc, giving $H_0 = 77.1_{-7.1}^{+7.3}$ km s$^{-1}$ Mpc$^{-1}$ - for a flat $\Lambda$ cold dark matter cosmology. The error budget accounts for all uncertainties, including the MSD inherent to the lens mass profile and the line-of-sight effects, and those related to the mass-anisotropy degeneracy and projection effects. Our new measurement is in excellent agreement with those obtained in the past using standard simply parametrized mass profiles for this single system ($H_0 = 78.3^{+3.4}_{-3.3}$ km s$^{-1}$ Mpc$^{-1}$) and for seven lenses ($H_0 = 74.2_{-1.6}^{+1.6}$ km s$^{-1}$ Mpc$^{-1}$), or for seven lenses using single-aperture kinematics and the same maximally flexible models used by us ($H_0 = 73.3^{+5.8}_{-5.8}$ km s$^{-1}$ Mpc$^{-1}$). This agreement corroborates the methodology of time-delay cosmography., Comment: 21 pages, 22 figures, 1 table. Accepted by A&A (this version: language improvements only)

The mm-Wave Interferometric Survey of Dark Object Masses (WISDOM) is probing supermassive black holes (SMBHs) in galaxies across the Hubble sequence via molecular gas dynamics. We present the first WISDOM study of a luminous infrared galaxy with an active galactic nuclei (AGN): Fairall 49. We use new ALMA observations of the CO(2-1) line with a spatial resolution of about 80 pc together with ancillary HST imaging. We reach the following results: (1) The CO kinematics are well described by a regularly rotating gas disk with a radial inflow motion, suggesting weak feedback on the cold gas from both AGN and starburst activity; (2) The dynamically inferred SMBH mass is 1.6 +/- 0.4 (rnd) +/- 0.8 (sys) x 10^8 Msun, assuming that we have accurately subtracted the AGN and starburst light contributions, which have a luminosity of about 10^9 Lsun; (3) The SMBH mass agrees with the SMBH-stellar mass relation but is about 50 times higher than previous estimates from X-ray variability; (4) The dynamically inferred molecular gas mass is 30 times smaller than that inferred from adopting the Galactic CO-to-H_2 conversion factor (X_CO) for thermalised gas, suggesting low values of X_CO; (5) the molecular gas inflow rate increases steadily with radius and may be as high as 5 Msun/yr. This work highlights the potential of using high-resolution CO data to estimate, in addition to SMBH masses, the X_CO factor and gas inflow rates in nearby galaxies., Comment: 19 pages, 13 figures, accepted for publication in MNRAS

I introduce some improvements to the pPXF method, which measures the stellar and gas kinematics, star formation history (SFH) and chemical composition of galaxies. I describe the new optimization algorithm that pPXF uses and the changes I made to fit both spectra and photometry simultaneously. I apply the updated pPXF method to a sample of 3200 galaxies at redshift $0.63\times10^{10}$ M$_\odot$), using spectroscopy from the LEGA-C survey (DR3) and 28-bands photometry from two different sources. I compare the masses from new JAM dynamical models with the pPXF stellar population $M_\ast$ and show the latter are more reliable than previous estimates. I use three different stellar population synthesis (SPS) models in pPXF and both photometric sources. I confirm the main trend of the galaxies' global ages and metallicity $[M/H]$ with stellar velocity dispersion $\sigma_\ast$ (or central density), but I also find that $[M/H]$ depends on age at fixed $\sigma_\ast$. The SFHs reveal a sharp transition from star formation to quenching for galaxies with $\lg(\sigma_\ast/\mathrm{km\, s^{-1}})>2.3$, or average mass density within 1 kpc $\lg(\Sigma_1^{\rm JAM}/\mathrm{M_\odot kpc^{-2}})>9.9$, or with $[M/H]>-0.1$, or with Sersic index $\lg n_{\rm Ser}>0.5$. However, the transition is smoother as a function of $M_\ast$. These results are consistent for two SPS models and both photometric sources, but they differ significantly from the third SPS model, which demonstrates the importance of comparing model assumptions. The pPXF software is available from https://pypi.org/project/ppxf/., Comment: 28 pages, 15 figures, LaTeX. Published in MNRAS. The method described in this paper is included in the pPXF software package at https://pypi.org/project/ppxf/

The ``Spectroscopy and H-band Imaging of Virgo cluster galaxies'' (SHIVir) survey is an optical and near-infrared survey which combines SDSS photometry, deep H-band photometry, and long-slit optical spectroscopy for 190 Virgo cluster galaxies (VCGs) covering all morphological types over the stellar mass range log (M_*/M_Sun) = 7.8-11.5$. We present the spectroscopic sample selection, data reduction, and analysis for this SHIVir sample. We have used and optimised the \texttt{pPXF} routine to extract stellar kinematics from our data. Ultimately, resolved kinematic profiles (rotation curves and velocity dispersion profiles) are available for 133 SHIVir galaxies. A comprehensive database of photometric and kinematic parameters for the SHIVir sample is presented with: grizH magnitudes, effective surface brightnesses, effective and isophotal radii, rotational velocities, velocity dispersions, and stellar and dynamical masses. Parameter distributions highlight some bimodal distributions and possible sample biases. A qualitative study of resolved extended velocity dispersion profiles suggests a link between the so-called ``sigma-drop'' kinematic profile and the presence of rings in lenticular S0 galaxies. Rising dispersion profiles are linked to early-type spirals or dwarf ellipticals for which a rotational component is significant, whereas peaked profiles are tied to featureless giant ellipticals., Comment: 21 pages, 16 figures, accepted for publication in MNRAS (doi: 10.1093/mnras/stac1347). Supplemental material available at https://www.queensu.ca/academia/courteau/links

We present a detailed study of a galaxy merger taking place at $z=1.89$ in the GOODS-S field. Here we analyze Keck/MOSFIRE spectroscopic observations from the MOSFIRE Deep Evolution Field (MOSDEF) survey along with multi-wavelength photometry assembled by the 3D-HST survey. The combined dataset is modeled to infer the past star-formation histories (SFHs) of both merging galaxies. They are found to be massive, with log$_{10}(M_{\ast}/M_{\odot}) > 11$, with a close mass ratio satisfying the typical major-merger definition. Additionally, in the context of delayed-$\tau$ models, GOODS-S 43114 and GOODS-S 43683 have similar SFHs and low star-formation rates (log$_{10}$(SFR(SED)/$M_{\odot}/\rm{yr}^{-1}$) $<$ 1.0) compared to their past averages. The best-fit model SEDs show elevated H$\delta_{\rm{A}}$ values for both galaxies, indicating that their stellar spectra are dominated by A-type stars, and that star formation peaked $\sim0.5-1$ Gyr ago and has recently declined. Additionally, based on SED fitting both merging galaxies turned on and shut off star formation within a few hundred Myr of each other, suggesting that their bursts of star formation may be linked. Combining the SFHs and H$\delta_{\rm{A}}$ results with recent galaxy merger simulations, we infer that these galaxies have recently completed their first pericentric passage and are moving apart. Finally, the relatively low second velocity moment of GOODS-S 43114 given its stellar mass, suggests a disk-like structure. However, including the geometry of the galaxy in the modeling does not completely resolve the discrepancy between the dynamical and stellar masses. Future work is needed to resolve this inconsistency in mass., Comment: 12 pages, 7 figures, 1 table, Accepted for publication in MNRAS

The star formation efficiency (SFE) has been shown to vary across different environments, particularly within galactic starbursts and deep within the bulges of galaxies. Various quenching mechanisms may be responsible, ranging from galactic dynamics to feedback from active galactic nuclei (AGN). Here, we use spatially-resolved observations of warm ionised gas emission lines from the imaging Fourier transform spectrograph SITELLE at the Canada-France-Hawaii Telescope (CFHT) and cold molecular gas (CO(2-1)) from the Atacama Large Millimeter/sub-millimeter Array (ALMA) to study the SFE in the bulge of the AGN-host galaxy NGC 3169. After distinguishing star-forming regions from AGN-ionised regions using emission-line ratio diagnostics, we measure spatially-resolved molecular gas depletion times (\tau_dep = 1/SFE) with a spatial resolution of \approx 100 pc within a galactocentric radius of 1.8 kpc. We identify a star-forming ring located at radii 1.25 \pm 0.6 kpc with an average \tau_dep of 0.3 Gyr. At radii < 0.9 kpc, however, the molecular gas surface densities and depletion times increase with decreasing radius, the latter reaching approximately 2.3 Gyr at a radius \approx 500 pc. Based on analyses of the gas kinematics and comparisons with simulations, we identify AGN feedback, bulge morphology and dynamics as the possible causes of the radial profile of SFE observed in the central region of NGC 3169., Comment: 25 pages, 15 figures, accepted for publication in MNRAS

We investigate the presence of ionised gas outflows in a sample of 141 main-sequence star-forming galaxies at $1.2 10.8$, AGN-dominated galaxies, suggesting that AGNs may be the primary drivers of these gas flows. Surprisingly, at $\log(M_\star/M_{\odot})\leq 9.6$, the observed line profiles are fully consistent with a rotating disc model, indicating that ionised gas outflows in dwarf galaxies might play a negligible role even during the peak of cosmic star-formation activity. Finally, we find that the observed mass loading factor scales with stellar mass as expected from the TNG50 cosmological simulation, but the ionised gas mass accounts for only 2$\%$ of the predicted value. This suggests that either the bulk of the outflowing mass is in other gaseous phases or the current feedback models implemented in cosmological simulations need to be revised., Comment: Accepted for publication in MNRAS, 28 pages, 20 figures

We use high-resolution maps of the molecular interstellar medium (ISM) in the centres of eighty-six nearby galaxies from the millimetre-Wave Interferometric Survey of Dark Object Masses (WISDOM) and Physics at High Angular Resolution in Nearby GalaxieS (PHANGS) surveys to investigate the physical mechanisms setting the morphology of the ISM at molecular cloud scales. We show that early-type galaxies tend to have smooth, regular molecular gas morphologies, while the ISM in spiral galaxy bulges is much more asymmetric and clumpy when observed at the same spatial scales. We quantify these differences using non-parametric morphology measures (Asymmetry, Smoothness and Gini), and compare these measurements with those extracted from idealised galaxy simulations. We show that the morphology of the molecular ISM changes systematically as a function of various large-scale galaxy parameters, including galaxy morphological type, stellar mass, stellar velocity dispersion, effective stellar mass surface density, molecular gas surface density, star formation efficiency and the presence of a bar. We perform a statistical analysis to determine which of these correlated parameters best predicts the morphology of the ISM. We find the effective stellar mass surface (or volume) density to be the strongest predictor of the morphology of the molecular gas, while star formation and bars maybe be important secondary drivers. We find that gas self-gravity is not the dominant process shaping the morphology of the molecular gas in galaxy centres. Instead effects caused by the depth of the potential well such as shear, suppression of stellar spiral density waves and/or inflow affect the ability of the gas to fragment., Comment: 20 pages, 9 figures. Accepted for publication in MNRAS

The Sloan Digital Sky Survey MaNGA program has now obtained integral field spectroscopy for over 10,000 galaxies in the nearby universe. We use the final MaNGA data release DR17 to study the correlation between ionized gas velocity dispersion and galactic star formation rate, finding a tight correlation in which sigma_Ha from galactic HII regions increases significantly from ~ 18-30 km/s broadly in keeping with previous studies. In contrast, sigma_Ha from diffuse ionized gas (DIG) increases more rapidly from 20-60 km/s. Using the statistical power of MaNGA, we investigate these correlations in greater detail using multiple emission lines and determine that the observed correlation of sigma_Ha with local star formation rate surface density is driven primarily by the global relation of increasing velocity dispersion at higher total SFR, as are apparent correlations with stellar mass. Assuming HII region models consistent with our finding that sigma_[O III] < sigma_Ha < sigma_[O I], we estimate the velocity dispersion of the molecular gas in which individual HII regions are embedded, finding values sigma_Mol = 5-30 km/s consistent with ALMA observations in a similar mass range. Finally, we use variations in the relation with inclination and disk azimuthal angle to constrain the velocity dispersion ellipsoid of the ionized gas sigma_z/sigma_r = 0.84 +- 0.03 and sigma_phi/sigma_r = 0.91 +- 0.03, similar to that of young stars in the Galactic disk. Our results are most consistent with theoretical models in which turbulence in modern galactic disks is driven primarily by star formation feedback., Comment: 23 pages, 18 figures. Revised version accepted for publication in ApJ

This paper documents the seventeenth data release (DR17) from the Sloan Digital Sky Surveys; the fifth and final release from the fourth phase (SDSS-IV). DR17 contains the complete release of the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey, which reached its goal of surveying over 10,000 nearby galaxies. The complete release of the MaNGA Stellar Library (MaStar) accompanies this data, providing observations of almost 30,000 stars through the MaNGA instrument during bright time. DR17 also contains the complete release of the Apache Point Observatory Galactic Evolution Experiment 2 (APOGEE-2) survey which publicly releases infra-red spectra of over 650,000 stars. The main sample from the Extended Baryon Oscillation Spectroscopic Survey (eBOSS), as well as the sub-survey Time Domain Spectroscopic Survey (TDSS) data were fully released in DR16. New single-fiber optical spectroscopy released in DR17 is from the SPectroscipic IDentification of ERosita Survey (SPIDERS) sub-survey and the eBOSS-RM program. Along with the primary data sets, DR17 includes 25 new or updated Value Added Catalogs (VACs). This paper concludes the release of SDSS-IV survey data. SDSS continues into its fifth phase with observations already underway for the Milky Way Mapper (MWM), Local Volume Mapper (LVM) and Black Hole Mapper (BHM) surveys., Comment: 40 pages, 8 figures, 6 tables. In press at ApJSS (arxiv v2 corrects some minor typos and updates references)

We calculate the fundamental stellar parameters effective temperature, surface gravity and iron abundance - T$_{\rm eff}$, log g, [Fe/H] - for the final release of the Mapping Nearby Galaxies at APO (MaNGA) Stellar Library (MaStar), containing 59,266 per-visit-spectra for 24,290 unique stars at intermediate resolution ($R\sim1800$) and high S/N (median = 96). We fit theoretical spectra from model atmospheres by both MARCS and BOSZ-ATLAS9 to the observed MaStar spectra, using the full spectral fitting code pPXF. We further employ a Bayesian approach, using a Markov Chain Monte Carlo (MCMC) technique to map the parameter space and obtain uncertainties. Originally in this paper, we cross match MaStar observations with Gaia photometry, which enable us to set reliable priors and identify outliers according to stellar evolution. In parallel to the parameter determination, we calculate corresponding stellar population models to test the reliability of the parameters for each stellar evolutionary phase. We further assess our procedure by determining parameters for standard stars such as the Sun and Vega and by comparing our parameters with those determined in the literature from high-resolution spectroscopy (APOGEE and SEGUE) and from lower-resolution matching template (LAMOST). The comparisons, considering the different methodologies and S/N of the literature surveys, are favourable in all cases. Our final parameter catalogue for MaStar cover the following ranges: $2592 \leq $ T$_{\rm eff} \leq 32983\;$K; $-0.7 \leq $ log g $ \leq 5.4\;$dex; $-2.9 \leq $ [Fe/H] $\leq 1.0\;$dex and will be available with the last SDSS-IV Data Release, in December 2021., Comment: 22 pages, 25 figures. MNRAS in press. Stellar parameters will be made available with the Data Release 17 of the SDSS, planned for December 2021

Supermassive black hole masses (MBH) can dynamically be estimated with various methods and using different kinematic tracers. Different methods have only been cross-checked for a small number of galaxies and often show discrepancies. To understand these discrepancies, detailed cross-comparisons of additional galaxies are needed. We present the first part of our cross-comparison between stellar- and gas-based MBH estimates in the nearby fast-rotating early-type galaxy NGC 6958. The measurements presented here are based on ground-layer adaptive optics-assisted Multi-Unit Spectroscopic Explorer (MUSE) science verification data at around 0.6 arcsec spatial resolution. The spatial resolution is a key ingredient for the measurement and we provide a Gaussian parametrisation of the adaptive optics-assisted point spread function (PSF) for various wavelengths. From the MUSE data, we extracted the stellar kinematics and constructed dynamical models. Using an axisymmetric Schwarzschild technique, we measured an MBH of (3.6+2.7-2.4)\times 10^8 Msun at 3\sigma significance taking kinematical and dynamical systematics (e.g.,radially-varying mass-to-light ratio) into account. We also added a dark halo, but our data does not allow to constrain the dark matter fraction. Adding dark matter with an abundance matching prior results in a 25 per cent more massive black hole. Jeans anisotropic models return MBH of (4.6+2.5-2.7) \times 10^8 Msun and (8.6+0.8-0.8) \times 10^8 Msun at 3\sigma confidence for spherical and cylindrical alignment of the velocity ellipsoid, respectively. In a follow-up study, we will compare the stellar-based MBH with those from cold and warm gas tracers, which will provide additional constraints for the MBH for NGC 6958, and insights into assumptions that lead to potential systematic uncertainty., Comment: Accepted for publication in MNRAS, 21 pages, 18 Figures (Including Appendix). Supplementary information is added at the end of the file (3 pages)

We present a comparison of the nitrogen-to-oxygen ratio (N/O) in 37 high-redshift galaxies at $z\sim$2 taken from the KMOS Lensed Emission Lines and VElocity Review (KLEVER) Survey with a comparison sample of local galaxies, taken from the Sloan Digital Sky Survey (SDSS). The KLEVER sample shows only a mild enrichment in N/O of $+$0.1 dex when compared to local galaxies at a given gas-phase metallicity (O/H), but shows a depletion in N/O of $-$0.36 dex when compared at a fixed stellar mass (M$_*$). We find a strong anti-correlation in local galaxies between N/O and SFR in the M$_*$-N/O plane, similar to the anti-correlation between O/H and SFR found in the mass-metallicity relation (MZR). We use this anti-correlation to construct a fundamental nitrogen relation (FNR), analogous to the fundamental metallicity relation (FMR). We find that KLEVER galaxies are consistent with both the FMR and the FNR. This suggests that the depletion of N/O in high-$z$ galaxies when considered at a fixed M$_*$ is driven by the redshift-evolution of the mass-metallicity relation in combination with a near redshift-invariant N/O-O/H relation. Furthermore, the existence of an fundamental nitrogen relation suggests that the mechanisms governing the fundamental metallicity relation must be probed by not only O/H, but also N/O, suggesting pure-pristine gas inflows are not the primary driver of the FMR, and other properties such as variations in galaxy age and star formation efficiency must be important., Comment: 23 pages, 11 figures, Accepted for publication in Monthly Notices of the Royal Astronomical Society

The stellar mass-to-light ratio gradient in SDSS $r-$band $\nabla (M_*/L_r)$ of a galaxy depends on its mass assembly history, which is imprinted in its morphology and gradients of age, metallicity, and stellar initial mass function (IMF). Taking a MaNGA sample of 2051 galaxies with stellar masses ranging from $10^9$ to $10^{12}M_\odot$ released in SDSS DR15, we focus on face-on galaxies, without merger and bar signatures, and investigate the dependence of the 2D $\nabla (M_*/L_r)$ on other galaxy properties, including $M_*/L_r$-colour relationships by assuming a fixed Salpeter IMF as the mass normalization reference. The median gradient is $\nabla M_*/L_r\sim -0.1$ (i.e., the $M_*/L_r$ is larger at the centre) for massive galaxies, becomes flat around $M_*\sim 10^{10} M_{\odot}$ and change sign to $\nabla M_*/L_r\sim 0.1$ at the lowest masses. The $M_*/L_r$ inside a half light radius increases with increasing galaxy stellar mass; in each mass bin, early-type galaxies have the highest value, while pure-disk late-type galaxies have the smallest. Correlation analyses suggest that the mass-weighted stellar age is the dominant parameter influencing the $M_*/L_r$ profile, since a luminosity-weighted age is easily affected by star formation when the specific star formation rate (sSFR) inside the half light radius is higher than $10^{-3} {\rm Gyr}^{-1}$. With increased sSFR gradient, one can obtain a steeper negative $\nabla (M_*/L_r)$. The scatter in the slopes of $M_*/L$-colour relations increases with increasing sSFR, for example, the slope for post-starburst galaxies can be flattened to $0.45$ from the global value $0.87$ in the $M_*/L$ vs. $g-r$ diagram. Hence converting galaxy colours to $M_*/L$ should be done carefully, especially for those galaxies with young luminosity-weighted stellar ages, which can have quite different star formation histories., Comment: 12 pages, 10 figures. Accepted for publication in MNRAS

We present the integral-field kinematics and stellar population properties of 64 galaxies (61 are Early-Type galaxies, ETGs) with Counter-Rotating stellar Disks (CRD) selected from about 4000 galaxies in the MaNGA survey, based on evidence of counter-rotation or two velocity dispersion peaks in the kinematic maps. For 17 CRDs, the counter-rotating components can also be separated spectroscopically. The frequency of CRDs in MaNGA is <5% for ellipticals, <3% for lenticulars and <1% for spirals (at 95% confidence level), consistent with previous estimates. We produced age and metallicity maps, and compared the stellar population properties to those of the general ETGs population. We found that CRDs have similar trends in age and metallicity to ETGs, but are less metallic at low masses, and show flatter age and steeper metallicity gradients, on average. A comparison of the velocity fields of the ionized gas and the stars reveals that in 33 cases the gas corotates with either the inner (15 cases) or outer (18 cases) stellar disk, in 9 cases it is misaligned. In most cases the gas corotates with the younger disk. Evidence of multimodality in the stellar population is found in 31 galaxies, while the 14 youngest and least massive galaxies show ongoing star formation; 14 galaxies, instead, exhibit unimodality, and are the oldest and most massive. As a general result, our work indicates that CRDs form primarily via gas accretion in retrograde rotation with respect to a pre-existing stellar disk., Comment: 26 pages, 22 figures, 1 table. Accepted for publication on December 16, 2021; published as MNRAS 511, 139 (2022)

We present parsec-scale kinematics of eleven nearby galactic nuclei, derived from adaptive-optics assisted integral-field spectroscopy at (near-infrared) CO band-head wavelengths. We focus our analysis on the balance between ordered rotation and random motions, which can provide insights into the dominant formation mechanism of nuclear star clusters (NSCs). We divide our target sample into late- and early-type galaxies, and discuss the nuclear kinematics of the two sub-samples, aiming at probing any link between NSC formation and host galaxy evolution. The results suggest that the dominant formation mechanism of NSCs is indeed affected by the different evolutionary paths of their hosts across the Hubble sequence. More specifically, nuclear regions in late-type galaxies are on average more rotation dominated, and the formation of nuclear stellar structures is potentially linked to the presence of gas funnelled to the center. Early-type galaxies, in contrast, tend to display slower-rotating NSCs with lower ellipticity. However, some exceptions suggest that in specific cases, early-type hosts can form NSCs in a way similar to spirals., Comment: 46 pages, 21 figures, 3 tables, published in ApJ

We construct a dynamical model of the Milky Way disk from a data set, which combines Gaia EDR3 and APOGEE data throughout Galactocentric radii between $5.0\leq R\leq19.5$ kpc. We make use of the spherically-aligned Jeans Anisotropic Method to model the stellar velocities and their velocity dispersions. Building upon our previous work, our model now is fitted to kinematic maps that have been extended to larger Galactocentric radii due to the expansion of our data set, probing the outer regions of the Galactic disk. Our best-fitting dynamical model suggests a logarithmic density slope of $\alpha_{\rm DM}=-1.602\pm0.079_{\rm syst}$ for the dark matter halo and a dark matter density of $\rho_{\rm DM}(R_{\odot})=(8.92\pm0.56_{\rm syst})\times 10^{-3}$ M$_{\odot}$ pc$^{-3}$ ($0.339\pm0.022_{\rm syst}$ GeV cm$^{3}$). We estimate a circular velocity at the solar radius of $v_{\rm circ}=(234.7\pm1.7_{\rm syst})$ km s$^{-1}$ with a decline towards larger radii. The total mass density is $\rho_{\rm tot}(R_{\odot})$=$(0.0672\pm0.0015_{\rm syst})$ M$_{\odot}$ pc$^{-3}$ with a slope of $\alpha_{\rm tot}$=$-2.367\pm0.047_{\rm syst}$ for $5\leq R\leq19.5$ kpc and the total surface density is $\Sigma(R{_\odot}, |z|\leq$ 1.1 kpc)=$(55.5\pm1.7_{\rm syst})$ M$_{\odot}$ pc$^{-2}$. While the statistical errors are small, the error budget of the derived quantities is dominated by the 3 to 7 times larger systematic uncertainties. These values are consistent with our previous determination, but systematic uncertainties are reduced due to the extended data set covering a larger spatial extent of the Milky Way disk. Furthermore, we test the influence of non-axisymmetric features on our resulting model and analyze how a flaring disk model would change our findings., Comment: 19 pages, 11 figures, accepted for publication in ApJ

We present high-resolution (synthesised beam size 0."088x0."083 or 25x23 pc$^2$) Atacama Large Millimetre/submillimetre Array (ALMA) $^{12}$CO(2-1) line and 236 GHz continuum observations, as well as 5 GHz enhanced Multi-Element Radio Linked Interferometer Network (e-MERLIN) continuum observations, of NGC 0708; the brightest galaxy in the low-mass galaxy cluster Abell 262. The line observations reveal a turbulent, rotating disc of molecular gas in the core of the galaxy, and a high-velocity, blue-shifted feature ~0."4 (~113 pc) from its centre. The sub-millimetre continuum emission peaks at the nucleus, but extends towards this anomalous CO emission feature. No corresponding elongation is found on the same spatial scales at 5 GHz with e-MERLIN. We discuss potential causes for the anomalous blue-shifted emission detected in this source, and conclude that it is most likely to be a low-mass in-falling filament of material condensing from the hot intra-cluster medium via chaotic cold accretion, but it is also possible that it is a jet-driven molecular outflow. We estimate the physical properties this structure has in these two scenarios, and show that either explanation is viable. We suggest future observations with integral field spectrographs will be able to determine the true cause of this anomalous emission, and provide further evidence for interaction between quenched cooling flows and mechanical feedback on both small and large scales in this source., Comment: 16 pages, 10 figures, MNRAS accepted

Supermassive black hole (SMBH) masses can be measured by resolving the dynamical influences of the SMBHs on spatially-resolved tracers of the central potentials. Modern long-baseline interferometers have enabled the use of molecular gas as such a tracer. We present here Atacama Large Millimeter/submillimeter Array observations of the elliptical galaxy NGC 7052 at 0.11 arcseconds (37 pc) resolution in the 12CO(2-1) line and 1.3mm continuum emission. This resolution is sufficient to resolve the region in which the potential is dominated by the SMBH. We forward model these observations, using a multi-Gaussian expansion of a Hubble Space Telescope F814W image and spatially-constant mass-to-light ratio to model the stellar mass distribution. We infer a SMBH mass of $2.5\pm0.3\times10^9\,\mathrm{M_\odot}$ and a stellar I-band mass-to-light ratio of $4.6\pm 0.2\,\mathrm{M_\odot/L_{\odot,I}}$ ($3\sigma$ confidence intervals). This SMBH mass is significantly larger than that derived using ionised gas kinematics, which however appear significantly more kinematically disturbed than the molecular gas. We also show that a central molecular gas deficit is likely to be the result of tidal disruption of molecular gas clouds due to the strong gradient in the central gravitational potential., Comment: 14 pages, 8 figures, accepted for publication in MNRAS

We use the statistical power of the MaNGA integral-field spectroscopic galaxy survey to improve the definition of strong line diagnostic boundaries used to classify gas ionization properties in galaxies. We detect line emission from 3.6 million spaxels distributed across 7400 individual galaxies spanning a wide range of stellar masses, star formation rates, and morphological types, and find that the gas-phase velocity dispersion sigma_HAlpha correlates strongly with traditional optical emission line ratios such as [S II]/HAlpha, [N II]/HAlpha, [O I]/HAlpha, and [O III]/HBeta. Spaxels whose line ratios are most consistent with ionization by galactic HII regions exhibit a narrow range of dynamically cold line of sight velocity distributions (LOSVDs) peaked around 25 km/s corresponding to a galactic thin disk, while those consistent with ionization by active galactic nuclei (AGN) and low-ionization emission-line regions (LI(N)ERs) have significantly broader LOSVDs extending to 200 km/s. Star-forming, AGN, and LI(N)ER regions are additionally well separated from each other in terms of their stellar velocity dispersion, stellar population age, HAlpha equivalent width, and typical radius within a given galaxy. We use our observations to revise the traditional emission line diagnostic classifications so that they reliably identify distinct dynamical samples both in two-dimensional representations of the diagnostic line ratio space and in a multi-dimensional space that accounts for the complex folding of the star forming model surface. By comparing the MaNGA observations to the SDSS single-fiber galaxy sample we note that the latter is systematically biased against young, low metallicity star-forming regions that lie outside of the 3 arcsec fiber footprint., Comment: 23 pages, 16 figures. Revised version accepted for publication in ApJ, updates all results to MaNGA final data release DR17

The SDSS-IV Mapping Nearby Galaxies at APO (MaNGA) program has been operating from 2014-2020, and has now observed a sample of 9,269 galaxies in the low redshift universe (z ~ 0.05) with integral-field spectroscopy. With rest-optical (\lambda\lambda 0.36 - 1.0 um) spectral resolution R ~ 2000 the instrumental spectral line-spread function (LSF) typically has 1sigma width of about 70 km/s, which poses a challenge for the study of the typically 20-30 km/s velocity dispersion of the ionized gas in present-day disk galaxies. In this contribution, we present a major revision of the MaNGA data pipeline architecture, focusing particularly on a variety of factors impacting the effective LSF (e.g., undersampling, spectral rectification, and data cube construction). Through comparison with external assessments of the MaNGA data provided by substantially higher-resolution R ~ 10,000 instruments we demonstrate that the revised MPL-10 pipeline measures the instrumental line spread function sufficiently accurately (<= 0.6% systematic, 2% random around the wavelength of Halpha) that it enables reliable measurements of astrophysical velocity dispersions sigma_Halpha ~ 20 km/s for spaxels with emission lines detected at SNR > 50. Velocity dispersions derived from [O II], Hbeta, [O III], [N II], and [S II] are consistent with those derived from Halpha to within about 2% at sigma_Halpha > 30 km/s. Although the impact of these changes to the estimated LSF will be minimal at velocity dispersions greater than about 100 km/s, scientific results from previous data releases that are based on dispersions far below the instrumental resolution should be reevaulated., Comment: 26 pages, 23 figures. Accepted for publication in AJ

Empirical correlations between the masses of supermassive black holes (SMBHs) and properties of their host galaxies are well-established. Among these is the correlation with the flat rotation velocity of each galaxy measured either at a large radius in its rotation curve or via a spatially-integrated emission line width. We propose here the use of the de-projected integrated CO emission line width as an alternative tracer of this rotation velocity, that has already been shown useful for the Tully-Fisher (luminosity-rotation velocity) relation. We investigate the correlation between CO line widths and SMBH masses for two samples of galaxies with dynamical SMBH mass measurements, with respectively spatially-resolved and unresolved CO observations. The tightest correlation is found using the resolved sample of 24 galaxies as $\log (M_\mathrm{BH}/\mathrm{M_\odot})=(7.5\pm0.1)+(8.5\pm0.9)[\log(W_\mathrm{50}/\sin i \,\mathrm{km\,s}^{-1})-2.7]$, where $M_\mathrm{BH}$ is the central SMBH mass, $W_{50}$ the full-width at half-maximum of a double-horned emission line profile, and $i$ the inclination of the CO disc. This relation has a total scatter of $0.6\,$dex, comparable to those of other SMBH mass correlations, and dominated by the intrinsic scatter of $0.5\,$dex. A tight correlation is also found between the de-projected CO line widths and the stellar velocity dispersions averaged within one effective radius. We apply our correlation to the COLD GASS sample to estimate the local SMBH mass function., Comment: 20 pages, 8 figures, accepted for publication in MNRAS. Supplemental material (17 pages, 4 figures) included. Underlying data are available via Zenodo at https://dx.doi.org/10.5281/zenodo.4067034

In the $(\lambda_{\rm R}, \varepsilon)$ and $(V/\sigma, \varepsilon)$ diagrams for characterizing dynamical states, the fast-rotator galaxies (both early-type and spirals) are distributed within a well-defined leaf-shaped envelope. This was explained as due to an upper limit to the orbital anisotropy increasing with galaxy intrinsic flattening. However, a physical explanation for this empirical trend was missing. Here we construct Jeans Anisotropic Models (JAM), with either cylindrically or spherically aligned velocity ellipsoid (two extreme assumptions), and each with either spatially-constant or -variable anisotropy. We use JAM to build mock samples of axisymmetric galaxies, assuming on average an oblate shape for the velocity ellipsoid (as required to reproduce the rotation of real galaxies), and limiting the radial anisotropy $\beta$ to the range allowed by physical solutions. We find that all four mock samples naturally predict the observed galaxy distribution on the $(\lambda_{\rm R}, \varepsilon)$ and $(V/\sigma, \varepsilon)$ diagrams, without further assumptions. Given the similarity of the results from quite different models, we conclude that the empirical anisotropy upper limit in real galaxies, and the corresponding observed distributions in the $(\lambda_{\rm R}, \varepsilon)$ and $(V/\sigma, \varepsilon)$ diagrams, are due to the lack of physical axisymmetric equilibrium solutions at high $\beta$ anisotropy when the velocity ellipsoid is close to oblate., Comment: 5 pages, 4 figures, accepted for publication on MNRAS Letters