k_f^2 = \frac{k^2+1}{1+k_s \left(\frac{\sqrt{k^2+k_p^2}+\sqrt{k^2+1}}{k_p+1}-1\right)}