Exploring the SMEFT landscape: Bayesian Model Selection for indirect discovery

TL;DR

This paper introduces a Bayesian model selection framework for SMEFT, enabling structured hypothesis testing over operator subsets using advanced algorithms and applying it to current experimental data.

Contribution

It presents a novel Bayesian approach to SMEFT analysis that treats operator subsets as hypotheses, improving the interpretation of experimental constraints.

Findings

No significant evidence for deviations from the Standard Model.

Bayesian model averaging enhances Wilson coefficient characterization.

Operator correlation matrix reveals relational structure and flat directions.

Abstract

We develop a framework for indirect discovery in the Standard Model Effective Field Theory (SMEFT) based on Bayesian model selection over operator subsets. We argue that SMEFT should be understood as a structured space of competing hypotheses rather than a single high-dimensional model, with each operator subset corresponding to a physically distinct low-energy realisation of new dynamics. Bayesian inference is applied at the level of model space itself, assigning posterior probabilities to operator subsets and marginal inclusion probabilities to individual operators. A genetic algorithm efficiently navigates the high-dimensional discrete model space, concentrating evaluations in the high-posterior region, while the Bayesian Information Criterion provides a tractable approximation to the Bayesian evidence. We apply this framework to a dataset comprising electroweak precision observables from LEP and Higgs, top-quark, and diboson measurements from LHC Run 2, at both linear and quadratic order in the Wilson coefficients, with one-loop renormalisation group evolution systematically included. The analysis finds no statistically significant evidence for any departure from the SM, and demonstrates that Bayesian Model Average posteriors on Wilson coefficients carry substantially improved characterisation potential compared to traditional global fits. The operator correlation matrix encodes the relational structure of the model posterior, identifying operator pairs that co-appear in high-posterior models and flat directions where additional measurements would be most valuable. The sensitivity of all results to the choice of matching scale $\mu_0$ is assessed, and its promotion to a continuous parameter of inference is identified as a natural extension of the framework.