Numerical Study of Order in a Gauge Glass Model

The XY model with quenched random phase shifts is studied by a T=0 finite size defect energy scaling method in 2d and 3d. The defect energy is defined by a change in the boundary conditions from those compatible with the true ground state configuration for a given realization of disorder. A numerical technique, which is exact in principle, is used to evaluate this energy and to estimate the stiffness exponent $\theta$. This method gives $\theta = -0.36\pm0.013$ in 2d and $\theta = +0.31\pm 0.015$ in 3d, which are considerably larger than previous estimates, strongly suggesting that the lower critical dimension is less than three. Some arguments in favor of these new estimates are given.