Scalable and sequential inference of the neutron star equation of state with the Einstein Telescope

Future gravitational-wave observatories such as the Einstein Telescope will detect tens of thousands of binary neutron star mergers per year, making gravitational waves a dominant probe of the neutron-star equation of state in the coming decades. However, extracting this information requires hierarchical inference across a rapidly growing catalog of events, and existing methods must restart from scratch whenever a new event arrives, making them impractical at this scale. In this Letter, we introduce a hybrid sequential Monte Carlo algorithm that combines data and likelihood tempering to adaptively update the equation-of-state posterior in batches, reusing previous inference results instead of restarting from a wide prior. Accelerated by GPU hardware, our method infers the equation of state from a simulated month of around 1500 binary neutron star mergers in a few hours on a single GPU, and we project that a full year of detections could be processed in a few days on current hardware. This establishes sequential Monte Carlo as a practical route for real-time hierarchical inference with future gravitational-wave detectors.