G. V. C. Allen, S. Pegwal, D. I. L. de Villiers, D. Anstey, K. Artuc, H. T. J. Bevins, G. Bernardi, M. Bucher, S. Carey, J. Cavillot, R. Chiello, A. S. Chu, W. Croukamp, J. Cumner, A. K. Dash, S. Dasgupta, E. de Lera Acedo, J. Dhandha, A. Dragovic, J. A. Ely, A. Fialkov, T. Gessey-Jones, W. J. Handley, C. Kirkham, G. Kulkarny, S. A. K. Leeney, A. Magro, P. Daan Meerburg, S. Mittal, D. Molnar, R. S. Patel, J. H. N. Pattison, C. M. Pieterse, J. R. Pritchard, G. Rajpoot, N. Razavi-Ghods, D. Robins, I. L. V. Roque, A. Saxena, K. H. Scheutwinkel, P. Scott, E. Shen, P. H. Sims, M. Spinelli, J. Zhu
9 min
Abstract
Recent experiments in cosmology, particularly those aimed at detecting the faint, redshifted, global 21 cm hydrogen line (depth < ~200 mK, z > 7.5), have imposed stringent new requirements on radiometer calibration. In this work, we present a framework for circuit modeling and parameter inference to strengthen these calibration pipelines. This new approach enables in situ characterization of otherwise immeasurable systematics using physically motivated models. A combination of frequentist and Bayesian techniques are employed in a pipeline that supports iterative modeling, robust parameter estimation, and detailed uncertainty quantification. The framework is applied to the REACH telescope, where the precise correction of variations in the radio signal paths arising from component aging or environmental effects is critical. Circuit models of REACH's calibration sources are developed, with the goal of predicting source temperature corrections that are conventionally obtained from laboratory measurements. By fitting the models to measured data using a convolutional cost function, a strong agreement with RMS residuals no worse than -37 dB is obtained. However, Bayesian inference reveals that the resulting temperature corrections can have uncertainties on the order of 1 to 2 K, caused by reflection coefficient degeneracies, measurement noise, and errors in the models. To combat this, posteriors obtained from laboratory measurements are employed as updated priors, reducing correction uncertainties down to 75 mK. Ultimately, the framework provides a means of dynamically accounting for drift in system non-idealities over time, addressing the increasing precision demands of global 21 cm radio astronomy.
Sam: Mismatches cause reflections that distort signals and add unwanted noise. The approach uses four numbers to describe how the receiver's own noise bounces around depending on the source connected—uncorrelated noise, cosine part, sine part, and a noise source term. These noise wave parameters let them solve for the true input temperature from power measurements. Losses in cables transform the physical temperature of a calibration load into an effective temperature at the receiver—like how a leaky hose cools hot water before it arrives.
Alex: Right, so the circuit models step in there to predict those cable details from what they can measure.
Sam: Yes. By modeling cables as transmission lines with resistance, inductance, conductance, and capacitance per length—like sections of pipe with different widths and frictions—they fit all 12 source reflections together. This breaks ambiguities and predicts the hidden details, enabling real-time temperature corrections that cut uncertainties from 1-2 kelvin down to 75 millikelvin—a clear match for the faint signal.
Alex: That logic ties it together neatly. So how do they actually tweak the model to match the on-site reflections?
Sam: They start with a trial-and-error process to find the best parameter values that minimize the difference between measured reflections and model predictions. Traditional ways sum up errors across frequencies and sources, but that often gets stuck in poor solutions because reflections near zero dominate less than those near one, and shared cable parts create overlaps. To fix this, they invented a convolutional cost function: it stacks error sizes for magnitude and phase separately across all 12 sources, then blends them smoothly like mixing layered tracks in audio to balance everything automatically.
Alex: Huh—like treating the errors as layered signals to average out biases. And that helps convergence on the shared parameters?
Sam: Yes, it outperforms standard sums, especially for joint fits. They built an open-source tool called ParamRF to run this efficiently in Python. But frequentist methods have limits: they risk local traps, ignore uncertainties, and might bias predictions without quantifying errors—crucial for cosmology.
Alex: So they turn to Bayesian inference for that robustness. Walk me through how it updates beliefs from data.
Sam: Bayesian inference starts with prior beliefs about parameters—like ranges from engineering know-how—multiplies by how well they match data via a likelihood, and gets a posterior distribution of possibilities. The likelihood assumes measurements are true values plus Gaussian noise—bell-shaped scatter independent in real and imaginary parts. For priors, they train normalizing flows: neural networks that reshape simple distributions into complex ones matching lab posterior samples, capturing correlations like how cable length ties to impedance.
Alex: Okay, so flows turn lab uncertainties into smart starting points for field fits. Does that sharpen the predictions?
Sam: Precisely—it shrinks field uncertainties by leveraging those correlations, enabling the temperature corrections we discussed. The paper suggests this combo provides rigorous error bars, vital since frequentist fits alone could mislead on drifts.
Alex: Overall advance makes sense. With all those parameters, how well do the models actually match real on-site reflections across the 12 sources?
Sam: They test this with joint fits on assembled system data—measuring reflections right where the telescope sits. Using the convolutional cost function, the model predictions hug the measurements closely, with root-mean-square errors no worse than minus 37 decibels across all sources. That's strong agreement in circuit terms, like nailing a puzzle where pieces interlock without forcing.
Alex: So it captures the full set without cherry-picking. But with intertwined parameters, aren't there spots where multiple combos fit the data equally?
Sam: Yes, Bayesian fits on lab data reveal degeneracies—situations where parameters trade off without changing reflections much. For a cable, length and its dielectric constant, which affects signal speed like insulation thickness on a wire, curve together in linked ways. Reflection data alone biases estimates, while lab priors avoid those traps.
Alex: Huh, so reflections don't fully pin down cable traits by themselves. How do they decide which model details to keep or trim?
Sam: They compare models using Bayesian evidence—a score favoring simpler ones that still explain data well, like picking the shortest clear path over a winding detour. For loads, tests show a single transmission line plus key parasitics—like extra capacitors and inductors mimicking tiny leaks—works well. They drop some extras, shrinking total parameters from 123 to 104 without losing fit quality.
Alex: Streamlining to dodge overcomplication. And the payoff for temperatures?
Sam: In simulations mimicking drifts, like boosting cable losses five percent, frequentist fits bias effective source temperatures by tens to hundreds of millikelvin. Bayesian ones with lab-trained priors narrow that to match the faint signal's scale, quantifying uncertainties properly—the paper suggests this rigor is key for real cosmic data.
Alex: So putting it all together, those temperature predictions—what do the simulations and real measurements actually show?
Sam: In simulations of drifts, even a perfect model with basic datasheet starting points leads to hundreds of millikelvin errors from noise and trade-offs. But using lab-updated priors shrinks that to tens of millikelvin. For the hot source measured on site, datasheet priors give 1 to 2 kelvin uncertainties, too wide for the signal. Updated lab priors cut that to around 75 millikelvin.
Alex: A clear drop. Biases meaning the predictions drift off true values?
Sam: Yes—systematic offsets the paper ties to cable deformations not captured, assuming straight uniform paths instead. This flags priorities for better modeling. It enables real-time tracking of drifts from on-site reflections alone, cutting downtime for remote telescopes like REACH. The paper notes it's the first such demo here, paving for continuous cosmic observations with sub-100 millikelvin precision.
Alex: That's a meaningful step for these experiments. Thanks, Sam—this has been a clear look at tackling instrument drifts head-on. Thanks for listening to ResearchPod.