A determination of $$\alpha _s(m_Z)$$ at $${{\textrm{aN}}}^3{{\textrm{LO}}}_{{\textrm{QCD}}}\otimes {{\textrm{NLO}}}_{{\textrm{QED}}}$$ accuracy from a global PDF analysis

We present a determination of the strong coupling αs(mZ) from a global dataset including both fixed-target and collider data from deep-inelastic scattering and a variety of hadronic processes, with a simultaneous determination of parton distribution functions (PDFs) based on the NNPDF4.0 methodology. This determination is performed at NNLO and approximate N3LO (aN3LO) perturbative QCD accuracy, including QED corrections and a photon PDF up to NLO accuracy. We extract αs using two independent methodologies, both of which take into account the cross-correlation between αs and the PDFs. The two methodologies are validated by closure tests that allow us to detect and remove or correct for several sources of bias, and lead to mutually consistent results. We account for all correlated experimental uncertainties, as well as correlated theoretical uncertainties related to missing higher order perturbative corrections (MHOUs). We study the perturbative convergence of our results and the impact of QED corrections. We assess individual sources of uncertainty, specifically MHOUs and the value of the top quark mass. We provide a detailed appraisal of methodological choices, including the choice of input dataset, the form of solution of evolution equation, the treatment of the experimental covariance matrix, and the details of Monte Carlo data generation. We find αs(mZ)=0.1194-0.0014+0.0007 at aN3LOQCD⊗NLOQED accuracy, consistent with the latest PDG average and with recent lattice results.