Numerical studies of the phase diagram of layered type II   superconductors in a magnetic field
  We report on simulations of layered superconductors using the Lawrence-Doniach model in the framework of the lowest Landau level approximation. We find a first order phase transition with a $B(T)$ dependence which agrees very well with the experimental ``melting'' line in YBaCuO. The transition is not associated with vortex lattice melting, but separates two vortex liquid states characterised by different degrees of short-range crystalline order and different length scales of correlations between vortices in different layers. The transition line ends at a critical end-point at low fields. We find the magnetization discontinuity and the location of the lower critical magnetic field to be in good agreement with experiments in YBaCuO. Length scales of order parameter correlations parallel and perpendicular to the magnetic field increase exponentially as 1/T at low temperatures. The dominant relaxation time scales grow roughly exponentially with these correlation lengths. We find that the first order phase transition persists in the presence of weak random point disorder but can be suppressed entirely by strong disorder. No vortex glass or Bragg glass state is found in the presence of disorder. The consistency of our numerical results with various experimental features in YBaCuO, including the dependence on anisotropy, and the temperature dependence of the structure factor at the Bragg peaks in neutron scattering experiments is demonstrated. 
A. K. Kienappel and M. A. Moore