docs(research): rendering equation as observerless observer

The rendering equation (Kajiya 1986) is the continuous limit of the
16D chiral observerless observer framework.

Mapping:
- BRDF f_r(ω_i, ω_o) = chiral coupling (braid crossing σ_i^ε)
- Irradiance cosine (ω_i · n) = q-profile (L₁/L₀ = poloidal/toroidal)
- Hemisphere integral ∫_Ω = CRT sum over n/2 channels
- Neumann series L_o = Σ Kᵏ[L_e] = eigensolid convergence
- Fixed-point recursion (L_o on both sides) = observerless observer

The Sidon property = discrete Nyquist criterion: channels must be
sufficiently separated to avoid aliasing in the directional integral.

Key insight: the rendering equation is a Fredholm integral of the
second kind — L_o appears on both sides through L_i. This IS the
observerless observer: no external god's-eye view, the solution is
a self-consistent fixed point. The eigensolid convergence
(BraidEigensolid.lean) is the discrete Neumann series.

The q-profile determines the BRDF shape:
- q >> 1: diffuse (many orthogonal channels, low coupling)
- q < 1: specular (few dominant channels, high coupling)
- q = 1: degenerate (single channel, no diversity)

This explains the q-profile sweep result: q > 1 = 100% Sidon because
low coupling = channels don't interfere (BRDF-orthogonal).

docs/research/RENDERING_EQUATION_OBSERVERLESS.md

@@ -0,0 +1,212 @@
+# The Rendering Equation as Observerless Observer
+
+**Status:** THEORETICAL — connects rendering equation to 16D chiral framework
+**Date:** 2026-07-04
+**Depends on:** `CHIRAL_CRT_MULTIPLEXING.md`, `HCMR_CRT_MULTIPLEXER.md`,
+`INVARIANT_COMPUTATION_GEOMETRY.md`, `OCTAGON_PRINCIPLE.md`
+**Source equation:** Kajiya (1986), "The Rendering Equation"
+
+---
+
+## 1. The Rendering Equation
+
+$$L_o(\mathbf{x}, \omega_o) = L_e(\mathbf{x}, \omega_o) + \int_{\Omega} f_r(\mathbf{x}, \omega_i, \omega_o) L_i(\mathbf{x}, \omega_i) (\omega_i \cdot \mathbf{n}) d\omega_i$$
+
+where:
+- `L_o(x, ω_o)` = outgoing radiance at point x in direction ω_o
+- `L_e(x, ω_o)` = emitted radiance (self-illumination)
+- `f_r(x, ω_i, ω_o)` = BRDF (bidirectional reflectance distribution function)
+- `L_i(x, ω_i)` = incoming radiance from direction ω_i
+- `(ω_i · n)` = irradiance factor (cosine with surface normal)
+- `Ω` = unit hemisphere above the surface
+
+## 2. Why This Is the Observerless Observer
+
+The rendering equation is a **Fredholm integral equation of the second kind**:
+`L_o` appears on both sides. The incoming radiance `L_i(x, ω_i)` is itself
+the outgoing radiance `L_o(x', ω_o)` at another point `x'` visible along
+direction `ω_i`. The equation is recursive:
+
+    L_o = L_e + K[L_o]
+
+where `K[·]` is the integral operator (the light transport kernel).
+
+This IS the observerless observer:
+- No external "god's-eye" reference frame
+- The observer (viewer at ω_o) and the observed (scene via L_i) are coupled
+- The solution is a **fixed point**: L_o = (I - K)⁻¹ L_e (Neumann series)
+- The observation emerges from self-consistency, not from an external frame
+
+In the 16D chiral framework, this is exactly the structure:
+- The 8-strand braid is a fixed point under crossing operations
+- The eigensolid convergence (BraidEigensolid.lean) IS the Neumann series
+  convergence: repeated application of the light transport operator
+- The "observerless observer" = no preferred direction = all directions
+  are treated equally in the hemisphere integral
+
+## 3. The Mapping: Rendering Equation ↔ 16D Chiral
+
+### 3.1 Component Map
+
+| Rendering equation | 16D chiral framework | Meaning |
+|---|---|---|
+| `L_o(x, ω_o)` | Strand output | What the observer strand "sees" |
+| `L_e(x, ω_o)` | Identity component (a mod L₀) | Intrinsic emission (poloidal) |
+| `f_r(x, ω_i, ω_o)` | Braid crossing σ_i | Chiral coupling (how i→o) |
+| `L_i(x, ω_i)` | Reflection component (S-a mod Lᵢ) | Incoming from environment (toroidal) |
+| `(ω_i · n)` | q-profile (L₁/L₀ ratio) | Angle-dependent irradiance factor |
+| `∫_Ω dω_i` | CRT sum over all channels | Hemisphere = all chiral channels |
+| Fixed-point (L_o = L_e + K[L_o]) | Observerless observer | No external reference frame |
+
+### 3.2 The BRDF as Chiral Coupling
+
+The BRDF `f_r(x, ω_i, ω_o)` encodes how light from direction ω_i reflects
+into direction ω_o. This is DIRECTIONAL — it depends on both angles.
+
+In the chiral framework:
+- Each braid crossing σ_i has chirality εᵢ ∈ {+1, -1}
+- σ_i⁺¹ = over-crossing = light reflects "over" (positive BRDF lobe)
+- σ_i⁻¹ = under-crossing = light reflects "under" (negative BRDF lobe)
+- The BRDF IS the chiral coupling: f_r(ω_i, ω_o) = f(σ_i^ε)
+
+A specular surface (mirror) has a sharp BRDF lobe = single chiral crossing.
+A diffuse surface (Lambertian) has uniform BRDF = all chiral configurations
+equally likely. The q-profile determines the BRDF shape:
+- q >> 1 (translation-dominated): diffuse-like (all channels active)
+- q < 1 (rotation-dominated): specular-like (few channels dominate)
+- q = 1: degenerate (single channel, no diversity)
+
+### 3.3 The Irradiance Factor as q-Profile
+
+The `(ω_i · n)` term is the cosine of the angle between incoming light and
+the surface normal. This is the "efficiency" of energy transfer.
+
+In the chiral framework:
+- `n` = the identity axis L₀ (the "normal" = the intrinsic direction)
+- `ω_i` = the reflection axis L₁ (the "incoming" = the toroidal direction)
+- `(ω_i · n)` = cos(angle between L₀ and L₁) ≈ L₁/L₀ = q
+
+When q < 1 (L₁ < L₀): the reflection axis is "aligned" with the identity
+(normal-like) → high irradiance → high coupling
+When q > 1 (L₁ > L₀): the reflection axis is "perpendicular" → low irradiance
+→ low coupling but more channels
+
+This explains the q-profile sweep result: q > 1 has 100% Sidon rate because
+low irradiance = low coupling = channels don't interfere (orthogonal).
+q < 1 has lower Sidon rate because high irradiance = high coupling = channels
+interfere (collisions).
+
+### 3.4 The Hemisphere Integral as CRT Sum
+
+The integral `∫_Ω dω_i` sums over all incoming directions in the hemisphere.
+This is the continuous version of summing over all chiral channels.
+
+In the discrete (CRT) framework:
+- The hemisphere Ω is discretized into n/2 chiral channels
+- Each channel = one (identity, reflection) pair
+- The integral becomes: Σ_{j=1}^{n/2} f_r(j) L_i(j) q_j
+- The Sidon property ensures channels are orthogonal (non-interfering)
+- Without Sidon: channels collide → the integral has aliasing artifacts
+
+## 4. The Neumann Series = Eigensolid Convergence
+
+### 4.1 Continuous Case (Rendering Equation)
+
+The rendering equation's solution is the Neumann series:
+
+    L_o = L_e + K[L_e] + K²[L_e] + K³[L_e] + ...
+    L_o = (I - K)⁻¹ L_e = Σ_{k=0}^∞ Kᵏ[L_e]
+
+This converges when the operator norm `||K|| < 1` (physically: energy is
+lost at each bounce, no perfect mirrors in a closed room).
+
+### 4.2 Discrete Case (BraidEigensolid)
+
+The eigensolid convergence (BraidEigensolid.lean) is the SAME series:
+
+    BraidState_final = Σ_{k=0}^∞ crossStepᵏ(BraidState_initial)
+
+where `crossStep` is the braid crossing operator (the discrete analog of
+the light transport kernel K).
+
+Convergence condition: the spectral radius of crossStep < 1.
+In HCMR terms: self_loop_prob < 1 (not fully contended).
+In rendering terms: ||K|| < 1 (energy lost per bounce).
+
+### 4.3 The Connection
+
+The eigensolid IS the rendering equation's solution in the discrete chiral
+framework:
+- Each braid crossing = one light bounce
+- The Sidon labels = the radiance values at each point
+- The crossStep operator = the light transport kernel K
+- The fixed point (eigensolid) = the steady-state radiance distribution
+- The "observerless observer" = the recursive fixed-point structure
+
+## 5. Implications for the Multiplexer
+
+### 5.1 The BRDF Determines Channel Quality
+
+In the CRT multiplexer, each channel's quality depends on the BRDF:
+- High BRDF lobe (specular) = strong coupling = one dominant channel
+- Low BRDF lobe (diffuse) = weak coupling = many channels, low each
+- The q-profile controls the BRDF shape
+
+### 5.2 The Rendering Equation Is the Continuous Limit
+
+The CRT multiplexer is the DISCRETE version of the rendering equation:
+- n/2 channels = n/2 directional samples of the hemisphere
+- CRT sum = discrete hemisphere integral
+- Sidon orthogonality = channels don't alias (Nyquist criterion)
+- The Neumann series = eigensolid convergence
+
+As n → ∞, the CRT multiplexer approaches the rendering equation.
+The Sidon property is the discrete Nyquist criterion: channels must be
+sufficiently separated to avoid aliasing.
+
+### 5.3 The Observerless Observer Is the Fixed Point
+
+The "observerless observer" from INVARIANT_COMPUTATION_GEOMETRY.md is
+the rendering equation's fixed point:
+- No external observer (L_o is defined self-consistently)
+- The observation emerges from the integral structure
+- The frame-independent invariants are the BRDF's symmetries
+
+In the chiral framework:
+- The braid's fixed point (eigensolid) = the steady-state radiance
+- The Sidon property = the BRDF's directional orthogonality
+- The q-profile = the BRDF's angular distribution
+
+## 6. Practical Implication: BRDF-Guided Channel Selection
+
+If the rendering equation is the continuous limit, then:
+1. The BRDF of a physical surface determines the optimal q-profile
+2. Specular surfaces → q < 1 (few dominant channels, high coupling)
+3. Diffuse surfaces → q > 1 (many channels, low coupling, orthogonal)
+4. The Sidon filter selects channels that are "BRDF-orthogonal"
+
+This means: for a given physical system (surface, network, workload),
+the BRDF (directional response function) determines which chiral
+configurations are useful. The Sidon filter selects exactly those.
+
+## 7. claim_boundary
+
+```
+rendering-equation-observerless:theoretical-connection:continuous-limit
+```
+
+The rendering equation (Kajiya 1986) is the continuous limit of the
+16D chiral observerless observer framework. The mapping:
+
+- BRDF = chiral coupling (braid crossing with chirality)
+- Irradiance cosine = q-profile (poloidal/toroidal ratio)
+- Hemisphere integral = CRT sum over channels
+- Neumann series = eigensolid convergence
+- Fixed-point recursion = observerless observer
+
+The Sidon property is the discrete Nyquist criterion: channels must be
+sufficiently separated to avoid aliasing in the directional integral.
+As n → ∞, the CRT multiplexer approaches the rendering equation.
+
+**OPEN:** Can the BRDF of a physical surface be used to predict the
+optimal q-profile for the CRT multiplexer?