---
title: Gravitational Wave Polarizations in f(R) Gravity
url: https://www.emergentmind.com/topics/gravitational-wave-polarizations-in-metric-f-r-gravity
type: topic
---

# Gravitational Wave Polarizations in f(R) Gravity

Metric $f(R)$ gravity is a class of modifications of General Relativity (GR) in which the Einstein–Hilbert Lagrangian $R$ is replaced by a general function $f(R)$ of the Ricci scalar. Gravitational wave propagation in metric $f(R)$ gravity exhibits essential differences from the pure tensorial wave content of GR, most notably in the emergence of extra polarization modes sourced by the theory's scalar degree of freedom. These supplementary polarizations lead to new experimental signatures and are key discriminants for $f(R)$ gravity versus GR, especially with the era of multi-messenger gravitational wave astronomy. This article presents a comprehensive technical analysis of all aspects governing gravitational wave polarizations in metric $f(R)$ gravity—including field equations, linear mode decomposition, explicit polarization tensors and Newman–Penrose (NP) scalars, propagation speeds, mathematical classification, and observational prospects.

## 1. Linearized Field Equations and Mode Structure

The metric $f(R)$ action is
\[
S = \frac{1}{2\kappa}\int d^4x \sqrt{-g}\,f(R),
\]
with $\kappa=8\pi G$. Varying with respect to $g_{\mu\nu}$ yields the fourth-order field equations,
\[
f'(R)R_{\mu\nu}-\frac{1}{2}f(R)g_{\mu\nu} + (g_{\mu\nu}\Box - \nabla_\mu\nabla_\nu)f'(R) = 0,
\]
where $f'(R)\equiv df/dR$.

In the weak-field limit around a constant-curvature background $R_0$ (often $R_0=0$ for the wave zone), perturbations are expanded as $g_{\mu\nu} = \eta_{\mu\nu} + h_{\mu\nu}$, $|h_{\mu\nu}|\ll1$, and $f(R)$ is Taylor-expanded to second order. The field equations decouple into:
- A transverse-traceless (TT) tensor sector governed by $\Box \bar h_{ij}^{TT}=0$ ($i,j=1,2,3$), describing two massless spin-2 polarizations ($h_+, h_\times$).
- A scalar sector emerging from the trace,
  \[
  3\Box f'(R) + R f'(R) - 2f(R) = 0,
  \]
  linearizing to the Klein–Gordon equation for the scalaron field $\phi$,
  \[
  (\Box - m_s^2)\phi = 0,
  \]
  where
  \[
  m_s^2 = \frac{1}{3}\bigg(\frac{f'_0}{f''_0} - R_0\bigg),
  \]
  with $f'_0 \equiv f'(R_0), f''_0 \equiv f''(R_0)$.

For specific forms, e.g., $f(R) = R^{1+\delta}/R_c^\delta$ ($\delta\ll1$), $m_s^2$ is typically positive and small for small $R_0$, but the scalar can be effectively massive in high-curvature regions, leading to chameleon-type screening [2401.07664, 1901.11277, 1606.07000].

## 2. Polarization Tensors and Physical Effects

The general plane-wave solution for a GW propagating in the $+z$ direction is decomposed as
\[
\bar h_{\mu\nu}(t,z) = h_+(t-z)\,e_{\mu\nu}^{(+)} + h_\times(t-z)\,e_{\mu\nu}^{(\times)} + \phi(\upsilon)\,e_{\mu\nu}^{(s)},
\]
where
- $e_{\mu\nu}^{(+)}$: $e_{xx}= -e_{yy}=1$, all other components zero;
- $e_{\mu\nu}^{(\times)}$: $e_{xy}=e_{yx}=1$;
- $e_{\mu\nu}^{(s)}$: $e_{zz}=1$ (and, in traceless representations, a breathing component $e_{xx}=e_{yy}$).

The scalar field $\phi$ excites both a transverse breathing mode (isotropic deformation in the $xy$ plane) and a longitudinal mode (stretch/compression along $z$) [2401.07664, 1906.01899, 1701.05998, 1608.01764]. The explicit effect on a ring of test particles is:
\begin{align*}
&\text{Tensor $+$}: && \delta x = +\tfrac12 h_+ x_0, \quad \delta y = -\tfrac12 h_+ y_0, \\
&\text{Tensor $\times$}: && \delta x = +\tfrac12 h_\times y_0, \quad \delta y = +\tfrac12 h_\times x_0, \\
&\text{Breathing (scalar)}: && \delta x, \delta y \propto \phi\, x_0, y_0, \\
&\text{Longitudinal (scalar)}: && \delta z \propto \phi\, z_0,
\end{align*}
where the scalar $\phi$ controls both breathing and longitudinal deformations as a single massive degree of freedom [1906.01899].

## 3. Newman–Penrose Scalars and E(2) Classification

The polarization content is rigorously characterized using the NP formalism, which computes projections of the Riemann tensor onto a null tetrad aligned with the GW direction:
- $\Psi_4$: two spin-2 transverse polarizations ($+$ and $\times$);
- $\Phi_{22}$: scalar breathing (transverse scalar);
- $\Psi_2$: scalar longitudinal.

For $f(R)$ gravity, 
\[
\Psi_4\neq0 \quad (\text{tensor}), \quad \Phi_{22}\neq0,~\Psi_2\neq0 \quad (\text{scalar }) 
\]
with $\Psi_{3} = 0$ (no vector polarisations). The non-vanishing of both $\Phi_{22}$ and $\Psi_2$ for a massive scalar mode implies a "mixed" breathing–longitudinal scalar polarization [1711.06220, 0908.0861]. In the $E(2)$ classification, these correspond to class II$_6$, i.e., up to four physical polarizations in principle.

However, the NP/E(2) analysis strictly applies only to null, plane waves. For massive scalar modes ($m_s\neq0$), some NP components may lack physical significance due to subluminal group velocity, and the number of independent dynamical degrees of freedom is three—not four: two TT tensor modes, and a single massive scalar mode whose observable effects project onto both breathing and longitudinal patterns [1701.05998, 1806.04027]. In the massless scalar limit ($m_s \rightarrow 0$), the longitudinal part vanishes, leaving a pure breathing mode [2409.19535, 1901.11277].

## 4. Propagation Speeds, Dispersion Relations, and Mass Dependence

The massless tensor modes satisfy the standard massless dispersion relation,
\[
\omega^2 = k^2, \quad v_g = 1,
\]
and thus propagate at the speed of light in vacuum. The scalar mode propagates according to
\[
\omega^2 = k^2 + m_s^2 \quad \Rightarrow \quad v_g = \frac{k}{\omega} = \sqrt{1-\frac{m_s^2}{\omega^2}},
\]
so its group velocity is strictly subluminal for $m_s>0$, tending to unity for $m_s \ll \omega$ [2401.07664, 1606.07000, 1901.11277, 1906.01899]. Consequently, a frequency-dependent arrival time difference between gravitational wave tensor and scalar components is a robust observable signature for $f(R)$-type modifications.

The Compton wavelength of the scalar mode $\lambda_s = m_s^{-1}$ determines the effective propagation range: for $m_s$ large (short $\lambda_s$), the scalar is strongly suppressed at astrophysical distances ("chameleon" effect). Only for sufficiently small $m_s$ does the scalar contribution persist over interferometric and astrophysical scales [1608.01764, 2006.04011].

In $f(R)$ models embedded as a subclass of Horndeski theory ($G_5=0$, $G_{4,X}=0$), the tensor propagation speed is forced to $c_g=1$ by the structure of the field equations [2401.07664].

## 5. Massless Limit, Stable Massless Scalar, and Model Discrimination

Typically, stability of cosmological perturbations requires $f_{RR}(R_0)>0$ (to avoid Dolgov–Kawasaki instability), so $m_s^2>0$. However, recent work demonstrates that certain $f(R)$ models with higher-order potential inflection points in the scalar sector can allow a strictly massless scalar mode ($m_s=0$) without violating stability, provided $V'(\Phi_d)=V''(\Phi_d)=\ldots=V^{(2K-1)}(\Phi_d)=0$, $V^{(2K)}(\Phi_d)>0$ (for some $K\geq2$) at the de Sitter background [2409.19535]. This yields a pure breathing scalar wave with no longitudinal component. Among popular dark energy $f(R)$ models, only the Gogoi–Dev construction can accommodate this, whereas Hu–Sawicki and Starobinsky types cannot. The detection of a pure breathing mode would thus directly discriminate among viable $f(R)$ models.

| Model          | $m_s=0$ Allowed? | Scalar Observability                                  |
|----------------|:----------------:|------------------------------------------------------|
| Hu–Sawicki     | No               | Only mixed mode if $m_s>0$                           |
| Starobinsky    | No               | Only mixed mode if $m_s>0$                           |
| Gogoi–Dev      | Yes              | Stable pure breathing mode if $m_s=0$                |

## 6. Observational Implications and Prospects

Metric $f(R)$ gravity predicts:
- Always two massless tensor modes ($h_+, h_\times$), as in GR.
- One extra scalar mode, appearing as a mixed massive breathing-longitudinal wave if $m_s >0$, reducing to a pure breathing mode in the $m_s\to 0$ limit [2401.07664, 1906.01899, 1701.05998].

Distinct consequences arise:
- The scalar polarization signal is generically suppressed in high-density environments and in the high-mass regime.
- A frequency-dependent phase and group velocity for the scalar mode leads to arrival time differences, unique angular correlation patterns in interferometers, and distinctive cross-correlations in pulsar timing arrays [1709.03313, 2006.04011].
- Laser interferometers (LIGO, Virgo, KAGRA, LISA) are orders of magnitude less sensitive to longitudinal than breathing polarizations; PTA are much better probes of the longitudinal scalar.
- Bounds on the scalar mass, currently $m_s \lesssim 10^{-22}$ eV, translate into lower limits on $f_{RR}(R_0)$, which constrain the functional form of $f(R)$ [1711.06220, 1606.07000].
- The absence of new polarizations or frequency-dependent dispersion in current GW data can exclude broad regions of $f(R)$ model space; detection of a nonstandard polarization would be a "smoking gun" for modified gravity.

## 7. Synthesis and Theoretical Context

The polarization content of gravitational waves in metric $f(R)$ gravity is fundamentally controlled by the dynamical scalaron degree of freedom, inherited from the higher derivative nature of the theory. The scalar sector modifies gravitational radiation via a single massive spin-0 mode that projects onto, but does not double, the set of observable non-tensor polarizations. Whereas the E(2)/NP formalism classifies up to six polarizations for general metric theories, the actual number of physical GW d.o.f. in metric $f(R)$ is three: two tensor and one scalar, with the scalar manifesting as a mixed (breathing + longitudinal) response unless $m_s=0$, in which case a pure breathing mode remains [2401.07664, 1701.05998, 1608.01764, 2409.19535, 1901.11277, 2006.04011]. Standard GR is recovered as the special case $f(R)=R$, where the scalar mode is infinitely massive and decouples.

Model-dependent predictions for the scalar mass tie GW phenomenology to cosmological and Solar-System constraints, positioning GW polarization measurements as uniquely powerful probes of gravitational physics beyond Einstein [2601.01028]. Future high-precision polarization measurements in both the interferometric and pulsar timing regimes have the potential not only to confirm or falsify $f(R)$ gravity in general, but, via the nature of the scalar mode, to directly distinguish among alternative $f(R)$ models.

**References:**  
- [2401.07664], [1606.07000], [1901.11277], [1906.01899], [1711.06220], [1806.04027], [1709.03313], [2409.19535], [2601.01028], [1912.12469], [1701.05998], [1608.01764], [2006.04011], [0908.0861].

Source: https://www.emergentmind.com/topics/gravitational-wave-polarizations-in-metric-f-r-gravity