Vectorized multigrid Poisson solver for the CDC CYBER 205

The full multigrid (FMG) method is applied to the two dimensional Poisson equation with Dirichlet boundary conditions. This has been chosen as a relatively simple test case for examining the efficiency of fully vectorizing of the multigrid method. Data structure and programming considerations and techniques are discussed, accompanied by performance details

Ginzburg-Landau Vortex Lattice in Superconductor Films of Finite Thickness

The Ginzburg-Landau equations are solved for ideally periodic vortex lattices in superconducting films of arbitrary thickness in a perpendicular magnetic field. The order parameter, current density, magnetic moment, and the 3-dimensional magnetic field inside and outside the film are obtained in the entire ranges of the applied magnetic field, Ginzburg Landau parameter kappa, and film thickness. The superconducting order parameter varies very little near the surface (by about 0.01) and the energy of the film surface is small. The shear modulus c66 of the triangular vortex lattice in thin films coincides with the bulk c66 taken at large kappa. In thin type-I superconductor films with kappa < 0.707, c66 can be positive at low fields and negative at high fields.Comment: 12 pages including 14 Figures, corrected, Fig.14 added, appears in Phys. Rev. B 71, issue 1 (2005

Equivalence of the Falicov-Kimball and Brandt-Mielsch forms for the free energy of the infinite-dimensional Falicov-Kimball model

Falicov and Kimball proposed a real-axis form for the free energy of the Falicov-Kimball model that was modified for the coherent potential approximation by Plischke. Brandt and Mielsch proposed an imaginary-axis form for the free energy of the dynamical mean field theory solution of the Falicov-Kimball model. It has long been known that these two formulae are numerically equal to each other; an explicit derivation showing this equivalence is presented here.Comment: 4 pages, 1 figure, typeset with ReVTe

Theory of Type-II Superconductors with Finite London Penetration Depth

Previous continuum theory of type-II superconductors of various shapes with and without vortex pinning in an applied magnetic field and with transport current, is generalized to account for a finite London penetration depth lambda. This extension is particularly important at low inductions B, where the transition to the Meissner state is now described correctly, and for films with thickness comparable to or smaller than lambda. The finite width of the surface layer with screening currents and the correct dc and ac responses in various geometries follow naturally from an equation of motion for the current density in which the integral kernel now accounts for finite lambda. New geometries considered here are thick and thin strips with applied current, and `washers', i.e. thin film squares with a slot and central hole as used for SQUIDs.Comment: 14 pages, including 15 high-resolution figure

Analytic Solution for the Critical State in Superconducting Elliptic Films

A thin superconductor platelet with elliptic shape in a perpendicular magnetic field is considered. Using a method originally applied to circular disks, we obtain an approximate analytic solution for the two-dimensional critical state of this ellipse. In the limits of the circular disk and the long strip this solution is exact, i.e. the current density is constant in the region penetrated by flux. For ellipses with arbitrary axis ratio the obtained current density is constant to typically 0.001, and the magnetic moment deviates by less than 0.001 from the exact value. This analytic solution is thus very accurate. In increasing applied magnetic field, the penetrating flux fronts are approximately concentric ellipses whose axis ratio b/a < 1 decreases and shrinks to zero when the flux front reaches the center, the long axis staying finite in the fully penetrated state. Analytic expressions for these axes, the sheet current, the magnetic moment, and the perpendicular magnetic field are presented and discussed. This solution applies also to superconductors with anisotropic critical current if the anisotropy has a particular, rather realistic form.Comment: Revtex file and 13 postscript figures, gives 10 pages of text with figures built i

Critical State in Thin Anisotropic Superconductors of Arbitrary Shape

A thin flat superconductor of arbitrary shape and with arbitrary in-plane and out-of-plane anisotropy of flux-line pinning is considered, in an external magnetic field normal to its plane. It is shown that the general three-dimensional critical state problem for this superconductor reduces to the two-dimensional problem of an infinitely thin sample of the same shape but with a modified induction dependence of the critical sheet current. The methods of solving the latter problem are well known. This finding thus enables one to study the critical states in realistic samples of high-Tc superconductors with various types of anisotropic flux-line pinning. As examples, we investigate the critical states of long strips and rectangular platelets of high-Tc superconductors with pinning either by the ab-planes or by extended defects aligned with the c-axis.Comment: 13 pages including 13 figure files in the tex

Properties of the Ideal Ginzburg-Landau Vortex Lattice

The magnetization curves M(H) for ideal type-II superconductors and the maximum, minimum, and saddle point magnetic fields of the vortex lattice are calculated from Ginzburg-Landau theory for the entire ranges of applied magnetic fields Hc1 <= H < Hc2 or inductions 0 <= B < Hc2 and Ginzburg-Landau parameters sqrt(1/2) <= kappa <= 1000. Results for the triangular and square flux-line lattices are compared with the results of the circular cell approximation. The exact magnetic field B(x,y) and magnetization M(H, kappa) are compared with often used approximate expressions, some of which deviate considerably or have limited validity. Useful limiting expressions and analytical interpolation formulas are presented.Comment: 11 pages, 8 figure

Vector screening masses in the quark-gluon plasma and their physical significance

Static and non-static thermal screening states that couple to the conserved vector current are investigated in the high-temperature phase of QCD. Their masses and couplings to the current are determined at weak coupling, as well as using two-flavor lattice QCD simulations. A consistent picture emerges from the comparison, providing evidence that non-static Matsubara modes can indeed be treated perturbatively. We elaborate on the physical significance of the screening masses.Comment: 4 pages, 3 figures. Submitted as a contribution to the proceedings of the Quark Matter 2014 conference (talk given by H. Meyer