function lambdaMinusScaled = lambdaMinusScaledFun(rhoi,rhoj,sigmaEffScaling,sigma)

lambdaMinusScaled = ((rhoi + rhoj)-((rhoi - rhoj)^2 + 4 * (sigmaEffScaling(1) * sigmaEffScaling(3) * sigma)^2)^0.5)/2;