\(\rho_{\text{hnf}}=\left(1-\phi_{2}\right)\left[\left(1-\phi_{1}\right)\rho_{\text{bf}}+\phi_{1}\rho_{np,1}\right]+\phi_{2}\rho_{np,2}\) (1)