Preprint
Article

This version is not peer-reviewed.

Palatini–PT-Even Torsion: Uniqueness, Three-Chain Equivalence, and Exact Luminal Tensors

Submitted:

20 September 2025

Posted:

29 September 2025

Read the latest preprint version here

Abstract
What is new. Within a first–order (Palatini) setting we impose a scalar PT projector on all observable scalar densities, restricting to at most one derivative per building block. Under this setting the framework implies: (i) a unique, pure–trace torsion aligned with ∂ϵ (C1); (ii) a bulk equivalence of three ostensibly different rank–one constructions with the same quadratic coefficient A⋆ = λ2/8 (C2); and (iii) a unique weight ratio between the rank–one determinant and closed–metric routes that removes TT–nonTT mixing and enforces the equal–coefficient identity K=G (C3), hence cT=1 at quadratic order. Setup. All observable scalars are mapped by a scalar PT projector ΠPT( · ) to a real, even sector; the internal phase ϵ(x) is a nondynamical spurion that enters only through ∂ϵ. With global assumptions A1–A6 (PT–invariant domain/measure, [PT , ∗]=0, projection–variation commutation, boundary/topology posture, Nieh–Yan as boundary counterterm, and a trace lock), we work at quadratic order with at most one derivative per building block. The trace lock is enforced algebraically so that Tμ = 3η ∂μϵ (no new canonical pairs under A4). Results. (C1) Palatini algebraicity plus the scalar projector force Sμ=0, qλμν=0, and TA BC = 2η δA [B∂C]ϵ, yielding the invariant IT = −6η2 ΠPT[(∂ϵ)2]. (C2) After C1 and projection, the three routes—(i) a rank–one determinant built from the canonical traceless matrix, (ii) a closed–metric rank–one deformation, and (iii) the PT–even CS/Nieh–Yan shadow—collapse to the same bulk quadratic piece A⋆ √−g (σϵ IT) up to improvements, with A⋆ = λ2/8. (Our “ROD” is a parity–even volume deformation, not Born–Infeld gravity.) (C3) On flat FRW (unitary gauge) the two independent TT–nonTT mixing entries are proportional and fix the locked ratio wROD : wCM = 2 : 3 . On generic admissible patches the two mixing equations are non–collinear, selecting the same unique ratio; with this locking we prove K(w)=G(w) (a total– divergence identity on the admissible variational domain), hence c2 T = 1 and a GR TT sector with no extra propagating degrees of freedom. Reproducibility and scope. All figures and coefficients are reproduced by the accompanying open–source notebooks and scripts; the 2:3 ratio is extracted by an automated mixing–matrix check on flat FRW and its generic uniqueness follows from the non–collinearity argument. Matter couplings respect the posture (axial channel vanishes; the trace coupling is removable up to a Nieh–Yan boundary shift). The leading PT–even NLO correction predicts δc2 T(k) = b k2/Λ2, offering band–limited constraints on Λ. Detailed proofs, improvement currents, and coefficient tables appear in Apps. A–E.
Keywords: 
;  ;  ;  ;  ;  ;  ;  ;  ;  ;  ;  

1. Introduction

Motivation.

Multimessenger bounds enforce luminal gravitational waves near the present epoch, sharply constraining modifications of GR. We ask whether a first–order (Palatini) framework can be formulated that (i) remains strictly P T –even at the level of observables, (ii) closes at quadratic order with at most one derivative per building block, (iii) preserves exact tensor luminality, and (iv) introduces no extra propagating degrees of freedom.

Posture and scope.

We work with independent vierbein e A μ and metric–compatible spin connection ω A B μ , and enforce a scalar  P T projector Π PT ( · ) on all observable scalar densities. An internal phase ϵ ( x ) enters the observed sector only through its gradient and as a spurion (not varied). At quadratic order with at most one derivative per factor, the projected, real scalar
Σ ϵ Π PT ( ϵ ) 2
serves as the unique spurionic invariant; we also use the shorthand I T 1 4 Π PT [ T A B C T A B C ] . Global assumptions (A1–A6)— P T -invariant domain/measure, orientation and [ PT , * ] = 0 , projection/variation commutation, boundary/topology posture, Nieh–Yan as boundary counterterm, and the trace-lock posture—are fixed once and for all (Section 2).

Positioning and falsifiable claims.

In contrast to EC/MAG lines where axial and traceless torsion may propagate in the bulk, our scalar P T posture plus Palatini algebraicity removes both irreps at the analyzed order. The framework is falsified if any of the following fails on admissible patches: (i) the pure–trace map (C1), (ii) the bulk equality of the three rank–one routes with A = λ 2 / 8 (C2), or (iii) the locked equal–coefficient identity K = G (C3), stated as a density total divergence K G = a 3 μ ( a 3 J Δ μ ) and proven in Section 7 (with a weak–field TT representative of J Δ μ recorded in Section 7.1). Paper-level null tests include: (a) recovering A = λ 2 / 8 from the P T -even CS/Nieh–Yan shadow on FRW/Minkowski; (b) reproducing the non-collinearity of the two mixing entries and the locked ratio w ROD : w CM ; and (c) verifying K G = a 3 μ ( a 3 J Δ μ ) with the listed J Δ μ .
  • Main results (at a glance). 
(C1) Palatini– P T  uniqueness under trace-lock (A6). Algebraic torsion consistent with the scalar projector is uniquely pure trace and aligned with ϵ :
T A B C = 2 η δ A [ B C ] ϵ , η > 0 , S μ = 0 , q λ μ ν = 0
and the surviving quadratic invariant is
I T = 6 η 2 Σ ϵ .
(C2) Three–chain equivalence (quadratic, modulo boundary). A rank–one determinant (ROD), a closed–metric rank–one deformation, and the parity–even projected CS/Nieh–Yan piece collapse after projection and C1 to the same bulk invariant:
δ 2 L ROD δ 2 L CM δ 2 L CS + A g I T ( mod μ J μ ) , A = λ 2 8 .
(C3) Coefficient lockingexact luminality and a GR TT sector. Eliminating TT–nonTT mixing fixes a unique weight ratio between ROD and closed–metric routes; the locked quadratic TT block satisfies the equal–coefficient identity
K G = a 3 μ a 3 J Δ μ
(Section 7), hence
c T 2 = 1 , δ 2 S TT = M Pl 2 8 d 4 x a 3 h ˙ i j TT h ˙ i j TT ( k h i j TT ) 2 a 2 , no new DoF .
A TT–gauge weak–field representative of J Δ μ on Minkowski with slowly varying ϵ is provided in Section 7.1, and a covariant FRW representative is listed in Appendix D.

Related work (concise map). 

In EC/MAG lines (e.g., [4,37]), axial and traceless torsion may propagate; by contrast, our scalar– P T Palatini posture plus algebraicity forces pure trace aligned with  ϵ at the analyzed order (C1), removing those irreps in the bulk. Parity–odd dCS gravity (e.g., [14,15]) ties tensor speed to a pseudo–scalar; here the P T  projected sector keeps only even scalars and fixes c T = 1 after locking (C3), modulo improvements. Nieh–Yan terms [9] in our posture serve as boundary counterterms (A5) and do not alter bulk equations. Operationally, the “three–chain’’ collapse (C2) is a parity–even, rank–one statement: a rank–one determinant (ROD), a closed–metric deformation, and the P T –even CS/Nieh–Yan shadow share the same bulk quadratic piece A g sgn ( Σ ϵ ) I T with A = λ 2 / 8 , differing only by improvement currents. We view our results as an even–parity, first–order complement to these traditional lines, tailored to tensor–sector falsifiability.

Falsifiability, stability, and organization. 

Any failure of the pure–trace map (C1), of the three–chain bulk collapse with A = λ 2 / 8 (C2), or of exact luminality after locking (C3) falsifies the framework at quadratic order. Locking is stable on generic patches: the 2 × 2 mixing determinant scales as Δ ( k , Σ ϵ ) = α 0 η 4 Σ ϵ 2 k 2 + O ( k 4 , Σ ϵ 3 ) with α 0 0 , so rank loss occurs only on the measure–zero set { Σ ϵ = 0 } { k = 0 } (Section 6, Appendix D). Beyond leading order, we record an order–of–magnitude data–facing NLO estimate in Section 9 (boxed remark; GW band). The paper is organized as follows: Section 2 fixes A1–A6 and the scalar projector (with an admissible–variation clarification); Section 3 and Section 4 establish closure and prove C1 (including a 3–line block–nondegeneracy proof sketch); Section 5 proves C2 and introduces flux–ratio diagnostics; Section 6 and Section 7 implement locking and prove c T = 1 , providing the weak–field TT representative of J Δ μ (Section 7); matter couplings and NLO/data–facing aspects appear in Section 8 and Section 9. Complete proofs and coefficient tables are collected in Apps. A–F.

2. Global Assumptions & the Scalar PT Projector

We work on an oriented ( 3 + 1 ) -dimensional spacetime with independent vierbein e A μ and metric-compatible spin connection ω A B μ = ω B A μ (Palatini posture). Greek indices μ , ν , are spacetime; Latin A , B , are Lorentz-frame indices with η A B = diag ( , + , + , + ) . Torsion and curvature are
T A D e A = d e A + ω A B e B , R A B d ω A B + ω A C ω C B .
Observable scalar densities are mapped to a real, P T -even sector by a projector Π PT ( · ) defined below. The internal phase ϵ ( x ) is a spurion entering observables only via ϵ and is not varied.
Assumption Posture at a Glance (A1–A6).
A1 
P T -invariant domain/measure. For any scalar density X, X = PT [ X ] . Used in: self-adjointness of Π PT and projector identities (Theorem 1); route-flux diagnostics in Section V.
A2 
Orientation and Hodge dual. The combined P T preserves the chosen orientation and [ PT , * ] = 0 on forms. Used in: parity of Levi–Civita density and commuting with *; boundary accounting.
A3 
Projection commutes with variation. For scalar densities, δ Π PT [ O ] = Π PT [ δ O ] . Used in: Palatini variation and irrep block-diagonalization (Section IV). Checkable proof: Lemma 1 below.
A4 
Topology/boundary posture. Work on trivial patches or impose P T -invariant boundary flux so that improvement currents are pure gauge (no extra canonical pairs). Used in: null tests and flux ratios (Section V); DoF count (Section VII).
A5 
Nieh–Yan as boundary counterterm.  NY d ( e A T A ) affects only boundary conventions; no bulk Euler–Lagrange effect. Used in: three-chain equivalence modulo total derivatives (Section V).
A6 
Trace-lock posture. We enforce T μ = 3 η μ ϵ via an algebraic Lagrange current. Under A4 this introduces no new canonical pairs and only removes the P T -even bilinear independence at quadratic order.

Sufficient conditions for A1. 

On oriented Lorentzian patches with the standard measure d 4 x g and P T -invariant boundaries (compact or AF/FRW fall-offs), the pullback by P T preserves both the integration domain and the measure, hence X = P T [ X ] for scalar densities.

Sufficient conditions for A2. 

If P T preserves the chosen orientation and the metric used for index operations, then P T commutes with the Hodge dual on forms: [ PT , * ] = 0 .

Action of P , T , P T on densities

Let P , T act as standard discrete isometries; T is anti-linear (complex conjugation). On oriented ( 3 + 1 ) -manifolds, the combined P T preserves orientation, so for the Levi–Civita density  ϵ μ ν ρ σ (weight + 1 ) we have
P : ϵ ϵ , T : ϵ ϵ , P T : ϵ + ϵ .
Assumption A2 implies [ PT , * ] = 0 . Consequently, a scalar of the form ϵ μ ν ρ σ X μ ν ρ σ is P T -even iff X is P T -even.

2.1. Scalar P T Projector: Definition And Basic Properties

The combined P T acts anti-linearly on complex-valued densities (complex conjugation accompanies time reversal). On scalar densities we define
Π PT [ O ] 1 2 O + PT [ O ] | real scalar , Π PT 2 = Π PT .
Self-adjointness on real scalars. Using A1 and PT 2 = id ,
Π PT [ O 1 ] Π PT [ O 2 ] = Π PT [ O 1 O 2 ] .
Lemma 1
(Project–vary commutation on scalar densities). Let Π PT [ O ] 1 2 ( O + P T [ O ] ) | real act on complex-valued scalar densities O L 2 ( D ) with A1 (PT-invariant domain/measure) and A4 (fall-offs) imposed. If the spurion ϵ is non-variational and variations act only on ( e , ω ) with compact support or admissible fall-offs, then δ Π PT [ O ] = Π PT [ δ O ] and Π PT is self-adjoint on real scalars.
Proof. 
Anti-linearity of T affects only i i ; variations here are real-linear on ( e , ω ) , and ϵ is fixed. With A1, X = P T [ X ] implies f Π PT [ g ] = Π PT [ f g ] on real scalars (self-adjointness). Then δ Π = 1 2 ( δ O + δ P T [ O ] ) = 1 2 ( δ O + P T [ δ O ] ) = Π [ δ O ] , where boundary terms vanish by A4.    □
δ Π PT [ O ] = Π PT [ δ O ] .

Function-analytic hypotheses and a minimal counterexample.

Assume O L 2 H 1 (scalar densities), ( δ e , δ ω ) H cpt 1 or with A4 fall-offs, and the spurion ϵ is non-variational. Since T is anti-linear (complex conjugation), P T acts as i i . Under these conditions δ is real-linear and δ Π PT [ O ] = Π PT [ δ O ] .
If, however, ϵ is varied, the commutation may fail. Let O = i μ ϵ V μ with a real V μ . Then Π PT [ O ] = 0 (purely P T -odd), but Π PT [ δ O ] = 1 2 ( δ O + P T [ δ O ] ) = δ O 0 when δ ϵ 0 . Hence δ Π PT [ O ] = 0 Π PT [ δ O ] . This shows why we hold ϵ fixed (spurion posture).

Admissible variation space.

Variations act on ( e , ω ) with compact support or with A4 fall-offs: on AF patches ( δ e , δ ω ) H δ 1 with δ < 3 / 2 ; on FRW slices, δ e = O ( r 1 σ ) , δ ω = O ( r 2 σ ) for some σ > 0 . The spurion ϵ is nondynamical (no variation) and T’s anti-linearity only sends i i . Under these conditions δ is real-linear and boundary fluxes vanish, so δ Π PT [ O ] = Π PT [ δ O ] holds term-by-term.

Variational domain for K G identity.

We take variations with compact support on spatial slices or with FRW/AF fall-offs: h i j TT = O ( r 1 σ ) , δ N = O ( r 2 σ ) , i N i = O ( r 2 σ ) with σ > 0 . Then J Δ μ = O ( r 3 σ ) and
δ d 4 x a 3 μ J Δ μ = d Σ μ a 3 δ J Δ μ = 0 .

2.2. P t Quick Tables (Projection-Ready)

Objects.

(Signs indicate intrinsic P/T parities; T is anti-linear.)
Object P T Note
e A μ + + vierbein (metric from its square)
ω A B μ + spin connection
T A μ ν + torsion 2-form components
R A B μ ν + + curvature 2-form components
ϵ (spurion) + enters observables via ϵ only
μ ϵ + flips under P; T anti-linear
ϵ μ ν ρ σ pseudo-density (weight + 1 ); P T : + ; [ PT , * ] = 0 by A2

Scalar monomials and the projector.

Monomial P T Π PT ( · )
( ϵ ) 2 even survives; real by definition
T A B C T A B C even survives; real under projection
ϵ μ ν ρ σ X μ ν ρ σ even X even kept iff X is P T -even; else projected out
T μ μ ϵ even survives pre-lock; reduces to 3 η Σ ϵ after C1/trace-lock
μ ( ) μ (vector density) — (total deriv.) boundary-only under A4; unaffected except taking real part

Notation Quick Reference (Projected, Real Scalars; σ ϵ Scheme)

Σ ϵ denotes the projected, real scalar
Σ ϵ Π PT ( ϵ ) 2
We also set
sgn ( Σ ϵ ) sgn ( Σ ϵ ) , n ^ μ μ ϵ | Σ ϵ | , T μ ν sgn ( Σ ϵ ) n ^ μ n ^ ν 1 4 δ μ ν , τ n ^ μ T μ = 3 η sgn ( Σ ϵ ) | Σ ϵ | .
These choices make orientation flips ( Σ ϵ 0 ) harmless and ensure Tr T = 0 , Tr ( T 2 ) = 3 4 . For later use we adopt the sign-compensated convention
I T 1 4 Π PT T A B C T A B C , sgn ( Σ ϵ ) I T sgn ( Σ ϵ ) I T = 6 η 2 | Σ ϵ | ,
and the key relation (used in Section V)
τ 2 = 9 η 2 | Σ ϵ | = 3 2 sgn ( Σ ϵ ) I T .
Regularity at Σ ϵ = 0 . All rank-one objects are defined on patches with Σ ϵ 0 and extended by continuity via a smooth regulator Σ ϵ ( Σ ϵ 2 + ε 2 ) 1 / 2 with ε 0 + . Quadratic densities and improvement currents remain finite in this limit.

2.3. Selection Rules Under the Scalar Projector

Theorem   1
(Selection rules). Under A1–A5 and [ PT , * ] = 0 (A2), for any complex scalar density O ,
Π PT [ O ] = O ( up to taking the real part ) iff O is P T - even ,
Π PT [ O ] = 0 iff O is P T - odd .
In particular, for admissible tensors X,
Π PT ϵ μ ν ρ σ X μ ν ρ σ = 0 iff X is P T - odd ; Π PT T A B C T A B C R , Π PT ( ϵ ) 2 R , Π PT T μ μ ϵ is P T - even and not annihilated by projection .
Sketch. By definition (2), Π PT averages an object with its P T image and then takes the real part. With A1, X = P T [ X ] ensures self-adjointness (3) on real scalars. A2 implies P T preserves orientation and commutes with *. Since the Levi–Civita density is P T -even while T is anti-linear, ϵ μ ν ρ σ X μ ν ρ σ inherits the P T parity of X, yielding the conditional statement above. Products such as T μ μ ϵ are P T -even (see the quick table), hence survive the projector; their later reduction to Σ ϵ uses the trace-lock/C1 map (Sections III–IV).

Roadmap and Where Each Assumption Enters

The definitions and rules above are the only projectors/boundary tools used later. In particular:
  • Section III (allowed sector/closure): we enumerate all P T -even quadratic monomials with at most one derivative per building block. Pre-lock the basis includes I 1 = Π [ T A B C T A B C ] , I 2 = Σ ϵ , and I 3 = Π [ T μ μ ϵ ] . Post-lock (trace-lock or C1) maps I 3 3 η Σ ϵ .
  • Section IV (C1): A3 enables project-then-vary in the Palatini connection equation; A2/A5 prevent hidden pseudo-scalar contaminations; A4 controls improvements.
  • Section V (C2, σ ϵ scheme): A5 ( NY boundary) and A1/A4 underwrite the equality of route-wise bulk pieces and flux-ratio diagnostics, all written with the sign-compensated invariant sgn ( Σ ϵ ) I T .
  • Sections VI–VII (C3): A1/A4 guarantee that improvement currents do not alter kinetic/gradient coefficients; A3 is used implicitly in all quadratic variations; the equal-coefficient identity is a total divergence on the admissible variational domain.

Boundary/Topology Posture and Fall-Offs (Pointer)

Assumption A4 is realized either by compact P T -invariant domains with vanishing flux or by standard asymptotically flat/FRW fall-offs; the explicit statements and the symplectic-flux check are compiled in Appendix A (boundary notes therein). We exclude torsion defects and multi-valued ϵ patches that would violate the projector’s scalar posture.

Spurion Posture, Ward Identity, and Field Redefinition (Sketch)

Treat ϵ as a spurion transforming as a scalar under diffeos but not varied. Then
δ ξ S = d 4 x g ξ ν μ T μ ν + E ( e , ω ) μ ν · δ ξ ( e , ω ) + μ L ( μ ϵ ) ν ϵ = 0 .
Using the Palatini EL equations E ( e , ω ) = 0 and C1 ( T μ = 3 η μ ϵ ), the spurion current becomes an improvement, hence μ T μ ν = 0 in the bulk (A4). For Dirac matter, the trace coupling T μ J μ is removed by ψ e i α ϵ ψ ; the Jacobian shifts the Nieh–Yan boundary term (A5) but does not affect bulk equations (cf. Section 8).
  • Summary of Section II. We fixed the global posture (A1–A6), defined the scalar P T projector with a checkable commutation lemma (A3), clarified the P T parities of basic objects and of the Levi–Civita density (A2), and stated conditional selection rules: Π PT removes P T -odd scalars and keeps P T -even ones (taking the real part). We emphasized the spurion posture and provided the admissible variational domain that makes the equal–coefficient identity integrable. A compact notation is set in which Σ ϵ is by definition a projected, real scalar and the sign-compensated invariant sgn ( Σ ϵ ) I T σ ϵ I T = 6 η 2 | Σ ϵ | will be used uniformly later. These ingredients feed directly into the Palatini analysis (Section III), the uniqueness theorem (Section IV), and the three-chain equivalence (Section V).

3. Palatini Setting and the Allowed PT-Even Scalar Sector

We work in the Palatini posture with independent vierbein e A μ and spin connection ω A B μ (metric compatible), and we impose the scalar  P T projector Π PT [ · ] on all observable scalar densities (Section 2). Throughout, Assumptions A1–A5 and the selection rules of Theorem 1 are in force. In particular, all statements below are understood after applying Π PT [ · ] , so the projected scalars are real.

3.1. Normalized Spurion Direction and Canonical Rank-One Traceless Tensor

The spurion enters the observed sector only through its gradient. We set, once and for all,
Σ ϵ Π PT ( ϵ ) 2 , sgn ( Σ ϵ ) sgn ( Σ ϵ ) , n ^ μ μ ϵ | Σ ϵ | ,
and define the canonical, dimensionless, traceless rank-one tensor built from n ^ ,
T μ ν sgn ( Σ ϵ ) n ^ μ n ^ ν 1 4 δ μ ν , Tr T = 0 , Tr ( T 2 ) = 3 4 .
The (dimension-one) trace-torsion scale along n ^ is
τ n ^ μ T μ = 3 η sgn ( Σ ϵ ) | Σ ϵ | ,
which will serve as the unique dimensionful spurion scale at this order.1

3.2. What the Projector Allows (Pre- Vs. Post-Lock)

By Theorem 1 and the parity assignments clarified in Section 2, the mixed bilinear T μ μ ϵ is P T -even. Hence it is not removed by the projector: pre-lock it survives as an independent scalar.2 Once the trace lock (introduced in Section 4, Equation (17)) or equivalently the C1 map is enforced, T μ = 3 η μ ϵ and the bilinear collapses to a multiple of Σ ϵ .
For later reference we list the fate of quadratic monomials (“≤ one derivative per building block”):
Monomial P T Projected fate (pre-lock → post-lock)
T A B C T A B C even survives → survives
( ϵ ) 2 even survives → survives
T μ μ ϵ even survives, independent → 3 η Σ ϵ
· ( ) improvement → improvement

3.3. Two-Stage Closure At “One Derivative Per Building Block”

We now state the restricted operator basis used throughout: allow at most one derivative per building block and work to quadratic order in fields. With the projector rules and standard integration-by-parts (IBP), the projected, real  P T -even scalar sector closes in two stages:
Pre - lock basis : I 1 Π PT T A B C T A B C , I 2 Σ ϵ = Π PT [ ( ϵ ) 2 ] , I 3 Π PT [ T μ μ ϵ ] ; Post - lock ( C 1 / trace - lock ) : I 3 T μ = 3 η μ ϵ 3 η Σ ϵ , so the basis reduces to { I 1 , Σ ϵ } .
Lemma 2
(Two-stage closure). Under A1–A5, for any P T -even scalar density O ( 2 ) built from ( T , ϵ ) with at most one derivative per factor,
Π PT [ O ( 2 ) ] = c 1 I 1 + c 2 I 2 + c 3 I 3 + μ J μ ( pre - lock ) ,
for some real constants c 1 , 2 , 3 and an improvement current J μ . After enforcing the trace lock (Equation (17)) or equivalently C1,
Π PT [ O ( 2 ) ] = c 1 I 1 + c ˜ 2 Σ ϵ + μ J ˜ μ ( post - lock ) ,
with c ˜ 2 c 2 + 3 η c 3 and a (possibly shifted) improvement current J ˜ μ .
Sketch. Enumerate quadratic monomials with ≤ one derivative per factor. Theorem 1 eliminates only P T -odd combinations; the three P T -even monomials { T 2 , ( ϵ ) 2 , T · ϵ } span the pre-lock space up to improvements. Self-adjointness of Π PT on real scalars (Equation (3)) justifies moving the projector through products; A4 removes boundary flux from IBP. Imposing T μ = 3 η μ ϵ gives I 3 3 η Σ ϵ , yielding the post-lock reduction.

3.4. Action Skeleton (Pre- And Post-Lock) And Notation For Later Sections

It is convenient to record a minimal bulk skeleton that respects the two-stage closure:
S bulk = d 4 x g M Pl 2 2 R ( e , ω ) + α 1 I 1 + α 2 I 2 + α 3 I 3 ,
where α 1 , 2 , 3 R . After the trace lock (or C1), (11) becomes
S bulk ( post lock ) = d 4 x g M Pl 2 2 R ( e , ω ) + α 1 I 1 + α 2 + 3 η α 3 Σ ϵ ,
so the effect of I 3 is a renormalization of the Σ ϵ coefficient at this order.
For compactness in later sections we will also use the invariant I T 1 4 Π PT T A B C T A B C = 6 η 2 Σ ϵ (which follows after C1 in Section 4), alongside the canonical traceless matrix T and scale τ defined in (8)–(9).

Null Test B (stated; paper-checkable).

On any admissible background with the fall-offs of A4, every projected, P T -even quadratic density with at most one derivative per building block reduces pre-lock to c 1 I 1 + c 2 I 2 + c 3 I 3 up to a total derivative, and post-lock to c 1 I 1 + c ˜ 2 Σ ϵ up to a total derivative. A constructive IBP reduction appears in the Appendix.

Order of Operations (Projection ↔ Variation; Lock vs. Closure)

The scalar projector is applied at the level of densities; by Lemma 1 (Section 2) project-then-vary and vary-then-project commute on scalar densities. The trace lock (Equation (17)) is an algebraic constraint introduced at the level of equations of motion or via a Lagrange current; it is not a P T projection. Accordingly, the correct sequence for building the quadratic sector is:
( i ) enumerate & project pre-lock basis { I 1 , I 2 , I 3 } ( ii ) impose trace lock / C 1 I 3 3 η Σ ϵ .
This sequencing will be used verbatim in Section 4 and Section 5.

4. Uniqueness Theorem (C1)

Under the posture of Section 2 and Section 3 (A1–A5, scalar P T projector, “one derivative per building block”), algebraic torsion consistent with the projector is uniquely pure trace and aligned with the spurion gradient:
T A B C = 2 η δ A [ B C ] ϵ , η > 0 ,
so that the axial and traceless irreps vanish, S μ = 0 and q λ μ ν = 0 . Equivalently,
Π PT T A B C T A B C = 2 3 T μ T μ = 6 η 2 Σ ϵ I T 1 4 Π PT [ T A B C T A B C ] = 6 η 2 Σ ϵ .
Here Σ ϵ Π PT ( ϵ ) 2 is the projected, real scalar, and we use the shorthands from Section 3.1: sgn ( Σ ϵ ) sgn ( Σ ϵ ) , n ^ μ μ ϵ / | Σ ϵ | , T μ ν sgn ( Σ ϵ ) n ^ μ n ^ ν 1 4 δ μ ν , and τ n ^ μ T μ = 3 η sgn ( Σ ϵ ) | Σ ϵ | .
Figure 1. Irreducible torsion content at quadratic order (log–log view). Ratios of projected scalar strengths comparing the pure-trace block against the axial and traceless blocks, shown as trace/axial and trace/tensor on logarithmic axes. The Palatini algebraicity (Section 4.2) and the P T projector drive axial and traceless pieces to zero, leaving the pure-trace map (13) as the unique survivor. [nb: fig_c1_pure_trace.py]
Figure 1. Irreducible torsion content at quadratic order (log–log view). Ratios of projected scalar strengths comparing the pure-trace block against the axial and traceless blocks, shown as trace/axial and trace/tensor on logarithmic axes. The Palatini algebraicity (Section 4.2) and the P T projector drive axial and traceless pieces to zero, leaving the pure-trace map (13) as the unique survivor. [nb: fig_c1_pure_trace.py]
Preprints 177575 g001

4.1. Most General Local Linear Ansatz (One Derivative)

At the derivative order relevant for the quadratic analysis, the only covector available is μ ϵ and the invariant tensors are δ A B and ϵ A B C D . The most general Lorentz-covariant linear ansatz is therefore
T A B C = A δ A [ B C ] ϵ + B ϵ A B C D D ϵ + ( P · ϵ ) A B C ,
with real A , B and where ( P · ϵ ) A B C denotes any attempted traceless mixed-symmetry piece built from a single vector.
  • Proposition B.1 (single-vector no-go; Appendix B). From one vector v μ one cannot construct a nonzero traceless torsion irrep q λ μ ν obeying q ν μ ν = q λ μ μ = 0 and q λ μ ν = q λ ν μ . Any such attempt reduces to the span of δ A [ B v C ] and ϵ A B C D v D .
By Prop. B.1 the last term in (15) is trivial, and the ansatz reduces to
T A B C = A δ A [ B C ] ϵ + B ϵ A B C D D ϵ .
The scalar projector (Section 2) removes all P T -odd scalars, but it does not by itself force B = 0 ; the Palatini equations will.

4.2. Palatini Equations: Algebraic, Irrep Blocks, and Alignment

We augment the pre-lock bulk skeleton (11) with a Lagrange current that enforces alignment of the torsion trace with ϵ ,
L Λ = g Λ μ T μ 3 η μ ϵ , T μ T ν μ ν .
By Lemma 1 (Section 2), project-then-vary and vary-then-project commute on scalar densities. Varying w.r.t. the independent connection then yields algebraic equations in the three irreps { T μ , S μ , q λ μ ν } :
δ ω S = d 4 x g 4 α 1 3 T μ δ T μ + α 1 12 S μ δ S μ + 2 α 1 q λ μ ν δ q λ μ ν + Λ μ δ T μ ,
which is block diagonal because the map from the connection variation to torsion irreps is non-degenerate:
Lemma 3
(Non-degeneracy of δ ω ( δ T μ , δ S μ , δ q λ μ ν ) ). In metric-compatible Palatini, the linear map from δ ω A B μ to the variations of the three torsion irreps is blockwise non-degenerate. Consequently the quadratic form in (18) splits into the three irreps with independent algebraic equations.
  • Proof sketch (3 lines). Varying only the spin connection, δ T A μ ν = δ ω A B [ μ e B ν ] . Projecting to irreps,
    δ T μ δ T ν μ ν = e A ν e B μ δ ω A B ν , δ S μ ϵ μ ν ρ σ δ T ν ρ σ = ϵ μ ν ρ σ e A ν e B ρ δ ω A B σ ,
    and δ q λ μ ν is the traceless remainder of δ T λ μ ν after subtracting the vector/axial projections. These linear maps are surjective and mutually orthogonal with respect to d 4 x g ( 4 3 δ T 2 + 1 12 δ S 2 + 2 δ q 2 ) , hence the quadratic form in (18) splits blockwise and the blocks do not interfere. □
The Euler–Lagrange equations read
( vector ) 4 α 1 3 T μ + Λ μ = 0 ,
( axial ) α 1 12 S μ = 0 S μ = 0 ,
( traceless ) 2 α 1 q λ μ ν = 0 q λ μ ν = 0 ,
and variation w.r.t. Λ μ enforces the lock T μ = 3 η μ ϵ . Substituting back into (16) fixes A = 2 η and forces B = 0 , i.e., alignment with ϵ and the absence of any axial piece, as claimed in (13).

Why the Trace Lock is Not an Extra Assumption

Palatini’s connection equations are algebraic in the three torsion irreps. In the vector block, the Euler–Lagrange equation necessarily constrains T μ to be collinear with the only available covector at this order, μ ϵ . Implementing the collinearity via a Lagrange current Λ μ ( T μ 3 η μ ϵ ) is an equivalent way to enforce this outcome; it introduces no extra canonical pairs under A4 and does not add assumptions beyond the algebraic Palatini dynamics.

On degrees of freedom.

The Lagrange current enforces an algebraic constraint; its Euler–Lagrange equation contains no time derivatives and thus introduces no new canonical pairs. Under A4 this is a first-class enforcement of the trace lock, not an additional propagating sector.

4.3. Positivity, Sign Choice, and the Invariant I T

Using the standard irrep split T A B C T A B C = 2 3 T μ T μ + 1 24 S μ S μ + 2 q λ μ ν q λ μ ν and (20)–(21), we obtain Π PT [ T A B C T A B C ] = 6 η 2 Σ ϵ , which gives (14) for I T 1 4 Π PT [ T A B C T A B C ] . Positivity of the TT kinetic coefficient after locking (Section 5) fixes the physical branch η > 0 .

4.4. Theorem and Three-Step Proof

Theorem   2
(Palatini– P T uniqueness (C1)). Under A1–A5 and the scalar P T projector, algebraic torsion solving the Palatini connection equation with the trace lock (17) is uniquely T A B C = 2 η δ A [ B C ] ϵ with η > 0 . In particular, S μ = 0 and q λ μ ν = 0 , and I T = 6 η 2 Σ ϵ .
Proof (three steps). Step 1 (Ansatz). The most general linear ansatz is (15); by Prop. B.1 (Appendix B) the single-vector traceless attempt vanishes, giving (16).
Step 2 (Palatini blocks). Using Lemma 1 to commute projection with variation and Lemma 3 for blockwise non-degeneracy, varying w.r.t. ω and Λ μ yields (19)–(21) and the lock T μ = 3 η μ ϵ , which fix A = 2 η and B = 0 .
Step 3 (Invariant/positivity). With S = q = 0 and the lock, the invariant reduces to (14). The branch η > 0 follows from positivity of the locked TT sector.    □

4.5. Frw Paper-Level Check And A Geometric Diagnostic

On flat FRW d s 2 = d τ 2 + a 2 ( τ ) d x 2 with homogeneous ϵ ( τ ) , the vector, axial and tensor blocks yield
4 α 1 3 T 0 + Λ 0 = 0 , α 1 12 S μ = 0 , 2 α 1 q λ μ ν = 0 ,
together with the lock T 0 = 3 η ϵ ˙ . Hence S μ = q λ μ ν = 0 and I T = 6 η 2 Σ ϵ , in agreement with (13)–(14).
Figure 2. Alignment-angle diagnostic. Distribution of the alignment angle θ between the torsion trace T μ and μ ϵ , cos θ T μ μ ϵ / T 2 ( ϵ ) 2 . The uniqueness map (13) predicts θ 0 up to finite-domain/boundary remainders. [nb: fig_c1_alignment.py]
Figure 2. Alignment-angle diagnostic. Distribution of the alignment angle θ between the torsion trace T μ and μ ϵ , cos θ T μ μ ϵ / T 2 ( ϵ ) 2 . The uniqueness map (13) predicts θ 0 up to finite-domain/boundary remainders. [nb: fig_c1_alignment.py]
Preprints 177575 g002

Corollaries, Scope, and Order-of-Operations Reminder

Corollary (basis reduction).

With S = q = 0 and T μ = 3 η μ ϵ , the first-order closure basis of Section 3.3 collapses to { I T , Σ ϵ } modulo a total derivative.

Scope and failure modes.

Violations of A2 (orientation/ [ PT , * ] = 0 ) or A4 (boundary/topology) can obstruct projector selection rules or boundary improvements globally; see Appendix B for representation-theoretic caveats and Appendix A for boundary notes. None affect the main (C1) statement on admissible patches.

Order of operations.

The trace lock (17) is an algebraic constraint introduced at the level of equations of motion (or via L Λ ); it is not a P T projection and introduces no new canonical pairs. The correct pipeline, used throughout Sections 5–9, is:
( i ) enumerate & project pre - lock basis { I 1 , I 2 , I 3 } ( ii ) impose trace lock / C 1 I 3 3 η Σ ϵ , use { I T , Σ ϵ } .

5. Three-Chain Equivalence (C2)

We adopt the sign-compensated (“ σ ϵ ”) scheme throughout C2: all rank-one bulk reductions are written in terms of sgn ( Σ ϵ ) I T σ ϵ I T so that patches with Σ ϵ < 0 are handled uniformly, sgn ( Σ ϵ ) I T = 6 η 2 | Σ ϵ | .
We prove that three ostensibly different quadratic routes— (i) a rank-one determinant route built out of the canonical traceless matrix T μ ν , (ii) a closed-metric rank-one deformation, and (iii) the P T -even shadow extracted from the CS/Nieh–Yan chain— collapse, after C1 and projection, to the same bulk invariant A g sgn ( Σ ϵ ) I T up to a total derivative. Differences are carried entirely by improvement currents μ J μ (closed forms on FRW/weak-field are in Appendix C).3

5.1. Preliminaries: Canonical Rank-One Objects And Normalization

Recall Σ ϵ Π PT ( ϵ ) 2 , n ^ μ μ ϵ | Σ ϵ | , and the dimensionless, traceless rank-one matrix
T μ ν sgn ( Σ ϵ ) n ^ μ n ^ ν 1 4 δ μ ν , Tr T = 0 , Tr ( T 2 ) = 3 4 ,
together with the dimension-one spurion scale (after C1),
τ n ^ μ T μ = 3 η sgn ( Σ ϵ ) | Σ ϵ | .
Key relation.   
τ 2 = 9 η 2 | Σ ϵ | = 3 2 ( sgn ( Σ ϵ ) I T ) , sgn ( Σ ϵ ) I T sgn ( Σ ϵ ) I T = 6 η 2 | Σ ϵ | .

5.2. Three One-Line Propositions (NY Split → * & Projection → Coefficient Match)

Proposition 1
(NY split).  T A T A = NY + e A e B R A B with NY d ( e A T A ) .
Proposition 2
( P T -even shadow after * and projection). Under A2, the combined P T preserves orientation and satisfies [ PT , * ] = 0 . Applying * to Prop. 1 and projecting to scalars, the term * NY contributes only as a boundary convention (A5), while the remaining P T -even bulk piece reduces, at quadratic order and after C1, to
δ 2 L CS + = A g sgn ( Σ ϵ ) I T + μ J CS + μ for some A R .
Proposition 3
(Coefficient match under the σ ϵ scheme). With Tr T = 0 , Tr ( T 2 ) = 3 4 , and τ = 3 η sgn ( Σ ϵ ) | Σ ϵ | , both (i) the rank-one determinant route det [ 1 + ( 2 3 ) λ τ T ] and (ii) the closed-metric deformation g ˜ μ ν = g μ ν + ( 2 3 ) λ τ n ^ μ n ^ ν yield the same quadratic bulk coefficient δ 2 L = A g sgn ( Σ ϵ ) I T with A = λ 2 / 8 .
  • Proof (paper-level). Expand det ( 1 + X ) = 1 + 1 2 Tr X + 1 8 [ ( Tr X ) 2 2 Tr X 2 ] + O ( X 3 ) ; with X = ( 2 3 ) λ τ T , Tr T = 0 gives δ 2 det ( 1 + X ) = 1 4 Tr ( X 2 ) = 1 4 ( 4 9 ) λ 2 τ 2 Tr ( T 2 ) = λ 2 τ 2 / 12 . Using (24), δ 2 L = g ( λ 2 / 8 ) sgn ( Σ ϵ ) I T . The closed-metric route has the same Jacobian g ˜ = g det [ 1 + ( 2 3 ) λ τ T ] , hence the identical coefficient.

5.3. Quick Derivations (ROD/CM/CS+)

Rank-one determinant (rank-one determinant route).

L ROD = g det 1 + 2 3 λ τ T 1 δ 2 L ROD = A g sgn ( Σ ϵ ) I T , A = λ 2 / 8 .

Closed-metric route.

g ˜ μ ν = g μ ν + 2 3 λ τ n ^ μ n ^ ν g ˜ = g det [ 1 + 2 3 λ τ T ] , so
δ 2 L CM = A g sgn ( Σ ϵ ) I T , A = λ 2 / 8 .

P T -even CS/Nieh–Yan shadow.

By Prop. 2,
δ 2 L CS + = A g sgn ( Σ ϵ ) I T + μ J CS + μ , A = λ 2 / 8 .

Bulk identity (summary).

δ 2 L ROD δ 2 L CM δ 2 L CS + A g sgn ( Σ ϵ ) I T ( mod μ J μ ) , A = λ 2 / 8 .
Figure 3. Residual scan for the three-chain reduction ( σ ϵ scheme). Quadratic reductions of the ROD/CM/CS+ routes are compared against the target bulk line δ 2 L = bulk A g sgn ( Σ ϵ ) I T with A = λ 2 / 8 . The vertical axis shows the residual after subtracting A g sgn ( Σ ϵ ) I T . All three routes saturate the target within tolerance. [nb: fig_c2_coeff_compare.py]
Figure 3. Residual scan for the three-chain reduction ( σ ϵ scheme). Quadratic reductions of the ROD/CM/CS+ routes are compared against the target bulk line δ 2 L = bulk A g sgn ( Σ ϵ ) I T with A = λ 2 / 8 . The vertical axis shows the residual after subtracting A g sgn ( Σ ϵ ) I T . All three routes saturate the target within tolerance. [nb: fig_c2_coeff_compare.py]
Preprints 177575 g003

5.4. Flux-Ratio Diagnostics and Convergence

For any two routes X , Y { ROD , CM , CS + } define improvement currents by
δ 2 L X A g sgn ( Σ ϵ ) I T = μ J X μ , δ 2 L Y A g sgn ( Σ ϵ ) I T = μ J Y μ .
Assumptions A1–A5 imply
R X / Y [ D ] = 1 ( quadratic order ; posture A 1 - - A 5 ) .
Figure 4. Flux-ratio diagnostic. Boundary flux ratios R X / Y ( R ) for representative pairs X / Y { DBI , CM , CS + } on finite FRW balls converge to 1 as the radius R grows, in agreement with Equation (30). The fit window used to extract the asymptote is annotated as R R min . Error bars reflect the IBP tolerance propagated to boundary terms. [nb: fig_flux_ratio.py]
Figure 4. Flux-ratio diagnostic. Boundary flux ratios R X / Y ( R ) for representative pairs X / Y { DBI , CM , CS + } on finite FRW balls converge to 1 as the radius R grows, in agreement with Equation (30). The fit window used to extract the asymptote is annotated as R R min . Error bars reflect the IBP tolerance propagated to boundary terms. [nb: fig_flux_ratio.py]
Preprints 177575 g004
  • Summary of Section V. At quadratic order and under A1–A5, the rank-one determinant route, closed-metric, and P T -even CS/Nieh–Yan routes share the same bulk coefficient A = λ 2 / 8 multiplying sgn ( Σ ϵ ) I T = 6 η 2 | Σ ϵ | , differing only by improvement currents (Appendix C). Boundary flux ratios R X / Y equal 1 within finite-domain tolerances.

6. Coefficient Locking (C3)

We now lock the relative weight of two bulk–equivalent routes so that the tensor (TT) sector (i) has no kinetic mixing with nonpropagating variables and (ii) is exactly luminal at quadratic order. Throughout we keep A1–A5, the scalar projector Π PT ( · ) , and the C1 pure–trace map T μ = 3 η μ ϵ (Sections 2–4).

6.1. Setup and Locking Posture

By Section 5, the rank–one determinant (ROD) and closed–metric (CM) routes share the same quadratic bulk reduction modulo a total derivative:
δ 2 L ROD δ 2 L CM = A g sgn ( Σ ϵ ) I T ( mod μ J μ ) , A = λ 2 8 .
We therefore consider the linear family
L chain ( w ) w ROD L ROD + w CM L CM , w ROD , w CM R ,
and determine the ratio w ROD : w CM by eliminating TT–nonTT mixing.

6.2. Quadratic ADM Block and Locking Conditions

Expanding L EH ( e , ω ) + L chain ( w ) to second order in ADM variables,
δ 2 L = M Pl 2 8 a 3 K ( w ) h ˙ i j TT h ˙ i j TT G ( w ) ( k h i j TT ) 2 a 2 + a 3 h TT · M ( w ) · Φ + a 3 2 Φ · P · Φ ,
where Φ collects nonpropagating fields (schematically Φ = { δ N , i N i , } ). The mixing block M ( w ) is linear in w and proportional to η 2 Σ ϵ . Two independent entries suffice to enforce (L1) no kinetic mixing:
M ( w ) h TT δ N = μ ROD ( N ) w ROD + μ CM ( N ) w CM η 2 Σ ϵ , M ( w ) h TT i N i = μ ROD ( N ) w ROD + μ CM ( N ) w CM η 2 Σ ϵ ,
with real, dimensionless coefficients μ ( · ) extracted from (i) the contorsion part of R ( e , ω ) under C1 and (ii) the rank–one volume variations. Explicit μ ’s and background dependence are tabulated in Appendix D.
We also impose (L2) exact luminality  K ( w ) = G ( w ) and (L3) GR normalization at ϵ = 0 so that the locked TT action reduces to GR with Planck mass M Pl .

6.3. The 2 × 2 Locking System and Non-Collinearity

Condition (L1) yields the linear system
μ ROD ( N ) μ CM ( N ) μ ROD ( N ) μ CM ( N ) w ROD w CM = 0 0 , ( Σ ϵ 0 ) .
On generic admissible backgrounds the two row vectors are not collinear; the determinant
Δ μ ROD ( N ) μ CM ( N ) μ ROD ( N ) μ CM ( N ) 0 ,
so the solution space is one-dimensional.

6.4. Locked Ratio, Normalization, and Exact Luminality

With Δ 0 there is a unique (up to scale) weight vector w solving (35):
w ROD w CM = μ CM ( N ) μ ROD ( N ) = μ CM ( N ) μ ROD ( N )
(the equality follows from Δ 0 ). The overall scale is fixed by (L3). For this choice, the kinetic and gradient coefficients obey the equal–coefficient identity
K ( w ) G ( w ) = a 3 μ a 3 J Δ μ ,
with a representative J Δ μ listed in Appendix D (TT gauge on weak–field and a covariant FRW form). Since the right-hand side is a total divergence on the admissible variational domain (A4), we obtain
K ( w ) = G ( w ) c T 2 = G K = 1 .
Theorem   3
(Coefficient Locking ⇒ Exact Luminality (C3)). Under A1–A5, the scalar projector, and C1, there exists a unique ratio w ROD : w CM (up to overall normalization) such that TT–nonTT mixing vanishes and K ( w ) G ( w ) = a 3 μ ( a 3 J Δ μ ) . Consequently the TT action is exactly the GR one at quadratic order:
δ 2 S TT = M Pl 2 8 d 4 x a 3 h ˙ i j TT h ˙ i j TT ( k h i j TT ) 2 a 2 , c T 2 = 1 ,
with no additional propagating degrees of freedom.
  • Proof sketch. (i) Palatini algebraicity and the rank–one structure yield the mixing form (34); non-collinearity (36) fixes w up to scale. (ii) For w , the difference K G integrates to a boundary term (38) by A1/A4 (self-adjoint projector; vanishing symplectic flux). (iii) The GR normalization (L3) fixes the overall scale and yields (40).

6.5. Rank Stability And Measure–Zero Degeneracies

Near small k and Σ ϵ ,
Δ ( k , Σ ϵ ; foliation ) = α 0 η 4 Σ ϵ 2 k 2 + O ( k 4 , Σ ϵ 3 ) , α 0 0 ,
so loss of rank occurs only on the measure–zero set { Σ ϵ = 0 } { k = 0 } . This does not affect locking on generic admissible patches.

6.6. Data Companion And Reproducibility (Pointer)

For FRW backgrounds used in figures, the locked ratio w and the two independent mixing entries are extracted directly from the companion data (see figs/data/mixing_matrix.csv with metadata in figs/data/mixing_matrix_meta.json). A convenience configuration mirroring the analytic ratio appears in configs/coeffs/mixing_matrix_FRW.json. These files are diagnostic only; the paper definition of w is Equation (37).
Figure 5. Heatmap of the tensor-speed deviation. Representative scan of | c T 1 | over the ( w ROD , w CM ) plane on an admissible background (A1–A5, C1 in force). The locking curve from (35) is overlaid; along it, TT–nonTT mixing vanishes and (39) holds. [nb: fig_c3_cT_heatmap.py]
Figure 5. Heatmap of the tensor-speed deviation. Representative scan of | c T 1 | over the ( w ROD , w CM ) plane on an admissible background (A1–A5, C1 in force). The locking curve from (35) is overlaid; along it, TT–nonTT mixing vanishes and (39) holds. [nb: fig_c3_cT_heatmap.py]
Preprints 177575 g005
Figure 6. Tensor dispersion c T ( k ) : locked vs. unlocked. Comparison of c T ( k ) for the locked ratio w from (37) (solid) and representative unlocked choices (dashed). [nb: fig_c3_dispersion.py]
Figure 6. Tensor dispersion c T ( k ) : locked vs. unlocked. Comparison of c T ( k ) for the locked ratio w from (37) (solid) and representative unlocked choices (dashed). [nb: fig_c3_dispersion.py]
Preprints 177575 g006
  • Summary of Section VI. Eliminating TT–nonTT mixing reduces to the full–rank system (35). Its nonzero determinant fixes a unique weight ratio w ROD : w CM (up to normalization). With this ratio, the equal–coefficient identity (38) gives K = G and hence c T = 1 .

7. Quadratic Action, Exact Luminality, and Degrees of Freedom

We now state and use the identity that fixes the tensor speed to be exactly luminal at quadratic order once the coefficients are locked, and we record the degrees–of–freedom (DoF) count. We assume A1–A5, the scalar P T projector, C1 with T μ = 3 η μ ϵ , and the locking of Section 6. All projected scalars are real.

Variational domain for the  K G  identity.

We take variations with compact support on spatial slices or with FRW/AF fall-offs:
h i j TT = O ( r 1 σ ) , δ N = O ( r 2 σ ) , i N i = O ( r 2 σ ) ( σ > 0 ) .
Then J Δ μ = O ( r 3 σ ) and δ d 4 x a 3 μ J Δ μ = d Σ μ a 3 δ J Δ μ = 0 , so the total-divergence terms do not feed back into the Euler–Lagrange equations.

7.1. Equal–Coefficient Identity And c T = 1

Expanding L tot = L EH ( e , ω ) + L chain ( w ) to quadratic order (about any admissible background with A4) yields
δ 2 L = M Pl 2 8 a 3 K ( w ) h ˙ i j TT h ˙ i j TT G ( w ) ( k h i j TT ) 2 a 2 + ( TT nonTT mix + constraints ) .
Using the contorsion decomposition of the Palatini curvature under C1, the bulk equivalence of Section 5, and TT transversality/tracelessness, one finds
K ( w ) G ( w ) = a 3 μ a 3 J Δ μ ,
with J Δ μ a quadratic improvement current fixed by the rank–one normalization (Appendix D).
Evaluated at the locked weights w = w (Section 6), the right-hand side is a total divergence. Under A4,
K ( w ) = G ( w ) c T 2 = G K = 1 ( exact at quadratic order ) .
On flat FRW this reduces to K G = a 3 τ ( a 3 J Δ 0 ) .

Weak-Field Check (Minkowski + Slowly Varying ϵ )

On η μ ν with a = 1 and a slowly varying spurion (keep α ϵ = const + O ( 2 ϵ ) ),
J Δ 0 = A 2 η 2 Σ ϵ h i j TT h ˙ i j TT , J Δ k = A 2 η 2 Σ ϵ h i j TT k h i j TT ,
and μ J Δ μ = A 2 η 2 Σ ϵ h ˙ i j TT h ˙ i j TT ( h i j TT ) ( h i j TT ) + μ ν X [ μ ν ] . Under an infinitesimal linearized diffeomorphism, δ ξ J Δ μ = ν X [ μ ν ] [ ξ ] , hence K G = μ J Δ μ is gauge-independent upon integration.

Mode-level check (FRW, finite ball).

For h i j TT ( τ , x ) = { H i j ( τ ) e i k · x } on a ball of radius R, the spatial flux is O ( R 1 ) and the temporal piece does not modify the Euler–Lagrange equation; hence ω 2 = k phys 2 and c T 2 = 1 .

7.2. Locked TT Action (Paper–Level Form)

Imposing the GR normalization (L3) and using (44), the TT action is
δ 2 S TT = M Pl 2 8 d 4 x a 3 h ˙ i j TT h ˙ i j TT ( k h i j TT ) 2 a 2 , c T 2 = 1 .
Figure 7. GW waveform overlay (GR vs. locked). Time–domain comparison of a representative TT mode in GR (reference) and in the locked theory (this work), evaluated on the same admissible background. The shaded band indicates the common numerical tolerance. The root–mean–square error (RMSE) and the best–fit phase offset between the two traces are annotated; both sit within the tolerance when coefficients are locked, consistent with c T = 1 from Equation (44). [nb: fig_gw_waveform_overlay.py]
Figure 7. GW waveform overlay (GR vs. locked). Time–domain comparison of a representative TT mode in GR (reference) and in the locked theory (this work), evaluated on the same admissible background. The shaded band indicates the common numerical tolerance. The root–mean–square error (RMSE) and the best–fit phase offset between the two traces are annotated; both sit within the tolerance when coefficients are locked, consistent with c T = 1 from Equation (44). [nb: fig_gw_waveform_overlay.py]
Preprints 177575 g007

7.3. Degrees Of Freedom (Dirac–Bergmann Count)

The spurion/posture introduces only algebraic Lagrange currents; no new canonical pairs appear on trivial domains with A4 fall–offs. The lapse and shift remain nonpropagating, the Lagrange sector is first class, and no secondary constraints revive extra modes. The propagating content is therefore exactly 2 (the TT tensor), independently of isotropy or slicing: # phys . DoF = 2 at quadratic order. Details and a side–by–side posture comparison appear in Appendix D.
Figure 8. DoF spectrum (eigenvalue stem plot). Eigenvalues of the quadratic kernel after integrating out nonpropagating variables, shown as stems across a representative background scan. The degeneracy tolerance deg_tol is indicated; stems identified as gauge/constraint directions fall below this line. The count of eigenvalues above deg_tol tracks # phys 2 across the scan, confirming the absence of extra propagating modes at quadratic order. [nb: fig_c3_degeneracy.py]
Figure 8. DoF spectrum (eigenvalue stem plot). Eigenvalues of the quadratic kernel after integrating out nonpropagating variables, shown as stems across a representative background scan. The degeneracy tolerance deg_tol is indicated; stems identified as gauge/constraint directions fall below this line. The count of eigenvalues above deg_tol tracks # phys 2 across the scan, confirming the absence of extra propagating modes at quadratic order. [nb: fig_c3_degeneracy.py]
Preprints 177575 g008

8. Coupling to Dirac Matter

We record how minimally coupled Dirac matter interacts with the Palatini– P T posture under A1–A5 and the uniqueness map (C1). Throughout, projected scalars are real (Section 2), and torsion is purely trace and aligned with the spurion gradient,
T μ = 3 η μ ϵ , S μ = 0 , q λ μ ν = 0 ,
as established in Theorem 2. Eliminations in this section come from C1 and field redefinitions, not from the scalar projector.
Two-line summary.
(i)
No axial channel: C1 forces S μ = 0 , so the axial coupling S μ J 5 μ is absent at tree level.
(ii)
Trace channel is removable: with T μ = 3 η μ ϵ , the linear vector coupling T μ J μ is removed by a vector phase redefinition ψ e i α ϵ ψ and reduces to a boundary improvement under A4/A5.

8.1. Setup And Conventions

Consider a Dirac spinor minimally coupled to Riemann–Cartan geometry,
L ψ = e ψ ¯ i e A μ γ A D μ m ψ , D μ ψ μ ψ + 1 4 ω A B μ γ A B ψ ,
with e g , γ A B 1 2 [ γ A , γ B ] , and metric-compatible ω A B μ . Splitting ω = ω ( LC ) + K , the torsion irreps { T μ , S μ , q λ μ ν } couple to J μ ψ ¯ γ μ ψ and J 5 μ ψ ¯ γ μ γ 5 ψ via
L ψ , T = e c A S μ J 5 μ + c V T μ J μ + μ ( ) ,
with real c A , V fixed by conventions (Appendix E). Improvements do not affect bulk Euler–Lagrange equations under A4.

8.2. Axial Channel: Null By C1

Theorem 2 gives S μ = 0 and q λ μ ν = 0 . Hence e c A S μ J 5 μ 0 .

8.3. Trace Channel: Removal By A Local Vector Rephasing

With T μ = 3 η μ ϵ , L ψ , T ( trace ) = 3 c V η e ( μ ϵ ) J μ . Perform ψ e i α ϵ ψ , ψ ¯ ψ ¯ e i α ϵ , which shifts e ψ ¯ i e A μ γ A μ ψ e ψ ¯ i e A μ γ A μ ψ α e ( μ ϵ ) J μ . Choosing α = 3 c V η cancels the trace channel. Up to a divergence, e ( μ ϵ ) J μ = μ ( e ϵ J μ ) ϵ μ ( e J μ ) , and μ ( e J μ ) 0 on the Dirac EOM, so the difference is an improvement fixed by J μ .

Parity remark and measure.

The projector does not remove T μ J μ (it is P T -even); its elimination uses C1 and a vector rephasing. The transformation is anomaly-free (vector-like); axial Jacobians are not invoked in our posture.

With U ( 1 ) .

If ψ carries charge q, the rephasing is equivalent to A μ A μ α q μ ϵ , leaving F μ ν invariant.

NLO Tensor Dispersion and a Band-Limited Estimate

NLO operator and estimate.

At the next order in the “one-derivative-per-building-block” posture, the leading P T -even correction in the TT block is
Δ L NLO = b Λ 2 ( h i j TT ) 2 sgn ( Σ ϵ ) I T δ c T 2 ( k ) = b k 2 Λ 2 .
Taking the multimessenger bound | c T 1 | 10 15 in the LIGO/Virgo band and using c T 1 + 1 2 δ c T 2 gives | b | k 2 / Λ 2 2 × 10 15 . At f = 100 Hz ( k 4.1 × 10 13 eV ),
Λ | b | k 2 × 10 15 9 × 10 6 eV ( length scale 2 cm ) ,
and the bound strengthens by × 10 at f = 1000 Hz . Thus current ground-based bands only weakly constrain k 2 / Λ 2 -type dispersion; lower-frequency probes (e.g., PTA) can improve this if b = O ( 1 ) .
  • Summary of Section VIII. Under the Palatini– P T posture and C1, the axial channel vanishes and the trace channel is removable by a local, anomaly-free vector rephasing, up to a boundary improvement controlled by A4/A5. In addition, the leading P T -even NLO operator yields δ c T 2 = b k 2 / Λ 2 ; a simple LIGO/Virgo estimate places band-limited lower bounds on Λ / | b | .

9. Next-to-Leading Order and Data-facing Remarks

We organize the leading, P T -even corrections beyond the quadratic, one-derivative-per-building-block closure and state their data-facing implication for tensor propagation. Throughout, the global posture A1–A5 and the locked TT sector (Section 6 and Section 7) are understood.

9.1. Minimal NLO Operator and Dispersion

At next-to-leading order (NLO) the unique, parity-even contribution that affects the TT dispersion at quadratic order is a four-derivative tensor operator.4 A convenient parameterization is
Δ L NLO = M Pl 2 8 a 3 b Λ 2 ( 2 h i j TT ) 2 a 4 ,
with a real, dimensionless coefficient b and a heavy scale Λ . For a Fourier mode with physical wavenumber k phys = k / a (today a = 1 ), the quadratic equation of motion gives
ω 2 = k phys 2 1 + b k phys 2 Λ 2 c T 2 ( k ) ω 2 k phys 2 = 1 + b k phys 2 Λ 2 ( k phys Λ ) ,
so that the leading deviation is quadratic in frequency.
Figure 9. NLO offsets and slope fit. Measured tensor-speed offset δ c T 2 ( k ) from the locked LO value c T 2 = 1 as a function of physical wavenumber k (log–log axes). The best-fit slope n ^ 2 is annotated; the vertical offset fixes b / Λ 2 in Equation (54) (a band is shown if multiple backgrounds are included). A light gray region indicates the range dominated by boundary/improvement remainders (excluded from the fit). [nb: fig_nlo_offsets.py]
Figure 9. NLO offsets and slope fit. Measured tensor-speed offset δ c T 2 ( k ) from the locked LO value c T 2 = 1 as a function of physical wavenumber k (log–log axes). The best-fit slope n ^ 2 is annotated; the vertical offset fixes b / Λ 2 in Equation (54) (a band is shown if multiple backgrounds are included). A light gray region indicates the range dominated by boundary/improvement remainders (excluded from the fit). [nb: fig_nlo_offsets.py]
Preprints 177575 g009

9.2. Dimensional Check And Normalization Of b

  • Dimensional check (today a = 1 ). With k = 2 π f the physical wavenumber at a given observed frequency band and Λ a heavy scale, the NLO prediction reads
    δ c T 2 ( k ) = ω 2 k 2 1 = b k 2 Λ 2 ,
    which is manifestly dimensionless. Our normalization of the rank-one tensor T and of τ (Section 5) singles out | b | = O ( 1 )  unless additional heavy operators are tuned to cancel each other at this order.
Remark (NLO estimate at GW band). With δ c T 2 ( k ) = b k 2 / Λ 2 and a representative multimessenger tolerance | c T 1 | 10 15 near f GW 100 Hz ,
| δ c T 2 | 2 | c T 1 | 2 × 10 15 , k 2 π f GW .
Assuming | b | 1 , this gives a lower bound
Λ k | δ c T 2 | 1.4 × 10 10 s 1 9.3 × 10 6 eV Λ c Λ 2.1 cm .
At f GW 1 kHz the bound scales linearly with f, Λ 9.3 × 10 5 eV ( Λ 2.1 mm ). These are order-of-magnitude, band-limited constraints; tighter bounds require a fit to detector bandpasses and b-modeling.

9.3. EFT Validity and Conservative Use of Bounds

Equation (53) is an EFT statement valid for
k phys Λ or equivalently , 2 π f Λ .
In this regime higher-order terms O ( k phys 4 / Λ 4 ) are negligible and the approximation c T 2 ( k ) 1 + b k phys 2 / Λ 2 is self-consistent. When interpreting band-limited constraints, we therefore adopt the conservative rule: only frequencies well below the inferred Λ should be used to quote limits on b / Λ 2 .

9.4. What To Report (Band-Limited, Paper-Level Recipe)

Given a measurement (or bound) on c T ( k ) in a finite band centered at f , report
δ c T 2 ( f ) c T 2 ( f ) 1 = b ( 2 π f ) 2 Λ 2 , k 2 π f , k Λ .
This is the only NLO, P T -even, projector-compatible modification of the locked TT sector at quadratic order that survives as a bulk effect. All other admissible NLO pieces either (i) reshuffle into improvements under A4/A5, or (ii) renormalize K and G equally and thus do not shift c T at this order.

Null Test (NLO).

Across disjoint frequency bands { f i } that all satisfy 2 π f i Λ , the ratio δ c T 2 ( f i ) / δ c T 2 ( f j ) = ( f i / f j ) 2 must hold up to experimental/systematic uncertainties; any systematic departure falsifies the NLO closure embodied by (52).
  • Summary of Section IX. At NLO the locked, P T -even Palatini posture predicts a single, band-limited correction to tensor propagation, δ c T 2 ( k ) = b k 2 / Λ 2 , with | b | = O ( 1 ) under our normalization of T . The EFT is valid for k Λ , under which higher-order terms are suppressed and bounds should be interpreted conservatively.

10. Supplement R: Reproducibility (Lean, Repository-Backed)

All code, figure generators, configs, and tests are open-sourced at:
This supplement gives a minimal, repository-backed map to rebuild the paper figures and to validate C1/C2/C3. We avoid code dumps; the repo already pins the environment (environment.yml), packaging (pyproject.toml), and task runners (Makefile, Snakefile).

Terminology vs. filenames.

The paper uses the term ROD for the rank-one determinant route. Repository filenames keep the legacy tokendbi(e.g., configs/coeffs/dbi.json); they refer to the same route.

R.0 Layout (Pointer)

Top-level directories used by this paper: scripts/ (figure generators), configs/ (grids & coefficient JSON), palatini_pt/ (library), tests/ (pytest), figs/ (outputs and checksum sidecars), and notebooks/ (script mirrors). A one-shot driver scripts/make_all_figs.py rebuilds all paper figures.

R.1 Figure Map (Generator → Artifact)

Generators live in scripts/ and write PDFs to figs/pdf/. Coefficients are in configs/coeffs/{closed,cspp,dbi}.json. Table 1 lists the mapping.
Table 1. Figure file map (repository-backed). 
Table 1. Figure file map (repository-backed). 
ID Generator (scripts/) Inputs (configs/) Output (PDF under figs/pdf/)
Figure 1 fig_c1_pure_trace.py figs/pdf/fig1c1puretrace.pdf
Figure 2 fig_c1_alignment.py figs/pdf/fig2c1alignment.pdf
Figure 3 fig_c2_coeff_compare.py coeffs/dbi.json, coeffs/closed.json, coeffs/cspp.json figs/pdf/fig3c2coeffcompare.pdf
Figure 4 fig_c3_cT_heatmap.py optional: papergrids.yaml figs/pdf/fig4c3cTheatmap.pdf
Figure 5 fig_c3_dispersion.py optional: papergrids.yaml figs/pdf/fig5c3dispersion.pdf
Figure 6 fig_c3_degeneracy.py figs/pdf/fig6c3degeneracy.pdf
Figure 7 fig_gw_waveform_overlay.py figs/pdf/fig7gwwaveformoverlay.pdf
Figure 8 fig_nlo_offsets.py optional: papergrids.yaml figs/pdf/fig8nlooffsets.pdf
Figure 9 fig_flux_ratio.py optional: papergrids.yaml figs/pdf/fig9fluxratio.pdf

R.2 Rebuild & Validate (Three Lines)

  • Environment (conda/mamba).conda env create -f environment.yml;   conda activate palpt
  • Rebuild all figures.python scripts/make_all_figs.py (writes PDFs to figs/pdf/)
  • Validate claims (C1/C2/C3 & diagnostics).pytest -q (covers tests/test_c1_torsion.py, test_c2_equivalence.py, test_c3_tensor.py, test_flux_ratio.py, test_nlo.py, ...)

R.3 Checksums (Sidecars)

Every PDF and data artifact ships a .md5 sidecar (e.g., figs/pdf/fig1c1puretrace.pdf.md5, figs/data/c2residuals.csv.md5). Verification: md5sum -c figs/pdf/*.md5; md5sum -c figs/data/*.md5. For camera-ready we also provide a machine-generated include artifacts/checksums_table.tex; if present, the paper auto-includes it.

R.4 Version Pin

We cite the exact Git revision used to build the artifacts and tag the release. The repository ships figs.tar.gz and notebooks.tar.gz snapshots matching committed artifacts. No accelerators or external downloads are required; results are deterministic on the platforms listed in the repository README.md.
  • Summary. Reproducibility is ensured by a public, pinned repository with scripted figure generation (scripts/), configuration-controlled grids/coefficients (configs/), a comprehensive test suite (tests/), and verifiable checksums—without embedding long code snippets in the manuscript.

11. Related Work & Operational Disambiguation

This section situates our Palatini, scalar- P T -projected posture within the established torsion/gravity lines and records the operational meaning of each global assumption (A1–A5) as used in Sections 2–7. Our goal is twofold: (i) anchor the assumptions in mainstream formalisms and (ii) make the boundary/improvement conventions that underlie the three-chain equivalence and coefficient locking fully checkable.
Table 2. Operational mapping of A1–A5 to mainstream gravitational frameworks (indicative refs).
Table 2. Operational mapping of A1–A5 to mainstream gravitational frameworks (indicative refs).
Assumption Mainstreamanalogue Indicativerefs
A1 measure/domain invariance under discrete symm.; AF/FRW fall-offs Regge–Teitelboim; Wald
A2 [ PT , * ] = 0 orientation-preserving duals; Holst/NY parity split Holst; Shapiro
A3 commute w/ variation Palatini eqs. & self-adjoint projectors in covariant phase space Hehl et al.; Iyer–Wald
A4 boundary posture covariant phase space flux; no extra canonical pairs Iyer–Wald; Regge–Teitelboim
A5 NY boundary NY = d ( e T ) exactness; topology caveats Nieh–Yan; Mercuri; Chandía–Zanelli
1 Remark. Non-Hermitian P T -symmetric QFT is cited only as an analogy for discrete projections on scalars, not as a structural input to gravity here.

11.1. Anchors In The Torsion/Gravity Literature

Our use of independent ( e , ω ) with metric compatibility and the irrep split { T μ , S μ , q λ μ ν } follows the Einstein–Cartan/metric-affine (EC/MAG) tradition and classic reviews [6,8]. The Holst term and its relation to the Nieh–Yan 4-form (exact on appropriate patches) provide the standard parity bookkeeping [9,10,11,12]. Boundary conditions and improvement currents are treated within the covariant phase-space/charge frameworks of Iyer–Wald and Regge–Teitelboim [21,22]. On the parity-odd side, Chern–Simons modified gravity (Jackiw–Pi; Alexander–Yunes) furnishes the canonical reference line [14,15]; our P T posture projects to a parity-even scalar piece accompanied by an improvement current (Section 5), rather than introducing a new propagating pseudo-scalar.

11.2. Operational Disambiguation Of A1–A5

A1 (domain/measure).

X = P T [ X ] ” means: work on oriented Lorentzian patches with the standard measure d 4 x g and P T -invariant boundaries (compact or AF/FRW fall-offs). The pullback by P T preserves both domain and measure for scalar densities, enabling self-adjointness of Π PT on real scalars (Equation (3)) and route-wise flux comparisons (Section 5).

A2 (orientation and Hodge dual).

We assume the combined P T preserves the chosen orientation and the metric used for index gymnastics; hence [ PT , * ] = 0 on forms. This eliminates P T -odd pseudo-scalars by projection (Theorem 1) and fixes parity assignments in the Holst/Nieh–Yan sector.

A3 (project–vary commutation).

For scalar densities and a linear, self-adjoint projector acting as Π PT [ O ] = 1 2 ( O + P T [ O ] ) | real , we use δ Π PT [ O ] = Π PT [ δ O ] (Equation (4)): the projector may be moved through the Palatini variation on ( e , ω ) because the spurion is nondynamical. The formal boundary-aware proof is in Appendix A.

A4 (boundary posture).

We adopt either compact P T -invariant boundaries with vanishing flux or standard AF/FRW fall-offs so that improvement currents contribute no canonical pairs (Iyer–Wald/Regge–Teitelboim). This makes statements like K G = μ ( a 3 J Δ μ a 3 )  operationally slicing/gauge independent at quadratic order (Section 7).

A5 (Nieh–Yan as boundary).

NY d ( e A T A ) is treated as a boundary counterterm on admissible patches: it affects only conventions for improvement currents and does not alter bulk Euler–Lagrange equations. Topological torsion defects or multi-valued ϵ that would spoil exactness are excluded in our posture (see “Exclusions” below).

11.3. “Improvement Current” And Route Equivalence

Throughout, “improvement current” means a total divergence μ J μ whose flux is fixed by A4. The three quadratic routes (DBI, closed metric, parity-even projected CS/Nieh–Yan piece) share the same bulk coefficient A = λ 2 / 8 multiplying I T (Equations (25)–(27)), and differ only by such improvements; on spatially flat FRW a convenient representative is J CS + 0 = λ 2 16 a 3 τ 2 Tr ( T 2 ) , J CS + i = 0 (Equation (45)). The boundary flux ratio diagnostic R X / Y [ D ] = 1 (Equation (30)) implements this equivalence in practice.

11.4. “Closed Metric” Vs. DBI and Rank-One Normalization

The DBI-type determinant det ( 1 + λ τ T ) and the closed-metric rank-one deformation g ˜ μ ν = g μ ν + λ τ n ^ μ n ^ ν have the same quadratic bulk reduction because both expand in traces of the canonical traceless matrix T μ ν with Tr T = 0 and Tr ( T 2 ) = 3 / 4 by construction (Section 5.1). This normalization is the only place where route-dependent choices enter; once fixed, A is universal at quadratic order.

11.5. Exclusions and Scope

We exclude (i) torsion defects or nontrivial bundles that render NY non-exact on the working patch, (ii) multi-valued ϵ or branch cuts, and (iii) boundary conditions that inject new canonical pairs. Within this scope, C1–C3 and the equal-coefficient identity are paper-level, checkable statements.

11.6. Terminology and Naming

To avoid confusion with parity-odd constructions, we refrain from phrases like “PT-even CS/Nieh–Yan shadow.” We instead use parity-even projected piece with improvement current for the scalar contribution that survives Π PT and differs by a total divergence from the DBI/closed-metric routes.

11.7. Reader/Referee Checklist (Operational)

  • Verify A1/A2 on the chosen background (orientation, measure, [ PT , * ] = 0 ).
  • Use the projector rules (Theorem 1) to confirm the first-order closure { I 1 = Π PT [ T 2 ] , I 2 = Σ ϵ } and that Π PT [ T · ϵ ] = 0 .
  • Apply C1 ( T μ = 3 η μ ϵ ) and reproduce A = λ 2 / 8 for DBI/closed-metric/CS+ up to μ J μ .
  • Check the flux ratio R X / Y 1 on growing FRW balls (Section 5.4).
  • Build the 2 × 2 mixing system (Equation (35)), confirm non-collinearity, extract w .
  • Confirm K ( w ) G ( w ) = μ ( a 3 J Δ μ a 3 ) and hence c T = 1 at quadratic order.

11.8. Representative Citations (For Bibliography)

Core anchors for the posture and disambiguation include: EC/MAG and torsion reviews [6,8]; Holst and the Holst–Nieh–Yan relation [9,10,11,12]; covariant phase space and boundary charges [21,22]; and Chern–Simons modified gravity [14,15]. These provide the intended “main gravitational line” against which our scalar- P T Palatini posture is operationally mapped.

Acknowledgments

The author is grateful to the anonymous referees for comments that improved the manuscript. Limited use of generative language tools was made for stylistic refinement; all scientific reasoning, derivations, and conclusions remain solely the responsibility of the author.

Appendix A. Projection–Variation Commutation (A3) and Variational Identities

This appendix establishes Assumption A3 in full generality for scalar densities and compiles the variational identities for g , ϵ μ ν ρ σ , and the Hodge star * under the Palatini posture with the scalar P T projector of Section 2. We also make explicit the boundary/topology posture (A4) used when trading improvements for boundary fluxes, and we record the P T selection rules used throughout the paper.

Appendix A.1. Setup and Conventions

We work with independent vierbein e A μ and a metric-compatible spin connection ω A B μ = ω B A μ (Palatini posture). The observable scalar densities are mapped to a real, P T -even sector by the projector
Π PT [ O ] 1 2 O + PT [ O ] | real scalar , Π PT 2 = Π PT ,
with the combined P T acting anti-linearly (complex conjugation accompanies T) and preserving the chosen orientation (A2), so that [ PT , * ] = 0 on forms. The internal phase ϵ ( x ) is a spurion: it enters observables only via ϵ and is not varied. All statements below are thus variations with respect to ( e , ω , matter ) while keeping ϵ fixed.

Appendix A.2. Projector Properties: Idempotence, Self-Adjointness, and Selection Rules

Proposition A1
(Self-adjointness of Π PT on real scalars). Under A1 (domain/measure P T -invariance) one has, for any scalar densities X , Y ,
Π PT [ X ] Π PT [ Y ] = Π PT [ X Y ] .
Proof. 
Expand the left-hand side using Equation (A1) and the fact that PT 2 = id :
1 4 X + PT [ X ] Y + PT [ Y ] = 1 4 X Y + X PT [ Y ] + PT [ X ] Y + PT [ X ] PT [ Y ] .
Using A1, Z = PT [ Z ] , the cross terms rearrange into 1 2 ( X Y + PT [ X Y ] ) = Π PT [ X Y ] .    □
Proposition A2
(Selection rules for the scalar projector). With A1–A2 and metric compatibility, the following hold for any admissible tensors X:
Π PT ϵ μ ν ρ σ X μ ν ρ σ = 0 , Π PT T μ μ ϵ = 0 ,
Π PT T A B C T A B C R , Π PT ( ϵ ) 2 R .
Proof (sign count). Under A2 the chosen orientation is preserved and [ PT , * ] = 0 ; ϵ μ ν ρ σ is P T -odd, hence any pseudo-scalar density built from it flips sign and is annihilated by Π PT . The spurion gradient μ ϵ picks opposite P / T parities relative to T μ T ν μ ν (table in Section 2.2), so T μ μ ϵ is P T -odd and is projected out. Quadratic contractions T 2 and ( ϵ ) 2 are P T -even and the projector returns their real parts, hence ().    □

Appendix A.3. Commutation of Projection with Variation (A3)

(A3)).Theorem A1 (Projection–variation commutation Let O be any local scalar density built from ( e , ω , ϵ , ) . Then, for variations with respect to ( e , ω , matter ) at fixed ϵ,
δ Π PT [ O ] = Π PT [ δ O ] .
Proof. 
By definition, δ Π PT [ O ] = 1 2 δ O + δ PT [ O ] | real . It suffices to show δ PT [ O ] = PT [ δ O ] . The P T action on fields is an involutive automorphism on the local functional algebra, and it is anti-linear only through global complex conjugation (time reversal). For any complex functional F one has δ F ¯ = δ F ¯ because the variation acts linearly on fields and does not act on the numerical i. Therefore, with Φ denoting collectively the fields that are varied and Φ PT their P T image,
δ PT [ O ( Φ ) ] = δ O ( Φ PT ) ¯ = δ O ( Φ PT ) ¯ = PT [ δ O ( Φ ) ] ,
where we used that PT does not touch the spurion (fixed) and commutes with derivatives and index operations under A2. Substituting back and using linearity of the “real” operation yields Equation (A5).    □

Remarks.

(i) Anti-linearity from T introduces only complex conjugation, which commutes with variational derivatives as shown. (ii) The assumption that ϵ is not varied (spurion posture) is essential; if one promotes ϵ to a dynamical field, additional boundary terms appear but Equation (A5) continues to hold for the scalar projector provided the same posture (A1–A2) is kept for the extended field space.

Appendix A.4. Variational Identities For -g, ϵ μνρσ , And The Hodge Star

We collect formulas used repeatedly in Sections 2–7. We write h μ ν δ g μ ν and h g μ ν h μ ν .

Metric and vierbein.

With g μ ν = η A B e A μ e B ν ,
δ g μ ν = η A B e A μ δ e B ν + e A ν δ e B μ , δ g μ ν = g μ α g ν β δ g α β .

Determinant and Levi–Civita tensor.

δ g = 1 2 g h , δ ϵ μ ν ρ σ = 1 2 ϵ μ ν ρ σ h , δ ϵ μ ν ρ σ = 1 2 ϵ μ ν ρ σ h .
These follow from ϵ μ ν ρ σ = g ε μ ν ρ σ and ϵ μ ν ρ σ = ε μ ν ρ σ / g , with ε μ ν ρ σ the (constant) Levi–Civita symbol.

Hodge star.

Let α be a p-form and h μ ν g μ λ h λ ν . Then the variation of * with respect to h is
δ ( * α ) = * δ α + 1 2 h ( * α ) 1 2 * ( H · α ) , ( H · α ) μ 1 μ p i = 1 p h μ i ν α μ 1 ν μ p .
In particular, for 2-forms F (frequent in the Palatini curvature/torsion algebra),
δ ( * F ) μ ν = ( * δ F ) μ ν + 1 2 h ( * F ) μ ν 1 2 h μ ρ ( * F ) ρ ν + h ν ρ ( * F ) μ ρ .

[PT,*]=0.

Because the metric is P T -even and the chosen orientation is preserved (A2), the Hodge map built from ( g , ϵ μ ν ρ σ ) commutes with P T :
PT [ * α ] = * PT [ α ] for any form α .
This identity is used both in the selection rules and in the projector proofs that involve p-form duals.

Appendix A.5. Boundary/Topology Posture and Improvement Currents

Assumption A4 is realized in either of the following equivalent ways:
(i)
Compact, P T -invariant domains with vanishing boundary flux: for any improvement current J μ arising from integration by parts, D d Σ μ J μ = 0 .
(ii)
Standard fall-offs on asymptotically flat or spatially flat FRW patches, for which d 4 x μ J μ reduces to a surface integral that vanishes in the R limit. A sufficient set is
h μ ν = O ( r 1 δ ) , K λ μ ν = O ( r 2 δ ) , μ ϵ = O ( r δ ) , δ > 0 ,
which ensures J r = O ( r 2 δ ) so that the flux through a sphere of radius R decays as R δ .
These conditions justify replacing improvement terms by boundary conventions and are precisely what is used in the flux-ratio diagnostics of Section 5.4.

Appendix A.6. Consequences Used in the Main Text

(C1) Palatini block-diagonalization.

Theorem A1 (A3) allows us to project then vary in the Palatini equations, so that the connection variation is algebraic and block-diagonal in the torsion irreps. Together with the selection rules (Prop. A2) this yields S μ = q λ μ ν = 0 and the pure-trace map quoted in Section 4.

(C2) Route equivalence modulo boundary.

Self-adjointness (Prop. A1) and the boundary posture (Appendix A.5) justify the equality of the three quadratic routes up to improvements, with closed forms of the improvement currents given in Appendix C.

(C3) Equal-coefficient identity and c T =1.

The star-variation identities (A8)–(A9) are used inside the ADM expansion behind the equal-coefficient identity K G = μ ( a 3 J Δ μ a 3 ) proven in Appendix D. The boundary posture then enforces K = G and c T = 1 at quadratic order.
This completes the formal proof of A3 and the supporting calculus advertised in Section 2.

Appendix B. Irrep Projectors & No-Go for q λμν (v)

This appendix collects the group-theoretic ingredients used in Section 4: (i) the irreducible decomposition of the torsion tensor under the Lorentz group, (ii) explicit, idempotent projectors onto the trace, axial, and traceless sectors, (iii) the quadratic identity for T A B C T A B C in our conventions, and (iv) the single-vector no-go that underlies the statement quoted in the main text as “Proposition B.1” for the q λ μ ν irrep. All statements are purely algebraic and hold before/after applying the scalar projector Π PT [ · ] ; after projection all scalar contractions are real (Section 2).

Appendix B.1. Torsion as a Lorentz Representation and Its Algebra

In index language (spacetime indices), torsion is a rank-3 tensor antisymmetric in its last two indices, T λ μ ν = T λ ν μ , with 4 × 6 = 24 independent components in d = 4 . The Lorentz-covariant irreducible content splits into
24 4 trace T μ 4 axial S μ 16 traceless tensor q λ μ ν , T μ T ν μ ν , S μ ϵ μ ν ρ σ T ν ρ σ , q ν μ ν = q λ μ μ = ϵ μ ν ρ σ q ν ρ σ = 0 .
We use ϵ 0123 = + 1 and the metric signature ( , + , + , + ) , and we adopt the standard scalar product ( X , Y ) X λ μ ν Y λ μ ν on this space.5

Appendix B.2. Idempotent Projectors

Define three linear maps P ( T ) , P ( S ) , P ( q ) on the torsion space by
P ( T ) T λ μ ν 1 3 g λ μ T ν g λ ν T μ ,
P ( S ) T λ μ ν 1 6 ϵ λ μ ν ρ S ρ = 1 6 ϵ λ μ ν ρ ϵ ρ α β γ T α β γ ,
P ( q ) T λ μ ν T λ μ ν P ( T ) T λ μ ν P ( S ) T λ μ ν .
These are the unique Lorentz-covariant, algebraic (derivative-free) projectors onto the three irreps in Equation (A12). A direct computation shows:
P ( X ) 2 = P ( X ) , P ( X ) P ( Y ) = 0 ( X Y ) , P ( T ) + P ( S ) + P ( q ) = 1 ,
and the images obey by construction
P ( T ) T ν μ ν = T μ , ϵ μ ν ρ σ P ( T ) T ν ρ σ = 0 , P ( S ) T ν μ ν = 0 , ϵ μ ν ρ σ P ( S ) T ν ρ σ = S μ , P ( q ) T ν μ ν = 0 , ϵ μ ν ρ σ P ( q ) T ν ρ σ = 0 .

Orthogonality and quadratic split.

With the scalar product ( · , · ) ,
P ( X ) T , P ( Y ) T = 0 ( X Y ) , T , T = P ( T ) T , P ( T ) T + P ( S ) T , P ( S ) T + P ( q ) T , P ( q ) T .
Using (A13)–(A15) one finds the standard quadratic identity
T A B C T A B C = 2 3 T μ T μ + 1 24 S μ S μ + q ˜ λ μ ν q ˜ λ μ ν ,
where we have denoted q ˜ P ( q ) T for the standard normalization of the traceless piece.

Normalization used in the main text.

For later convenience—and to match the coefficient choice used in Section 4—we rescale the traceless irrep by a constant factor and define
q λ μ ν 1 2 q ˜ λ μ ν , T A B C T A B C = 2 3 T μ T μ + 1 24 S μ S μ + 2 q λ μ ν q λ μ ν .
The projector formulas (A13)–(A15) are unchanged; only the bookkeeping name “q” for the traceless image carries the fixed 2 factor.6 After applying Π PT [ · ] to either side of (A20), the scalar is manifestly real (Theorem 1).

Appendix B.3. Compatibility With The Scalar Projector

The projectors P ( X ) are algebraic and commute with Π PT [ · ] at the scalar level: for any two torsions T , T ,
Π PT ( P ( X ) T ) λ μ ν ( P ( Y ) T ) λ μ ν = δ X Y Π PT ( P ( X ) T ) λ μ ν ( P ( X ) T ) λ μ ν .
Moreover, the mixed scalar T μ S μ is P T -odd and is annihilated by Π PT [ · ] (Section 2.3). Thus the orthogonal split (A18) remains valid as an identity between projected, real scalars.

Appendix B.4. Proposition B.1: Single-Vector No-Go For The Traceless Irrep

We now formalize the statement invoked in Section 4.1.
Proposition A3
(single-vector no-go). Let v μ be any nonzero covector. There is no nonvanishing tensor of the form q λ μ ν ( v ) , linear in v μ , that (i) is antisymmetric in μ ν , (ii) obeys q ν μ ν = q λ μ μ = 0 , and (iii) satisfies ϵ μ ν ρ σ q ν ρ σ = 0 . Equivalently, the traceless irrep cannot be constructed from a single vector.
Proof. 
The most general Lorentz-covariant tensor built linearly from a single v μ and antisymmetric in its last two indices is a linear combination of the two rank-one seeds
Y A B C ( v ) = α δ A [ B v C ] + β ϵ A B C D v D ,
with real α , β . Compute its traces and axial contraction:
Y ν μ ν ( v ) = α δ ν [ μ v ν ] = α 2 δ ν μ v ν δ ν ν v μ = 3 α 2 v μ ,
Y λ μ μ ( v ) = α δ λ [ μ v ν ] g μ ν = 0 ,
ϵ μ ν ρ σ Y ν ρ σ ( v ) = β ϵ μ ν ρ σ ϵ ν ρ σ D v D = 3 ! β δ μ D v D = 6 β v μ ,
where we used ϵ μ ν ρ σ ϵ α ν ρ σ = 3 ! δ μ α and d = 4 . Requiring the trace constraints Y ν μ ν = Y λ μ μ = 0 forces α = 0 , and the axial constraint forces β = 0 . Therefore the only admissible linear combination is the trivial one, Y A B C 0 , proving the claim.    □
Corollary A1.
For any single covector v μ the P ( q ) projector annihilates the two rank-one seeds: P ( q ) [ δ A [ B v C ] ] = 0 = P ( q ) [ ϵ A B C D v D ] .

Appendix B.5. Consequences For The C1 Ansatz

Applying Prop. A3 to the most general linear ansatz with one derivative (Section 4.1),
T A B C = A δ A [ B C ] ϵ + B ϵ A B C D D ϵ + P · ϵ A B C ,
shows that the attempted traceless piece P · ϵ A B C necessarily vanishes: it is a linear, single-vector construct and is thus killed by P ( q ) (Cor. A1). The ansatz collapses to
T A B C = A δ A [ B C ] ϵ + B ϵ A B C D D ϵ ,
recovering Equation (4. 2) of the main text. The scalar projector Π PT [ · ] removes P T -odd scalars built from the axial seed (Section 2.3), and the Palatini connection equation then sets B = 0 while fixing A = 2 η under the trace lock T μ = 3 η μ ϵ (Section 4.2). This yields the uniqueness map T A B C = 2 η δ A [ B C ] ϵ quoted in Theorem 2.

Appendix B.6. Consistency Check With The Quadratic Invariant

With the normalization (A20) and the C1 map (so S μ = 0 and q λ μ ν = 0 ),
Π PT T A B C T A B C = 2 3 T μ T μ = 2 3 ( 3 η ) 2 Σ ϵ = 6 η 2 Σ ϵ , I T 1 4 Π PT T A B C T A B C = 6 η 2 Σ ϵ ,
as used throughout Section 4 and Section 5. Here Σ ϵ Π PT ( ϵ ) 2 is the projected, real scalar, and the sign bookkeeping is carried by sgn ( Σ ϵ ) sgn ( Σ ϵ ) ; the unit 1-form n ^ μ μ ϵ / | Σ ϵ | , the canonical traceless rank-one matrix T μ ν sgn ( Σ ϵ ) ( n ^ μ n ^ ν 1 4 δ μ ν ) , and the trace scale τ n ^ μ T μ = 3 η sgn ( Σ ϵ ) | Σ ϵ | are recalled from Section 3.

Appendix B.7. Edge Cases and Patches

On loci where Σ ϵ = 0 the normalized direction n ^ is defined patchwise (or by continuity); all algebraic projector statements remain valid, and the conclusions above hold on any patch with Σ ϵ 0 . Global/topological subtleties (multi-valued ϵ , nontrivial bundles) lie outside the posture A1–A5 (Section 2).
  • Summary of Appendix B. We have given explicit, idempotent projectors onto the three torsion irreps, fixed the quadratic identity in the normalization used in the main text, and proved the single-vector no-go: from one covector v μ no nonzero traceless torsion irrep can be built. This reduces the most general one-derivative ansatz to the { δ A [ B C ] ϵ , ϵ A B C D D ϵ } span, after which the Palatini equations and the scalar projector select the pure-trace map used in the C1 uniqueness theorem.

Appendix C. Appendix C: Three–Chain Reductions & Improvement Currents (σ ϵ Scheme)

  • Scope ( σ ϵ scheme and naming). This appendix provides the paper–checkable reductions behind Section 5 under the sign-compensated convention sgn ( Σ ϵ ) I T sgn ( Σ ϵ ) I T = 6 η 2 | Σ ϵ | and the rank-one determinant route naming (formerly “DBI”-type; not Born–Infeld gravity): (i) a rank-one determinant route built out of the canonical traceless matrix T μ ν , (ii) a closed–metric rank–one deformation, and (iii) the P T –even CS/Nieh–Yan shadow. At quadratic order, each route reduces in the bulk to the same invariant line
    δ 2 L X = A g sgn ( Σ ϵ ) I T + μ J X μ , A = λ 2 / 8 , X { r a n k o n e d e t e r m i n a n t r o u t e , CM , CS + } ,
    with improvements μ J X μ differing by boundary choices (A4/A5). Closed representatives for J μ are given on FRW/weak–field backgrounds.
  • Notation and key relation. We use the preamble shorthands Σ ϵ Π PT ( ϵ ) 2 , sgn ( Σ ϵ ) sgn ( Σ ϵ ) , n ^ μ μ ϵ | Σ ϵ | , T μ ν sgn ( Σ ϵ ) n ^ μ n ^ ν 1 4 δ μ ν , τ n ^ μ T μ = 3 η sgn ( Σ ϵ ) | Σ ϵ | so that Tr T = 0 , Tr ( T 2 ) = 3 4 , and (after C1)
    τ 2 = 9 η 2 | Σ ϵ | = 3 2 ( sgn ( Σ ϵ ) I T ) , sgn ( Σ ϵ ) I T sgn ( Σ ϵ ) I T = 6 η 2 | Σ ϵ | .

Appendix C.1. rank-one determinant route: Determinant Algebra with the 2 3 Normalization

Consider the rank-one determinant route Lagrangian
L ROD = g det 1 + 2 3 λ τ T 1 , λ R .
Using det ( 1 + X ) = 1 + 1 2 Tr X + 1 8 [ ( Tr X ) 2 2 Tr X 2 ] + O ( X 3 ) and Tr T = 0 , the quadratic piece is
δ 2 L ROD = g 1 8 2 Tr ( 2 3 λ τ ) 2 T 2 = 1 4 2 3 2 λ 2 g τ 2 Tr ( T 2 ) = 1 4 · 4 9 · 3 4 λ 2 g τ 2 = 1 12 λ 2 g τ 2 .
With τ 2 = 9 η 2 | Σ ϵ | and sgn ( Σ ϵ ) I T = 6 η 2 | Σ ϵ | ,
1 12 λ 2 τ 2 = 1 12 λ 2 ( 9 η 2 | Σ ϵ | ) = 3 4 λ 2 η 2 | Σ ϵ | = λ 2 8 sgn ( Σ ϵ ) I T ,
so that
δ 2 L ROD = A g sgn ( Σ ϵ ) I T + μ J ROD μ , A = λ 2 8 .
At the bulk-density level one may take J ROD μ 0 . For unified boundary diagnostics we adopt a common canonical representative J can μ for all three routes (Appendix C.4).

Appendix C.2. Closed–Metric Route: Rank–One Deformation Equals rank-one determinant route To O(T 2 )

Take the rank–one deformation g ˜ μ ν = g μ ν + 2 3 λ τ n ^ μ n ^ ν so that g ˜ = g det [ 1 + 2 3 λ τ T ] and define L CM g ˜ g . The same algebra gives
δ 2 L CM = A g sgn ( Σ ϵ ) I T + μ J CM μ , A = λ 2 8 .
Thus rank-one determinant route and CM have the same bulk coefficient A and differ only by improvements.

Appendix C.3. Pt–Even CS/Nieh–Yan Shadow: Quadratic Reduction ( σ ϵ )

Using T A T A = d ( e A T A ) e A e B R A B , applying * and the scalar P T projector (A2), and evaluating after C1, the P T –even piece reduces at quadratic order to
δ 2 L CS + = A g sgn ( Σ ϵ ) I T + μ J CS + μ , A = λ 2 8 ,
where the Nieh–Yan assignment (A5) reshuffles only boundary conventions inside the P T –even sector.

Appendix C.4. A Universal Canonical Improvement at Quadratic Order

For route–by–route flux comparisons it is convenient to select the same improvement representative for all routes:
J can μ [ h TT ] A 2 a 3 h i j TT μ h i j TT δ μ 0 h i j TT h ˙ i j TT ,
so that we adopt the convention
δ 2 L X = A g sgn ( Σ ϵ ) I T + μ J can μ , X { r a n k o n e d e t e r m i n a n t r o u t e , CM , CS + } .
Different representatives differ by ν X [ μ ν ] and yield identical integrated fluxes under A4.

Check (FRW).

On spatially flat FRW in TT gauge,
μ J can μ = A a 3 2 h ˙ i j TT h ˙ i j TT k h i j TT k h i j TT a 2 + τ 3 2 A a 2 a ˙ h i j TT h i j TT ,
i.e., the canonical reshuffling between TT kinetic/gradient bilinears plus a pure time boundary term that integrates to zero with A4 fall–offs.

Appendix C.5. Closed Forms for FRW and Weak Field

FRW.

For d s 2 = a ( τ ) 2 ( d τ 2 + d x 2 ) and homogeneous ϵ ( τ ) ,
J can 0 = A 2 a 3 h i j TT h ˙ i j TT , J can i = A 2 a h j k TT i h j k TT ,
which satisfies (A36).

Weak field (AF).

At a 1 ,
J can 0 = A 2 h i j TT h ˙ i j TT , J can i = A 2 h j k TT i h j k TT .

Appendix C.6. Flux–Ratio Identity & Finite–Domain Convergence

With the unified choice (A35), D d Σ μ J ROD μ = D d Σ μ J CM μ = D d Σ μ J CS + μ , hence
R X / Y [ D ] D d Σ μ J X μ D d Σ μ J Y μ = 1 , X , Y { r a n k o n e d e t e r m i n a n t r o u t e , CM , CS + } .
On finite FRW balls (or AF shells) residuals scale away with the radius R, agreeing with Section V.

Summary of Appendix C

At quadratic order and under A1–A5 plus C1, the rank-one determinant route, closed–metric, and P T –even CS/Nieh–Yan routes share the same bulk reduction δ 2 L X = A g sgn ( Σ ϵ ) I T ( mod μ J μ ) , A = λ 2 / 8 . A single canonical improvement J can μ (Equations (A34)–(A36)) is used for all routes and underlies the flux–ratio plots in Section V. [ nb : fig _ c 2 _ coeff _ compare . py ; fig _ flux _ ratio . py ]

Appendix D. Appendix D: Mixing Matrix and the Equal–Coefficient Identity

This appendix contains (i) extraction rules and tables for the 2 × 2 mixing matrix used in Section 6, including a non-collinearity proof of its two row vectors on admissible backgrounds, and (ii) a covariant derivation of the equal–coefficient identity  K ( w ) G ( w ) = a 3 μ ( a 3 J Δ μ ) quoted in Section 7. We assume A1–A5, the scalar P T projector (Section 2), and the C1 map T μ = 3 η μ ϵ . Projected scalars are real by construction; the σ ϵ scheme only affects the bulk line through sgn ( Σ ϵ ) I T = sgn ( Σ ϵ ) I T = 6 η 2 | Σ ϵ | , not the kinematical identity K G = J .

Variational domain (used below).

We take variations with compact support on spatial slices or with FRW/AF fall-offs: h i j TT = O ( r 1 σ ) , δ N = O ( r 2 σ ) , i N i = O ( r 2 σ ) with σ > 0 . Then J Δ μ = O ( r 3 σ ) and δ d 4 x a 3 μ J Δ μ = d Σ μ a 3 δ J Δ μ = 0 .

Appendix D.1. Adm Conventions and Extraction of Mixing Entries

With d s 2 = N 2 d τ 2 + γ i j ( d x i + N i d τ ) ( d x j + N j d τ ) , N = 1 + δ N , N i = i χ + N i T ( i N i T = 0 ), and γ i j = a 2 ( τ ) ( δ i j + h i j ) , we project h i j h i j TT with i h i j TT = 0 = δ i j h i j TT . The quadratic Lagrangian takes the block form
δ 2 L = M Pl 2 8 a 3 K ( w ) h ˙ i j TT h ˙ i j TT G ( w ) ( k h i j TT ) 2 a 2 + a 3 h TT · M ( w ) · Φ + a 3 2 Φ · P · Φ ,
where Φ = { δ N , i N i , } collects nonpropagating pieces. Define the dimensionless mixing entries by
M ( w ) h TT δ N 8 M Pl 2 1 a 3 ( δ N ) δ 2 L | bilinear in h TT , δ N ,
M ( w ) h TT i N i 8 M Pl 2 1 a 3 ( i N i ) δ 2 L | bilinear in h TT , i N i .
For L chain ( w ) = w ROD L ROD + w CM L CM , the mixing block is linear in w and proportional to η 2 Σ ϵ :
M ( w ) h TT δ N = μ ROD ( N ) w ROD + μ CM ( N ) w CM η 2 Σ ϵ , M ( w ) h TT i N i = μ ROD ( N ) w ROD + μ CM ( N ) w CM η 2 Σ ϵ ,
defining the four dimensionless coefficients μ ROD ( N ) , μ CM ( N ) , μ ROD ( N ) , μ CM ( N ) .

Appendix D.2. Background Invariants And Compact Parametrization

Introduce
H a a , Υ τ ln τ / a = τ τ H ,
(prime is conformal-time derivative). Each μ admits a linear decomposition
μ X ( Y ) = c X , 0 ( Y ) + c X , H ( Y ) H + c X , Υ ( Y ) Υ , X { ROD , CM } , Y { N , N } ,
with c’s real O ( 1 ) numbers fixed by the quadratic expansion rules (route Jacobians plus contorsion under C1).

Appendix D.3. Coefficient Tables (FRW and Weak Field)

Table A1. FRW coefficients with ϵ = ϵ ( τ ) . Entries are the dimensionless μ ’s of (A43), written as μ = c 0 + c H H + c Υ Υ [Equation (A45)]. Overall factor η 2 Σ ϵ multiplies the mixing (not shown here).
Table A1. FRW coefficients with ϵ = ϵ ( τ ) . Entries are the dimensionless μ ’s of (A43), written as μ = c 0 + c H H + c Υ Υ [Equation (A45)]. Overall factor η 2 Σ ϵ multiplies the mixing (not shown here).
Coefficient c 0 c H c Υ
μ ROD ( N ) c ROD , 0 ( N ) c ROD , H ( N ) c ROD , Υ ( N )
μ ROD ( N ) c ROD , 0 ( N ) c ROD , H ( N ) c ROD , Υ ( N )
μ CM ( N ) c CM , 0 ( N ) c CM , H ( N ) c CM , Υ ( N )
μ CM ( N ) c CM , 0 ( N ) c CM , H ( N ) c CM , Υ ( N )
  Regeneration recipe (FRW). (1) Write δ g μ ν in terms of ( δ N , i N i , h i j TT ) ; (2) expand T , τ to linear order using n ^ μ = μ ϵ / | Σ ϵ | and C1; (3) insert into L ROD and L CM with the 2 3 scaling; (4) keep only h TT δ N and h TT i N i bilinears to read off c 0 , c H , c Υ .
Table A2. Weak field (AF) coefficients with slowly varying ϵ ( t , x ) . Set a 1 ( H = 0 ) and define Υ t ln τ . Overall factor η 2 Σ ϵ multiplies the mixing (not shown here).
Table A2. Weak field (AF) coefficients with slowly varying ϵ ( t , x ) . Set a 1 ( H = 0 ) and define Υ t ln τ . Overall factor η 2 Σ ϵ multiplies the mixing (not shown here).
Coefficient c 0 c Υ
μ ROD ( N ) c ˜ ROD , 0 ( N ) c ˜ ROD , Υ ( N )
μ ROD ( N ) c ˜ ROD , 0 ( N ) c ˜ ROD , Υ ( N )
μ CM ( N ) c ˜ CM , 0 ( N ) c ˜ CM , Υ ( N )
μ CM ( N ) c ˜ CM , 0 ( N ) c ˜ CM , Υ ( N )
  Regeneration recipe (weak field). Repeat FRW steps with a 1 , replace conformal by physical time, and expand consistently in spatial gradients so that Σ ϵ remains the projected, real scalar.

Closed-form labels (for Section VI cross-references).

For later citation we record the compact forms as Equations (D.12)–(D.15):
μ ROD ( N ) = c ROD , 0 ( N ) + c ROD , H ( N ) H + c ROD , Υ ( N ) Υ ,
μ CM ( N ) = c CM , 0 ( N ) + c CM , H ( N ) H + c CM , Υ ( N ) Υ ,
μ ROD ( N ) = c ROD , 0 ( N ) + c ROD , H ( N ) H + c ROD , Υ ( N ) Υ ,
μ CM ( N ) = c CM , 0 ( N ) + c CM , H ( N ) H + c CM , Υ ( N ) Υ .

Appendix D.4. Non-Collinearity of the Two Locking Equations

Let the two rows be r 1 = ( μ ROD ( N ) , μ CM ( N ) ) , r 2 = ( μ ROD ( N ) , μ CM ( N ) ) .
Proposition A4
(Non-collinearity). Under A1–A5, C1, and the one-derivative-per-building-block posture, r 1 ¬ r 2 on any admissible background with at least one of { H , Υ } 0 and μ n ^ ν 0 . Equivalently,
Δ μ ROD ( N ) μ CM ( N ) μ ROD ( N ) μ CM ( N ) 0 .
  • Proof (sketch). Using (A45), proportionality would require ρ X μ X ( N ) / μ X ( N ) to be route–independent and constant under independent shifts of H and Υ . But H ρ X and Υ ρ X cannot both vanish for X = ROD , CM unless c X , H ( Y ) = c X , Υ ( Y ) = 0 (static trivial branch H = Υ = 0 , n ^ = 0 ). Hence r 1 ¬ r 2 generically. □

Small-k structure (symbolic).

On weak-field patches one finds
Δ ( k , Σ ϵ ) = α 0 η 4 Σ ϵ 2 k 2 + O ( k 4 , Σ ϵ 3 ) , α 0 0 ,
so rank loss occurs only on the measure-zero set { k = 0 } { Σ ϵ = 0 } or special foliations.

Appendix D.5. Equal–Coefficient Identity: Covariant Derivation And Gauge Shift

Define the quadratic current (indices contracted with the background spatial metric)
J Δ μ = M Pl 2 8 J Δ 0 = Π i j h i j TT J Δ k = Ξ i j k h i j TT ,
with Π i j δ ( δ 2 L ) / δ h ˙ i j TT and Ξ i j δ ( δ 2 L ) / δ ( 2 h i j TT ) . Using TT conditions, the rank–one traceless normalization, and the bulk equality of the routes (Section V),
K ( w ) G ( w ) = 1 a 3 μ a 3 J Δ μ .
Under δ ξ h i j TT = 0 , δ ξ h 0 μ = μ ξ 0 + , the representative shifts as δ ξ J Δ μ = ν X [ μ ν ] [ ξ ] , hence the integrated identity is gauge independent.

FRW and weak-field representatives.

On spatially flat FRW ( d s 2 = a 2 ( d τ 2 + d x 2 ) , homogeneous ϵ ),
J Δ 0 = A η 2 2 Σ ϵ h i j TT h ˙ i j TT + i ( ) ,
J Δ k = A η 2 2 Σ ϵ 1 a 2 h i j TT k h i j TT + ( ) [ k ] ,
so K G = a 3 τ ( a 3 J Δ 0 ) + a 3 k ( a 3 J Δ k ) is a total divergence. In weak field (AF, a 1 ),
J Δ 0 = A η 2 2 Σ ϵ h i j TT h ˙ i j TT + i ( ) i , J Δ k = A η 2 2 Σ ϵ h i j TT k h i j TT + ( ) [ k ] .

Boundary consequence and luminality.

With the variational domain above (A4), a 3 ( K G ) = μ ( a 3 J Δ μ ) = 0 . At the locked weights w (Section VI; no TT–nonTT mixing), K ( w ) = G ( w ) c T 2 = 1 at quadratic order, slicing–independently.

Appendix D.6. Locked Weights and Summary Box

Non-collinearity with the two mixing equations yields a unique ratio w ROD / w CM (up to the GR normalization). At these weights, K G = a 3 μ ( a 3 J Δ μ ) implies exact luminality.
Appendix D (at a glance).
  • Mixing matrix. Four dimensionless coefficients μ X ( Y ) ( X = ROD , CM ; Y = N , N ) control the TT–constraint mixing with an overall η 2 Σ ϵ factor.
  • Non-collinearity. The two locking rows are not proportional on evolving admissible backgrounds; the determinant Δ is generically nonzero ( Δ α 0 η 4 Σ ϵ 2 k 2 + ).
  • Equal–coefficient identity.  K G = a 3 μ ( a 3 J Δ μ ) ; under A4 this yields K = G and c T 2 = 1 at the locked weights w .

Reproducibility Note

Table A1Table A2 come from a single symbolic pipeline that: (i) expands T and τ to linear order in ADM variables with the 2 3 normalization, (ii) forms the rank-one determinant route/CM quadratic densities, (iii) projects onto the bilinears of Equations (A41)–(A42). Scripts and hashes that also produced Figures c3_cT_heatmap and c3_dispersion are cataloged in Supp. R. [ nb : mixing _ matrix _ extract . py ]

Appendix E. Dirac Coupling, Field Redefinition, and the Nieh–Yan Scope

This appendix supplies the details promised in Section 8: (i) the torsion–fermion couplings in our conventions, (ii) the precise field redefinition that removes the trace channel once the pure-trace map (C1) is imposed, (iii) the role of the Nieh–Yan density as a boundary convention (A5), and (iv) the scope limitations on nontrivial topology (A4).
Throughout we work under the global posture A1–A5 of Section 2, use metric signature η A B = diag ( , + , + , + ) , and adopt the gamma-matrix conventions
{ γ A , γ B } = 2 η A B , γ 5 i γ 0 γ 1 γ 2 γ 3 , ϵ 0123 = + 1 ,
so that γ A B 1 2 [ γ A , γ B ] and γ A B C γ [ A γ B γ C ] = i ϵ A B C D γ 5 γ D .

Appendix E.1. Minimal Dirac Coupling In Riemann–Cartan Space

The minimally coupled Dirac Lagrangian is
L ψ = e ψ ¯ i e A μ γ A D μ m ψ , D μ ψ = μ ψ + 1 4 ω A B μ γ A B ψ , e g .
Splitting the spin connection into Levi–Civita plus contorsion, ω A B μ = ω A B μ ( LC ) + K A B μ , the torsion-dependent piece of (A51) is
Δ L ψ [ K ] = i 4 e ψ ¯ e A μ γ A K B C μ γ B C ψ = i 4 e ψ ¯ γ A K B C A γ B C ψ ,
where tangent indices are moved with η A B , and K A B C e A μ K B C μ .
It is convenient to decompose contorsion into the standard torsion irreps (vector trace T A , axial S A , and traceless q A B C ):
K A B C = 1 3 η A C T B η A B T C 1 6 ϵ A B C D S D + q A B C , T A T B A B , S A ϵ A B C D T B C D ,
with q B A B = 0 and ϵ A B C D q B C D = 0 . Using the gamma identities listed below Equation (A50) and the antisymmetry of K B C A in ( B , C ) , Equation (A52) reduces to
Δ L ψ [ K ] = e c A S μ J 5 μ + c V T μ J μ + e C μ ν ρ ( q ) ψ ¯ γ [ μ γ ν γ ρ ] ψ ,
where J μ ψ ¯ γ μ ψ and J 5 μ ψ ¯ γ μ γ 5 ψ . In our normalization (A50) the numerical coefficients are
c A = 3 8 , c V = 1 4
while the q-channel is proportional to the totally antisymmetric part of q, which vanishes by definition; equivalently, the q-coupling can be re-expressed in terms of S μ and drops out when S μ = 0 . Therefore the torsion-induced Dirac channels in our conventions are precisely
Δ L ψ [ T , S ] = e 3 8 S μ J 5 μ 1 4 T μ J μ ,
up to P T -even improvement terms annihilated by A1/A4 in the bulk.

Cross-check.

Equation (A56) reproduces the two-channel bookkeeping used in Section 8 with c A = 3 8 and c V = 1 4 . Any alternative gamma/Hodge convention simply rescales (A55) by an overall sign; our later conclusions (vanishing axial channel by C1 and removability of the trace channel) are insensitive to such simultaneous flips.

Appendix E.2. C1 Implies no Axial Channel; The Trace Channel is Removable

By the uniqueness theorem (C1), Equation (13), the axial and traceless torsion irreps vanish, S μ = 0 , q λ μ ν = 0 , and the trace aligns with the spurion gradient, T μ = 3 η μ ϵ . Equation (A56) therefore reduces to the single trace channel
L ψ , T ( trace ) = e c V T μ J μ = 3 c V η e ( μ ϵ ) J μ .
Proposition A5
(Vector rephasing removes the trace channel). Consider the local vector phase redefinition ψ ψ = e i α ϵ ψ , ψ ¯ ψ ¯ = ψ ¯ e i α ϵ with a constant α. Then
L ψ [ ψ , ψ ¯ ] = L ψ [ ψ , ψ ¯ ] α e ( μ ϵ ) J μ + μ ,
so choosing α = 3 c V η ( = 3 4 η in our normalization ) cancels (A57) pointwise. Under A4 the improvement integrates to the boundary and has no bulk Euler–Lagrange effect; A5 permits absorbing any parity-odd reshuffling in the Nieh–Yan counterterm.
Sketch 
Insert ψ = e i α ϵ ψ into the kinetic term e ψ ¯ i e A μ γ A μ ψ and use the Leibniz rule. The derivative acting on e i α ϵ generates the shift α e ( μ ϵ ) J μ in (A58); the mass and spin-connection pieces are invariant under a vector (not axial) phase. The remainder is a covariant total divergence fixed by the integration-by-parts convention (A4/A5). □
Combining (A57) with Prop. A5 and the coefficient assignment (A55) yields precisely the two-line summary in Section 8: no axial channel and a removable trace channel.

Path-integral measure (anomaly) check.

The field redefinition in Prop. A5 is a vector  U ( 1 ) rotation, hence the fermionic measure is invariant (Jacobian equal to one). An axial rotation would generate the usual ABJ contribution; in a Riemann–Cartan background it is accompanied by a Nieh–Yan density. We do not perform an axial rotation anywhere in this work.

Appendix E.3. Nieh–Yan as a Boundary Convention (A5)

The Nieh–Yan 4-form NY d ( e A T A ) is an exact form. In our posture (A5), any explicit use of NY enters only as a boundary convention after applying the scalar projector Π PT : it does not modify Euler–Lagrange equations in the bulk and is indistinguishable from a choice of improvement current. This applies equally to (i) the three-chain equivalence (where the P T -even CS/Nieh–Yan shadow differs from the DBI/CM routes by an improvement current) and (ii) the Dirac sector (where vector rephasing reshuffles boundary terms that can be absorbed into the chosen Nieh–Yan convention). No physical statement in Sections 5–8 depends on a specific Nieh–Yan choice.

Appendix E.4. Scope and Caveats: Topology and Boundary Flux

Our conclusions rely on the posture A4: either compact P T -invariant domains with vanishing boundary flux or standard asymptotic fall-offs (FRW/flat) so that improvement currents integrate to zero. They also assume a single-valued spurion phase ϵ ( x ) so that the rephasing in Prop. A5 is a globally well-defined map ψ e i α ϵ ψ .
  • Nontrivial topology (not covered). If ϵ is multi-valued or admits nontrivial holonomy (e.g., on manifolds with nontrivial π 1 or H 1 ), the map e i α ϵ may fail to be single-valued. In such cases the local cancellation in (A58) can leave a global remnant proportional to the winding; our bulk equivalences and removability statements are not asserted in that setting.
  • Nonvanishing boundary flux (not covered). On domains where the P T -invariant boundary flux of the relevant improvement currents does not vanish, boundary terms may carry physical information (e.g., in explicitly finite boxes with prescribed inflow). Our null tests and cancellations are presented only under A4 fall-offs.

Appendix E.5. Bookkeeping Table And Quick References

For convenience we summarize the conventions and coefficients used in the Dirac sector:
Object / Convention Value / Definition
Gamma matrices { γ A , γ B } = 2 η A B , γ 5 = i γ 0 γ 1 γ 2 γ 3
Hodge / Levi-Civita ϵ 0123 = + 1 , [ PT , * ] = 0 (A2)
Contorsion split K A B C = 1 3 ( η A C T B η A B T C ) 1 6 ϵ A B C D S D + q A B C
Dirac currents J μ = ψ ¯ γ μ ψ , J 5 μ = ψ ¯ γ μ γ 5 ψ
Torsion→Dirac Δ L ψ = e 3 8 S μ J 5 μ 1 4 T μ J μ
C1 (pure trace) S μ = 0 , q λ μ ν = 0 , T μ = 3 η μ ϵ
Vector rephasing ψ e i α ϵ ψ with α = 3 c V η = 3 4 η
Outcome Trace channel removed up to a total derivative (A4/A5)
Nieh–Yan Boundary convention only (A5); no bulk Euler–Lagrange effect

Appendix E.6. Corollary: No LO Four-Fermion Contact from Torsion

At the order analyzed in this paper (quadratic in fields; at most one derivative per building block), torsion is algebraic and fixed by (C1), and the only linear Dirac–torsion channel that survives projection is removed by the vector redefinition above. Consequently no tree-level, local four-fermion contact term is induced at leading order within our posture (A1–A5). Any such effect would require (i) an axial channel ( S μ 0 ), (ii) higher-derivative completions beyond our closure basis, or (iii) loop corrections outside the present scope.
  • Summary of Appendix E. In our conventions the minimally coupled Dirac field interacts with torsion through e ( 3 8 S μ J 5 μ 1 4 T μ J μ ) . Under the pure-trace map (C1) the axial channel vanishes and the remaining trace channel is removed by a local vector rephasing ψ e i ( 3 c V η ) ϵ ψ , up to a total derivative consistent with A4/A5. Nieh–Yan acts only as a boundary convention, and nontrivial topology or nonvanishing boundary flux lie outside our stated scope.

Appendix F. Operational Diagnostics (Extended)

This appendix collects a practitioner-oriented checklist and an expanded comparison table. Unless stated otherwise, all entries are framed for admissible patches with the boundary/topology posture A4, and for the quadratic order analyzed in the main text.

Appendix F.1. Full Decision Tree (D1–D8)

  • (D1) Parity channels. Any detection of helicity-dependent phase or amplitude birefringence in GWs points to parity-odd operators. Our scalar- P T posture forbids such effects at LO (Theorem 1; Section 2.2).7
  • (D2) Axial/tensor torsion irreps. Nonzero axial ( S μ ) or traceless ( q λ μ ν ) torsion signals are incompatible with C1. We predict the pure-trace map (Theorem 2); Figure 1, Figure 2 provide diagnostics.
  • (D3) Three-route bulk collapse. Reconstruct quadratic kernels from simulations/perturbation theory: the DBI, closed-metric, and CS + routes must share the same bulk piece A g sgn ( Σ ϵ ) I T with A = λ 2 / 8 (Section 5); residuals should enter only through improvements (Figure 3).
  • (D4) Boundary equivalence. Boundary flux ratios R X / Y [ D ] converge to 1 on growing domains (A4) (Section 5.4; Figure 4). Persistent deviations falsify C2.
  • (D5) Exact luminality by identity. After eliminating TT–nonTT mixing with the full-rank 2 × 2 system (Section 6.3), the equal-coefficient identity K G = μ ( a 3 J Δ μ a 3 ) enforces c T = 1 without tuning (Section 7; Figure 7).
  • (D6) DoF count. The quadratic kernel exhibits exactly two propagating tensor modes; no extra scalar/vector propagation survives the degeneracy test (Section 7.3; Figure 8). Any additional mode falsifies our posture at this order.
  • (D7) NLO slope. At next order, the leading P T -even dispersion predicts δ c T 2 ( k ) = b k 2 / Λ 2 (Equation (54)): a log–log slope 2 in clean frequency windows is characteristic (Figure 9; Section 9.1).8
  • (D8) Fermion channel null. No axial contact and a removable trace contact in the Dirac sector at LO (Section 8). Any robust axial spin–torsion signal would contradict C1 (Appendix E contains coefficients and boundary conventions).

Appendix F.2. Expanded Disambiguation Table

Entries are generic at LO on admissible backgrounds with A4 fall-offs; tuned subclasses and specific parameter choices may alter individual cells. See the notes in Appendix F.3 for caveats and representative citations.

Appendix F.3. Notes, Caveats, And Representative Anchors

Parity-odd lines.

Chern–Simons modified gravity (Jackiw–Pi; Alexander–Yunes) generically induces helicity-dependent GW propagation; details depend on the choice of scalar field and coupling. Our scalar- P T projector eliminates parity-odd scalars at LO; see Theorem 1. Representative anchors: [14,15].

Horndeski/DHOST.

Post-GW170817 constraints force c T 1 by parameter tuning or by restricting to subclasses; however, no equal-coefficient identity of the type K G = μ ( ) is generally present, and an extra scalar DoF is typical. See the multimessenger bounds and reviews [27,28,29,30,35,36].

Teleparallel/f(T).

Many models yield c T = 1 at LO; torsion irreps and matter couplings are model-dependent, and our route-equality/flux-ratio diagnostics (C2) do not directly apply. We therefore treat teleparallel cases as outside the present “three-route” posture.

EC/ECSK-like.

With propagating torsion or nonminimal fermion couplings, axial/traceless torsion irreps can be present and Dirac axial contacts typically survive. Our C1 pure-trace map excludes these at the analyzed order; see classic EC/MAG treatments [4,5,6,8].

NLO dispersion.

The δ c T 2 k 2 slope follows from the projected, P T -even closure and the normalization of T (Section 5); in practice, slopes may be obscured by finite-band systematics, so conservative fits are recommended (Section 9).

Boundary posture.

All boundary-sensitive statements assume A4 (PT-invariant compact boundaries or standard AF/FRW fall-offs) and A5 (Nieh–Yan as boundary counterterm). Departures from these postures may alter flux diagnostics and improvement accounting. Operational anchors: covariant phase space and surface charges [20,21,22], and the Holst/Nieh–Yan parity bookkeeping [9,10,11,12].

Representative anchors (one-line map).

Metric-affine/Einstein–Cartan: [4,6,8]. Holst/Nieh–Yan: [9,10,11,12]. Covariant phase space & boundary charges: [20,21,22]. CS-modified gravity: [14,15]. GW170817 speed bounds & implications: [23,24,25,26,27,28,29,30].

References

  1. Palatini, A. Deduzione invariantiva delle equazioni gravitazionali dal principio di Hamilton. Rend. Circ. Mat. Palermo 1919, 43, 203–212. [Google Scholar] [CrossRef]
  2. Kibble, T.W.B. Lorentz invariance and the gravitational field. J. Math. Phys. 1961, 2, 212–221. [Google Scholar] [CrossRef]
  3. Sciama, D.W. On the analogy between charge and spin in general relativity. in Recent Developments in General Relativity (Pergamon, 1962), pp. 415–439.
  4. Hehl, F.W.; von de r Heyde, P.; Kerlick, G.D.; Nester, J.M. General Relativity with Spin and Torsion: Foundations and Prospects. Rev. Mod. Phys. 1976, 48, 393–416. [Google Scholar] [CrossRef]
  5. Hehl, F.W.; Datta, B.K. Nonlinear spinor equation and asymmetric connection in general relativity. J. Math. Phys. 1971, 12, 1334–1339. [Google Scholar] [CrossRef]
  6. Hehl, F.W.; McCrea, J.D.; Mielke, E.W.; Ne’eman, Y. Metric-affine gauge theory of gravity: Field equations, Noether identities, world spinors, and breaking of dilaton invariance. Phys. Rept. 1995, 258, 1–171. [Google Scholar] [CrossRef]
  7. Hammond, R.T. Torsion gravity. Rep. Prog. Phys. 2002, 65, 599–649. [Google Scholar] [CrossRef]
  8. Shapiro, I.L. Physical aspects of the space–time torsion. Phys. Rept. 2002, 357, 113–213. [Google Scholar] [CrossRef]
  9. Nieh, H.T.; Yan, M.L. An identity in Riemann–Cartan geometry. J. Math. Phys. 1982, 23, 373–374. [Google Scholar] [CrossRef]
  10. Chandia, O.; Zanelli, J. Topological invariants, instantons and the chiral anomaly on spaces with torsion. Phys. Rev. D 1997, 55, 7580–7585. [Google Scholar] [CrossRef]
  11. Holst, S. Barbero’s Hamiltonian derived from a generalized Hilbert–Palatini action. Phys. Rev. D 1996, 53, 5966–5969. [Google Scholar] [CrossRef]
  12. Mercuri, S. , Fermions in the Ashtekar–Barbero connection formalism: The Nieh–Yan invariant as a source of the Immirzi parameter. Phys. Rev. D 2006, 73, 084016. [Google Scholar] [CrossRef]
  13. Sotiriou, T.P.; Faraoni, V. “f(R) theories of gravity,” Rev. Mod. Phys. 2010, 82, 451–497. [Google Scholar] [CrossRef]
  14. Jackiw, R.; Pi, S.-Y. Chern–Simons modification of general relativity. Phys. Rev. D 2003, 68, 104012. [Google Scholar] [CrossRef]
  15. Alexander, S.; Yunes, N. Chern–Simons Modified Gravity. Phys. Rept. 2009, 480, 1–55. [Google Scholar] [CrossRef]
  16. Bekenstein, J.D. The Relation between physical and gravitational geometry. Phys. Rev. D 1993, 48, 3641–3647. [Google Scholar] [CrossRef] [PubMed]
  17. Deser, S.; Gibbons, G.W. Born–Infeld–Einstein actions? Class. Quantum Grav. 1998, 15, L35–L39. [Google Scholar] [CrossRef]
  18. Ba nados, M.; Ferreira, P.G. Eddington’s theory of gravity and its progeny. Phys. Rev. Lett. 2010, 105, 011101. [Google Scholar] [CrossRef]
  19. Arnowitt, R.; Deser, S.; Misner, C.W. The Dynamics of General Relativity. in Gravitation: An Introduction to Current Research, ed. L. Witten (Wiley, 1962), pp. 227–265; reprinted in Gen. Rel. Grav. 2008, 40, 1997–2027. [Google Scholar] [CrossRef]
  20. Wald, R.M. General Relativity (University of Chicago Press, 1984).
  21. Iyer, V.; Wald, R.M. Some properties of Noether charge and a proposal for dynamical black hole entropy. Phys. Rev. D 1994, 50, 846–864. [Google Scholar] [CrossRef]
  22. Regge, T.; Teitelboim, C. Role of surface integrals in the Hamiltonian formulation of general relativity. Annals Phys. 1974, 88, 286–318. [Google Scholar] [CrossRef]
  23. B. P. Abbott et al. (LIGO Scientific Collaboration and Virgo Collaboration), GW170817: Observation of Gravitational Waves from a Binary Neutron Star Inspiral. Phys. Rev. Lett. 2017, 119, 161101. [CrossRef]
  24. Abbott, B.P.; et al. (LIGO Scientific Collaboration and Virgo Collaboration), Multi-messenger Observations of a Binary Neutron Star Merger. Astrophys. J. Lett. 2017, 848, L12. [Google Scholar] [CrossRef]
  25. Goldstein, A.; et al. , An Ordinary Short Gamma-Ray Burst with Extraordinary Implications: Fermi-GBM Detection of GRB 170817A. Astrophys. J. Lett. 2017, 848, L14. [Google Scholar] [CrossRef]
  26. Savchenko, V.; et al. INTEGRAL Detection of the First Prompt Gamma-Ray Signal Coincident with the Gravitational-wave Event GW170817. Astrophys. J. Lett. 2017, 848, L15. [Google Scholar] [CrossRef]
  27. Creminelli, P.; Vernizzi, F. Dark Energy after GW170817 and GRB170817A. Phys. Rev. Lett. 2017, 119, 251302. [Google Scholar] [CrossRef]
  28. Baker, T.; Bellini, E.; Ferreira, P.G.; Lagos, M.; Noller, J.; Sawicki, I. Strong constraints on cosmological gravity from GW170817 and GRB170817A. Phys. Rev. Lett. 2017, 119, 251301. [Google Scholar] [CrossRef]
  29. Ezquiaga, J.M.; Zumalacárregui, M. Dark Energy After GW170817: Dead Ends and the Road Ahead. Phys. Rev. Lett. 2017, 119, 251304. [Google Scholar] [CrossRef]
  30. Sakstein, J.; Jain, B. Implications of the Neutron Star Merger GW170817 for Cosmological Scalar–Tensor Theories. Phys. Rev. Lett. 2017, 119, 251303. [Google Scholar] [CrossRef]
  31. Agazie, N.V.; et al. (NANOGrav Collaboration), The NANOGrav 15-year Data Set: Evidence for a Gravitational-Wave Background. Astrophys. J. Lett. 2023, 951, L8. [Google Scholar] [CrossRef]
  32. Antoniadis, J.; et al. (EPTA Collaboration and InPTA Collaboration). The second data release from the European Pulsar Timing Array: Search for signals from new physics. Astron. Astrophys. 2023, 678, A50. [Google Scholar]
  33. Reardon, D.J.; et al. (PPTA Collaboration), Search for an Isotropic Gravitational-Wave Background with the Parkes Pulsar Timing Array. Astrophys. J. Lett. 2023, 951, L7. [Google Scholar] [CrossRef]
  34. Liu, F.A.A.; et al. (CPTA Collaboration), Searching for the nano-Hertz stochastic gravitational-wave background with the Chinese Pulsar Timing Array Data Release I. Res. Astron. Astrophys. 2023, 23, 075024. [Google Scholar]
  35. Kase, R.; Tsujikawa, S. Dark energy in Horndeski theories after GW170817: A review. Int. J. Mod. Phys. D 2019, 28, 1942005. [Google Scholar] [CrossRef]
  36. Langlois, D. Dark Energy and Modified Gravity in Degenerate Higher-Order Scalar–Tensor (DHOST) theories. Int. J. Mod. Phys. D 2019, 28, 1942006. [Google Scholar] [CrossRef]
  37. Obukhov, Y.N. Poincaré gauge gravity: Selected topics. Int. J. Geom. Methods Mod. Phys. 2006, 3, 95–138. [Google Scholar] [CrossRef]
  38. Obukhov, Y.N.; Hehl, F.W. Rotation, acceleration, and gravity in the framework of classical electrodynamics. Phys. Lett. A 2008, 372, 3946–3952. [Google Scholar] [CrossRef]
1
When Σ ϵ = 0 the constructions are understood on patches with Σ ϵ 0 and extended by continuity. The overall sign bookkeeping is carried by sgn ( Σ ϵ ) , so orientation flips are harmless.
2
Independence follows because, before imposing any relation between T μ and μ ϵ , no integration-by-parts identity reduces the linear bilinear T · ϵ to a combination of T 2 and ( ϵ ) 2 up to improvements.
3
We use “rank-one determinant route’’ (formerly “DBI’’-type) to denote the volume deformation det [ 1 + ( 2 3 ) λ τ T ] in our posture; this is not Born–Infeld gravity.
4
Other NLO scalars either renormalize K and G equally (leaving c T unchanged at this order) or reduce to improvements under A4/A5.
5
Frame and spacetime contractions are equivalent once e A μ is inserted; all formulas below hold with either T A B C or T λ μ ν , and we freely move between the two notations.
6
This harmless convention matches Equation (4. 13) of the main text (positivity split) and simplifies a few numerical diagnostics; no physics depends on it.
7
Chern–Simons modified gravity is the canonical benchmark for helicity-split propagation; see e.g., Jackiw–Pi; Alexander–Yunes.
8
Caveat: finite-window fits may be biased by instrumental systematics or astrophysical priors; interpret bounds conservatively.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.
Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.
Prerpints.org logo

Preprints.org is a free preprint server supported by MDPI in Basel, Switzerland.

Subscribe

Disclaimer

Terms of Use

Privacy Policy

Privacy Settings

© 2025 MDPI (Basel, Switzerland) unless otherwise stated