Showing 1–50 of 299 results for author: Cappellari, M

WISDOM Project -- XXI. Giant molecular clouds in the central region of the barred spiral galaxy NGC 613: a steep size -- linewidth relation

Authors: Woorak Choi, Martin Bureau, Lijie Liu, Michele Cappellari, Timothy A. Davis, Jindra Gensior, Fu-Heng Liang, Anan Lu, Sanghyuk Moon, Ilaria Ruffa, Thomas G. Williams, Aeree Chung

Abstract: 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… ▽ More 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 ($σ_{\mathrm{obs,los}}$) and gas mass surface densities are larger. The GMC size -- linewidth relation ($σ_{\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α_{\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. △ Less

Submitted 30 May, 2024; v1 submitted 30 May, 2024; originally announced May 2024.

Comments: 15 pages, 8 figures, accepted for publication in MNRAS. arXiv admin note: text overlap with arXiv:2304.10471

Asymmetric drift in MaNGA: Mass and radially-dependent stratification rates in galaxy disks

Authors: Matthew A. Bershady, Kyle B. Westfall, Shravan Shetty, David R. Law, Michele Cappellari, Niv Drory, Kevin Bundy, Renbin Yan

Abstract: 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 qualitati… ▽ More 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<b_r<0.5 for MaNGA galaxies. Overall b_r increases with galaxy mass, decreases with radius for galaxies above 10.4 dex (M_solar) in stellar mass, but is ~constant with radius at lower mass. The measurement scatter indicates the stratification model is too simple to capture the complexity seen in the data, unsurprising given the many possible astrophysical processes that may lead to stellar population dynamical stratification. Nonetheless, the data show dynamical stratification is broadly present in the galaxy population, with systematic trends in mass and density. The amplitude of the asymmetric drift signal is larger for the MaNGA sample than the MW, and better represented in the mean by what is observed in the disks of M31 and M33. Either typical disks have higher surface-density or, more likely, are dynamically hotter (hence thicker) than the MW. △ Less

Submitted 3 May, 2024; originally announced May 2024.

Comments: 21 pages, 18 figures. Accepted for publication in MNRAS

WISDOM Project -- XIX. Figures of merit for supermassive black hole mass measurements using molecular gas and/or megamaser kinematics

Authors: Hengyue Zhang, Martin Bureau, Mark D. Smith, Michele Cappellari, Timothy A. Davis, Pandora Dominiak, Jacob S. Elford, Fu-Heng Liang, Ilaria Ruffa, Thomas G. Williams

Abstract: 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 t… ▽ More 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}/σ_\mathrm{e}^2$, where $σ_\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. △ Less

Submitted 29 April, 2024; v1 submitted 25 April, 2024; originally announced April 2024.

Comments: 13 pages, 9 figures. Accepted by MNRAS

WISDOM Project -- XXIV. Cross-checking supermassive black hole mass estimates from ALMA CO gas kinematics and SINFONI stellar kinematics in the galaxy NGC 4751

Authors: Pandora Dominiak, Michele Cappellari, Martin Bureau, Timothy A. Davis, Marc Sarzi, Ilaria Ruffa, Satoru Iguchi, Thomas G. Williams, Hengyue Zhang

Abstract: 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 cl… ▽ More 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σ$ 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. △ Less

Submitted 17 April, 2024; originally announced April 2024.

Comments: 15 pages, 13 figures, submitted to MNRAS. arXiv admin note: text overlap with arXiv:2401.16376

SDSS-IV MaNGA: Calibration of astrophysical line-widths in the Hα region using HexPak observations

Authors: Sabyasachi Chattopadhyay, Matthew A. Bershady, David R. Law, Kyle Westfall, Shravan Shetty, Camilo Machuca, Michele Cappellari, Kate H. R. Rubin, Kevin Bundy, Samantha Penny

Abstract: 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 ion… ▽ More 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. △ Less

Submitted 2 February, 2024; originally announced February 2024.

INSPIRE: INvestigating Stellar Population In RElics VI -- The low-mass end slope of the stellar Initial Mass Function and chemical composition

Authors: Michalina Maksymowicz-Maciata, Chiara Spiniello, Ignacio Martín-Navarro, Anna Ferré-Mateu, Davide Bevacqua, Michele Cappellari, Giuseppe D'Ago, Crescenzo Tortora, Magda Arnaboldi, Johanna Hartke, Paolo Saracco, Diana Scognamiglio

Abstract: The INSPIRE project has built the largest sample of ultra-compact massive galaxies (UCMGs) at 0.1<z<0.4 and obtained their star formation histories (SFHs). Due to their preserved very old stellar populations, relics are the perfect systems to constrain the earliest epochs of mass assembly in the Universe and the formation of massive early-type galaxies. The goal of this work is to investigate whet… ▽ More The INSPIRE project has built the largest sample of ultra-compact massive galaxies (UCMGs) at 0.1<z<0.4 and obtained their star formation histories (SFHs). Due to their preserved very old stellar populations, relics are the perfect systems to constrain the earliest epochs of mass assembly in the Universe and the formation of massive early-type galaxies. The goal of this work is to investigate whether a correlation exists between the degree of relicness (DoR), quantifying the fraction of stellar mass formed at z>2, 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. △ Less

Submitted 28 January, 2024; originally announced January 2024.

Comments: 16 pages, 11 figures, submitted to MNRAS

The stellar Fundamental Metallicity Relation: the correlation between stellar mass, star-formation rate and stellar metallicity

Authors: Tobias J. Looser, Francesco D'Eugenio, Joanna M. Piotrowska, Francesco Belfiore, Roberto Maiolino, Michele Cappellari, William M. Baker, Sandro Tacchella

Abstract: We present observational evidence for a stellar Fundamental Metallicity Relation (FMR), a smooth relation between stellar mass, star-formation rate (SFR) and the light-weighted stellar metallicity of galaxies (analogous to the well-established gas-phase FMR). We use the flexible, non-parametric software pPXF to reconstruct simultaneously the star-formation and chemical-enrichment history of a repr… ▽ More We present observational evidence for a stellar Fundamental Metallicity Relation (FMR), a smooth relation between stellar mass, star-formation rate (SFR) and the light-weighted stellar metallicity of galaxies (analogous to the well-established gas-phase FMR). We use the flexible, non-parametric software pPXF to reconstruct simultaneously the star-formation and chemical-enrichment history of a representative sample of galaxies from the local MaNGA survey. We find that (i) the metallicity of individual galaxies increases with cosmic time and (ii) at all stellar masses, the metallicity of galaxies is progressively higher, moving from the star-burst region above the main sequence (MS) towards the passive galaxies below the MS, manifesting the stellar FMR. These findings are in qualitative agreement with theoretical expectations from IllustrisTNG, where we find a mass-weighted stellar FMR. The scatter is reduced when replacing the stellar mass $M_{*}$ with $M_{*}/R_{\rm e}$ (with $R_{\rm e}$ being the effective radius), in agreement with previous results using the velocity dispersion $σ_{\rm e}$, which correlates with $M_{*}/R_{\rm e}$. Our results point to starvation as the main physical process through which galaxies quench, showing that metal-poor gas accretion from the intergalactic/circumgalactic medium -- or the lack thereof -- plays an important role in galaxy evolution by simultaneously shaping both their star-formation and their metallicity evolutions, while outflows play a subordinate role. This interpretation is further supported by the additional finding of a young stellar FMR, tracing only the stellar populations formed in the last 300 Myr. This suggests a tight co-evolution of the chemical composition of both the gaseous interstellar medium and the stellar populations, where the gas-phase FMR is continuously imprinted onto the stars over cosmic times. △ Less

Submitted 16 January, 2024; originally announced January 2024.

WISDOM Project -- XVI. The link between circumnuclear molecular gas reservoirs and active galactic nucleus fuelling

Authors: Jacob S. Elford, Timothy A. Davis, Ilaria Ruffa, Martin Bureau, Michele Cappellari, Jindra Gensior, Satoru Iguchi, Fu-Heng Liang, Lijie Liu, Anan Lu, Thomas G. Williams

Abstract: 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… ▽ More 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. △ Less

Submitted 24 December, 2023; v1 submitted 29 November, 2023; originally announced November 2023.

Comments: 15 pages plus 3 in the appendix, 8 figures plus 1 in the appendix, 3 tables plus 4 in the appendix

WISDOM project -- XVIII. Molecular gas distributions and kinematics of three megamaser galaxies

Authors: Fu-Heng Liang, Mark D. Smith, Martin Bureau, Feng Gao, Timothy A. Davis, Michele Cappellari, Jacob S. Elford, Jenny E. Greene, Satoru Iguchi, Federico Lelli, Anan Lu, Ilaria Ruffa, Thomas G. Williams, Hengyue Zhang

Abstract: 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 megamas… ▽ More 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. △ Less

Submitted 10 December, 2023; v1 submitted 26 November, 2023; originally announced November 2023.

Comments: 17 pages, 5 figures, accepted by MNRAS

MaNGA DynPop -- VI. Matter density slopes from dynamical models of 6000 galaxies versus cosmological simulations: the interplay between baryonic and dark matter

Authors: Shubo Li, Ran Li, Kai Zhu, Shengdong Lu, Michele Cappellari, Shude Mao, Chunxiang Wang, Liang Gao

Abstract: 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… ▽ More 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γ_*$, dark $\overlineγ_{\rm DM}$, and total $\overlineγ_{\rm T}$ mass from the MaNGA DynPop project. The $\overlineγ_{\rm T}$ approaches a constant value of 2.2 for high $σ_{\rm e}$ galaxies, and flattens for lg$(σ_{\rm e}/{\rm km\ s^{-1}})\lesssim2.3$, reaching 1.5 for lg$(σ_{\rm e}/{\rm km\ s^{-1}})\approx1.8$. The total and stellar slopes track each other tightly, with $\overlineγ_{\rm T}\approx\overlineγ_*-0.174$ over the full $σ_{\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γ_*$ and $\overlineγ_{\rm DM}$ do not vary inversely within the $σ_{\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γ_{\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. △ Less

Submitted 1 April, 2024; v1 submitted 20 October, 2023; originally announced October 2023.

Comments: 17 pages, 14 figures; accepted by MNRAS

The WISDOM of power spectra: how the galactic gravitational potential impacts a galaxy's central gas reservoir in simulations and observations

Authors: Jindra Gensior, Timothy A. Davis, Martin Bureau, J. M. Diederik Kruijssen, Michele Cappellari, Ilaria Ruffa, Thomas G. Williams

Abstract: 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 stru… ▽ More 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. △ Less

Submitted 10 October, 2023; originally announced October 2023.

Comments: 14 pages, 7 figures + appendices (7 pages, 9 figures); accepted by MNRAS (October 09th, 2023)

INSPIRE: INvestigating Stellar Population In RElics V. A catalogue of ultra-compact massive galaxies outside the local Universe and their degree of relicness

Authors: C. Spiniello, G. D'Ago, L. Coccato, J. Hartke, C. Tortora, A. Ferré-Mateu, C. Pulsoni, M. Cappellari, M. Maksymowicz-Maciata, M. Arnaboldi, D. Bevacqua, A. Gallazzi, L. K. Hunt, F. La Barbera, I. Martín-Navarro, N. R. Napolitano, M. Radovich, P. Saracco, D. Scognamiglio, M. Spavone, S. Zibetti

Abstract: This paper presents the third data release of the INvestigating Stellar Population In RElics (INSPIRE) project, comprising 52 ultra-compact massive galaxies (UCMGs) observed with the X-Shooter spectrograph. We measure integrated stellar velocity dispersion, [Mg/Fe] abundances, ages, and metallicities for all the INSPIRE objects. We thus infer star formation histories and confirm the existence of a… ▽ More This paper presents the third data release of the INvestigating Stellar Population In RElics (INSPIRE) project, comprising 52 ultra-compact massive galaxies (UCMGs) observed with the X-Shooter spectrograph. We measure integrated stellar velocity dispersion, [Mg/Fe] abundances, ages, and metallicities for all the INSPIRE objects. We thus infer star formation histories and confirm the existence of a degree of relicness (DoR), defined in terms of the fraction of stellar mass formed by $z=2$, the time at which a galaxy has assembled 75\% of its mass, and the final assembly time. Objects with a high DoR assembled their stellar mass at early epochs, while low-DoR objects show a non-negligible fraction of later-formed populations and hence a spread in ages and metallicities. A higher DoR correlates with larger [Mg/Fe], super-solar metallicity, and larger velocity dispersion values. The 52 UMCGs span a large range of DoR from 0.83 to 0.06, with 38 of them having formed more than 75\% of their mass by $z=2$. Of these, nine are extreme relics (DoR$>0.7$), since they formed the totality ($>99\%$) of their stellar mass by redshift $z=2$. The remaining 14 UCMGs cannot be considered relics, as they are characterised by more extended star formation histories. With INSPIRE, we built the first sizeable sample of relics outside the local Universe, up to $z\sim0.4$, increasing the number of confirmed relics by a factor of $>10$, and opening up an important window to explain the mass assembly of massive galaxies in the high-z Universe. △ Less

Submitted 14 November, 2023; v1 submitted 22 September, 2023; originally announced September 2023.

Comments: accepted for publication on MNRAS, 20 pages, 15 figures, 3 tables

MaNGA DynPop -- V. The dark-matter fraction versus stellar velocity dispersion relation and stellar initial mass function variations in galaxies: dynamical models and full spectrum fitting of integral-field spectroscopy

Authors: Shengdong Lu, Kai Zhu, Michele Cappellari, Ran Li, Shude Mao, Dandan Xu

Abstract: 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 stella… ▽ More 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 $σ_{\rm e}$ and low $f_{\rm DM}$ at larger $σ_{\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 $α_{\rm IMF}$ increases with $σ_{\rm e}$, consistent with previous studies that reported a transition from Chabrier-like to Salpeter IMF among galaxies. We show that the $α_{\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 $α_{\rm IMF}$ and age, but no correlations with metallicity. We stack galaxy spectra according to their $α_{\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. △ Less

Submitted 23 April, 2024; v1 submitted 21 September, 2023; originally announced September 2023.

Comments: 20 pages, 12 figures, 2 tables; Accepted for publication in MNRAS on 22 April 2024

WISDOM Project -- XVII. Beam-by-beam Properties of the Molecular Gas in Early-type Galaxies

Authors: Thomas G. Williams, Martin Bureau, Timothy A. Davis, Michele Cappellari, Woorak Choi, Jacob S. Elford, Satoru Iguchi, Jindra Gensior, Fu-Heng Liang, Anan Lu, Ilaria Ruffa, Hengyue Zhang

Abstract: We present a study of the molecular gas of seven early-type galaxies with high angular resolution data obtained as part of the mm-Wave Interferometric Survey of Dark Object Masses (WISDOM) project with the Atacama Large Millimeter/submillimeter Array. Using a fixed spatial scale approach, we study the mass surface density ($Σ$) and velocity dispersion ($σ$) of the molecular gas on spatial scales r… ▽ More 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 ($Σ$) and velocity dispersion ($σ$) 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 $α_\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 $α_\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. △ Less

Submitted 9 August, 2023; originally announced August 2023.

Comments: 32 pages (16 of Appendices), 39 Figures (27 in Appendices). Accepted for publication in MNRAS

A fundamental plane of black hole accretion at millimetre wavelengths

Authors: Ilaria Ruffa, Timothy A. Davis, Jacob S. Elford, Martin Bureau, Michele Cappellari, Jindra Gensior, Daryl Haggard, Satoru Iguchi, Federico Lelli, Fu-Heng Liang, Lijie Liu, Marc Sarzi, Thomas G. Williams, Hengyue Zhang

Abstract: 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 ν, 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 ($σ_{\rm int}$) of $0.40$ dex. The plane is found for a sample of 48 nearby galaxies… ▽ More 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 ν, 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 ($σ_{\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 ν, mm}$ at a highly significant level, although such correlation is less tight than the mm fundamental plane ($σ_{\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. △ Less

Submitted 6 November, 2023; v1 submitted 25 July, 2023; originally announced July 2023.

Comments: 9 pages, 5 figures

Velocity Dispersion $σ_{\rm aper}$ Aperture Corrections as a Function of Galaxy Properties from Integral-field Stellar Kinematics of 10,000 MaNGA Galaxies

Authors: Kai Zhu, Ran Li, Xiaoyue Cao, Shengdong Lu, Michele Cappellari, Shude Mao

Abstract: The second moment of the stellar velocity within the effective radius, denoted by $σ_{\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 $σ_{\rm e}$ directly, instead providing $σ_{\rm aper}$, the second moment of the stellar velocity within a fixed fiber apertu… ▽ More The second moment of the stellar velocity within the effective radius, denoted by $σ_{\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 $σ_{\rm e}$ directly, instead providing $σ_{\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 $σ_{\rm aper}/σ_{\rm e}=(R_{\rm aper}/R_{\rm e})^α$, 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 $α$ on the $r$-band absolute magnitude $M_{\rm r}$, $g-i$ color, and Sersic index $n_{\rm Ser}$, where $α$ 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 $α$ values for different galaxy types, with parameters in the ranges of $-18>M_{\rm r}>-24$, $0.4<g-i<1.6$, and $0<n_{\rm Ser}<8$. A Python script is provided to obtain the correction factors from the lookup table. △ Less

Submitted 23 July, 2023; originally announced July 2023.

Comments: 12 pages, 10 figures, 1 table, published in Research in Astronomy and Astrophysics

Journal ref: Research in Astronomy and Astrophysics, Volume 23, Issue 8, id.085001, 12 pp. (2023)

A complete catalogue of broad-line AGNs and double-peaked emission lines from MaNGA integral-field spectroscopy of 10K galaxies: stellar population of AGNs, supermassive black holes, and dual AGNs

Authors: Youquan Fu, Michele Cappellari, Shude Mao, Shengdong Lu, Kai Zhu, Ran Li

Abstract: 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α$ and [OIII] $λ$5007 lines, (2) 101 galaxies with broad $Hα$ lines but no broad [OIII] $λ$5007 lines, a… ▽ More 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α$ and [OIII] $λ$5007 lines, (2) 101 galaxies with broad $Hα$ lines but no broad [OIII] $λ$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α$ lines. We use the broad-$Hα$ line width and luminosity to estimate masses of black hole in these galaxies, and test the $M_{BH}-σ_{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. △ Less

Submitted 20 September, 2023; v1 submitted 4 May, 2023; originally announced May 2023.

Comments: 21 pages, 17 figures, LaTeX. Accepted by MNRAS

MaNGA DynPop -- IV. Stacked total density profile of galaxy groups and clusters from combining dynamical models of integral-field stellar kinematics and galaxy-galaxy lensing

Authors: Chunxiang Wang, Ran Li, Kai Zhu, Huanyuan Shan, Weiwei Xu, Michele Cappellari, Liang Gao, Nan Li, Shengdong Lu, Shude Mao, Ji Yao, Yushan Xie

Abstract: 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 severa… ▽ More 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 $γ_{\rm gnfw}$= $1.82_{-0.25}^{+0.15}$ and $γ_{\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σ$ level. △ Less

Submitted 18 December, 2023; v1 submitted 23 April, 2023; originally announced April 2023.

Comments: 12 pages, 9 figures, accepted by MNRAS

MaNGA DynPop -- III. Stellar dynamics versus stellar population relations in 6000 early-type and spiral galaxies: Fundamental Plane, mass-to-light ratios, total density slopes, and dark matter fractions

Authors: Kai Zhu, Shengdong Lu, Michele Cappellari, Ran Li, Shude Mao, Liang Gao, Junqiang Ge

Abstract: 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… ▽ More 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 $σ_{\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σ_{\rm e}$, with the $M/L$ increasing with $σ_{\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{γ_{\rm T}}$ and $σ_{\rm e}$: $\overline{γ_{\rm T}}$ become steeper with increasing $σ_{\rm e}$ until $\lg(σ_{\rm e}/{\rm km\,s^{-1}})\approx 2.2$ and then remain constant around $\overline{γ_{\rm T}}\approx -2.2$. The $\overline{γ_{\rm T}}-σ_{\rm e}$ variation is larger for LTGs than ETGs. At fixed $σ_{\rm e}$ the total density profiles steepen with galaxy age and for ETGs. We find generally low dark matter fractions, median $f_{\rm DM}(<R_{\rm e})=8$ per cent, within a sphere of radius $R_{\rm e}$. However, we find that $f_{\rm DM}(<R_{\rm e})$ depends on $σ_{\rm e}$ better than stellar mass: dark matter increases to a median $f_{\rm DM}(<R_{\rm e})=33$ percent for galaxies with $σ_{\rm e}\lesssim100{\rm km\,s^{-1}}$. The increased $f_{\rm DM}(<R_{\rm e})$ at low $σ_{\rm e}$ explains the parabolic $\lg(M/L)_{\rm e}-\lgσ_{\rm e}$ relation. △ Less

Submitted 4 November, 2023; v1 submitted 23 April, 2023; originally announced April 2023.

Comments: 26 pages, 21 figures, 1 table; Accepted by MNRAS

MaNGA DynPop -- II. Global stellar population, gradients, and star-formation histories from integral-field spectroscopy of 10K galaxies: link with galaxy rotation, shape, and total-density gradients

Authors: Shengdong Lu, Kai Zhu, Michele Cappellari, Ran Li, Shude Mao, Dandan Xu

Abstract: 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… ▽ More 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 $σ_{\rm e}$, but also find that younger galaxies are more metal-poor at fixed $σ_{\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 $(σ_{\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 $σ_{\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. △ Less

Submitted 27 September, 2023; v1 submitted 23 April, 2023; originally announced April 2023.

Comments: 26 pages, 20 figures, 1 table. Accepted for publication in MNRAS. The catalogue is available at https://manga-dynpop.github.io/

MaNGA DynPop -- I. Quality-assessed stellar dynamical modelling from integral-field spectroscopy of 10K nearby galaxies: a catalogue of masses, mass-to-light ratios, density profiles and dark matter

Authors: Kai Zhu, Shengdong Lu, Michele Cappellari, Ran Li, Shude Mao, Liang Gao

Abstract: 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 Je… ▽ More 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. △ Less

Submitted 4 June, 2023; v1 submitted 23 April, 2023; originally announced April 2023.

Comments: 30 pages, 15 figures, 4 tables, accepted by MNRAS, the catalogue is available here https://manga-dynpop.github.io/

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 522, Issue 4, pp.6326-6353, (2023)

WISDOM Project -- XV. Giant Molecular Clouds in the Central Region of the Barred Spiral Galaxy NGC 5806

Authors: Woorak Choi, Lijie Liu, Martin Bureau, Michele Cappellari, Timothy A. Davis, Jindra Gensior, Fu-Heng Liang, Anan Lu, Thomas G. Williams, Aeree Chung

Abstract: 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). The… ▽ More 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 ($σ_{\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 ($σ_{\mathrm{obs,los}}\propto R_{\mathrm{c}}^{1.20}$), the clouds are on average only marginally bound (with a mean virial parameter $\langleα_{\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. △ Less

Submitted 21 April, 2023; v1 submitted 20 April, 2023; originally announced April 2023.

Comments: Accepted for publication in MNRAS, 20 pages, 16 figures

WISDOM project -- XIV. SMBH mass in the early-type galaxies NGC0612, NGC1574, and NGC4261 from CO dynamical modelling

Authors: Ilaria Ruffa, Timothy A. Davis, Michele Cappellari, Martin Bureau, Jacob S. Elford, Satoru Iguchi, Federico Lelli, Fu-Heng Liang, Lijie Liu, Anan Lu, Marc Sarzi, Thomas G. Williams

Abstract: 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… ▽ More 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$σ$ 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σ$ 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}-σ_{\star}$ relation. △ Less

Submitted 6 November, 2023; v1 submitted 12 April, 2023; originally announced April 2023.

Comments: Main text: 20 pages, 14 Figures; Appendix: 7 pages, 6 Figures. Accepted for publication in MNRAS on 2023 March 28

Supermassive black hole mass in the massive elliptical galaxy M87 from integral-field stellar dynamics using OASIS and MUSE with adaptive optics: assessing systematic uncertainties

Authors: David A. Simon, Michele Cappellari, Johanna Hartke

Abstract: 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 observatio… ▽ More 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 [\text{random}] \pm1.3 [\text{systematic}]) \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 find that this is the primary cause of the difference between our results and previous work, in addition to smaller contributions due to kinematics and modelling method. 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. △ Less

Submitted 28 November, 2023; v1 submitted 31 March, 2023; originally announced March 2023.

Comments: 21 pages, 17 figures, 4 tables, accepted by MNRAS

MaNGA integral-field stellar kinematics of LoTSS radio galaxies: Luminous radio galaxies tend to be slow rotators

Authors: Xuechen Zheng, Huub Röttgering, Arjen van der Wel, Michele Cappellari

Abstract: 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… ▽ More 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 $λ_{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 $λ_{Re}$, and the $F_{radio}$ at a fixed stellar mass decreases with $λ_{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. △ Less

Submitted 25 March, 2023; originally announced March 2023.

Comments: 10 pages,7 figures. Accepted in A&A

Journal ref: A&A 673, A12 (2023)

The massive relic galaxy NGC 1277 is dark matter deficient. From dynamical models of integral-field stellar kinematics out to five effective radii

Authors: Sébastien Comerón, Ignacio Trujillo, Michele Cappellari, Fernando Buitrago, Luis E. Garduño, Javier Zaragoza-Cardiel, Igor A. Zinchenko, Maritza A. Lara-López, Anna Ferré-Mateu, Sami Dib

Abstract: According to the $Λ$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… ▽ More According to the $Λ$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. △ Less

Submitted 8 June, 2023; v1 submitted 20 March, 2023; originally announced March 2023.

Comments: Accepted for publication in A&A

Journal ref: A&A 675, A143 (2023)

Simulating supermassive black hole mass measurements for a sample of ultra massive galaxies using ELT/HARMONI high spatial resolution integral-field stellar kinematics

Authors: Dieu D. Nguyen, Michele Cappellari, Miguel Pereira-Santaella

Abstract: 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$, $|δ+24^{\circ}|<45$°, $|b|>8$°), includi… ▽ More 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$, $|δ+24^{\circ}|<45$°, $|b|>8$°), 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 ($Δ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. △ Less

Submitted 17 September, 2023; v1 submitted 20 February, 2023; originally announced February 2023.

Comments: Accepted to MNRAS

TDCOSMO. XII. Improved Hubble constant measurement from lensing time delays using spatially resolved stellar kinematics of the lens galaxy

Authors: Anowar J. Shajib, Pritom Mozumdar, Geoff C. -F. Chen, Tommaso Treu, Michele Cappellari, Shawn Knabel, Sherry H. Suyu, Vardha N. Bennert, Joshua A. Frieman, Dominique Sluse, Simon Birrer, Frederic Courbin, Christopher D. Fassnacht, Lizvette Villafaña, Peter R. Williams

Abstract: 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… ▽ More 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_{Δt} = 2180_{-271}^{+472}$ Mpc, giving $H_0 = 77.1_{-7.1}^{+7.3}$ km s$^{-1}$ Mpc$^{-1}$ - for a flat $Λ$ 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. △ Less

Submitted 19 November, 2023; v1 submitted 6 January, 2023; originally announced January 2023.

Comments: 21 pages, 22 figures, 1 table. Accepted by A&A (this version: language improvements only)

Journal ref: A&A 673, A9 (2023)

WISDOM Project -- XIII. Feeding molecular gas to the supermassive black hole in the starburst AGN-host galaxy Fairall 49

Authors: Federico Lelli, Timothy A. Davis, Martin Bureau, Michele Cappellari, Lijie Liu, Ilaria Ruffa, Mark D. Smith, Thomas G. Williams

Abstract: 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 ancil… ▽ More 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. △ Less

Submitted 1 September, 2022; originally announced September 2022.

Comments: 19 pages, 13 figures, accepted for publication in MNRAS

Full spectrum fitting with photometry in pPXF: stellar population versus dynamical masses, non-parametric star formation history and metallicity for 3200 LEGA-C galaxies at redshift z~0.8

Authors: Michele Cappellari

Abstract: 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.6<z<1$ (median $z=0.76$, stel… ▽ More 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.6<z<1$ (median $z=0.76$, stellar mass $M_\ast>3\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 $σ_\ast$ (or central density), but I also find that $[M/H]$ depends on age at fixed $σ_\ast$. The SFHs reveal a sharp transition from star formation to quenching for galaxies with $\lg(σ_\ast/\mathrm{km\, s^{-1}})>2.3$, or average mass density within 1 kpc $\lg(Σ_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/. △ Less

Submitted 7 September, 2023; v1 submitted 31 August, 2022; originally announced August 2022.

Comments: 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: Data Catalogue and Kinematic Profiles

Authors: Nathalie N. -Q. Ouellette, Stéphane Courteau, Jon A. Holtzman, Michael McDonald, Michele Cappellari, Joel C. Roediger, Patrick Côté, Julianne J. Dalcanton, Elena Dalla Bontà, Laura Ferrarese, R. Brent Tully, Connor Stone, Eric W. Peng

Abstract: 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 reductio… ▽ More 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. △ Less

Submitted 24 June, 2022; originally announced June 2022.

Comments: 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

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 514, Issue 2, August 2022, Pages 2356-2375

The MOSDEF Survey: A New View of a Remarkable z=1.89 Merger

Authors: Jordan N. Runco, Alice E. Shapley, Mariska Kriek, Michele Cappellari, Michael W. Topping, Ryan L. Sanders, Vasily I. Kokorev, Sedona H. Price, Naveen A. Reddy, Alison L. Coil, Bahram Mobasher, Brian Siana, Tom Zick, Georgios E. Magdis, Gabriel Brammer, James Aird

Abstract: 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… ▽ More 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-$τ$ 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$δ_{\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$δ_{\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. △ Less

Submitted 6 October, 2022; v1 submitted 10 June, 2022; originally announced June 2022.

Comments: 12 pages, 7 figures, 1 table, Accepted for publication in MNRAS

WISDOM project -- XI. Star Formation Efficiency in the Bulge of the AGN-host Galaxy NGC 3169 with SITELLE and ALMA

Authors: Anan Lu, Hope Boyce, Daryl Haggard, Martin Bureau, Fu-Heng Liang, Lijie Liu, Woorak Choi, Michele Cappellari, Laurent Chemin, Mélanie Chevance, Timothy A. Davis, Laurent Drissen, Jacob S. Elford, Jindra Gensior, J. M. Diederik Kruijssen, Thomas Martin, Etienne Massé, Carmelle Robert, Ilaria Ruffa, Laurie Rousseau-Nepton, Marc Sarzi, Gabriel Savard Thomas G. Williams

Abstract: 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 Fourie… ▽ More 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 (τ_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 τ_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. △ Less

Submitted 7 June, 2022; originally announced June 2022.

Comments: 25 pages, 15 figures, accepted for publication in MNRAS

Being KLEVER at cosmic noon: ionised gas outflows are inconspicuous in low-mass star-forming galaxies but prominent in massive AGN hosts

Authors: Alice Concas, Roberto Maiolino, Mirko Curti, Connor Hayden-Pawson, Michele Cirasuolo, Gareth C. Jones, Amata Mercurio, Francesco Belfiore, Giovanni Cresci, Fergus Cullen, Filippo Mannucci, Alessandro Marconi, Michele Cappellari, Claudia Cicone, Yingjie Peng, Paulina Troncoso

Abstract: We investigate the presence of ionised gas outflows in a sample of 141 main-sequence star-forming galaxies at $1.2<z<2.6$ from the KLEVER (KMOS Lensed Emission Lines and VElocity Review) survey. Our sample covers an exceptionally wide range of stellar masses, $8.1<\log(M_\star/M_{\odot})<11.3$, pushing outflow studies into the dwarf regime thanks to gravitationally lensed objects. We stack optical… ▽ More We investigate the presence of ionised gas outflows in a sample of 141 main-sequence star-forming galaxies at $1.2<z<2.6$ from the KLEVER (KMOS Lensed Emission Lines and VElocity Review) survey. Our sample covers an exceptionally wide range of stellar masses, $8.1<\log(M_\star/M_{\odot})<11.3$, pushing outflow studies into the dwarf regime thanks to gravitationally lensed objects. We stack optical rest-frame emission lines (H$β$, [OIII], H$α$ and [NII]) in different mass bins and seek for tracers of gas outflows by using a novel, physically motivated method that improves over the widely used, simplistic double Gaussian fitting. We compare the observed emission lines with the expectations from a rotating disc (disc+bulge for the most massive galaxies) model, whereby significant deviations are interpreted as a signature of outflows. We find clear evidence for outflows in the most massive, $\log(M_\star/M_{\odot}) > 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. △ Less

Submitted 22 March, 2022; originally announced March 2022.

Comments: Accepted for publication in MNRAS, 28 pages, 20 figures

WISDOM Project -- X. The morphology of the molecular ISM in galaxy centres and its dependence on galaxy structure

Authors: Timothy A. Davis, Jindra Gensior, Martin Bureau, Michele Cappellari, Woorak Choi, Jacob S. Elford, J. M. Diederik Kruijssen, Federico Lelli, Fu-Heng Liang, Lijie Liu, Ilaria Ruffa, Toshiki Saito, Marc Sarzi, Andreas Schruba, Thomas G. Williams

Abstract: 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 galaxi… ▽ More 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. △ Less

Submitted 7 March, 2022; v1 submitted 2 March, 2022; originally announced March 2022.

Comments: 20 pages, 9 figures. Accepted for publication in MNRAS

The Westerbork Coma Survey: A blind, deep, high-resolution HI survey of the Coma cluster

Authors: D. Cs. Molnar, P. Serra, T. van der Hulst, T. H. Jarrett, A. Boselli, L. Cortese, J. Healy, E. de Blok, M. Cappellari, K. M. Hess, G. I. G. Jozsa, R. M. McDermid, T. A. Oosterloo, M. A. W. Verheijen

Abstract: We present the blind Westerbork Coma Survey probing the HI content of the Coma galaxy cluster with the Westerbork Synthesis Radio Telescope. The survey covers the inner $\sim$ 1 Mpc around the cluster centre, extending out to 1.5 Mpc towards the south-western NGC 4839 group. The survey probes the atomic gas in the entire Coma volume down to a sensitivity of $\sim$ 10$^{19}$ cm$^{-2}$ and 10$^8$ M… ▽ More We present the blind Westerbork Coma Survey probing the HI content of the Coma galaxy cluster with the Westerbork Synthesis Radio Telescope. The survey covers the inner $\sim$ 1 Mpc around the cluster centre, extending out to 1.5 Mpc towards the south-western NGC 4839 group. The survey probes the atomic gas in the entire Coma volume down to a sensitivity of $\sim$ 10$^{19}$ cm$^{-2}$ and 10$^8$ M$_{\odot}$. Combining automated source finding with source extraction at optical redshifts and visual verification, we obtained 40 HI detections of which 24 are new. Over half of the sample displays perturbed HI morphologies indicative of an ongoing interaction with the cluster environment. With the use of ancillary UV and mid-IR, data we measured their stellar masses and star formation rates and compared the HI properties to a set of field galaxies spanning a similar stellar mass and star formation rate range. We find that $\sim$ 75 % of HI-selected Coma galaxies have simultaneously enhanced star formation rates (by $\sim$ 0.2 dex) and are HI deficient (by $\sim$ 0.5 dex) compared to field galaxies of the same stellar mass. According to our toy model, the simultaneous HI deficiency and enhanced star formation activity can be attributed to either HI stripping of already highly star forming galaxies on a very short timescale, while their H$_2$ content remains largely unaffected, or to HI stripping coupled to a temporary boost of the HI-to-H$_2$ conversion, causing a brief starburst phase triggered by ram pressure before eventually quenching the galaxy. △ Less

Submitted 6 January, 2022; v1 submitted 22 December, 2021; originally announced December 2021.

Comments: Accepted for publication in A&A, 27 pages, 21 figures

Journal ref: A&A 659, A94 (2022)

SDSS-IV MaNGA: Understanding Ionized Gas Turbulence using Integral Field Spectroscopy of 4500 Star-Forming Disk Galaxies

Authors: David R. Law, Francesco Belfiore, Matthew A. Bershady, Michele Cappellari, Niv Drory, Karen L. Masters, Kyle B. Westfall, Dmitry Bizyaev, Kevin Bundy, Kaike Pan, Renbin Yan

Abstract: 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 i… ▽ More 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. △ Less

Submitted 4 March, 2022; v1 submitted 21 December, 2021; originally announced December 2021.

Comments: 23 pages, 18 figures. Revised version accepted for publication in ApJ

The Seventeenth Data Release of the Sloan Digital Sky Surveys: Complete Release of MaNGA, MaStar and APOGEE-2 Data

Authors: Abdurro'uf, Katherine Accetta, Conny Aerts, Victor Silva Aguirre, Romina Ahumada, Nikhil Ajgaonkar, N. Filiz Ak, Shadab Alam, Carlos Allende Prieto, Andres Almeida, Friedrich Anders, Scott F. Anderson, Brett H. Andrews, Borja Anguiano, Erik Aquino-Ortiz, Alfonso Aragon-Salamanca, Maria Argudo-Fernandez, Metin Ata, Marie Aubert, Vladimir Avila-Reese, Carles Badenes, Rodolfo H. Barba, Kat Barger, Jorge K. Barrera-Ballesteros, Rachael L. Beaton , et al. (316 additional authors not shown)

Abstract: 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… ▽ More 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. △ Less

Submitted 13 January, 2022; v1 submitted 3 December, 2021; originally announced December 2021.

Comments: 40 pages, 8 figures, 6 tables. In press at ApJSS (arxiv v2 corrects some minor typos and updates references)

SDSS-IV MaStar: Theoretical Atmospheric Parameters for the MaNGA Stellar Library

Authors: Lewis Hill, Daniel Thomas, Claudia Maraston, Renbin Yan, Justus Neumann, Andrew Lundgren, Daniel Lazarz, Yan-Ping Chen, Michele Cappellari, Jon A. Holtzman, Julie Imig, Katia Cunha, Guy Stringfellow, Dmitry Bizyaev, David R. Law, Keivan G. Stassun, Niv Drory, Michael Merrifield, Timothy C. Beers

Abstract: 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… ▽ More 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. △ Less

Submitted 9 November, 2021; originally announced November 2021.

Comments: 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

Cross-checking SMBH mass estimates in NGC 6958 -- I: Stellar dynamics from adaptive optics-assisted MUSE observations

Authors: Sabine Thater, Davor Krajnović, Peter M. Weilbacher, Dieu D. Nguyen, Martin Bureau, Michele Cappellari, Timothy A. Davis, Satoru Iguchi, Richard McDermid, Kyoko Onishi, Marc Sarzi, Glenn van de Ven

Abstract: 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- an… ▽ More 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σ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σ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. △ Less

Submitted 2 November, 2021; originally announced November 2021.

Comments: Accepted for publication in MNRAS, 21 pages, 18 Figures (Including Appendix). Supplementary information is added at the end of the file (3 pages)

The MBHBM$^{\star}$ Project -- II. Molecular Gas Kinematics in the Lenticular Galaxy NGC 3593 Reveal a Supermassive Black Hole

Authors: Dieu D. Nguyen, Martin Bureau, Sabine Thater, Kristina Nyland, Mark den Brok, Michelle Cappellari, Timothy A. Davis, Jenny E. Greene, Nadine Neumayer, Masatoshi Imanishi, Takuma Izumi, Taiki Kawamuro, Shunsuke Baba, Phuong M. Nguyen, Satoru Iguchi, Takafumi Tsukui, Lam N. T., Than Ho

Abstract: As part of the Measuring Black Holes in Below Milky Way-mass (M$^\star$) galaxies (MBHBM$^\star$) Project, we present a dynamical measurement of the supermassive black hole (SMBH) mass in the nearby lenticular galaxy NGC 3593, using cold molecular gas $^{12}$CO(2-1) emission observed at an angular resolution of $\approx0''.3$ ($\approx10$ pc) with the Atacama Large Millimeter/submillimeter Array (… ▽ More As part of the Measuring Black Holes in Below Milky Way-mass (M$^\star$) galaxies (MBHBM$^\star$) Project, we present a dynamical measurement of the supermassive black hole (SMBH) mass in the nearby lenticular galaxy NGC 3593, using cold molecular gas $^{12}$CO(2-1) emission observed at an angular resolution of $\approx0''.3$ ($\approx10$ pc) with the Atacama Large Millimeter/submillimeter Array (ALMA). Our ALMA observations reveal a circumnuclear molecular gas disc (CND) elongated along the galaxy major axis and rotating around the SMBH. Using dynamical modelling, the molecular gas kinematics allow us to infer a SMBH mass $M_{\rm BH}=2.40_{-1.05}^{+1.87}\times10^6$ M$_\odot$ (only statistical uncertainties at the $3σ$ level). We also detect a massive core of cold molecular gas (CMC) of mass $M_{\rm CMC}=(5.4\pm1.2)\times10^6$ M$_\odot$ and effective (half-mass) radius $r_{\rm CMC,e}=11.2\pm2.8$ pc, co-spatial with a nuclear star cluster (NSC) of mass $M_{\rm NSC}=(1.67\pm0.48)\times10^7$ M$_\odot$ and effective radius $r_{\rm NSC,e}=5.0\pm1.0$~pc (or $0''.15\pm0''.03$). The mass profiles of the CMC and NSC are well described by Sérsic functions with indices $1-1.4$. Our $M_{\rm BH}$ and $M_{\rm NSC}$ estimates for NGC 3593 agree well with the recently compiled $M_{\rm BH}$-$M_{\rm NSC}$ scaling relation. Although the $M_{\rm NSC}$ uncertainty is twice the inferred $M_{\rm BH}$, the rapid central rise of the rotation velocities of the CND (as the radius decreases) clearly suggests a SMBH. Indeed, our dynamical models show that even if $M_{\rm NSC}$ is at the upper end of its allowed range, the evidence for a black hole does not vanish, but remains with a lower limit of $M_{\rm BH}>3\times10^5$ M$_\odot$. △ Less

Submitted 16 October, 2021; originally announced October 2021.

Comments: 27 pages, 18 Figures, 8 Tables. Accepted to MNRAS

The KLEVER survey: Nitrogen abundances at $z\sim$2 and probing the existence of a fundamental nitrogen relation

Authors: Connor Hayden-Pawson, Mirko Curti, Roberto Maiolino, Michele Cirasuolo, Francesco Belfiore, Michele Cappellari, Alice Concas, Giovanni Cresci, Fergus Cullen, Chiaki Kobayashi, Filippo Mannucci, Alessandro Marconi, Massimo Meneghetti, Amata Mercurio, Yingjie Peng, Mark Swinbank, Fiorenzo Vincenzo

Abstract: 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 metall… ▽ More 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. △ Less

Submitted 9 March, 2022; v1 submitted 30 September, 2021; originally announced October 2021.

Comments: 23 pages, 11 figures, Accepted for publication in Monthly Notices of the Royal Astronomical Society

SDSS-IV MaNGA: Stellar M/L gradients and the M/L-colour relation in galaxies

Authors: Junqiang Ge, Shude Mao, Youjun Lu, Michele Cappellari, Richard J. Long, Renbin Yan

Abstract: 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… ▽ More 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. △ Less

Submitted 11 August, 2021; originally announced August 2021.

Comments: 12 pages, 10 figures. Accepted for publication in MNRAS

SDSS-IV MaNGA: Integral-field kinematics and stellar population of a sample of galaxies with counter-rotating stellar disks selected from about 4000 galaxies

Authors: Davide Bevacqua, Michele Cappellari, Silvia Pellegrini

Abstract: 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… ▽ More 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. △ Less

Submitted 16 February, 2022; v1 submitted 20 July, 2021; originally announced July 2021.

Comments: 26 pages, 22 figures, 1 table. Accepted for publication on December 16, 2021; published as MNRAS 511, 139 (2022)

Journal ref: MNRAS 511, 139 (2022)

Resolved nuclear kinematics link the formation and growth of nuclear star clusters with the evolution of their early and late-type hosts

Authors: Francesca Pinna, Nadine Neumayer, Anil Seth, Eric Emsellem, Dieu D. Nguyen, Torsten Boeker, Michele Cappellari, Richard M. McDermid, Karina Voggel, C. Jakob Walcher

Abstract: 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- a… ▽ More 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. △ Less

Submitted 26 October, 2021; v1 submitted 19 July, 2021; originally announced July 2021.

Comments: 46 pages, 21 figures, 3 tables, published in ApJ

Journal ref: ApJ 921 (2021) 8

Dynamical model of the Milky Way using APOGEE and Gaia data

Authors: Maria Selina Nitschai, Anna-Christina Eilers, Nadine Neumayer, Michele Cappellari, Hans-Walter Rix

Abstract: 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 extend… ▽ More 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 $α_{\rm DM}=-1.602\pm0.079_{\rm syst}$ for the dark matter halo and a dark matter density of $ρ_{\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 $ρ_{\rm tot}(R_{\odot})$=$(0.0672\pm0.0015_{\rm syst})$ M$_{\odot}$ pc$^{-3}$ with a slope of $α_{\rm tot}$=$-2.367\pm0.047_{\rm syst}$ for $5\leq R\leq19.5$ kpc and the total surface density is $Σ(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. △ Less

Submitted 9 June, 2021; originally announced June 2021.

Comments: 19 pages, 11 figures, accepted for publication in ApJ

WISDOM project -- VIII. Multi-scale feedback cycles in the brightest cluster galaxy NGC 0708

Authors: Eve V. North, Timothy A. Davis, Martin Bureau, Massimo Gaspari, Michele Cappellari, Satoru Iguchi, Lijie Liu, Kyoko Onishi, Marc Sarzi, Mark D. Smith, Thomas G. Williams

Abstract: 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 re… ▽ More 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. △ Less

Submitted 16 March, 2021; originally announced March 2021.

Comments: 16 pages, 10 figures, MNRAS accepted

WISDOM Project -- VII. Molecular gas measurement of the supermassive black hole mass in NGC 7052

Authors: Mark D. Smith, Martin Bureau, Timothy A. Davis, Michele Cappellari, Lijie Liu, Kyoko Onishi, Satoru Iguchi, Eve V. North, Marc Sarzi, Thomas G. Williams

Abstract: 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… ▽ More 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σ$ 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. △ Less

Submitted 16 March, 2021; originally announced March 2021.

Comments: 14 pages, 8 figures, accepted for publication in MNRAS

SDSS-IV MaNGA: Refining Strong Line Diagnostic Classifications Using Spatially Resolved Gas Dynamics

Authors: David R. Law, Xihan Ji, Francesco Belfiore, Matthew A. Bershady, Michele Cappellari, Kyle B. Westfall, Renbin Yan, Dmitry Bizyaev, Joel R. Brownstein, Niv Drory, Brett H. Andrews

Abstract: 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 t… ▽ More 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. △ Less

Submitted 3 May, 2021; v1 submitted 11 November, 2020; originally announced November 2020.

Comments: 23 pages, 16 figures. Revised version accepted for publication in ApJ, updates all results to MaNGA final data release DR17

SDSS-IV MaNGA: Modeling the Spectral Line Spread Function to Sub-Percent Accuracy

Authors: David R. Law, Kyle B. Westfall, Matthew A. Bershady, Michele Cappellari, Renbin Yan, Francesco Belfiore, Dmitry Bizyaev, Joel R. Brownstein, Yanping Chen, Brian Cherinka, Niv Drory, Daniel Lazarz, Shravan Shetty

Abstract: 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 (λλ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 challen… ▽ More 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 (λλ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. △ Less

Submitted 13 November, 2020; v1 submitted 9 November, 2020; originally announced November 2020.

Comments: 26 pages, 23 figures. Accepted for publication in AJ