Axions from the anomalous global Peccei–Quinn symmetry resolve the strong CP problem but are susceptible to Planck-scale suppressed operators. A gauged abelian flavor symmetry U(1)_F, introduced to account for mass hierarchies via the Froggatt–Nielsen mechanism, shields the axion, yielding an accidentally high-quality flavored axion with a domain wall number of one. The model predicts two complementary signals: flavor-changing currents in K→π a decays (typical for flavor scales Λ_FN ≳ f_a) and stochastic gravitational waves from the evolution and decay of networks of gauged flavor and axion cosmic strings. Global axion strings efficiently radiate axions, potentially explaining the observed dark matter density. The gravitational-wave spectrum acquires a “plateau–valley” structure—a powerful and characteristic probe for high-quality flavored axion dark matter models, supplementing low-energy flavor experiments.
The Standard Model of particle physics, despite its successes, leaves two fundamental puzzles unanswered. The strong CP problem—the negligibly small CP violation parameter in quantum chromodynamics—requires fine-tuning. The Peccei–Quinn mechanism turns this parameter into a dynamical field, the axion, but the axion's quality suffers from Planck-scale corrections. Simultaneously, the hierarchy of quark and lepton masses and mixings—the flavor puzzle—hints at a hidden symmetry. Vera Rubin convincingly showed that dark matter dominates galaxies, and the axion remains one of the best candidates. This work weaves these threads together: a U(1)F gauge flavor symmetry protects the axion from quantum gravity effects, and the cosmic strings born when it breaks become sources of gravitational waves with a unique spectral structure.
The authors construct a model extending the standard DFSZ axion scheme with a flavon field X and a U(1)F gauge symmetry. The fermion mass hierarchy is generated via the Froggatt–Nielsen mechanism: Yukawa couplings are suppressed by powers of a small parameter ε~0.22. Spontaneous breaking of U(1)F forms a network of two types of cosmic strings—'axionic' (S field) and 'flavonic' (X field). Their evolution in the expanding Universe is traced from Planckian scales down to the QCD epoch. Gravitational-wave emission is computed by summing contributions from oscillating string loops; for axionic strings, axion emission dominates, but after the QCD phase transition, flavonic strings reorganize into purely gauge strings and start radiating gravitational waves efficiently. Key parameters: the flavor scale ΛFN~10^12 GeV, the vacuum expectation value ratio r=vS/vX, and the axion decay constant fa.
The gravitational-wave spectrum exhibits a characteristic 'plateau–dip' in the millihertz frequency range. The plateau arises from radiation by purely gauge string loops, while the dip results from the abrupt transition from global to gauge strings at the QCD epoch. The cutoff frequency fcut≈10^(−3) Hz is directly related to the flavor scale ΛFN. For axion dark matter with fa~1.26×10^10−10^11 GeV (the relic density window), the amplitude ΩGWh^2 reaches 10^(−11)−10^(−9), falling within the sensitivity band of future detectors LISA, BBO, DECIGO. For certain r, the axion becomes neutronphobic (Can=0) or electrophobic (Cae=0), relaxing astrophysical bounds and narrowing the parameter space. Combined analysis with NA62 data on rare K-meson decays highlights a preferred band ΛFN~10^12 GeV, where solutions to all five puzzles intersect.
The proposed mechanism links for the first time five fundamental problems into one coherent picture: neutrino masses, baryon asymmetry (leptogenesis), the flavor hierarchy, the strong CP problem, and dark matter. The characteristic plateau–dip structure in the gravitational-wave spectrum becomes a smoking-gun signature of a high-quality flavored axion. Detecting such a signal would not only confirm the gauge origin of flavors but also rule out many alternative axion models, providing a crucial test for quantum gravity.
The subject will advance along two fronts. On one hand, precision measurement of the frequency and shape of the gravitational-wave dip with LISA, DECIGO, and ground-based facilities (ET, CE) will allow extraction of the flavor scale and ratio r, turning GW astronomy into a spectroscope of high-energy physics. On the other, improved limits on rare K→πa decays from HIKE experiments will narrow the allowed ΛFN range, while direct axion detection (e.g., ADMX) could complete the chain of evidence.
The results will deeply impact gravitational-wave astronomy, particle physics, and cosmology, as they tie the origin of dark matter to the dynamics of the early Universe. The methods for analyzing string networks will be applicable to a broad class of models beyond the Standard Model.
Numerical simulations of the reorganization of string-wall networks at the QCD epoch taking axion emission into account are needed. In parallel, targeted searches for a correlated signal in NA62 data and preparations for HIKE are necessary.
The model simultaneously addresses: the strong CP problem (axion), the flavor hierarchy (FN mechanism), dark matter (axion), neutrino masses (type I seesaw), and baryon asymmetry (resonant leptogenesis). Cosmic strings and gravitational waves become the experimental bridge between these challenges.
🎯 Within the model, an axionic string can decay into several unit strings that annihilate, avoiding catastrophic domain wall domination. And the characteristic dip in the gravitational-wave spectrum serves as a 'cosmic thermometer,' pinning down the temperature of the QCD phase transition.