Zel'dovich smearing approximation of the BAO feature for model-agnostic cosmological inference

A model-agnostic description of the baryon acoustic oscillation (BAO) feature in redshift space requires a number of ingredients. Physically, one must describe the impact of cosmological bulk flows which progressively and anisotropically smear out the feature over time. One must also model the effects of the scale dependence of tracer bias and the mode coupling between short and long scales. All of these can be incorporated using the Zel'dovich approximation alone, without reference to any particular cosmological model. On the technical front, one needs a robust, complete and cosmology-independent basis to describe the shape of the real space BAO feature in linear theory, which can then be propagated to the nonlinearly evolved, measured feature in redshift space. In this work, we describe how these ingredients -- which we have systematically constructed in recent work -- come together in an accurate framework capable of describing the BAO-scale pairwise measurements of state-of-the-art galaxy surveys. Using mock observations and $N$-body simulations, we show that our template-free framework can produce unbiased and precise cosmological constraints for samples with realistic levels of nonlinearity. This work represents one of the final steps in constructing a data-ready analysis framework for model-agnostic cosmological inference from the BAO feature.