deletions | additions       diff --git a/bacoso2.tex b/bacoso2.tex  index 4a47e9e..83614eb 100644  --- a/bacoso2.tex  +++ b/bacoso2.tex 
...
In order to not miss crucial seed structures which ultimately led to the experimental structure, we found that the randomly generated initial population of structures must be sufficiently large. An initial population of size 300 was sufficient with a single generation size of 60. In total $\sim 700$ metastable structures were produced in 8 generations. In addition, spin polarization is crucial for the local relaxations performed within each USPEX generation in order to find the experimental structure. When non-spin-polarized DFT was used, we could not find some of the lowest energy materials (including the observed structure)  What causes the relative shift in energies as function of U? We first roughly grouped the structures by the behavior of the $E$ vs. $U-J$ curve, which we have indicated by colors in Fig.~\ref{fig:reordering}. For example, in red are a small group of structures which are found to be the lowest energy structures in paramagnetic LDA calculations, but which quickly shoot upwards in energies as $U-J$ is increased.  In order to rationalize this behavior, we  Roughly, the correction depends on the occupation of the 3d orbitals to which U is applied, namely the energy difference between two states with occupations n1 and n2 is roughly $∆E \sim U(n1-1/2)(n2-n1)$. [Gabi, Ran, this is just my hypothesis. Ran, could you check this from the data?]