Quantization and simulation of Born-Infeld non-linear electrodynamics on a lattice

Born-Infeld non-linear electrodynamics arises naturally as a field theory description of the dynamics of strings and branes. Most analyses of this theory have been limited to studying it as a classical field theory. We quantize this theory on a Euclidean 4-dimensional space-time lattice and determine its properties using Monte-Carlo simulations. The electromagnetic field around a static point charge is measured using Luscher-Weisz methods to overcome the sign problem associated with the introduction of this charge. The D field appears identical to that of Maxwell QED. However, the E field is enhanced by quantum fluctuations, while still showing the short distance screening observed in the classical theory. In addition, whereas for the classical theory, the screening increases without bound as the non-linearity increases, the quantum theory approaches a limiting conformal field theory.Comment: 24 pages, 10 figures. Latex with postscript figure

Evidence for O(2) universality at the finite temperature transition for lattice QCD with 2 flavours of massless staggered quarks

We simulate lattice QCD with 2 flavours of massless quarks on lattices of temporal extent N_t=8, to study the finite temperature transition from hadronic matter to a quark-gluon plasma. A modified action which incorporates an irrelevant chiral 4-fermion interaction is used, which allows simulations at zero quark mass. We obtain excellent fits of the chiral condensates to the magnetizations of a 3-dimensional O(2) spin model on lattices small enough to model the finite size effects. This gives predictions for correlation lengths and chiral susceptibilities from the corresponding spin-model quantities. These are in good agreement with our measurements over the relevant range of parameters. Binder cumulants are measured, but the errors are too large to draw definite conclusions. From the properties of the O(2) spin model on the relatively small lattices with which we fit our `data', we can see why earlier attempts to fit staggered lattice data to leading-order infinite-volume scaling functions, as well as finite size scaling studies, failed and led to erroneous conclusions.Comment: 27 pages, Latex with 10 postscript figures. Some of the discussions have been expanded to satisfy a referee. Typographical errors were correcte

The RHMC algorithm for theories with unknown spectral bounds

The Rational Hybrid Monte Carlo (RHMC) algorithm extends the Hybrid Monte Carlo algorithm for lattice QCD simulations to situations involving fractional powers of the determinant of the quadratic Dirac operator. This avoids the updating increment ($dt$) dependence of observables which plagues the Hybrid Molecular-dynamics (HMD) method. The RHMC algorithm uses rational approximations to fractional powers of the quadratic Dirac operator. Such approximations are only available when positive upper and lower bounds to the operator's spectrum are known. We apply the RHMC algorithm to simulations of 2 theories for which a positive lower spectral bound is unknown: lattice QCD with staggered quarks at finite isospin chemical potential and lattice QCD with massless staggered quarks and chiral 4-fermion interactions ($\chi$QCD). A choice of lower bound is made in each case, and the properties of the RHMC simulations these define are studied. Justification of our choices of lower bounds is made by comparing measurements with those from HMD simulations, and by comparing different choices of lower bounds.Comment: Latex(Revtex 4) 25 pages, 8 postscript figure

Simulations of a Scintillator Compton Gamma Imager for Safety and Security

We are designing an all-scintillator Compton gamma imager for use in security investigations and remediation actions involving radioactive threat material. To satisfy requirements for a rugged and portable instrument, we have chosen solid scintillator for the active volumes of both the scatter and absorber detectors. Using the BEAMnrc/EGSnrc Monte Carlo simulation package, we have constructed models using four different materials for the scatter detector: LaBr_3, NaI, CaF_2 and PVT. We have compared the detector performances using angular resolution, efficiency, and image resolution. We find that while PVT provides worse performance than that of the detectors based entirely on inorganic scintillators, all of the materials investigated for the scatter detector have the potential to provide performance adequate for our purposes.Comment: Revised text and figures, Presented at SORMA West 2008, Published in IEEE Transactions on Nuclear Scienc

Density of states and Fisher's zeros in compact U(1) pure gauge theory

We present high-accuracy calculations of the density of states using multicanonical methods for lattice gauge theory with a compact gauge group U(1) on 4^4, 6^4 and 8^4 lattices. We show that the results are consistent with weak and strong coupling expansions. We present methods based on Chebyshev interpolations and Cauchy theorem to find the (Fisher's) zeros of the partition function in the complex beta=1/g^2 plane. The results are consistent with reweighting methods whenever the latter are accurate. We discuss the volume dependence of the imaginary part of the Fisher's zeros, the width and depth of the plaquette distribution at the value of beta where the two peaks have equal height. We discuss strategies to discriminate between first and second order transitions and explore them with data at larger volume but lower statistics. Higher statistics and even larger lattices are necessary to draw strong conclusions regarding the order of the transition.Comment: 14 pages, 16 figure