pub fn dissociation_equilibrium_fraction_kernel<R>(
k_eq: EquilibriumConstant<R>,
nuclei_density: R,
) -> Result<DissociationFraction<R>, PhysicsError>where
R: RealField + FromPrimitive,Expand description
Equilibrium dissociation fraction of a diatomic pool A₂ ⇌ 2A at fixed
nuclei density: solve [A]²/[A₂] = K with [A₂] = (n_nuclei − [A])/2 for
the atom share x = [A]/n_nuclei. Closed form:
[A] = (−K + √(K² + 8·K·n_nuclei))/4. k_eq and nuclei_density must be
in one consistent concentration basis.