Event-driven Monte Carlo algorithm for general potentials

We extend the event-chain Monte Carlo algorithm from hard-sphere interactions to the micro-canonical ensemble (constant potential energy) for general potentials. This event-driven Monte Carlo algorithm is non-local, rejection-free, and allows for the breaking of detailed balance. The algorithm uses a discretized potential, but its running speed is asymptotically independent of the discretization. We implement the algorithm for the cut-off linear potential, and discuss its possible implementation directly in the continuum limit.Comment: 3 pages 4 figures, brief repor

Damage spreading and coupling in Markov chains

In this paper, we relate the coupling of Markov chains, at the basis of perfect sampling methods, with damage spreading, which captures the chaotic nature of stochastic dynamics. For two-dimensional spin glasses and hard spheres we point out that the obstacle to the application of perfect-sampling schemes is posed by damage spreading rather than by the survey problem of the entire configuration space. We find dynamical damage-spreading transitions deeply inside the paramagnetic and liquid phases, and show that critical values of the transition temperatures and densities depend on the coupling scheme. We discuss our findings in the light of a classic proof that for arbitrary Monte Carlo algorithms damage spreading can be avoided through non-Markovian coupling schemes.Comment: 6 pages, 8 figure

Hard-disk equation of state: First-order liquid-hexatic transition in two dimensions with three simulation methods

We report large-scale computer simulations of the hard-disk system at high densities in the region of the melting transition. Our simulations reproduce the equation of state, previously obtained using the event-chain Monte Carlo algorithm, with a massively parallel implementation of the local Monte Carlo method and with event-driven molecular dynamics. We analyze the relative performance of these simulation methods to sample configuration space and approach equilibrium. Our results confirm the first-order nature of the melting phase transition in hard disks. Phase coexistence is visualized for individual configurations via the orientational order parameter field. The analysis of positional order confirms the existence of the hexatic phase.Comment: 9 pages, 8 figures, 2 table

Two-step melting in two dimensions: First-order liquid-hexatic transition

Melting in two spatial dimensions, as realized in thin films or at interfaces, represents one of the most fascinating phase transitions in nature, but it remains poorly understood. Even for the fundamental hard-disk model, the melting mechanism has not been agreed on after fifty years of studies. A recent Monte Carlo algorithm allows us to thermalize systems large enough to access the thermodynamic regime. We show that melting in hard disks proceeds in two steps with a liquid phase, a hexatic phase, and a solid. The hexatic-solid transition is continuous while, surprisingly, the liquid-hexatic transition is of first-order. This melting scenario solves one of the fundamental statistical-physics models, which is at the root of a large body of theoretical, computational and experimental research.Comment: 5 pages, 4 figure

Risk factors for treatment failure in orthopedic device-related methicillin-resistant Staphylococcus aureus infection

The purpose of this study was to determine the clinical and microbiological risk factors for treatment failure of methicillin-resistant Staphylococcus aureus (MRSA) orthopedic device-related infection (ODRI). A retrospective cohort study of patients with MRSA ODRI who were treated at Geneva University Hospitals between 2000 and 2008 was undertaken. Stored MRSA isolates were retrieved for genetic characterization and determination of the vancomycin minimum inhibitory concentration (MIC). Fifty-two patients were included, of whom 23 (44%) had joint arthroplasty and 29 (56%) had osteosynthesis. All 41 of the retrieved MRSA isolates were susceptible to vancomycin (MIC ≤ 2mg/L) and 35 (85%) shared genetic characteristics of the South German clone (ST228). During a median follow-up of 391days (range, 4-2,922days), 18 patients (35%) experienced treatment failure involving MRSA persistence or recurrence. Microbiological factors such as infection with the predominant clone and a vancomycin MIC of 2mg/L were not associated with treatment failure. Using a Cox proportional hazards model, implant retention (hazard ratio [HR], 4.9; 95% confidence interval [CI], 1.3-18.2; P = 0.017) and single-agent antimicrobial therapy (HR, 4.4; 95% CI, 1.2-16.3; P = 0.025) were independent predictors of treatment failure after debridement. Therapy using a combination of antimicrobials should be considered for patients with MRSA ODRI, especially when implant removal is not feasibl

Internal avalanches in a pile of superconducting vortices

Using an array of miniature Hall probes, we monitored the spatiotemporal variation of the internal magnetic induction in a superconducting niobium sample during a slow sweep of external magnetic field. We found that a sizable fraction of the increase in the local vortex population occurs in abrupt jumps. The size distribution of these avalanches presents a power-law collapse on a limited range. In contrast, at low temperatures and low fields, huge avalanches with a typical size occur and the system does not display a well-defined macroscopic critical current.Comment: 5 pages including 5 figure