A quantitative understanding of cometary outbursts requires robust constraints on the size distribution of ejected particles, which governs outburst dynamics and underpins estimates of released gas and dust. In the absence of direct measurements of particle sizes, assumptions about the size distribution play a central role in modelling dust-trail formation, their dynamical evolution and observability, and the potential production of meteor showers following encounters with Earth. We analyse brightness amplitude variations associated with outbursts of comet 17P/Holmes from 1892 to 2021, with particular emphasis on the exceptional 2007 mega-outburst. During this event the comet underwent a rapid and substantial brightening; at its peak, the expanding coma reached a diameter larger than that of the Sun and briefly became the largest object in the Solar System visible to the naked eye. We constrain the size distribution and total mass of porous agglomerates composed of ice, organics, and dust ejected during the outburst. The inferred particle size distribution is consistent with a power law of index q, yielding effective particle sizes between 1.15 x 10^-6 m for q = 4 and 5 x 10^-3 m for q = 2. Accounting for effective particle size, sublimation flux, and bulk density, we find that the total number of ejected particles increases with both q and sublimation flux. These results place quantitative constraints on the physical properties of outburst ejecta and provide physically motivated initial conditions for long-term dust-trail evolution modelling, relevant to the origin of meteoroid streams and the interplanetary dust population.
We report the discovery of bimodal structure in the drift rate distribution of upward-drifting burst clusters from the hyperactive repeating fast radio burst FRB 20240114A. Using unsupervised machine learning (UMAP dimensionality reduction combined with HDBSCAN density-based clustering) applied to 233 upward-drifting burst clusters from the FAST telescope dataset, we identify a distinct subpopulation of 45 burst clusters (Cluster C1) with mean drift rates 2.5x higher than typical upward-drifting burst clusters (245.6 vs 98.1 MHz/ms). Gaussian mixture modeling reveals strong evidence for bimodality (delta-BIC = 296.6), with clearly separated modes (Ashman's D = 2.70 > 2) and a statistically significant gap in the distribution (11.3 sigma). Crucially, we demonstrate that this bimodality persists when restricting the analysis to single-component (U1) burst clusters only (delta-BIC = 19.9, Ashman's D = 2.71), confirming that the result is not an artifact of combining single- and multi-component burst clusters with different drift rate definitions. The extreme-drift subpopulation also exhibits systematically lower peak frequencies (-7%), shorter durations (-29%), and distinct clustering in multi-dimensional feature space. These findings are suggestive of two spatially separated emission regions in the magnetosphere, each producing upward-drifting burst clusters with distinct physical characteristics, although confirmation requires observations from additional epochs and sources.
We present Active Learning for Accelerated Bayesian Inference (\texttt{alabi}): an open-source Python package for performing Bayesian inference with computationally expensive models. Given a forward model and observational data to construct a likelihood and priors, \texttt{alabi}\ uses a Gaussian Process (GP) surrogate model trained to predict posterior probability as a function of input parameters, and employs active learning to iteratively improve GP predictive performance in high-likelihood regions where the GP is most uncertain. \texttt{alabi}\ provides a uniform interface for using Markov chain Monte Carlo (MCMC) with different packages, including the affine-invariant sampler \texttt{emcee}, and nested samplers \texttt{dynesty}, \texttt{multinest}, and \texttt{ultranest}. This approach facilitates accurate estimation of the desired posterior distribution, while reducing the number of computationally expensive model evaluations required by factors of thousands. We demonstrate the performance of \texttt{alabi}\ on a variety of test cases, including where inference is challenging due to complex posterior structure or high dimensionality. We show that \texttt{alabi}\ offers a substantial improvement for likelihood functions with evaluation times $\gtrsim 1$\,s, speeding up MCMC computations by a factor of $10-1000\times$ when tested on problems with up to 64 dimensions.