Extremely Large Telescope

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:

Abstract Ultra-short period lava worlds offer a unique window into the coupled evolution of planetary interior and atmospheres under extreme irradiation. In this study, we investigate the mantle dynamics, nightside volcanism, and volatile outgassing on lava world K2-141 b (1.54 R⊕, 5.31 M⊕) using two-dimensional convection models with tracer-based volatile tracking. Our simulations explore a range of interior configurations, including models with and without plastic yielding, basal versus mixed heating, core cooling, and melt intrusion. In models without plastic yielding (i.e. with a strong lithosphere), we find that mantle upwellings form at the substellar and antistellar points, while downwellings form near the day-night terminators at the boundary between the magma ocean and cold, solid nightside. These downwellings facilitate the recycling of crustal material, representing a form of asymmetric, single-lid tectonics. The resulting magma ocean thickness varies from 200 to 300 km depending on the model parameters, corresponding to about 2-3 % of the planet’s radius. Continuous nightside volcanism produces a basaltic crust and gradually depletes the mantle of volatiles. We find that over a billion years, volcanic eruptions can outgas tens of bars of CO2 and H2O. We show that even relatively large volcanic eruptions on the nightside produce thermal emission signals of no more than 1 ppm, remaining below the current detectability threshold in thermal phase curves. However, for most models, outgassing rates are increased near the day-night terminators and future studies should assess whether such localised outgassing could lead to atmospheric signatures in transmission spectroscopy.

Abstract:

Current spatially resolved kinematic measurements of supermassive black hole (SMBH) masses are largely confined to the local Universe (distances ≲ 100 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 1 ≲ z ≲ 2. 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 %. 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 z ≈ 1 (F814W, 20–20.5 mag) to 5 hours for galaxies at 1 < z ≲ 2 (F160W, 20.8 mag). These findings highlight the ELT’s capability to push the frontier of SMBH mass measurements to z ≈ 2, enabling crucial tests of SMBH-galaxy co-evolution at the top end of the galaxies mass function.