A semi-analytical surrogate model for giant planet evolution: bypassing ordinary differential equation solvers with localised thermodynamics, softplus asymptotes, and B-spline photometry

Abstract

Context. Evolutionary models translate the observable luminosity, temperature, and colours of giant planets and brown dwarfs into mass and age. Generating their cooling tracks normally requires integrating the internal energy over time with an ordinary differential equation (ODE) solver coupled to pre-computed atmospheric grids, which becomes numerically stiff at sharp transitions such as cloud condensation and the onset of electron degeneracy, and is fragile inside Bayesian retrievals. Aims. We…