Comment on "Critical Dynamics of a Vortex-Loop Model for the Superconducting Transition"

Recently, Aji and Goldenfeldt [Phys. Rev. Lett. 87, 197003 (2001), cond-mat/0105622] put forward an explanation for the value of the dynamic critical exponent z observed in certain Monte Carlo simulations of the superconducting phase transition in zero magnetic field. In this Comment, we point out that their analysis is based on incorrect assumptions regarding the scaling dimension of the vortex density.Comment: 1 page, no figure

Improving the efficiency of extended ensemble simulations: The accelerated weight histogram method

We propose a method for efficient simulations in extended ensembles, useful, e.g., for the study of problems with complex energy landscapes and for free energy calculations. The main difficulty in such simulations is the estimation of the a priori unknown weight parameters needed to produce flat histograms. The method combines several complementary techniques, namely, a Gibbs sampler for the parameter moves, a reweighting procedure to optimize data use, and a Bayesian update allowing for systematic refinement of the free energy estimate. In a certain limit the scheme reduces to the 1/t algorithm of B.E. Belardinelli and V.D. Pereyra [Phys. Rev. E 75, 046701 (2007)]. The performance of the method is studied on the two-dimensional Ising model, where comparison with the exact free energy is possible, and on an Ising spin glass.Comment: 5 page

Influence of vortices and phase fluctuations on thermoelectric transport properties of superconductors in a magnetic field

We study heat transport and thermoelectric effects in two-dimensional superconductors in a magnetic field. These are modeled as granular Josephson-junction arrays, forming either regular or random lattices. We employ two different models for the dynamics, relaxational model-A dynamics or resistively and capacitively shunted Josephson junction (RCSJ) dynamics. We derive expressions for the heat current in these models, which are then used in numerical simulations to calculate the heat conductivity and the Nernst coefficient for different temperatures and magnetic fields. At low temperatures and zero magnetic field the heat conductivity in the RCSJ model is calculated analytically from a spin wave approximation, and is seen to have an anomalous logarithmic dependence on the system size, and also to diverge in the completely overdamped limit C -> 0. From our simulations we find at low magnetic fields that the Nernst signal displays a characteristic "tilted hill" profile similar to experiments and a non-monotonic temperature dependence of the heat conductivity. We also investigate the effects of granularity and randomness, which become important for higher magnetic fields. In this regime geometric frustration strongly influences the results in both regular and random systems and leads to highly non-trivial magnetic field dependencies of the studied transport coefficients

Evidence of many thermodynamic states of the three-dimensional Ising spin glass

We present a large-scale simulation of the three-dimensional Ising spin glass with Gaussian disorder to low temperatures and large sizes using optimized population annealing Monte Carlo. Our primary focus is investigating the number of pure states regarding a controversial statistic, characterizing the fraction of centrally peaked disorder instances, of the overlap function order parameter. We observe that this statistic is subtly and sensitively influenced by the slight fluctuations of the integrated central weight of the disorder-averaged overlap function, making the asymptotic growth behaviour very difficult to identify. Modified statistics effectively reducing this correlation are studied and essentially monotonic growth trends are obtained. The effect of temperature is also studied, finding a larger growth rate at a higher temperature. Our state-of-the-art simulation and variance reduction data analysis suggest that the many pure state picture is most likely and coherent.Comment: 8 pages, 5 figure

Modeling and simulations of quantum phase slips in ultrathin superconducting wires

We study quantum phase slips (QPS) in ultrathin superconducting wires. Starting from an effective one-dimensional microscopic model, which includes electromagnetic fluctuations, we map the problem to a (1+1)-dimensional gas of interacting instantons. We introduce a method to calculate the tunneling amplitude of quantum phase slips directly from Monte Carlo simulations. This allows us to go beyond the dilute instanton gas approximation and study the problem without any limitations of the density of QPS. We find that the tunneling amplitude shows a characteristic scaling behavior near the superconductor-insulator transition. We also calculate the voltage-charge relation of the insulating state, which is the dual of the Josephson current-phase relation in ordinary superconducting weak links. This evolves from a sinusoidal form in the regime of dilute QPS to more exotic shapes for higher QPS densities, where interactions are important.Comment: 12 pages, 11 figure