Duality in Shearing Rheology Near the Athermal Jamming Transition

We consider the rheology of soft-core frictionless disks in two dimensions in the neighborhood of the athermal jamming transition. From numerical simulations of bidisperse, overdamped, particles, we argue that the divergence of the viscosity below jamming is characteristic of the hard-core limit, independent of the particular soft-core interaction. We develop a mapping from soft-core to hard-core particles that recovers all the critical behavior found in earlier scaling analyses. Using this mapping we derive a duality relation that gives the exponent of the non-linear Herschel-Bulkley rheology above jamming in terms of the exponent of the diverging viscosity below jamming.Comment: 5 pages, 4 figures. Manuscript revisions: new title, additional text concerning connections to experiment, revised Fig. 4, other minor changes and clarifications in text. Conclusions remain essentially unchanged. Accepted for publication in Phys. Rev. Let

Kink-Antikink Unbinding Transition in the Two Dimensional Fully Frustrated XY Model

We carry out the first numerical simulations to directly confirm the existence of a kink-antikink unbinding transition, T_w, along Ising-like domain walls in the two dimensional fully frustrated XY model. We comment on the possible implications of kink-antikink unbinding for the bulk phase transitions of the model.Comment: 8 pages, 10 figures - expanded version of original submission as accepted by Phys. Rev. B. Some correction to the analysis lead to an increase in the value of the kink unbinding transition temperatur

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

Melting of icosahedral gold nanoclusters from molecular dynamics simulations

Molecular dynamics simulations show that gold clusters with about 600--3000 atoms crystallize into a Mackay icosahedron upon cooling from the liquid. A detailed surface analysis shows that the facets on the surface of the Mackay icosahedral gold clusters soften but do not premelt below the bulk melting temperature. This softening is found to be due to the increasing mobility of vertex and edge atoms with temperature, which leads to inter-layer and intra-layer diffusion, and a shrinkage of the average facet size, so that the average shape of the cluster is nearly spherical at melting.Comment: 40 pages, 27 figure

Thermodynamics of a Heavy Ion-Irradiated Superconductor: the Zero-Field Transition

Specific heat measurements show that the introduction of amorphous columnar defects considerably affects the transition from the normal to the superconducting state in zero magnetic field. Experimental results are compared to numerical simulations of the 3D XY model for both the pure system and the system containing random columnar disorder. The numerics reproduce the salient features of experiment, showing in particular that the specific heat peak changes from cusp-like to smoothly rounded when columnar defects are added. By considering the specific heat critical exponent alpha, we argue that such behavior is consistent with recent numerical work [Vestergren et al., PRB 70, 054508 (2004)] showing that the introduction of columnar defects changes the universality class of the transition.Comment: 4 pages, 2 figure

Observation of Ising-like critical fluctuations in frustrated Josephson junction arrays with modulated coupling energies

We report the results of ac sheet conductance measurements performed on fully frustrated square arrays of Josephson junctions whose coupling energy is periodically modulated in one of the principal lattice directions. Such systems are predicted to exhibit two distinct transitions: a low-temperature Ising-like transition triggered by the proliferation of domain walls and a high-temperature transition driven by the vortex unbinding mechanism of the Beresinskii-Kosterlitz-Thouless (BKT) theory. Both the superfluid and dissipative components of the conductance are found to exhibit features which unambiguously demonstrate the existence of a double transition whose properties are consistent with the Ising-BKT scenario.Comment: To be published in Physica C (Proceedings of the 2nd European Conference in School Format 'Vortex Matter in Superconductors'

Finite-Size-Scaling at the Jamming Transition: Corrections to Scaling and the Correlation Length Critical Exponent

We carry out a finite size scaling analysis of the jamming transition in frictionless bi-disperse soft core disks in two dimensions. We consider two different jamming protocols: (i) quench from random initial positions, and (ii) quasistatic shearing. By considering the fraction of jammed states as a function of packing fraction for systems with different numbers of particles, we determine the spatial correlation length critical exponent $\nu\approx 1$, and show that corrections to scaling are crucial for analyzing the data. We show that earlier numerical results yielding $\nu<1$ are due to the improper neglect of these corrections.Comment: 5 pages, 4 figures -- slightly revised version as accepted for Phys. Rev. E Rapid Communication

Patterns of Striped order in the Classical Lattice Coulomb Gas

We obtain via Monte Carlo simulations the low temperature charge configurations in the lattice Coulomb gas on square lattices for charge filling ratio $f$ in the range $1/3 < f < 1/2$. We find a simple regularity in the low temperature charge configurations which consist of a suitable periodic combination of a few basic striped patterns characterized by the existence of partially filled diagonal channels. In general there exist two separate transitions where the lower temperature transition ($T_p$) corresponds to the freezing of charges within the partially filled channels. $T_p$ is found to be sensitively dependent on $f$ through the charge number density $\nu = p_{1}/q_{1}$ within the channels.Comment: 4 pages, 8 figure

Boundary effects on the scaling of the superfluid density

We study numerically the influence of the substrate (boundary conditions) on the finite--size scaling properties of the superfluid density $\rho_s$ in superfluid films of thickness $H$ within the XY model employing the Monte Carlo method. Our results suggest that the jump $\rho_s H/T_c$ at the Kosterlitz--Thouless transition temperature $T_c$ depends on the boundary conditions.Comment: 2 pages, 1 Latex file, 1 postscript figure, 2 style file

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