Lattice Gas Dynamics; Application to Driven Vortices in Two Dimensional Superconductors

A continuous time Monte Carlo lattice gas dynamics is developed to model driven steady states of vortices in two dimensional superconducting networks. Dramatic differences are found when compared to a simpler Metropolis dynamics. Subtle finite size effects are found at low temperature, with a moving smectic that becomes unstable to an anisotropic liquid on sufficiently large length scales.Comment: 5 pages, 4 figure

Glassiness Vs. Order in Densely Frustrated Josephson Arrays

We carry out extensive Monte Carlo simulations on the Coulomb gas dual to the uniformly frustrated two dimensional XY model, for a sequence of frustrations f converging to the irraltional (3-sqrt 5)/2. We find in these systems a sharp first order equilibrium phase transition to an ordered vortex structure at a T_c which varies only slightly with f. This ordered vortex structure remains in general phase incoherent until a lower pinning transition T_p(f) that varies with f. We argue that the glassy behaviors reported for this model in earlier simulations are dynamic effects.Comment: 4 pages, 4 eps figure

Spin and chiral orderings of frustrated quantum spin chains

Ordering of frustrated S=1/2 and 1 XY and Heisenberg spin chains with the competing nearest- and next-nearest-neighbor antiferromagnetic couplings is studied by exact diagonalization and density-matrix renormalization-group methods. It is found that the S=1 XY chain exhibits both gapless and gapped `chiral' phases characterized by the spontaneous breaking of parity, in which the long-range order parameter is a chirality, $\kappa_i = S_i^xS_{i+1}^y-S_i^yS_{i+1}^x$, whereas the spin correlation decays either algebraically or exponentially. Such chiral phases are not realized in the S=1/2 XY chain nor in the Heisenberg chains.Comment: 4 pages, 5 EPS-figures, LaTeX(RevTeX),to appear in J.Phys.Soc.Japa

The Effect of Columnar Disorder on the Superconducting Transition of a Type-II Superconductor in Zero Applied Magnetic Field

We investigate the effect of random columnar disorder on the superconducting phase transition of a type-II superconductor in zero applied magnetic field using numerical simulations of three dimensional XY and vortex loop models. We consider both an unscreened model, in which the bare magnetic penetration length is approximated as infinite, and a strongly screened model, in which the magnetic penetration length is of order the vortex core radius. We consider both equilibrium and dynamic critical exponents. We show that, as in the disorder free case, the equilibrium transitions of the unscreened and strongly screened models lie in the same universality class, however scaling is now anisotropic. We find for the correlation length exponent $\nu=1.2\pm 0.1$, and for the anisotropy exponent $\zeta=1.3\pm 0.1$. We find different dynamic critical exponents for the unscreened and strongly screened models.Comment: 30 pages 12 ps figure

Positional Disorder (Random Gaussian Phase Shifts) in the Fully Frustrated Josephson Junction Array (2D XY Model)

We consider the effect of positional disorder on a Josephson junction array with an applied magnetic field of f=1/2 flux quantum per unit cell. This is equivalent to the problem of random Gaussian phase shifts in the fully frustrated 2D XY model. Using simple analytical arguments and numerical simulations, we present evidence that the ground state vortex lattice of the pure model becomes disordered, in the thermodynamic limit, by any amount of positional disorder.Comment: 4 pages, 4 eps figures embedde

Continuous Time Monte Carlo and Spatial Ordering in Driven Lattice Gases: Application to Driven Vortices in Periodic Superconducting Networks

We consider the two dimensional (2D) classical lattice Coulomb gas as a model for magnetic field induced vortices in 2D superconducting networks. Two different dynamical rules are introduced to investigate driven diffusive steady states far from equilibrium as a function of temperature and driving force. The resulting steady states differ dramatically depending on which dynamical rule is used. We show that the commonly used driven diffusive Metropolis Monte Carlo dynamics contains unphysical intrinsic randomness that destroys the spatial ordering present in equilibrium (the vortex lattice) over most of the driven phase diagram. A continuous time Monte Carlo (CTMC) is then developed, which results in spatially ordered driven states at low temperature in finite sized systems. We show that CTMC is the natural discretization of continuum Langevin dynamics, and argue that it gives the correct physical behavior when the discrete grid represents the minima of a periodic potential. We use detailed finite size scaling methods to analyze the spatial structure of the steady states. We find that finite size effects can be subtle and that very long simulation times can be needed to arrive at the correct steady state. For particles moving on a triangular grid, we find that the ordered moving state is a transversely pinned smectic that becomes unstable to an anisotropic liquid on sufficiently large length scales. For particles moving on a square grid, the moving state is a similar smectic at large drives, but we find evidence for a possible moving solid at lower drives. We find that the driven liquid on the square grid has long range hexatic order, and we explain this as a specifically non-equilibrium effect. We show that, in the liquid, fluctuations are diffusive in both the transverse and longitudinal directions.Comment: 29 pages, 35 figure

Vortex lattce melting in 2D superconductors and Josephson arrays

Monte Carlo simulations of 2D vortex lattice melting in a thin superconducting film (or alternatively an array of Josephson junctions) are performed in the London limit. Finite size scaling analyses are used to make a detailed test of the dislocation mediated melting theory of KTNHY. We find that the melting transition is weakly first order, with a jump in the shear modulus very close to that predicted by the KTNHY theory. No hexatic liquid phase is found.Comment: 12 pages, 4 figures (available on request from [email protected]), REVTEX [we revise our conclusion on the order of the melting transition from second to first order - new figure 4 added

Heterogeneous Dynamics, Marginal Stability and Soft Modes in Hard Sphere Glasses

In a recent publication we established an analogy between the free energy of a hard sphere system and the energy of an elastic network [1]. This result enables one to study the free energy landscape of hard spheres, in particular to define normal modes. In this Letter we use these tools to analyze the activated transitions between meta-bassins, both in the aging regime deep in the glass phase and near the glass transition. We observe numerically that structural relaxation occurs mostly along a very small number of nearly-unstable extended modes. This number decays for denser packing and is significantly lowered as the system undergoes the glass transition. This observation supports that structural relaxation and marginal modes share common properties. In particular theoretical results [2, 3] show that these modes extend at least on some length scale $l^*\sim (\phi_c-\phi)^{-1/2}$ where $\phi_c$ corresponds to the maximum packing fraction, i.e. the jamming transition. This prediction is consistent with very recent numerical observations of sheared systems near the jamming threshold [4], where a similar exponent is found, and with the commonly observed growth of the rearranging regions with compression near the glass transition.Comment: 6 pages, improved versio

Current-voltage scaling of a Josephson-junction array at irrational frustration

Numerical simulations of the current-voltage characteristics of an ordered two-dimensional Josephson junction array at an irrational flux quantum per plaquette are presented. The results are consistent with an scaling analysis which assumes a zero temperature vortex glass transition. The thermal correlation length exponent characterizing this transition is found to be significantly different from the corresponding value for vortex-glass models in disordered two-dimensional superconductors. This leads to a current scale where nonlinearities appear in the current-voltage characteristics decreasing with temperature $T$ roughly as $T^2$ in contrast with the $T^3$ behavior expected for disordered models.Comment: RevTex 3.0, 12 pages with Latex figures, to appear in Phys. Rev. B 54, Rapid. Com

Kinetic Inductance of Josephson Junction Arrays: Dynamic and Equilibrium Calculations

We show analytically that the inverse kinetic inductance $L^{-1}$ of an overdamped junction array at low frequencies is proportional to the admittance of an inhomogeneous equivalent impedance network. The $ij^{th}$ bond in this equivalent network has an inverse inductance $J_{ij}\cos(\theta_i^0-\theta_j^0-A_{ij})$, where $J_{ij}$ is the Josephson coupling energy of the $ij^{th}$ bond, $\theta_i^0$ is the ground-state phase of the grain $i$, and $A_{ij}$ is the usual magnetic phase factor. We use this theorem to calculate $L^{-1}$ for square arrays as large as $180\times 180$. The calculated $L^{-1}$ is in very good agreement with the low-temperature limit of the helicity modulus $\gamma$ calculated by conventional equilibrium Monte Carlo techniques. However, the finite temperature structure of $\gamma$, as a function of magnetic field, is \underline{sharper} than the zero-temperature $L^{-1}$, which shows surprisingly weak structure. In triangular arrays, the equilibrium calculation of $\gamma$ yields a series of peaks at frustrations $f = \frac{1}{2}(1-1/N)$, where $N$ is an integer $\geq 2$, consistent with experiment.Comment: 14 pages + 6 postscript figures, 3.0 REVTe