Finite-Response Coupled Field Dynamics

Finite-Response Coupled Field Dynamics (FRCFD)

FRCFD Core Equations

∂²S/∂t² − c²∇²S + βS³ = σ(x,t)|Ψ|²
∂²Ψ/∂t² − v²∇²Ψ + μΨ + λ|Ψ|²Ψ = κSΨ

🔍 Equation 1: Medium Field S(x,t)

∂²S/∂t² − c²∇²S + βS³ = σ(x,t)|Ψ|²

Term-by-term structure:

1. Inertial term
∂²S/∂t² — second time derivative defining temporal acceleration.

2. Propagation term
−c²∇²S — Laplacian operator enforcing finite-speed spatial propagation.

3. Nonlinear saturation
+βS³ — cubic self-interaction dominant for large |S|, ensuring bounded solutions.

4. Source term
σ(x,t)|Ψ|² — local driving proportional to matter-field intensity.

Compact structure:

□₍c₎ S + βS³ = σ|Ψ|²

where:

□₍c₎ = ∂²/∂t² − c²∇²

🔍 Equation 2: Matter Field Ψ(x,t)

∂²Ψ/∂t² − v²∇²Ψ + μΨ + λ|Ψ|²Ψ = κSΨ

Term-by-term structure:

1. Inertial term
∂²Ψ/∂t²

2. Propagation term
−v²∇²Ψ

3. Linear term
+μΨ — sets characteristic frequency scale.

4. Nonlinear self-interaction
+λ|Ψ|²Ψ — cubic nonlinearity producing amplitude-dependent response.

5. Coupling term
κSΨ — bilinear interaction linking Ψ to the substrate field S.

Compact structure:

□₍v₎ Ψ + μΨ + λ|Ψ|²Ψ = κSΨ

🔁 Coupled System Structure

Forward coupling: Ψ → S : σ|Ψ|²
Back coupling: S → Ψ : κSΨ

⚙️ Mathematical Properties

1. Hyperbolic operators
Both equations contain ∂²/∂t² − c²∇² → finite propagation speed and causal structure.

2. Nonlinearity
S³ and |Ψ|²Ψ → amplitude-dependent dynamics and mode coupling.

3. Boundedness
βS³ dominates at large |S| → prevents divergence.

4. Energy exchange
σ|Ψ|² and κSΨ → bidirectional transfer between fields.

📊 Limiting Cases

Case 1: No coupling (σ = κ = 0)

□₍c₎ S + βS³ = 0
□₍v₎ Ψ + μΨ + λ|Ψ|²Ψ = 0

Case 2: Linear regime (S, Ψ → 0)

□₍c₎ S = 0
□₍v₎ Ψ + μΨ = 0

Case 3: Strong amplitude
βS³ and λ|Ψ|²Ψ dominate → saturation and bounded states.

📌 System Summary (Equation Form Only)

∂²S/∂t² − c²∇²S + βS³ = σ|Ψ|²
∂²Ψ/∂t² − v²∇²Ψ + μΨ + λ|Ψ|²Ψ = κSΨ

🔥 Key Takeaways (Strictly Mathematical)

• second-order hyperbolic system
• nonlinear cubic self-interactions
• bilinear coupling
• finite propagation speeds c and v
• bounded solutions for β, λ > 0
• supports mode coupling and energy transfer

1. Lagrangian Density

We seek a Lagrangian density L(S, Ψ) such that the Euler–Lagrange variations reproduce the field equations for S and Ψ:

δL/δS ⇒ Eq. (S),    δL/δΨ ⇒ Eq. (Ψ)

Consistent choice:

L = 1/2 (∂ₜS)² − c²/2 |∇S|² − β/4 S⁴
  + 1/2 (∂ₜΨ)² − v²/2 |∇Ψ|² − μ/2 Ψ² − λ/4 Ψ⁴ − κ/2 SΨ²

We have replaced |Ψ|² → Ψ² (real scalar field assumption). The source term σ(x,t) is treated as external and is not included in the conservative Lagrangian.

2. Euler–Lagrange Derivation

General form:

∂L/∂ϕ − ∂ₜ(∂L/∂(∂ₜϕ)) − ∇·(∂L/∂(∇ϕ)) = 0

For S:

∂L/∂S        = −βS³ − κ/2 Ψ²
∂L/∂(∂ₜS)    = ∂ₜS
∂L/∂(∇S)     = −c²∇S

Result:

∂²S/∂t² − c²∇²S + βS³ = −κ/2 Ψ²

For Ψ:

∂L/∂Ψ        = −μΨ − λΨ³ − κSΨ

Result:

∂²Ψ/∂t² − v²∇²Ψ + μΨ + λΨ³ = κSΨ

To match the original coupled system with a source term σ|Ψ|², we identify σ = κ/2, so the full dynamics are derivable from a single Lagrangian.

3. Energy (Hamiltonian Density)

Definition:

H = Σϕ (∂ₜϕ · πϕ) − L

with conjugate momenta:

π_S  = ∂L/∂(∂ₜS)  = ∂ₜS
π_Ψ  = ∂L/∂(∂ₜΨ) = ∂ₜΨ

Hamiltonian density:

H = 1/2 (∂ₜS)² + c²/2 |∇S|² + β/4 S⁴
  + 1/2 (∂ₜΨ)² + v²/2 |∇Ψ|² + μ/2 Ψ² + λ/4 Ψ⁴ + κ/2 SΨ²

4. Conserved Quantities

(A) Energy conservation

d/dt ∫ H d³x = 0

(B) Momentum conservation (spatially uniform case):

P = ∫ d³x (∂ₜS ∇S + ∂ₜΨ ∇Ψ)

(C) Noether symmetries: time translations → energy conservation; space translations → momentum conservation. There is no U(1) symmetry (real field + coupling), so no conserved charge.

5. Interpretation of Energy Terms

Kinetic: 1/2 (∂ₜS)², 1/2 (∂ₜΨ)²

Spatial (gradient): c²/2 |∇S|², v²/2 |∇Ψ|²

Potential: β/4 S⁴, μ/2 Ψ², λ/4 Ψ⁴

Coupling: κ/2 SΨ² → energy exchange channel between substrate and matter-field.

Final Compact Form

Lagrangian:

L = 1/2 (∂ₜS)² − c²/2 |∇S|² − β/4 S⁴
  + 1/2 (∂ₜΨ)² − v²/2 |∇Ψ|² − μ/2 Ψ² − λ/4 Ψ⁴ − κ/2 SΨ²

Hamiltonian:

H = 1/2 (∂ₜS)² + c²/2 |∇S|² + β/4 S⁴
  + 1/2 (∂ₜΨ)² + v²/2 |∇Ψ|² + μ/2 Ψ² + λ/4 Ψ⁴ + κ/2 SΨ²

1. Relativistic Formulation

Spacetime coordinates:

x^μ = (t, x),   μ = 0, 1, 2, 3

Metric (flat spacetime):

η_{μν} = diag(1, −1, −1, −1)

d'Alembert operator:

□ = ∂_μ ∂^μ = ∂_t^2 − ∇^2

Relativistic units: set c = v = 1 (can be restored later).

Relativistic Field Equations

□S + βS^3 = σ|Ψ|^2
□Ψ + μΨ + λΨ^3 = κSΨ

2. Covariant Lagrangian Density

L = 1/2 ∂_μS ∂^μS + 1/2 ∂_μΨ ∂^μΨ − V(S, Ψ)

Potential term:

V(S, Ψ) = β/4 S^4 + μ/2 Ψ^2 + λ/4 Ψ^4 + κ/2 SΨ^2

Euler–Lagrange equations:

∂_μ∂^μ S + ∂V/∂S = 0
∂_μ∂^μ Ψ + ∂V/∂Ψ = 0

3. Stress–Energy Tensor

Canonical (Noether) form:

T_{μν} = Σ_{ϕ = S, Ψ} (∂L/∂(∂^μϕ)) ∂_νϕ − η_{μν} L

Derivatives:

∂L/∂(∂^μS) = ∂_μS
∂L/∂(∂^μΨ) = ∂_μΨ

Final stress–energy tensor:

T_{μν} = ∂_μS ∂_νS + ∂_μΨ ∂_νΨ − η_{μν} L

4. Expanded Form

T_{μν} = ∂_μS ∂_νS + ∂_μΨ ∂_νΨ
        − η_{μν} [ 1/2 (∂_αS ∂^αS) + 1/2 (∂_αΨ ∂^αΨ) − V(S, Ψ) ]

5. Energy Density

Energy density: T₀₀

T_{00} = 1/2 (∂_tS)^2 + 1/2 |∇S|^2
       + 1/2 (∂_tΨ)^2 + 1/2 |∇Ψ|^2
       + V(S, Ψ)

6. Momentum Density

T_{0i} = ∂_tS ∂_iS + ∂_tΨ ∂_iΨ

7. Conservation Law

∂_μ T^{μν} = 0

Energy conservation:

d/dt ∫ T_{00} d^3x = 0

Momentum conservation:

d/dt ∫ T_{0i} d^3x = 0

8. Physical Structure of T_{μν}

• Energy density: T₀₀
• Energy flux / momentum: T_0i
• Stress tensor: T_ij

Field gradients (∂_μS ∂_νS, ∂_μΨ ∂_νΨ), potential energy V(S, Ψ), and the coupling term κ/2 SΨ² all contribute to the local energy–momentum distribution.

Final Covariant System

Lagrangian:

L = 1/2 ∂_μS ∂^μS + 1/2 ∂_μΨ ∂^μΨ
  − (β/4 S^4 + μ/2 Ψ^2 + λ/4 Ψ^4 + κ/2 SΨ^2)

Stress–energy tensor:

T_{μν} = ∂_μS ∂_νS + ∂_μΨ ∂_νΨ − η_{μν} L

1. Promote Flat → Curved Spacetime

Replace flat metric and partial derivatives with curved-space objects:

η_{μν} → g_{μν}(x)
∂_μ   → ∇_μ

where:

• g_{μν}(x) = spacetime metric
• ∇_μ = covariant derivative

Curved d'Alembertian:

□_g ϕ = ∇_μ ∇^μ ϕ
      = (1/√−g) ∂_μ (√−g g^{μν} ∂_ν ϕ)

2. Action (GR + Coupled Fields)

Total action:

S_total = ∫ d^4x √−g [ (1/16πG) R + L_fields ]

Components:

Einstein–Hilbert term:

(1/16πG) R

Field Lagrangian:

L_fields = 1/2 g^{μν} ∇_μS ∇_νS
         + 1/2 g^{μν} ∇_μΨ ∇_νΨ
         − V(S, Ψ)

Potential:

V(S, Ψ) = β/4 S^4 + μ/2 Ψ^2 + λ/4 Ψ^4 + κ/2 SΨ^2

3. Field Equations (Curved Spacetime)

Scalar S:

□_g S + βS^3 = −(κ/2) Ψ^2

Field Ψ:

□_g Ψ + μΨ + λΨ^3 = κSΨ

4. Einstein Field Equations

Varying the action with respect to g_{μν} gives:

G_{μν} = 8πG T_{μν}

Einstein tensor:

G_{μν} = R_{μν} − 1/2 g_{μν} R

5. Stress–Energy Tensor (Curved Form)

T_{μν} = ∇_μS ∇_νS + ∇_μΨ ∇_νΨ − g_{μν} L_fields

Expanded:

T_{μν} = ∇_μS ∇_νS + ∇_μΨ ∇_νΨ
       − g_{μν} [ 1/2 g^{αβ} ∇_αS ∇_βS
                 + 1/2 g^{αβ} ∇_αΨ ∇_βΨ
                 − V(S, Ψ) ]

6. Conservation Law

∇_μ T^{μν} = 0

This is required for consistency with General Relativity.

7. Full Coupled System

(A) Geometry:

G_{μν} = 8πG T_{μν}

(B) Fields:

□_g S + βS^3 = −(κ/2) Ψ^2
□_g Ψ + μΨ + λΨ^3 = κSΨ

8. Interpretation (Strictly Mathematical)

• Metric g_{μν} is determined by T_{μν}.
• T_{μν} is determined by S and Ψ.
• Fields evolve in the geometry they generate.
• Result: a fully coupled nonlinear PDE system.

9. Weak-Field Limit (Recover Newtonian Gravity)

Assume small perturbations:

g_{μν} = η_{μν} + h_{μν},   |h_{μν}| ≪ 1

Then, in the Newtonian limit:

∇^2 Φ ∼ T_{00}

with

T_{00} ≈ 1/2 (∂_tS)^2 + 1/2 (∂_tΨ)^2 + V(S, Ψ)

1. Promote Flat → Curved Spacetime

Replace flat metric and partial derivatives with curved-space objects:

η_{μν} → g_{μν}(x)
∂_μ   → ∇_μ

where:

• g_{μν}(x) = spacetime metric
• ∇_μ = covariant derivative

Curved d'Alembertian:

□_g ϕ = ∇_μ ∇^μ ϕ
      = (1/√−g) ∂_μ (√−g g^{μν} ∂_ν ϕ)

2. Action (GR + Coupled Fields)

Total action:

S_total = ∫ d^4x √−g [ (1/16πG) R + L_fields ]

Einstein–Hilbert term:

(1/16πG) R

Field Lagrangian:

L_fields = 1/2 g^{μν} ∇_μS ∇_νS
         + 1/2 g^{μν} ∇_μΨ ∇_νΨ
         − V(S, Ψ)

Potential:

V(S, Ψ) = β/4 S^4 + μ/2 Ψ^2 + λ/4 Ψ^4 + κ/2 SΨ^2

3. Field Equations (Curved Spacetime)

Substrate field S:

□_g S + βS^3 = −(κ/2) Ψ^2

Matter field Ψ:

□_g Ψ + μΨ + λΨ^3 = κSΨ

4. Einstein Field Equations

Variation with respect to g_{μν} gives:

G_{μν} = 8πG T_{μν}

Einstein tensor:

G_{μν} = R_{μν} − 1/2 g_{μν} R

5. Stress–Energy Tensor (Curved Form)

T_{μν} = ∇_μS ∇_νS + ∇_μΨ ∇_νΨ − g_{μν} L_fields

Expanded:

T_{μν} = ∇_μS ∇_νS + ∇_μΨ ∇_νΨ
       − g_{μν} [
           1/2 g^{αβ} ∇_αS ∇_βS
         + 1/2 g^{αβ} ∇_αΨ ∇_βΨ
         − V(S, Ψ)
       ]

6. Conservation Law

∇_μ T^{μν} = 0

This ensures consistency with General Relativity.

7. Full Coupled System

(A) Geometry:

G_{μν} = 8πG T_{μν}

(B) Fields:

□_g S + βS^3 = −(κ/2) Ψ^2
□_g Ψ + μΨ + λΨ^3 = κSΨ

8. Interpretation (Strictly Mathematical)

• Metric g_{μν} is determined by T_{μν}.
• T_{μν} is determined by S and Ψ.
• Fields evolve in the geometry they generate.
• Result: a fully coupled nonlinear PDE system.

9. Weak-Field Limit (Recover Newtonian Gravity)

Assume small perturbations:

g_{μν} = η_{μν} + h_{μν},   |h_{μν}| ≪ 1

Then:

∇^2 Φ ∼ T_{00}

with:

T_{00} ≈ 1/2 (∂_tS)^2 + 1/2 (∂_tΨ)^2 + V(S, Ψ)

---

PART I — Linear Perturbations

1. Background + Perturbations

Start from homogeneous background fields:

S(t),  Ψ(t),  a(t)

Introduce perturbations:

S(t,x)  = S(t)  + δS(t,x)
Ψ(t,x)  = Ψ(t)  + δΨ(t,x)
g_{μν}  = g^{FRW}_{μν} + δg_{μν}

2. Fourier Decomposition

δS(t,x)  = ∫ d^3k δS_k(t)  e^{i k·x}
δΨ(t,x)  = ∫ d^3k δΨ_k(t)  e^{i k·x}

3. Linearized Equations

Substrate perturbation:

δS̈_k + 3H δṠ_k + (k²/a² + 3βS²) δS_k = −κΨ δΨ_k

Matter-field perturbation:

δΨ̈_k + 3H δΨ̇_k + (k²/a² + μ + 3λΨ² − κS) δΨ_k = κΨ δS_k

4. Structure

These form a coupled oscillator system:

δS_k ↔ δΨ_k

• time-dependent effective masses
• Hubble damping
• cross-coupling

5. Power Spectrum

P(k) = |δΨ_k|²
P(k) = (k³ / 2π²) |δΨ_k|²

Encodes:

• distribution of fluctuations
• structure formation seeds
• CMB signatures

6. Expected Behavior

Large scales (k ≪ aH):

• modes freeze
• spectrum becomes constant

Small scales (k ≫ aH):

• oscillatory behavior
• wave-like propagation

---

PART II — Numerical Simulation (Python)

Minimal working simulation for background + one perturbation mode:

import numpy as np
import matplotlib.pyplot as plt
from scipy.integrate import solve_ivp

# Parameters
beta = 0.1
lambda_ = 0.1
mu = 1.0
kappa = 0.5
G = 1.0

k_mode = 0.1  # perturbation scale

# Initial conditions
S0, dS0 = 0.5, 0.0
Psi0, dPsi0 = 0.5, 0.0
a0 = 1.0

dS_k0 = 1e-5
dPsi_k0 = 1e-5

# Time span
t_span = (0, 50)
t_eval = np.linspace(*t_span, 1000)

def system(t, y):
    S, dS, Psi, dPsi, a, dS_k, ddS_k, dPsi_k, ddPsi_k = y

    # Energy density
    V = (beta/4)*S**4 + (mu/2)*Psi**2 + (lambda_/4)*Psi**4 + (kappa/2)*S*Psi**2
    rho = 0.5*dS**2 + 0.5*dPsi**2 + V

    H = np.sqrt((8*np.pi*G/3)*rho)

    # Background equations
    ddS = -3*H*dS - beta*S**3 - (kappa/2)*Psi**2
    ddPsi = -3*H*dPsi - mu*Psi - lambda_*Psi**3 + kappa*S*Psi

    # Perturbations
    k_phys = k_mode / a

    ddS_k_new = (
        -3*H*ddS_k
        - (k_phys**2 + 3*beta*S**2)*dS_k
        - kappa*Psi*dPsi_k
    )

    ddPsi_k_new = (
        -3*H*ddPsi_k
        - (k_phys**2 + mu + 3*lambda_*Psi**2 - kappa*S)*dPsi_k
        + kappa*Psi*dS_k
    )

    # Scale factor evolution
    da = H * a

    return [
        dS, ddS,
        dPsi, ddPsi,
        da,
        ddS_k, ddS_k_new,
        ddPsi_k, ddPsi_k_new
    ]

# Initial vector
y0 = [
    S0, dS0,
    Psi0, dPsi0,
    a0,
    dS_k0, 0.0,
    dPsi_k0, 0.0
]

# Solve
sol = solve_ivp(system, t_span, y0, t_eval=t_eval)

# Extract
t = sol.t
Psi_k = sol.y[7]

# Power spectrum
P_k = np.abs(Psi_k)**2

# Plot
plt.figure()
plt.plot(t, P_k)
plt.xlabel("Time")
plt.ylabel("Power Spectrum P(k)")
plt.title("Evolution of Perturbation Power")
plt.show()

Finite-Response Coupled Field Dynamics (FRCFD)

1. Mathematical Core

Substrate Equation (S):
Models the Reactive Medium using a nonlinear cubic interaction term βS³, which enforces bounded field amplitudes and removes mathematical singularities.

Matter-Field Equation (Ψ):
Describes localized excitations through a bilinear coupling κSΨ, creating a bidirectional feedback loop where matter modifies the substrate and the substrate modifies matter.

2. Spatio-Temporal Dynamics

Causal Constraints:
The V‑shaped light-cone structure in Ψ evolution demonstrates a finite-response substrate with a bounded propagation speed c, ruling out instantaneous information transfer.

Interference:
Ripple patterns and wavefront interactions reveal interference, dispersion, and nonlinear scattering as excitations propagate through the medium.

3. Physical Coupling and Stability

Correlated Oscillations:
Mid-time field profiles show strong correlation between S and Ψ, where matter-induced stress “dents” the substrate and guides propagation.

Emergent Coherence:
Regions of high-amplitude coupling remain coherent despite surrounding high-frequency jitter, demonstrating the mechanical origin of stable physical structures.

4. Irreversibility and Unification

Spectral Entropy:
Nonlinear mode-mixing transfers energy from coherent modes into a high-frequency noise floor, producing a mechanical Arrow of Time through monotonic entropy growth.

Systemic Synthesis:
Both classical (large-scale) and stochastic (quantum-scale) behaviors emerge from the same coupled system. This resolves apparent paradoxes of non-locality and singular density through finite-response realism and scale-dependent approximations.