Current-voltage characteristics of the two-dimensional XY model with Monte Carlo dynamics

Current-voltage characteristics and the linear resistance of the two-dimensional XY model with and without external uniform current driving are studied by Monte Carlo simulations. We apply the standard finite-size scaling analysis to get the dynamic critical exponent $z$ at various temperatures. From the comparison with the resistively-shunted junction dynamics, it is concluded that $z$ is universal in the sense that it does not depend on details of dynamics. This comparison also leads to the quantification of the time in the Monte Carlo dynamic simulation.Comment: 5 pages in two columns including 5 figures, to appear in PR

Nonperturbative QCD Vacuum Effects in Nonlocal Quark Dynamics

A straightforward calculation reveals the essentially nonlocal character of the leading heavy $Q\bar{Q}$ interaction arising from nonperturbative gluon field correlations in the model of a fluctuating QCD vacuum. In light of this quarkonium spin splitting ratio predictions which have supported the scalar confinement ansatz are reconsidered as a specific example of possible consequences for spectroscopy.Comment: Latex, 9 page

Collapse of the vortex-lattice inductance and shear modulus at the melting transition in untwinned YBa2Cu3O7\rm YBa_2Cu_3O_7

The complex resistivity $\hat{\rho}(\omega)$ of the vortex lattice in an untwinned crystal of 93-K $\rm YBa_2Cu_3O_7$ has been measured at frequencies $\omega/2\pi$ from 100 kHz to 20 MHz in a 2-Tesla field $\bf H\parallel c$, using a 4-probe RF transmission technique that enables continuous measurements versus $\omega$ and temperature $T$. As $T$ is increased, the inductance ${\cal L}_s(\omega) ={\rm Im} \hat{\rho}(\omega)/ \omega$ increases steeply to a cusp at the melting temperature $T_m$, and then undergoes a steep collapse consistent with vanishing of the shear modulus $c_{66}$. We discuss in detail the separation of the vortex-lattice inductance from the `volume' inductance, and other skin-depth effects. To analyze the spectra, we consider a weakly disordered lattice with a low pin density. Close fits are obtained to $\rho_1(\omega)$ over 2 decades in $\omega$. Values of the pinning parameter $\kappa$ and shear modulus $c_{66}$ obtained show that $c_{66}$ collapses by over 4 decades at $T_m$, whereas $\kappa$ remains finite.Comment: 11 pages, 8 figures, Phys. Rev. B, in pres

Dynamic Simulations of the Kosterlitz-Thouless Phase Transition

Based on the short-time dynamic scaling form, a novel dynamic approach is proposed to tackle numerically the Kosterlitz-Thouless phase transition. Taking the two-dimensional XY model as an example, the exponential divergence of the spatial correlation length, the transition temperature $T_{KT}$ and all critical exponents are computed. Compared with Monte Carlo simulations in equilibrium, we obtain data at temperatures nearer to $T_{KT}$.Comment: to appear in Phys. Rev. E in Rapid Communicatio

Numerical study of the strongly screened vortex glass model in an external field

The vortex glass model for a disordered high-T_c superconductor in an external magnetic field is studied in the strong screening limit. With exact ground state (i.e. T=0) calculations we show that 1) the ground state of the vortex configuration varies drastically with infinitesimal variations of the strength of the external field, 2) the minimum energy of global excitation loops of length scale L do not depend on the strength of the external field, however 3) the excitation loops themself depend sensibly on the field. From 2) we infer the absence of a true superconducting state at any finite temperature independent of the external field.Comment: 6 pages RevTeX, 5 eps-figures include

Application of a minimum cost flow algorithm to the three-dimensional gauge glass model with screening

We study the three-dimensional gauge glass model in the limit of strong screening by using a minimum cost flow algorithm, enabling us to obtain EXACT ground states for systems of linear size L<=48. By calculating the domain-wall energy, we obtain the stiffness exponent theta = -0.95+/-0.03, indicating the absence of a finite temperature phase transition, and the thermal exponent nu=1.05+/-0.03. We discuss the sensitivity of the ground state with respect to small perturbations of the disorder and determine the overlap length, which is characterized by the chaos exponent zeta=3.9+/-0.2, implying strong chaos.Comment: 4 pages RevTeX, 2 eps-figures include

Critical dynamics in the 2d classical XY-model: a spin dynamics study

Using spin-dynamics techniques we have performed large-scale computer simulations of the dynamic behavior of the classical three component XY-model (i.e. the anisotropic limit of an easy-plane Heisenberg ferromagnet), on square lattices of size up to 192^2, for several temperatures below, at, and above T_KT. The temporal evolution of spin configurations was determined numerically from coupled equations of motion for individual spins by a fourth order predictor-corrector method, with initial spin configurations generated by a hybrid Monte Carlo algorithm. The neutron scattering function S(q,omega) was calculated from the resultant space-time displaced spin-spin correlation function. Pronounced spin-wave peaks were found both in the in-plane and the out-of-plane scattering function over a wide range of temperatures. The in-plane scattering function S^xx also has a large number of clear but weak additional peaks, which we interpret to come from two-spin-wave scattering. In addition, we observed a small central peak in S^xx, even at temperatures well below the phase transition. We used dynamic finite size scaling theory to extract the dynamic critical exponent z. We find z=1.00(4) for all T <= T_KT, in excellent agreement with theoretical predictions, although the shape of S(q,omega) is not well described by current theory.Comment: 31 pages, LaTex, 13 figures (38 subfigures) included as eps-files, needs psfig, 260 K

Monte Carlo study of the Villain version of the fully frustrated XY model

The fully frustrated XY model with Villain interaction on a square lattice is studied by means of Monte Carlo simulations. On the basis of the universal jump condition it is argued that there are two distinct transitions in the model, corresponding to the loss of XY order and Z_2 order, respectively. The Kosterlitz-Thouless (KT) transition is analyzed by finite size scaling of the helicity modulus at lattices of size L = 32 through 128, giving T_KT = 0.8108(1). The vorticity-vorticity correlation function is used to determine two different characteristic lengths, the Z_2 correlation length \xi, and the screening length \lambda, associated with the KT transition and free vortices. The temperature dependence of \xi is examined in order to determine T_c and the correlation length exponent, \nu. The exponent is found to be consistent with the 2D Ising value, \nu = 1, and the obtained critical temperature is T_c = 0.8225(5). The determinations of both \xi and \nu are done carefully, first applying the techniques to the 2D Ising model, which serves as a convenient testing ground.Comment: 56 pages, 31 sub-figures, submitted to Phys Rev

Order parameter for two-dimensional critical systems with boundaries

Conformal transformations can be used to obtain the order parameter for two-dimensional systems at criticality in finite geometries with fixed boundary conditions on a connected boundary. To the known examples of this class (such as the disk and the infinite strip) we contribute the case of a rectangle. We show that the order parameter profile for simply connected boundaries can be represented as a universal function (independent of the criticality model) raised to the power eta/2. The universal function can be determined from the Gaussian model or equivalently a problem in two-dimensional electrostatics. We show that fitting the order parameter profile to the theoretical form gives an accurate route to the determination of eta. We perform numerical simulations for the Ising model and percolation for comparison with these analytic predictions, and apply this approach to the study of the planar rotor model.Comment: 10 pages, 14 figures. Revisions: Removed many typos, improved presentation of most of the figure

A Gauge Invariant Unitary Theory for Pion Photoproduction

A covariant, unitary and gauge invariant theory for pion photoproduction on a single nucleon is presented. To achieve gauge invariance at the operator level one needs to include both the $\pi N$ and $\gamma\pi N$ thresholds. The final amplitude can be written in terms of a distorted wave in the final $\pi N$ channel provided one includes additional diagrams to the standard Born term in which the photon is coupled to the final state pion and nucleon. These additional diagrams are required in order to satisfy gauge invariance.Comment: 4 pages, LaTeX, 1 figure as a separate uuencoded compressed tar fil