The questions of how the bulk of the Universe’s visible mass emerges and how it is manifest in the existence and properties of hadrons are profound, and probe the heart of strongly interacting matter. Paradoxically, the lightest pseudoscalar mesons appear to be key to a further understanding of the emergent mass and structure mechanisms. These mesons, namely, the pion and kaon, are the Nambu–Goldstone boson modes of quantum chromodynamics (QCD). Unravelling their partonic structure and the interplay between emergent and Higgs-boson mass mechanisms is a common goal of three interdependent approaches—continuum QCD phenomenology, lattice-regularised QCD, and the global analysis of parton distributions—linked to experimental measurements of hadron structure. Experimentally, the anticipated electron–ion collider will enable a revolution in our ability to study pion and kaon structures, accessed by scattering from the ‘meson cloud’ of the proton through the Sullivan process. With the goal of enabling a suite of measurements that can address these questions, we examine key reactions that identify the critical detector-system requirements needed to map tagged pion and kaon cross-sections over a wide range of kinematics. The excellent prospects for extracting pion structural, functional, and form-factor data are outlined, and similar prospects for kaon structures are discussed in the context of a worldwide programme. The successful completion of the programme outlined herein will deliver deep, far-reaching insights into the emergence of pions and kaons, their properties, and their role as QCD’s Goldstone boson modes.

Beginning with results for the leading-twist two-particle distribution amplitudes of $$\pi $$π - and K-mesons, each of which exhibits dilation driven by the mechanism responsible for the emergence of hadronic mass, we develop parameter-free predictions for the pointwise behaviour of all $$\pi $$π and K distribution functions (DFs), including glue and sea. The large-x behaviour of each DF meets expectations based on quantum chromodynamics; the valence-quark distributions match extractions from available data, including the pion case when threshold resummation effects are included; and at $$\zeta _5=5.2\,$$ζ 5 = 5.2 GeV, the scale of existing measurements, the light-front momentum of these hadrons is shared as follows: $$\langle x_{\mathrm{valence}} \rangle ^\pi = 0.41(4)$$⟨ x valence ⟩ π = 0.41 ( 4 ) , $$\langle x_{\mathrm{glue}} \rangle ^\pi = 0.45(2)$$⟨ x glue ⟩ π = 0.45 ( 2 ) , $$\langle x_{\mathrm{sea}} \rangle ^\pi = 0.14(2)$$⟨ x sea ⟩ π = 0.14 ( 2 ) ; and $$\langle x_{\mathrm{valence}} \rangle ^K = 0.42(3)$$⟨ x valence ⟩ K = 0.42 ( 3 ) , $$\langle x_{\mathrm{glue}} \rangle ^K = 0.44(2)$$⟨ x glue ⟩ K = 0.44 ( 2 ) , $$\langle x_{\mathrm{sea}} \rangle ^K = 0.14(2)$$⟨ x sea ⟩ K = 0.14 ( 2 ) . The kaon’s glue and sea distributions are similar to those in the pion, although the inclusion of mass-dependent splitting functions introduces some differences on the valence-quark domain. This study should stimulate improved analyses of existing data and motivate new experiments sensitive to all $$\pi $$π and K DFs. With little known empirically about the structure of the Standard Model’s (pseudo-) Nambu-Goldstone modes and analyses of existing, limited data being controversial, it is likely that new generation experiments at upgraded and anticipated facilities will provide the information needed to resolve the puzzles and complete the picture of these complex bound states.

Analyses of the pion valence-quark distribution function (DF), "Equation missing", which explicitly incorporate the behaviour of the pion wave function prescribed by quantum chromodynamics (QCD), predict "Equation missing", $$\beta (\zeta \gtrsim m_p)>2$$β ( ζ ≳ m p ) > 2 , where $$m_p$$m p is the proton mass. Nevertheless, more than forty years after the first experiment to collect data suitable for extracting the $$x\simeq 1$$x ≃ 1 behaviour of "Equation missing", the empirical status remains uncertain because some methods used to fit existing data return a result for "Equation missing" that violates this constraint. Such disagreement entails one of the following conclusions: the analysis concerned is incomplete; not all data being considered are a true expression of qualities intrinsic to the pion; or QCD, as it is currently understood, is not the theory of strong interactions. New, precise data are necessary before a final conclusion is possible. In developing these positions, we exploit a single proposition, viz. there is an effective charge which defines an evolution scheme for parton DFs that is all-orders exact. This proposition has numerous corollaries, which can be used to test the character of any DF, whether fitted or calculated.

A unified set of predictions for pion and kaon elastic electromagnetic and gravitational form factors is obtained using a symmetry-preserving truncation of each relevant quantum field equation. A key part of the study is a description of salient aspects of the dressed graviton + quark vertices. The calculations reveal that each meson’s mass radius is smaller than its charge radius, matching available empirical inferences; and meson core pressures are commensurate with those in neutron stars. The analysis described herein paves the way for a direct calculation of nucleon gravitational form factors.

Supposing only that there is an effective charge which defines an evolution scheme for parton distribution functions (DFs) that is all-orders exact, strict lower and upper bounds on all Mellin moments of the valence-quark DFs of pionlike systems are derived. Exploiting contemporary results from numerical simulations of lattice-regularized quantum chromodynamics (QCD) that are consistent with these bounds, parameter-free predictions for pion valence, glue, and sea DFs are obtained. The form of the valence-quark DF at large values of the light-front momentum fraction is consistent with predictions derived using the QCD-prescribed behavior of the pion wave function.