Phase transitions in the one-dimensional frustrated quantum XY model and Josephson-junction ladders

A one-dimensional quantum version of the frustrated XY (planar rotor) model is considered which can be physically realized as a ladder of Josephson-junctions at half a flux quantum per plaquette. This system undergoes a superconductor to insulator transition at zero temperature as a function of charging energy. The critical behavior is studied using a Monte Carlo transfer matrix applied to the path-integral representation of the model and a finite-size-scaling analysis of data on small system sizes. Depending on the ratio between the interchain and intrachain couplings the system can have single or double transitions which is consistent with the prediction that its critical behavior should be described by the two-dimensional classical XY-Ising model.Comment: 13 pages, Revtex, J. Appl. Phys. (to appear), Inpe-las-00

Phase transitions in a frustrated XY model with zig-zag couplings

We study a new generalized version of the square-lattice frustrated XY model where unequal ferromagnetic and antiferromagnetic couplings are arranged in a zig-zag pattern. The ratio between the couplings $\rho$ can be used to tune the system, continuously, from the isotropic square-lattice to the triangular-lattice frustrated XY model. The model can be physically realized as a Josephson-junction array with two different couplings, in a magnetic field corresponding to half-flux quanta per plaquette. Mean-field approximation, Ginzburg-Landau expansion and finite-size scaling of Monte Carlo simulations are used to study the phase diagram and critical behavior. Depending on the value of $\rho$, two separate transitions or a transition line in the universality class of the XY-Ising model, with combined $Z_2$ and U(1) symmetries, takes place. In particular, the phase transitions of the standard square-lattice and triangular-lattice frustrated XY models correspond to two different cuts through the same transition line. Estimates of the chiral ($Z_2$) critical exponents on this transition line deviate significantly from the pure Ising values, consistent with that along the critical line of the XY-Ising model. This suggests that a frustrated XY model or Josephson-junction array with a zig-zag coupling modulation can provide a physical realization of the XY-Ising model critical line.Comment: 11 pages, 9 figures, RevTex, to appear in Phys. Rev.

Phase-coherence threshold and vortex-glass state in diluted Josephson-junction arrays in a magnetic field

We study numerically the interplay of phase coherence and vortex-glass state in two-dimensional Josephson-junction arrays with average rational values of flux quantum per plaquette $f$ and random dilution of junctions. For $f=1/2$, we find evidence of a phase coherence threshold value $x_s$, below the percolation concentration of diluted junctions $x_p$, where the superconducting transition vanishes. For $x_s < x < x_p$ the array behaves as a zero-temperature vortex glass with nonzero linear resistance at finite temperatures. The zero-temperature critical currents are insensitive to variations in $f$ in the vortex glass region while they are strongly $f$ dependent in the phase coherent region.Comment: 6 pages, 4 figures, to appear in Phys. Rev.

Numerical Studies of the Two Dimensional XY Model with Symmetry Breaking Fields

We present results of numerical studies of the two dimensional XY model with four and eight fold symmetry breaking fields. This model has recently been shown to describe hydrogen induced reconstruction on the W(100) surface. Based on mean-field and renormalization group arguments,we first show how the interplay between the anisotropy fields can give rise to different phase transitions in the model. When the fields are compatible with each other there is a continuous phase transition when the fourth order field is varied from negative to positive values. This transition becomes discontinuous at low temperatures. These two regimes are separated by a multicritical point. In the case of competing four and eight fold fields, the first order transition at low temperatures opens up into two Ising transitions. We then use numerical methods to accurately locate the position of the multicritical point, and to verify the nature of the transitions. The different techniques used include Monte Carlo histogram methods combined with finite size scaling analysis, the real space Monte Carlo Renormalization Group method, and the Monte Carlo Transfer Matrix method. Our numerical results are in good agreement with the theoretical arguments.Comment: 29 pages, HU-TFT-94-36, to appear in Phys. Rev. B, Vol 50, November 1, 1994. A LaTeX file with no figure

Phase Transitions in the Two-Dimensional XY Model with Random Phases: a Monte Carlo Study

We study the two-dimensional XY model with quenched random phases by Monte Carlo simulation and finite-size scaling analysis. We determine the phase diagram of the model and study its critical behavior as a function of disorder and temperature. If the strength of the randomness is less than a critical value, $\sigma_{c}$, the system has a Kosterlitz-Thouless (KT) phase transition from the paramagnetic phase to a state with quasi-long-range order. Our data suggest that the latter exists down to T=0 in contradiction with theories that predict the appearance of a low-temperature reentrant phase. At the critical disorder $T_{KT}\rightarrow 0$ and for $\sigma > \sigma_{c}$ there is no quasi-ordered phase. At zero temperature there is a phase transition between two different glassy states at $\sigma_{c}$. The functional dependence of the correlation length on $\sigma$ suggests that this transition corresponds to the disorder-driven unbinding of vortex pairs.Comment: LaTex file and 18 figure

Nonlinear sliding friction of adsorbed overlayers on disordered substrates

We study the response of an adsorbed monolayer on a disordered substrate under a driving force using Brownian molecular-dynamics simulation. We find that the sharp longitudinal and transverse depinning transitions with hysteresis still persist in the presence of weak disorder. However, the transitions are smeared out in the strong disorder limit. The theoretical results here provide a natural explanation for the recent data for the depinning transition of Kr films on gold substrate.Comment: 8 pages, 8 figs, to appear in Phys. Rev.

Phase Diagram and Commensurate-Incommensurate Transitions in the Phase Field Crystal Model with an External Pinning Potential

We study the phase diagram and the commensurate-incommensurate transitions in a phase field model of a two-dimensional crystal lattice in the presence of an external pinning potential. The model allows for both elastic and plastic deformations and provides a continuum description of lattice systems, such as for adsorbed atomic layers or two-dimensional vortex lattices. Analytically, a mode expansion analysis is used to determine the ground states and the commensurate-incommensurate transitions in the model as a function of the strength of the pinning potential and the lattice mismatch parameter. Numerical minimization of the corresponding free energy shows good agreement with the analytical predictions and provides details on the topological defects in the transition region. We find that for small mismatch the transition is of first-order, and it remains so for the largest values of mismatch studied here. Our results are consistent with results of simulations for atomistic models of adsorbed overlayers

Numerical studies of the 2 and 3D gauge glass at low temperature

We report results from Monte Carlo simulations of the two- and three-dimensional gauge glass at low temperature using parallel tempering Monte Carlo. In two dimensions, we find strong evidence for a zero-temperature transition. By means of finite-size scaling, we determine the stiffness exponent theta = -0.39 +/- 0.03. In three dimensions, where a finite-temperature transition is well established, we find theta = 0.27 +/- 0.01, compatible with recent results from domain-wall renormalization group studies.Comment: 3 pages, 3 figures. Proceedings of "2002 MMM Conference", Tampa, F

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

Two phase transitions in the fully frustrated XYXY model

The fully frustrated $XY$ model on a square lattice is studied by means of Monte Carlo simulations. A Kosterlitz-Thouless transition is found at $T_{\rm KT} \approx 0.446$, followed by an ordinary Ising transition at a slightly higher temperature, $T_c \approx 0.452$. The non-Ising exponents reported by others, are explained as a failure of finite size scaling due to the screening length associated with the nearby Kosterlitz-Thouless transition.Comment: REVTEX file, 8 pages, 5 figures in uuencoded postscrip