2026-09-30 | | Total: 59
We investigate self-interacting dark matter in an $O(N)$-symmetric real scalar multiplet model with Higgs-portal interactions involving the Standard Model Higgs and a real singlet scalar. After spontaneous breaking of the singlet $\mathbb{Z}_2$ symmetry, the singlet scalar mixes with the SM Higgs boson. We focus on the regime where the resulting light scalar serves as a mediator, leading to sizable and velocity-dependent self-interactions. We impose constraints from Higgs measurements, rare $B$-meson decays, the dark matter relic abundance, spin-independent direct-detection searches, and scalar coupling perturbativity. We find viable parameter regions in which destructive interference between scalar-mediated scattering amplitudes suppresses the direct-detection cross section while allowing self-interactions relevant to dwarf-galaxy scales. Among these viable parameter regions, a subset also satisfies electroweak vacuum stability up to the Planck scale. Finally, we show that future long-lived-particle searches can probe part of the viable parameter space, providing complementary tests of Higgs-portal self-interacting dark matter with a light scalar mediator.
The LHCb Collaboration has observed a narrow state in the $\bar{B}_s^0π^0$ spectrum, about 74 MeV below the $B^-K^+$ threshold, which is the natural candidate for the bottom partner of the $D_{s0}^*(2317)$. Such a state had been predicted as a $\bar{B}K$ bound state by unitarized heavy-meson chiral perturbation theory with lattice-QCD input. We analyze the LHCb spectrum with these coupled-channel amplitudes, including isospin breaking through the physical meson masses, $π^0$-$η$ mixing, and $m_u\neq m_d$ in the chiral amplitudes, thereby allowing the $\bar{B}_{s0}^*$ to decay through its only open strong channel, $\bar{B}_s^0π^0$. With the parameters fixed by lattice QCD and no experimental input from the LHCb measurement, the predicted state lies within 1.1-1.7 standard deviations of the observed peak. Adjusting a single parameter, the spectrum is described as well as with the zero-width limit of the LHCb analysis. The new state thus emerges naturally from these amplitudes, with a large $\bar{B}K$ component and an isospin-violating width of tens of keV, and its $J^P=1^+$ heavy-quark spin partner is predicted with a precision comparable to that of the measured mass.
We consider a planar and a non-planar two-loop Feynman integral family arising in non-factorisable quantum chromodynamics corrections to Higgs-boson production via vector boson fusion. We construct basis integrals which fulfil canonical differential equations with respect to seven variables for the external kinematics and an internal mass. One of the integral families involves nested square roots. We present strategies for deriving canonical master integrals, $\mathrm{d}\log$ forms and analytic continuation for such cases. Including also families with crossed massless legs, we express the $\varepsilon$-expanded integrals in terms of an algebraically independent function basis. For their evaluation, we employ numerical integration of the differential equations.
We consider one-loop QCD corrections to Higgs-boson production via vector-boson fusion at hadron colliders. As an input for future exact NNLO computations, we derive the one-loop helicity amplitudes to higher orders in the dimensional regulator. We present a direct calculation of the amplitudes with spinor-helicity techniques and different $γ_5$ schemes. In addition, we compare this approach to a projector-based calculation. In order to derive analytical results to higher orders in $ε$, we calculate the pentagon master integrals by the method of $ε$-factorised differential equations and express them in terms of an algebraically independent set of basis functions.
By performing a fit to the merged BES 2008 and BESIII 2024 $e^+e^-\to D\bar D$ cross-section data, we find that the $ψ(2\mathrm S)$ subthreshold pole influences this process markedly. The statistical significance of the $ψ(2\mathrm S)$ as a subthreshold pole is 8$σ$. A fit assuming it and seven resonances excels with the fit quality characterized by $χ^2/\mathrm{DOF}=240.4/284$, which means a $p$-value of 97\%. The interference of resonances with the continuum provided by the subthreshold pole changes their shapes, positions, and widths.
We compute jet quenching observables--the nuclear modification factor $R_{AA}$ and elliptic anisotropy coefficient $v_2$--for open-heavy flavor mesons with high transverse momentum in an improved CUJET model, in which the strong coupling parameter as well as the temperature-dependent fraction of chromo-magneto monopoles are extracted through Bayesian inference using high-$p_T$ light-hadron $R_{\mathrm{AA}}$ and $v_2$ measurements from RHIC and the LHC. We show in this work that such a Bayesian-inference-improved CUJET model provides a good description of available high-$p_T$ charm- and bottom-hadron observables based on the same posterior parameters. The corresponding heavy-quark spatial diffusion parameter, $2πT D_s$, is also computed, which is broadly consistent with current phenomenological and lattice-QCD constraints. These results support a universal description for the medium modification of high-$p_T$ light- and heavy-flavor energy loss and indicate an important role played by the chromo-magnetic degrees of freedom near the confinement transition region of quantum chromodynamics (QCD).
We begin a systematic study of neutrinoless muon-to-positron conversion, $μ^-+(A,Z)\to e^++(A,Z-2)$, which violates lepton number by two units and changes charged-lepton flavor. This first paper in a series constructs nuclear operators for local dimension-nine interactions containing a scalar or pseudoscalar lepton bilinear and two contracted vector or axial quark currents. The covariant nucleon currents are built with all six vector and axial form factors. For the working nuclear operator we neglect second-class currents and adopt a pion-pole-improved reduction: non-pole recoil terms generated by the displayed one-body currents are retained through second order in inverse nucleon mass, whereas the pion-pole longitudinal axial current is kept intact before the two currents are contracted. We preserve the two independent momentum transfers at the nucleon vertices and the Coulomb-shifted spatial phase of the outgoing positron. Starting from the spacetime transition element, we show explicitly how the interaction-coordinate integral fixes the sum, but not the sharing, of the two current momenta. Unequal form-factor arguments then generate relative-coordinate multipoles absent in a common-momentum approximation. We give an exact recoupling into two-body nuclear tensors and specialize it to the physically important $0^+\to0^+$ transition. In this case the external nuclear rank is zero and all magnetic sums collapse, although finite-momentum internal orbital and spin tensors remain. The result is expressed in a form directly contractible with charge-changing two-body transition densities. General six-form-factor currents, arbitrary nuclear angular momenta, and power-counting qualifications are collected in appendices.
The relativized quark model of Godfrey and Isgur (1985) remains a reference for meson spectroscopy four decades after its publication, yet despite its impact the calculation behind it has never been publicly available. I present an independent reproduction of the whole paper, released as the open-source Julia package GIModel.jl with companion packages for transition observables and for the comparison with the original. The published parameters are used without refitting, and the paper's three-stage harmonic-oscillator algorithm is cross-checked with an independent finite-difference solver. All 209 levels of the paper's spectrum figures are reproduced within the precision of the printed diagrams. The strong-decay, radiative, annihilation, and electromagnetic tables are reproduced with median model-to-paper ratios between 0.97 and 1.06, and the $1P$ mixing angles agree with the values the authors revised after 1985. With the calculation open, the model's own wavefunctions become available for decay amplitudes and structure visualization. The package, its teaching material, and its agent-accessible interface make the model a baseline for experimental questions and theoretical extensions.
We study the interaction of $φ$ with $h_1(1380)$, considering the latter as a $K^*\bar K$ molecular state, and find the generation of a state with the same quantum numbers, together with a mass and width compatible with the experimental results, as those of $X(2370)$, recently advocated as a candidate for a glueball. Within our description, the suppressed decay modes of $X(2370)$ to $γω$, $ωφ$ and $K^*\bar K$, reported in the experimental paper as clues to its glueball nature, are naturally explained. Conversely, the non $K^* \bar K$ decay to $K \bar K π$ (the latter being the invariant mass where the resonance is observed) is shown to be largely enhanced due to the combination of a triangle singularity and the presence of $K(1630)$ in the $ρφ\to \bar K$ vertex involved in the triangle loop.
We determine the decay rate $Γ$ of a scalar particle bound to a light nucleus, numerically and analytically. As a model system, we consider a negatively charged kaon in the ground state of a hydrogen-like ion decaying into charged and neutral scalar pions, $K^-\to π^-π^0$. All particles are assumed point-like; the nucleus is infinitely heavy. The partial wave expansion converges slowly, as in a recent analysis of the bound-muon decay. To overcome this difficulty, we introduce a subtraction method: the free $π^-$ contribution is evaluated separately, while the partial-wave sum is used only for the Coulomb distortion. The ratio of bound-to-free decay rates, $Γ/Γ_0$, is found for light nuclei, $1\leq Z\leq 10$. The numerical calculation is complemented by a diagrammatic small-$αZ$ expansion. We find $ Γ/Γ_0 = 1-(αZ)^2/2 + 17(αZ)^4/24 +\mathcal{O}((αZ)^5)$, in agreement with the all-orders numerical calculation at low $Z$.
The scale evolution of parton distributions is governed by splitting functions. We compute the four-loop pure-singlet contribution to the quark-to-quark splitting function in perturbative QCD. Together with the known non-singlet contribution, our result completes the quark-to-quark splitting function at this order. We confirm previous partial results and obtain, for the first time, a fully analytic expression valid for all momentum fractions. We examine its behavior at small momentum fractions and provide precise numerical representations suitable for parton evolution.
The LUX-ZEPLIN (LZ) experiment has reported a single nuclear-recoil candidate at $248\pm23\,(\mathrm{stat})\pm23\,(\mathrm{sys})$ keV. We suggest that it comes from a GeV-scale dark state $χ_1$ that downscatters exothermically, $χ_1A\toχ_2A$, and is continuously replenished in the halo by the slow decay $χ_0\toχ_1χ_1$ of a long-lived parent. Because the recoil energy is supplied by the mass splitting rather than by the halo velocity, the event is a narrow line with no low-energy tail. Kinematics requires $m_{χ_1}\gtrsim247$ MeV and fixes the splitting for each $m_{χ_1}$. The same line should appear at 700$-$810 keV in argon for $m_{χ_1}=0.3-10$ GeV. Normalizing to one LZ event also fixes the local rate of the radiative decay $χ_1\toχ_2γγ$, up to the mediator's photon coupling. A rough $γ$-ray estimate then requires this coupling to be $\lesssim10^{-2}$ of a proton-loop reference at $m_{χ_1}=10$ GeV, but $\lesssim10^{-9}$ at 0.3 GeV, which favours daughters near 10 GeV.
Renormalisation-group evolution enters PDF determinations through the running coupling and the splitting functions. At fixed perturbative order, the solution of the evolution equations is not unique beyond the nominal accuracy, and this residual freedom can be probed through resummation-scale variations. This contribution summarises an NNLO study of this effect in inclusive HERA DIS fits performed with the xFitter fitting framework. The induced shifts are sizeable in a low-scale baseline fit, especially for the gluon and sea PDFs, and can affect LHC and FCC top-quark pair predictions at the few-percent level. Raising the PDF starting scale and applying a higher DIS virtuality cut considerably reduces the effect, a conclusion also supported by treating the resummation scale as a free parameter constrained directly by the data.
The LUX-ZEPLIN experiment has reported a nuclear-recoil candidate near 248 keV. We investigate a decay source for an endothermic interpretation: a long-lived dark parent produces energetic ground-state particles that upscatter on xenon. The parent mass fixes the injection energy, while its lifetime determines the absolute flux through the Galactic dark-matter column density. We include the Doppler broadening from parent and observer motion and use exact scattering kinematics with natural xenon isotopes. For a scalar contact interaction, a benchmark with a 10 MeV incident state and a 506 MeV parent gives an efficiency-weighted recoil spectrum peaked at 248.5 keV. One event in a specified high-energy window requires $τ_X \simeq 4.4 \times 10^{20}$ s $f_X [σ_{\text{ref}}/(10^{-36} \text{cm}^2)]$, with a normalization sensitive to the nuclear response. We quantify the accompanying elastic-scattering requirement and the target-dependent thresholds. The calculation identifies a decay realization of the boosted flux and the portal conditions needed to preserve the recoil signal; a detector-level fit and ultraviolet-specific constraints remain necessary for a complete experimental assessment.
Many experiments have provided possible hints of gluon saturation, including geometric scaling in deep inelastic scattering (DIS) at HERA, the suppression of forward particle yields in d+Au collisions relative to p+p collisions at RHIC, and the suppression of away-side peaks in di-hadron correlations at forward rapidities in proton-nucleus collisions at RHIC and the LHC. Previous studies within the Color Glass Condensate (CGC) framework have primarily focused on fits to HERA deep inelastic scattering data, using the resulting dipole amplitudes to predict hadronic observables in proton-nucleus and proton-proton collisions. A simultaneous fit of DIS and forward particle production data, and a systematic study of their interplay, has so far been absent. In the work presented here, we focus on two processes: inclusive DIS and forward single inclusive hadron production in p+p collisions, incorporating data from HERA, RHIC, and the LHC. We solve the Balitsky-Kovchegov (BK) small-x evolution equation and perform a global fit to the combined datasets. Our analysis accounts for uncertainties in fragmentation functions and provides final uncertainties for the fitted parameters. The outcomes demonstrate the potential of the CGC formalism to provide a universal description of these diverse processes at the level of the dipole amplitude. In these proceedings, after very briefly summarizing the setup and our results, we highlight one particular aspect of the analysis, namely the complete uncertainty treatment. We validated the final uncertainties with two independent methods: the Hessian method and the Bayesian inference method, finding very good agreement.
We compute the combinations $E_M=H^{u+d}+E^{u+d}$ and $H_E=H^{u-d}+tE^{u-d}/(4M_N^2)$ of unpolarized generalized parton distributions (GPDs) of the nucleon in the chiral quark soliton model at nonzero skewness $ξ$ and nonzero momentum transfer $t$. The combinations contain the GPDs $H^{u-d}$ and $E^{u+d}$, which are subleading in the expansion in the inverse number of colors. Both combinations are double sums over the quark states. A multipole expansion of the light-cone operators makes the sums computable at every $ξ$ and every $t\le0$. The contribution of the Dirac continuum is summed over the quark states, without an interpolation formula and without a gradient expansion. The first moments satisfy the form-factor sum rules, and the second moments satisfy polynomiality.
We present a complete analytic calculation of the next-to-leading-order (NLO) QCD cross sections for the production of an electroweak vector boson ($W^\pm$, $Z/γ^*$) at finite transverse momentum, decomposed into the full set of helicity cross sections that determine the angular coefficients $A_0$-$A_7$ of the lepton-pair angular distribution in the Collins-Soper frame. The calculation includes the unpolarised cross section and all eight helicity projections for both neutral- and charged-current Drell-Yan production, the axial quark-triangle contributions in the $m_t\to\infty$ limit and with full top-quark-mass dependence, and closed-form expressions for all coefficient functions. The results are implemented in the DYTurbo program and validated at the level of a few times $10^{-4}$ against its numerical implementation derived from MCFM, as well as through the analytic $q_T\to0$ limit. The analytic implementation achieves a speed-up of more than five orders of magnitude with respect to the numerical code, enabling fast and numerically stable predictions for precision electroweak measurements at the LHC. The calculation was carried out with extensive use of large language models (Anthropic Claude) as autonomous agents; we document the methodology and the associated computational and token budgets.
We investigate the structure of the pion in a magnetized nuclear medium at finite temperature by analyzing the in-medium modifications of its parton distribution functions (PDFs) and electromagnetic form factors (EMFFs). We use a unified framework combining the chiral SU(3) quark mean field (CQMF) model and the light-cone quark model (LCQM) where the in-medium quark masses are evaluated using the CQMF model and subsequently employed as inputs in the light-cone formalism. The study systematically explores the effects of magnetic field, baryon density, temperature, and isospin asymmetry on the pion structure. In addition, Dirac sea contributions are incorporated to account for vacuum polarization effects in the medium. Our results of PDFs and EMFFs provide valuable insights into the modification of pion properties under extreme conditions, highlighting the significant role of the surrounding medium in altering its internal structure.
We present nCTEQ26OHF, a global analysis of nuclear parton distribution functions (nPDFs) within the nCTEQ framework. In this analysis, constrasting the main nCTEQ26 release, we replace the Crystal-Ball function, an effective matrix element that is used in nCTEQ26 for open heavy-flavor (OHF) and quarkonium production, by an SACOT GM-VFNS calculation for the OHF predictions. The resulting nCTEQ26OHF PDFs are compatible with the nCTEQ26 PDFs within uncertainties, implying that the OHF data alone provide ample constraint.
We infer the initial condition for the Balitsky-Kovchegov (BK) evolution equation at next-to-leading order accuracy from HERA proton structure function data. A simultaneous description of both the total cross section and charm production measurements is obtained. We then show that inclusive hadron production in deep inelastic scattering can provide complementary constraints on the non-perturbative initial condition of the BK equation. As a more differential observable, it serves as a sensitive probe of gluon saturation phenomena. We demonstrate this by calculating predictions for the nuclear modification factor in hadron production in photon-nucleus scattering at leading order in the dipole picture. We find that saturation effects lead to significant nuclear suppression at small $x_\mathrm{Bj}$.
We introduce the basics of Effective Field Theories, using the Standard Model Effective Field Theory as our central example. The first part of this chapter gives a pedagogical introduction to EFTs - including matching, resummation and the construction of an operator basis - with explicit examples and provides a starting point for newcomers. The second part reviews the current status of SMEFT and collects suitable references for those who want to go deeper.
The electron-positron Future Circular Collider will provide an unprecedented opportunity for precision studies at the $Z$ pole. In addition to improving the electroweak measurements pioneered at LEP, Tera-$Z$ statistics and modern jet flavor reconstruction could make rare, flavor-resolved final states accessible as precision observables. We study the exclusive four-body processes $e^+e^-\to Z\to Q\bar Q\ell^+\ell^-$ with $ Q=c,\,b$ and $\ell=μ,\,τ$ as probes of flavor-specific new physics contact interactions. A fast analysis yields representative multi-TeV reaches for vector-current interactions involving bottom or charm quarks and muons or taus, providing a baseline for future detector-level studies. For charm-current interactions, the projected direct sensitivities are comparable to those obtained indirectly from electroweak precision observables. Two-dimensional fits further show that these processes retain quark-flavor information that can be poorly resolved by the renormalization-group-induced electroweak response. These results illustrate how flavor-resolved four-body $Z$ processes can bridge the electroweak and flavor programs at FCC-ee, opening a new class of complementary precision observables.
The $J^{P}=2^{+}$ $cc\bar c\bar s$ and $bb\bar b\bar s$ tetraquarks are studied in parallel within the two-point QCD Sum Rules framework using the same symmetric axial-vector$\otimes$axial-vector interpolating current. The current-specific OPE sectors retained through dimension seven within the stated background-field truncation are evaluated with a common symbolic construction, with the bottom channel obtained by the heavy-flavour replacement $c\to b$. The same pole-dominance and OPE-convergence criteria are imposed in both channels. The adopted charm region is $\sqrt{s_0}=5.7\GeV$ with $3.6\le M_B^2\le4.3\GeVsq$, whereas the bottom region is $\sqrt{s_0}=14.3\GeV$ with $9.8\le M_B^2\le10.8\GeVsq$. At representative Borel points, the extracted masses are $5.239\GeV$ and $13.832\GeV$, respectively. After the QCD-input, continuum-threshold, and Borel-window variations are combined, the final results are $m_{cc\bar c\bar s,\,2^+}=5.239\pm0.061\GeV$ and $m_{bb\bar b\bar s,\,2^+}=13.832\pm0.059\GeV$.
New physics beyond the Standard Model (SM) may induce new four-fermion interactions. We investigate the sensitivity of radiative muon decay, $μ^-\to e^-\barν_eν_μγ$, to such four-fermion interactions of mass dimension seven and eight. For each operator, we calculate its interference with the SM amplitude and obtain an analytic expression for the decay rate. We apply these to a log-likelihood comparison for a benchmark sample of $10^5$ events, and find that these operator contributions can be distinguished from the SM prediction at characteristic coefficient magnitudes between $10^{-9}$ and $10^{-8}$.
The LUX--ZEPLIN (LZ) Collaboration has reported a nuclear-recoil candidate at $248\pm23~({\rm stat})\pm23~({\rm sys})~{\rm keV}$ in an extended high-energy search window. We investigate whether this event can arise from the hyperfine excitation of hydrogen-like atomic dark matter. For a TeV-scale dark atom, a hyperfine splitting near $340~{\rm keV}$ places xenon scattering close to the observed value by naturally selecting the high-velocity tail of the Galactic halo while suppressing the leading low-energy elastic response. The resulting endothermic kinematics predict a pronounced target hierarchy: the benchmark transition is inaccessible on Ar and Ge, lies close to threshold on Xe, and remains open on heavier targets such as W, which is specific to our model parameters. Because xenon probes the extreme high-speed tail, the signal exhibits a large annual modulation and a strong dependence on the assumed halo distribution. We also examine how the infall of the Large Magellanic Cloud enhances the high-speed tail of the Milky Way's local dark matter distribution, extending the kinematic reach to larger hyperfine splittings.