Enhanced Stochastic Gravitational Waves signals from Wess-Zumino chiral superfield

Technical Summary: Enhanced Stochastic Gravitational Waves from Wess-Zumino Chiral Superfields

Problem Statement
The paper investigates the potential for supersymmetric (SUSY) structures to leave observable imprints on the stochastic gravitational-wave (GW) background generated during the reheating era of the early Universe. While previous studies have established that inflaton decay into massive particles accompanied by graviton bremsstrahlung (three-body decay) produces a stochastic GW background, the amplitude of these signals in standard Yukawa-coupling scenarios is generally modest. Specifically, prior work (e.g., Ref. [16]) indicated that the GW strain amplitude exhibits only a slight dependence on the spin of the final-state particles, with differences well below one order of magnitude. The authors pose the question: if these particles are organized into multiplets transforming under internal symmetries, such as chiral multiplets in SUSY, does the associated GW amplitude exhibit universal behavior, or are there significant deviations?

Methodology
The authors construct a phenomenological model based on the Wess-Zumino chiral superfield formalism.

1. Interaction Construction: They formulate an effective phenomenological Lagrangian describing the coupling between a single inflaton ($\phi$) and the D-term sectors of a pair of chiral ($\Phi$) and anti-chiral ($\bar{\Phi}$) superfields. The interaction is given by $S_{int} = y_D \kappa \int d^8z \, \phi [\bar{\Phi}\Phi]_D$, where $y_D$ is a dimensionless coupling and $\kappa \equiv \sqrt{16\pi}/M_P$.
2. Component Expansion: Using superspace expansion techniques, they extract the cubic interaction terms in component fields. This yields two distinct interaction vertices:
   • A coupling between the inflaton and a pair of Majorana fermions ($\psi$): $S^{(\phi\psi\psi)}_{int} \propto \phi(\psi^\dagger \bar{\sigma}^\mu \partial_\mu \psi - \partial_\mu \psi^\dagger \bar{\sigma}^\mu \psi)$.
   • A coupling between the inflaton and a pair of complex scalars ($\phi$): $S^{(\phi\phi\phi)}_{int} \propto \phi(2\eta^{\mu\nu}\partial_\mu \phi^* \partial_\nu \phi + \dots)$.
     Crucially, these interactions contain derivative terms arising from the D-term construction, unlike standard non-derivative Yukawa couplings.
3. Decay Calculations:
   • Two-Body Decay: They derive the Feynman rules and compute the two-body decay widths ($\Gamma^{(0)}_{\phi \to \psi\psi}$ and $\Gamma^{(0)}_{\phi \to \phi\phi^*}$) to establish the baseline inflaton decay rate.
   • Three-Body Decay (GW Source): They analyze the three-body decay processes where the inflaton decays into a pair of particles (scalars or Majorana fermions) accompanied by the emission of a single massless graviton ($\phi \to \psi\psi g$ or $\phi \to \phi\phi^* g$). This involves calculating the squared amplitudes for all relevant Feynman diagrams, including those where the graviton is emitted from the inflaton leg and those where it is emitted from the final-state particle legs.
   • Formalism: The calculations employ two-component spinor techniques for Majorana fermions, utilizing $\sigma$-matrix algebra and path integral methods to derive the necessary Feynman rules for graviton-fermion and graviton-scalar interactions.
4. Spectrum Evolution: The differential decay rates ($d\Gamma^{(1)}_\phi/dE_l$) are incorporated into Boltzmann equations governing the evolution of the inflaton energy density, the radiation bath, and the graviton energy density during reheating. This allows for the computation of the present-day GW spectrum $\Omega_{GW}h^2$ and characteristic strain $h_c(f)$.

Key Contributions and Results

• Derivation of Enhanced Amplitudes: The primary result is the calculation of the GW spectrum generated by inflaton decay within this SUSY framework. The authors find that the derivative nature of the interactions (arising from the D-term construction) leads to a substantial enhancement in the GW signal compared to standard Yukawa scenarios.
• Quantitative Enhancement: The dimensionless GW strain amplitude is enhanced by at least one order of magnitude compared to typical non-supersymmetric Yukawa coupling scenarios.
  • For the scalar channel, the numerical coefficient $C_{GW}$ is found to be $\approx 43.3 \times 10^{-8}$.
  • For the Majorana fermion channel, the enhancement is significant but dependent on the mass ratio $y = m/M$. For specific parameter choices (e.g., $y=0.1$), the result approaches the scalar enhancement.
• Spectral Characteristics: The paper presents the GW spectrum $\Omega_{GW}h^2$ and strain $h_c(f)$ for two benchmark scenarios involving different inflaton masses ($M$) and reheating temperatures ($T_{rh}$). The signals peak in the high-frequency range ($10^9 - 10^{12}$ Hz), potentially accessible to next-generation high-frequency detectors.
• Comparison with Benchmarks: The predicted signal amplitude ($\Omega_{GW}h^2 \sim 10^{-9}$ to $10^{-12}$ depending on the benchmark) exceeds typical inflaton-decay predictions by over an order of magnitude, placing the signal within the discovery reach of proposed detectors like BBO, DECIGO, and high-frequency resonant cavities.

Significance and Claims
The authors claim that their work demonstrates a direct link between the microscopic structure of supersymmetric physics and macroscopic cosmological observables. Specifically:

1. New Observational Window: The stochastic GW background produced during reheating may serve as a new probe for possible remnants of supersymmetry, particularly if the signal amplitude is found to be significantly enhanced.
2. Role of SUSY Structure: The enhancement is attributed to the specific structure of the chiral supermultiplet and the derivative couplings inherent in the D-term construction. This suggests that the "interplay between reheating-era stochastic GWs and supersymmetric phenomenology" is a fertile ground for further theoretical and phenomenological development.
3. Feasibility: Even within a minimal toy model (the simplest Wess-Zumino model), the calculated enhancements are large enough to be potentially detectable, motivating the extension of these calculations to more realistic SUSY models involving non-trivial Kähler potentials and gauge symmetries.

The paper concludes that the supersymmetric framework significantly improves the theoretical detectability of reheating-era GWs, warranting further exploration of extended Yukawa couplings and higher-derivative superfield structures.