paper Review Profile
Log-Periodic Signatures from Discrete Scale Invariance in the Stochastic Gravitational-Wave Background: Walking Technicolor as a Candidate Ultraviolet Completion
The paper shows that discrete scale invariance (DSI) in the anisotropic stress of a first-order cosmological phase transition produces a multiplicative log-periodic modulation of the stochastic gravitational-wave background, and proves a factorization theorem that this modulation survives the unequal-time integrals to percent-level accuracy under the short-correlation-time approximation. As an explicit ultraviolet completion, it demonstrates that walking technicolor can realize the required DSI and predicts a falsifiable parameter band (epsilon ~ 0.04–0.18, b ~ 1.7–2.8) that lies in the…
Full breakdown: https://theoryofeverything.ai/papers/log-periodic-signatures-from-discrete-scale-invariance-in-the-stochastic-gravitational-wave-background-walking-technicolor-as-a-candidate-ultraviolet-completion
This paper presents an interesting phenomenological proposal that discrete scale invariance (DSI) in the anisotropic stress tensor during first-order phase transitions can imprint log-periodic modulation on the stochastic gravitational-wave background. The core theoretical framework is mathematically sound: if the source UETC has a multiplicative k-only DSI modulation as assumed in Eq. (6), then the factorization theorem rigorously shows this modulation survives the unequal-time integration under short-correlation-time conditions. However, significant consistency issues emerge in connecting this phenomenological backbone to the proposed walking technicolor UV completion. The math specialists identified critical gaps where the WTC construction produces modulation through internal-momentum convolution (Eqs. 33-38), but the paper does not demonstrate this yields the same mathematical form as the external multiplicative factor C(k) assumed in the phenomenology. Additionally, the short-correlation approximation k^{-1} >> τ_corr conflicts with the paper's own parameter estimates k τ_corr ~ 1 at the spectral peak. The abstract claims 'percent-level' accuracy while Table 1 shows ~15% error at the stated threshold β/H_* = 10. Despite these internal inconsistencies, the work provides a concrete falsifiable prediction with well-defined observable signatures and explicit detectability analysis for LISA.
I agree with the higher-score assessments that the paper is unusually explicit about its hierarchy of claims and that, at the purely phenomenological level, the logic is coherent: if Eq. (6) is assumed as a k-only multiplicative modulation of the UETC, then Eq. (4) carries that factor into P_h and Ω_GW by linearity. I also acknowledge the strongest pro-consistency point raised by the 4/5 assessment: the paper reportedly distinguishes derived results, controlled approximations, UV-completion claims, and conjectural motivation, and it distinguishes potential-level parameters (ε_f,b0) from observable parameters (ε,b) at leading order. That transparency mitigates some possible objections. However, it does not remove the central consistency problem. The same modulation object is used in two logically different ways: first as the phenomenological UETC multiplier C(k) in Eq. (6), and later as the endpoint of a WTC chain from a periodic potential to propagators to a convolutional UETC. The paper does not show that the latter has the same mathematical form, amplitude, period, and phase as the former. This matters because the main conclusions include not only the conditional template but also the WTC parameter band. The strongest low-score objection is also persuasive: Eq. (9) invokes k^{-1} >> τcorr, but with the paper’s later characteristic estimates k ~ β/vw and τcorr ~ vw/β, one obtains kτcorr ~ 1, so the stated short-correlation control is not internally aligned with the parameter regime used for predictions. In addition, the claimed percent-level validity for β/H_* ≳ 10 appears inconsistent with the reported Table 1 value of roughly 15% at β/H_*=10, and the Sec. 5.4 ε=ε_f identification after a two-cross-term convolution is not a harmless normalization choice. Because these issues affect the core definition and quantitative use of ε,b,C(k), the central-definition-drift cap applies and the internal-consistency score cannot exceed 2. A consensus round resolved an earlier panel split before this score was finalized.
The phenomenological backbone (Sec. 3.2, factorization theorem) is solid: the short-correlation-time approximation is justified, the Riemann-Lebesgue-style bound on the Green's-function oscillations (Eq. 13) is correctly applied, the factorization of C(k) out of the η integral is algebraically exact (correctly noted in Sec. 5.1), and the numerical validation in Fig. 4 is appropriate as a consistency check on the algebraic claim. Eq. (15) and the observable template (16) follow rigorously from the stated assumptions. The matched-filter SNR derivation (Sec. 4.2) is standard and correct. However, the WTC UV-completion derivation (Sec. 5.3) contains two load-bearing unverified steps: Eq. (31), the holographic derivation of the multiplicative log-periodic 4D potential from a sinusoidal warp-factor perturbation, is asserted with citation but no derivation — the factor of 4, the cos(ln φ) structure, and ln b_0 = kL/n_p require an explicit Goldberger-Wise-style radion calculation that is not shown; and Eq. (32), the propagation of the modulation from φ to q via chain rule, glosses over the Δ_φ scaling-dimension factor that should modify the period. These steps are central to the paper's quantitative WTC prediction band (Eq. 39), which is the headline falsifiable result. The author's Table 2 epistemic hierarchy partially mitigates this by labelling the UV layer 'candidate' and 'conjectural', but the specific numerical band ε ∈ [0.04, 0.18], b ∈ [1.7, 2.8] is presented as a sharp prediction, and its derivation chain has unfilled steps. Score capped at 3 per the unverified_central_derivation rule.
The work is reasonably falsifiable. Its main phenomenological prediction is a specific spectral template: a sinusoidal modulation in ln f multiplying the SGWB baseline, parameterized by ε, b, and φ0. That is a concrete, differentiable signature rather than a generic statement that 'the spectrum is altered.' The paper also supplies a candidate parameter range, ε≈0.04–0.18 and b≈1.7–2.8, and gives an explicit detectability scaling for matched filtering in the LISA band. A null search for persistent log-periodic oscillations over multiple periods in the relevant frequency range would directly constrain or exclude the proposed parameter region, so the proposal is genuinely testable in principle and plausibly in practice. The main limitation is that the falsification criteria are not stated as sharply as they could be. The paper does not specify a detection pipeline threshold, model comparison criterion, or exclusion contour that would count as ruling out the framework, beyond qualitative statements that non-detection would bound ε as a function of b. In addition, much of the UV-completion discussion is only a candidate realization, so falsification of the WTC embedding is less clean than falsification of the phenomenological DSI template itself. Still, the observable template and parameter band are specific enough to merit a strong score.
The paper is generally organized well: the progression from general SGWB formalism to DSI ansatz, observable template, detectability estimate, and candidate UV completion is easy to follow structurally. Many assumptions are explicitly signposted, and Table 2 is particularly helpful in separating derived results from conjectural UV motivation. The author also does a good job of distinguishing the phenomenological backbone from the speculative embedding, which improves interpretability. However, clarity is reduced by a few important issues. First, there is material overclaim in the abstract relative to the body, especially concerning the WTC realization and the level of proof behind the factorization claim. Second, notation/terminology shifts require care: SNR_baseline is used in a specialized per-log-period sense that is easy to misread, and several modulation parameters (C, δ, ε, ε_f, b, b0) are related but not always introduced in the cleanest sequence. Third, some argumentative steps are compressed and would force a graduate-level reader to reread, especially where simulation validation is presented as support for a broader theorem even though the tested setup omits the non-separable corrections later emphasized. Because of the red-flag caps from term redefinition and abstract overclaim, clarity cannot exceed 3.
The paper's central novelty is the proposal that discrete scale invariance in the source anisotropic stress of a first-order phase transition yields a multiplicative log-periodic modulation in the SGWB, together with the claim that this structure survives the unequal-time integration under short-correlation assumptions. That is a nontrivial reinterpretation of SGWB source physics and, if correct, gives a distinctive template not standard in phase-transition GW phenomenology. The attempt to connect this to a candidate UV realization in walking technicolor further broadens the scope beyond a purely abstract signal model. The score is not a 5 because the work is partly a synthesis of known ingredients: DSI/log-periodicity, standard SGWB phase-transition machinery, matched filtering, and an existing WTC phase-transition context. The manuscript itself acknowledges prior discussion of log-periodic GW features in other contexts. The most original contribution is therefore the specific source-to-observable mechanism and the proposed WTC embedding, not an entirely new mathematical structure. Since the UV completion remains conjectural rather than derived, the novelty is substantial but not fully secured at the deepest mechanistic level.
The paper presents a well-developed argument from phenomenological framework to specific predictions. The factorization theorem is rigorously derived with explicit error bounds (Table 1). All key variables are defined, boundary conditions are addressed through the short-correlation-time approximation (τ_corr H_* << 1), and limitations are clearly stated. The WTC embedding includes a complete calculation chain from potential modulation to observable spectrum. Minor gaps include: (1) the holographic origin of the periodic warp factor (Eq. 30) is motivated but not fully derived from first principles, and (2) some intermediate steps in the convolution calculation could be more detailed. However, the core argument is complete and the main stated goals are fully addressed.
Strengths
- +Novel phenomenological mechanism linking discrete scale invariance to stochastic gravitational-wave background through multiplicative log-periodic modulation
- +Rigorous factorization theorem (Eqs. 10-15) showing DSI modulation survives unequal-time integration under stated approximations
- +Concrete falsifiable predictions with sharp parameter band ε ∈ [0.04,0.18], b ∈ [1.7,2.8] and explicit LISA detectability analysis
- +Exceptional epistemic transparency with Table 2 explicitly distinguishing derived results from controlled approximations and conjectural UV completion
- +Comprehensive matched-filter SNR analysis providing practical detectability scaling and enhancement factors
Areas for Improvement
- -Demonstrate equivalence between phenomenological UETC multiplier C(k) and WTC convolution output, or quantify the non-separable remainder
- -Resolve inconsistency between short-correlation condition k^{-1} >> τ_corr and parameter regime k τ_corr ~ 1
- -Reconcile abstract claim of 'percent-level' accuracy with Table 1 showing ~15% error at β/H_* = 10
- -Provide explicit derivation of holographic mapping from periodic warp factor to technidilaton potential (Eq. 31)
- -Justify factor-of-2 absorption and ε = ε_f identification in convolution calculation (Sec. 5.4)
- -Include numerical validation that tests non-separable corrections rather than only separable case
Log-Periodic Signatures from Discrete Scale Invariance in the Stochastic Gravitational-Wave Background Walking Technicolor as a Candidate Ultraviolet Completion Jill F. Rankin Independent Researcher jill.rankin@g.austincc.edu May 2026(preprint) Abstract We show that discrete scale invariance (DSI) in the anisotropic stress tensor during a first-order cosmological phase transition imprints a multiplicative log- periodic modulation on the stochastic gravitational-wave background (SGWB). Under the physically motivated short-correlation-time approximation (τ corr H ∗ ≪1, satisfied forβ/H ∗ ≳10), the DSI modulation factorizes from the source unequal-time correlator to the observable energy-density spectrum at the percent level, yielding Ω GW (f ) = Ω 0 GW (f ) 1 + ε cos 2π ln(f/f ∗ ) lnb
- φ 0 , with modulation amplitudeε≪1 and discrete scaling ratiob >1. Matched-filter de- tectability of the oscillatory component scales asSNR osc ≃(ε/ √ 2)SNR baseline p N periods , whereN periods =ln(f max /f min )/ lnbis the number of complete log-periods in the detector band andSNR baseline is the per-log-period baseline SNR, giving a useful enhancement over the naive ε suppression. As a candidate ultraviolet completion we explore the realization of the required DSI within walking technicolor (WTC), a strongly coupled hidden-sector gauge theory that (i) naturally provides approximate continuous scale invariance broken to DSI by a small periodic modulation of the technidilaton potential, and (ii) produces a strong first-order phase transition already known to generate LISA-detectable gravitational waves. A convolution calculation shows that, under the same short- correlation-time approximation, the DSI propagates from the technidilaton potential to the observable SGWB with errors≲1%. The WTC parameter space predicts ε ∈[0.04,0.18],b ∈[1.7,2.8], which occupies the high-SNR region of the LISA detectability plane, providing a sharp falsifiable target. All approximations are quantitatively bounded.
Contents 1 Introduction3 2 Gravitational-Wave Tensor Power Spectrum4 3 Discrete Scale Invariance in the Source UETC5 3.1 DSI ansatz . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .5 3.2 Factorization theorem . . . . . . . . . . . . . . . . . . . . . . . . . . . . .5 4 Observable Signatures7 4.1 DSI-modulated energy-density spectrum . . . . . . . . . . . . . . . . . . .7 4.2 Matched-filter detectability . . . . . . . . . . . . . . . . . . . . . . . . . . .7 5 Ultraviolet Completion: Walking Technicolor9 5.1 Numerical validation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .9 5.2 Phase-transition parameter space . . . . . . . . . . . . . . . . . . . . . . . 11 5.3 Engineering discrete scale invariance . . . . . . . . . . . . . . . . . . . . . 11 5.4 Convolution for the UETC . . . . . . . . . . . . . . . . . . . . . . . . . . . 14 5.5 WTC predictions . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15 6 Discussion15 7 Conclusions16 2
1 Introduction The stochastic gravitational-wave background (SGWB) from first-order cosmological phase transitions is among the most promising observational targets for current and next-generation gravitational-wave detectors. The Laser Interferometer Space Antenna (LISA) [1] will be sensitive to phase transitions occurring at temperaturesT ∗ ∼10–10 4 GeV , covering a broad class of beyond-Standard-Model (BSM) scenarios. Pulsar-timing arrays (PTAs) have now reported evidence for a gravitational-wave background at nano-hertz frequencies [3–5], with spectra consistent with — though not yet uniquely identified as — a cosmological phase-transition origin. In this environment, spectral features that go beyond the smooth envelope predicted by conventional calculations take on special importance: they carry direct information about the microphysics of the transition and the nature of any BSM sector responsible for it. Standard calculations of the SGWB from a first-order phase transition predict a broad-band spectrum shaped by three source contributions — bubble collisions [8], sound waves [6,7], and magneto-hydrodynamic turbulence [9] — each with a characteristic broken power-law profile. A variety of beyond-standard effects can modify this picture: strong supercooling can sharpen the bubble-collision peak [8]; non-runaway walls alter the sound-wave contribution [6]; and non-equilibrium dynamics can generate additional log contributions [7]. However, none of these mechanisms generically produces a coherent log-periodic oscillation superimposed on the spectrum. Discrete scale invariance (DSI) is the symmetry that does. A system is said to possess DSI with ratiob >1 if it is invariant only under the discrete rescalingx→ b n xfor integer n, rather than under all continuous dilations [10]. DSI arises in hierarchical lattice models, fractal structures, iterated-function-system attractors, and — crucially for our purposes — near-conformal gauge theories with explicit periodic modulations. Its universal observable consequence is a log-periodic correction to any power-law observable, F (x) = x D 0 1 + A cos 2π lnx lnb
- φ ,(1) arising from complex scaling dimensionsD n =D 0 ±2πin/ lnbin the spectrum of the dilatation operator [10]. DSI and its signatures have been studied in condensed-matter physics [10] and in financial time-series analysis [11], but its imprint on the SGWB has received comparatively little attention. Log-periodic features in the SGWB have been discussed in the context of non-standard inflationary scenarios and beyond-Einstein-gravity models [12]. In this paper we pursue a more direct route: we show that DSI in the anisotropic stress tensor of a first-order phase transition itself imprints a multiplicative log-periodic modulation on the observable SGWB. The mechanism operates at the level of the source unequal-time correlator (UETC) and is not specific to any particular BSM sector. The key technical result is a factorization theorem: under the physically well-motivated short- correlation-time approximation, valid for all realistic first-order phase transitions with β/H ∗ ≳10, the DSI modulation passes through the double time-integral of the tensor power spectrum unchanged, at the percent level. As a candidate ultraviolet (UV) completion we explore the realization of the required DSI within walking technicolor (WTC) [13], a strongly coupled hidden-sector gauge theory in the near-conformal regime. Two features motivate WTC as a host for DSI: (i) walking dynamics naturally provide approximate continuous scale invariance over a 3
wide range of energies, which can be broken to DSI by a small periodic modulation of the technidilaton effective potential — motivated (but not derived from first principles) by holographic models with periodic warp factors and by RG-group limit-cycle structure near the quasi-fixed point; and (ii) the WTC phase transition is already known to generate LISA-detectable gravitational waves [13], placing the DSI-modulated prediction squarely in the observable band without requiring any new tuning. We perform an explicit convolution calculation that traces the DSI modulation from the technidilaton potential through the UETC to the observable Ω GW (f ), with every approximation quantified. The resulting prediction is sharp: the WTC parameter space maps onto a specific bandε ∈[0.04,0.18],b ∈[1.7,2.8] in the DSI amplitude–ratio plane, which overlaps the high-SNR region of the LISA detectability forecast. A matched-filter search for the log-periodic template provides an optimal discriminant. A companion paper [18] demonstrates the log-periodic spectral imprinting mechanism in a controlled one-dimensional electromagnetic cavity using finite-difference time-domain (FDTD) simulations, providing a numerical proof of concept independent of gravitational- wave physics. The theoretical framework of dynamic mode-accessibility engineering that unifies both papers is developed in Ref. [19]. The paper is structured as follows. Section 2 reviews the tensor power spectrum and sets up notation. Section 3 states the DSI ansatz and derives the factorization theorem. Section 4 works out the observable signature, quadratic relic corrections, and matched-filter detectability. Section 5 develops the WTC UV completion. Section 6 discusses robustness, distinguishability, and extensions. Section 7 summarizes the main results. Throughout we use natural units c =ℏ = k B = 1 and metric signature (−, +, +, +). 2 Gravitational-Wave Tensor Power Spectrum Tensor metric perturbationsh ij in a flat Friedmann–Lemaˆıtre–Robertson–Walker (FLRW) background satisfy h ′′ ij (k,η) + 2Hh ′ ij (k,η) + k 2 h ij (k,η) = 16πGa 2 (η) Π TT ij (k,η),(2) where primes denote derivatives with respect to conformal timeη,H=a ′ /a,a(η) is the scale factor, and Π TT ij is the transverse-traceless projected anisotropic stress sourced by the phase transition. The two-point function of the source defines the unequal-time correlator,
Π TT ij (k,η) Π TT∗ ij (k ′ ,η ′ ) = (2π) 3 δ (3) (k−k ′ ) Π(k,η,η ′ ),(3) where statistical isotropy has been used to write Π as a function ofk=|k|. Solving Eq. (2) with the retarded Green’s function G k (η,η ′ ) gives the tensor power spectrum, P h (k,η) = (16πG) 2 Z dη 1 dη 2 G k (η,η 1 )G k (η,η 2 )a 2 (η 1 )a 2 (η 2 ) Π(k,η 1 ,η 2 ).(4) The fractional GW energy density per logarithmic frequency interval, referred to the critical density today, is [7] Ω GW (k,η)≃ k 3 12a 2 H 2 P h (k,η),(5) valid for sub-horizon modesk ≫ H. In what follows we work in terms of the observed frequency f = k/(2πa 0 ). 4
3 Discrete Scale Invariance in the Source UETC 3.1 DSI ansatz We assume that the source UETC carries a discrete scale invariance with ratiob >1 and amplitude ε≪ 1: Π(k,η,η ′ ) = Π 0 (k,η,η ′ ) 1 + ε cos 2π ln(k/k ∗ ) lnb
- φ 0 ,(6) where Π 0 is the smooth DSI-free UETC,k ∗ is a reference scale, andφ 0 is an overall phase. Equation (6) is the leading-order expression consistent with invariance underk → b n kfor integern; the log-periodic modulation is the real part of the complex power-law correction associated with complex scaling dimensions [10]. 3.2 Factorization theorem For a first-order phase transition the UETC naturally separates into macroscopic (slow) and microscopic (fast) parts, Π(k,η,η ′ ) = S(η,η ′ )F (k,η− η ′ ),(7) whereS(η,η ′ ) describes the macroscopic source evolution (slowly varying on the Hubble timeH −1 ∗ ) andF(k,∆η) encodes temporal correlations (decaying onτ corr ≪ H −1 ∗ ). This separation holds when the source is stationary on timescalesτ corr ≪∆η ≪ H −1 ∗ : a good approximation for the envelope and sound-shell contributions [6,8], for whichSis approximately constant while F decays rapidly. The phase-transition source decorrelates on the bubble radius/wall-speed timescale τ corr ∼ R ∗ ∼ v w /β, giving τ corr H ∗ ∼ v w β/H ∗ ≪ 1for β/H ∗ ≳ 10.(8) In this limitF(k,∆η) is sharply peaked at ∆η= 0. To bound the error, expandFabout ∆η= 0:F(k,∆η) =F(k,0)δ τ corr (∆η) +O(τ corr H ∗ ), whereδ τ corr is a nascent delta function of widthτ corr . Substituting into the double (η 1 ,η 2 ) integral of Eq. (4), theη 2 integral is dominated by the region|η 2 − η 1 |≲ τ corr . SinceG k (η,η 2 )a 2 (η 2 ) varies on timescales k −1 ≫ τ corr (for sub-horizon modes) andH −1 ∗ ≫ τ corr , it may be evaluated atη 2 =η 1 , giving F (k,η− η ′ )≃ F (k)δ(η− η ′ ) +O(τ corr H ∗ ),(9) whereF(k)≡ R F(k,∆η)d(∆η). The relative error isO(τ corr H ∗ ) =O(v w /(β/H ∗ )), bounded explicitly in Table 1. We decompose the spectral kernel asF(k) =C(k)F 0 (k), whereF 0 (k) is the smooth baseline kernel andC(k) = 1 +ε cos[2π ln(k/k ∗ )/ lnb+φ 0 ] carries the DSI modulation. Substituting into Eq. (4) and performing the η 2 integral using the delta function: P h (k,η)≃ (16πG) 2 C(k) Z dη 1 G 2 k (η,η 1 )a 4 (η 1 )F 0 (k)S(η 1 ,η 1 ).(10) C (k) depends only onk, not onη 1 , and therefore factors out of theη 1 integral exactly. The remaining integral, together with the prefactors, defines the smooth tensor power spectrum: P 0 h (k,η)≡ (16πG) 2 F 0 (k) Z dη 1 G 2 k (η,η 1 )a 4 (η 1 )S(η 1 ,η 1 ).(11) 5
We now show explicitly thatP 0 h acquires no log-periodic structure from the Green’s function. ExpandingG 2 k (η,η 1 ) =sin 2 [k(η − η 1 )]/k 2 = [1− cos(2k(η − η 1 ))]/(2k 2 ), the integral splits as Z dη 1 G 2 k a 4 S = I 0 (η) 2k 2 − e I(k,η) 2k 2 ,(12) whereI 0 (η) = R dη 1 a 4 (η 1 )S(η 1 ,η 1 ) is strictlyk-independent, and e I(k,η) = R dη 1 cos[2k(η− η 1 )]a 4 Sis an oscillatory Fourier transform of the slowly varying envelope. Integrating e I by parts once (boundary terms vanish since S = 0 outside the source epoch): e I(k,η) ≤ 1 2k Z ∂ η 1 [a 4 S] dη 1 ≲ H ∗ k I 0 (η),(13) sincea 4 Svaries on the Hubble timescale:|∂ η 1 [a 4 S]|≲ H ∗ a 4 S. For sub-horizon GW modes k ≫ H ∗ , the ratio | e I|/I 0 ≲ H ∗ /k. To boundH ∗ /kin terms of the factorization error, note that the characteristic GW wavenumber isk ∼ βa ∗ /(v w a 0 ), givingH ∗ /k ∼ H ∗ v w a 0 /(βa ∗ ) =v w /(β/H ∗ )≡ τ corr H ∗ . Thus | e I|/I 0 =O(τ corr H ∗ ) and Z dη 1 G 2 k a 4 S = I 0 (η) 2k 2 1 +O(τ corr H ∗ ) .(14) Thek-oscillations of e Ioccur on the linear scale ∆k ∼1/η ∗ (producing the sound-wave broken power-law spectral features [6]); in log-kspace this corresponds to ∆lnk ∼ 1/(kη ∗ )∼ τ corr H ∗ ≪ lnb. These sound-wave features are spectrally separated from the DSI modulation by the large factorlnb/(τ corr H ∗ )∼(β/H ∗ /v w )lnb≫1: there is no overlap in log-frequency space, andP 0 h acquires no log-periodic structure at periodlnb. Hence, combining Eq. (14) with the exact factorization of C(k): P h (k,η) = C(k)P 0 h (k,η) 1 +O(τ corr H ∗ ) ,(15) whereP 0 h is the tensor power spectrum evaluated with the smooth UETC Π 0 . This is the factorization theorem: the DSI modulation transfers multiplicatively from the source to the tensor power spectrum, with the Green’s-function contribution bounded explicitly by Eq. (13). Higher-order corrections are quantified in Table 1. Table 1: Relative error bound on the factorizationP h =C(k)P 0 h , forε= 0.1,v w = 1. Two independent corrections contribute. (i) Delta-function approximation, Eq. (9): relative error ≤ τ corr H ∗ ≡ v w /(β/H ∗ ), obtained by bounding the correction R dη 1 |G k (η,η 1 )|·O(τ corr )· |G k (η 2 ,η 1 )a 2 |relative to the leading term. (ii) Convolution factorization, Sec. 5.4: relative error≤(2π/ lnb)ε f τ corr H ∗ ≲13ε f τ corr H ∗ forb ≥1.5, using|dδ/d lnq| ≤ ε f ×2π/ lnb (the absolute derivative, uniformly bounded even at zeros ofδ). (iii) Green’s function oscillations, Eq. (13): relative error≤ H ∗ /k ≡ τ corr H ∗ from the Riemann–Lebesgue bound. The combined bound at ε = 0.1 is (1 + 2π/ lnb + 1)τ corr H ∗ ≈ 15τ corr H ∗ for b = 2. β/H ∗ τ corr H ∗ Combined bound Dominant term 100.10≲ 15%δ-fn 1000.01≲ 1.5%convolution 10000.001≲ 0.15%Green’s fn 6
Numerical validation. A numerical validation of the factorization theorem is presented in Sec. 5.1 (Fig. 4). 4 Observable Signatures 4.1 DSI-modulated energy-density spectrum Combining Eq. (6) with the factorization (15) and using Eq. (5), the observable GW energy-density spectrum is Ω GW (f ) = Ω 0 GW (f ) 1 + ε cos 2π ln(f/f ∗ ) lnb
- φ 0 .(16) The fractional residual R(f )≡ Ω GW (f )− Ω 0 GW (f ) Ω 0 GW (f ) = ε cos 2π ln(f/f ∗ ) lnb
- φ 0 (17) is exactly a sinusoid inlnfwith period ∆lnf=lnb, amplitudeε, and phaseφ 0 . By construction, its Pearson correlation coefficient with a fixed-period cosine template equals r= 1.00 for any bandwidth spanning complete log-periods. This analyticr= 1 should be distinguished from the empiricalr= 0.81±0.04 reported in the companion FDTD paper [18]: in that setting, finite time-series length, Hann-window spectral leakage, and imperfect power-law envelope subtraction all reduce the observed correlation below unity. The FDTD value quantifies detection efficiency in a realistic finite-bandwidth experiment; r= 1 is the property of the underlying physics, recovered in the limit of infinite bandwidth and exact envelope knowledge. For the smooth baseline Ω 0 GW we adopt the standard sound-wave contribution [6, 7], Ω sw (f )h 2 = 2.65× 10 −6 H ∗ β 2 κ sw α 1 + α 2 100 g ∗ 1/3 v w S sw (f ),(18) S sw (f ) = f f sw 3 7 4 + 3(f/f sw ) 2 7/2 ,(19) with peak frequency f sw = 1.9× 10 −5 Hz 1 v w β H ∗ T ∗ 100 GeV g ∗ 100 1/6 .(20) Hereαis the transition strength,κ sw is the fraction of the released latent heat converted to fluid bulk motion,g ∗ is the number of relativistic degrees of freedom atT ∗ , andv w is the wall velocity. We set the DSI reference scalef ∗ ∼ f sw . Figure 1 shows the spectrum, residual, and log-period spacing for representative parameter values. 4.2 Matched-filter detectability The oscillatory component of the signal is δΩ GW (f ) = ε Ω 0 GW (f ) cos 2π ln(f/f ∗ ) lnb
- φ 0 .(21) 7
10 3 10 2 10 1 10 0 10 1 Frequency f [Hz] 10 13 10 12 10 11 10 10 10 9 10 8 10 7 h 2 GW ( f ) (a) = 0.10,b = 2.0, 0 = 0 Smooth baseline 0 GW DSI-modulated GW 10 3 10 2 10 1 10 0 10 1 Frequency f [Hz] 0 + R ( f ) / 0 (b) f * Analytic: r = 1.00 (see text; cf. FDTD: r = 0.81 ± 0.04) ( GW 0 )/ 0 Fixed-period fit: cos(2ln(f/f * )/ln b) 10 3 10 2 10 1 10 0 10 1 Frequency f [Hz] ln f = ln b f * (c) log-period spacing Figure 1: Log-periodic modulation of the SGWB. (a) Power spectrumh 2 Ω GW (f) (orange, solid) and smooth baselineh 2 Ω 0 GW (f) (blue, dashed) versus frequency, forε= 0.1,b= 2, φ 0 = 0. (b) Fractional residualR(f)≡[Ω GW (f)−Ω 0 GW (f)]/Ω 0 GW (f), showing the clean sinusoidal oscillation inlnfpredicted by Eq. (16). The orange curve is the fixed-period cosine fit. (c) Log-period spacing: vertical ticks mark frequencies where the modulation peaks (cos = +1), equally spaced by ∆ lnf = lnb. The reference scale f ∗ is indicated. 8
The squared matched-filter signal-to-noise ratio for a search with fixed template parameters (b,φ 0 ) is SNR 2 osc
Z [δΩ GW (f )] 2 σ 2 (f ) d lnf = ε 2 Z [Ω 0 GW (f )] 2 σ 2 (f ) cos 2 2π ln(f/f ∗ ) lnb
- φ 0 d lnf.(22) Over N periods complete log-periods ⟨cos 2 ⟩ = 1/2, giving SNR 2 osc = ε 2 2 N periods SNR 2 bin ,(23) where SNR bin is the baseline SNR per log-period of width lnb: SNR 2 bin ≡ Z lnb [Ω 0 GW ] 2 σ 2 d lnf.(24) Hence SNR osc = ε √ 2 p N periods SNR bin ,(25) with N periods = ln(f max /f min ) lnb .(26) Throughout this paperSNR baseline ≡ SNR bin denotes the per-log-period baseline SNR; the total-band baseline SNR isSNR total = p N periods SNR bin . For LISA with effective band [f min ,f max ] = [10 −4 , 1]Hz(ln(f max /f min )≈9.21) the per-log-period baseline SNR at the WTC signal level isSNR bin ≈20–25; the factor p N periods / √ 2ranges from 2.6 atb= 2 to 1.7 atb= 5, providing meaningful amplification. We writeSNR osc ≈(ε/ √ 2)SNR bin p N periods in what follows; figures use SNR bin = 20 to set contours. The detectability plane (bvs.ε) is shown in Figs. 2 and 3, with SNR contours at {1, 5, 10, 20} and the WTC prediction band overlaid. 5 Ultraviolet Completion: Walking Technicolor Before developing the WTC embedding it is worth stating explicitly which claims in this paper rest on what kind of argument. The phenomenological backbone (Sec. 3–4) follows from the UETC ansatz and the controlled-approximation bounds of Sec. 3.2; the UV-completion layer developed in this section is a candidate realization motivated by holography and near-conformal dynamics but not derived from a complete microscopic model. Table 2 makes this hierarchy explicit. We emphasize that the cited literature [10, 13, 16, 17] motivates the individual ingre- dients — DSI in near-conformal systems, walking dynamics, holographic warp factors, radion potentials — but the complete chain from a microscopic WTC Lagrangian to the periodic technidilaton potential (29) is not, to our knowledge, established in the literature. The construction below should be read as a plausibility argument for the existence of a UV completion within the 5.1 Numerical validation As an independent check on the factorization argument we evaluate Eq. (4) numerically for a separable UETC Π(k,η 1 ,η 2 ) =F 0 (k)C(k)S(η 1 ,η 2 )g(η 1 −η 2 ) with a tophat macroscopic 9
23456 Discrete scaling factor b 0.01 0.05 0.10 0.50 Modulation amplitude SNR base = 20; LISA band [10 4 , 1] Hz Forecast SNR contours for DSI oscillations in the SGWB WTC [0.04, 0.18] b[1.7, 2.8] SNR=1 SNR=5 SNR=10 SNR=20 Figure 2: Forecast matched-filter SNR contours for the DSI oscillatory component in the (b,ε) plane, assumingSNR baseline = 20 and a LISA frequency band [10 −4 ,1]Hz. Contours are shown atSNR osc = 1,5,10,20. The orange shaded region is the WTC prediction band ε∈ [0.04, 0.18], b∈ [1.7, 2.8]. The model populates the high-SNR portion of the plane. 23456 Discrete scaling factor b 0.01 0.05 0.10 0.50 Modulation amplitude SNR base = 20; LISA band [10 4 , 1] Hz Forecast SNR contours with LISA 5 sensitivity and WTC prediction WTC [0.04, 0.18] b[1.7, 2.8] LISA 5 threshold (SNR osc = 5) LISA accessible (SNR osc 5) SNR=1 SNR=5 SNR=10 SNR=20 Figure 3: Same as Fig. 2, with the approximate LISA 5σdetection threshold (blue line, SNR osc = 5 forSNR baseline = 20) and LISA-accessible region (purple shading) overlaid. The WTC prediction band lies entirely within the LISA-accessible region. 10
Table 2: Scope and epistemic status of the principal claims of this paper. Rows 1–2 form the phenomenological backbone and are derived under explicit, quantified approximations. Rows 3–4 form the UV-completion layer and should be read as a candidate realization rather than a first-principles derivation. ClaimStatusSection DSI in UETC ⇒ log-periodic SGWB template derived3.2, 4.1 Factorization in short-correlation regimecontrolled approximation3.2, Table 1 WTC as DSI hostcandidate UV completion 5.3–5.5 Holographic origin of periodic warp factorconjectural motivation5.3 sourceSand a Gaussian temporal correlationg(∆η) of widthτ corr . Crucially the numerical computation uses the full Gaussian, not the δ-function limit invoked in Eq. (9). Figure 4 shows the resulting residualR(k) = (P h − P 0 h )/P 0 h alongside the analytic predictionε cos(2π ln(k/k ∗ )/ lnb+φ 0 ). Two features confirm the theorem: (i) atτ corr H ∗
0.05 the numerical residual matches the analytic template to within∼10 −16 , the floor of double-precision arithmetic; and (ii) repeating the calculation atτ corr H ∗ ∈{0.01,0.05,0.20} yields residuals that are pointwise identical to within numerical precision. Both observations follow from the fact that, whenC(k) depends only onkand is independent ofη 1 ,η 2 , it factors out of the double time integral algebraically — not merely up toO(τ corr H ∗ ). The error bounds in Table 1 arise instead from non-separable corrections (the convolution and Green’s-function terms): these are not probed by the present test and would require a more elaborate numerical setup. WTC framework, not as a derivation from first principles. 5.2 Phase-transition parameter space We adopt the benchmark large-N f QCD realization of walking technicolor [13]. The hidden sector is anSU(N c ) gauge theory withN f fundamental techniquarks in the near-conformal windowN f /N c ≳4–8. Near this window the gauge coupling walks — evolves slowly over many decades of energy scale — providing approximate scale invariance; the theory is attracted toward a quasi-fixed point (the Banks–Zaks fixed point) before condensing at Λ TC . Benchmark values areN c = 8,N f = 8, technidilaton decay constantF φ ≈1TeV, with an ultra-supercooled first-order phase transition (FOPT) characterized by [13] α≈ 0.73–0.83, β/H ∗ ≈ 100–1000, v w ≈ 1.(27) These give a sound-wave-dominated SGWB with h 2 Ω 0 GW (f peak )∼ 10 −9 –10 −8 at f peak ∼ 0.1–10 Hz,(28) comfortably within the LISA sensitivity band [1, 2], and satisfy τ corr H ∗ ≲ 0.01≪ 1. 5.3 Engineering discrete scale invariance Walking dynamics provide approximate continuous scale invariance: the technidilatonφis the pseudo-Nambu–Goldstone boson of the approximate scale symmetry, and its effective 11
10 0 10 1 wavenumber k (units with k * = 5) 0.15 0.10 0.05 0.00 0.05 0.10 0.15 R ( k ) ( P h P 0 h )/ P 0 h (a) = 0.10, b = 2.0, corr H * = 0.05, max dev. = 1.80e16 Factorization theorem: numerical residual vs.\ analytic template Analytic: cos(2ln(k/k * )/ln b + 0 ) Numerical: (P h P 0 h )/P 0 h 10 0 10 1 wavenumber k 0.15 0.10 0.05 0.00 0.05 0.10 0.15 R ( k ) (b) Independence of corr H * for separable UETCs: factorization is algebraic, not perturbative Analytic template Numerical ( corr H * = 0.01) Numerical ( corr H * = 0.05) Numerical ( corr H * = 0.2) Figure 4: Numerical validation of the factorization theorem (Sec. 3.2). The DSI-modulated tensor power spectrumP h (k) is computed by direct numerical integration of the double time integral (4) for a separable UETC Π =F 0 (k)C(k)S(η 1 ,η 2 )g(η 1 −η 2 ) with a Gaussian temporal correlationg(∆η) of widthτ corr (the full Gaussian; not theδ-function limit of Eq. (9)). (a) ResidualR(k) = (P h − P 0 h )/P 0 h (orange circles) plotted against the analytic predictionε cos(2π ln(k/k ∗ )/ lnb+φ 0 ) (gray line) atτ corr H ∗ = 0.05; maximum deviation ∼10 −16 (numerical floor). (b) The same residual evaluated atτ corr H ∗ ∈{0.01,0.05,0.20} collapses onto a single curve, confirming that the factorization is algebraic (independent ofτ corr H ∗ ) for separable UETCs. TheO(τ corr H ∗ ) corrections in Table 1 arise from non- separable structure (convolution and Green’s-function terms) that is beyond the scope of this clean test. 12
potential is of Coleman–Weinberg form [14]. DSI arises when this symmetry is broken from continuous to discrete. We realize this by adding a small explicit periodic modulation, V (φ) = V CW (φ) 1 + ε f cos 2π ln(φ/φ 0 ) lnb 0 , ε f ≪ 1, b 0
1,(29) whereV CW is the Coleman–Weinberg potential [14]. Such modulations are motivated by two independent sources. First, in AdS/CFT dual descriptions of near-conformal dynamics, periodic warp factors in the extra dimension generate exactly this type of potential modulation in the 4D effective theory. Second, near the quasi-fixed point the RGβ-function has no zero; instead the integrated RG flow over one cycle inφ-space is zero, corresponding to a limit cycle rather than a fixed point — the RG-flow realization of DSI [10]. Holographic origin ofε f andb 0 . Equation (29) arises naturally in the holographic dual of WTC without requiring fine-tuning of five-dimensional parameters. Modelling the technidilaton as the radion field in a Randall–Sundrum-type AdS 5 geometry [16], a small periodic modulation of the standard warp factor A(y) = ky of the form A(y) → ky + δA 0 sin(n p ky), δA 0 ≪ 1,(30) generates, at linear order inδA 0 , the multiplicative log-periodic correction to the 4D technidilaton potential [17] V 4 (φ)≈ V CW (φ) 1 + 4δA 0 cos 2π ln(φ/φ 0 ) lnb 0 +O(δA 2 0 ) ,(31) withε f = 4δA 0 andlnb 0 =kL/n p , whereLis the proper length of the extra dimension andn p is the number of warp-factor oscillation periods. For the ETC hierarchykL ≈ ln(Λ ETC /Λ TC )≈2–3 andn p = 3–4 (bothO(1) integers in AdS units), one obtains b 0 =e kL/n p ∈ [1.7,2.8]. The required amplitudeε f ∈[0.04,0.18] corresponds toδA 0 ∈ [0.01,0.045], a 1–5% warp-factor perturbation that is technically natural (protected by the approximate discreteφ→ φ+L/n p shift symmetry of the periodic modulation) and requires no independent fine-tuning. Both DSI parameters therefore emerge fromO(1) choices of the 5D geometry. The modulation in Eq. (29) induces a log-periodic correction to the gauge-field propa- gator at momentumq. We derive this at leading order inε f . In the near-conformal WTC regime, the technidilaton VEV⟨φ(q)⟩at renormalisation scaleqis related to its UV value by ⟨φ(q)⟩=⟨φ UV ⟩(q/q 0 ) −∆ φ , where ∆ φ is the technidilaton scaling dimension (∆ φ ≈1 near the quasi-fixed point). The gauge-boson mass is generated viam 2 V (φ) =y 2 ⟨φ⟩ 2 ; a modulation δV ∝ ε f cos(2π lnφ/ lnb 0 ) shifts the mass asδm 2 V /m 2 V =ε f cos(2π ln(q/q ∗ )/ lnb 0 ) +O(ε 2 f ) by the chain rule. Propagating to the full propagator at leading order in ε f : D(q; ∆η) = D 0 (q; ∆η) [1 + δ(q)], δ(q) = ε f cos 2π ln(q/q ∗ ) lnb 0
- φ 0 ,(32) withq ∗ ∼ q 0 . The modulation inherits the same log-periodb 0 as the potential, up to the conformal-dimension factor ∆ φ which isO(1) near the fixed point. Higher-order corrections enter at O(ε 2 f ). Note that |δ(q)|≤ ε f ≪ 1 so D is positive definite for all q. 13
5.4 Convolution for the UETC The transverse-traceless anisotropic stress is bilinear in the gauge fields, so the UETC is the convolution Π(k,η,η ′ )∝ Z d 3 p (2π) 3 P TT D(p; ∆η)D(|k−p|; ∆η),(33) whereP TT projects onto the transverse-traceless sector. Expanding to linear order inε f and retaining only the cross-term (the self-term is O(ε 2 f )), Π(k)⊃ Z d 3 p (2π) 3 P TT D 0 (p)D 0 (|k−p|) δ(p) + δ(|k−p|) .(34) We bound the two terms separately. Termδ(p). The baseline propagatorD 0 (p) is sharply peaked atp∼ q ∗ ∼ β/v w with relative half-width ∆p/p∼ τ corr H ∗ ≪1. The variation ofδacross this peak is bounded using the absolute derivative: dδ d lnq = ε f 2π lnb 0 sin
2π lnq/ lnb 0
- φ 0 ≤ ε f 2π lnb 0 ,(35) uniformly bounded for allq, including near the zeros ofδ(where the logarithmic derivative d lnδ/d lnqwould diverge, but the divergence is integrable since the weightD 0 is smooth and |δ|→ 0). Integrating over the support | ln(p/k)|≲ τ corr H ∗ : |δ(p)− δ(k)|≲ 2πε f lnb 0 τ corr H ∗ ≲ 13ε f τ corr H ∗ (b 0 ≥ 1.5),(36) giving δ(p) = δ(k)[1 +O(ε f τ corr H ∗ )] over the support of D 0 (p). Termδ(|k−p|). Nearp≈k,|k−p|→0 andδ(|k−p|) oscillates rapidly. However, WTC gauge bosons have a mass gapm V ∼Λ TC , and for|k−p|≪ m V the propagator satisfiesD 0 (|k −p|)≲ D 0 (0)∼1/m 2 V . At the signal scalek ∼ q ∗ ≫ m V ,D 0 (q ∗ )∼ 1/(q 2 ∗ − m 2 V )≪1/m 2 V , soD 0 (0)/D 0 (q ∗ )∼ m 2 V /q 2 ∗ ≡ ρ ≪1. The contribution of the |k−p| →0 region is therefore suppressed byρrelative to the peak contribution from |k−p|∼ q ∗ , where the same slow-variation argument as forδ(p) applies. The combined bound is δ(p)≈ δ(|k−p|)≈ δ(k) 1 +O(ε f τ corr H ∗ ) +O(ρ) ,(37) withρ≡ m 2 V /q 2 ∗ ≪ 1 in the WTC regime. The cross-term factorizes as 2δ(k)×Π 0 (k,η,η ′ ). The factor of 2 corresponds to the two linear cross-terms; absorbing it into a redefined baseline normalization setsε=ε f at leading order (the factor of 2 multiplies the baseline amplitude, not the modulation depth). Hence Π(k,η,η ′ ) = Π 0 (k,η,η ′ ) 1 + ε cos 2π ln(k/k ∗ ) lnb
- φ 0 1 +O(ετ corr H ∗ ) +O(ρ) , (38) withε=ε f andb=b 0 at leading order. For WTC benchmark parametersβ/H ∗ ≳100 and ρ∼ m 2 V v 2 w /β 2 ≲ 0.01, the combined relative correction is≲ 1.4% (Table 1). 14
5.5 WTC predictions Combining Eq. (38) with the short-correlation-time factorization theorem of Sec. 3.2, the DSI modulation propagates multiplicatively to the observable SGWB, recovering Eq. (16) up to controlledO(ετ corr H ∗ ) corrections (andO(ρ) mass-gap corrections from the WTC convolution, withρ ≡ m 2 V /q 2 ∗ ≲ 0.01). The WTC parameter space [13], spanned by F φ ≈ 1 TeV, Λ ETC ∼ 5–10 TeV, and soft masses m p ∼ 1–100 GeV, maps onto ε∈ [0.04, 0.18], b∈ [1.7, 2.8].(39) This band is shown in Figs. 2–3 and overlaps the high-SNR region of the LISA detectability forecast. The map from WTC parameters to (b,ε) follows from the holographic benchmark of Sec. 5.3:ε= 4δA 0 andb=exp(kL/n p ), withkL≈ ln(Λ ETC /Λ TC )≈ ln(5–10) = 1.6–2.3 andn p = 3–4 from the WTC benchmark [13];δA 0 ∈[0.01,0.045] from theO(1–5%) warp-factor perturbation range. For a baselineSNR bin = 20 the matched-filter SNR satisfies SNR osc ≳ 4 over most of the band (using Eq. (25) with the 1/ √ 2 factor). The chain from the WTC Lagrangian to the observable Ω GW (f) is now complete: every step has been individually justified and the cumulative relative error is below 2% for β/H ∗ ≳ 100. 6 Discussion Robustness of the factorization. The key approximation is the short-correlation-time limitτ corr H ∗ ≪1. Its validity requiresβ/H ∗ ≫1, i.e. a transition that completes rapidly compared to the Hubble time. This is satisfied for the WTC benchmark (β/H ∗ ∼100– 1000) and is a generic property of strong first-order transitions. Slow transitions with β/H ∗ ≲10 would require higher-order corrections, which can be computed systematically as an expansion inτ corr H ∗ . The separate factorization condition|dδ/d lnq|·(∆q/q)≪1 — using the absolute derivative|dδ/d lnq| ≤ ε f ×2π/ lnb 0 , which is uniformly bounded for allq(the logarithmic derivatived lnδ/d lnqwould diverge at the zeros ofδ, but the absolute derivative does not; see Sec. 5.4) — is equally well controlled and introduces no additional tuning. FDTD analogy and cosmological causal structure. The companion FDTD pa- per [18] demonstrates log-periodic spectral imprinting in a controlled electromagnetic cavity; the correspondence to the cosmological FOPT warrants explicit comment. In the cavity, rigid static boundaries enforce global mode selection via discrete standing-wave conditions (Dirichlet or absorbing boundary conditions at the walls): the mode spectrum is shaped by the entire geometry simultaneously. In a FOPT no global boundary condition exists: bubbles nucleate independently within their past light cones, and causal horizons preclude global mode coherence. The structural role of the geometric boundary is instead played by the characteristic bubble spacingR ∗ ∼ v w /β, which acts as a local, dynamic filter. Modes withk ≫ R −1 ∗ are exponentially suppressed by the decay of the temporal correlation functionF(k,∆η) at large separations; modes withk ≪ R −1 ∗ see a nearly homo- geneous source and are coherently accessible. In the language of the companion framework paper [19],R ∗ plays the role of the boundary-conditioned mode density cutoff, with the plasma mean free path providing the dynamic spectral participation filter. Crucially, the factorization theorem of Sec. 3.2 relies only on the local conditionτ corr H ∗ ≪1 — set by β/H ∗ ≫1, independent of any global causal horizon structure. The FDTD result therefore 15
validates the mathematical mechanism of DSI imprinting (that a log-periodically structured boundary parameter transfers its signature multiplicatively to the power spectrum), while the factorization theorem independently establishes the validity of that transfer in the cosmological context via purely local causal arguments. Distinguishability from other spectral features. The log-periodic modulation (16) produces a coherent, phase-stable sinusoid inlnf, persisting overN periods ∼6–13 full oscillations across the LISA band forb∈[1.7,2.8]. This is qualitatively distinct from other known spectral features: (i) The kink at the crossover from sound-wave to turbulence domination is a single discontinuity in the spectral slope, not a periodic oscillation. (ii) A sharp bubble-collision peak is a feature of limited frequency extent, not a multi-period sinusoid. (iii) Stochastic backgrounds from astrophysical sources produce spectra that are smooth inlnfto high accuracy. A likelihood-ratio test between the smooth template Ω 0 GW and the DSI-modulated template (16) provides the optimal discriminant. The three-parameter family (b,ε,φ 0 ) can be mapped from the data by standard matched-filter techniques [2]. Parameter degeneracies. The phaseφ 0 merely shifts the oscillation inlnfand does not affect detectability;εandbcan be independently constrained from the oscillation depth and period respectively. The frequency resolution needed to resolve individual oscillations is ∆f/f ∼ lnb/(2π); forb= 2 this is ∆f/f ≈0.11, well within LISA’s capabilities over its four-year nominal mission. Alternative UV completions. The factorization result and the observable template (16) are model-independent consequences of DSI in the UETC, requiring onlyτ corr H ∗ ≪1. Walking technicolor is one concrete realization; other BSM models with approximate conformal symmetry and explicit periodic modulations — extended Higgs sectors with Coleman–Weinberg potentials modified by threshold corrections, Randall–Sundrum–type models with periodic radion potentials, or clockwork models [15] — are equally valid candidates and will produce the same spectral template with different (b,ε) values. A detection of log-periodic oscillations in the SGWB would uniquely fixbandε, allowing discrimination among UV completions. Multi-messenger signatures. Beyond gravitational waves, the DSI in the WTC potential generates log-periodic modulations in the technidilaton production rate and hence in the energy density of any dark-radiation component coupled to the hidden sector, providing in principle an independent observational handle on the same (ε,b) parameters. 7 Conclusions We have demonstrated that discrete scale invariance in the anisotropic stress tensor of a first-order cosmological phase transition imprints a multiplicative log-periodic modulation on the stochastic gravitational-wave background. The main results are: 1. Factorization theorem. In the physically motivated short-correlation-time limit (τ corr H ∗ ≪1, satisfied forβ/H ∗ ≳10), the DSI modulation passes from the source UETC to the observable Ω GW (f) at the per-cent level:P h =C(k)P 0 h [1 +O(τ corr H ∗ )]. 16
2.Universal spectral template. The observable signature is Ω GW = Ω 0 GW [1 + ε cos(2π ln(f/f ∗ )/ lnb+φ 0 )] — a sinusoid inlnfsuperimposed on the smooth baseline, characterized at leading order in τ corr H ∗ by three parameters (ε,b,φ 0 ). 3.Matched-filter detectability.SNR osc ≃(ε/ √ 2)SNR baseline p N periods , withN periods
6–13 oscillations in the LISA band. 4.Walking technicolor UV completion. An explicit convolution calculation con- firms that the WTC potential modulation propagates to the SGWB with≲1% error. The WTC prediction bandε ∈[0.04,0.18],b ∈[1.7,2.8] sits in the high-SNR osc region of the LISA detectability plane. A non-detection by LISA would place sharp upper limits onεas a function ofb, directly constraining the allowed parameter space for near-conformal BSM phase transitions. A detection would simultaneously reveal the discrete scaling ratio, the DSI amplitude, and the phase of the modulation, providing a unique window into the self-similar structure of the hidden-sector dynamics. The log-periodic template (16) is simple, well-defined, and implementable in any LISA data-analysis pipeline via standard matched-filter methods. Acknowledgments The author thanks the gravitational-wave and beyond-Standard-Model communities for stimulating discussions. No external funding was received for this work. References [1]P. Amaro-Seoane et al. (LISA Collaboration), “Laser Interferometer Space Antenna,” (2017) [arXiv:1702.00786]. [2] C. Caprini et al., “Science with the space-based interferometer eLISA. II: Grav- itational waves from cosmological phase transitions,” JCAP 1604, 001 (2016) [arXiv:1512.06239]. [3]G. Agazie et al. (NANOGrav Collaboration), “The NANOGrav 15 yr Data Set: Evidence for a Gravitational-Wave Background,” Astrophys. J. Lett. 951, L8 (2023) [arXiv:2306.16213]. [4]D. J. Reardon et al. (PPTA Collaboration), “Search for an Isotropic Gravitational- Wave Background with the Parkes Pulsar Timing Array,” Astrophys. J. Lett. 951, L6 (2023) [arXiv:2306.16215]. [5]J. Antoniadis et al. (EPTA Collaboration), “The second data release from the European Pulsar Timing Array: V. Implications for massive black holes, dark matter and the early Universe,” Astron. Astrophys. 678, A50 (2023) [arXiv:2306.16227]. [6]M. Hindmarsh and M. Hijazi, “Gravitational waves from first-order cosmological phase transitions in the Sound Shell Model,” JCAP 12, 062 (2019) [arXiv:1909.10040]. [7] D. G. Figueroa, A. Florio, F. Guedes, and F. Torrenti, “Cosmological phase tran- sitions: From theory to gravitational wave phenomenology,” JCAP 03, 027 (2021) [arXiv:2010.00972]. 17
[8]J. R. Espinosa, T. Konstandin, J. M. No, and G. Servant, “Energy Budget of Cosmological First-Order Phase Transitions,” JCAP 06, 028 (2010) [arXiv:1004.0691]. [9]C. Caprini and R. Durrer, “Gravitational waves from stochastic relativistic sources: Primordial turbulence and magnetic fields,” Phys. Rev. D 74, 063521 (2006) [arXiv:astro-ph/0603476]. [10]D. Sornette, “Discrete scale invariance and complex dimensions,” Phys. Rep. 297, 239 (1998) [arXiv:cond-mat/9707012]. [11]D. Sornette, Critical Phenomena in Natural Sciences: Chaos, Fractals, Self- Organization and Disorder: Concepts and Tools, 2nd ed. (Springer, 2006). [12]G. Calcagni and S. Kuroyanagi, “Log-periodic gravitational-wave background beyond Einstein gravity,” Class. Quantum Grav. 41, 015031 (2024) [arXiv:2308.05904]. [13]M. Miura, K. Ohnishi, T. Sawanaka, and K. Yamawaki, “Gravitational waves from walking technicolor,” (2019) [arXiv:1811.05670]. [14] S. Coleman and E. Weinberg, “Radiative corrections as the origin of spontaneous symmetry breaking,” Phys. Rev. D 7, 1888 (1973). [15]G. F. Giudice, Y. Kats, M. McCullough, R. Torre, and A. Urbano, “Clockwork/linear dilaton: structure and phenomenology,” JHEP 06, 098 (2018) [arXiv:1711.08437]. [16]L. Randall and R. Sundrum, “A large mass hierarchy from a small extra dimension,” Phys. Rev. Lett. 83, 3370 (1999) [arXiv:hep-ph/9905221]. [17] W. D. Goldberger and M. B. Wise, “Moduli stabilization with bulk fields,” Phys. Rev. Lett. 83, 4922 (1999) [arXiv:hep-ph/9907447]. [18]J. F. Rankin, “Log-periodic spectral hierarchies in a boundary-driven electromagnetic cavity: evidence from FDTD simulations,” preprint (2026). [19]J. F. Rankin, “A phenomenological framework for mode-accessibility engineering in structured field environments,” preprint (2026). 18
This paper presents a thorough theoretical analysis of discrete scale invariance signatures in gravitational wave backgrounds. The core phenomenological framework is complete, with the factorization theorem rigorously derived and error bounds explicitly quantified. The authors clearly distinguish between the well-controlled theoretical backbone (Sections 3-4) and the candidate UV completion (Section 5), appropriately characterizing the epistemic status of different claims. While the holographic motivation for the WTC embedding involves some conjectural elements, the mathematical chain from assumed periodic potential to observable predictions is complete and well-bounded. The work successfully addresses its stated goals within its declared framework.
This paper is substantially more complete than a speculative note: it presents a followable core derivation for how a log-periodic modulation in the source UETC propagates into the SGWB spectrum, and it gives an explicit observable template plus a matched-filter scaling argument. The phenomenological heart of the paper is therefore reasonably self-contained and supported within its own approximation scheme.
The main limitations lie in the transition from that phenomenological core to the advertised UV completion and quantitative error control. The paper is commendably honest that the WTC realization is a candidate rather than a first-principles derivation, but that also means one of its headline claims is only partially realized. Likewise, the strongest numerical check validates a separable special case, not the more difficult corrections responsible for the quoted percent-level bounds. Overall, the submission is coherent and mostly complete at the phenomenological level, but only moderately complete as a full end-to-end supported paper because the UV-completion and error-validation layers remain partly schematic.
The paper's internal architecture is presented clearly and its claim hierarchy is well-flagged, but its central technical claims hinge on a definitional bridge between the Sec. 3 phenomenological ansatz (a k-only multiplicative C(k) on the UETC) and the Sec. 5 WTC construction (modulation entering via internal-momentum convolution). This bridge is asserted under a peaking/dominance argument that is not made quantitative, and is accompanied by a factor-of-2 absorption with simultaneous identification ε = ε_f that is not algebraically justified. Both are load-bearing for the paper's headline falsifiable target (ε ∈ [0.04, 0.18], b ∈ [1.7, 2.8] in LISA's high-SNR region). Additionally, the short-correlation approximation introduced at Eq. (9) with declared O(τ_corr H_) error is treated as exact in the observable spectrum and SNR analysis, and the abstract's percent-level claim for β/H_ ≳ 10 conflicts with the reported ~15% error at β/H_* = 10 in Table 1. After weighing the optimistic 5/5 and 4/5 assessments — which emphasize the surface coherence of the ansatz and notational hygiene — against the structural issues identified by the two 2/5 assessments, I conclude the lower scores are correct. The red-flag caps for central definition drift, approximation escalation, and unverified central derivation all apply.
⚑Derivation Flags (28)
- highEq. (29) ⇒ Eq. (32) — Claim that a periodic modulation of the technidilaton potential induces δmV^2/mV^2=εf cos(2π ln(q/q*)/ln b0)+O(εf^2) and hence a multiplicative propagator modulation D=D0(1+δ(q)) is not derived; it assumes a specific dependence of masses/propagators on φ(q) and a direct ‘chain rule’ transfer of the cosine.
If wrong: If δ(q) is not of the asserted form/amplitude, the subsequent convolution does not yield Eq. (38), and the claimed mapping from WTC/holographic parameters to the SGWB modulation (ε,b) is unsupported.
- highEq. (31), Sec. 5.3 — Multiplicative log-periodic 4D potential V_4(φ) ≈ V_CW(φ)[1 + 4δA_0 cos(2π ln(φ/φ_0)/ln b_0)] derived from a sinusoidally perturbed AdS_5 warp factor is stated with citation to [17] (Goldberger-Wise) but no explicit derivation. The factor of 4, the cosine-in-log-φ structure, and the identification ln b_0 = kL/n_p are presented as results.
If wrong: The mapping from WTC/RS holographic parameters (δA_0, kL, n_p) to the observable (ε, b) band would be unsupported. The central falsifiable WTC prediction ε ∈ [0.04, 0.18], b ∈ [1.7, 2.8] would lose its microphysical anchor, reducing the UV-completion section to dimensional estimates.
- highEq. (32) — The derivation from a log-periodic technidilaton potential to a gauge-field propagator modulation D(q)=D_0(q)[1+delta(q)] is highly compressed. The asserted same-amplitude transfer by the chain rule, the treatment of the scaling dimension Delta_phi, and the positivity claim for D are not demonstrated.
If wrong: The starting point for the WTC UETC convolution is unsupported; Eq. (38) and the WTC prediction band Eq. (39) would not follow.
- highEq. (33)–(37) — Convolution factorization relies on (a) D0(p) being sharply peaked with width Δp/p∼τcorr H*; (b) bounding δ(p) variation by |dδ/d ln q| times support width; (c) suppressing δ(|k−p|) near |k−p|→0 by a mass-gap ratio ρ; these are plausible but not established with explicit models for D0 and the TT projector weighting.
If wrong: If D0 is not sharply peaked or the TT projection weights regions where δ varies rapidly, the modulation may be smeared, phase-shifted, or reduced; then Eq. (38) and the ≲1% propagation claim fail.
- highEq. (37) — The mass-gap suppression estimate appears algebraically reversed. For a massive propagator of the form 1/(q^2+m_V^2), D_0(0)/D_0(q_*) is approximately q_*^2/m_V^2, not m_V^2/q_*^2. The text claims the small ratio rho=m_V^2/q_*^2 suppresses the |k-p|->0 region.
If wrong: The bound excluding rapidly oscillatory small-|k-p| contributions is invalid; the convolution factorization in Eq. (38) is not reliable.
- highEq. (38) — The statement that the factor of 2 from the two linear cross terms can be absorbed into the baseline normalization while keeping epsilon=epsilon_f is not mathematically justified. A multiplicative baseline rescaling does not remove a relative modulation-depth factor.
If wrong: The predicted modulation amplitude range Eq. (39) is at least quantitatively wrong, and the WTC detectability conclusions shift.
- highEq. (39) — The map from WTC parameters to epsilon in [0.04,0.18] and b in [1.7,2.8] depends on the unproven holographic mapping, the compressed propagator derivation, and the problematic convolution factorization.
If wrong: The paper's falsifiable WTC prediction band is unsupported, although the generic phenomenological template remains definable.
- highEq. (7) and transition to Eq. (9) — Replacement of F(k,Δη) by F(k)δ(Δη)+O(τcorr H*) is sketched via a ‘nascent delta function’ expansion; no explicit control norm or assumptions on F are given to justify the stated relative error scaling for the GW integral.
If wrong: If the δ-approximation does not control the double-time integral in Eq. (4) as claimed, then Eq. (10)–(15) (factorization theorem and error bounds) are unreliable, undermining the main phenomenological claim that the modulation survives to percent level.
- highEq. (9) — The replacement F(k,Delta eta) ≃ F(k) delta(eta-eta') + O(tau_corr H_*) is asserted via a nascent delta expansion, but the normalization and dimensional transition from F(k,0) delta_tau to F(k)=integral F dDelta eta are not fully derived. The estimate also requires control of variation of the rest of the integrand over tau_corr.
If wrong: The claimed O(tau_corr H_*) factorization error in the central theorem is not established; Eq. (15)'s accuracy bound would be unsupported.
- highEqs. (12)-(14) — The bound on the oscillatory Green-function contribution assumes |partial_eta(a^4 S)| less than about H_* a^4 S. For a first-order phase transition, the source envelope generally turns on/off on beta^{-1}, not necessarily H_*^{-1}. The boundary-term and smooth-envelope assumptions are not proved for the stated source models.
If wrong: The conclusion that the Green-function contribution is only O(H_*/k) and cannot contaminate the log-periodic modulation is not quantitatively established.
- highEqs. (34)-(36) — The replacement delta(p) ≈ delta(k) assumes the convolution support satisfies |ln(p/k)| less than about tau_corr H_*. The paper only states that D_0(p) is peaked at p~q_*; it does not derive that p tracks the external momentum k over the spectral range. In a three-dimensional convolution, p, |k-p|, and k are geometrically distinct.
If wrong: The convolution need not preserve the same log-periodic phase and amplitude as a function of external k; the claimed UETC ansatz Eq. (38) may fail.
- highSec. 3.2, statement before Eq. (9) — The paper states that G_k(eta,eta_2)a^2(eta_2) varies on timescale k^{-1} >> tau_corr for subhorizon modes. But for the characteristic phase-transition signal k is later taken to be of order beta/v_w, while tau_corr ~ v_w/beta, giving k tau_corr ~ 1 rather than k^{-1} >> tau_corr.
If wrong: The delta-function reduction of the unequal-time integral cannot be justified by slow variation of the Green function; the central short-correlation-time error estimate may fail.
- mediumEq. (10)–(11) — After inserting the δ-function, the resulting single-integral expression uses S(η1,η1) and defines F(k)=∫F(k,Δη)dΔη, but the mapping between the original Π0(k,η1,η2) and the new objects F0(k), S(η1,η1) is not derived; it is asserted as a decomposition.
If wrong: If Π0 does not admit this reduction, P0h(k,η) may inherit additional k-structure that can interfere with or mimic log-periodicity, weakening the separation claimed in Eq. (15)–(16).
- mediumEq. (12)–(14) — Integration-by-parts bound on Ĩ(k,η) assumes vanishing boundary terms ‘since S=0 outside the source epoch’ and uses |∂η1(a^4 S)|≲H* a^4 S without specifying conditions under which this inequality holds (e.g., sharp turn-on/off, differentiability, and whether H* should be conformal H).
If wrong: If boundary terms or envelope derivatives are not controlled, the estimate |Ĩ|/I0≲H*/k may fail, so the claim that Green’s-function oscillations cannot generate log-periodic structure at period ln b is not secured.
- mediumEq. (14), relation H_*/k ~ v_w/(beta/H_*) — The relation mixes conformal variables from Eq. (2), physical Hubble quantities used in beta/H_*, and observed/comoving wave numbers without a consistent scale-factor derivation.
If wrong: The numerical error estimates in Table 1 and the percent-level accuracy claim may be mis-scaled.
- mediumEq. (16)–(17) — Template written as exact multiplicative modulation and ‘exact sinusoid in ln f’ downstream, despite earlier O(τcorr H*) and later O(ε τcorr H*)+O(ρ) corrections.
If wrong: Data-analysis claims (perfect template correlation, exact sinusoidality) would be overstated; practical matched-filter performance could degrade or require additional nuisance parameters.
- mediumEq. (18) — The adopted sound-wave amplitude is given with (H_*/beta)^2. The standard formula usually contains a single power H_*/beta unless an additional finite-lifetime suppression is separately included and derived. No derivation of the squared dependence is provided.
If wrong: The baseline amplitude, LISA SNR estimates, and placement of the WTC band in the detectability plots could be quantitatively wrong.
- mediumEq. (32), Sec. 5.3 — Propagation of the potential modulation to the gauge-boson propagator via δm_V^2/m_V^2 = ε_f cos(2π ln(q/q_*)/ln b_0). The chain-rule argument substitutes the running VEV ⟨φ(q)⟩ ∝ q^{-Δ_φ} into the modulation, but the conversion of a cos(ln φ) modulation into a cos(ln q) modulation with the same period b_0 (up to an O(1) Δ_φ factor) is asserted rather than carried through with the dimensional factor made explicit.
If wrong: The identification b = b_0 at leading order in the WTC convolution (Eq. 38) would be modified by a Δ_φ-dependent factor, shifting the WTC b-band away from [1.7, 2.8]. The qualitative existence of the log-periodic modulation in the propagator would survive.
- mediumEq. (38) and statement 'absorbing [factor 2] into a redefined baseline normalization sets ε=εf at leading order' — The identification of the UETC modulation depth ε with εf is asserted. In a generic structure Π=Π0+εf Π1, the relative modulation depth depends on Π1/Π0; the factor 2 alone does not guarantee ε=εf.
If wrong: The predicted ε range in Eq. (39) could shift systematically (e.g., by an O(1) factor), affecting the claimed occupancy of the high-SNR region.
- mediumEq. (6) — DSI introduced as a purely k-dependent multiplicative modulation of the full UETC Π(k,η,η′) without showing that such factorization is compatible with general UETC structure for the relevant phase-transition sources.
If wrong: If the DSI modulation depends on (η,η′) as well as k, it will not factor out of the unequal-time integrals as used in Eq. (10)–(15), and the clean template Eq. (16) may be distorted or suppressed.
- mediumEqs. (30)-(31) — The mapping from a periodic 5D warp-factor perturbation A(y)->ky+delta A_0 sin(n_p k y) to the 4D technidilaton potential modulation with coefficient 4 delta A_0 and period ln b_0=kL/n_p is stated without derivation.
If wrong: The WTC/holographic origin of the predicted b and epsilon ranges is not established, though the phenomenological DSI template would remain intact.
- mediumSec. 5.1 / Fig. 4 — The numerical validation uses a separable UETC with C(k) independent of both times, for which factorization is algebraic. It explicitly does not test the non-separable convolution and Green-function corrections that determine the claimed percent-level error bound.
If wrong: The numerical experiment cannot validate the central approximation bounds; it only confirms the trivial separable case.
- mediumSec. 5.4, the claim 'the baseline propagator D_0(p) is sharply peaked at p ~ q_* with relative half-width Δp/p ~ τ_corr H_*' — The narrow-width property of the propagator in momentum space is asserted without derivation from the WTC dynamics. The statement appears to conflate temporal correlation length (τ_corr) with spatial momentum width; a rigorous justification of Δp/p is missing.
If wrong: If the peak is significantly broader, the bound on |δ(p)−δ(k)| could be larger than estimated, potentially invalidating the ≤1.4% error claim for the WTC convolution. The predicted ε, b band for the WTC model would be less reliable, although the model-independent DSI imprinting mechanism (factorization theorem) remains unaffected.
- lowEq. (13) and identification H*/k ≡ τcorr H* — The step relating H*/k to v_w/(β/H*) uses an asserted characteristic wavenumber k∼β a*/(v_w a0) but mixes quantities at different epochs and does not show the correct dimensionless ratio in conformal variables.
If wrong: Quantitative error bars in Table 1 and the claimed percent-level bounds could shift; the qualitative suppression with k≫H* may still hold, but the numerical mapping to β/H* could be off.
- lowEq. (13), Sec. 3.2 — Integration-by-parts bound |Ĩ(k,η)|/I_0 ≲ H_*/k uses |∂_{η_1}[a^4 S]| ≲ H_* a^4 S, which is plausible for a slowly varying envelope but not explicitly justified for the source envelope shape during a FOPT (where S can have sharp turn-on/turn-off features).
If wrong: The O(τ_corr H_*) bound on the Green's-function-induced log-periodic contamination could be weaker by an O(1) factor, slightly degrading the percent-level factorization accuracy but not invalidating the theorem.
- lowEq. (22)–(26) — Derivation of SNRosc scaling uses ⟨cos^2⟩=1/2 over Nperiods complete periods but does not address weighting by σ(f) and Ω0GW(f) variations across the band; Eq. (23)–(25) effectively assumes the integral decomposes into equal-contribution periods.
If wrong: The √Nperiods enhancement may be weaker/stronger depending on band-shape; this affects quantitative detectability contours but not the existence of an oscillatory matched-filter statistic.
- lowEq. (36), Sec. 5.4 — Numerical bound |δ(p) - δ(k)| ≲ 13 ε_f τ_corr H_* for b_0 ≥ 1.5 uses the absolute derivative bound 2π/ln b_0. The relative half-width ∆p/p ~ τ_corr H_* of the propagator peak is asserted but not derived from the propagator structure.
If wrong: Numerical coefficient in the convolution error bound could change by an O(1) factor; the qualitative O(ε_f τ_corr H_*) scaling would survive.
- lowTable 1 combined bound and Sec. 5.1 numerical validation — Numerical test uses a separable Π=F0(k)C(k)S(η1,η2)g(η1−η2), which algebraically factorizes C(k); it does not probe the non-separable corrections that are said to dominate the error budget (convolution, Green’s-function term).
If wrong: The presented ‘validation’ does not validate the stated percent-level error bounds for realistic non-separable sources; thus the empirical support for Table 1 is weaker than implied.
This is a scientifically interesting and reasonably original paper whose strongest contribution is phenomenological: it proposes that discrete scale invariance in the source of a first-order phase transition yields a log-periodic SGWB modulation that could be searched for with matched filtering. That signature is specific, observationally meaningful, and different enough from standard smooth spectral templates to make the work legitimately testable. The manuscript also does better than many speculative papers at marking which claims are derived and which are conjectural.
The main weakness is communication discipline around the UV completion and validation claims. The body carefully says the walking-technicolor construction is only a candidate realization, but the abstract presents it more definitively. Likewise, the numerical check confirms a simplified separable setup rather than the hardest part of the claimed theorem. So the paper's core idea appears novel and potentially useful, but its strongest rhetoric should be softened to match what is actually established. Overall: worthwhile phenomenological proposal with good testability, solid but not maximal novelty, and moderate clarity limited by some overclaim and notation ambiguity.
This is a scientifically well-crafted submission that combines a clean factorization theorem (derived under explicit, quantified approximations valid for realistic β/H* ≳ 10 transitions) with a candidate UV completion (walking technicolor with holographically motivated periodic warp factor) to produce a sharp, falsifiable observational prediction for LISA. The author is unusually careful to distinguish what is derived from what is conjectured (Table 2), and to bound every approximation explicitly (Table 1). The observable template — a log-periodic sinusoid in ln(f) with amplitude ε ∈ [0.04, 0.18] and ratio b ∈ [1.7, 2.8] — is detectable by matched-filter analysis with SNR_osc ≳ 4 in the LISA band, providing a clean discriminant against competing spectral features (sound-wave kinks, bubble-collision peaks, astrophysical backgrounds).
The principal limitations are appropriately self-flagged: the microscopic derivation of the periodic technidilaton potential from a full WTC Lagrangian is not established, and the numerical validation tests only the separable-UETC limit. Neither limitation undermines the central phenomenological backbone, because the factorization theorem and the observable template are model-independent consequences of DSI in the source UETC — WTC is offered as a representative, not unique, UV completion (clockwork and Randall–Sundrum alternatives are explicitly noted). The submission exemplifies the discipline of stating what one has derived, what one has approximated, and what one has conjectured, and it deserves to be evaluated as a serious theoretical-phenomenological proposal with a clear experimental target.
The paper is exceptionally well-structured and self-contained. It defines the DSI ansatz for the source UETC, rigorously derives the factorization theorem under the short-correlation-time approximation, bounds all errors, and calculates the observable log-periodic spectrum and its matched-filter detectability. A walking technicolor UV completion is presented as a candidate realization, with the convolution from the technidilaton potential to the SGWB explicitly worked out and its error estimated. Every approximation is quantified, and the epistemic status of each component is clearly labelled. The work fully addresses its stated objectives and includes discussions of robustness, distinguishability, and alternative completions. No missing variables, skipped steps, or unaddressed goals compromise the completeness.
Internally, the paper is coherent in its phenomenological core: if one assumes the Sec. 3 DSI ansatz Π=Π0·C(k), then the subsequent use of linear response and short-correlation-time reasoning can consistently motivate an approximately multiplicative log-periodic modulation in the observable spectrum.
The main internal-consistency failure arises when the paper attempts to weld the UV-completion construction onto that ansatz. Introducing DSI at the potential/propagator level and then claiming the resulting anisotropic-stress UETC is modulated by the same external C(k) requires an explicit equivalence argument or a quantified remainder term; without it, the narrative implicitly shifts the meaning/realization of the central modulation object while continuing to apply conclusions derived under the earlier meaning. This is central (not peripheral) because it underwrites the claimed WTC-derived falsifiable (ε,b) band.
⚑Derivation Flags (28)
- highEq. (29) ⇒ Eq. (32) — Claim that a periodic modulation of the technidilaton potential induces δmV^2/mV^2=εf cos(2π ln(q/q*)/ln b0)+O(εf^2) and hence a multiplicative propagator modulation D=D0(1+δ(q)) is not derived; it assumes a specific dependence of masses/propagators on φ(q) and a direct ‘chain rule’ transfer of the cosine.
If wrong: If δ(q) is not of the asserted form/amplitude, the subsequent convolution does not yield Eq. (38), and the claimed mapping from WTC/holographic parameters to the SGWB modulation (ε,b) is unsupported.
- highEq. (31), Sec. 5.3 — Multiplicative log-periodic 4D potential V_4(φ) ≈ V_CW(φ)[1 + 4δA_0 cos(2π ln(φ/φ_0)/ln b_0)] derived from a sinusoidally perturbed AdS_5 warp factor is stated with citation to [17] (Goldberger-Wise) but no explicit derivation. The factor of 4, the cosine-in-log-φ structure, and the identification ln b_0 = kL/n_p are presented as results.
If wrong: The mapping from WTC/RS holographic parameters (δA_0, kL, n_p) to the observable (ε, b) band would be unsupported. The central falsifiable WTC prediction ε ∈ [0.04, 0.18], b ∈ [1.7, 2.8] would lose its microphysical anchor, reducing the UV-completion section to dimensional estimates.
- highEq. (32) — The derivation from a log-periodic technidilaton potential to a gauge-field propagator modulation D(q)=D_0(q)[1+delta(q)] is highly compressed. The asserted same-amplitude transfer by the chain rule, the treatment of the scaling dimension Delta_phi, and the positivity claim for D are not demonstrated.
If wrong: The starting point for the WTC UETC convolution is unsupported; Eq. (38) and the WTC prediction band Eq. (39) would not follow.
- highEq. (33)–(37) — Convolution factorization relies on (a) D0(p) being sharply peaked with width Δp/p∼τcorr H*; (b) bounding δ(p) variation by |dδ/d ln q| times support width; (c) suppressing δ(|k−p|) near |k−p|→0 by a mass-gap ratio ρ; these are plausible but not established with explicit models for D0 and the TT projector weighting.
If wrong: If D0 is not sharply peaked or the TT projection weights regions where δ varies rapidly, the modulation may be smeared, phase-shifted, or reduced; then Eq. (38) and the ≲1% propagation claim fail.
- highEq. (37) — The mass-gap suppression estimate appears algebraically reversed. For a massive propagator of the form 1/(q^2+m_V^2), D_0(0)/D_0(q_*) is approximately q_*^2/m_V^2, not m_V^2/q_*^2. The text claims the small ratio rho=m_V^2/q_*^2 suppresses the |k-p|->0 region.
If wrong: The bound excluding rapidly oscillatory small-|k-p| contributions is invalid; the convolution factorization in Eq. (38) is not reliable.
- highEq. (38) — The statement that the factor of 2 from the two linear cross terms can be absorbed into the baseline normalization while keeping epsilon=epsilon_f is not mathematically justified. A multiplicative baseline rescaling does not remove a relative modulation-depth factor.
If wrong: The predicted modulation amplitude range Eq. (39) is at least quantitatively wrong, and the WTC detectability conclusions shift.
- highEq. (39) — The map from WTC parameters to epsilon in [0.04,0.18] and b in [1.7,2.8] depends on the unproven holographic mapping, the compressed propagator derivation, and the problematic convolution factorization.
If wrong: The paper's falsifiable WTC prediction band is unsupported, although the generic phenomenological template remains definable.
- highEq. (7) and transition to Eq. (9) — Replacement of F(k,Δη) by F(k)δ(Δη)+O(τcorr H*) is sketched via a ‘nascent delta function’ expansion; no explicit control norm or assumptions on F are given to justify the stated relative error scaling for the GW integral.
If wrong: If the δ-approximation does not control the double-time integral in Eq. (4) as claimed, then Eq. (10)–(15) (factorization theorem and error bounds) are unreliable, undermining the main phenomenological claim that the modulation survives to percent level.
- highEq. (9) — The replacement F(k,Delta eta) ≃ F(k) delta(eta-eta') + O(tau_corr H_*) is asserted via a nascent delta expansion, but the normalization and dimensional transition from F(k,0) delta_tau to F(k)=integral F dDelta eta are not fully derived. The estimate also requires control of variation of the rest of the integrand over tau_corr.
If wrong: The claimed O(tau_corr H_*) factorization error in the central theorem is not established; Eq. (15)'s accuracy bound would be unsupported.
- highEqs. (12)-(14) — The bound on the oscillatory Green-function contribution assumes |partial_eta(a^4 S)| less than about H_* a^4 S. For a first-order phase transition, the source envelope generally turns on/off on beta^{-1}, not necessarily H_*^{-1}. The boundary-term and smooth-envelope assumptions are not proved for the stated source models.
If wrong: The conclusion that the Green-function contribution is only O(H_*/k) and cannot contaminate the log-periodic modulation is not quantitatively established.
- highEqs. (34)-(36) — The replacement delta(p) ≈ delta(k) assumes the convolution support satisfies |ln(p/k)| less than about tau_corr H_*. The paper only states that D_0(p) is peaked at p~q_*; it does not derive that p tracks the external momentum k over the spectral range. In a three-dimensional convolution, p, |k-p|, and k are geometrically distinct.
If wrong: The convolution need not preserve the same log-periodic phase and amplitude as a function of external k; the claimed UETC ansatz Eq. (38) may fail.
- highSec. 3.2, statement before Eq. (9) — The paper states that G_k(eta,eta_2)a^2(eta_2) varies on timescale k^{-1} >> tau_corr for subhorizon modes. But for the characteristic phase-transition signal k is later taken to be of order beta/v_w, while tau_corr ~ v_w/beta, giving k tau_corr ~ 1 rather than k^{-1} >> tau_corr.
If wrong: The delta-function reduction of the unequal-time integral cannot be justified by slow variation of the Green function; the central short-correlation-time error estimate may fail.
- mediumEq. (10)–(11) — After inserting the δ-function, the resulting single-integral expression uses S(η1,η1) and defines F(k)=∫F(k,Δη)dΔη, but the mapping between the original Π0(k,η1,η2) and the new objects F0(k), S(η1,η1) is not derived; it is asserted as a decomposition.
If wrong: If Π0 does not admit this reduction, P0h(k,η) may inherit additional k-structure that can interfere with or mimic log-periodicity, weakening the separation claimed in Eq. (15)–(16).
- mediumEq. (12)–(14) — Integration-by-parts bound on Ĩ(k,η) assumes vanishing boundary terms ‘since S=0 outside the source epoch’ and uses |∂η1(a^4 S)|≲H* a^4 S without specifying conditions under which this inequality holds (e.g., sharp turn-on/off, differentiability, and whether H* should be conformal H).
If wrong: If boundary terms or envelope derivatives are not controlled, the estimate |Ĩ|/I0≲H*/k may fail, so the claim that Green’s-function oscillations cannot generate log-periodic structure at period ln b is not secured.
- mediumEq. (14), relation H_*/k ~ v_w/(beta/H_*) — The relation mixes conformal variables from Eq. (2), physical Hubble quantities used in beta/H_*, and observed/comoving wave numbers without a consistent scale-factor derivation.
If wrong: The numerical error estimates in Table 1 and the percent-level accuracy claim may be mis-scaled.
- mediumEq. (16)–(17) — Template written as exact multiplicative modulation and ‘exact sinusoid in ln f’ downstream, despite earlier O(τcorr H*) and later O(ε τcorr H*)+O(ρ) corrections.
If wrong: Data-analysis claims (perfect template correlation, exact sinusoidality) would be overstated; practical matched-filter performance could degrade or require additional nuisance parameters.
- mediumEq. (18) — The adopted sound-wave amplitude is given with (H_*/beta)^2. The standard formula usually contains a single power H_*/beta unless an additional finite-lifetime suppression is separately included and derived. No derivation of the squared dependence is provided.
If wrong: The baseline amplitude, LISA SNR estimates, and placement of the WTC band in the detectability plots could be quantitatively wrong.
- mediumEq. (32), Sec. 5.3 — Propagation of the potential modulation to the gauge-boson propagator via δm_V^2/m_V^2 = ε_f cos(2π ln(q/q_*)/ln b_0). The chain-rule argument substitutes the running VEV ⟨φ(q)⟩ ∝ q^{-Δ_φ} into the modulation, but the conversion of a cos(ln φ) modulation into a cos(ln q) modulation with the same period b_0 (up to an O(1) Δ_φ factor) is asserted rather than carried through with the dimensional factor made explicit.
If wrong: The identification b = b_0 at leading order in the WTC convolution (Eq. 38) would be modified by a Δ_φ-dependent factor, shifting the WTC b-band away from [1.7, 2.8]. The qualitative existence of the log-periodic modulation in the propagator would survive.
- mediumEq. (38) and statement 'absorbing [factor 2] into a redefined baseline normalization sets ε=εf at leading order' — The identification of the UETC modulation depth ε with εf is asserted. In a generic structure Π=Π0+εf Π1, the relative modulation depth depends on Π1/Π0; the factor 2 alone does not guarantee ε=εf.
If wrong: The predicted ε range in Eq. (39) could shift systematically (e.g., by an O(1) factor), affecting the claimed occupancy of the high-SNR region.
- mediumEq. (6) — DSI introduced as a purely k-dependent multiplicative modulation of the full UETC Π(k,η,η′) without showing that such factorization is compatible with general UETC structure for the relevant phase-transition sources.
If wrong: If the DSI modulation depends on (η,η′) as well as k, it will not factor out of the unequal-time integrals as used in Eq. (10)–(15), and the clean template Eq. (16) may be distorted or suppressed.
- mediumEqs. (30)-(31) — The mapping from a periodic 5D warp-factor perturbation A(y)->ky+delta A_0 sin(n_p k y) to the 4D technidilaton potential modulation with coefficient 4 delta A_0 and period ln b_0=kL/n_p is stated without derivation.
If wrong: The WTC/holographic origin of the predicted b and epsilon ranges is not established, though the phenomenological DSI template would remain intact.
- mediumSec. 5.1 / Fig. 4 — The numerical validation uses a separable UETC with C(k) independent of both times, for which factorization is algebraic. It explicitly does not test the non-separable convolution and Green-function corrections that determine the claimed percent-level error bound.
If wrong: The numerical experiment cannot validate the central approximation bounds; it only confirms the trivial separable case.
- mediumSec. 5.4, the claim 'the baseline propagator D_0(p) is sharply peaked at p ~ q_* with relative half-width Δp/p ~ τ_corr H_*' — The narrow-width property of the propagator in momentum space is asserted without derivation from the WTC dynamics. The statement appears to conflate temporal correlation length (τ_corr) with spatial momentum width; a rigorous justification of Δp/p is missing.
If wrong: If the peak is significantly broader, the bound on |δ(p)−δ(k)| could be larger than estimated, potentially invalidating the ≤1.4% error claim for the WTC convolution. The predicted ε, b band for the WTC model would be less reliable, although the model-independent DSI imprinting mechanism (factorization theorem) remains unaffected.
- lowEq. (13) and identification H*/k ≡ τcorr H* — The step relating H*/k to v_w/(β/H*) uses an asserted characteristic wavenumber k∼β a*/(v_w a0) but mixes quantities at different epochs and does not show the correct dimensionless ratio in conformal variables.
If wrong: Quantitative error bars in Table 1 and the claimed percent-level bounds could shift; the qualitative suppression with k≫H* may still hold, but the numerical mapping to β/H* could be off.
- lowEq. (13), Sec. 3.2 — Integration-by-parts bound |Ĩ(k,η)|/I_0 ≲ H_*/k uses |∂_{η_1}[a^4 S]| ≲ H_* a^4 S, which is plausible for a slowly varying envelope but not explicitly justified for the source envelope shape during a FOPT (where S can have sharp turn-on/turn-off features).
If wrong: The O(τ_corr H_*) bound on the Green's-function-induced log-periodic contamination could be weaker by an O(1) factor, slightly degrading the percent-level factorization accuracy but not invalidating the theorem.
- lowEq. (22)–(26) — Derivation of SNRosc scaling uses ⟨cos^2⟩=1/2 over Nperiods complete periods but does not address weighting by σ(f) and Ω0GW(f) variations across the band; Eq. (23)–(25) effectively assumes the integral decomposes into equal-contribution periods.
If wrong: The √Nperiods enhancement may be weaker/stronger depending on band-shape; this affects quantitative detectability contours but not the existence of an oscillatory matched-filter statistic.
- lowEq. (36), Sec. 5.4 — Numerical bound |δ(p) - δ(k)| ≲ 13 ε_f τ_corr H_* for b_0 ≥ 1.5 uses the absolute derivative bound 2π/ln b_0. The relative half-width ∆p/p ~ τ_corr H_* of the propagator peak is asserted but not derived from the propagator structure.
If wrong: Numerical coefficient in the convolution error bound could change by an O(1) factor; the qualitative O(ε_f τ_corr H_*) scaling would survive.
- lowTable 1 combined bound and Sec. 5.1 numerical validation — Numerical test uses a separable Π=F0(k)C(k)S(η1,η2)g(η1−η2), which algebraically factorizes C(k); it does not probe the non-separable corrections that are said to dominate the error budget (convolution, Green’s-function term).
If wrong: The presented ‘validation’ does not validate the stated percent-level error bounds for realistic non-separable sources; thus the empirical support for Table 1 is weaker than implied.
The paper presents a logically coherent framework in which discrete scale invariance in the source UETC imprints a log-periodic modulation on the SGWB. The core argument—ansatz to factorization to observable template—is structurally sound and internally consistent in its broad strokes. However, the mathematical rigor of the central approximation that enables factorization is compromised by a mismatch between the approximation's justification and the frequency range where it is applied: the condition k^{-1} >> tau_corr is violated at the peak of the GW spectrum, where k*tau_corr ~ 1. This means the error bound presented in the derivation does not fully capture the leading correction, and the claimed percent-level accuracy for the peak modes is not adequately supported. The UV-completion section adds further gaps, with the mapping from model parameters to observable parameters sketched rather than derived. These issues are not fatal—the qualitative prediction of a log-periodic modulation remains plausible—but they lower confidence in the quantitative predictions. Internal consistency is rated 3/5, and mathematical validity 3/5, reflecting a paper with a solid conceptual skeleton but insufficiently rigorous flesh on the bones of its most consequential derivations.
⚑Derivation Flags (28)
- highEq. (29) ⇒ Eq. (32) — Claim that a periodic modulation of the technidilaton potential induces δmV^2/mV^2=εf cos(2π ln(q/q*)/ln b0)+O(εf^2) and hence a multiplicative propagator modulation D=D0(1+δ(q)) is not derived; it assumes a specific dependence of masses/propagators on φ(q) and a direct ‘chain rule’ transfer of the cosine.
If wrong: If δ(q) is not of the asserted form/amplitude, the subsequent convolution does not yield Eq. (38), and the claimed mapping from WTC/holographic parameters to the SGWB modulation (ε,b) is unsupported.
- highEq. (31), Sec. 5.3 — Multiplicative log-periodic 4D potential V_4(φ) ≈ V_CW(φ)[1 + 4δA_0 cos(2π ln(φ/φ_0)/ln b_0)] derived from a sinusoidally perturbed AdS_5 warp factor is stated with citation to [17] (Goldberger-Wise) but no explicit derivation. The factor of 4, the cosine-in-log-φ structure, and the identification ln b_0 = kL/n_p are presented as results.
If wrong: The mapping from WTC/RS holographic parameters (δA_0, kL, n_p) to the observable (ε, b) band would be unsupported. The central falsifiable WTC prediction ε ∈ [0.04, 0.18], b ∈ [1.7, 2.8] would lose its microphysical anchor, reducing the UV-completion section to dimensional estimates.
- highEq. (32) — The derivation from a log-periodic technidilaton potential to a gauge-field propagator modulation D(q)=D_0(q)[1+delta(q)] is highly compressed. The asserted same-amplitude transfer by the chain rule, the treatment of the scaling dimension Delta_phi, and the positivity claim for D are not demonstrated.
If wrong: The starting point for the WTC UETC convolution is unsupported; Eq. (38) and the WTC prediction band Eq. (39) would not follow.
- highEq. (33)–(37) — Convolution factorization relies on (a) D0(p) being sharply peaked with width Δp/p∼τcorr H*; (b) bounding δ(p) variation by |dδ/d ln q| times support width; (c) suppressing δ(|k−p|) near |k−p|→0 by a mass-gap ratio ρ; these are plausible but not established with explicit models for D0 and the TT projector weighting.
If wrong: If D0 is not sharply peaked or the TT projection weights regions where δ varies rapidly, the modulation may be smeared, phase-shifted, or reduced; then Eq. (38) and the ≲1% propagation claim fail.
- highEq. (37) — The mass-gap suppression estimate appears algebraically reversed. For a massive propagator of the form 1/(q^2+m_V^2), D_0(0)/D_0(q_*) is approximately q_*^2/m_V^2, not m_V^2/q_*^2. The text claims the small ratio rho=m_V^2/q_*^2 suppresses the |k-p|->0 region.
If wrong: The bound excluding rapidly oscillatory small-|k-p| contributions is invalid; the convolution factorization in Eq. (38) is not reliable.
- highEq. (38) — The statement that the factor of 2 from the two linear cross terms can be absorbed into the baseline normalization while keeping epsilon=epsilon_f is not mathematically justified. A multiplicative baseline rescaling does not remove a relative modulation-depth factor.
If wrong: The predicted modulation amplitude range Eq. (39) is at least quantitatively wrong, and the WTC detectability conclusions shift.
- highEq. (39) — The map from WTC parameters to epsilon in [0.04,0.18] and b in [1.7,2.8] depends on the unproven holographic mapping, the compressed propagator derivation, and the problematic convolution factorization.
If wrong: The paper's falsifiable WTC prediction band is unsupported, although the generic phenomenological template remains definable.
- highEq. (7) and transition to Eq. (9) — Replacement of F(k,Δη) by F(k)δ(Δη)+O(τcorr H*) is sketched via a ‘nascent delta function’ expansion; no explicit control norm or assumptions on F are given to justify the stated relative error scaling for the GW integral.
If wrong: If the δ-approximation does not control the double-time integral in Eq. (4) as claimed, then Eq. (10)–(15) (factorization theorem and error bounds) are unreliable, undermining the main phenomenological claim that the modulation survives to percent level.
- highEq. (9) — The replacement F(k,Delta eta) ≃ F(k) delta(eta-eta') + O(tau_corr H_*) is asserted via a nascent delta expansion, but the normalization and dimensional transition from F(k,0) delta_tau to F(k)=integral F dDelta eta are not fully derived. The estimate also requires control of variation of the rest of the integrand over tau_corr.
If wrong: The claimed O(tau_corr H_*) factorization error in the central theorem is not established; Eq. (15)'s accuracy bound would be unsupported.
- highEqs. (12)-(14) — The bound on the oscillatory Green-function contribution assumes |partial_eta(a^4 S)| less than about H_* a^4 S. For a first-order phase transition, the source envelope generally turns on/off on beta^{-1}, not necessarily H_*^{-1}. The boundary-term and smooth-envelope assumptions are not proved for the stated source models.
If wrong: The conclusion that the Green-function contribution is only O(H_*/k) and cannot contaminate the log-periodic modulation is not quantitatively established.
- highEqs. (34)-(36) — The replacement delta(p) ≈ delta(k) assumes the convolution support satisfies |ln(p/k)| less than about tau_corr H_*. The paper only states that D_0(p) is peaked at p~q_*; it does not derive that p tracks the external momentum k over the spectral range. In a three-dimensional convolution, p, |k-p|, and k are geometrically distinct.
If wrong: The convolution need not preserve the same log-periodic phase and amplitude as a function of external k; the claimed UETC ansatz Eq. (38) may fail.
- highSec. 3.2, statement before Eq. (9) — The paper states that G_k(eta,eta_2)a^2(eta_2) varies on timescale k^{-1} >> tau_corr for subhorizon modes. But for the characteristic phase-transition signal k is later taken to be of order beta/v_w, while tau_corr ~ v_w/beta, giving k tau_corr ~ 1 rather than k^{-1} >> tau_corr.
If wrong: The delta-function reduction of the unequal-time integral cannot be justified by slow variation of the Green function; the central short-correlation-time error estimate may fail.
- mediumEq. (10)–(11) — After inserting the δ-function, the resulting single-integral expression uses S(η1,η1) and defines F(k)=∫F(k,Δη)dΔη, but the mapping between the original Π0(k,η1,η2) and the new objects F0(k), S(η1,η1) is not derived; it is asserted as a decomposition.
If wrong: If Π0 does not admit this reduction, P0h(k,η) may inherit additional k-structure that can interfere with or mimic log-periodicity, weakening the separation claimed in Eq. (15)–(16).
- mediumEq. (12)–(14) — Integration-by-parts bound on Ĩ(k,η) assumes vanishing boundary terms ‘since S=0 outside the source epoch’ and uses |∂η1(a^4 S)|≲H* a^4 S without specifying conditions under which this inequality holds (e.g., sharp turn-on/off, differentiability, and whether H* should be conformal H).
If wrong: If boundary terms or envelope derivatives are not controlled, the estimate |Ĩ|/I0≲H*/k may fail, so the claim that Green’s-function oscillations cannot generate log-periodic structure at period ln b is not secured.
- mediumEq. (14), relation H_*/k ~ v_w/(beta/H_*) — The relation mixes conformal variables from Eq. (2), physical Hubble quantities used in beta/H_*, and observed/comoving wave numbers without a consistent scale-factor derivation.
If wrong: The numerical error estimates in Table 1 and the percent-level accuracy claim may be mis-scaled.
- mediumEq. (16)–(17) — Template written as exact multiplicative modulation and ‘exact sinusoid in ln f’ downstream, despite earlier O(τcorr H*) and later O(ε τcorr H*)+O(ρ) corrections.
If wrong: Data-analysis claims (perfect template correlation, exact sinusoidality) would be overstated; practical matched-filter performance could degrade or require additional nuisance parameters.
- mediumEq. (18) — The adopted sound-wave amplitude is given with (H_*/beta)^2. The standard formula usually contains a single power H_*/beta unless an additional finite-lifetime suppression is separately included and derived. No derivation of the squared dependence is provided.
If wrong: The baseline amplitude, LISA SNR estimates, and placement of the WTC band in the detectability plots could be quantitatively wrong.
- mediumEq. (32), Sec. 5.3 — Propagation of the potential modulation to the gauge-boson propagator via δm_V^2/m_V^2 = ε_f cos(2π ln(q/q_*)/ln b_0). The chain-rule argument substitutes the running VEV ⟨φ(q)⟩ ∝ q^{-Δ_φ} into the modulation, but the conversion of a cos(ln φ) modulation into a cos(ln q) modulation with the same period b_0 (up to an O(1) Δ_φ factor) is asserted rather than carried through with the dimensional factor made explicit.
If wrong: The identification b = b_0 at leading order in the WTC convolution (Eq. 38) would be modified by a Δ_φ-dependent factor, shifting the WTC b-band away from [1.7, 2.8]. The qualitative existence of the log-periodic modulation in the propagator would survive.
- mediumEq. (38) and statement 'absorbing [factor 2] into a redefined baseline normalization sets ε=εf at leading order' — The identification of the UETC modulation depth ε with εf is asserted. In a generic structure Π=Π0+εf Π1, the relative modulation depth depends on Π1/Π0; the factor 2 alone does not guarantee ε=εf.
If wrong: The predicted ε range in Eq. (39) could shift systematically (e.g., by an O(1) factor), affecting the claimed occupancy of the high-SNR region.
- mediumEq. (6) — DSI introduced as a purely k-dependent multiplicative modulation of the full UETC Π(k,η,η′) without showing that such factorization is compatible with general UETC structure for the relevant phase-transition sources.
If wrong: If the DSI modulation depends on (η,η′) as well as k, it will not factor out of the unequal-time integrals as used in Eq. (10)–(15), and the clean template Eq. (16) may be distorted or suppressed.
- mediumEqs. (30)-(31) — The mapping from a periodic 5D warp-factor perturbation A(y)->ky+delta A_0 sin(n_p k y) to the 4D technidilaton potential modulation with coefficient 4 delta A_0 and period ln b_0=kL/n_p is stated without derivation.
If wrong: The WTC/holographic origin of the predicted b and epsilon ranges is not established, though the phenomenological DSI template would remain intact.
- mediumSec. 5.1 / Fig. 4 — The numerical validation uses a separable UETC with C(k) independent of both times, for which factorization is algebraic. It explicitly does not test the non-separable convolution and Green-function corrections that determine the claimed percent-level error bound.
If wrong: The numerical experiment cannot validate the central approximation bounds; it only confirms the trivial separable case.
- mediumSec. 5.4, the claim 'the baseline propagator D_0(p) is sharply peaked at p ~ q_* with relative half-width Δp/p ~ τ_corr H_*' — The narrow-width property of the propagator in momentum space is asserted without derivation from the WTC dynamics. The statement appears to conflate temporal correlation length (τ_corr) with spatial momentum width; a rigorous justification of Δp/p is missing.
If wrong: If the peak is significantly broader, the bound on |δ(p)−δ(k)| could be larger than estimated, potentially invalidating the ≤1.4% error claim for the WTC convolution. The predicted ε, b band for the WTC model would be less reliable, although the model-independent DSI imprinting mechanism (factorization theorem) remains unaffected.
- lowEq. (13) and identification H*/k ≡ τcorr H* — The step relating H*/k to v_w/(β/H*) uses an asserted characteristic wavenumber k∼β a*/(v_w a0) but mixes quantities at different epochs and does not show the correct dimensionless ratio in conformal variables.
If wrong: Quantitative error bars in Table 1 and the claimed percent-level bounds could shift; the qualitative suppression with k≫H* may still hold, but the numerical mapping to β/H* could be off.
- lowEq. (13), Sec. 3.2 — Integration-by-parts bound |Ĩ(k,η)|/I_0 ≲ H_*/k uses |∂_{η_1}[a^4 S]| ≲ H_* a^4 S, which is plausible for a slowly varying envelope but not explicitly justified for the source envelope shape during a FOPT (where S can have sharp turn-on/turn-off features).
If wrong: The O(τ_corr H_*) bound on the Green's-function-induced log-periodic contamination could be weaker by an O(1) factor, slightly degrading the percent-level factorization accuracy but not invalidating the theorem.
- lowEq. (22)–(26) — Derivation of SNRosc scaling uses ⟨cos^2⟩=1/2 over Nperiods complete periods but does not address weighting by σ(f) and Ω0GW(f) variations across the band; Eq. (23)–(25) effectively assumes the integral decomposes into equal-contribution periods.
If wrong: The √Nperiods enhancement may be weaker/stronger depending on band-shape; this affects quantitative detectability contours but not the existence of an oscillatory matched-filter statistic.
- lowEq. (36), Sec. 5.4 — Numerical bound |δ(p) - δ(k)| ≲ 13 ε_f τ_corr H_* for b_0 ≥ 1.5 uses the absolute derivative bound 2π/ln b_0. The relative half-width ∆p/p ~ τ_corr H_* of the propagator peak is asserted but not derived from the propagator structure.
If wrong: Numerical coefficient in the convolution error bound could change by an O(1) factor; the qualitative O(ε_f τ_corr H_*) scaling would survive.
- lowTable 1 combined bound and Sec. 5.1 numerical validation — Numerical test uses a separable Π=F0(k)C(k)S(η1,η2)g(η1−η2), which algebraically factorizes C(k); it does not probe the non-separable corrections that are said to dominate the error budget (convolution, Green’s-function term).
If wrong: The presented ‘validation’ does not validate the stated percent-level error bounds for realistic non-separable sources; thus the empirical support for Table 1 is weaker than implied.
The paper is coherent as a conditional phenomenological statement: if the UETC has the multiplicative k-only DSI form in Eq. (6), then the observable spectrum inherits that modulation. That part of the argument is not internally contradictory and is mathematically simple.
The broader submission, however, makes stronger claims: that the short-correlation approximation controls realistic UETCs at the percent level and that walking technicolor predicts a specific observable ε,b band. Those claims require consistent identification of the modulation across potential, propagator, UETC, and observable levels. The provided derivations do not close that chain, and several quantitative statements conflict with the assumptions used to justify them. The strongest high-score argument about the paper’s transparent hierarchy is acknowledged, but it does not cure the central parameter-definition and approximation-regime inconsistencies. Therefore the internal-consistency score is 2/5.
⚑Derivation Flags (28)
- highEq. (29) ⇒ Eq. (32) — Claim that a periodic modulation of the technidilaton potential induces δmV^2/mV^2=εf cos(2π ln(q/q*)/ln b0)+O(εf^2) and hence a multiplicative propagator modulation D=D0(1+δ(q)) is not derived; it assumes a specific dependence of masses/propagators on φ(q) and a direct ‘chain rule’ transfer of the cosine.
If wrong: If δ(q) is not of the asserted form/amplitude, the subsequent convolution does not yield Eq. (38), and the claimed mapping from WTC/holographic parameters to the SGWB modulation (ε,b) is unsupported.
- highEq. (31), Sec. 5.3 — Multiplicative log-periodic 4D potential V_4(φ) ≈ V_CW(φ)[1 + 4δA_0 cos(2π ln(φ/φ_0)/ln b_0)] derived from a sinusoidally perturbed AdS_5 warp factor is stated with citation to [17] (Goldberger-Wise) but no explicit derivation. The factor of 4, the cosine-in-log-φ structure, and the identification ln b_0 = kL/n_p are presented as results.
If wrong: The mapping from WTC/RS holographic parameters (δA_0, kL, n_p) to the observable (ε, b) band would be unsupported. The central falsifiable WTC prediction ε ∈ [0.04, 0.18], b ∈ [1.7, 2.8] would lose its microphysical anchor, reducing the UV-completion section to dimensional estimates.
- highEq. (32) — The derivation from a log-periodic technidilaton potential to a gauge-field propagator modulation D(q)=D_0(q)[1+delta(q)] is highly compressed. The asserted same-amplitude transfer by the chain rule, the treatment of the scaling dimension Delta_phi, and the positivity claim for D are not demonstrated.
If wrong: The starting point for the WTC UETC convolution is unsupported; Eq. (38) and the WTC prediction band Eq. (39) would not follow.
- highEq. (33)–(37) — Convolution factorization relies on (a) D0(p) being sharply peaked with width Δp/p∼τcorr H*; (b) bounding δ(p) variation by |dδ/d ln q| times support width; (c) suppressing δ(|k−p|) near |k−p|→0 by a mass-gap ratio ρ; these are plausible but not established with explicit models for D0 and the TT projector weighting.
If wrong: If D0 is not sharply peaked or the TT projection weights regions where δ varies rapidly, the modulation may be smeared, phase-shifted, or reduced; then Eq. (38) and the ≲1% propagation claim fail.
- highEq. (37) — The mass-gap suppression estimate appears algebraically reversed. For a massive propagator of the form 1/(q^2+m_V^2), D_0(0)/D_0(q_*) is approximately q_*^2/m_V^2, not m_V^2/q_*^2. The text claims the small ratio rho=m_V^2/q_*^2 suppresses the |k-p|->0 region.
If wrong: The bound excluding rapidly oscillatory small-|k-p| contributions is invalid; the convolution factorization in Eq. (38) is not reliable.
- highEq. (38) — The statement that the factor of 2 from the two linear cross terms can be absorbed into the baseline normalization while keeping epsilon=epsilon_f is not mathematically justified. A multiplicative baseline rescaling does not remove a relative modulation-depth factor.
If wrong: The predicted modulation amplitude range Eq. (39) is at least quantitatively wrong, and the WTC detectability conclusions shift.
- highEq. (39) — The map from WTC parameters to epsilon in [0.04,0.18] and b in [1.7,2.8] depends on the unproven holographic mapping, the compressed propagator derivation, and the problematic convolution factorization.
If wrong: The paper's falsifiable WTC prediction band is unsupported, although the generic phenomenological template remains definable.
- highEq. (7) and transition to Eq. (9) — Replacement of F(k,Δη) by F(k)δ(Δη)+O(τcorr H*) is sketched via a ‘nascent delta function’ expansion; no explicit control norm or assumptions on F are given to justify the stated relative error scaling for the GW integral.
If wrong: If the δ-approximation does not control the double-time integral in Eq. (4) as claimed, then Eq. (10)–(15) (factorization theorem and error bounds) are unreliable, undermining the main phenomenological claim that the modulation survives to percent level.
- highEq. (9) — The replacement F(k,Delta eta) ≃ F(k) delta(eta-eta') + O(tau_corr H_*) is asserted via a nascent delta expansion, but the normalization and dimensional transition from F(k,0) delta_tau to F(k)=integral F dDelta eta are not fully derived. The estimate also requires control of variation of the rest of the integrand over tau_corr.
If wrong: The claimed O(tau_corr H_*) factorization error in the central theorem is not established; Eq. (15)'s accuracy bound would be unsupported.
- highEqs. (12)-(14) — The bound on the oscillatory Green-function contribution assumes |partial_eta(a^4 S)| less than about H_* a^4 S. For a first-order phase transition, the source envelope generally turns on/off on beta^{-1}, not necessarily H_*^{-1}. The boundary-term and smooth-envelope assumptions are not proved for the stated source models.
If wrong: The conclusion that the Green-function contribution is only O(H_*/k) and cannot contaminate the log-periodic modulation is not quantitatively established.
- highEqs. (34)-(36) — The replacement delta(p) ≈ delta(k) assumes the convolution support satisfies |ln(p/k)| less than about tau_corr H_*. The paper only states that D_0(p) is peaked at p~q_*; it does not derive that p tracks the external momentum k over the spectral range. In a three-dimensional convolution, p, |k-p|, and k are geometrically distinct.
If wrong: The convolution need not preserve the same log-periodic phase and amplitude as a function of external k; the claimed UETC ansatz Eq. (38) may fail.
- highSec. 3.2, statement before Eq. (9) — The paper states that G_k(eta,eta_2)a^2(eta_2) varies on timescale k^{-1} >> tau_corr for subhorizon modes. But for the characteristic phase-transition signal k is later taken to be of order beta/v_w, while tau_corr ~ v_w/beta, giving k tau_corr ~ 1 rather than k^{-1} >> tau_corr.
If wrong: The delta-function reduction of the unequal-time integral cannot be justified by slow variation of the Green function; the central short-correlation-time error estimate may fail.
- mediumEq. (10)–(11) — After inserting the δ-function, the resulting single-integral expression uses S(η1,η1) and defines F(k)=∫F(k,Δη)dΔη, but the mapping between the original Π0(k,η1,η2) and the new objects F0(k), S(η1,η1) is not derived; it is asserted as a decomposition.
If wrong: If Π0 does not admit this reduction, P0h(k,η) may inherit additional k-structure that can interfere with or mimic log-periodicity, weakening the separation claimed in Eq. (15)–(16).
- mediumEq. (12)–(14) — Integration-by-parts bound on Ĩ(k,η) assumes vanishing boundary terms ‘since S=0 outside the source epoch’ and uses |∂η1(a^4 S)|≲H* a^4 S without specifying conditions under which this inequality holds (e.g., sharp turn-on/off, differentiability, and whether H* should be conformal H).
If wrong: If boundary terms or envelope derivatives are not controlled, the estimate |Ĩ|/I0≲H*/k may fail, so the claim that Green’s-function oscillations cannot generate log-periodic structure at period ln b is not secured.
- mediumEq. (14), relation H_*/k ~ v_w/(beta/H_*) — The relation mixes conformal variables from Eq. (2), physical Hubble quantities used in beta/H_*, and observed/comoving wave numbers without a consistent scale-factor derivation.
If wrong: The numerical error estimates in Table 1 and the percent-level accuracy claim may be mis-scaled.
- mediumEq. (16)–(17) — Template written as exact multiplicative modulation and ‘exact sinusoid in ln f’ downstream, despite earlier O(τcorr H*) and later O(ε τcorr H*)+O(ρ) corrections.
If wrong: Data-analysis claims (perfect template correlation, exact sinusoidality) would be overstated; practical matched-filter performance could degrade or require additional nuisance parameters.
- mediumEq. (18) — The adopted sound-wave amplitude is given with (H_*/beta)^2. The standard formula usually contains a single power H_*/beta unless an additional finite-lifetime suppression is separately included and derived. No derivation of the squared dependence is provided.
If wrong: The baseline amplitude, LISA SNR estimates, and placement of the WTC band in the detectability plots could be quantitatively wrong.
- mediumEq. (32), Sec. 5.3 — Propagation of the potential modulation to the gauge-boson propagator via δm_V^2/m_V^2 = ε_f cos(2π ln(q/q_*)/ln b_0). The chain-rule argument substitutes the running VEV ⟨φ(q)⟩ ∝ q^{-Δ_φ} into the modulation, but the conversion of a cos(ln φ) modulation into a cos(ln q) modulation with the same period b_0 (up to an O(1) Δ_φ factor) is asserted rather than carried through with the dimensional factor made explicit.
If wrong: The identification b = b_0 at leading order in the WTC convolution (Eq. 38) would be modified by a Δ_φ-dependent factor, shifting the WTC b-band away from [1.7, 2.8]. The qualitative existence of the log-periodic modulation in the propagator would survive.
- mediumEq. (38) and statement 'absorbing [factor 2] into a redefined baseline normalization sets ε=εf at leading order' — The identification of the UETC modulation depth ε with εf is asserted. In a generic structure Π=Π0+εf Π1, the relative modulation depth depends on Π1/Π0; the factor 2 alone does not guarantee ε=εf.
If wrong: The predicted ε range in Eq. (39) could shift systematically (e.g., by an O(1) factor), affecting the claimed occupancy of the high-SNR region.
- mediumEq. (6) — DSI introduced as a purely k-dependent multiplicative modulation of the full UETC Π(k,η,η′) without showing that such factorization is compatible with general UETC structure for the relevant phase-transition sources.
If wrong: If the DSI modulation depends on (η,η′) as well as k, it will not factor out of the unequal-time integrals as used in Eq. (10)–(15), and the clean template Eq. (16) may be distorted or suppressed.
- mediumEqs. (30)-(31) — The mapping from a periodic 5D warp-factor perturbation A(y)->ky+delta A_0 sin(n_p k y) to the 4D technidilaton potential modulation with coefficient 4 delta A_0 and period ln b_0=kL/n_p is stated without derivation.
If wrong: The WTC/holographic origin of the predicted b and epsilon ranges is not established, though the phenomenological DSI template would remain intact.
- mediumSec. 5.1 / Fig. 4 — The numerical validation uses a separable UETC with C(k) independent of both times, for which factorization is algebraic. It explicitly does not test the non-separable convolution and Green-function corrections that determine the claimed percent-level error bound.
If wrong: The numerical experiment cannot validate the central approximation bounds; it only confirms the trivial separable case.
- mediumSec. 5.4, the claim 'the baseline propagator D_0(p) is sharply peaked at p ~ q_* with relative half-width Δp/p ~ τ_corr H_*' — The narrow-width property of the propagator in momentum space is asserted without derivation from the WTC dynamics. The statement appears to conflate temporal correlation length (τ_corr) with spatial momentum width; a rigorous justification of Δp/p is missing.
If wrong: If the peak is significantly broader, the bound on |δ(p)−δ(k)| could be larger than estimated, potentially invalidating the ≤1.4% error claim for the WTC convolution. The predicted ε, b band for the WTC model would be less reliable, although the model-independent DSI imprinting mechanism (factorization theorem) remains unaffected.
- lowEq. (13) and identification H*/k ≡ τcorr H* — The step relating H*/k to v_w/(β/H*) uses an asserted characteristic wavenumber k∼β a*/(v_w a0) but mixes quantities at different epochs and does not show the correct dimensionless ratio in conformal variables.
If wrong: Quantitative error bars in Table 1 and the claimed percent-level bounds could shift; the qualitative suppression with k≫H* may still hold, but the numerical mapping to β/H* could be off.
- lowEq. (13), Sec. 3.2 — Integration-by-parts bound |Ĩ(k,η)|/I_0 ≲ H_*/k uses |∂_{η_1}[a^4 S]| ≲ H_* a^4 S, which is plausible for a slowly varying envelope but not explicitly justified for the source envelope shape during a FOPT (where S can have sharp turn-on/turn-off features).
If wrong: The O(τ_corr H_*) bound on the Green's-function-induced log-periodic contamination could be weaker by an O(1) factor, slightly degrading the percent-level factorization accuracy but not invalidating the theorem.
- lowEq. (22)–(26) — Derivation of SNRosc scaling uses ⟨cos^2⟩=1/2 over Nperiods complete periods but does not address weighting by σ(f) and Ω0GW(f) variations across the band; Eq. (23)–(25) effectively assumes the integral decomposes into equal-contribution periods.
If wrong: The √Nperiods enhancement may be weaker/stronger depending on band-shape; this affects quantitative detectability contours but not the existence of an oscillatory matched-filter statistic.
- lowEq. (36), Sec. 5.4 — Numerical bound |δ(p) - δ(k)| ≲ 13 ε_f τ_corr H_* for b_0 ≥ 1.5 uses the absolute derivative bound 2π/ln b_0. The relative half-width ∆p/p ~ τ_corr H_* of the propagator peak is asserted but not derived from the propagator structure.
If wrong: Numerical coefficient in the convolution error bound could change by an O(1) factor; the qualitative O(ε_f τ_corr H_*) scaling would survive.
- lowTable 1 combined bound and Sec. 5.1 numerical validation — Numerical test uses a separable Π=F0(k)C(k)S(η1,η2)g(η1−η2), which algebraically factorizes C(k); it does not probe the non-separable corrections that are said to dominate the error budget (convolution, Green’s-function term).
If wrong: The presented ‘validation’ does not validate the stated percent-level error bounds for realistic non-separable sources; thus the empirical support for Table 1 is weaker than implied.
Factorization theorem: under the short-correlation-time approximation the DSI modulation C(k) factors multiplicatively out of the unequal-time integrals to percent-level accuracy.
Observable SGWB energy-density spectrum: the predicted universal log-periodic template (sinusoid in ln f) superimposed on the smooth baseline spectrum.
Matched-filter detectability scaling: the oscillatory component SNR scales with the modulation amplitude ε, the baseline per-log-period SNR, and the number of complete log-periods in the detector band.
The SGWB from a first-order cosmological phase transition with DSI in the source UETC will exhibit a multiplicative log-periodic modulation described by Ω_GW(f)=Ω0_GW(f)[1+ε cos(2π ln(f/f_*)/ln b + φ0)].
Falsifiable if: A sufficiently sensitive search (e.g. matched-filter search on LISA data) that rules out log-periodic modulations with amplitude ε at the predicted frequencies and periods (for given b and φ0) at the claimed SNR thresholds would falsify the claim for those parameter values.
Walking technicolor provides a viable UV completion realizing the DSI modulation, predicting ε in [0.04,0.18] and b in [1.7,2.8], which lie in the high-SNR region for LISA.
Falsifiable if: A null result by LISA (or other sufficiently sensitive SGWB searches) that excludes the predicted modulation amplitudes across the b∈[1.7,2.8] band at the forecasted SNR would rule out the WTC parameter band as the source of an observable DSI-modulated SGWB.
Under the short-correlation-time approximation (τ_corr H_* ≪ 1, e.g. β/H_* ≳ 10) the DSI modulation factorizes through the unequal-time integrals so that the modulation survives to the observable spectrum with percent-level relative error.
Falsifiable if: High-precision numerical evaluation of the full double-time integral (without the short-correlation-time approximation) for realistic UETCs that finds deviations from multiplicative factorization exceeding the stated percent-level bounds for β/H_* ≳ 10 would falsify the factorization claim in that regime.
Share this Review
Post your AI review credential to social media, or copy the link to share anywhere.
theoryofeverything.ai/review-profile/paper/2b3f4e99-df22-4c37-b7f2-83c659258b30Embed a live badge
Paste into a GitHub README, preprint page, or personal site. It reads the current review, so it updates on its own if the score changes or the work is withdrawn.
[](https://theoryofeverything.ai/review-profile/paper/2b3f4e99-df22-4c37-b7f2-83c659258b30?utm_source=toeshare_badge)How to use it: Markdown goes in a GitHub README or any Markdown page. HTML goes in your own web page, wherever you want the badge to appear. Both show the live badge and link back to this review profile.
Zenodo descriptions strip images; there, link the text TOE-Share 3.3 / 5 to your review profile instead. Show us you shared your work on social media — tag @TOE_Share in your post, or reply to our newsletter with the link. Either works — we’ll add one free review credit to your account.
Share by Email
Email clients cannot render the full review profile page. We send a branded HTML summary plus a link to the live credential.
Sign in as the submission owner to send a branded HTML email from TOE-Share. Anyone can still copy the text or open their email app.
This review was conducted by TOE-Share's multi-agent AI specialist pipeline. Each dimension is independently evaluated by specialist agents (Math/Logic, Sources/Evidence, Science/Novelty), then synthesized by a coordinator agent. This methodology is aligned with the multi-model AI feedback approach validated in Thakkar et al., Nature Machine Intelligence 2026.
TOE-Share — theoryofeverything.ai