30
Figure 3.7: Corresponding density profiles of the atomic cloud for N = 1083 at Λ = 50
tigation using grand-canonical mode is going on and we hope to find out at what particle
number the densities in radial direction begin to separate.
3.4 Future Research Plan
3.4.1 Improving the algorithm
One of the most urgent objective for us is to improve the convergence efficiency of our
self-consistent calculation and prevent the divergence in the iteration process. We are look-
ing to borrow some experiences from the studies on DFT calculation [34, 35], in which
a lot of efforts have been devoted to obtain the self-consistency in density functional cal-
culations and some novel techniques have been introduced which are considered as the
better candidates for self-consistent calculations than the Broyden method we are currently
employing.