Fastest or Significant: A Systematic Framework for Validating Global Minimum Variability Timescale Measurements of Gamma-ray Bursts

TL;DR

GRBs exhibit extreme variability, but reliably measuring the Minimum Variability Timescale (MVT) is challenging due to Poisson noise and analysis choices. The authors extend Haar-based MVT estimation with a comprehensive simulation suite across Gaussian, triangular, and Norris pulse shapes to derive an MVT Validation Curve and a practical workflow for robust interpretation. They apply the framework to Fermi-GBM GRBs and reclassify several published MVT values as upper limits, highlighting the need to account for signal-to-noise when constraining emission-region sizes and jet dynamics, since shorter timescales require higher significance to be resolvable. The work provides a standard, openly available toolset for MVT analyses and emphasizes caution in drawing physical inferences from single MVT measurements in complex bursts.

Abstract

The minimum variability timescale (MVT) is a key observable used to probe the central engines of Gamma-Ray Bursts (GRBs) by constraining the emission region size and the outflow Lorentz factor. However, its interpretation is often ambiguous: statistical noise and analysis choices can bias measurements, making it difficult to distinguish genuine source variability from artifacts. Here we perform a comprehensive suite of simulations to establish a quantitative framework for validating Haar-based MVT measurements. We show that in multi--component light curves, the MVT returns the most statistically significant structure in the interval, which is not necessarily the fastest intrinsic timescale, and can therefore converge to intermediate values. Reliability is found to depend jointly on the MVT value and its signal-to-noise ratio ($\mathrm{SNR}_{\mathrm{MVT}}$), with shorter intrinsic timescales requiring proportionally higher $\mathrm{SNR}_{\mathrm{MVT}}$ to be resolved. We use this relation to define an empirical MVT Validation Curve, and provide a practical workflow to classify measurements as robust detections or upper limits. Applying this procedure to a sample of Fermi-GBM bursts shows that several published MVT values are better interpreted as upper limits. These results provide a path toward standardizing MVT analyses and highlight the caution required when inferring physical constraints from a single MVT measurement in complex events.

Figures (15)

• Figure 1: The median Minimum Variability Timescale (MVT) as a function of analysis bin width (BW) for high--SNR Gaussian pulses of varying intrinsic widths ($\sigma$). The plot shows two distinct regimes: a bin--limited regime (red region) at large BW where the MVT is systematically overestimated, and a source--dominated regime at small BW where the MVT converges to a stable plateau that reflects the true $\sigma$. This demonstrates that a BW significantly smaller than the timescale of interest is required for a reliable measurement.
• Figure 2: The median MVT as a function of peak count–rate amplitude for Gaussian pulses with varying intrinsic widths ($\sigma$). Each marker style represents a different $\sigma$, as shown in the legend. At low amplitudes (red region), the MVT is in the noise–dominated regime, yielding highly scattered and systematically overestimated values. As the amplitude increases, the MVT converges toward the true intrinsic width of the pulse.
• Figure 3: The median MVT as a function of peak amplitude, identical to Figure \ref{['fig:gaussian_mvt_vs_amp']}. Here, the data points are color-coded by the measurement success rate, defined as the percentage of the 300 Monte Carlo realizations that yielded a valid MVT. This figure represents the measurement's reliability, showing a clear transition from unreliable (dark blue, low success rate) at low amplitudes to highly reliable (bright yellow, $\approx$100% success rate) at high amplitudes.
• Figure 4: The median MVT as a function of the $\mathrm{SNR}_{\mathrm{MVT}}$ for simulated Gaussian pulses. Color indicates the intrinsic width ($\sigma$) and marker shape indicates the BW. The plot shows that MVT measurements converge from a noise-dominated, overestimated regime at low SNR to the true intrinsic timescale at high SNR. The specific SNR required for this convergence is dependent on $\sigma$, with shorter timescales requiring a higher SNR for a reliable measurement.
• Figure 5: The median MVT as a function of BW for high-amplitude triangular pulses of varying total widths (colors) and rise-to-fall time ratios (marker shapes). The plot confirms the bin-limited behavior at large BWs and convergence to a stable plateau at small BWs.