Abstract
This paper generalizesthe method of using the magnetic Bloch states to obtain flux-latticesolutions of the Bogoliubov-de Gennes equations.Developed by several groups, this method has had a limited applicability to the s-wave pairingand to high-κ materials where thespatial variation of the magnetic field may be neglected.It is extended here so that:(i) an arbitrary anisotropic pairing can be handled;(ii) the effect of the screening current can be taken into account.Hence the formalism will enable a fully microscopicdescription of various Abrikosov lattices in clean type-II superconductorsover Hc1 ≤ H ≤ Hc2.