2026-07-16 | | Total: 26
We present fixed-order perturbative calculations up to four-loop order for a generic non-global QCD observable in hadron-hadron collisions. Specifically, we study the invariant-mass distribution of the highest-$p_t$ jet produced in association with a vector boson or a Higgs boson, where jets are defined using the Cambridge/Aachen sequential recombination algorithm. This work is part of a series of papers~\cite{Khelifa-Kerfa:2015mma, Khelifa-Kerfa:2024hwx, Khelifa-Kerfa:2024dut, Khelifa-Kerfa:2025cdn, Khelifa-Kerfa:2024udm, Khelifa-Kerfa:2025jev} examining the impact of various jet algorithms on the perturbative structure of non-global observables across different collision environments. We find that at fixed order the Cambridge/Aachen algorithm suppresses the large non-global logarithms more effectively than the anti-$k_t$ and $k_t$ algorithms, while at all orders it performs comparably to the $k_t$ algorithm. Furthermore, comparisons of all-orders resummed form factors reveal that finite-$N_c$ corrections remain at the percent level, consistent with previous findings.
We show that the concept of primary dimensions, first introduced in [1] as an organizing principle for chiral Lagrangians, is inconsistent. Although this had been pointed out already in [2], the notion of primary dimensions has re-appeared in recent literature. We briefly comment on the proper power counting for such effective field theories, including the electroweak chiral Lagrangian with a light Higgs, which is based on chiral dimensions, equivalent to the counting of loop orders.
Parton Distribution Functions (PDFs) at large Bjorken-$x$ are mostly constrained by high-energy measurements, and thus risk to absorb the energy-growing effects of unaccounted new physics (NP) if it were present in the high-energy tails of hadron-collider observables used in global PDF fits. If undetected, such contamination biases the resulting PDFs and can potentially hide the absorbed NP signals from subsequent searches. We illustrate this effect in a Standard Model Effective Field Theory (SMEFT) risk-assessment scenario affecting high-mass Drell-Yan production at the High-Luminosity LHC, and compare two strategies to obtain robust PDFs and SMEFT bounds: conservative fits that exclude data above an energy cut, and simultaneous fits of PDFs and Wilson coefficients performed with the {\tt SIMUnet} tool. Both approaches are shown to successfully recover the injected new physics that is otherwise absorbed and hidden by the PDFs.
We investigate the production of high-energy neutrinos from microscopic black holes formed through the gravitational collapse of asymmetric dark matter accumulated inside neutron stars. When Hawking evaporation dominates over accretion, long-lived, feebly interacting particles beyond the Standard Model escape the neutron star and subsequently decay into high-energy neutrinos. We analyze the repeated cycle of dark matter capture, black hole formation, and evaporation, identifying two distinct regimes determined by the competition between the dark matter thermalization time and the collapse cycle. In particular, we identify a partially thermalized regime in which the dark matter cloud evolves toward a quasi-stationary state with a temperature significantly exceeding that of the neutron star core. We derive the time-integrated Hawking emission, the resulting secondary neutrino spectra, and the expected Galactic and diffuse extragalactic neutrino fluxes. The predicted signal exhibits two distinctive signatures: a broad neutrino spectrum with a characteristic energy scale set by the initial Hawking temperature of the evaporating black hole, whose spectral peak naturally lies above $\mathcal{O}(10\,{\rm TeV})$, and an extended Galactic component strongly concentrated toward the Galactic Center. Although the predicted event rates are generally small, the resulting signal may contribute at the percent level to the observed Galactic high-energy neutrino flux under favorable microscopic and astrophysical conditions. The proposed mechanism provides a new observational window on Hawking evaporation through microscopic black holes continuously produced inside neutron stars, linking dark matter, compact objects, black hole thermodynamics and high-energy neutrino astronomy.
We study polarization observables in $τ^- \to V_1^- V_2^0 ν_τ$ decays for $(V_1^-, V_2^0)=(ρ^-, ρ^0), (ρ^-, ω), (K^{*-}, ρ^0), (K^{*-}, ω), (ρ^-, \bar{K}^{*0})$ using the helicity formalism within the generalized hidden local symmetry framework. The predicted branching fractions are comparable with those from the effective chiral model, but smaller than those from the angular momentum algebra model. Our results for the longitudinal polarization fractions $f_L$ indicate that nonnegligible corrections to the endpoint limit ($f_L = 1/3$), a consequence of Lorentz symmetry, are present at the scale of tau mass. Future measurements of $f_L$ at Belle II for $τ^- \to V_1^- V_2^0 ν_τ$ decays will provide crucial information on polarization patterns observed in nonleptonic charm decays.
We study heavy right-handed neutrino production induced by a derivative coupling to an axion-like inflaton and its implications for non-thermal leptogenesis. We develop a unified treatment of Majorana-fermion production during inflation and preheating. During inflation, the rolling inflaton background generates a helicity-asymmetric spectrum that can be obtained analytically in the slow-roll regime, while during preheating the oscillating background drives repeated non-adiabatic production events with a helicity-dependent effective momentum and a natural momentum cutoff. We compute the resulting abundance in both the delayed-decay and prompt-decay limits, including the effects of Pauli blocking and successive production stages. We further clarify the relation between the fermion basis used to identify non-adiabatic production and the Hamiltonian-diagonal basis used to define instantaneous occupation numbers. Applying these results to non-thermal leptogenesis, we identify parameter regions consistent with the observed baryon asymmetry.
Electron-positron annihilation is an ideal place to study the hadronization process which is described by fragmentation functions. In contrast to conventional approaches based on single-inclusive or semi-inclusive annihilation, the use of fully inclusive jets to study fragmentation functions is more advantageous. This approach provides novel insight into the underlying hadronization mechanisms. For example, the mass of the jet, which breaks chiral symmetry, can be effectively used to probe twist-3 chiral-odd effects. In this paper, we present a systematic theoretical framework for calculating the differential cross-section up to twist-3 for the unpolarized annihilation process. Calculations take into account both the electromagnetic interaction and weak interaction. We find that leading-twist fragmentation functions are convoluted into the twist-3 effects. This probably provides a new perspective for studying fragmentation functions. To this end, we introduce two twist-3 azimuthal asymmetries and give the corresponding numerical estimates. We notice that the twist-3 effects are not significantly suppressed in high-energy regions due to the presence of the weak interaction. In other words, it is feasible to measure twist-3 effects at future high-energy electron-positron colliders to study hadronization mechanisms.
The development of next-generation, high-luminosity, and high-precision charged lepton beam facilities at JLab and MAMI has opened, in recent years, a new frontier in the exploration of weak interaction processes induced by electrons and positrons in the neutral current sector, which can also be used to study weak interaction processes induced by charged currents. In particular, these processes in the intermediate energy regime, spanning from a few hundred MeV to a few GeV, play a crucial role in understanding electroweak dynamics, nucleon structure, and hadronic response functions. This review presents a comprehensive theoretical study of weak charged-current interactions of electrons and positrons with free protons, encompassing quasielastic scattering in both the strangeness conserving and strangeness changing channels, together with inelastic production of the $P_{33}$(1232), $P_{11}$(1440), $S_{11}$(1535) resonances, $η$ and $K$ mesons, and associated production of strange particles. We analyse differential and total cross sections, polarization observables of the final baryons, and spin asymmetries of the proton target, demonstrating their sensitivity to the underlying weak interaction dynamics and to possible second class currents, thereby enabling stringent tests of G- and T- invariance. The explored kinematic region also offers a unique and independent opportunity to constrain the axial vector sector of the weak interaction, and it provides a discussion of alternative ways to determine the axial dipole mass in the quasielastic scattering region, a fundamental parameter that is in debate for nearly two decades. It also focuses on the study of the axial-vector form factors associated with the excitation of the $P_{33}(1232)$ resonance in a manner that is free from the uncertainties inherent in their determination from studies of (anti)neutrino-induced weak processes.
In this talk, I summarize the highlights of our recent and ongoing work on exotic hadrons involving explicitly exotic heavy-flavor mesons. By heavy flavors, I refer in particular to strange and charm quarks. As discussed in this manuscript, several puzzle-like features appearing in recent experimental data can be understood in terms of coupled-channel interactions. The main result is the prediction of an exotic meson carrying both charm and strangeness quantum numbers, whose properties cannot be explained as those of a simple quark-antiquark bound state. I also discuss the future prospects of our research program on this topic.
Using a sequence of effective field theories (EFTs), we calculate the rate for the leptonic decay $B^-\toτ^-\barν_τ(γ)$ including real and virtual QED corrections, with a cut $E_{\rm cut}\llΛ_{\rm QCD}$ imposed on electromagnetic radiation in the $B$-meson rest frame. We establish a factorization theorem for the rate and evaluate it at $\mathcal{O}(α)$, resumming the leading logarithmic corrections to all orders in perturbation theory. The large mass $m_τ$ allows us to treat the tau lepton as a heavy fermion below $μ\sim m_B\sim m_τ$, leading to an EFT construction that is structurally different and noticeably simpler than that for the muon case. In particular, hadron-structure dependent QED corrections can be described in terms of QED-induced $B^-\toτ^-$ form factors, which should be calculable on the lattice. Our analysis includes the leading $Λ_{\rm QCD}/m_τ$ corrections as well as the leading corrections in $E_{\rm cut}/Λ_{\rm QCD}$. We show that contributions of the form $Λ_{\rm QCD}m_B/m_τ^2$, which are naively sizable, cancel among each other. Contrary to the muon channel, structure-dependent corrections involving $BB^\astγ$ transitions are numerically negligible for the tau case. The remaining logarithmic dependence on $E_\mathrm{cut}$ is mild. We estimate that the present uncertainty in the calculation of the structure-dependent QED corrections is about 0.5\% of the rate. As a byproduct, we present a state-of-the-art prediction for the lepton flavor universality ratio of the tau and muon channels.
We develop a semi-analytical framework to map ultra-high-energy neutrino fluxes onto event rates at neutrino telescopes. The formulation is based on the Boltzmann equation for the distribution of secondary muons produced by neutrino interactions in matter, and naturally accounts for the effective target volume relevant for through-going tracks. This "soft volume" is controlled by muon propagation and stochastic energy losses, and can be significantly larger than the instrumented detector volume. Exploiting the dominance of soft energy losses, we derive a controlled second-order expansion of the collision operator, reducing the transport problem to a drift-diffusion equation in energy space, with rare hard scatterings treated perturbatively. The resulting master formula provides a fast alternative to full Monte Carlo simulations while retaining a direct connection to the microscopic muon energy-loss processes. We apply the formalism to IceCube through-going muon data, marginalizing over theoretical uncertainties in the transport coefficients, and obtain a diffuse-flux fit compatible with the experimental result. We also revisit the interpretation of the ultra-high-energy track event reported by KM3NeT, assessing its consistency with IceCube non-observation in the same energy range. Our results provide a first-principles bridge between neutrino-flux models and track-event observables, with direct applications to precision neutrino astronomy and searches for physics beyond the Standard Model.
We investigate the contributions of Lepton Number Violating (LNV) effective operators to the rare decays $B\to K^{(\ast)}νν$ and $K\to πνν$. Such operators can modify the kinematic distributions of these processes, providing distinctive probes of physics beyond the Standard Model. Through a renormalization-group analysis, we show that the Standard Model Effective Field Theory (SMEFT) operators responsible for these effects are subject to stringent indirect constraints from neutrino physics. In particular, we find that the mild excess reported by Belle-II in the $B^+\to K^+νν$ channel cannot be explained by LNV SMEFT operators without introducing significant fine-tuning in neutrino masses. We then show that these constraints can be evaded in the SMEFT minimally extended by a light right-handed neutrino, allowing for sizable effects in rare meson decays. Finally, we explore the implications of this viable scenario for low-energy processes, including neutrinoless double-beta decays.
Self-interacting dark matter turns the structure of galactic halos into a direct requirement on a low-energy scattering amplitude. We show that, for weakly coupled scalar dark matter, this requirement implies a much stronger mass bound on the dark matter particle than partial-wave unitarity alone. Using analyticity, crossing symmetry, locality and partial-wave unitarity, we compute the maximal allowed threshold amplitude with a dispersive primal S-matrix bootstrap, assuming only a weakly coupled EFT below a scale $Λ$ and allowing arbitrary UV particle content above $Λ$. For the benchmark self-interaction cross section $σ_{\rm self}=10^{-24}(M/\mathrm{GeV})\mathrm{cm}^2$, the mass of a generic weakly coupled scalar satisfies $M\lesssim 0.3\,\mathrm{GeV}$ in the controlled EFT regime. If dark matter is a derivative-dominated pseudo-Nambu-Goldstone boson, the mass bound is lowered to the MeV scale or below, depending on the hierarchy $M/Λ$.
CKM unitarity triangle constraints are insensitive to New Physics scenarios in which the $b$-column is uniformly rescaled. This is because B physics data over-constrain the angles of the (bd) unitarity triangle, but only calibrate its normalization. We identify this flat direction and derive model-independent constraints on the rescaling. We find that $|\varepsilon_K|$ already constrains the rescaling at the same level as the direct measurement of the $b$-column normalization, with current data allowing for ${\cal O}(4\%)$ deviations. Projected inputs promote the kaon sector to be the leading probe. A toy model employing mixing with a vector-like up-type singlet provides a proof of principle for a possible UV realization.
Dynamical dark energy offers an alternative to a cosmological constant with distinct observational signatures. However, the small energy density scale, Hubble-sized mass, and Planckian excursions make simple models fine-tuned and unnatural. In this work, we show that a weak version of the axion, identified with the phase field of the anomalous $U(1)_{B+L}$ of the Standard Model, can generate the scale hierarchies expected for dark energy. The axion potential is controlled by sources of explicit baryon and lepton number violation and is radiatively stable. We show that the leading contribution comes from two inequivalent Weinberg operators, one $B+L$-conserving and one $B+L$-violating, which generate the axion potential. We propose a flavor selection rule based on a spontaneously broken $S_3$ permutation symmetry in the lepton sector that simultaneously removes the quadratic divergence and dominant temperature-dependent contributions. The resulting potential first appears at quartic order in the axion-dependent neutrino masses and, for the observed departure from tribimaximal mixing, its amplitude is parametrically close to the dark-energy density. The dominant uncertainty comes from $δ_{\rm CP}$ and $θ_{23}$, so experiments like Hyper-K and DUNE can directly test the model in the future. Cosmologically, the field behaves as thawing quintessence for $f$ close to $M_{pl}$, stays frozen by Hubble friction until late times and never enters the adiabatic regime.
The domain wall is a kind of topological defect that can appear when a discrete symmetry is broken. If the discrete symmetry appears as an intermediate symmetry during a $U(1)$ symmetry breaking, the domain walls are connected to cosmic strings, forming walls bounded by strings. Intuitively, the domain wall disappears if the breaking scale of the discrete symmetry is comparable to that of the $U(1)$ symmetry. In this paper, relying on a 3:1 $U(1)$ model, we show the detailed processes of the disappearance of the domain wall. Due to the existence of the non-negligible ``bias angle'' $β$, the relevance of the ``$Z_3$ symmetry'' and the domain wall is blurred, and thereby the evaluations of the string profiles in a hybrid wall-string network should be revised. We also made some preliminary calculations of the gravitational waves generated by the wall-string network created in the early universe.
The Standard Model prediction for the rare $K_L\toμ^+μ^-$ decay depends critically on the long-distance contribution coming from the exchange of two photons. Such a contribution can be computed using lattice QCD and an effective three-flavor theory including only the $u$, $d$ and $s$ quarks, provided terms falling as the inverse square of the omitted charm quark mass, $1/m_c^2$ are neglected. Because of the missing Glashow-Iliopoulos-Maiani cancelation, this three-flavor theory contains additional low-energy constants that depend on $m_c$. Here we show how these constants can be determined from a practical four-flavor lattice QCD calculation performed on a small volume with $u$ and $d$ quark masses that are heavier than physical.
We study the spatial diffusion coefficient and the momentum transport coefficients of D mesons through a dense nuclear medium in the relaxation time approximation of the kinetic theory. The in medium modifications of the D meson transport properties are computed in the chiral SU(3) hadronic model. Relaxation time is estimated using dilute and degenerate gas approximations for low and high baryonic densities, respectively. We have noticed that relaxation time and spatial diffusion of D meson decrease rapidly in the low density dilute gas domain and mildly in the high density degenerate gas domain. The detailed result of the present work on D meson diffusion is quite contemporary and important towards the compressed baryonic matter densities which can be assessed in future heavy ion collision experiments.
We investigate the origin of the ultra-high-energy $γ$-ray emission detected by the Large High Altitude Air Shower Observatory (LHAASO) from the source LHAASO J0341+5258, which has not yet been associated with any known astrophysical object within the detector's field of view. The observed UHE emission is modelled within two independent frameworks: initially with a time-dependent, source-independent hadronic scenario implemented with the numerical package GAMERA, where particles propagate through the interstellar medium and subsequently interact with the molecular gas observed in the region and finally with an analytical description of the interaction between the accelerated cosmic-ray population from a supernova remnant (SNR) and the molecular gas around it. The relevant parameter space for the hypothetical past SNR is explored using the observed TeV $γ$-rays. We show that, for physically plausible source-cloud separations and propagation timescales, GAMERA provides a first-order approximation to the spectral modifications induced by particle transport and yields an adequate fit to the ultra-high-energy $γ$-ray data. Within the framework of our analytical approach, we demonstrated that the observed GeV emission can plausibly originate from the SNR itself, while the TeV emission detected by LHAASO can be consistently interpreted as arising from particles that have escaped from the SNR and are illuminating nearby molecular clouds. We invoke a spatio-temporally evolved SNR-molecular cloud interaction scenario to account self-consistently for the entire $γ$-ray spectrum from GeV to TeV energies. Despite the remaining uncertainty regarding the nature of the acceleration source, we conclude that the TeV emission detected by LHAASO can be consistently interpreted within the framework of an illumination scenario.
Cosmologically coupled black holes (CCBHs), whose masses evolve in response to the cosmological expansion, have recently attracted significant theoretical and observational interest. Existing studies have treated CCBHs as purely classical objects, neglecting the effect of Hawking radiation (HR), which competes with the cosmological coupling (CC) mechanism. We take a first step towards studying evaporating CCBHs, adopting a quasi-adiabatic approximation in which the HR rate is evaluated at the instantaneous CCBH mass, and modeling the CC mechanism through the phenomenological scaling of the mass with the scale factor, $M \propto a^k$. We show that, depending on the coupling strength $k$, even late-time CC activation can significantly delay Hawking evaporation, or lead to asymptotic CC-dominated mass growth, with important implications. We set limits on the abundance of primordial CCBHs from $γ$-ray observations, finding limits which are weaker than their uncoupled counterparts, as CCBHs are kept farther from the endpoint of evaporation for a longer time. Unlike standard primordial black holes, the CCBH formation and present-day masses no longer approximately coincide, even for formation masses $M_{\text{form}} \gtrsim 10^{15}\,{\text{g}}$. Therefore, the same population of primordial CCBHs may be subject to evaporation limits through its past emission history, as well as to other limits (such as microlensing) through its present-day mass.
Axion-Like-Particles (ALPs) are pseudo-scalar particles actively searched as light dark matter candidates. ALPs can couple to photons which give rise to the possibility of oscillations with photons in an external magnetic field. If sufficiently strong, this coupling can imprint distinctive spectral irregularities in the gamma ray spectrum of astrophysical sources. We present a prospective study on the sensitivity of probing ALP-photon interactions using stacked observations of selected active galactic nuclei (AGNs) located behind galaxy clusters. The ALP-photon conversion in cluster magnetic fields produces absorption-like features in AGN spectra that are difficult to predict for individual sources. To address this, we apply a stacking analysis of multiple AGN-cluster pairs, yielding a controlled prediction of the expected ALP induced spectral patterns and enhancing the sensitivity to such irregularities. Using simulated data for selected hard-spectrum Fermi/LAT AGNs that can be observed by Imaging Atmospheric Cherenkov Telescopes such as H.E.S.S., we evaluate the performance of this method. The combination of mock IACT observations with our stacking approach enable exploration of the previously uncharted ALP dark matter parameter space in the neV mass range.
Space-based gravitational-wave (GW) detectors, including LISA and Taiji, offer unprecedented access to regimes where alternative theories of gravity may deviate from General Relativity (GR). Gravitational Quantum Field Theory (GQFT) provides a novel framework in which the Poincaré-type inhomogeneous spin symmetry of Weyl-type fermions in the Standard Model is elevated to a gauge symmetry. Within this construction, the fundamental gravitational field is identified with a gravigauge field which behaves as a Goldstone-type bi-covariant vector field. Unlike GR, GQFT predicts additional polarization states: one transverse scalar (breathing) mode and two vector modes. In this work, we focus on the transverse, isotropic scalar mode and investigate its detectability with Taiji. To isolate this mode, we employ the null-response channel (NRC), a specific interferometric combination designed to suppress contributions from other polarizations. We implement an analytical, dynamic orbital model to realistically simulate a triangular constellation. We compute the response functions and sensitivity curves for various interferometric channels, compare them with the standard Michelson channel, and demonstrate the effectiveness of the NRC approach. Our results show that the NRC provides a reliable, waveform-independent criterion for testing non-GR polarizations, and we anticipate that it will serve as a valuable tool for probing gravitational theories in future space-based GW missions.
Cosmology stands at a crossroads. The Hubble tension has reached a nominal significance above $7σ$, while analyses combining DESI BAO and Type Ia supernova data show emerging hints of departures from $Λ$CDM. Meanwhile, high-precision CMB measurements from ACT and SPT enable a timely and more stringent reassessment of proposed solutions to the tension. In this paper, we revisit the $H_0$ Olympics, a systematic contest comparing proposed alternatives to $Λ$CDM using common datasets, likelihoods, and statistical criteria. In this updated edition, the $H_0$ World Cup, we subject fourteen representative solutions to a common analysis of current CMB, BAO, and SN data. The contenders span four broad mechanisms: late-time modifications of the expansion history, modified recombination, additional pre-recombination radiation, and early non-radiative energy injection. Relative to the original analysis, the present competition includes models and mechanisms proposed in the intervening years and evaluates all contenders using both Bayesian and Frequentist tests of tension and model performance, letting the neutrino mass sum vary. We further test if late-time extensions through curvature or the Chevallier-Polarski-Linder (CPL) dark energy parametrization can aid the success of the models. Finally, we subject the leading contenders to dedicated robustness tests involving alternative CMB likelihoods and multipole cuts, supernova samples, large-scale-structure information, and big-bang nucleosynthesis constraints. This framework assesses both the ability of each mechanism to ease the Hubble tension and the robustness of our conclusions to datasets and analysis choices.
The Hubble tension has reached a nominal significance above $7σ$, while new high-precision measurements of the cosmic microwave background (CMB) and baryon acoustic oscillations (BAO) sharpen the test of proposed solutions. Using a common framework, we compare fourteen representative alternatives to the standard $Λ$ Cold Dark Matter ($Λ$CDM) model in light of up-to-date CMB, BAO and supernovae data to gauge their ability to resolve the tension. The models span late-time modifications, modified recombination, and exotic pre-recombination expansion histories driven by additional radiation or a localized dark energy injection. We evaluate each proposal with complementary frequentist and Bayesian measures of the residual calibration tension and of the improvement in the joint fit. Both approaches identify the same broad hierarchy. Early dark energy and early modified gravity models perform best, shifting the $H_0$ inference without local measurement priors toward $70\,\mathrm{km\,s^{-1}\,Mpc^{-1}}$ and reducing the residual discrepancy with SH0ES to approximately $2.5-3.6σ$, depending on the model and statistic, while receiving strong support over $Λ$CDM in the combined fit. Varying the electron mass at recombination yields an intermediate improvement, whereas the enhanced-radiation and late-time scenarios do not improve over $Λ$CDM. This Letter summarizes the group stage of the competition; in a companion paper (Paper II) we present the results of an exhaustive set of analyses and assess their robustness to variations in modeling assumptions and datasets.
With current advances in gravitational wave (GW) detection made by the worldwide LIGO-Virgo-KAGRA (LVK) network of detectors, ever-more sensitive tests of gravity in the strong-field regime are now possible. This enables one to test gauge theories beyond Einstein-Hilbert action, such as Weyl gauge theories of gravity. The only anomaly-free (quantum) gauge theory of a space-time symmetry beyond Poincaré is based on Weyl gauge group (of dilatations and Poincaré symmetry) with Weyl conformal geometry as its natural underlying geometry. This gauge theory has spontaneous breaking of Weyl gauge symmetry to Einstein-Hilbert and Proca actions, plus a positive cosmological constant. We investigate the GW polarisation modes of Weyl (quadratic) gauge theory of gravity in Weyl geometry and compare our findings to the most recent experimental data. We show how the geodesic deviation equation from Riemannian geometry translates to Weyl geometry, and explain why it is crucial to perform the analysis around de Sitter background, which is the correct low-energy limit of Weyl quadratic gravity, to not alter the GW content, and then compute the polarisation modes. In addition to the two transverse-traceless tensor modes predicted by Einstein-Hilbert action, we find two additional vector modes induced by the transverse fluctuations of the Weyl gauge field. If detected, these vectors modes would be important evidence for Weyl gauge symmetry.