2026-08-18 | | Total: 7
Uncertainty decompositions in profile-likelihood fits are commonly reported through nuisance-parameter impacts, although shifting a fitted parameter and fluctuating the observation that constrains it answer different questions. Recently, Pinto et al. provided an explicit construction for uncertainty decomposition based on fluctuating the observations and argued that such a construction allows for a cleaner interpretation of systematic uncertainties in terms of physical sources. We provide a geometric description of the distinction between the two methods by using a coupled pair of fiber bundles. The geometrical approach clarifies the relationship between several methods traditionally used to estimate systematic uncertainties in high-energy physics. By describing the profiling as an information-orthogonal horizontal lift and the Schur complement as the induced metric on the parameter-of-interest manifold, physical-source uncertainties arise by mapping observation fluctuations to score covectors, raising them with the inverse total information, and pushing the resulting estimator covariance forward to the parameters of interest. This construction clarifies why nuisance-parameter impacts do not generally coincide with repeated-experiment source variances.
Localized detector inefficiencies produce incomplete calorimeter data that limit the ability to perform precision measurements. We address this problem in relativistic heavy-ion collisions from a Bayesian perspective using pretrained calorimeter diffusion models as priors to reconstruct the missing signal conditioned on surrounding measurements. In this work, we conduct a systematic comparison of several diffusion-based inpainting algorithms, whose performance is evaluated using Bayesian posterior diagnostics of energy response, spatial bias, and uncertainty calibration. The reconstruction fidelity is also analyzed across collision centralities and masked region sizes. This study establishes a general validation strategy for probabilistic reconstruction of missing detector information.
Track seeding strongly affects both the quality and computational cost of charged-particle reconstruction, yet its many configuration parameters are commonly tuned through expert intuition and repeated trial and error. ACTS reduces this burden with an Optuna Tree-structured Parzen Estimator auto-tuner, but expensive evaluations, a restricted search space, and a scalarized objective can limit evaluation efficiency, exclude promising configurations, and obscure performance trade-offs. We investigate whether Bayesian optimization can address these limitations using ACTS with the Open Data Detector (ODD). Under identical search ranges and a common 100-trial budget, we compare Expected Improvement and Upper Confidence Bound with TPE and random search on the existing eight-parameter problem, extend the best-performing Bayesian method to fifteen parameters, and apply Expected Hypervolume Improvement to optimize efficiency, fake rate, duplicate rate, and runtime without fixed scalar weights. Candidate configurations are evaluated through the full ACTS reconstruction chain and validated on disjoint held-out events. The Bayesian acquisition methods identify strong configurations earlier than TPE, and their advantage persists in held-out validation. Expanding the search further improves performance, while multi-objective optimization reveals competitive non-dominated solutions spanning distinct trade-offs. These results indicate that Bayesian optimization can strengthen ACTS auto-tuning through efficient evaluation, broader parameter searches, and post-hoc expert selection among non-dominated alternatives.
Thermal-aware optimization of multi-die 3D integrated circuits evaluates many designs, each a costly heat-equation solve. Operator-learning surrogates replace this solve with a fast forward pass, ideally trained from physics alone, without labeled data. DeepOHeat-v1 made such surrogates fast and trustworthy, but only on low-contrast geometries. High-contrast multi-die stacks break it in two ways: discontinuous conductivities make the continuous physics loss ill-defined at material interfaces, and ill-conditioning ($κ_2(A_h) \approx 6 \times 10^4$) puts the discretized strong-form loss beyond first-order optimization. We propose DeepOHeat-v2 to overcome both. First, we train on a discretized physics loss that handles the discontinuities natively; its energy form reduces the prediction-space loss-Hessian conditioning from $κ^2$ to $κ$, and a matrix-preconditioned optimizer cuts the mean peak temperature error from over 30 K to 0.55 K. Second, because optimization leaves the training distribution, we propose a self-improving framework: a hotspot trust gate sends flagged placements to a reference solver, and the surrogate incrementally retrains on the refined solutions, keeping an update only when it improves held-out validation error. On a multi-die benchmark, the surrogate-true peak gap on the returned design falls from 1.12 K to 0.11 K, matching a solve-at-every-step optimizer while running $56\times$ faster.
Spatial early warning signals (EWSs) seek evidence of an approaching tipping point from a single snapshot of many interacting elements. Existing theory largely assumes spatial homogeneity, whereas networks introduce systematic differences among nodes that may obscure fluctuation-based warning signals. We develop a mathematical framework for spatial EWSs in stochastic dynamical systems on networks. We find that the expected spatial variance, a popular spatial EWS, decomposes exactly into a structural contribution from heterogeneity in the equilibrium state and a fluctuation contribution determined by the stationary covariance. Near a simple steady-state bifurcation, the potentially divergent covariance concentrates along the critical eigendirection: the left eigenvector determines how strongly noise excites the critical fluctuation, while the right eigenvector determines its spatial pattern. Consequently, the spatial variance has a divergent fluctuation contribution when the limiting critical eigendirection is noise-excited and spatially nonuniform after centering. In contrast, the spatial coefficient of variation generally saturates, while skewness, kurtosis, and Moran's $I$ approach network-dependent limits without a universal warning direction. We also derive results for homogeneous networks, node-wise baseline subtraction as preprocessing, and Hopf bifurcations, for which the limiting distributions are qualitatively different. These results clarify when spatial EWSs provide reliable warnings and why their performance depends on network structure, noise, and preprocessing.
Acoustic mode frequencies in the Sun and Sun-like stars change due to magnetic activity, on time-scales much larger than the star's rotation and much smaller than its evolution. Given the poor S/N of the observed stellar p-modes, it is challenging to measure the changes of individual mode frequencies. Typically, power spectra of different time series segments are cross-correlated to estimate a mean p-mode frequency change, which ends up averaging over the individual mode contributions. We seek to enhance the cross-correlation method, by introducing a novel and computationally cheap method, thus enabling us to disentangle p-mode frequency changes for different spherical harmonic degree $\ell$. Assuming that the inclination angle and rotation rate are already measured, filters are designed, which enable the isolation of $δω_\ell$, frequency changes of modes with a given $\ell$, while preventing bias creeping in from neighbouring modes. Monte-Carlo simulations are performed to quantify uncertainty in the estimation of $δω_\ell$. We validate our method against well-studied solar data (SOHO/VIRGO and BiSON) and demonstrate its applicability to the solar-like Kepler star KIC 8006161.
We consider observables $X$ whose realizations in a sampled dataset are restricted, for example by measurement resolution, to a finite set of distinguishable categories within their possibly infinite theoretical domain. Given the probabilities of observable categories, we study the coarse-grained probability mass of families of possible datasets generated through a forward sampling process. The combinatorial construction induces an intrinsically discrete $p$-value defined directly from the sampling process rather than through additional probabilistic structure on observables. Specifying the sample means of $d+k$ arbitrary functions $g_α(X)$ defines a linear family of datasets whose probability mass is obtained as a weighted sum over integer lattice points contained within the associated polyhedron. To overcome the intractable large-$N$ combinatorics, we derive via saddle-point techniques a density approximating these probability masses in the continuum limit of forward sampling within the multinomial universality class. As a demonstration, we consider conditional sampling, where $d$ structural means are fixed while $k$ means vary over admissible datasets. The information geometry emerging from the saddle-point density, together with the spherical symmetry arising at large $N$ from the intrinsic $p$-value construction, enables efficient computation of the $p$-value in the Laplace approximation via the $χ^2_k$ distribution. The resulting statistic is given by the semi-analytic expression $2N$ times the Kullback-Leibler divergence between the information projections associated with the corresponding structural and observed linear families. These projections can be computed efficiently via standard numerical routines converging for sufficiently well-behaved sample means.