From vortex molecules to the Abrikosov lattice in thin mesoscopic   superconducting disks
  Stable vortex states are studied in large superconducting thin disks (for numerical purposes we considered with radius R = 50 \xi). Configurations containing more than 700 vortices were obtained using two different approaches: the nonlinear Ginzburg-Landau (GL) theory and the London approximation. To obtain better agreement with results from the GL theory we generalized the London theory by including the spatial variation of the order parameter following Clem's ansatz. We find that configurations calculated in the London limit are also stable within the Ginzburg-Landau theory for up to ~ 230 vortices. For large values of the vorticity (typically, L > 100), the vortices are arranged in an Abrikosov lattice in the center of the disk, which is surrounded by at least two circular shells of vortices. A Voronoi construction is used to identify the defects present in the ground state vortex configurations. Such defects cluster near the edge of the disk, but for large L also grain boundaries are found which extend up to the center of the disk. 
L. R. E. Cabral, B. J. Baelus and F. M. Peeters