paper Review Profile
Log-Periodic Signatures from Discrete Scale Invariance in the Stochastic Gravitational-Wave Background: Walking Technicolor as a Candidate Ultraviolet Completion
Discrete scale invariance (DSI) in the anisotropic stress of a first-order cosmological phase transition imprints a multiplicative log-periodic modulation on the stochastic gravitational-wave background, and under a short-correlation-time factorization theorem this modulation propagates to the observable spectrum at the percent level. As a concrete UV completion, walking technicolor can produce the required DSI and predicts ε∈[0.04,0.18], b∈[1.7,2.8], placing the signal in the high-SNR region for LISA and enabling enhanced matched-filter detectability.
Read the Full BreakdownFull 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-mp7qc1gm
The conditional phenomenological core is mostly internally coherent: if the UETC has the multiplicative form Π = Π_0 C(k) as in Eq. (6), then the linear Green-function integral Eq. (4) indeed carries C(k) into P_h and hence into Omega_GW. The matched-filter notation is also clarified by defining SNR_baseline as a per-log-period quantity. However, there are moderate internal consistency problems. The paper alternates between describing the factorization as an approximate short-correlation-time theorem and as an exact algebraic consequence of separability; these are not equivalent scopes. It also presents percent-level accuracy in the abstract and conclusions while Table 1 gives 15% at beta/H_* = 10 and Section 5.3 gives about 10% WTC mass-gap corrections. Most importantly, the WTC parameter band is presented as sharp in the abstract and figures, but Section 5.4 admits O(1) model-dependent rescalings and an order-of-magnitude status. These inconsistencies undermine the precision of the WTC/detectability conclusions, though they do not destroy the narrower conditional statement from Eq. (6) to Eq. (16).
Many equations are standard (Eqs. (2)–(5) are conventional, and the step from Eq. (22) to (25) is mathematically correct under the stated averaging assumptions). However, the core ‘factorization theorem’ depends on compressed and partially self-undermined approximations: Eq. (9) is asserted without a rigorous error bound propagated through Eq. (4), and the Green’s-function evaluation at η2=η1 is used even where kτ_corr∼1 (peak), where the paper notes possible O(1) corrections. The integration-by-parts bound (Eq. (13)) is mathematically fine given smoothness assumptions, but the paper itself notes those assumptions can fail for β-timescale features, which removes smallness. The qualitative ‘spectral separation’ argument is plausible but not a proof that no log-periodic component is generated by the remaining integrals. In the UV section, Eqs. (32) and (33)–(38) involve heuristic scaling arguments and unspecified O(1) factors (c_geom), so the quantitative WTC ranges in Eq. (39) are not derived. Because the unverified steps are central to the claimed percent-level propagation into Ω_GW, mathematical validity cannot exceed 3 under the rubric.
The work makes a clear observational prediction: a sinusoidal residual in ln f superimposed on the SGWB spectrum, with measurable parameters (ε, b, φ0), and it identifies a concrete detection channel via LISA matched filtering. This is a good falsifiability structure because it proposes a feature not generic to standard smooth broken-power-law SGWB templates, so the theory can be distinguished observationally rather than merely accommodated. The paper also states a candidate parameter region ε∈[0.04,0.18], b∈[1.7,2.8], which is substantially more testable than a purely qualitative claim. The main limitation is that the strongest quantitative detectability claims depend on assumed baseline SNR values and on a UV mapping that the paper itself admits is only order-of-magnitude and model-dependent. The paper does not give a crisp statement of what specific observed non-detection contour in (ε,b) would falsify the WTC completion versus merely constrain it, and some claimed error control in the abstract is stronger than warranted by the body. Still, the phenomenological template itself is clearly testable with near-term planned instrumentation, so the paper deserves a high but not top falsifiability score.
The paper is generally organized well, with clear sectioning, a helpful hierarchy-of-claims table, and repeated efforts to distinguish derived phenomenology from conjectural UV completion. A scientifically literate reader can follow the intended narrative: DSI ansatz at the source level, factorization claim, observational template, then candidate realization in WTC. That communicative structure is a real strength. However, clarity is reduced by several issues. First, the paper oscillates between exact and approximate statements about factorization, and the reader must work to disentangle what is proven algebraically for separable toy UETCs from what is only argued heuristically for physical sources. Second, the abstract materially overclaims the status of the results relative to the caveats later admitted in Sec. 5.4, which weakens communicative trust. Third, some passages are dense with caveats, nested assumptions, and competing error estimates, making the central claim harder to parse than necessary. Because the overclaim is material and because important limitations are only clarified later, clarity cannot be scored above 3.
The application of discrete scale invariance to the anisotropic stress UETC of a first-order phase transition, with a derived factorization theorem propagating the log-periodic modulation to the observable SGWB at the percent level, appears to be a genuinely novel synthesis. Log-periodic GW signatures have been discussed in inflationary/beyond-Einstein contexts (Calcagni & Kuroyanagi cited), but the route via DSI in the UETC of an FOPT, and the explicit WTC realization with a holographic warp-factor origin (Eq. 30–31), is a new combination. The author appropriately situates the work relative to prior literature on DSI in condensed matter and finance and on FOPT GW signatures. Novelty is somewhat moderated by the fact that the WTC UV-completion chain is acknowledged as 'plausibility argument' rather than first-principles derivation, so the truly novel contribution is concentrated in the model-independent factorization plus template, rather than in the UV-completion itself.
The paper is reasonably complete on its core phenomenological claim but only partially complete on its claimed UV completion. On the positive side, the main objects in the SGWB calculation are mostly defined before use, assumptions are repeatedly stated, and the author does a good job separating derived phenomenology from conjectural model-building. The factorization argument is laid out in a followable sequence, approximation regimes are discussed, and limitations are explicitly acknowledged in Secs. 5.4 and 6. The paper also addresses its principal observable goal by providing the template, detectability scaling, and a falsifiability discussion. The main incompleteness lies in the WTC-specific chain from potential modulation to a concrete predicted parameter band. Several key steps are only sketched: the periodic technidilaton potential is motivated but not microscopically derived; the transfer from V(φ) to D(q) is heuristic; the convolution introduces an uncomputed geometric factor c_geom; and the final mapping to ε ∈ [0.04, 0.18], b ∈ [1.7, 2.8] retains acknowledged O(1) uncertainty. The paper does explicitly admit these weaknesses, which helps internal honesty, but they still limit completeness because the title and abstract present WTC as a concrete candidate UV completion with quantitative predictions. There is also a notable internal tension in Sec. 3.2: the discussion alternates between a short-correlation-time/δ-function approximation, a broader separability argument that makes factorization algebraic, and caveats that the peak-scale slow-variation approximation can be marginal with O(1) residuals. The author partially resolves this by distinguishing separable from non-separable corrections, but the treatment of boundary cases is not fully unified. Overall, the central phenomenology is followable, but the UV completion and quantitative prediction layer remains structurally incomplete.
This paper presents a genuinely novel and scientifically interesting proposal: that discrete scale invariance (DSI) in the anisotropic stress tensor of a first-order cosmological phase transition could imprint log-periodic modulations on the stochastic gravitational-wave background, with walking technicolor (WTC) as a candidate ultraviolet completion. The phenomenological core is mathematically sound and internally coherent, providing a concrete, falsifiable template for LISA observations. However, significant mathematical gaps emerge in both the factorization error control and the WTC ultraviolet completion chain. The factorization theorem itself is correct as an algebraic statement: if the unequal-time correlator (UETC) has the exact separable form Π(k,η,η') = S(η,η')F(k)g(Δη) with a purely k-dependent modulation C(k), then this modulation factors through the tensor power spectrum integral exactly. The numerical validation in Figure 1 confirms this algebraic property at machine precision. However, the extension to realistic non-separable phase transition sources involves approximations that are not fully controlled. The specialists identify several critical gaps: Equation (9)'s δ-function approximation lacks rigorous error bounds, particularly near the spectral peak where kτ_corr ~ 1 makes the slow-variation assumption marginal with acknowledged O(1) corrections. The integration-by-parts bound (Equation 13) assumes smooth source envelopes, but realistic transition boundaries may vary on β^(-1) timescales, weakening the bound to O(1). More seriously, the WTC ultraviolet completion contains multiple unverified steps that undermine the quantitative predictions. Equation (31) presents the holographic mapping ε_f = 4δA_0 without derivation from the cited Goldberger-Wise reference. Equation (32)'s chain-rule transfer from potential modulation to propagator modulation is sketched but not rigorously established. Most critically, the convolution calculation in Equations (35)-(38) that produces the final WTC parameter band relies on replacing internal momenta p and |k-p| with external momentum k inside log-periodic functions - a step that is not generally valid for broad convolution integrals. The geometric factor c_geom remains uncomputed, introducing unconstrained O(1) uncertainty into the predicted ranges ε ∈ [0.04, 0.18] and b ∈ [1.7, 2.8]. Despite these mathematical limitations, the work makes valuable contributions. The phenomenological template (Equation 16) is well-defined and provides a specific, testable signature that could distinguish DSI-induced signals from standard SGWB features. The matched-filter detectability analysis gives concrete SNR scaling and places the predicted signal in LISA's high-sensitivity region. The authors demonstrate commendable epistemic discipline by clearly separating derived results from conjectural elements (Table 2) and acknowledging the order-of-magnitude status of their WTC predictions in Section 5.4, though this honest assessment somewhat contradicts the sharper claims in the abstract.
Strengths
- +Proposes a genuinely novel observable signature: coherent log-periodic oscillations in the SGWB spectrum that could be distinguished from standard phase transition features
- +Provides a concrete, falsifiable template Ω_GW(f) = Ω^0_GW(f)[1 + ε cos(2π ln(f/f*)/ln b + φ_0)] with explicit parameter ranges testable by LISA
- +Demonstrates rigorous algebraic factorization for separable UETCs with machine-precision numerical validation
- +Exhibits exceptional epistemic discipline by clearly distinguishing derived phenomenological results from conjectural UV completion elements (Table 2)
- +Connects to a specific BSM framework (walking technicolor) that already predicts LISA-detectable gravitational waves
- +Provides quantitative error bounds and approximation structure (Table 1) with transparent discussion of limitations
Areas for Improvement
- -Derive rigorous error bounds for the δ-function approximation (Eq. 9), particularly addressing the acknowledged marginal regime where kτ_corr ~ 1 at the spectral peak
- -Provide the missing holographic derivation of Eq. (31) linking warp-factor perturbations to technidilaton potential modulation, or acknowledge this as a phenomenological ansatz
- -Rigorously establish the chain-rule transfer from V(φ) modulation to propagator modulation D(q) in Eq. (32), including proper treatment of scaling dimensions and phase factors
- -Justify or correct the convolution approximation in Eqs. (35)-(38) where internal momenta are replaced by external momentum inside log-periodic functions
- -Compute the geometric factor c_geom and other O(1) model-dependent coefficients to support the quoted numerical prediction ranges
- -Reconcile the tension between 'percent-level' claims in the abstract and the 15% error bounds for β/H* = 10 shown in Table 1
- -Address potential confusion with detector systematics or calibration-induced oscillations that could mimic log-periodic signals in matched-filter searches
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 jillfarleyrankin@gmail.com 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 . . . . . . . . . . . . . . . . . . . . . . . . . . .9 5 Ultraviolet Completion: Walking Technicolor11 5.1 Phase-transition parameter space . . . . . . . . . . . . . . . . . . . . . . . 13 5.2 Engineering discrete scale invariance . . . . . . . . . . . . . . . . . . . . . 13 5.3 Convolution for the UETC . . . . . . . . . . . . . . . . . . . . . . . . . . . 14 5.4 WTC predictions . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15 6 Discussion16 7 Conclusions18 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 . The Green’s functionG k (η,η 2 )a 2 (η 2 ) varies on timescalek −1 . For sub-horizon modes well above the peak (k ≫ β/v w ),k −1 ≪ τ corr is automatic; at the characteristic peak (k ∼ β/v w ),kτ corr ∼1 and the slow-variation approximation is marginal, withO(1) residual corrections that we do not compute explicitly. In either case we evaluate G k (η,η 2 ) at η 2 = η 1 , giving F (k,η− η ′ )≃ F (k)δ(η− η ′ ) +O(τ corr H ∗ ),(9) whereF(k)≡ R F(k,∆η)d(∆η). The relativeO(τ corr H ∗ ) error is bounded in Table 1; we emphasize, however, that the factorization theorem proved below rests ultimately on separability of the UETC in (k,∆η), not on theδ-function limit per se — for any separable Π =S(η,η ′ )F(k)g(∆η) the modulationC(k) factors out of the time integral algebraically, regardless of how rapidly the Green’s function varies (Fig. 1 validates this with the full g(∆η), not itsδ-function limit). Theδ-function limit is the simplest path to the result; the deeper property is separability. 5
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) 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. This bound assumes the macroscopic envelope is smooth (differentiable onH −1 ∗ ); ifShas sharper features at the transition boundaries varying onβ −1 , the bound is weakened to| e I|/I 0 ≲ β/k ∼ v w , which isO(1) rather than small. The deeper argument that no log-periodic contamination is generated rests on the spectral-separation property ( e Ioscillates in lineark, not inlnk) discussed below, which holds regardless of the envelope smoothness. For sub-horizon GW modes k ≫ H ∗ , the ratio | e I|/I 0 ≲ H ∗ /k for smooth envelopes. 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. 6
Table 1: Relative error bound on the factorizationP h =C(k)P 0 h , forε= 0.1,v w = 1. The factorization-step (phenomenological) corrections are: (i) Delta-function approximation, Eq. (9): relative error≤ τ corr H ∗ ≡ v w /(β/H ∗ ) for smooth envelopes. (ii) Bulk convolution slow variation, Sec. 5.3: relative error≤(2π/ lnb)ε f τ corr H ∗ forb≥1.5, using the uniformly bounded absolute derivative|dδ/d lnq|≤ ε f ×2π/ lnb. (iii) Green’s function oscillations, Eq. (13): relative error≤ H ∗ /k ≡ τ corr H ∗ (Riemann–Lebesgue, smooth envelopes). Errors (i)–(iii) are quantitatively controlled and shown below. In addition, the WTC convolution carries a separateO(m V /q ∗ ) mass-gap correction (Sec. 5.3) that is model-dependent and is the dominant error for the WTC parameter band (∼10% form V /q ∗ ∼0.1); it is not included in the phenomenological budget below. β/H ∗ τ corr H ∗ Combined bound Dominant term 100.10≲ 15%δ-fn 1000.01≲ 1.5%convolution 10000.001≲ 0.15%Green’s fn 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 sourceSand a Gaussian temporal correlationg(∆η) of widthτ corr . This is a toy validation of the factorization structure rather than a full WTC simulation: it tests whether the Green’s-function convolution and time integration in Eq. (4) preserve ak-dependent modulation imposed at the UETC level, and is not intended to validate the microphysics of any specific BSM source. Crucially the numerical computation uses the full Gaussian, not theδ-function limit invoked in Eq. (9). Figure 1 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. 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) 7
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 1: 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. 8
is a sinusoid inlnfwith period ∆lnf=lnb, amplitudeε, and phaseφ 0 at leading order inτ corr H ∗ . By construction, its Pearson correlation coefficient with a fixed-period cosine template equalsr= 1.00 for any bandwidth spanning complete log-periods, modulo the O(τ corr H ∗ ) factorization corrections quantified in Table 1. 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 leading-order 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 2 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) 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) 9
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 2: 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. 10
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. 3 and 4, with SNR contours at {1, 5, 10, 20} and the WTC prediction band overlaid. 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 3: 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. 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 ingredi- ents — DSI in near-conformal systems, walking dynamics, holographic warp factors, radion 11
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 4: Same as Fig. 3, 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. 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.2–5.4 Holographic origin of periodic warp factorconjectural motivation5.2 12
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 therefore be read as a plausibility argument for a WTC-style UV completion, not as a derivation from first principles. 5.1 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.2 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 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) 13
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. 5.3 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) 14
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. For a massive gauge-boson propagatorD 0 (q) = 1/(q 2 +m 2 V ), this region hasD 0 (|k−p|)∼1/m 2 V , which is in fact larger than the bulk valueD 0 (q ∗ )∼1/q 2 ∗ by the factorq 2 ∗ /m 2 V ≫ 1 when q ∗ ≫ m V . The integrated contribution of the small-|k−p|region is nevertheless suppressed by the three-dimensional phase-space measurep 2 dp: contributions from|k−p|≲ m V scale asm 3 V ×D 0 (0)×D 0 (q ∗ )∼ m V /q 2 ∗ , compared to the bulk contribution at|k−p|∼ q ∗ scaling asq 3 ∗ × D 0 (q ∗ ) 2 ∼1/q ∗ . The ratio is (boundary)/(bulk)∼ m V /q ∗ ≡ ρ 1/2 , with ρ≡ m 2 V /q 2 ∗ ≪ 1. The combined bound is δ(p)≈ δ(|k−p|)≈ δ(k) 1 +O(ε f τ corr H ∗ ) +O(m V /q ∗ ) ,(37) withm V /q ∗ ≪1 in the WTC regime. For typical WTC benchmark parametersm V /q ∗ ∼ 0.1, the boundary contribution sets a residual relative error of order 10% on the convolution factorization, exceeding theO(τ corr H ∗ ) bulk correction. We retainρ 1/2 ≡ m V /q ∗ in error budgets below. The cross-term contributionsδ(p) andδ(|k−p|) each factor asδ(k)×Π 0 (k,η,η ′ ) over the dominant support, giving a total cross-term contribution of 2δ(k)Π 0 from the two linear terms. This means the UETC modulation depthεinherited from the propagator modulationδ=ε f cos(···) satisfiesε= 2ε f × c geom at leading order, wherec geom =O(1) is a geometric coefficient arising from the angular average ofP TT over the convolution support. We do not computec geom explicitly; for the WTC band of Sec. 5.4 we takeε∼ ε f as an order-of-magnitude estimate, with the understanding that the precise mapping carries an O(1) uncertainty. Hence Π(k,η,η ′ ) = Π 0 (k,η,η ′ ) 1 + ε cos 2π ln(k/k ∗ ) lnb
- φ 0 1 +O(ετ corr H ∗ ) +O(m V /q ∗ ) , (38) withε∼ ε f andb=b 0 at leading order (up toO(1) model-dependent factors). For WTC benchmark parametersβ/H ∗ ≳100 andm V /q ∗ ∼0.1, the combined relative correction is dominated by the mass-gap term at ∼ 10%, not the ∼ 1% bulk error. 5.4 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 from the factorization step andO(m V /q ∗ ) mass- gap corrections from the WTC convolution. 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. 3–4 and overlaps the high-SNR region of the LISA detectability forecast. The map from WTC parameters to (b,ε) follows from the holographic benchmark 15
of Sec. 5.2:ε= 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). Approximation hierarchy and model-dependent factors. The chain from a WTC- type Lagrangian to the observable Ω GW (f) runs through several steps that we have sketched but not derived from first principles; each carries anO(1) model-dependent multiplicative factor. In particular: (i) The Goldberger–Wise-type relationε f = 4δA 0 between the holographic warp-factor perturbation and the 4D technidilaton potential modulation [Eq. (31)] depends on the specific radion-dilaton identification in the 5D dual. (ii) The chain-rule transfer of the periodic modulation fromV(φ) to the gauge propagator D(q) [Eq. (32)] involves the technidilaton scaling dimension ∆ φ , which isO(1) near the quasi-fixed point but is not exactly unity; this introduces anO(1) shift in the inherited log-periodbrelative tob 0 . (iii) The cross-term factor of 2 from the convolution Eq. (34) and the angular average of the TT projectorP TT give anO(1) geometric factorc geom in the relationε∼ c geom ε f , which we have not computed. (iv) The mass-gap suppression of the small-|k−p|region of the convolution scales asm V /q ∗ rather than (m V /q ∗ ) 2 once the three-dimensional phase-space measure is included; for typical WTC parameters this is ∼10%. Taken together, these factors mean the WTC prediction band Eq. (39) should be understood as an order-of-magnitude forecast rather than a precise mapping: the central values are correct up toO(1) model-dependent rescalings, and a precise determination of the (ε,b) band requires explicit microphysical computation that is beyond the scope of this work. The phenomenological chain of Sec. 3–4 — DSI in the UETC implies a log-periodic SGWB template at theτ corr H ∗ level — is independent of these model-dependent factors, and is what the matched-filter observable in Eq. (16) tests. 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.3) — 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 16
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 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. 17
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 ∗ )]. 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. The phenomenological chain connecting a WTC-inspired DSI modulation to the observable Ω GW (f) is internally consistent: a convolution calculation confirms that the potential modulation propagates to the SGWB up to controlledO(τ corr H ∗ ) corrections andO(1) model-dependent factors detailed in Sec. 5.4. The WTC prediction bandε ∈[0.04,0.18],b ∈[1.7,2.8] — understood as an order-of-magnitude forecast rather than a precise mapping — 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]. 18
[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]. [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). 19
⚑Derivation Flags (23)
- highEq. (15) — The main factorization theorem is exact only when the UETC modulation is a k-only multiplicative factor. The extension to realistic nonseparable corrections is asserted through error estimates rather than proved in a general norm or with explicit hypotheses on the UETC.
If wrong: The universal claim that DSI in the source passes to the observable spectrum at percent level would be reduced to a conditional statement valid only for the assumed separable ansatz Eq. (6).
- highEq. (31) — The relation V_4(phi) approximately equals V_CW(phi)[1 + 4 delta A_0 cos(...)] with epsilon_f = 4 delta A_0 and ln b_0 = kL/n_p is stated as a result of a holographic warp-factor modulation, but no derivation from the 5D action or radion effective potential is provided.
If wrong: The claimed WTC/holographic origin of the numerical epsilon_f and b_0 ranges would not be established, and Eq. (39)'s parameter band would lose its stated basis.
- highEq. (32) — The transfer from a periodic technidilaton potential modulation to a multiplicative gauge-field propagator modulation is justified only by a brief chain-rule argument. The dependence of the propagator on the full effective action, the scaling dimension Delta_phi, and possible phase/amplitude renormalization is not derived.
If wrong: The DSI modulation may not appear in D(q) with the same amplitude or period, so the subsequent UETC convolution and the observable WTC prediction would not follow.
- highEq. (32) (propagator modulation D=D0[1+δ(q)] derived from potential modulation) — ‘Chain rule’ argument from δV(φ) to δm_V^2 to δ(q) is sketched but not shown; in particular, how a modulation of the effective potential induces a log-periodic multiplicative modulation of the momentum-space propagator at finite temperature is not derived.
If wrong: The UV-completion claim that WTC produces a DSI-modulated UETC (Eq. (38)) is not established; predicted ε,b ranges may not follow.
- highEq. (33)–(38) (convolution factorization and boundary-region estimate ∼ m_V/q_*) — Convolution estimates are parametric and partially heuristic: (i) assumes D0(p) sharply peaked at p∼q_* with width Δp/p∼τ_corr H_*; (ii) replaces δ(p) by δ(k) over support using a derivative bound; (iii) estimates the small-|k−p| contribution scaling as m_V/q_* without a detailed integral or treatment of TT projector angular structure; c_geom is left uncomputed while later used to set ε∼ε_f.
If wrong: Eq. (38) may not hold even approximately; the modulation depth ε could be much smaller, phase-shifted, or not purely multiplicative in k, directly impacting the WTC ‘candidate UV completion’ quantitative predictions.
- highEq. (39) — The map from WTC parameters to epsilon in [0.04, 0.18] and b in [1.7, 2.8] depends on uncomputed O(1) factors: the radion-dilaton mapping, Delta_phi, the TT angular average c_geom, and the mass-gap correction. The paper later calls the band order-of-magnitude, but presents it earlier as a sharp prediction.
If wrong: The claimed falsifiable WTC target region in the detectability plane may shift by O(1), possibly changing the detectability and model-discrimination conclusions.
- highEq. (9) (F(k,Δη) ≃ F(k) δ(η−η′) + O(τ_corr H_*)) — Replacement of a narrow correlation function by a delta function is stated with an error term but without a quantitative bound propagated through Eq. (4) (e.g., using moments of g(Δη) and bounds on Green’s function derivatives).
If wrong: Eq. (10)–(15) may not hold even approximately; the modulation could be distorted in amplitude/phase or mixed with additional k-dependence, undermining the central template Eq. (16).
- highEqs. (35)-(38) — The convolution argument assumes that the internal-momentum support allows delta(p) and delta(|k-p|) to be replaced by delta(k). This is not demonstrated. The stated propagator D_0(q)=1/(q^2+m_V^2) is not sharply peaked at q_*; moreover, a convolution over p generally samples p and |k-p| over a broad angular domain, so a log-periodic function of internal momentum need not reproduce the same log-periodic function of external k.
If wrong: The central claim that WTC microphysical DSI propagates to the UETC with epsilon approximately epsilon_f and b approximately b_0 fails. The WTC prediction band and its placement in the LISA SNR plots become unsupported.
- highSec. 3.2 statement: 'evaluate G_k(η,η2) at η2=η1' and discussion around kτ_corr∼1 — The paper acknowledges that at k∼β/v_w one has kτ_corr∼1 and the slow-variation approximation is ‘marginal’ with O(1) corrections, but does not compute them; nonetheless it continues to use the factorized form with percent-level claims for β/H_*≳10.
If wrong: At/near the spectral peak—most important for detectability—the factorization amplitude could be misestimated by order unity, invalidating the stated percent-level transfer and detectability forecasts.
- mediumEq. (13)–(14) (integration by parts bound |Ĩ|/I0 ≲ H_*/k) — Bound assumes differentiability and Hubble-timescale variation of a^4 S; the paper notes that if S varies on β^{-1} the bound weakens to O(1), but no uniform bound is provided for realistic source turn-on/turn-off profiles.
If wrong: The claim that Green’s-function oscillatory term is a small correction may fail, affecting the clean multiplicative mapping from C(k) to P_h.
- mediumEq. (13), bound on Ĩ via integration by parts — Bound |Ĩ| ≤ (1/2k) ∫|∂_η1 [a^4 S]| dη_1 ≲ (H*/k) I_0 assumes boundary terms vanish and that |∂_η1[a^4 S]| ≲ H* a^4 S. The latter holds for smooth envelopes only; the paper acknowledges the bound degrades to O(1) for sharper features but then appeals to a 'spectral-separation' argument that is asserted rather than rigorously proven.
If wrong: If the spectral-separation argument fails for non-smooth macroscopic envelopes (e.g., sharp transition boundaries varying on β^{-1}), the no-log-periodic-contamination claim for P_h^0 weakens, and the cleanliness of the factorization result Eq. (15) for realistic phase transitions becomes uncertain.
- mediumEq. (29)–(31) (periodic modulation of V(φ) from warp-factor modulation) — Key mapping ε_f=4δA0 and ln b0=kL/n_p is stated as arising ‘at linear order’ but without derivation; depends on a specific radion/technidilaton identification.
If wrong: The claimed WTC parameter-to-(ε,b) mapping and the numerical ranges in Eq. (39) become unsupported (though the phenomenological DSI→template claim could still stand independently).
- mediumEq. (32), chain-rule propagation of modulation to the propagator — The transfer of the periodic modulation from V(φ) to the gauge propagator D(q) is sketched via m_V^2 ∝ ⟨φ⟩^2 and 'chain rule', but the relationship between the φ-dependence of V and the q-dependence of D(q) at scale q is not rigorously established. The identification q* ~ q_0 and the inheritance of b_0 are asserted with O(1) caveat.
If wrong: The inherited log-period in the UETC may differ from b_0 by an O(1) factor (the conformal dimension Δ_φ), shifting the WTC predicted band for b. The paper acknowledges this in Sec. 5.4(ii).
- mediumEq. (7) and Sec. 3.2 (UETC separability Π=S(η,η′)F(k,Δη)) — Separation of UETC into macroscopic S and microscopic F is asserted as ‘natural’ for phase transitions; no explicit conditions (stationarity, mixing bounds, or derivation from a known model UETC) are given.
If wrong: If Π is not approximately separable in this way, the later factorization of C(k) may fail, and the claimed multiplicative transfer to P_h and Ω_GW is not justified.
- mediumEq. (9) — The replacement F(k, Delta eta) approximately equals F(k) delta(eta - eta') plus O(tau_corr H_*) is sketched rather than derived. The error estimate does not clearly control dependence on k tau_corr, and the text itself notes that near the peak k tau_corr is about 1, where slow-variation is marginal.
If wrong: The quantitative error bounds attached to the short-correlation-time factorization would not follow. The exact algebraic factorization for a strictly separable k-only C(k) could remain valid, but the claimed O(tau_corr H_*) control for realistic non-delta UETCs would be unsupported.
- mediumEqs. (12)-(14) — The integration-by-parts bound on the oscillatory Green-function contribution assumes a smooth source envelope with vanishing boundary terms. The paper acknowledges that sharper transition-boundary features weaken the bound to O(1), but then continues to use the small O(tau_corr H_*) estimate for the factorization theorem.
If wrong: The claim that Green-function oscillations cannot contaminate or alter the log-periodic residual at the stated small level would be unproven for nonsmooth or rapidly varying macroscopic source profiles.
- mediumSec. 3.2 ‘spectral separation’ argument (Δ ln k ∼ 1/(kη_*) ≪ ln b ⇒ no log-periodic contamination) — Qualitative argument that linear-k oscillations cannot generate log-periodic features at period ln b is asserted without a proof that windowing/convolution/integration cannot produce aliasing or beat phenomena in ln k.
If wrong: Even if C(k) factors approximately, additional structure in ln k could mimic or bias the extracted (b,φ0), weakening falsifiability claims.
- mediumSec. 5.3, mass-gap suppression argument (boundary contribution ~ m_V/q*) — The scaling (boundary)/(bulk) ~ m_V/q* from phase-space measure is presented as a heuristic estimate without explicit integration. The crossover between the bulk and boundary regions and the precise coefficient are not computed.
If wrong: The ~10% error budget for the WTC convolution could be larger; would primarily affect the precision of the WTC band, not the existence of the log-periodic template.
- lowEq. (22)→(25) (SNR_osc scaling with ⟨cos^2⟩=1/2 and N_periods) — Uses averaging ⟨cos^2⟩=1/2 over ‘complete log-periods’ and defines SNR_bin as an integral over width ln b; does not address mismatch when Ω0/σ varies significantly within a period, nor correlations between periods due to colored noise and overlapping window functions used in practice.
If wrong: Detectability scaling could deviate by O(1) factors; the qualitative √N_periods enhancement likely remains but the contour values in Figs. 3–4 could shift.
- lowEq. (25) — The matched-filter scaling SNR_osc = (epsilon/sqrt(2)) sqrt(N_periods) SNR_bin assumes complete log-periods, a fixed known template phase and b, and effectively comparable baseline SNR per log-period. The derivation does not include look-elsewhere penalties or variation of noise and baseline shape across the band.
If wrong: The plotted SNR contours would be optimistic or require reinterpretation, but the existence of a log-periodic template would not be affected.
- lowEq. (31), holographic derivation of V_4(φ) — The 4D potential modulation with ε_f = 4δA_0 and ln b_0 = kL/n_p is stated as following from the warp-factor perturbation Eq. (30) 'at linear order' but the derivation is not shown; cited Ref. [17] (Goldberger-Wise) does not contain this specific result for periodic warp factors.
If wrong: The specific numerical mapping of WTC parameters to (ε, b) in Eq. (39) would carry additional uncertainty, but this is already disclosed as 'order-of-magnitude' in Sec. 5.4. The phenomenological backbone is unaffected.
- lowEq. (36), bound |δ(p) - δ(k)| ≲ (2π ε_f/ ln b_0) τ_corr H* — The bound combines the uniform absolute derivative |dδ/d ln q| ≤ ε_f · 2π/ln b_0 with support width τ_corr H*, but the support width is identified with relative half-width ∆p/p of D_0 without justification linking it directly to τ_corr H*.
If wrong: Bound coefficient may be off by O(1); does not affect the structural result.
- lowNumerical validation, Fig. 1, claim of 10^-16 agreement — The pointwise identical residual across τ_corr H* ∈ {0.01, 0.05, 0.20} is consistent with the algebraic factorization property only when C(k) depends purely on k. This is essentially a check of algebra, not of the physical approximations bounded in Table 1 — as the paper itself notes. The numerical check thus validates only the trivial factorization step.
If wrong: The numerical demonstration is honest about its limited scope; no central claim is at risk.
Within the paper’s own framework, the cleanest mathematically valid statement is: if the source UETC can be written as Π(k,η1,η2)=C(k)Π0(k,η1,η2) with C(k) independent of time arguments, then the same multiplicative factor C(k) appears in P_h(k,η) and hence in Ω_GW(f), modulo any later processing that does not mix k. That algebraic point is correct and is nicely separated from the microphysics. What is not established with the same rigor is the stronger, quantitative claim that generic first-order phase transition sources with τ_corr H_*≪1 (or β/H_*≳10) satisfy the required separability to percent accuracy near the most relevant frequencies, especially given the paper’s own admission that kτ_corr∼1 at the peak makes the key approximation marginal and potentially O(1). The UV-completion section is explicitly presented as a candidate, and mathematically it remains at the level of parametric scaling estimates; the mapping to the stated (ε,b) band is therefore not a derived consequence of the equations shown.
⚑Derivation Flags (23)
- highEq. (15) — The main factorization theorem is exact only when the UETC modulation is a k-only multiplicative factor. The extension to realistic nonseparable corrections is asserted through error estimates rather than proved in a general norm or with explicit hypotheses on the UETC.
If wrong: The universal claim that DSI in the source passes to the observable spectrum at percent level would be reduced to a conditional statement valid only for the assumed separable ansatz Eq. (6).
- highEq. (31) — The relation V_4(phi) approximately equals V_CW(phi)[1 + 4 delta A_0 cos(...)] with epsilon_f = 4 delta A_0 and ln b_0 = kL/n_p is stated as a result of a holographic warp-factor modulation, but no derivation from the 5D action or radion effective potential is provided.
If wrong: The claimed WTC/holographic origin of the numerical epsilon_f and b_0 ranges would not be established, and Eq. (39)'s parameter band would lose its stated basis.
- highEq. (32) — The transfer from a periodic technidilaton potential modulation to a multiplicative gauge-field propagator modulation is justified only by a brief chain-rule argument. The dependence of the propagator on the full effective action, the scaling dimension Delta_phi, and possible phase/amplitude renormalization is not derived.
If wrong: The DSI modulation may not appear in D(q) with the same amplitude or period, so the subsequent UETC convolution and the observable WTC prediction would not follow.
- highEq. (32) (propagator modulation D=D0[1+δ(q)] derived from potential modulation) — ‘Chain rule’ argument from δV(φ) to δm_V^2 to δ(q) is sketched but not shown; in particular, how a modulation of the effective potential induces a log-periodic multiplicative modulation of the momentum-space propagator at finite temperature is not derived.
If wrong: The UV-completion claim that WTC produces a DSI-modulated UETC (Eq. (38)) is not established; predicted ε,b ranges may not follow.
- highEq. (33)–(38) (convolution factorization and boundary-region estimate ∼ m_V/q_*) — Convolution estimates are parametric and partially heuristic: (i) assumes D0(p) sharply peaked at p∼q_* with width Δp/p∼τ_corr H_*; (ii) replaces δ(p) by δ(k) over support using a derivative bound; (iii) estimates the small-|k−p| contribution scaling as m_V/q_* without a detailed integral or treatment of TT projector angular structure; c_geom is left uncomputed while later used to set ε∼ε_f.
If wrong: Eq. (38) may not hold even approximately; the modulation depth ε could be much smaller, phase-shifted, or not purely multiplicative in k, directly impacting the WTC ‘candidate UV completion’ quantitative predictions.
- highEq. (39) — The map from WTC parameters to epsilon in [0.04, 0.18] and b in [1.7, 2.8] depends on uncomputed O(1) factors: the radion-dilaton mapping, Delta_phi, the TT angular average c_geom, and the mass-gap correction. The paper later calls the band order-of-magnitude, but presents it earlier as a sharp prediction.
If wrong: The claimed falsifiable WTC target region in the detectability plane may shift by O(1), possibly changing the detectability and model-discrimination conclusions.
- highEq. (9) (F(k,Δη) ≃ F(k) δ(η−η′) + O(τ_corr H_*)) — Replacement of a narrow correlation function by a delta function is stated with an error term but without a quantitative bound propagated through Eq. (4) (e.g., using moments of g(Δη) and bounds on Green’s function derivatives).
If wrong: Eq. (10)–(15) may not hold even approximately; the modulation could be distorted in amplitude/phase or mixed with additional k-dependence, undermining the central template Eq. (16).
- highEqs. (35)-(38) — The convolution argument assumes that the internal-momentum support allows delta(p) and delta(|k-p|) to be replaced by delta(k). This is not demonstrated. The stated propagator D_0(q)=1/(q^2+m_V^2) is not sharply peaked at q_*; moreover, a convolution over p generally samples p and |k-p| over a broad angular domain, so a log-periodic function of internal momentum need not reproduce the same log-periodic function of external k.
If wrong: The central claim that WTC microphysical DSI propagates to the UETC with epsilon approximately epsilon_f and b approximately b_0 fails. The WTC prediction band and its placement in the LISA SNR plots become unsupported.
- highSec. 3.2 statement: 'evaluate G_k(η,η2) at η2=η1' and discussion around kτ_corr∼1 — The paper acknowledges that at k∼β/v_w one has kτ_corr∼1 and the slow-variation approximation is ‘marginal’ with O(1) corrections, but does not compute them; nonetheless it continues to use the factorized form with percent-level claims for β/H_*≳10.
If wrong: At/near the spectral peak—most important for detectability—the factorization amplitude could be misestimated by order unity, invalidating the stated percent-level transfer and detectability forecasts.
- mediumEq. (13)–(14) (integration by parts bound |Ĩ|/I0 ≲ H_*/k) — Bound assumes differentiability and Hubble-timescale variation of a^4 S; the paper notes that if S varies on β^{-1} the bound weakens to O(1), but no uniform bound is provided for realistic source turn-on/turn-off profiles.
If wrong: The claim that Green’s-function oscillatory term is a small correction may fail, affecting the clean multiplicative mapping from C(k) to P_h.
- mediumEq. (13), bound on Ĩ via integration by parts — Bound |Ĩ| ≤ (1/2k) ∫|∂_η1 [a^4 S]| dη_1 ≲ (H*/k) I_0 assumes boundary terms vanish and that |∂_η1[a^4 S]| ≲ H* a^4 S. The latter holds for smooth envelopes only; the paper acknowledges the bound degrades to O(1) for sharper features but then appeals to a 'spectral-separation' argument that is asserted rather than rigorously proven.
If wrong: If the spectral-separation argument fails for non-smooth macroscopic envelopes (e.g., sharp transition boundaries varying on β^{-1}), the no-log-periodic-contamination claim for P_h^0 weakens, and the cleanliness of the factorization result Eq. (15) for realistic phase transitions becomes uncertain.
- mediumEq. (29)–(31) (periodic modulation of V(φ) from warp-factor modulation) — Key mapping ε_f=4δA0 and ln b0=kL/n_p is stated as arising ‘at linear order’ but without derivation; depends on a specific radion/technidilaton identification.
If wrong: The claimed WTC parameter-to-(ε,b) mapping and the numerical ranges in Eq. (39) become unsupported (though the phenomenological DSI→template claim could still stand independently).
- mediumEq. (32), chain-rule propagation of modulation to the propagator — The transfer of the periodic modulation from V(φ) to the gauge propagator D(q) is sketched via m_V^2 ∝ ⟨φ⟩^2 and 'chain rule', but the relationship between the φ-dependence of V and the q-dependence of D(q) at scale q is not rigorously established. The identification q* ~ q_0 and the inheritance of b_0 are asserted with O(1) caveat.
If wrong: The inherited log-period in the UETC may differ from b_0 by an O(1) factor (the conformal dimension Δ_φ), shifting the WTC predicted band for b. The paper acknowledges this in Sec. 5.4(ii).
- mediumEq. (7) and Sec. 3.2 (UETC separability Π=S(η,η′)F(k,Δη)) — Separation of UETC into macroscopic S and microscopic F is asserted as ‘natural’ for phase transitions; no explicit conditions (stationarity, mixing bounds, or derivation from a known model UETC) are given.
If wrong: If Π is not approximately separable in this way, the later factorization of C(k) may fail, and the claimed multiplicative transfer to P_h and Ω_GW is not justified.
- mediumEq. (9) — The replacement F(k, Delta eta) approximately equals F(k) delta(eta - eta') plus O(tau_corr H_*) is sketched rather than derived. The error estimate does not clearly control dependence on k tau_corr, and the text itself notes that near the peak k tau_corr is about 1, where slow-variation is marginal.
If wrong: The quantitative error bounds attached to the short-correlation-time factorization would not follow. The exact algebraic factorization for a strictly separable k-only C(k) could remain valid, but the claimed O(tau_corr H_*) control for realistic non-delta UETCs would be unsupported.
- mediumEqs. (12)-(14) — The integration-by-parts bound on the oscillatory Green-function contribution assumes a smooth source envelope with vanishing boundary terms. The paper acknowledges that sharper transition-boundary features weaken the bound to O(1), but then continues to use the small O(tau_corr H_*) estimate for the factorization theorem.
If wrong: The claim that Green-function oscillations cannot contaminate or alter the log-periodic residual at the stated small level would be unproven for nonsmooth or rapidly varying macroscopic source profiles.
- mediumSec. 3.2 ‘spectral separation’ argument (Δ ln k ∼ 1/(kη_*) ≪ ln b ⇒ no log-periodic contamination) — Qualitative argument that linear-k oscillations cannot generate log-periodic features at period ln b is asserted without a proof that windowing/convolution/integration cannot produce aliasing or beat phenomena in ln k.
If wrong: Even if C(k) factors approximately, additional structure in ln k could mimic or bias the extracted (b,φ0), weakening falsifiability claims.
- mediumSec. 5.3, mass-gap suppression argument (boundary contribution ~ m_V/q*) — The scaling (boundary)/(bulk) ~ m_V/q* from phase-space measure is presented as a heuristic estimate without explicit integration. The crossover between the bulk and boundary regions and the precise coefficient are not computed.
If wrong: The ~10% error budget for the WTC convolution could be larger; would primarily affect the precision of the WTC band, not the existence of the log-periodic template.
- lowEq. (22)→(25) (SNR_osc scaling with ⟨cos^2⟩=1/2 and N_periods) — Uses averaging ⟨cos^2⟩=1/2 over ‘complete log-periods’ and defines SNR_bin as an integral over width ln b; does not address mismatch when Ω0/σ varies significantly within a period, nor correlations between periods due to colored noise and overlapping window functions used in practice.
If wrong: Detectability scaling could deviate by O(1) factors; the qualitative √N_periods enhancement likely remains but the contour values in Figs. 3–4 could shift.
- lowEq. (25) — The matched-filter scaling SNR_osc = (epsilon/sqrt(2)) sqrt(N_periods) SNR_bin assumes complete log-periods, a fixed known template phase and b, and effectively comparable baseline SNR per log-period. The derivation does not include look-elsewhere penalties or variation of noise and baseline shape across the band.
If wrong: The plotted SNR contours would be optimistic or require reinterpretation, but the existence of a log-periodic template would not be affected.
- lowEq. (31), holographic derivation of V_4(φ) — The 4D potential modulation with ε_f = 4δA_0 and ln b_0 = kL/n_p is stated as following from the warp-factor perturbation Eq. (30) 'at linear order' but the derivation is not shown; cited Ref. [17] (Goldberger-Wise) does not contain this specific result for periodic warp factors.
If wrong: The specific numerical mapping of WTC parameters to (ε, b) in Eq. (39) would carry additional uncertainty, but this is already disclosed as 'order-of-magnitude' in Sec. 5.4. The phenomenological backbone is unaffected.
- lowEq. (36), bound |δ(p) - δ(k)| ≲ (2π ε_f/ ln b_0) τ_corr H* — The bound combines the uniform absolute derivative |dδ/d ln q| ≤ ε_f · 2π/ln b_0 with support width τ_corr H*, but the support width is identified with relative half-width ∆p/p of D_0 without justification linking it directly to τ_corr H*.
If wrong: Bound coefficient may be off by O(1); does not affect the structural result.
- lowNumerical validation, Fig. 1, claim of 10^-16 agreement — The pointwise identical residual across τ_corr H* ∈ {0.01, 0.05, 0.20} is consistent with the algebraic factorization property only when C(k) depends purely on k. This is essentially a check of algebra, not of the physical approximations bounded in Table 1 — as the paper itself notes. The numerical check thus validates only the trivial factorization step.
If wrong: The numerical demonstration is honest about its limited scope; no central claim is at risk.
The phenomenological core is valid as a conditional result: if the source UETC contains a multiplicative, time-independent log-periodic factor C(k), then the linear gravitational-wave power-spectrum integral transmits that same factor to the observable spectrum. In that restricted sense, the route from Eq. (6) to Eq. (16) is mathematically coherent. The broader claims are much less secure. The short-correlation-time error estimates are not fully derived, and the WTC ultraviolet-completion chain contains load-bearing gaps in the holographic potential mapping, the potential-to-propagator transfer, and the anisotropic-stress convolution. These gaps directly affect the claimed WTC parameter band and its LISA detectability forecast. Thus the paper is internally plausible as a conditional phenomenological template, but it does not mathematically establish the advertised WTC prediction with the stated precision.
⚑Derivation Flags (23)
- highEq. (15) — The main factorization theorem is exact only when the UETC modulation is a k-only multiplicative factor. The extension to realistic nonseparable corrections is asserted through error estimates rather than proved in a general norm or with explicit hypotheses on the UETC.
If wrong: The universal claim that DSI in the source passes to the observable spectrum at percent level would be reduced to a conditional statement valid only for the assumed separable ansatz Eq. (6).
- highEq. (31) — The relation V_4(phi) approximately equals V_CW(phi)[1 + 4 delta A_0 cos(...)] with epsilon_f = 4 delta A_0 and ln b_0 = kL/n_p is stated as a result of a holographic warp-factor modulation, but no derivation from the 5D action or radion effective potential is provided.
If wrong: The claimed WTC/holographic origin of the numerical epsilon_f and b_0 ranges would not be established, and Eq. (39)'s parameter band would lose its stated basis.
- highEq. (32) — The transfer from a periodic technidilaton potential modulation to a multiplicative gauge-field propagator modulation is justified only by a brief chain-rule argument. The dependence of the propagator on the full effective action, the scaling dimension Delta_phi, and possible phase/amplitude renormalization is not derived.
If wrong: The DSI modulation may not appear in D(q) with the same amplitude or period, so the subsequent UETC convolution and the observable WTC prediction would not follow.
- highEq. (32) (propagator modulation D=D0[1+δ(q)] derived from potential modulation) — ‘Chain rule’ argument from δV(φ) to δm_V^2 to δ(q) is sketched but not shown; in particular, how a modulation of the effective potential induces a log-periodic multiplicative modulation of the momentum-space propagator at finite temperature is not derived.
If wrong: The UV-completion claim that WTC produces a DSI-modulated UETC (Eq. (38)) is not established; predicted ε,b ranges may not follow.
- highEq. (33)–(38) (convolution factorization and boundary-region estimate ∼ m_V/q_*) — Convolution estimates are parametric and partially heuristic: (i) assumes D0(p) sharply peaked at p∼q_* with width Δp/p∼τ_corr H_*; (ii) replaces δ(p) by δ(k) over support using a derivative bound; (iii) estimates the small-|k−p| contribution scaling as m_V/q_* without a detailed integral or treatment of TT projector angular structure; c_geom is left uncomputed while later used to set ε∼ε_f.
If wrong: Eq. (38) may not hold even approximately; the modulation depth ε could be much smaller, phase-shifted, or not purely multiplicative in k, directly impacting the WTC ‘candidate UV completion’ quantitative predictions.
- highEq. (39) — The map from WTC parameters to epsilon in [0.04, 0.18] and b in [1.7, 2.8] depends on uncomputed O(1) factors: the radion-dilaton mapping, Delta_phi, the TT angular average c_geom, and the mass-gap correction. The paper later calls the band order-of-magnitude, but presents it earlier as a sharp prediction.
If wrong: The claimed falsifiable WTC target region in the detectability plane may shift by O(1), possibly changing the detectability and model-discrimination conclusions.
- highEq. (9) (F(k,Δη) ≃ F(k) δ(η−η′) + O(τ_corr H_*)) — Replacement of a narrow correlation function by a delta function is stated with an error term but without a quantitative bound propagated through Eq. (4) (e.g., using moments of g(Δη) and bounds on Green’s function derivatives).
If wrong: Eq. (10)–(15) may not hold even approximately; the modulation could be distorted in amplitude/phase or mixed with additional k-dependence, undermining the central template Eq. (16).
- highEqs. (35)-(38) — The convolution argument assumes that the internal-momentum support allows delta(p) and delta(|k-p|) to be replaced by delta(k). This is not demonstrated. The stated propagator D_0(q)=1/(q^2+m_V^2) is not sharply peaked at q_*; moreover, a convolution over p generally samples p and |k-p| over a broad angular domain, so a log-periodic function of internal momentum need not reproduce the same log-periodic function of external k.
If wrong: The central claim that WTC microphysical DSI propagates to the UETC with epsilon approximately epsilon_f and b approximately b_0 fails. The WTC prediction band and its placement in the LISA SNR plots become unsupported.
- highSec. 3.2 statement: 'evaluate G_k(η,η2) at η2=η1' and discussion around kτ_corr∼1 — The paper acknowledges that at k∼β/v_w one has kτ_corr∼1 and the slow-variation approximation is ‘marginal’ with O(1) corrections, but does not compute them; nonetheless it continues to use the factorized form with percent-level claims for β/H_*≳10.
If wrong: At/near the spectral peak—most important for detectability—the factorization amplitude could be misestimated by order unity, invalidating the stated percent-level transfer and detectability forecasts.
- mediumEq. (13)–(14) (integration by parts bound |Ĩ|/I0 ≲ H_*/k) — Bound assumes differentiability and Hubble-timescale variation of a^4 S; the paper notes that if S varies on β^{-1} the bound weakens to O(1), but no uniform bound is provided for realistic source turn-on/turn-off profiles.
If wrong: The claim that Green’s-function oscillatory term is a small correction may fail, affecting the clean multiplicative mapping from C(k) to P_h.
- mediumEq. (13), bound on Ĩ via integration by parts — Bound |Ĩ| ≤ (1/2k) ∫|∂_η1 [a^4 S]| dη_1 ≲ (H*/k) I_0 assumes boundary terms vanish and that |∂_η1[a^4 S]| ≲ H* a^4 S. The latter holds for smooth envelopes only; the paper acknowledges the bound degrades to O(1) for sharper features but then appeals to a 'spectral-separation' argument that is asserted rather than rigorously proven.
If wrong: If the spectral-separation argument fails for non-smooth macroscopic envelopes (e.g., sharp transition boundaries varying on β^{-1}), the no-log-periodic-contamination claim for P_h^0 weakens, and the cleanliness of the factorization result Eq. (15) for realistic phase transitions becomes uncertain.
- mediumEq. (29)–(31) (periodic modulation of V(φ) from warp-factor modulation) — Key mapping ε_f=4δA0 and ln b0=kL/n_p is stated as arising ‘at linear order’ but without derivation; depends on a specific radion/technidilaton identification.
If wrong: The claimed WTC parameter-to-(ε,b) mapping and the numerical ranges in Eq. (39) become unsupported (though the phenomenological DSI→template claim could still stand independently).
- mediumEq. (32), chain-rule propagation of modulation to the propagator — The transfer of the periodic modulation from V(φ) to the gauge propagator D(q) is sketched via m_V^2 ∝ ⟨φ⟩^2 and 'chain rule', but the relationship between the φ-dependence of V and the q-dependence of D(q) at scale q is not rigorously established. The identification q* ~ q_0 and the inheritance of b_0 are asserted with O(1) caveat.
If wrong: The inherited log-period in the UETC may differ from b_0 by an O(1) factor (the conformal dimension Δ_φ), shifting the WTC predicted band for b. The paper acknowledges this in Sec. 5.4(ii).
- mediumEq. (7) and Sec. 3.2 (UETC separability Π=S(η,η′)F(k,Δη)) — Separation of UETC into macroscopic S and microscopic F is asserted as ‘natural’ for phase transitions; no explicit conditions (stationarity, mixing bounds, or derivation from a known model UETC) are given.
If wrong: If Π is not approximately separable in this way, the later factorization of C(k) may fail, and the claimed multiplicative transfer to P_h and Ω_GW is not justified.
- mediumEq. (9) — The replacement F(k, Delta eta) approximately equals F(k) delta(eta - eta') plus O(tau_corr H_*) is sketched rather than derived. The error estimate does not clearly control dependence on k tau_corr, and the text itself notes that near the peak k tau_corr is about 1, where slow-variation is marginal.
If wrong: The quantitative error bounds attached to the short-correlation-time factorization would not follow. The exact algebraic factorization for a strictly separable k-only C(k) could remain valid, but the claimed O(tau_corr H_*) control for realistic non-delta UETCs would be unsupported.
- mediumEqs. (12)-(14) — The integration-by-parts bound on the oscillatory Green-function contribution assumes a smooth source envelope with vanishing boundary terms. The paper acknowledges that sharper transition-boundary features weaken the bound to O(1), but then continues to use the small O(tau_corr H_*) estimate for the factorization theorem.
If wrong: The claim that Green-function oscillations cannot contaminate or alter the log-periodic residual at the stated small level would be unproven for nonsmooth or rapidly varying macroscopic source profiles.
- mediumSec. 3.2 ‘spectral separation’ argument (Δ ln k ∼ 1/(kη_*) ≪ ln b ⇒ no log-periodic contamination) — Qualitative argument that linear-k oscillations cannot generate log-periodic features at period ln b is asserted without a proof that windowing/convolution/integration cannot produce aliasing or beat phenomena in ln k.
If wrong: Even if C(k) factors approximately, additional structure in ln k could mimic or bias the extracted (b,φ0), weakening falsifiability claims.
- mediumSec. 5.3, mass-gap suppression argument (boundary contribution ~ m_V/q*) — The scaling (boundary)/(bulk) ~ m_V/q* from phase-space measure is presented as a heuristic estimate without explicit integration. The crossover between the bulk and boundary regions and the precise coefficient are not computed.
If wrong: The ~10% error budget for the WTC convolution could be larger; would primarily affect the precision of the WTC band, not the existence of the log-periodic template.
- lowEq. (22)→(25) (SNR_osc scaling with ⟨cos^2⟩=1/2 and N_periods) — Uses averaging ⟨cos^2⟩=1/2 over ‘complete log-periods’ and defines SNR_bin as an integral over width ln b; does not address mismatch when Ω0/σ varies significantly within a period, nor correlations between periods due to colored noise and overlapping window functions used in practice.
If wrong: Detectability scaling could deviate by O(1) factors; the qualitative √N_periods enhancement likely remains but the contour values in Figs. 3–4 could shift.
- lowEq. (25) — The matched-filter scaling SNR_osc = (epsilon/sqrt(2)) sqrt(N_periods) SNR_bin assumes complete log-periods, a fixed known template phase and b, and effectively comparable baseline SNR per log-period. The derivation does not include look-elsewhere penalties or variation of noise and baseline shape across the band.
If wrong: The plotted SNR contours would be optimistic or require reinterpretation, but the existence of a log-periodic template would not be affected.
- lowEq. (31), holographic derivation of V_4(φ) — The 4D potential modulation with ε_f = 4δA_0 and ln b_0 = kL/n_p is stated as following from the warp-factor perturbation Eq. (30) 'at linear order' but the derivation is not shown; cited Ref. [17] (Goldberger-Wise) does not contain this specific result for periodic warp factors.
If wrong: The specific numerical mapping of WTC parameters to (ε, b) in Eq. (39) would carry additional uncertainty, but this is already disclosed as 'order-of-magnitude' in Sec. 5.4. The phenomenological backbone is unaffected.
- lowEq. (36), bound |δ(p) - δ(k)| ≲ (2π ε_f/ ln b_0) τ_corr H* — The bound combines the uniform absolute derivative |dδ/d ln q| ≤ ε_f · 2π/ln b_0 with support width τ_corr H*, but the support width is identified with relative half-width ∆p/p of D_0 without justification linking it directly to τ_corr H*.
If wrong: Bound coefficient may be off by O(1); does not affect the structural result.
- lowNumerical validation, Fig. 1, claim of 10^-16 agreement — The pointwise identical residual across τ_corr H* ∈ {0.01, 0.05, 0.20} is consistent with the algebraic factorization property only when C(k) depends purely on k. This is essentially a check of algebra, not of the physical approximations bounded in Table 1 — as the paper itself notes. The numerical check thus validates only the trivial factorization step.
If wrong: The numerical demonstration is honest about its limited scope; no central claim is at risk.
This paper presents a two-tier theoretical framework with contrasting levels of completeness. The phenomenological core (discrete scale invariance in phase transition sources leading to log-periodic gravitational wave signatures) is mathematically complete with controlled approximations and explicit error bounds. The factorization theorem is properly derived and numerically validated. However, the walking technicolor ultraviolet completion represents a candidate realization rather than a complete microscopic theory. While the authors are appropriately transparent about these limitations, the incomplete UV completion prevents the work from achieving full theoretical completeness. The phenomenological framework alone provides a testable prediction for gravitational wave observations, but the specific parameter predictions require the incompletely developed WTC model.
This paper is strongest as a phenomenological analysis of what follows if a first-order phase-transition source UETC carries a k-dependent DSI modulation. Within that stated setup, the argument is fairly complete: the modulation template is explicit, the detectability scaling is concrete, and the approximation structure is more transparently discussed than in many speculative papers. The author is also commendably explicit about which parts are derived and which are conjectural. Where the submission is less complete is in the advertised UV completion. The walking-technicolor section reads as a plausible embedding rather than a quantitatively established realization, and the paper itself ultimately says so. That honesty prevents a harsher score, but the result is that the WTC-specific predictions are not fully supported at the level the abstract and introduction suggest. In short: the phenomenological core is substantially developed, but the concrete UV-completion layer remains partial and should be presented more consistently as suggestive rather than sharply derived.
This paper presents a phenomenological mechanism for imprinting DSI-induced log-periodic modulations on the SGWB, with a carefully argued factorization theorem and a concrete, though model-dependent, UV completion candidate. The derivation is logically coherent and well-documented within its own axioms, with approximations quantitatively bounded. The main limitations lie in the assumptions of UETC separability and the incompletely derived WTC-microphysics chain, but the paper acknowledges these explicitly and separates the more solid phenomenological backbone from the speculative UV layer. Overall, the work is substantially complete for its stated goals and provides a clear, testable prediction for LISA.
This is a scientifically interesting submission with a genuinely testable phenomenological idea: if a first-order cosmological phase transition source exhibits discrete scale invariance, the SGWB may carry a log-periodic residual across detector frequency space. That is a useful contribution because it suggests a specific spectral template and a detection strategy that could discriminate this framework from smoother SGWB models. The paper's strongest contribution is therefore not the walking-technicolor embedding, but the proposal that DSI could be searched for directly in gravitational-wave data as a structured oscillatory residual in ln f. The main weakness is communicative and epistemic calibration. The paper is more careful in the body than in the abstract: the UV completion is openly acknowledged as a candidate realization with O(1) uncertainties, yet the headline language presents it as sharper than the derivation supports. Likewise, the factorization argument is exact in a separable toy setting but only approximately justified for the physical case, and the numerical check does not fully address that gap. Overall, the paper is novel and meaningfully falsifiable, but its presentation would be stronger if it more cleanly separated the robust phenomenological template from the speculative UV-specific mapping.
This paper proposes that discrete scale invariance in the anisotropic stress of a first-order cosmological phase transition imprints a multiplicative log-periodic modulation on the SGWB, derives a factorization theorem under the short-correlation-time approximation, and offers walking technicolor as a candidate UV completion with a specific predicted band (ε ∈ [0.04, 0.18], b ∈ [1.7, 2.8]) that overlaps the high-SNR region of LISA. The work is communicatively strong: it is well-organized, explicitly distinguishes derived results from conjectural motivation (Table 2), quantifies its approximations (Table 1), and provides an independent numerical check. The observable template is concrete, falsifiable with a planned near-term observatory, and comes with a ready-to-implement matched-filter scheme — these are unambiguous merits. The principal scientific concerns are concentrated in the UV-completion layer rather than the phenomenological core. The chain from a WTC Lagrangian to the predicted (ε, b) band involves several uncomputed O(1) factors and a ~10% mass-gap correction that the author acknowledges but does not integrate into the headline error budget. The phenomenological backbone (DSI in UETC → log-periodic SGWB template) is the more robust result and is itself a novel, well-defined, and testable claim. Overall this is a carefully scoped contribution whose central observable is sharply falsifiable and whose author has been notably disciplined about labeling the epistemic status of each link in the argument.
Factorization theorem: the DSI modulation C(k) factors multiplicatively from the UETC into the tensor power spectrum up to controlled short-correlation-time corrections (Eq. 15).
Observable SGWB energy-density spectrum: smooth baseline multiplied by a log-periodic oscillatory modulation in ln f (main observable template; Eq. 16).
Matched-filter detectability scaling: oscillatory-component SNR in terms of modulation amplitude ε, number of log-periods in band, and baseline per-period SNR (Eq. 25).
The SGWB from a first-order cosmological phase transition with DSI in the source anisotropic stress will show a multiplicative log-periodic modulation Ω_GW(f)=Ω_0_GW(f)[1+ε cos(2π ln(f/f_*)/ln b + φ_0)] across the detector band.
Falsifiable if: A LISA (or comparable) measurement of the cosmological SGWB that, with matched-filter sensitivity accounting for baseline SNR, excludes a coherent sinusoidal modulation in ln f with amplitude ε above the instrument-specific upper limit (e.g. rules out ε values down to the projected sensitivity level in the relevant b-range).
Walking technicolor (WTC) can realize the DSI necessary to produce the SGWB modulation and predicts a parameter band ε∈[0.04,0.18], b∈[1.7,2.8] that falls in the high-SNR region for LISA.
Falsifiable if: A non-detection by LISA that, at the expected matched-filter sensitivity (e.g. SNR_osc threshold chosen for discovery), excludes the entire band ε∈[0.04,0.18] for b∈[1.7,2.8]; or model-specific microphysical calculations demonstrating that WTC microphysics cannot produce the required modulation amplitude or periodicity.
Under the short-correlation-time condition (τ_corr H_* ≪ 1, satisfied for β/H_* ≳ 10), the DSI modulation factors algebraically through the UETC to the observable spectrum with percent-level systematic error.
Falsifiable if: Empirical or numerical computation of the full double-time integral for realistic non-separable UETCs (with β/H_* ≳ 10) showing residuals or contamination of the log-periodic template exceeding the stated percent-level bound, or experimental detection of a mismatch between source-level DSI and observed modulation inconsistent with factorization bounds.
Share this Review
Post your AI review credential to social media, or copy the link to share anywhere.
theoryofeverything.ai/review-profile/paper/2d57a959-b0b0-427f-98c9-568d2873b23fThis 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