2026-05-29 | | Total: 26
In this paper, we present HyperPrecision, a Mathematica package for high-precision numerical evaluation of general Horn-type multivariate hypergeometric functions and their Laurent expansions in a small parameter $ε$. Such functions appear widely in physics and mathematics, with applications ranging from quantum field theory and string theory to number theory and statistics. Their high-precision numerical evaluation, however, remains challenging, since their defining series converge only in restricted domains and analytic continuation beyond these domains is, in general, non-trivial. HyperPrecision addresses this problem by automatically constructing the Pfaffian system of partial differential equations for a given hypergeometric function and restricting it to a one-dimensional contour in the space of variables connecting the starting to the target point. The resulting ordinary differential equation is then solved by the Frobenius method, with boundary conditions fixed analytically by the defining series. We illustrate the use of the package by evaluating commonly occurring multivariate hypergeometric functions, including the Appell $F_1$, $F_2$, $F_3$, and $F_4$ functions, the Horn $G$- and $H$-series, and the Lauricella $F_A$, $F_B$, $F_C$, and $F_D$ functions, as well as by considering applications to angular integrals, Feynman integrals, and cosmological and holographic correlators.
We study the thermodynamics of the symmetry-improved Cornwall-Jackiw-Tomboulis (SICJT) formalism and apply it to a low-energy effective theory of QCD. In the symmetry-improved formulation, Ward-Takahashi identities are restored by auxiliary sources whose values are fixed self-consistently by the equilibrium state. While this construction improves the symmetry properties of the loop-wise truncated two-particle-irreducible (2PI) theory, it also makes the thermodynamic interpretation of the pressure nontrivial. We formulate several pressure prescriptions, including the conventional vacuum-subtracted pressure, a source-matched subtraction, and a pulled-back pressure in which the explicit source-induced energy shift is removed. Using the three-flavor linear sigma model with quarks, we analyze the equation of state, isentropic trajectories, adiabatic sound velocity, and trace anomaly across the chiral transition. We find that the global thermodynamic structure is stable under the different prescriptions, while quantitative differences are concentrated near the crossover and first-order transition region. These results establish a practical framework for constructing thermodynamically consistent observables in symmetry-improved 2PI approaches.
We present a quantum-field-theoretic treatment of massive chiral fields in which particles possess well-defined chirality and helicity. This framework reproduces the chiral oscillation formula previously obtained in first-quantized approaches and provides a consistent description of weak-interaction processes. We further derive corresponding chiral-energy uncertainty relations.
We study the one-loop Polyakov-loop effective potential in pure gluonic matter under constant acceleration. We perform the computation in both the Euclidean Rindler spacetime and the optical spacetime, which are related via a conformal transformation. The results from the two formulations correspond to physically different observables, and we clarify their connection to specific components of the energy-momentum tensor. This identification resolves a discrepancy previously noted for fields on conical backgrounds. For the Polyakov-loop expectation value, we should minimize the effective potential computed in the optical metric formulation, which concludes that real acceleration strengthens deconfining properties. We also discuss analytic continuation from real to imaginary acceleration and find a perturbatively confined phase. We point out some suggestive similarities and differences between systems under imaginary acceleration and imaginary rotation.
In this study we apply the Remler formalism to charmonium production at SPS and GSI/FAIR energies in order to investigate the effects of baryon-rich matter on charmonium production and dissociation in heavy-ion collisions within the Parton-Hadron-String Dynamics (PHSD). As a first step the Remler formalism is tested in p+p collisions and then applied to p+A collisions in order to extract the nuclear absorption cross section of charmonium, which is then utilized in heavy-ion collisions. We find that the Remler formalism successfully describes charmonium production in heavy-ion collisions at SPS energies when an in-medium heavy quark potential is implemented, in which $J/ψ$ dissociates near $T_c$. Finally the same formalism is applied to the GSI/FAIR energies, where we estimate charmonium production.
Axion dark matter in high-scale inflation is subject to the isocurvature constraint, since quantum fluctuations of the axion field during inflation may exceed the current CMB bound. One conventional way to suppress these fluctuations is to assume that the Peccei-Quinn field has a large expectation value during inflation. However, this mechanism becomes ineffective when the axion domain wall number is large. In this work, we point out that a high-quality axion protected by a discrete gauge symmetry can naturally evade this problem. A Peccei-Quinn-violating but gauge-invariant operator induces a large effective axion mass during inflation, thereby suppressing the axion fluctuation. The same setup can address both the axion quality problem and the isocurvature problem, while leading to a prediction for the axion parameter space to be verified in future experiments.
Contrary to conventional expectations, we identify a class of extended scalar-sector scenarios in which final states with two, three, or four Higgs bosons constitute the leading discovery channels for new physics at the LHC. In these scenarios, higher-dimensional interactions, together with suppressed Higgs-scalar mixing near the alignment limit, reorganize the decay patterns of new scalar states, suppressing conventional modes while enhancing multi-Higgs final states. We illustrate the emergence of dominant triple- and quadruple-Higgs signatures in two representative realizations: a single-scalar extension of the Standard Model, where higher-dimensional operators suppress conventional two-body decays while preserving couplings to higher-multiplicity Higgs final states; and a two-singlet scenario, where similar signatures arise through cascade decays with a simpler operator structure. In both cases, the new scalar states can be produced via gluon fusion, yielding potentially observable rates for multi-Higgs production at the LHC. Although both realizations lead to identical final states, they exhibit distinct kinematic features reflecting their underlying topologies, providing a direct handle on the dynamics.
The CMS collaboration has reported a measurement of $Υ$ decays to ditaus using $61.9~{\rm fb}^{-1}$ of scouting data. If interpreted as the decay of $Υ(1S,2S,3S)$, the measured ditau rate is more than ten times that seen in the dimuon final states at the $\sim 3 σ$ level, and is likewise inconsistent with the branching ratios measured at $B$-factories. If confirmed with more data and at higher significance, such a violation of lepton flavor universality would necessitate new physics. In this Letter, we present a simple model with a light pseudoscalar coincidentally near the $Υ(1S)$ mass, which mixes with two Higgs doublets in the alignment limit. Such a particle naturally decays primarily to taus and evades all existing experimental constraints, while implying a number of predictions that can be tested in the near future.
In this work we generalize the thermal decoupling, or freeze-out, of weakly interacting massive particle dark matter within the Tsallis nonextensive formalism. The generalization is implemented through $q$-deformed distribution functions obtained from the maximum entropy principle with Curado-Tsallis constraints. The Tsallis parameter $q$, which measures deviations from extensivity with respect to the limit $q=1$, is promoted to a dynamical quantity depending on the dimensionless variable $x=m_χ/T$, where $m_χ$ is the dark matter mass. This dynamical evolution is characterized by a relaxation toward extensivity, while requiring that the nonextensive deformation is not completely erased before freeze-out. We solve the Boltzmann equation assuming a sectorial deformation, where only the dark matter equilibrium abundance is generalized and the radiation background remains extensive. The relic abundance is computed for different dark matter masses and initial values of the Tsallis parameter. From this evolution, we extract the residual value $q_χ^{\rm fo}$ at freeze-out, which is then used as the initial input for a phenomenological memory channel. This channel propagates the residual nonextensivity into the radiation sector, specifically into the electron-positron plasma and neutrinos, while photons are kept extensive in order to avoid direct tensions with CMB physics. The resulting deformation modifies the neutrino energy density and the photon reheating contribution, producing a correction to $N_{\rm eff}$. We compare the predicted values with the compressed CMB+BAO constraint on $N_{\rm eff}$ and find that the residual-memory scenario can remain phenomenologically compatible with current bounds.
Astrophysical neutrinos provide a unique probe of neutrino flavor changes over cosmological baselines. While the tau component of the neutrino flux is expected to arise almost entirely from mixing, current measurements rely primarily on rare double-cascade signatures. We investigate a complementary method to measure the tau fraction using the visible inelasticity of starting track events in neutrino telescopes. Muonic decays of tau leptons produce tracks with systematically larger visible inelasticity than those from muon neutrino interactions, potentially enabling statistical separation of the two flavors. Using realistic IceCube exposures and detector performance, we show that this observable already yields competitive sensitivity to the tau-to-muon flux ratio, $R_{τμ}$, achievable with existing data. This approach may further enable flavor measurements of individual sources and the selection of tau-enhanced source catalogs. Starting-track inelasticity thus provides a powerful and immediately accessible probe of astrophysical neutrino flavor and of potential physics beyond standard neutrino mixing.
We present a comprehensive study of dark matter phenomenology in standard model extensions featuring an electroweak triplet scalar or fermion with hypercharge $Y = 0$ or $Y = 2$. These minimal triplet extensions provide well-motivated dark matter candidates stabilised by the $Z_2$ discrete symmetry. We perform a detailed analysis of the parameter space consistent with current cosmological and experimental constraints, including the relic abundance, direct detection limits, and indirect detection bounds. We find that the scalar triplet with $Y=0$ is ruled out by a combination of relic density, direct detection and indirect detection constraints. On the other hand, the scalar and fermionic triplets with $Y=2$ are both excluded by current direct detection experiments due to their large spin-independent scattering cross-sections. The viable parameter space of the remaining $Y=0$ fermion triplet dark matter lies within the projected sensitivity of near-future experiments, particularly those targeting indirect detection signatures. Collider prospects for these triplet extensions are also discussed in the Appendix.
We propose a new framework for solving the strong CP problem via a heavy axion, using mirror symmetry and grand unification. The mirror GUT sector remains unbroken and dynamically generates a calculable heavy mass scale via confinement without fine tuning. Models in this class feature a heavy axion, whose potential is less sensitive to Planck scale corrections, as well as a rich hidden sector from the confined mirror GUT. The solution to the strong CP problem remains unspoiled by the presence of additional phases in the GUT Yukawas, yet allowing the possibility of electric dipole moments within the reach of future experiments. Our proposal offers new directions in GUT model building, axion phenomenology, dark matter and cosmology.
Recently observed empirical scaling laws describe the performance of foundation-type models as three independent key quantities -- dataset size, compute, and model parameters -- are modified. Extracting these scaling laws informs the training of large complex models for which the tuning of hyperparameters in traditional ways is not feasible. This work for the first time explores if scaling laws can also be observed for the task of particle jet generation -- both relevant as a pre-training objective for foundation models and as in-situ simulation by itself. We indeed replicate the key logarithmic scaling law behavior for model-size scaling. Beyond studying the next token prediction validation loss of the generative model, we also study the sliced Wasserstein distance of five physical quantities that are not immediately available to the model during training. Our study shows that this quantity is monotonically related to the next token prediction validation loss, meaning that this loss is indeed a good proxy for the physics performance. For the scaling with dataset size and compute, we observe substantially weaker scaling behavior of both the loss and the sliced Wasserstein distance. We analyze this behavior by introducing the concept of a learnable window, and argue that autoregressive next token prediction on jet constituents exhibits comparatively rapid saturation relative to language-model studies. We discuss possible origins of this behavior, including the stochastic nature of QCD radiation and differences between generative and supervised learning tasks in collider physics.
We introduce FunKit, a Mathematica package for the derivation and tracing of functional equations from arbitrary master equations. FunKit provides an expression vocabulary and a set of rules that allow for derivations in any given field theory and master equation. It also allows users to add extensions for more specific equation systems. Therefore, it can be used in a wide range of situations, for example Dyson--Schwinger or functional RG equations, flowing reparametrisations, nPI equations, (modified) STIs and WTIs, functional Polchinski and Wegner flows, functional master equations with sources, and many others. Besides interfacing with the \FORM language to trace large tensor expressions efficiently, FunKit also provides facilities to export arbitrary Mathematica expressions to C++, Julia or Fortran code, including the results of derivations, which can then be evaluated numerically. Both the tracing and code generation can also be used independently and in combination with other packages.
Energy correlators exhibit well-understood scaling behavior in the collinear limit, governed by perturbative QCD dynamics. We explore how this scaling regime is broken by new physics, converting precise energy correlator measurements into a broadband search for new physics. Under generic assumptions, unitarity and positivity are sufficient to classify and constrain the relevant signatures, which imprint an angular resonance on top of this smoothly scaling background. This converts the search into bump hunting within jets. As a proof of principle, we derive projected LHC sensitivity for a light hadrophilic $Z'$, finding competitive constraints with existing searches.
We propose a gauge-symmetry origin of a rank-two Dirac neutrino mass matrix that enforces one exactly massless neutrino, while being consistent with the oscillation data, as well as cosmological constraints. The mechanism relies on a minimal dark $SU(2)_D$ gauge symmetry under which one right-handed-neutrino-like Weyl fermion is charged, thereby forbidding its Standard Model Yukawa coupling. Quantum consistency then fixes the minimal dark-sector completion: Cancellation of the Witten anomaly requires a second fermionic $SU(2)_D$ doublet, while a discrete $Z_4$ symmetry that forbids Majorana masses allows the two dark doublets to form a vectorlike pair. This anomaly-free completion gives rise to a secluded, confining dark sector with a viable dark matter candidate, linking the protected neutrino texture to dark infrared dynamics.
The existence of a landscape of metastable vacua raises the possibility that our Universe may have undergone quantum vacuum decay at late times. This work explores how such a transition can be tested with cosmological observables, focusing on precision distance measurements and cosmic microwave background (CMB) anisotropies. A set of phenomenological models is constructed in which late-time quantum tunneling changes the vacuum energy and may convert a subcomponent of dark matter into dark radiation, possibly accompanied by domain-wall production. The resulting expansion histories are compared with DESI DR2 baryon acoustic oscillation data; supernova distance measurements from DES-Dovekie, Pantheon+, and Union3; and a compressed CMB likelihood. For quantum-tunneling models, current cosmological distance measurements still allow a 50% decrease in the total vacuum energy for a transition redshift $z_t<1$. The model with dark-matter conversion and domain-wall production provides a good fit to resolve the tension between cosmological observables and the $Λ$CDM model, with a preferred transition around $z_t \sim 7$ and about 10% of dark matter participating in the transition. Additionally, CMB anisotropy constraints from bubble nucleation and the associated domain-wall network are derived and shown to strongly restrict slow or sparse late transitions. Applied to the minimal quantum-tunneling model, these constraints allow an $\mathcal{O}(10\%)$ decrease in the total vacuum energy for a transition redshift $z_t$ of order unity. For nonminimal models, dark-matter-density-dependent tunneling can proceed rapidly enough to evade such bounds. These results demonstrate that late-time quantum vacuum decay is a testable cosmological phenomenon and provide a concrete observational handle on metastable-vacuum physics motivated by landscape scenarios.
Low threshold cryogenic calorimeters are a key technology for the advancement of rare-event searches. However, since a few years their sensitivity reach is challenged by the presence of a rising spectrum at low energies named low-energy excess (LEE), ascribed to an unknown background. In this work, we describe the LEE as absorber events induced by the relative thermal-contraction coefficient mismatch between the absorber and the SiO$_2$ amorphous layer underneath the transition-edge sensors (TESs), present for example in the case of CRESST detectors. The relative contraction in processes with temperature changes, such as during sensor fabrication and cooldown from room temperature to the temperature of operation, can induce surface dislocation nucleation. Other interfaced materials with thermal-expansion mismatch can also generate dislocations during temperature-variation processes. We formulate a simple elastic model to bridge this solid-state effect and the LEE observations. Double-TES modules have been designed to provide surface background rejection. We highlight that the presence of the LEE in the coincident event band of double-TES modules does not exclude the explanation given in this work. Exemplary, we discuss the role of the thermal boundary resistance between absorber and sensor as explanation for the presence of the LEE in the coincident-event band. We propose detector designs to test these hypotheses and mitigate the LEE.
We present a lattice QCD determination of the Mellin moments of the unpolarized gluon parton distribution function in the proton. The analysis is based on matrix elements of nonlocal gluon operators coupled to momentum-boosted proton states. The calculation relies on an $N_f=2+1+1$ ensemble of maximally twisted mass fermions with clover improvement and the Iwasaki-improved gauge action, at a pion mass of approximately 260 MeV. Working within the short-distance operator product expansion (OPE) of the reduced gluon Ioffe-time distribution, we extract ratios of higher-order gluon moments, $\langle x^n\rangle$ with $n>1$, to the gluon momentum fraction, $\langle x\rangle$. We investigate systematic effects associated with the truncation of the order of moment in the OPE, the choice of minimum and maximum Wilson-line separations entering the analysis, and the treatment of mixing with the quark-singlet under perturbative matching. The stability of the extracted moments is further studied under scale evolution using DGLAP equations, allowing us to assess uncertainties related to perturbative truncation by varying the scale. Our work provides a determination of the ratio $\langle x^3\rangle_g/\langle x\rangle_g$ at a scale of 2 GeV, with uncertainties that account for both statistical and the dominant theoretical systematic uncertainties.
We present results on the Mellin moments of the unpolarized parton distribution function (PDF) of the pion and kaon up to the fourth order. The computation is done using one $N_f=2+1+1$ gauge ensemble of twisted mass fermions with quark masses tuned to approximately their physical values. We reconstruct the valence pion and kaon PDFs using the connected contributions to the three Mellin moments. We compare our results on the Mellin moments and the reconstructed PDFs with other lattice QCD and phenomenological determinations.
We study the quantum chaos bound in the photon ring region surrounding black holes. By evaluating the Lyapunov exponent associated with unstable null geodesics in a broad class of generalized Kerr geometries, as well as the temperature induced by a string probe, we show that the quantum chaos bound is exactly saturated on equatorial circular orbits of the photon ring. We confirm our result by deriving the same exponent from out-of-time-order correlators in the near ring region. As a byproduct, we show that the photon ring saturation of the quantum chaos bound implies the saturation of the Bekenstein bound on the rate of information emission from the ringdown phase through the quasi-normal modes in the eikonal limit. Our results extend the known correspondence between black hole thermodynamics and chaotic dynamics, highlighting the role of the photon ring as a probe of the fundamental limits on thermalization and information scrambling in black holes.
We present a lattice QCD calculation of the charged-neutral pion mass splitting $M_{π^+} - M_{π^0}$ at $\mathcal{O}(α_\mathrm{em})$ using a recently proposed framework based on a Pauli-Villars (PV) regulated photon propagator defined in the continuum and infinite-volume limit, with $Λ$ acting as an additional UV cutoff scale. The use of this propagator avoids power-law finite-volume effects, allowing for a straightforward treatment of the infinite-volume limit. We perform the calculation using CLS ensembles, studying finite-volume effects, the continuum limit and the extrapolation to the physical point for several values of the scale $Λ$. By means of the Cottingham formula, we further decompose the result into elastic and inelastic contributions at fixed $Λ$. Our final result, after removing the cutoff scale $Λ$, is $M_{π^+} - M_{π^0} = 4.56(22)$ MeV, in good agreement with the experimental measurement. This calculation serves as a validation of the formalism in a well-controlled setting and offers useful insights into the application of electromagnetic corrections to other observables.
Magic relations are a class of integration-by-parts identities where all integrals in the generating sector drop out. Since their presence causes several otherwise successful methods in the Feynman-integral computational pipeline to break down, they are important to detect and understand. In this paper, we take a first step toward a systematic characterization of such identities. Specifically, we observe and argue that the occurrence of magic relations always coincides with the presence of higher-dimensional critical varieties in the generating sector. This provides a practical computational test to check if a family of Feynman integrals can contain magic relations and to find them, which we implement in the ancillary Mathematica file Magic-Test.m. Additionally, we discuss how to count the number of master integrals in the presence of higher-dimensional critical varieties, classify the behavior of magic relations under symmetries, and we discuss their interplay with cuts.
Bayesian analyses in the context of relativistic heavy-ion collisions have so far relied almost exclusively on bulk hadronic observables constructed from momentum degrees of freedom to constrain the transport properties of the quark-gluon plasma. In this work, we perform the Bayesian inference after incorporating the longitudinal spin polarization of $Λ$ hyperons alongside conventional bulk measurements in Pb+Pb collisions at $\sqrt{s_{NN}}=5.02$ TeV to constrain the shear and bulk viscosity to entropy density ratios, $η/s$ and $ζ/s$. We demonstrate that the inclusion of spin polarization, which provides complementary sensitivity to the space-time structure and vorticity of the medium, shifts the posterior distribution of $ζ/s$ toward larger values, although current uncertainties do not allow a statistically significant separation at the 68% credibility level. Nevertheless, the results establish spin polarization as a valuable probe in quantitative studies of QGP transport properties and indicate that it should be incorporated in comprehensive and systematically constrained Bayesian extractions of the medium's dynamical parameters.
We consider systems characterized by the presence of a rapidly oscillating force. A general method is presented for the construction of the effective action governing the large-scale nonlinear dynamics of such systems order by order in inverse powers of the oscillation frequency $ω$. The explicit expression for the effective Lagrangian is derived up to ${\cal O}(1/ω^6)$ next-to-next-to-leading approximation. The general structure of the high-frequency expansion reveals a broad class of nonlinear systems whose transition curves are identical to those of the linear Mathieu equation, which enables a fully nonperturbative stability analysis in the case of strong driving and nonlinearity. The method is generalized to velocity-dependent forces and configuration space with curvature, characteristic to systems with constraints. Several applications are discussed in detail, including the dynamical magnetic trapping of electric charges.