We establish an approach to analyze the free hadron and transition (nonperturbative) regions of near-side energy-energy correlators (EECs) based on dihadron fragmentation functions (DiFFs). We introduce a (nonperturbative) function we call the "EEC DiFF" and explicitly show that expanding it for large relative transverse momentum between the two hadrons gives the O(α_{s}) expression for the "EEC jet" function used in the quark-gluon (perturbative) region. This connection indicates that a formal theoretical matching will be able to bridge the free-hadron region, transition, and quark-gluon regions and allow all of them to be analyzed simultaneously. We further derive a result valid for near-side EECs in the free hadron and transition regions of e^{+}e^{-} annihilation in terms of the EEC DiFF. Using a simple model for the function, we perform the first fit within the dihadron framework to experimental data in this regime. We find reasonable agreement with the measurements and reproduce the salient features of near-side EECs in the free hadron and transition regions.
Three-body forces, a type of nonperturbative strong interactions, are widely studied in nuclear physics. However, whether their inclusion is necessary in nuclear systems remains a topic of intense debate. In this Letter, we propose that the existence of three-body forces in certain three-body hadronic systems with definite C parity is certain. Such systems consist of two components whose interactions are mediated by three-body forces-a mechanism not easily realized in conventional three-nucleon systems. We investigate two specific three-body hadronic systems, D[over ¯]_{s}DK and D[over ¯]^{*}Dη, using contact-range potentials. The two-body hadron-hadron interactions are constrained by reproducing their scattering lengths, while the three-body couplings are constrained by charge symmetry. Our results indicate that three-body forces play a minor role in binding the I(J^{PC})=0(0^{--})D[over ¯]_{s}DK system, but a crucial one in binding the I(J^{PC})=0(1^{-+})D[over ¯]^{*}Dη system. In fact, three-body forces determine whether D[over ¯]^{*}Dη forms a bound state, making this system a promising candidate for exploring three-body forces in hadronic physics.
Carbon ion radiotherapy (CIRT) combines the physical precision of charged particles with high-linear energy transfer (LET) biological effectiveness, enabling superior control of radioresistant or anatomically complex tumors. Initially developed in Japan, CIRT is slowly expanding globally with increasing institutional and clinical maturity. To summarize current clinical evidence, operational strategies, and expansion plans for carbon ion facilities worldwide, as discussed during the "Hadrontherapy-for-Life" symposium held in Caen, France (March 2025). This white paper compiles presentations and institutional reports from major CIRT centers in Japan (QST, Gunma), Europe (HIT, CNAO, MedAustron, Cyclhad), and other emerging sites. Data sources include institutional registries, peer-reviewed studies, and national frameworks for hadrontherapy development. Over 50,000 patients have received CIRT worldwide. Japan remains the leading contributor through the J-CROS network, emphasizing registry-based research, hypofractionation, and treatment of mobile tumors. European centers adopt translational programs aimed at further exploiting biological properties of carbon and helium ions. Major barriers include heterogeneous dose-modeling (local effect model vs modified microdosimetric kinetic model), limited image guidance, lack of standardized indications, and funding disparities. Ongoing efforts emphasize registry harmonization, pragmatic trials, and cost-effectiveness modeling. New multi-ion facilities such as Cyclhad (Caen) aim to propose multi-ion therapy, increase the level of evidence for clinical use, and promote research. CIRT is transitioning from pioneering programs to a coordinated global network. Collaborative initiatives and shared data platforms are essential to establish evidence-based indications, optimize biological modeling, and ensure economic sustainability. The "Hadrontherapy-for-Life" initiative calls for a strategic step toward internationally standardized, clinically integrated heavy-ion therapy.
We present a systematic study of identified-hadron production (protons, [Formula: see text], [Formula: see text]) in [Formula: see text] at [Formula: see text] over [Formula: see text]-[Formula: see text] using PYTHIA 8.315 with the Angantyr extension. To match the NA61/SHINE production-event definition, we apply a generator-level selection that suppresses events dominated by a single forward, beam-momentum proton at very small angle (elastic-/quasi-elastic-like "beam survivors"), while leaving hadron production in accepted inelastic events unchanged. We find that both the default Angantyr baseline and the non-suppressed fixed-radii geometry underestimate forward proton yields for [Formula: see text], despite describing mid-angle production reasonably well. To address this, we introduce an impact-parameter-dependent suppression applied after the fixed-radii ND/SD/DD/EL classification, modifying only the relative fractions of ND/SD/DD/EL NN encounter types while leaving particle production within accepted encounters unchanged. Relative to a b-invariant ("flat") alternative, the b-dependent scheme hardens the forward tail without introducing unphysical distortions outside the forward region. Including modest hadronic rescattering enhances soft yields at larger angles while preserving the very-forward high-p component; with a common strangeness retune, the combined setup gives the most consistent description across protons, pions, and kaons. These results show that a minimal impact-parameter-dependent change in the mixture of NN encounter types can improve forward yields at SPS energies without distorting large-angle production, once the NA61/SHINE production-event class is reproduced.
The "Hadrontherapy for Life" symposium in Caen, France, highlighted that a new era of radiobiology is fundamental for advancing particle therapy to the next level. A radiobiology capable of integrating molecular biology and omics technologies is needed to deeply analyze treatment responses and underlying mechanisms. Key challenges discussed at the symposium included tumor hypoxia, which remains only partially mitigated by high-LET radiation, and the specificity of carbon ions, or more broadly, high-LET particles, considered as "new drugs" capable of providing systemic benefits beyond local tumor control, including their potential to promote immunogenicity. Moreover, emerging modalities, such as Ultra High Dose Rate irradiation and spatial fractionated beams, were also discussed, with consensus that all require dedicated and coordinated radiobiological investigations. Infrastructure presentations highlighted the international capabilities of leading centers in Europe and Asia, emphasizing the importance of integrating radiobiology into clinical programs, advancing multi-ion experimentation, and adopting innovative experimental models, such as organoids and/or 3D cell cultures. Participants also stressed the need for greater access to animal experimentation facilities, which are essential for accelerating progress in the field. Furthermore, the meeting underscored translational endpoints such as biomarker development, a hot topic in current radiotherapy. The C400 accelerator enables Caen to incorporate radiobiology from its very inception, establishing a European hub for collaborative research. Round-table discussions emphasized the importance of harmonized protocols, dedicated in vivo irradiation rooms, international training programs with exchange of students and researchers, and comprehensive patient biobanking. In summary, the symposium reinforced the essential role of radiobiology in advancing hadron therapy (HT), providing strategic directions for translational research, infrastructure development, and international collaborations to accelerate personalized and effective particle therapy.
For 50 years, the standard model of particle physics has been very successful in describing subatomic phenomena. In the past quarter of a century, this was challenged by a mismatch between its predictions and precision measurements of the anomalous magnetic moment of the muon, aμ. This disagreement was eventually reconciled, first through a determination in an ab initio lattice calculation1 of the most uncertain theoretical contribution, the leading-order hadronic vacuum polarization (LO-HVP), a μ LO-HVP , and subsequently by experimental results2 and updates of the reference standard-model predictions using lattice results for a μ LO-HVP (ref. 3). Here we present a new calculation for this crucial quantity, obtaining a μ LO-HVP = 715.1 ( 2.5 ) ( 2.3 ) [ 3.4 ] × 1 0 - 10 . This reduces the uncertainty by a factor of 1.6 compared with our earlier computation1. We use a hybrid approach that includes a small, long-distance contribution from experiments in a low-energy regime in which they all agree. Our approach combines the strengths of experimental and lattice data in different energy ranges, achieving better precision than with either alone. Our lattice quantum chromodynamics (QCD) simulations are performed on finer lattices than in ref. 1, allowing for an even more accurate continuum extrapolation. Combined with the calculations of the other standard-model contributions summarized in ref. 3, our result leads to a prediction that differs from the recent measurement of aμ (ref. 4) by only 0.5 standard deviations. This provides a notable validation of the standard model to 11 digits.
Multiturn charge-exchange injection is the primary method of creating high-intensity hadron beams in circular accelerators, and phase space painting during injection enables tailoring of the accumulated phase space distribution. A technique we call eigenpainting allows injection of particles into a single mode of a coupled ring, providing full four-dimensional control of the phase space distribution. Under ideal conditions, uniform eigenpainting generates a linear-force equilibrium distribution in the transverse plane, with zero volume in four-dimensional transverse phase space, even including space charge. We have implemented eigenpainting for the first time in the spallation neutron source accumulator ring. Injecting 8.8  μC of an 800 MeV beam, we obtain a final ratio of intrinsic transverse emittances of ≈2.4. We analyze the effect of space charge on the final distribution through comparison of the reconstructed phase space to particle-in-cell simulations.
We study the helicity correlation of two Λ hyperons produced in unpolarized semi-inclusive deep inelastic scatterings (SIDIS), with one hyperon detected in the current fragmentation region and the other in the target fragmentation region. This observable provides direct access to the spin-dependent fragmentation function G 1 L q and the spin-dependent fracture function l 1 q L even in unpolarized lepton nucleon collisions. Utilizing the perturbative matching of the fracture function, we present numerical predictions for the helicity correlation, revealing significant variations with flavor and kinematic regions. This observable offers a unique way to investigate the spin-dependent hadronization mechanism across both the current and target fragmentation regions. It also provides new insights into the spin transfer effects in SIDIS processes.
Quantum Chromodynamics predicts a phase transition from hadronic matter to quark-gluon plasma (QGP) at high temperatures and energy densities, where quarks and gluons (partons) are no longer confined within hadrons. The QGP forms in ultrarelativistic heavy-ion collisions. Anisotropic flow coefficients, quantifying the azimuthal expansion of produced matter, probe QGP properties. Flow measurements in high-energy heavy-ion collisions show a distinctive grouping of anisotropic flow for baryons and mesons at intermediate transverse momentum - a feature associated with flow imparted at the quark level, confirming QGP existence. The observation of QGP-like features in proton-proton and proton-ion collisions has sparked debate about QGP formation in smaller systems. For the first time, we demonstrate the distinctive grouping of anisotropic flow for baryons and mesons in high-multiplicity proton-lead and proton-proton collisions at the Large Hadron Collider (LHC). These results are described by a model including hydrodynamic flow followed by hadron formation via quark coalescence, consistent with the formation of partonic flowing systems in these collisions.
The outcomes of modern particle physics experiments, such as proton-proton collisions at the Large Hadron Collider at CERN (European Organization for Nuclear Research), depend crucially on the precise description of the scattering processes in terms of the fundamental forces. Among all the known forces that contribute, the limited understanding of the strong nuclear force is a key source of inaccuracy. At the fundamental level, the strong force is described by quantum chromodynamics, the theory of quarks and gluons. Their coupling, αs, becomes weaker at high energies (asymptotic freedom), enabling power series expansions in αs, but the confinement of quarks in hadronic bound states usually requires additional model assumptions. Consequently, determinations of αs from experiment mostly remain with large systematic theory errors1,2. Here we report the model-free determination of αs with unprecedented precision from low-energy experimental input combined with large-scale numerical simulations of the first-principles formulation of quantum chromodynamics on a space-time lattice. The uncertainty, about half that of all other results combined3, originates predominantly from the statistical Monte Carlo evaluation and has a clear probabilistic interpretation. The result for αs describes both low-energy hadronic physics with the help of lattice quantum chromodynamics and high-energy scattering using the perturbative expansion. By removing a source of theoretical uncertainty, our estimate of αs could enable markedly improved analyses of many high-energy experiments4. This will contribute to the likelihood that small effects of yet unknown physics are uncovered, as well as enable stringent precision tests of the Standard Model.
This Letter reports the first calculation of the gravitational form factors (GFFs) of the scalar glueball, performed via lattice field theory in Yang-Mills theory at a single lattice spacing. The glueball GFFs are compared with those of other hadrons as determined in previous lattice calculations, providing strong indications that glueballs have a different gluonic structure than typical hadronic states. A mass radius of 0.263(31) fm is predicted, supporting previous suggestions that the scalar glueball is significantly smaller than other hadrons.
Flat hadron beam collisions, though expected to enhance peak luminosity by about an order of magnitude, have not yet been demonstrated. We identify a new mechanism in which the hourglass-induced synchrobetatron resonance, when coupled to unavoidable machine fluctuations, drives emittance transfer and degrades the flat-beam condition. This effect becomes a critical limitation for maintaining flat hadron beams on operational timescales. Using Electron-Ion Collider (EIC) parameters, we show that it imposes stringent tolerances on machine noise. We further demonstrate that a dynamic focusing scheme, implemented by combining sextupoles with crab cavities, suppresses the resonance and relaxes these constraints, providing a practical path toward robust flat-beam collisions in future lepton-hadron colliders.
In phenomenological studies of multiparticle production, transverse-momentum spectra measured in experiments frequently display an approximately power-law falloff, for which the Tsallis-type functional form is commonly employed as an effective parametrization. Within this framework, the emergence of such spectra is interpreted as a manifestation of nonextensive statistical behavior. An analogous power-law structure, however, can be reproduced without explicitly postulating Tsallis statistics by assuming the presence of intrinsic fluctuations of the local temperature (T) in the hadronizing medium; in that case, the observed deviations from a purely exponential spectrum are encapsulated by the nonextensivity index (q). We show that temperature fluctuation mechanisms capable of generating Tsallis-like power-law distributions in multiparticle production necessarily induce nontrivial inter-particle correlations among the emitted hadrons. Building on this observation, we outline a strategy to discriminate fluctuations realized on an event-by-event basis from those arising predominantly through event-to-event variability. Such a separation may be particularly pertinent for the characterization of high-multiplicity (high-density) final states produced at the Large Hadron Collider.
The mass of the lightest hadron in nature, the pion, is one seventh of that of the nucleon and one tenth of the mass of its first excited state, the π(1300). This enormous energy difference opens an interesting window into the confinement of quarks and the structure of the lightest hadrons. In this Letter, we provide the first calculation of resonance parameters of the π(1300) from lattice quantum chromodynamics (QCD). For this purpose, recently derived state-of-the-art tools are adapted and applied both in the construction of three-hadron operators and for mapping finite-volume spectra to infinite-volume amplitudes, subsequently analytically continuing these to complex energies. For our heavy pion mass ensembles, we find a clear signal of the resonance. Making a simple assumption of vanishing pion mass dependence for the three-body force, but incorporating constraints from Chiral Perturbation Theory for all the two-body channels, enables a robust extrapolation to the physical point. Applying model averaging, we extract a pole position of M_{π(1300)}=(1169±46)-i(62_{-62}^{+168})  MeV supporting values from phenomenology.
We study the QCD scaling behavior of the small-angle energy-energy correlator (EEC), focusing on the transition between its perturbative preconfinement and nonperturbative postconfinement regimes. Applying the light-ray operator product expansion (OPE), we develop a formalism that describes the scaling of the EEC with the input energy Q in the transition and the postconfinement region, where the latter quantum scaling is determined by the J=5 DGLAP anomalous dimension. A key result of our Letter is a novel connection between the light-ray OPE and the dihadron fragmentation function (DFF), where we show that the nonperturbative OPE coefficients correspond to moments of the DFF. This finding establishes a new paradigm for studying hadronization. Our theoretical predictions are validated against Monte Carlo simulations for both e^{+}e^{-} and pp collisions, showing excellent agreement. The potential role of the quantum scaling in the precision determination of α_{s} is also discussed.
Supernova remnants (SNRs) have been considered as the primary contributors to cosmic rays (CRs) in our Galaxy. However, the maximum energy of particles that can be accelerated by shocks of SNRs is uncertain, and SNRs' contribution to CRs around PeV energies is unclear. In this Letter, we present observations of high-energy γ-ray emission from the SNR IC 443 using the Large High Altitude Air Shower Observatory (LHAASO). The morphological analysis reveals a pointlike source whose location and spectrum are consistent with those of the Fermi-LAT-detected compact source with π^{0}-decay signature, and a more extended source that is associated with a newly discovered Fermi source. The spectrum of the point source can be described by a power-law function with an index of ∼3.0, extending beyond ∼30  TeV without apparent cutoff. Assuming a hadronic origin of the γ-ray emission, the 95% lower limit of accelerated protons reaches about 300 TeV. The extended source might be associated with IC 443, SNR G189.6+3.3, or the putative pulsar wind nebula CXOU J061705.3+222127, and can be explained by either a hadronic or a leptonic model with particles reaching hundreds of TeV. These LHAASO results provide compelling evidence that sub-PeV CRs can be accelerated efficiently by shocks of SNRs.
A reliable calculation of radiative corrections to τ→ππν_{τ} decays is an important prerequisite for using hadronic τ decays for a data-driven evaluation of the hadronic-vacuum-polarization (HVP) contribution to the anomalous magnetic moment of the muon, a_{μ}^{HVP,LO}[ππ,τ]. In this Letter, we present an improved model-independent analysis of these radiative corrections, including, for the first time, effects beyond pointlike pions in the evaluation of the loop diagrams. These structure-dependent corrections, implemented via a dispersive representation of the pion form factor, lead to significant changes compared to previous calculations due to enhancements near the ρ(770) resonance. We also devise strategies for the matching to chiral perturbation theory and a stable implementation of the real corrections down to the two-pion threshold, which shows that some higher-order isospin-breaking corrections need to be kept due to a strong threshold enhancement. Finally, we perform dispersive fits to the currently available τ→ππν_{τ} spectra and discuss the consequences for isospin-breaking corrections in the evaluation of a_{μ}^{HVP,LO}[ππ,τ].
The vacuum is now understood to have a rich and complex structure, characterized by fluctuating energy fields1 and a condensate of virtual quark-antiquark pairs. The spontaneous breaking of the approximate chiral symmetry2, signalled by the nonvanishing quark condensate ⟨ q q ¯ ⟩ , is dynamically generated through topologically nontrivial gauge configurations such as instantons3. The precise mechanism linking the chiral symmetry breaking to the mass generation associated with quark confinement4 remains a profound open question in quantum chromodynamics (QCD)-the fundamental theory of strong interaction. High-energy proton-proton collisions could liberate virtual quark-antiquark pairs from the vacuum that subsequently undergo confinement to form hadrons, whose properties could serve as probes into QCD confinement and the quark condensate. Here we report evidence of spin correlations in Λ Λ ¯ hyperon pairs inherited from spin-correlated strange quark-antiquark virtual pairs. Measurements by the STAR experiment at the Relativistic Heavy Ion Collider (RHIC) at Brookhaven National Laboratory reveal a relative polarization signal of (18 ± 4)% that links the virtual spin-correlated quark pairs from the QCD vacuum to their final-state hadron counterparts. Crucially, this correlation vanishes when the hyperon pairs are widely separated in angle, consistent with the decoherence of the quantum system. Our findings provide a new experimental model for exploring the dynamics and interplay of quark confinement and entanglement.
Regge trajectories provide a simple geometric picture of hadron spectra, but particles at higher-spin and missing states often scatter from linear fits, raising questions about both experimental completeness and theoretical limits of the model. Here, we develop a unified data framework that integrates standard particle listings with hypergraph-based decay features, enabling systematic comparison across baryons and mesons. We employed orthogonal distance regression with bootstrap resampling to quantify uncertainties in slope and intercept estimates, while hypergraph-derived structural invariants (community purity, motif z-scores, and product entropy) serve as quantitative predictors of spectroscopic regularity, establishing decay topology as a microscopic determinant of macroscopic Regge behavior. Applying this hybrid approach to 20 [Formula: see text] baryon resonances, we obtain strong linear correlation (R2 = 0.90) with slope [Formula: see text] GeV-2, though elevated scatter ([Formula: see text]) correlates strongly with resonance width (r = 0.88, p < 0.001). While hypergraph features did not significantly improve explanatory power beyond quality controls in this pilot dataset ([Formula: see text], p = 0.42), we introduce a hypergraph-informed confidence framework for missing resonance predictions, where structural decay coherence provides quantitative reliability metrics beyond traditional trajectory extrapolation. Together, these results demonstrate how combining trajectory analysis with network-inspired methods can improve hadron classification and provide concrete predictions for future experimental searches.
In this Research Perspective, we briefly review the diffusion wake, a distinctive consequence of the Mach-cone wake induced by the supersonic jets in ultrarelativistic heavy-ion collisions. The diffusion wake depletes soft hadrons in the direction opposite to the propagating jet. According to coupled transport and hydrodynamic simulations, a valley in the 2-dimensional jet-hadron correlation in azimuthal angle and rapidity arises on the top of the multiple-parton interaction ridge as an unambiguous signal of the diffusion wake induced by γ -jets in heavy-ion collisions. In dijet events with a finite rapidity gap, the rapidity asymmetry of the jet-hadron correlation has been shown to be a robust signal of the diffusion wake. The same rapidity asymmetry can also be applied to γ -jet events, and both are background-free. Experimental measurements of these signals can provide valuable insights into the properties of the quark-gluon plasma formed in high-energy heavy-ion collisions.