trust our results. Unfortunately, since $\<\Phi\>\neq 0$ and it is no longer a simple matter to determine the Goldstone directions, the calculation of the masses must be done for all fields and becomes significantly more involved [^5]. As a result of this, we discovered non-trivial contributions to not only $\Phi_0$’s potential, but $\phi_0$ and $\Phi_1$’s.

Our approach was to start at a point close to the global ISS minimum, compute the one-loop potential to second order in the fields at this point and then minimise the tree-level + one-loop potential [^6]. This new point was then used to re-compute the one-loop potential, allowing us to minimise yet again. This process was repeated until the solution converged and the eigenvalues of the matrix of second derivatives of the full potential were consistent with the expected number of Goldstones (which would not be the case if the effective potential were computed away from the minimum). Regrettably, this is not sufficient to make the computation tractable for large numbers of flavours, as the computation time grows rapidly with the number of flavours and $N_f=4$, $N_c=3$ is already takes prohibitively long. However, for $N_f=3$, $N_c=2$ the one-loop potential can be computed in a reasonable 