Computation of Equilibrium Vortex Structures for Type-II Superconductors

Abstract

A model of the vortex structure and configuration for homogeneous type-II superconductors can be deter mined by minimization of the Ginzburg-Landau free energy functional. A generalization of this formulation, the Lawrence-Doniach model, can be used for layered systems. A finite-dimensional approximation of this functional leads to a very large, sparse optimization problem. A number of algorithms have been used to attempt to solve this problem. The most successful of these is a damped, inexact Newton algorithm. Its com putational kernel is the iterative solution of a large, sparse, linear system that is nearly singular. We present computational results for solving these three- dimensional problems on the Intel DELTA. Using a general-purpose, scalable, iterative solver based on the preconditioned conjugate gradient method, we ob tained sustained computational rates of up to 4.26 GFLOPS on 512 processors. We find improvements of over a factor of 100 in the total execution time when compared with the same problem run on the CRAY-2. This computational tool has given us the first compu tational view of three-dimensional vortex phenomena such as vortex locking.

Get full access to this article

View all access options for this article.

References

Abrikosov, A. 1988. Fundamentals of the theory of metals. Amsterdam: North-Holland.

Bednorz, G. , and Müller, A. 1986. Possible high Tc superconductivity in Ba-La-Cu-O system. Z. Phys. B 64:189-193.

Coleman, T.F. , and Moré, J.J. 1983. Estimation of sparse Jacobian matrices and graph coloring problems. SIAM J. Numer. Anal. 20:187-209.

Doria, M.M. , Gubernatis, J.E. , and Rainer, D. 1990. On solving the Ginzburg-Landau equations by simulated annealing. Phys. Rev. B 41:6335-6340.

Garner, J. 1992. Unpublished information. Material Sciences Division. Argonne, Ill.: Argonne National Laboratory.

Garner, J. , Spanbauer, M. , Benedek, R. , Strand-burg, K. , Wright, S. , and Plassmann, P. 1992. Critical fields of Josephsoncoupled superconducting multilayers . Phys. Rev. B 45:7973-7983.

Golub, G.H. , and Van Loan, C.F. 1983. Matrix computations. Baltimore : Johns Hopkins University Press.

Gropp, W. 1992. Unpublished information. Mathematics and Computer Sciences Division. Argonne, Ill.: Argonne National Laboratory .

Jones, M.T. , and Plassmann, P.E. 1991a. An efficient parallel iterative solver for large sparse linear systems. In Proceedings of the IMA Workshop on Sparse Matrix Computations: Graph Theory Issues & Algorithms Minneapolis: University of Minnesota .

10.

—. 1991b. Scalable iterative solution of sparse linear systems . Preprint MCS-P277-1191. Mathematics and Computer Science Division. Argonne, Ill.: Argonne National Laboratory .

11.

—. 1992a. BlockSolve v1.0: scalable library software for the parallel solution of sparse linear systems. ANL Report ANL-92/46. Mathematics and Computer Science Division. Argonne, Ill.: Argonne National Laboratory.

12.

—. 1992b. The effect of many-color orderings on the convergence of iterative methods. In Proceedings of the Copper Mountain Conference on Iterative Methods Philadelphia: SIAM LA-SIG.

13.

—. 1993. A parallel graph coloring heuristic. SIAM J. Sci. Statist. Comput. 14: 654-669.

14.

Lawrence, W.E. , and Doniach, S. 1971. Theory of layer structure superconductors. In Proc. 12th Int'l Conf. on Low Temp. Physics, edited by E. Kanda. Kyoto: Academic Press of Japan , pp. 361— 362.

15.

Liu, D.C. , and Nocedal, J. 1990. On the limited-memory BFGS method for large scale optimization . Math, Programming 45:503-528.

16.

Liu, L. , Kouvel, J. , and Brun, T. 1991. Directional vortex-flux bifurcation in polycrystalline YBa2Cu3O7 rotated in a magnetic field. Phys. Rev. B 43:7859-7865.

17.

Liu, L. , Kouvel, J. , and Brun, T. 1992. Anisotropic-vortex-flux production and retention in grain-oriented YBa2Cu3O7. Phys. Rev. B 45:3054-3058.

18.

Manteuffel, T.A. 1980. An incomplete factorization technique for positive definite linear systems. Math. Comp. 34:473-497.

19.

Pothen, A. , Simon, H. , and Liou, K.-P. 1990. Partitioning sparse matrices with eigenvectors of graphs . SIAM J. Matrix Anal. 11:430-452.

20.

Schreiber, R. , and Tang, W.-P. 1982. Vectorizing the conjugate gradient method. Unpublished information. Department of Computer Science. Stanford, Calif.: Stanford University.

21.

Vavasis, S. 1991. Automatic domain partitioning in three dimensions . SIAM J. Sci. Statist. Comput. 12:950-970.

22.

Wang, Z.D. , and Hu, C.R. 1991. Numerical relaxation approach for solving the general Ginzburg-Landau equations for type-II superconductors. Phys. Rev. B 44: 11918-11923.

23.

Wu, M. , Ashburn, J. , Torng, C. , Hor, P. , Meng, R. , Gao, L. , Huang, Z. , Wang, Y. , and Chu, C. 1987. Superconductivity at 93 K in a new mixed-phase Y-Ba-Cu-O compound system at ambient pressure. Phys. Rev. Lett. 58:908-910.