2026-11-10 | | Total: 82
At future high-energy colliders, such as a multi-TeV muon collider or the high-energy stages of an $e^+e^-$ linear collider, collinear radiation of electroweak gauge bosons becomes a dominant phenomenon. In this regime, many key processes are primarily initiated by photons and gauge bosons originating from lepton splittings. Electroweak Parton Distribution Functions (EW PDFs) provide an appropriate framework to describe these interactions by resumming large collinear logarithms arising from all possible Standard Model splittings, thereby enabling a partonic interpretation of the hard scattering. The formalism is particularly suitable for studying new physics beyond the Standard Model in the electroweak sector at higher scales. In this paper, we present the implementation of the EW PDF formalism within the Monte Carlo event generator \texttt{Whizard}. We demonstrate several applications of this framework and analyze its interplay with full fixed-order calculations as well as the purely kinematic Effective Weak Boson Approximation. We highlight scenarios where the EW PDF approach is especially useful for phenomenology and physics-potential studies for future lepton colliders. The EW PDF formalism is a vital step towards the ``electroweak jet regime'' where electroweak particles behave similarly to QCD jets.
Higher-form symmetries can protect extra-dimensional axions, but additional symmetry breaking can still spoil the solution to the strong CP problem. We focus on electrically charged matter and magnetic gauging in a five-dimensional circle compactification where the axion arises from an Abelian gauge field and the radius is dynamical. We discuss two mechanisms for stabilising the radius/radion using Casimir energy from fields neutral under the Abelian gauge field, with and without an additional periodic scalar winding around the compact circle. Both admit Minkowski and de Sitter minima. At the Minkowski minima (with radion masses of $100$--$170\,\mathrm{TeV}$), electrically charged matter induces a residual strong CP angle of order $10^{-17}$ and shifts the radius negligibly. Magnetic gauging remains compatible after including its additional degree of freedom in the Casimir spectrum, provided its coupling is sufficiently small. We observe that increasing the radius suppresses charged-matter breaking relative to QCD but increases magnetic breaking contribution. We also determine the radion decay channels, including the graviphoton, and study reheating from a possible radion-dominated era without scalar winding or magnetic gauging. For a non-minimal Higgs-curvature coupling $ξ=7$, the visible-sector reheating temperature is about $85\,\mathrm{MeV}$ with a dark-radiation contribution of $ΔN_{\rm eff}\approx0.02$. The entropy released during radion decay dilutes the axion abundance, allowing axions with $f_a=10^{14}\,\mathrm{GeV}$ and an initial misalignment angle $\vartheta_i\approx0.43$ to account for the entire observed dark matter abundance. For the same decay constant, reproducing this abundance in radiation domination requires $\vartheta_i\approx0.065$.
We examine a possible signature of the neutrino angular-momentum contribution to the muon magnetic anomaly measurement in the J-PARC E34 experiment. In previous papers we argued that the angular momentum carried by the unobserved neutrino pair may leave an uncertainty in the anomaly inferred from the time modulation of detected positrons. At Fermilab the highly relativistic muons make the forward-emission approximation suitable for this estimate. The lower muon momentum at J-PARC requires the decay kinematics to be treated explicitly. We integrate the longitudinal momentum of the neutrino pair using the positron energy distribution and asymmetry derived from the theoretical weak-decay spectrum, and the same event-counting weight $N(E)A(E)$ is employed for both experiments. For positron energies between $1.5$ and $2.9$ GeV at Fermilab and between $200$ and $275$ MeV at J-PARC, the proposed correction is smaller at J-PARC by $0.415\pm0.176$ ppm, when the most optimistic E34 precision is assumed. The expected difference amounts to $2.36$ standard deviations if the future measurement is centred on this value. The calculation thus identifies an experimental comparison which may test the proposed effect independently of the hadronic contribution to the Standard Model prediction. We also examine the possibility of testing the effect by comparing measurements at two different muon momenta. With a fixed positron energy threshold of 200 MeV, the maximum sensitivity is obtained at a muon beam energy of 245 MeV. The predicted difference relative to the nominal J-PARC configuration has significance of $4.88σ$ for the most optimistic experimental sensitivity.
We extract collinear transversity distributions and dihadron fragmentation functions in a joint neural-network analysis of dihadron production in electron-positron annihilation, semi-inclusive deep-inelastic scattering, and polarized proton-proton collisions. We incorporate lattice QCD calculation of isovector transversity, which constrains the up-down flavor differences for quarks and antiquarks over broad momentum fraction range complementary to data. Matched fits with and without the lattice input describe the experimental data comparably well. Including the lattice input shifts the down-valence distribution toward more negative values and enhances up-valence transversity at large $x$. The magnitude of the central down-valence result exceeds the reference Soffer inequalities, motivating further test of positivity constraints.
We extract the relative phase between strong and electromagnetic amplitudes of the decays of a vector charmonium into spin-$1/2$ baryon-antibaryon pairs. Starting from a flavour-SU(3) effective Lagrangian with vector and tensor currents, we construct a parametrization in terms of partial waves coefficients. The electric and magnetic phases, the phase of the integrated complex overlap and the effective rate angle are derived from the same source amplitudes. We retain two possible solutions and both of them describe the scalar $Σ^+\barΣ^-$ scan phase within its uncertainty, while presenting tension with several electric-magnetic production phases.
We present a basis of special functions and an efficient numerical implementation for the planar two-loop master integrals that enter the leading-colour amplitudes for $pp\to t\bar{t}H$ with a closed light-quark loop. We provide analytic results for the special functions at transcendental weight one and two in terms of logarithms and dilogarithms, while higher weights are evaluated as one-fold integrals. We focus on the nested roots that arise, studying their analytic structure on the physical region and constructing numerically stable representations of letters and $\mathrm{d}\log$ forms that involve them. We find that a number of $\mathrm{d}\log$ forms have an antiderivative that cannot be represented by a single logarithm over the physical region and tackle them numerically with adaptive Chebyshev approximation. We implement the functions in \texttt{C++} as an extension of the \texttt{PentagonFunctions++} library. We study the performance of our implementation over physical phase-space points and find that the evaluation is accurate and fast enough for phenomenological application. This provides an important building block for NNLO studies of $t\bar{t}H$ production at the LHC.
An angular analysis of the purely baryonic decays $Λ_b^0\to np\bar p$ and $Λ_b^0\toΛp\bar p$ is presented within a phenomenological generalized-factorization model. Within this model, two current contributions leading to the same final state are combined coherently. Pair-angle distributions, forward-backward hemisphere asymmetries, and Legendre moments are calculated in both channels, together with the polarization of the produced $Λ$ for an unpolarized parent. For the neutron channel, the full-region forward-backward asymmetries are $4.46\times10^{-4}$ in the $n\bar p$ frame and $0.856$ in the $p\bar p$ frame at $N_c^{\mathrm{eff}}=3$ and zero additional relative phase. From mass-bin integration, the small first asymmetry is found to arise from compensation between signed contributions from different regions. In the $Λ$ channel, the angular distributions are modified by the selection $m(p\bar p)<2.85\,\mathrm{GeV}$, and the rate-weighted longitudinal polarization along the $Λ$ flight direction is shifted from $-0.803$ to $-0.954$. These results show how angular and spin observables depend on the relative current contributions and the selected phase-space region. The model-dependent analysis is intended as a first step toward understanding the dynamics of these decays and as a basis for further theoretical studies and comparisons with measurements.
We study a generalized NMSSM in which a gauge singlet links the Higgs bilinear to the type-II seesaw scale. The singlet vacuum expectation value contributes to $μ_{\rm eff}$ and generates the supersymmetric mass of a vectorlike triplet pair. Integrating out the triplets gives the Weinberg operator, and their gauge charges admit a possible messenger interpretation. Large cancellation between the bare $μ_0$ term and the singlet VEV generated bilinear Higgs term that keeps $μ_{\rm eff}$ near the soft scale is still natural in the supersymmetric limit (by the supersymmetric non-renormalization theorem), although the concrete value of $μ_{\rm eff}$ also depends slightly on the soft SUSY breaking parameters and finite threshold matching. The studied bino-like and singlino-like neutralino solutions accommodate triplet input scales approaching $10^{15}$ GeV together with electroweak-to-TeV dark matter. Representative spectra consistent with heavy-Higgs funnels or coannihilation with gluinos or winos yield allowed relic abundances while satisfying direct-detection constraints.
The observation of suppression and finite high-$p_\perp$ elliptic anisotropy in O+O and Ne+Ne raises the question of whether a common quark--gluon plasma (QGP) description can simultaneously account for both observables across systems from oxygen to lead. We present predictions using event-by-event bulk evolution constrained by soft-sector observables and coupled to DREENA radiative and collisional energy loss, with the same microscopic parameter set for O+O, Ne+Ne, Xe+Xe, and Pb+Pb. The calculation reproduces the observed ordering of suppression and predicts elliptic anisotropies of comparable magnitude across the four systems in 10--20$\%$ collisions. We predict a modest centrality dependence of $v_2^{\mathrm{SP}}$ and a more pronounced decrease of $v_3^{\mathrm{SP}}$ toward peripheral light-ion collisions, in contrast to the corresponding trends in larger systems. In light ions, we predict that $v_2^{\mathrm{SP}}/(1-R_{AA})$ is nearly independent of momentum at high $p_\perp$. We further predict that the differential ratio $v_2\{4\}/v_2\{2\}$ depends only weakly on momentum at high $p_\perp$, with magnitudes and centrality trends broadly comparable to those of the measured soft-sector ratios, suggesting that the soft-sector fluctuation pattern persists in hard--soft correlations. Predictions for D and B mesons extend the comparison to probes of different flavors. Taken together, these system-size, harmonic, fluctuation, and flavor dependences provide experimental tests of how far a common QGP description extends.
We compute the magnetic moments of the $Σ_c$ and $Σ_b$ isotriplets with external-field QCD sum rules. The heavy-quark propagator is kept with its full mass dependence and with its couplings to the photon, to soft gluons and to both simultaneously, and the magnetic sum rule is normalized by the mass sum rule of the same interpolating field, which removes the residue and gives the baryon masses within a few percent. The quark-condensate magnetic susceptibility $χ$, which otherwise dominates the uncertainty, is fixed by the proton magnetic moment computed with the same interpolating field; the neutron moment then follows to about $5\%$. Uncertainties from the QCD inputs are evaluated with a Monte Carlo analysis, and those from the interpolating field, the Borel window, the scale of the heavy-quark mass, the calibration and the normalization are added. With $\overline{\rm MS}$ heavy-quark masses we obtain $μ_{Σ_c^{++}}=2.37\pm0.65$, $μ_{Σ_c^{+}}=0.50\pm0.15$, $μ_{Σ_c^{0}}=-1.39\pm0.37$, $μ_{Σ_b^{+}}=2.16\pm0.65$, $μ_{Σ_b^{0}}=0.55\pm0.16$ and $μ_{Σ_b^{-}}=-1.06\pm0.32$ in nuclear magnetons; at this order the results depend sizably on the heavy-quark mass scheme. The moments are carried by the light diquark; the heavy-quark contribution, whose separation depends on the interpolating field, is $-0.13\pm0.02\,μ_N$ for $Σ_c$ and $+0.013\pm0.003\,μ_N$ for $Σ_b$.
We consider the possibility that electroweak naturalness may be epoch dependent. In no-scale supergravity, a rolling quintessence modulus can control the supersymmetry-breaking scale, allowing superpartners to be much lighter in the early universe while remaining heavy today. Fifth-force and varying-constant effects can be suppressed via correlated modular weights in the visible-sector Kähler metric and modulus-dependent superpotential Yukawa couplings, while radiative stability of the quintessence potential requires additional tuning of the field-dependent vacuum energy. The resulting effective theory has distinctive implications for early-universe phase transitions, relic abundances, and gravitational waves.
We study negative-parity hidden-flavor pentaquarks with three to five heavy quarks, $QQ\bar{Q}qq$, $QQQ\bar{Q}q$, and $QQQQ\bar{Q}$ with $Q=c,b$ and $q=u,d$, in QCD sum rules, using the complete set of S-wave baryon--meson currents with $J^P=1/2^-$, $3/2^-$, and $5/2^-$ and the operator product expansion up to dimension 10, 9, and 6. We focus on the splittings between states of the same baryon--meson configuration, which, unlike the absolute masses, are only weakly affected by the choice of the heavy-quark mass scheme and by a finite width. Configurations whose light degrees of freedom have a unique spin are nearly degenerate or show splittings that decrease from charm to bottom by a factor compatible with $\bar m_c/\bar m_b$, whereas configurations with two light quarks in different hadrons show splittings that are not reduced, as expected from heavy-quark spin symmetry. A decomposition of the correlators into factorizable two-hadron and connected quark-exchange contributions shows that the correlators are dominated by the two-hadron part and that this pattern holds separately for both parts: the quark-exchange contribution decreases from charm to bottom when a heavy quark is exchanged and is not reduced when a light quark is exchanged. In the fully heavy systems the $1/2^-$ and $3/2^-$ states are degenerate within about $20$~MeV, while the $5/2^-$ state lies $126$--$139$~MeV ($cccc\bar c$) and $183$--$228$~MeV ($bbbb\bar b$) above the $3/2^-$ state. This splitting is generated by the exchange terms between the identical heavy quarks and does not decrease from charm to bottom, which points to the Pauli principle rather than to a hyperfine effect. Currents whose baryon and meson share no quark flavor are exactly factorizable at this order and cannot distinguish bound states from two-hadron states.
Recent JWST observations motivate mechanisms for rapid black-hole growth in the early Universe. Extending our previous Schwarzschild study, we show that critical relativistic Bondi accretion of repulsively self-interacting dark matter remains an efficient super-Eddington channel for Kerr black holes with rotating inflow. Within a stationary equatorial quasi-spherical framework, smooth critical solutions exist over a finite range of inflow angular momentum set by the Kerr marginally bound capture values. Rotation changes the sonic structure and accretion efficiency, while the critical density cancellation and dark-sector/halo matching retain their Schwarzschild form. The resulting rapid growth persists until the finite halo supply becomes limiting. Because the inflow also carries angular momentum, stellar-mass seeds can grow and spin up simultaneously; illustrative population evolutions therefore provide a joint mass--spin route to early supermassive Kerr black holes.
The next-to-leading-order (NLO) QCD corrections have been investigated for the production of doubly heavy baryons via the inclusive $W^{+}$-boson decays $W^{+}\toΞ_{cc}+\bar{c}+\bar{s}+X$ and $W^{+}\toΞ_{bc}+\bar{b}+\bar{s}+X$ within the nonrelativistic QCD framework. The decay widths, the kinematic distributions, and the uncertainties induced by the heavy-quark masses and the renormalization scale have been analyzed for both the color-antitriplet ($\bar{\mathbf{3}}$) and the color-sextet ($\mathbf{6}$) $S$-wave states. With the central input parameters $m_{c}=1.8~\mathrm{GeV}$ and $m_{b}=5.1~\mathrm{GeV}$, the $K$ factors are $2.27$ for $Ξ_{cc}$ and $2.95$ for $Ξ_{bc}$, i.e., the NLO corrections increase the corresponding LO widths by about $127\%$ and $195\%$, which shows that the NLO corrections are indispensable for quantitative predictions. The individual-channel $K$ factors range from $1.94$ to $5.53$, with the largest one occurring for $Ξ_{bc}$ in the $[{}^{3}S_{1}]_{\mathbf{6}}$ channel. The distributions show that the NLO corrections enhance the low invariant-mass region and shift the baryon energy spectrum toward smaller energy fractions. With $Γ_{W}=2.085~\mathrm{GeV}$, the inclusive branching ratios are $8.86\times10^{-5}$ and $3.63\times10^{-6}$ for $Ξ_{cc}$ and $Ξ_{bc}$, respectively, which provide an NLO benchmark for future experimental searches.
Taking the real singlet scalar dark matter model as a representative test case, we examine the constraints imposed on the model's parameter space by Higgs trilinear coupling measurements and perturbative unitarity requirements. For the low dark matter mass regime, the constraining strength delivered by the current highest-precision $κ_3$ measurement is notably inferior to the exclusion limit derived from perturbative unitarity requirements, while the perturbative unitarity bound itself remains considerably less restrictive than the constraints extracted from the extant experimental upper limit on the Higgs invisible decay branching ratio $\mathrm{BR}_{\mathrm{inv}} \leq 0.107$. Both the coupling modifiers $κ_3$ and $κ_4$ exhibit a well-defined decoupling behavior, with their respective magnitudes monotonically decreasing as the dark matter mass $m_S$ increases. As a direct consequence, the corrections to both couplings asymptotically approach zero for dark matter masses in excess of 60 TeV, when the parameter $a_2$ remains within its perturbative validity domain. Neither existing collider measurements of the Higgs trilinear coupling nor measurements targeted for planned future collider facilities possess the requisite sensitivity to impose more stringent constraints on the model's viable parameter space than the well-established bounds derived from experimental probes of dark matter relic density, invisible Higgs decay signatures, dark matter direct detection, and dark matter indirect detection. Nevertheless, under the scenario that the experimental measurement uncertainty of $κ_3$ reaches the $\mathcal{O}(0.1\%)$ scale or lower, non-negligible constraints acting on a subset of the large mass parameter space associated with the WIMP dark matter scenario can be expected.
In these proceedings we review the features of the recently-released version 5.0 of the FORM computer algebra system, which is one of the core software packages for the large-scale symbolic manipulations required by calculations in theoretical high-energy physics. We also discuss developments since the release of version 5.0, focussing in particular on recent performance improvements.
We investigate how the late-time decays of unstable particles into dark matter and neutrinos can leave measurable imprints on both cosmological observables and collider experiments. We study this mechanism within the scotogenic model, where the lightest singlet fermion $N_1$ serves as the dark matter candidate and is produced via the late decay of the lightest inert doublet scalar via the superWIMP mechanism. By incorporating a kination-dominated era in the early Universe, the altered expansion rate allows to match the observed $N_1$ relic density in parameter regions that would be underabundant in standard cosmology. The energy exchange between the injected neutrinos and photons resulting from the late-time decay can produce CMB spectral $μ$-distortions, detectable shifts in the effective number of relativistic species $N_{\rm eff}$, or both simultaneously. Concentrating on the lightest component of inert doublet scalars with masses between $200$ and $500~\text{GeV}$, we find that these cosmological signals are accompanied by complementary signatures accessible at the LHC. Although our analysis is derived within the scotogenic framework, our findings apply more broadly to scenarios in which the unstable particle undergoes late-time out-of-equilibrium decay to produce dark matter and non-thermal neutrinos, offering how future CMB measurements and collider searches can jointly probe dark matter production and the early Universe's thermal history.
A fundamental requirement for the consistency of quantum field theories is that their Higgs potential is lower-bounded, such that no direction in the $Φ$ field space leads to $V(Φ)\to \, -\infty$. This work presents a comprehensive analysis of the various theoretical constraints on potentials in extended Higgs sectors, highlighting their relative impact and restrictiveness. We illustrate our study on the potential parameters with two particular Higgs sector extensions: the Real Singlet Extension (RSE) and the Two-Higgs-Doublet Model (2HDM), separating $\mathbb{Z}_{2}$-conserving and $\mathbb{Z}_{2}$-violating cases. In order to do this, we derive and compare these conditions to other studies found in the literature. Additionally, we developed a publicly available code capable of testing positivity conditions in complex potentials, including particularly challenging cases such as the general $\mathbb{Z}_2$-violating 2HDM. Finally, we apply those restrictions to the physical parameters and deduce an upper and a lower bound on the masses of new heavy scalars.
We provide a continuously verified diagnostic that distinguishes a compatible single-pole fit from a mass region required by QCD sum-rule moments. For specified glueball sum-rule inputs, the moments can require contributions in broad mass bands even when narrower candidate bands are empty. We compare fourteen two- and three-gluon currents at central fixed-scale inputs, undeformed continua and 1\% matching tolerance. All twelve nonempty inputs permit an empty candidate band, yet nine still require a broader contribution. This separates two statements that a pole fit conflates: a successful pole fit does not make its mass region necessary, and excluding every single pole does not remove the information that a broader region must contribute. The five diagnostic classes of the survey organize these outcomes by whether a pole is compatible and whether the 30\% band is required or can be empty; two inputs admit no nonnegative spectrum. Four of the nine positive bounds persist on subwindows satisfying the fixed-scale input criteria. For the $GG,0^{-+}$ current, closest to the nominal $X(2370)$ mass, the 20\% band containing $X(2370)$ can be empty while the 30\% band requires about 44.2\% of each spectrum's own low-energy Borel integral. Determining the physical nature of $X(2370)$ requires further theoretical and experimental progress.
We investigate photon-ALP conversion as a simultaneous explanation of the ultra-high-energy photons observed from GRB 221009A by LHAASO and of the $\sim300$ TeV photon reported by Carpet-3. We model the magnetic field of the host galaxy as a set of domains with randomly oriented field directions and variable domain lengths, and propagate the photon-ALP system through the host galaxy, the intergalactic medium and the Milky Way. Fitting the combined LHAASO and Carpet-3 spectra, we find that a broad distribution of domain lengths raises the photon survival probability at the highest energies and accommodates the Carpet-3 point with photon-ALP conversion alone. The best fit, $m_a=3.85\times10^{-7}$~eV and $g_{aγ}=3.07\times10^{-11}\,\mathrm{GeV}^{-1}$, lies below the CAST bound, and parts of the $1σ$ and $2σ$ regions are also compatible with it.
The spin correlation of $Λ\barΛ$ hyperon pairs produced in heavy-ion collisions is investigated through scalar and pseudo-scalar mesonic channels that mediate the creation of $s\bar{s}$ pairs. Because a $J=0$ meson forces $L=S$ while parity requires $(-1)^{L+1}=+/-1$, the two vertices project the pair onto distinct total-spin states: a scalar ($0^+$) vertex selects the $L=S=1$ configuration with a spin weight of $+1$, whereas a pseudo-scalar ($0^-$) vertex selects the $L=S=0$ configuration with a weight of $-3$; these weights are the traces of the correlation tensor, corresponding to $P_{s\bar{s}}=+1/3$ and $-1$, respectively. Implementing this framework within the three-flavor Nambu-Jona-Lasinio model at the Random-Phase-Approximation level, we compute the scalar and pseudo-scalar spectral weights. Our results show that the pseudo-scalar channel dominates, resulting in a negative net correlation fraction throughout. And the spin correlation also depends on the opening angle: strongest spin correlation at around $θ=0$ while weakest at $θ=π$. Furthermore, we show that the collision-energy dependence of the momentum-integrated spin correlation along the freeze-out line is sensitive to the trajectory's path relative to the chiral phase boundary. Trajectories passing near the crossover region exhibit a non-monotonic structure, providing a sensitive probe for the location of the chiral phase transition.
We investigate the stability of homogeneous chiral condensates in magnetized quark matter within the two-flavor Nambu-Jona-Lasinio model. By expanding the effective action to quadratic order, we derive the full momentum-dependent bosonic two-point function in a constant magnetic field. At zero temperature, our analytical results show that each occupied Landau level produces a sharp logarithmic divergence at its corresponding Fermi momentum. This severe Peierls instability is driven by one-dimensional Fermi surface nesting within the effectively 1+1-dimensional subsystems created by Landau quantization. Numerically, we map out the momentum-space topography of the two-point function across varying magnetic fields and temperatures. We find that strong magnetic fields induce deep negative canyons. While thermal fluctuations smooth these canyons into shallower valleys as the temperature increases, the negative regions remain robust at intermediate temperatures and only vanish entirely at sufficiently high temperatures. Intriguingly, at larger chemical potentials where the system is strictly stable, a broad positive basin with a flat curvature at the origin persists at non-zero momentum. This non-monotonic dispersion is the hallmark of a moat regime, acting as a precursor state that strongly signals a residual tendency toward inhomogeneous chiral symmetry breaking. Our findings provide valuable insights into the phase structure of dense, magnetized quark matter.
The microscopic origin of the deuteron tensor polarized structure function $b_1$ remains unresolved. Conventional nuclear models attribute a nonzero $b_1$ mainly to the $D$-wave component of the deuteron wave function, but generally fail to reproduce the magnitude and momentum-dependence observed by the HERMES experiment. Using a QCD-based light-front framework with degrees of freedom beyond the conventional nucleonic picture, in particular including hidden-color components and relativistic spin--momentum correlations, we find that a small but nonzero estimated hidden-color admixture can be sufficient to generate a substantial tensor response ($b_1(x)$), while it is essentially insensitive to the $D$-wave component. The physical $D$-wave probability of $0.85\%$, fixed by the empirical quadrupole moment, changes $b_1$ by only $1.7\%$, with the dominant contribution arising from the colored-cluster $S$-wave sector. These results demonstrate that the tensor-polarized structure of the deuteron is driven primarily by hidden-color dynamics rather than by the $D$-wave component.
We study the connection between asymmetric dark matter and nanohertz gravitational waves from a supercooled first-order phase transition in a classically conformal dark sector. The dark sector contains a real scalar, a complex scalar and a Dirac fermion, whose mass arises from radiative symmetry breaking. An additional Affleck--Dine condensate generates a fermion asymmetry and heats the dark bath before nucleation. The conserved dark charge fixes the asymmetry between fermions and antifermions, while entropy production during delayed visible-sector reheating dilutes the relic and determines the redshift of the gravitational-wave spectrum. We calculate the thermal nucleation rate in a three-dimensional effective theory with finite-density corrections, next-to-leading-order contributions and fluctuation determinants (NLOdet). We follow the coupled thermal evolution of the dark and visible sectors. For a benchmark chosen to match the observed dark-matter abundance, the transition completes at MeV-scale temperatures and produces a nanohertz signal under the adopted dynamical and thermal assumptions. We compare the predicted collision-dominated gravitational-wave spectra with free-spectrum posteriors from pulsar timing arrays. For a benchmark point and variations of the dark-sector parameters, the gravitational-wave spectra and dark matter relic abundances are compared, and compatibility with Big Bang nucleosynthesis is assessed.
Sequential kinetic freeze-out in high-energy collisions probes the dynamics governing the final stages of strongly interacting matter. Inspired by the recently established quadruple freeze-out scenario, we present a thermodynamic analysis of transverse momentum spectra for a set of identified particles produced in proton-proton collisions at $\sqrt{s} = 2.76$ and $7$ TeV. Using a modified Tsallis-Pareto function, we extract the effective temperature and the non-extensivity parameter for each species. We observe a robust five-layer hierarchy in the effective temperature. At both energies, the layers are: non-strange hadrons, singly-strange hadrons, multi-strange hadrons, light nuclei, and the $φ(1020)$ meson. At $7$~TeV, the $η$ meson joins the $φ(1020)$ meson in the highest-temperature layer. At both energies, each pair of neighboring layers is separated by more than $ 3σ$, and the separation at the top exceeds $7σ$. The grouping persists at both energies, while the temperatures shift upward at 7 TeV. The observed pattern shows that the decoupling sequence is driven by an interplay of strangeness, cross-sections, and the OZI rule rather than by particle mass alone. Notably, the $φ(1020)$ meson has a mass close to the proton, but a significantly higher temperature due to its OZI-suppressed coupling to the hadronic medium. The unexpected coincidence of $η$ and $φ(1020)$ at 7 TeV presents a puzzle that may point to a common early production mechanism or in-medium effects, which are not yet understood. The five-layer hierarchy provides new constraints on theoretical models of hadronization. It suggests that the sequential decoupling of particles from the evolving medium is even more sensitive to their internal structure and interaction than previously assumed, with implications for the QGP-like matter formation in small systems.