Michele Cappellari

Abstract:

Abstract Axisymmetric Jeans modelling is widely used to infer galaxy mass profiles from integral-field kinematics, but existing implementations maintain tractability by adopting highly restricted anisotropy prescriptions. I present a new spectral method that solves the axisymmetric Jeans equations as a two-dimensional boundary-value problem. Remarkably, this breaks the traditional trade-off between model flexibility and computational cost, accommodating completely general anisotropy distributions β(r, θ) while executing significantly faster than standard restrictive techniques. The method relies on three key choices: (i) solving for the intrinsic dispersion $\overline{v_r^2}$ rather than the rapidly varying pressure $\nu \overline{v_r^2}$ to improve numerical conditioning; (ii) working in logarithmic radius to efficiently resolve the large dynamic range of galaxies, uniquely matching scale-free (power-law) regimes; and (iii) imposing a Robin outer boundary condition that enforces the correct asymptotic decay on a finite computational domain. Orbit integrations in realistic galaxy potentials motivate spherical alignment of the velocity ellipsoid as a physically plausible default, though the framework easily adapts to other alignments. Validated against exact analytic benchmarks—including new analytic Jeans solutions derived herein—the solver recovers intrinsic second moments with high accuracy, showing radially uniform residuals for power-law tests. In practice, it delivers orders-of-magnitude speed-ups over high-accuracy quadrature schemes and is naturally suited to massive GPU parallelization. Released in the public JamPy package, this enables the routine application of highly general Jeans models to large surveys and the extensive parameter-space exploration required for rigorous uncertainty quantification.

Abstract:

We present a new stellar-dynamical measurement of the supermassive black hole (SMBH) mass in the nearby spiral galaxy NGC 4258 (M106), a critical benchmark for extragalactic mass measurements. We use archival James Webb Space Telescope (JWST) Near-Infrared Spectrograph (NIRSpec) integral field unit data (G235H/F170LP grating) to extract high-resolution two-dimensional stellar kinematics from the CO bandhead absorption features within the central 3″ × 3″. We extract the stellar kinematics after correcting for instrumental artifacts and separating the stellar light from the nonthermal active galactic nucleus (AGN) continuum. We employ Jeans anisotropic models to fit the observed kinematics, exploring a grid of 12 models to systematically test the impact of different assumptions for the point-spread function, stellar mass-to-light ratio profile, and orbital anisotropy. All 12 models provide broadly acceptable fits, albeit with minor differences. The ensemble median and 68% (1σ) bootstrap confidence interval of our 12 models yield a black hole mass of MBH=(4.08−0.33+0.19)×107 M⊙. This paper showcases the utility of using the full model ensemble to robustly account for systematic uncertainties, rather than relying on formal errors from a single preferred model, as has been common practice. Our result is just 5% larger than, and consistent with, the benchmark SMBH mass derived from water-maser dynamics, validating the use of NIRSpec stellar kinematics for robust SMBH mass determination. Our analysis demonstrates JWST’s ability to resolve the SMBH’s sphere of influence and deliver precise dynamical masses, even in the presence of significant AGN continuum emission.

Abstract:

Current spatially resolved kinematic measurements of supermassive black hole (SMBH) masses are largely confined to the local Universe (distances Mpc). We investigate the potential of the Extremely Large Telescope’s (ELT) first-light instruments, MICADO and HARMONI, to extend these dynamical measurements to galaxies at redshift . We select a sample of five bright, massive, quiescent galaxies at these redshifts, adopting their Sérsic profiles, from HST photometry, as their intrinsic surface brightness distributions. Based on these intrinsic models, we generate mock MICADO images using SimCADO and mock HARMONI integral-field spectroscopic data cubes using hsim. The HARMONI simulations utilize input stellar kinematics derived from Jeans Anisotropic Models (JAM). We then process these mock observations: the simulated MICADO images are fitted with Multi-Gaussian Expansion (MGE) to derive stellar mass models, and stellar kinematics are extracted from mock HARMONI cubes with pPXF. Finally, these derived stellar mass models and kinematics are used to constrain JAM dynamical models within a Bayesian framework. Our analysis demonstrates that SMBH masses can be recovered with an accuracy of 10 per cent. We find that MICADO can provide detailed stellar mass models with 1 hour of on-source exposure. HARMONI requires longer minimum integrations for reliable stellar kinematic measurements of SMBHs. The required on-source time scales with apparent brightness, ranging from 5–7.5 hours for galaxies at (F814W, 20–20.5 mag) to 5 hours for galaxies at (F160W, 20.8 mag). These findings highlight the ELT’s capability to push the frontier of SMBH mass measurements to , enabling crucial tests of SMBH-galaxy co-evolution at the top end of the galaxies mass function.

Abstract:

Spatially resolved stellar kinematics has become a key ingredient in time-delay cosmography to break the mass-sheet degeneracy in the mass profile and in turn provide a precise constraint on the Hubble constant and other cosmological parameters. In this paper, we present the first measurements of 2D resolved stellar kinematics for the lens galaxy in the quadruply lensed quasar system łensname using integral field spectroscopy from JWST's Near-Infrared Spectrograph (NIRSpec), marking the first such measurement conducted with JWST. In extracting robust kinematic measurements from this first-of-its-kind dataset, we have made methodological improvements both in the data reduction and kinematic extraction. In our kinematic extraction procedure, we performed joint modeling of the lens galaxy, the quasar, and its host galaxy's contributions in the spectra to deblend the lens galaxy component and robustly constrain its stellar kinematics. Our improved methodological frameworks are released as software pipelines for future use: squirrel , for extracting stellar kinematics, and , for JWST-NIRSpec data reduction. We incorporated additional artifact cleaning beyond the standard JWST pipeline. We compared our measured stellar kinematics from the JWST NIRSpec with previously obtained ground-based measurements from the Keck Cosmic Web Imager integral field unit and find that the two datasets are statistically consistent at a ∼1.1σ confidence level. Our measured kinematics will be used in a future study to improve the precision of the Hubble constant measurement. RegalJumper