2025ApJ...993L..29G

We present the Transiting Exoplanet Survey Satellite (TESS) discovery of only the second system of transiting exocomets with a sufficient number of events to measure the size distribution in the RZ Psc system, enabling comparisons with the β Pictoris and solar system size distributions. Twenty-four transits with absorption depths (ADs) of 1%-20% were observed across three TESS sectors of the 20-50 Myr K0V star, detected as part of our TESS survey of extreme debris disks identified by their IR excess. We discover that the ADs (and hence exocomet radii) follow a broken power-law cumulative frequency distribution not previously seen in extrasolar contexts but similar to that observed in solar system Kuiper Belt object sizes, with power-law slopes above and below the break of γ_AD>break​​ = 2.32 ± 0.12 and γ_AD<break = 0.11 ± 0.04, respectively. We derive size distributions of 1-7 km from two independent lines of evidence. We use the RZ Psc exocomet rate to predict exocomet yields for the Early eVolution Explorer (EVE) NASA astrophysics Small Explorer mission concept to obtain simultaneous photometry of 10⁴ young stars in near-ultraviolet, optical, and near-infrared bands. Assuming occurrence rates scaled from RZ Psc, EVE would detect 590 exocomets from ≈70 young systems in the optical band, with ≈120 simultaneous 5σ detections in all three bands. These data would enable grain sizes of 200-700 nm and graphite-olivine compositions of dozens of events to be distinguished at 2.5σ-3σ, as well as a 4σ determination of the accuracy of the Herschel-derived M debris disk fraction.