Casimir Interactions between Magnetic Flux Tubes in a Dense Lattice

We use the worldline numerics technique to study a cylindrically symmetric model of magnetic flux tubes in a dense lattice and the non-local Casimir forces acting between regions of magnetic flux. Within a superconductor the magnetic field is constrained within magnetic flux tubes and if the background magnetic field is on the order the quantum critical field strength, $B_k = frac{m^2}{e} = 4.4 times 10^{13}$ Gauss, the magnetic field is likely to vary rapidly on the scales where acs{QED} effects are important. In this paper, we construct a cylindrically symmetric toy model of a flux tube lattice in which the non-local influence of acs{QED} on neighbouring flux tubes is taken into account. We compute the effective action densities using the worldline numerics technique. The numerics predict a greater effective energy density in the region of the flux tube, but a smaller energy density in the regions between the flux tubes compared to a locally-constant-field approximation. We also compute the interaction energy between a flux tube and its neighbours as the lattice spacing is reduced from infinity. Because our flux tubes exhibit compact support, this energy is entirely non-local and predicted to be zero in local approximations such as the derivative expansion. This Casimir-Polder energy can take positive or negative values depending on the distance between the flux tubes, and it may cause the flux tubes in neutron stars to form bunches. In addition to the above results we also discuss two important subtleties of determining the statistical uncertainties within the worldline numerics technique and recommend a form of jackknife analysis.

Parallel Worldline Numerics: Implementation and Error Analysis

We give an overview of the worldline numerics technique, and discuss the parallel CUDA implementation of a worldline numerics algorithm. In the worldline numerics technique, we wish to generate an ensemble of representative closed-loop particle trajectories, and use these to compute an approximate average value for Wilson loops. We show how this can be done with a specific emphasis on cylindrically symmetric magnetic fields. The fine-grained, massive parallelism provided by the GPU architecture results in considerable speedup in computing Wilson loop averages. Furthermore, we give a brief overview of uncertainty analysis in the worldline numerics method. There are uncertainties from discretizing each loop, and from using a statistical ensemble of representative loops. The former can be minimized so that the latter dominates. However, determining the statistical uncertainties is complicated by two subtleties. Firstly, the distributions generated by the worldline ensembles are highly non-Gaussian, and so the standard error in the mean is not a good measure of the statistical uncertainty. Secondly, because the same ensemble of worldlines is used to compute the Wilson loops at different values of $T$ and $x_mathrm{ cm}$, the uncertainties associated with each computed value of the integrand are strongly correlated. We recommend a form of jackknife analysis which deals with both of these problems.