Dynamical Modelling of Galactic Kinematics Using Neural Networks

Chapter in Astrophysics and Space Science Proceedings, Springer Nature Switzerland (2026) 117-123

Authors:

David A Simon, Michele Cappellari, Shude Mao, Jiani Chu, Dandan Xu

TDCOSMO. XXIII. Measurement of the Hubble constant from the doubly lensed quasarHE1104-1805

Astronomy & Astrophysics EDP Sciences (2025)

Authors:

Eric Paic, Frédéric Courbin, D Christopher Fassnacht, Aymeric Galan, Martin Millon, Dominique Sluse, M Devon Williams, Simon Birrer, J Elizabeth Buckley-Geer, Michele Cappellari, Frédéric Dux, Xiang-Yu Huang, Shawn Knabel, Cameron Lemon, J Anowar Shajib, H Sherry Suyu, Tommaso Treu, C Kenneth Wong, Lise Christensen, Veronica Motta, Alessandro Sonnenfeld

Abstract:

Time-delay cosmography leverages strongly lensed quasars to measure the Universe's current expansion rate, _ independently from other methods. The latest TDCOSMO milestone measurement primarily used quadruply lensed quasars for their mass profile constraints. However, doubly lensed quasars, being more abundant and offering precise time delays, could expand the sample by a factor of 5, significantly advancing towards a 1% precision measurement of We present the first TDCOSMO analysis of a doubly imaged source, ̋Eonze, including the measurement of the four necessary ingredients. First, by combining 17 years of data from the SMARTS, Euler, and WFI telescopes, we measured a time delay of 176.3 +11.4 -10.3 days. Second, using MUSE data, we extracted stellar velocity dispersion measurements in three radial bins with 5% to 13% precision. Third, employing F160W HST imaging for lens modelling and marginalising over various modelling choices, we measured the Fermat potential difference between the images. Fourth, using wide-field imaging, we measured the convergence added by objects not included in the lens modelling. By combining these four ingredients, we measured the time delay distance and the angular diameter distance to the deflector, favouring a power-law mass model over a baryonic and dark matter composite model. The measurement was performed blindly to prevent experimenter bias and resulted in a Hubble constant of hc = 64.2^ +5.8 _ -5.0 times łint ̨msmpc, where łint is the internal mass sheet degeneracy parameter. This is in agreement with the TDCOSMO-2025 milestone and its precision for łint=1 is comparable to that obtained with the best-observed quadruply lensed quasars (4-6%). This work is a stepping stone towards a precise measurement of using a large sample of doubly lensed quasars, supplementing the current sample. The next TDCOSMO milestone paper will include this system in its hierarchical analysis, constraining łint and jointly with multiple lenses.

TDCOSMO. XXIV. Measurement of the Hubble constant from the doubly lensed quasar HE1104-1805

(2025)

Authors:

Eric Paic, Frà dà ric Courbin, Christopher D Fassnacht, Aymeric Galan, Martin Millon, Dominique Sluse, Devon M Williams, Simon Birrer, Elizabeth J Buckley-Geer, Michele Cappellari, Frà dà ric Dux, Xiang-Yu Huang, Shawn Knabel, Cameron Lemon, Anowar J Shajib, Sherry H Suyu, Tommaso~Treu, Kenneth C Wong, Lise Christensen, Veronica Motta, Alessandro Sonnenfeld

Multimodal atmospheric characterization of β Pictoris b

Astronomy & Astrophysics EDP Sciences 704 (2025) a325

Authors:

M Ravet, M Bonnefoy, G Chauvin, S Lacour, M Nowak, B Charnay, P Tremblin, D Homeier, C Morley, J Fortney, A Denis, S Petrus, P Palma-Bifani, R Landman, LT Parker, M Houllé, A Chomez, K Worthen, F Kiefer, G-D Marleau, Z Zhang, JL Birkby, F Millour, A-M Lagrange, A Vigan, GPPL Otten, J Shangguan

Abstract:

Context. Characterizations of giant exoplanets such as β Pictoris b (hereafter β Pic b) are now routinely performed with multiple spectrographs and imagers exploring different spectral bandwidths and resolutions, allowing for atmospheric retrieval of spectra with or without the conservation of the planet spectral continuum. The accounting of data multimodality in the analysis could provide a more comprehensive determination of the planets physical and chemical properties and inform on their formation history. Aims. We present the first VLTI observations at R λ ∼4000 of β Pic b obtained for an exoplanet with GRAVITY at such a high resolution. We upgraded the forward modelling code ForMoSA to account for the data multimodality, including low-, medium-, and high-resolution spectroscopy based on both a direct model-data comparison and an analysis of cross-correlation signals. We used the ForMoSA code to refine the constraints on the atmospheric properties of the exoplanet and evaluated the sensitivity of the retrieved values to the input dataset. Methods. We obtained four high-signal-to-noise (S/N ∼ 20) spectra of β Pic b in the K band with GRAVITY at R λ ∼4000 conserving both the pseudo-continuum and the pattern of molecular absorptions. We used ForMoSA with four grids of self-consistent forward models (Exo-REM, ATMO, BT-Settl, and Sonora) to explore different T e ff , log(g), metallicity, C/O, and 12 CO/ 13 CO ratio values. We then combined the GRAVITY spectra with published 1–5 µm photometry (NaCo, VisAO, NICI, and SPHERE), low-to-mediumresolution ( R λ ≤ 700 broadband, 0.9–7 µm) spectra, and echelle spectra covering narrower bandwidths ( R λ ∼ 100 000, 2.1–5.2 µm). Results. Sonora and Exo-REM are statistically preferred among all four models, regardless of the dataset used. Exo-REM predicts T eff  = 1607.45 −6.20 +4.85 K and log(g) = 4.46 −0.04 +0.02 dex when using only the GRAVITY epochs, whereas we have T eff  = 1502.74 −2.14 +2.32 K log(g) = 4.00 ± 0.01 dex when incorporating all available datasets. The inclusion of archival data significantly affects all retrieved posteriors. When using all datasets, C/O mostly remains solar (0.552 −0.002 +0.003 ), while [M/H] reaches super-solar values (0.50 ± 0.01). We report the first tentative constraint on the isotopic ratio log( 12 CO/ 13 CO) = 1.12 −0.08 +0.11 in β Pic b’s atmosphere; however, we note that this detection remains inconclusive due to telluric residuals affecting both the GRAVITY and SINFONI data. Additionally, we estimated the bolometric luminosity as log(L/L ⊙ ) = −4.01 −0.05 +0.04 dex. Using a system age of 23 ± 3 Myr, along with this bolometric luminosity and the constraints on the dynamical mass of β Pic b, we were able to constrain the maximum of heavy element content of the planet to be on the order of 5% (20–80 M Earth ). Conclusions. The joint access to the pseudo-continuum and molecular lines in the K band provided by GRAVITY have a significant impact on the retrieved metallicity, possibly owing to the collision-induced absorption driving the continuum shape of the K band. The echelle spectra do not dominate the final fit with respect to lower resolution data covering a broader portion of the spectral energy distribution and the latter keeps encapsulating more robust information on T eff . Future multimodal frameworks should include a weighting scheme to account for the bandwidth and central wavelength of the observations.

TDCOSMO 2025: Cosmological constraints from strong lensing time delays

Astronomy & Astrophysics EDP Sciences 704 (2025) a63

Authors:

Simon Birrer, Elizabeth J Buckley-Geer, Michele Cappellari, Frédéric Courbin, Frédéric Dux, Christopher D Fassnacht, Joshua A Frieman, Aymeric Galan, Daniel Gilman, Xiang-Yu Huang, Shawn Knabel, Danial Langeroodi, Huan Lin, Martin Millon, Takahiro Morishita, Veronica Motta, Pritom Mozumdar, Eric Paic, Anowar J Shajib, William Sheu, Dominique Sluse, Alessandro Sonnenfeld, Chiara Spiniello, Massimo Stiavelli, Sherry H Suyu, Chin Yi Tan, Tommaso Treu, Lyne Van de Vyvere, Han Wang, Patrick Wells, Devon M Williams, Kenneth C Wong

Abstract:

We present cosmological constraints from eight strongly lensed quasars (hereafter, the TDCOSMO-2025 sample). Building on previous work, our analysis incorporated new deflector stellar velocity dispersions measured from spectra obtained with the James Webb Space Telescope (JWST), the Keck Telescopes, and the Very Large Telescope (VLT), utilizing improved methods. We used integrated JWST stellar kinematics for five lenses, VLT-MUSE for 2, and resolved kinematics from Keck and JWST for RX J1131−1231. We also considered two samples of non-time-delay lenses: 11 from the Sloan Lens ACS (SLACS) sample with Keck-KCWI resolved kinematics; and four from the Strong Lenses in the Legacy Survey (SL2S) sample. We improved our analysis of line-of-sight effects, the surface brightness profile of the lens galaxies, and orbital anisotropy, and corrected for projection effects in the dynamics. Our uncertainties are maximally conservative by accounting for the mass-sheet degeneracy in the deflectors’ mass density profiles. The analysis was blinded to prevent experimenter bias. Our primary result is based on the TDCOSMO-2025 sample, in combination with Ω m constraints from the Pantheon+ Type Ia supernovae (SN) dataset. In the flat Λ cold dark matter (CDM), we find H 0 = 71.6 +3.9 −3.3 km s −1 Mpc −1 . The SLACS and SL2S samples are in excellent agreement with the TDCOSMO-2025 sample, improving the precision on H 0 in flat ΛCDM to 4.6%. Using the Dark Energy Survey SN Year-5 dataset (DES-SN5YR) or DESI-DR2 baryonic acoustic oscillations (BAO) likelihoods instead of Pantheon+ yields very similar results. We also present constraints in the open ΛCDM, w CDM, w 0 w a CDM, and w ϕ CDM cosmologies. The TDCOSMO H 0 inference is robust and consistent across all presented cosmological models, and our cosmological constraints in them agree with those from the BAO and SN.