SOA/EDFA Carrier-Photon Rate Equations (dim=4)
PARITYS2 · dim 4No clear winner. The survival gap is under 10 percentage points and the balanced-score gap is under 0.05, so neither SolvSRK nor the best baseline clears the win threshold. Either works - choose on cost, licensing, or integration effort. All verdicts →
4-state semiconductor optical amplifier model. States: carrier density N, co-propagating signal photon density, ASE photon density, gain saturation integral h. Stiff from ~ps carrier lifetime vs ~ns photon cavity time scale separation. ND-3 phase reference.
Problem definition
Connelly (2001) IEEE J Quantum Electron; Agrawal & Olsson (1989) IEEE J Quantum Electron
Canonical RHS excerpt from the registered callable used for this benchmark cell. Expand it to verify the state equations; it is not a standalone runnable fixture.
Show canonical RHS excerpt
def _rhs_soa(t, y):
del t
n = max(float(y[0]), 1.0e-9)
p = max(float(y[1]), _PHOTON_FLOOR)
p_ase = max(float(y[2]), _PHOTON_FLOOR)
h = float(y[3])
g = _G_ND * (n - 1.0)
# SSB-validated 2-state carrier-photon core.
dn_tau = _P_ND - n * _DECAY - g * p
dp_tau = g * p - p + _BETA_ND * n
# Passive auxiliaries for dim=4 catalog.
dpa_tau = (p - p_ase) / _TAU_ASE_ND
dh_tau = (n - 1.0 - h) / _TAU_H_ND
if y[1] <= _PHOTON_FLOOR and dp_tau < 0.0:
dp_tau = 0.0
inv_tau_p = 1.0 / _TAU_P
return np.array(
[dn_tau * inv_tau_p, dp_tau * inv_tau_p, dpa_tau * inv_tau_p, dh_tau * inv_tau_p],
dtype=np.float64,
)- Parameters
- _BETA_ND = 5e-07
- _DECAY = 0.0166666666667
- _G_ND = 3.85714285714
- _PHOTON_FLOOR = 1e-12
- _P_ND = 0.0330687830688
- _TAU_ASE_ND = 10
- _TAU_H_ND = 1000
- _TAU_P = 5e-12
- Initial condition
- y(0) = [1.25925925926, 0.0120811287478, 0.0120811287478, 0.259259259259]
- Horizon
- t ∈ [0, 1e-08]
Canonical RHS excerpt captured from the same registered callable used for the published benchmark. Frozen closure values are summarized below; helper imports and solver settings are intentionally omitted.
Fingerprint
Spread: low
Default noise: none
Recommendation snapshot
Clean best: SolvSRK
Noisy best: SciPy RK45
Coverage
14 solver arms · clean + 5 noise levels
Ranked on survival, precision, and speed
Versions & freeze
Methodology →- Freeze
- 2026-08-13
- libsolvsrk
- 2.3.0
- SciPy
- 1.14
- SUNDIALS
- CVODE (bundled backend)
20 seeds/cell default · 14 arms · TRL 4–5 · simulation-lab validated · this page: SOA/EDFA Carrier-Photon Rate Equations (dim=4) (soa-edfa-carrier-photon-rate-equations-dim-4)
Governed SolvTune benchmark freeze; per-arm medians only. RHS definitions and raw trial rows are not published.
Self-reported by Resonix Labs · not independently verified
Results matrix
Pick an objective and a noise level to rank all arms on survival, median SCD, median nfev, and median wall time. Medians across seeds.
Objective
Best overall trade-off of survival, precision, and speed.
Noise level
| # | Solver | Survival | SCD | nfev | Wall | Score |
|---|---|---|---|---|---|---|
| 1 | SolvSRK | 100% | 11.2 | 972 | 21 ms | 0.885 |
| 2 | SciPy RadauSciPy | 100% | 10.3 | 295 | 6 ms | 0.865 |
| 3 | SciPy DOP853SciPy | 100% | 9.3 | 1,142 | 7 ms | 0.840 |
| 4 | SciPy BDFSciPy | 100% | 8.1 | 186 | 7 ms | 0.812 |
| 5 | SciPy LSODASciPy | 100% | 8.1 | 257 | 1 ms | 0.811 |
| 6 | CVODE BDFexternal | 100% | 7.6 | 552 | 10 ms | 0.800 |
| 7 | CVODE Adamsexternal | 100% | 7.6 | 432 | 9 ms | 0.799 |
| 8 | SciPy RK23SciPy | 100% | 7.6 | 845 | 7 ms | 0.799 |
| 9 | SciPy RK45SciPy | 100% | 7.3 | 992 | 7 ms | 0.793 |
| 10 | Tsit5external | 100% | 6.7 | 1,038 | 862 ms | 0.779 |
At Clean, best balanced arm is SolvSRK.
Values are medians across seeds, measured by Resonix Labs on Resonix hardware and not independently verified; nfev and wall are on reference lab hardware (indicative). Under injected noise only SolvSRK and the SciPy arms are run. How we measure accuracy → · Verification status →
Cite this page
Replace the access date. Pin the freeze ID and library versions when comparing against a later export. Cite it as what it is - a self-reported vendor benchmark, not an independently verified result. The note field says so; please keep it.
@misc{resonix_evidence_soa_edfa_carrier_photon_rate_equations_dim_4_2026,
title = {Resonix Evidence Portal: SOA/EDFA Carrier-Photon Rate Equations (dim=4)},
author = {{Resonix Labs (Canada) Inc.}},
year = {2026},
howpublished = {\url{https://resonixusa.com/evidence/problems/soa-edfa-carrier-photon-rate-equations-dim-4}},
note = {Self-reported vendor benchmark; internally generated by Resonix Labs and not independently verified. Accessed YYYY-MM-DD. Freeze 2026-08-13; libsolvsrk 2.3.0; SciPy 1.14.}
}Related
TRL 4–5 · simulation-lab validated · 398 problems · 14 solver arms · clean + 5 noise levels
Freeze: 2026-08-13 · scipy 1.14 · libsolvsrk 2.3.0 · Methodology
Self-reported by Resonix Labs · not independently verified · Verification status