The JWST Ultramassive Galaxy Sample -- I. Overview, Photometry, and Stellar Kinematics: Selecting Optimal Template Subsets via Non-Negative Matrix Factorization

ArXiv 2609.21979 (2026)

Authors:

Michele Cappellari, Dieu D Nguyen, Susan A Kassin

The JWST Ultramassive Galaxy Sample -- II. Supermassive Black Hole Masses for 8 Extreme Early-Type Galaxies via Jeans Anisotropic Modelling of High-Resolution Integral-Field NIRSpec Stellar Kinematics

ArXiv 2609.2198 (2026)

Authors:

Michele Cappellari, Dieu D Nguyen, Susan A Kassin

Hydrodynamic Simulations of Eccentric Black Hole Encounters in AGN Discs

(2026)

Authors:

Andrew Imai, Henry Whitehead, Bence Kocsis

Eccentric stellar-mass binary black holes: Population, detectability, and waveform analysis in the LISA and LIGO era

Physical Review D American Physical Society (APS) 114:6 (2026) 063036

Authors:

Zeyuan Xuan, Smadar Naoz, Kyle Kremer, Michael L Katz, Bence Kocsis, Erez Michaely

Abstract:

Eccentric binary black holes (BBHs) formed through dynamical interactions can significantly contribute to gravitational wave (GW) detections. In this work, we present a simulated catalog of dynamically formed, stellar-mass BBHs in the local universe, incorporating contributions from the Galactic field (flyby interactions), Galactic nucleus (eccentric Kozai-Lidov evolution), and globular clusters (N-body interactions). Our results predict a wide, highly eccentric BBH population in the Milky Way, with source counts of ∼25, 9, 3, 1, and <1 for a 4-year LISA mission (∼36, 13, 5, 2, and 1 for a 10-year extended mission), for SNR>1, 3, 8, 20, and 50, respectively. Extending this model to cosmological populations, we show that different dynamical channels can produce distinct eccentricity distributions in the LVK band and can contribute hundreds of additional low-SNR mHz sources. Specifically, our model yields a merger rate of Γ∼9 Gpc-3 yr-1 and ∼150 extragalactic mHz BBHs with SNR>1. However, due to the lower mass and weaker GW signals of stellar-mass BBHs, this number declines sharply at higher detection thresholds (e.g., ∼1 for SNR>8). We further highlight the impact of eccentric BBH signals on the LISA global fit, showing that their individual harmonics can be independently detected in the Milky Way, and may mimic circular binaries with systematically biased chirp masses. Lastly, we show that post-Newtonian waveforms converge reliably for eccentric BBHs with masses of ≲103M⊙ in the mHz band. Overall, eccentric BBHs represent a prevalent and promising target for future space-based GW observatories.

TDCOSMO XXV: Measuring H0 to 6.5% precision with quasar strong lensing and maximally flexible mass-sheet dynamics

Physical Review D American Physical Society (APS) 114:6 (2026) 063509

Authors:

William Sheu, Tommaso Treu, Martin Millon, Frédéric Dux, Devon Williams, Shawn Knabel, Simon Birrer, Pritom Mozumdar, Giacomo Queirolo, Anowar J Shajib, Michele Cappellari, Kenneth C Wong, Ildar M Asfandiyarov, Otabek A Burkhonov, Frédéric Courbin, Shuhrat A Ehgamberdiev, Sofía Rojas-Ruiz, Asadulla M Shaymanov, Talat A Akhunov

Abstract:

We present a blind time-delay cosmography measurement of the Hubble-Lemaître constant H0 based on the quadruply imaged quasar SDSSJ1433 + 6007. Our analysis combines deep Hubble Space Telescope imaging, extended time-delay monitoring from the Wendelstein and Maidanak Observatories, and spatially resolved stellar kinematics from the Keck Cosmic Web Imager and Reionization Mapper. We build a robust lens model to reconstruct the mass distribution and high-signal-to-noise kinematic maps to break the mass-sheet degeneracy, explicitly accounting for the lens galaxy’s oblateness, rotation, and anisotropy. Furthermore, we constrain the external convergence (κext) by characterizing the line-of-sight environment using wide-field photometry from the Dark Energy Spectroscopic Instrument (DESI) Legacy Survey data release 10. We incorporate these constraints into our joint lensing and dynamical model, running multiple iterations to estimate random and systematic uncertainties. Accounting for maximal flexibility of the mass-sheet transformation, and assuming a flat Λ cold dark matter cosmology and an Ωm,0 prior from DESI data release 2, we infer H0=73.2-4.7+4.8 km s-1 Mpc-1 (a 6.5% precision) and an internal mass-sheet parameter λint=1.12-0.06+0.05. Notably, λint is 2σ away from unity for this system, highlighting the importance of treating it as a free parameter. Our H0 measurement is consistent with the result from our 2025 milestone paper, and it will be included in our next hierarchical analysis to improve the overall precision. Moving forward, the comprehensive pipeline demonstrated herein establishes a robust framework that can be readily applied to future strongly lensed systems to further refine cosmological constraints.