“…One may gain some physical insight into this difference by considering the different types of quartic trap, which is isotropic in [30], while extreme anisotropic in our model. Thus, compared with the results reported in [30], although the boundary between (b) and (c) is more sensitive to μ, our numerical results are, to some extent, similar to the results of [30], that is, the area of phase (c) increases with the increasing of the strength of μ. Finally, it is noteworthy that in both regions (b) and (c), the density distributions are stretched along y-axis direction, and we attribute this to the special external potential considered here.…”