“…These structural ordering compete with each other and the resulted vortex arrangement is determined by the interplay between the geometry-induced confining potential and the mutual repulsive interaction between vortices [1,2,3,4]. For large vorticities, a piece(s) of the triangular vortex lattice remains near the disk center because of densely packed vortices [3,5], whereas for small vorticities, vortices are arranged themselves into regular polygons situated on concentric circular rings from the disk center [1,2,4]. This is called "vortex shells", of which configurations have been visualized by some experimental studies on mesoscopic superconducting disks of Nb and amorphuous MoGe thin films [6,7].…”