Search

We have studied numerically the appearance of multiscaling behavior in the three-dimensional ferromagnetic Ising site diluted model, in the form of a multifractal distribution of the decay exponents for the spatial correlation functions at the critical temperature. We have computed the exponents of the long-distance decay of higher moments of the correlation function, up to the 10th power, by studying three different quantities: global susceptibilities, local susceptibilities and correlation functions. We have found very clear evidences for multiscaling behavior., Comment: 19 pages and 5 figures. Final version of the paper

Rejuvenation and memory, long considered the distinguishing features of spin glasses, have recently been proven to result from the growth of multiple length scales. This insight, enabled by simulations on the Janus~II supercomputer, has opened the door to a quantitative analysis. We combine numerical simulations with comparable experiments to introduce two coefficients that quantify memory. A third coefficient has been recently presented by Freedberg et al. We show that these coefficients are physically equivalent by studying their temperature and waiting-time dependence., Comment: 10 pages, 8 figures

We unveil the multifractal behavior of Ising spin glasses in their low-temperature phase. Using the Janus II custom-built supercomputer, the spin-glass correlation function is studied locally. Dramatic fluctuations are found when pairs of sites at the same distance are compared. The scaling of these fluctuations, as the spin-glass coherence length grows with time, is characterized through the computation of the singularity spectrum and its corresponding Legendre transform. A comparatively small number of site pairs controls the average correlation that governs the response to a magnetic field. We explain how this scenario of dramatic fluctuations (at length scales smaller than the coherence length) can be reconciled with the smooth, self-averaging behavior that has long been considered to describe spin-glass dynamics., Comment: 15 pages, 15 figures

Genome assembly, the process of reconstructing a long genetic sequence by aligning and merging short fragments, or reads, is known to be NP-hard, either as a version of the shortest common superstring problem or in a Hamiltonian-cycle formulation. That is, the computing time is believed to grow exponentially with the the problem size in the worst case. Despite this fact, high-throughput technologies and modern algorithms currently allow bioinformaticians to handle datasets of billions of reads. Using methods from statistical mechanics, we address this conundrum by demonstrating the existence of a phase transition in the computational complexity of the problem and showing that practical instances always fall in the 'easy' phase (solvable by polynomial-time algorithms). In addition, we propose a Markov-chain Monte Carlo method that outperforms common deterministic algorithms in the hard regime., Comment: 9 pages, 5 figures. Version accepted for publication in Phys. Rev. E

The extended principle of superposition has been a touchstone of spin glass dynamics for almost thirty years. The Uppsala group has demonstrated its validity for the metallic spin glass, CuMn, for magnetic fields $H$ up to 10 Oe at the reduced temperature $T_\mathrm{r}=T/T_\mathrm{g} = 0.95$, where $T_\mathrm{g}$ is the spin glass condensation temperature. For $H > 10$ Oe, they observe a departure from linear response which they ascribe to the development of non-linear dynamics. The thrust of this paper is to develop a microscopic origin for this behavior by focusing on the time development of the spin glass correlation length, $\xi(t,t_\mathrm{w};H)$. Here, $t$ is the time after $H$ changes, and $t_\mathrm{w}$ is the time from the quench for $T>T_\mathrm{g}$ to the working temperature $T$ until $H$ changes. We connect the growth of $\xi(t,t_\mathrm{w};H)$ to the barrier heights $\Delta(t_\mathrm{w})$ that set the dynamics. The effect of $H$ on the magnitude of $\Delta(t_\mathrm{w})$ is responsible for affecting differently the two dynamical protocols associated with turning $H$ off (TRM, or thermoremanent magnetization) or on (ZFC, or zero field-cooled magnetization). In this paper, we display the difference between the zero-field cooled $\xi_{\text {ZFC}}(t,t_\mathrm{w};H)$ and the thermoremanent magnetization $\xi_{\text {TRM}}(t,t_\mathrm{w};H)$ correlation lengths as $H$ increases, both experimentally and through numerical simulations, corresponding to the violation of the extended principle of superposition in line with the finding of the Uppsala Group., Comment: 23 pages and 18 figures

Memory and rejuvenation effects in the magnetic response of off-equilibrium spin glasses have been widely regarded as the doorway into the experimental exploration of ultrametricity and temperature chaos (maybe the most exotic features in glassy free-energy landscapes). Unfortunately, despite more than twenty years of theoretical efforts following the experimental discovery of memory and rejuvenation, these effects have thus far been impossible to simulate reliably. Yet, three recent developments convinced us to accept this challenge: first, the custom-built Janus II supercomputer makes it possible to carry out "numerical experiments" in which the very same quantities that can be measured in single crystals of CuMn are computed from the simulation, allowing for parallel analysis of the simulation/experiment data. Second, Janus II simulations have taught us how numerical and experimental length scales should be compared. Third, we have recently understood how temperature chaos materializes in aging dynamics. All three aspects have proved crucial for reliably reproducing rejuvenation and memory effects on the computer. Our analysis shows that (at least) three different length scales play a key role in aging dynamics, while essentially all theoretical analyses of the aging dynamics emphasize the presence and the crucial role of a single glassy correlation length., Comment: Main text:13 pages, 4 figures SM: 7 pages, 8 figures

Use of dedicated computers in spin glass simulations allows one to equilibrate very large samples (of size as large as $L=32$) and to carry out "computer experiments" that can be compared to (and analyzed in combination with) laboratory experiments on spin-glass samples. In the absence of a magnetic field, the most economic conclusion of the combined analysis of equilibrium and non-equilibrium simulations is that an RSB spin glass phase is present in three spatial dimensions. However, in the presence of a field, the lower critical dimension for the de Almeida-Thouless transition seems to be larger than three., Comment: 28 pages and 11 figures. We have added an additional scenario in Sect. 4. To appear as a contribution to the edited volume "Spin Glass Theory & Far Beyond - Replica Symmetry Breaking after 40 Years", World Scientific

A growing body of evidence indicates that the sluggish low-temperature dynamics of glass formers (e.g. supercooled liquids, colloids or spin glasses) is due to a growing correlation length. Which is the effective field theory that describes these correlations? The natural field theory was drastically simplified by Bray and Roberts in 1980. More than forty years later, we confirm the tenets of Bray and Roberts theory by studying the Ising spin glass in an externally applied magnetic field, both in four spatial dimensions (data obtained from the Janus collaboration) and on the Bethe lattice., Comment: 8 pages, 4 figures, accepted in PRE

A recent 3-XORSAT challenge required to minimize a very complex and rough energy function, typical of glassy models with a random first order transition and a golf course like energy landscape. We present the ideas beyond the quasi-greedy algorithm and its very efficient implementation on GPUs that are allowing us to rank first in such a competition. We suggest a better protocol to compare algorithmic performances and we also provide analytical predictions about the exponential growth of the times to find the solution in terms of free-energy barriers., Comment: 7 pages, 7 figure, EPL format + SM (2 pages)

The synergy between experiment, theory, and simulations enables a microscopic analysis of spin-glass dynamics in a magnetic field in the vicinity of and below the spin-glass transition temperature $T_\mathrm{g}$. The spin-glass correlation length, $\xi(t,t_\mathrm{w};T)$, is analysed both in experiments and in simulations in terms of the waiting time $t_\mathrm{w}$ after the spin glass has been cooled down to a stabilised measuring temperature $T

We find a dynamic effect in the non-equilibrium dynamics of a spin glass that closely parallels equilibrium temperature chaos. This effect, that we name dynamic temperature chaos, is spatially heterogeneous to a large degree. The key controlling quantity is the time-growing spin-glass coherence length. Our detailed characterization of dynamic temperature chaos paves the way for the analysis of recent and forthcoming experiments. This work has been made possible thanks to the most massive simulation to date of non-equilibrium dynamics, carried out on the Janus~II custom-built supercomputer., Comment: Version accepted for publication in Communication Physics 10 pages, 9 figures

The correlation length $\xi$, a key quantity in glassy dynamics, can now be precisely measured for spin glasses both in experiments and in simulations. However, known analysis methods lead to discrepancies either for large external fields or close to the glass temperature. We solve this problem by introducing a scaling law that takes into account both the magnetic field and the time-dependent spin-glass correlation length. The scaling law is successfully tested against experimental measurements in a CuMn single crystal and against large-scale simulations on the Janus II dedicated computer., Comment: Revised version, including supplemental material

The growth of the spin-glass correlation length has been measured as a function of the waiting time $t_{\mathrm{w}}$ on a single crystal of CuMn (6 at.\%), reaching values $\xi\sim 150$ nm, larger than any other glassy correlation-length measured to date. We find an aging rate $\mathrm{d}\ln\,t_{\mathrm{w}}/\mathrm{d}\ln\,\xi$ larger than found in previous measurements, which evinces a dynamic slowing-down as $\xi$ grows. Our measured aging rate is compared with simulation results by the Janus collaboration. After critical effects are taken into account, we find excellent agreement with the Janus data., Comment: Version accepted in PRB (see https://journals.aps.org/prb/abstract/10.1103/PhysRevB.100.094202)

Motivated by recent experiments of exceptional accuracy, we study numerically the spin-glass dynamics in a film geometry. We cover all the relevant time regimes, from picoseconds to equilibrium, at temperatures at and below the 3D critical point. The dimensional crossover from 3D to 2D dynamics, that starts when the correlation length becomes comparable to the film thickness, consists of four dynamical regimes. Our analysis, based on a Renormalization Group transformation, finds consistent the overall physical picture employed by Orbach et al. in the interpretation of their experiments., Comment: 8 pages, 8 figures. Version accepted for publication in Phys. Rev. B

We study the problem of glassy relaxations in the presence of an external field in the highly controlled context of a spin-glass simulation. We consider a small spin glass in three dimensions (specifically, a lattice of size L=8, small enough to be equilibrated through a Parallel Tempering simulations at low temperatures, deep in the spin glass phase). After equilibrating the sample, an external field is switched on, and the subsequent dynamics is studied. The field turns out to reduce the relaxation time, but huge statistical fluctuations are found when different samples are compared. After taking care of these fluctuations we find that the expected linear regime is very narrow. Nevertheless, when regarded as a purely numerical method, we find that the external field is extremely effective in reducing the relaxation times., Comment: 22 pages, 10 figures; Published version

We consider the spatial correlation function of the two-dimensional Ising spin glass under out-equilibrium conditions. We pay special attention to the scaling limit reached upon approaching zero temperature. The field-theory of a non-interacting field makes a surprisingly good job at describing the spatial shape of the correlation function of the out-equilibrium Edwards-Anderson Ising model in two dimensions., Comment: 20 pages + 5 Figures

Recent high precision experimental results on spin-glass films ask for a detailed understanding of the domain-growth dynamics of two-dimensional spin glasses. To achieve this goal, we numerically simulate the out-equilibrium dynamics of the Ising spin glass for a time that spans close to twelve orders of magnitude (from picoseconds to order of a second), in systems large enough to avoid finite-size effects. We find that the time-growth of the size of the glassy domains is excellently described by a single scaling function. A single time-scale $\tau(T)$ controls the dynamics. $\tau(T)$ diverges upon approaching the $T=0$ critical point. The divergence of $\tau(T\to 0)$ is Arrhenius-like, with a barrier height that depends very mildly on temperature. The growth of this barrier-height is best described by critical dynamics. As a side product we obtain an impressive confirmation of universality of the equilibrium behavior of two-dimensional spin-glasses., Comment: 21 pages, 9 figures. Updated references. Added DOI and Journal ref

The Mpemba effect occurs when a hot system cools faster than an initially colder one, when both are refrigerated in the same thermal reservoir. Using the custom built supercomputer Janus II, we study the Mpemba effect in spin glasses and show that it is a non-equilibrium process, governed by the coherence length \xi of the system. The effect occurs when the bath temperature lies in the glassy phase, but it is not necessary for the thermal protocol to cross the critical temperature. In fact, the Mpemba effect follows from a strong relationship between the internal energy and \xi that turns out to be a sure-tell sign of being in the glassy phase. Thus, the Mpemba effect presents itself as an intriguing new avenue for the experimental study of the coherence length in supercooled liquids and other glass formers., Comment: Version accepted for publication in PNAS. 6 pages, 7 figures

Experiments on spin glasses can now make precise measurements of the exponent $z(T)$ governing the growth of glassy domains, while our computational capabilities allow us to make quantitative predictions for experimental scales. However, experimental and numerical values for $z(T)$ have differed. We use new simulations on the Janus II computer to resolve this discrepancy, finding a time-dependent $z(T, t_w)$, which leads to the experimental value through mild extrapolations. Furthermore, theoretical insight is gained by studying a crossover between the $T = T_c$ and $T = 0$ fixed points., Comment: Version accepted for publication in PRL. 12 pages, 9 figures

We have introduced a variational method to improve the computation of integrated correlation times in the Parallel Tempering Dynamics, obtaining a better estimate (a lower bound, at least) of the exponential correlation time. Using this determination of the correlation times, we revisited the problem of the characterization of the chaos in temperature in finite dimensional spin glasses by means of the study of correlations between different chaos indicators computed in the static and the correlation times of the Parallel Tempering dynamics. The sample-distribution of the characteristic time for the Parallel Tempering dynamics turns out to be fat-tailed and it obeys finite-size scaling., Comment: New appendix added. Version to appear in JSTAT

We first reproduce on the Janus and Janus II computers a milestone experiment that measures the spin-glass coherence length through the lowering of free-energy barriers induced by the Zeeman effect. Secondly we determine the scaling behavior that allows a quantitative analysis of a new experiment reported in the companion Letter [S. Guchhait and R. Orbach, Phys. Rev. Lett. 118, 157203 (2017)]. The value of the coherence length estimated through the analysis of microscopic correlation functions turns out to be quantitatively consistent with its measurement through macroscopic response functions. Further, non-linear susceptibilities, recently measured in glass-forming liquids, scale as powers of the same microscopic length., Comment: 6 pages, 4 figures

Chaotic size dependence makes it extremely difficult to take the thermodynamic limit in disordered systems. Instead, the metastate, which is a distribution over thermodynamic states, might have a smooth limit. So far, studies of the metastate have been mostly mathematical. We present a numerical construction of the metastate for the d=3 Ising spin glass. We work in equilibrium, below the critical temperature. Leveraging recent rigorous results, our numerical analysis gives evidence for a "dispersed" metastate, supported on many thermodynamic states., Comment: 5 pages, 5 figures

We report a high-precision numerical estimation of the critical exponent $\alpha$ of the specific heat of the random-field Ising model in four dimensions. Our result $\alpha = 0.12(1)$ indicates a diverging specific-heat behavior and is consistent with the estimation coming from the modified hyperscaling relation using our estimate of $\theta$ via the anomalous dimensions $\eta$ and $\bar{\eta}$. Our analysis benefited form a high-statistics zero-temperature numerical simulation of the model for two distributions of the random fields, namely a Gaussian and Poissonian distribution, as well as recent advances in finite-size scaling and reweighting methods for disordered systems. An original estimate of the critical slowing down exponent $z$ of the maximum-flow algorithm used is also provided., Comment: 11 pages, 3 figures, 1 table. arXiv admin note: text overlap with arXiv:1512.06571

The unifying feature of glass formers (such as polymers, supercooled liquids, colloids, granulars, spin glasses, superconductors, ...) is a sluggish dynamics at low temperatures. Indeed, their dynamics is so slow that thermal equilibrium is never reached in macroscopic samples: in analogy with living beings, glasses are said to age. Here, we show how to relate experimentally relevant quantities with the experimentally unreachable low-temperature equilibrium phase. We have performed a very accurate computation of the non-equilibrium fluctuation-dissipation ratio for the three-dimensional Edwards-Anderson Ising spin glass, by means of large-scale simulations on the special-purpose computers Janus and Janus II. This ratio (computed for finite times on very large, effectively infinite, systems) is compared with the equilibrium probability distribution of the spin overlap for finite sizes. The resulting quantitative statics-dynamics dictionary, based on observables that can be measured with current experimental methods, could allow the experimental exploration of important features of the spin-glass phase without uncontrollable extrapolations to infinite times or system sizes., Comment: Version accepted for publication in PNAS. Reuploaded to fix a typo in the list of authors

Temperature chaos plays a role in important effects, like for example memory and rejuvenation, in spin glasses, colloids, polymers. We numerically investigate temperature chaos in spin glasses, exploiting its recent characterization as a rare-event driven phenomenon. The peculiarities of the transformation from periodic to anti-periodic boundary conditions in spin glasses allow us to conclude that temperature chaos is non-local: no bounded region of the system causes it. We precise the statistical relationship between temperature chaos and the free-energy changes upon varying boundary conditions., Comment: 15 pages, 8 figures. Version accepted for publication in JSTAT

We use finite size scaling to study Ising spin glasses in two spatial dimensions. The issue of universality is addressed by comparing discrete and continuous probability distributions for the quenched random couplings. The sophisticated temperature dependency of the scaling fields is identified as the major obstacle that has impeded a complete analysis. Once temperature is relinquished in favor of the correlation length as the basic variable, we obtain a reliable estimation of the anomalous dimension and of the thermal critical exponent. Universality among binary and Gaussian couplings is confirmed to a high numerical accuracy., Comment: 11 pages,7 figures

We use a sample-dependent analysis, based on medians and quantiles, to analyze the behavior of the overlap probability distribution of the Sherrington-Kirkpatrick and 3D Edwards-Anderson models of Ising spin glasses. We find that this approach is an effective tool to distinguish between RSB-like and droplet-like behavior of the spin-glass phase. Our results are in agreement with a RSB-like behavior for the 3D Edwards-Anderson model., Comment: Version accepted in PRB. 12 pages, 10 figures

We perform equilibrium parallel-tempering simulations of the 3D Ising Edwards-Anderson spin glass in a field. A traditional analysis shows no signs of a phase transition. Yet, we encounter dramatic fluctuations in the behaviour of the model: Averages over all the data only describe the behaviour of a small fraction of it. Therefore we develop a new approach to study the equilibrium behaviour of the system, by classifying the measurements as a function of a conditioning variate. We propose a finite-size scaling analysis based on the probability distribution function of the conditioning variate, which may accelerate the convergence to the thermodynamic limit. In this way, we find a non-trivial spectrum of behaviours, where a part of the measurements behaves as the average, while the majority of them shows signs of scale invariance. As a result, we can estimate the temperature interval where the phase transition in a field ought to lie, if it exists. Although this would-be critical regime is unreachable with present resources, the numerical challenge is finally well posed., Comment: 42 pages, 19 figures. Minor changes and added figure (results unchanged)

We report a high-precision finite-size scaling study of the critical behavior of the three-dimensional Ising Edwards-Anderson model (the Ising spin glass). We have thermalized lattices up to L=40 using the Janus dedicated computer. Our analysis takes into account leading-order corrections to scaling. We obtain Tc = 1.1019(29) for the critical temperature, \nu = 2.562(42) for the thermal exponent, \eta = -0.3900(36) for the anomalous dimension and \omega = 1.12(10) for the exponent of the leading corrections to scaling. Standard (hyper)scaling relations yield \alpha = -5.69(13), \beta = 0.782(10) and \gamma = 6.13(11). We also compute several universal quantities at Tc., Comment: 9 pages, 5 figures

This paper describes the architecture, the development and the implementation of Janus II, a new generation application-driven number cruncher optimized for Monte Carlo simulations of spin systems (mainly spin glasses). This domain of computational physics is a recognized grand challenge of high-performance computing: the resources necessary to study in detail theoretical models that can make contact with experimental data are by far beyond those available using commodity computer systems. On the other hand, several specific features of the associated algorithms suggest that unconventional computer architectures, which can be implemented with available electronics technologies, may lead to order of magnitude increases in performance, reducing to acceptable values on human scales the time needed to carry out simulation campaigns that would take centuries on commercially available machines. Janus II is one such machine, recently developed and commissioned, that builds upon and improves on the successful JANUS machine, which has been used for physics since 2008 and is still in operation today. This paper describes in detail the motivations behind the project, the computational requirements, the architecture and the implementation of this new machine and compares its expected performances with those of currently available commercial systems., Comment: 28 pages, 6 figures

We characterize the phase diagram of anisotropic Heisenberg spin glasses, finding both the spin and the chiral glass transition. We remark the presence of strong finite-size effects on the chiral sector. We find a unique phase transition for the chiral and spin glass sector, in the Universality class of Ising spin glasses. We focus on keeping finite-size effects under control, and we stress that they are important to understand experiments. Thanks to large GPU clusters we have been able to thermalize cubic lattices with up to 64x64x64 spins, over a vast range of temperatures., Comment: 12 pages, 6 figures. New simulations added for L=64, D=1 (results unchanged)

We study the off-equilibrium dynamics of the three-dimensional Ising spin glass in the presence of an external magnetic field. We have performed simulations both at fixed temperature and with an annealing protocol. Thanks to the Janus special-purpose computer, based on FPGAs, we have been able to reach times equivalent to 0.01 seconds in experiments. We have studied the system relaxation both for high and for low temperatures, clearly identifying a dynamical transition point. This dynamical temperature is strictly positive and depends on the external applied magnetic field. We discuss different possibilities for the underlying physics, which include a thermodynamical spin-glass transition, a mode-coupling crossover or an interpretation reminiscent of the random first-order picture of structural glasses., Comment: Version accepted in PRE. 23 pages, 23 figures

Temperature chaos has often been reported in literature as a rare-event driven phenomenon. However, this fact has always been ignored in the data analysis, thus erasing the signal of the chaotic behavior (still rare in the sizes achieved) and leading to an overall picture of a weak and gradual phenomenon. On the contrary, our analysis relies on a large-deviations functional that allows to discuss the size dependencies. In addition, we had at our disposal unprecedentedly large configurations equilibrated at low temperatures, thanks to the Janus computer. According to our results, when temperature chaos occurs its effects are strong and can be felt even at short distances., Comment: 5 pages, 5 figures

A recent interesting paper [Yucesoy et al. Phys. Rev. Lett. 109, 177204 (2012), arXiv:1206:0783] compares the low-temperature phase of the 3D Edwards-Anderson (EA) model to its mean-field counterpart, the Sherrington-Kirkpatrick (SK) model. The authors study the overlap distributions P_J(q) and conclude that the two models behave differently. Here we notice that a similar analysis using state-of-the-art, larger data sets for the EA model (generated with the Janus computer) leads to a very clear interpretation of the results of Yucesoy et al., showing that the EA model behaves as predicted by the replica symmetry breaking (RSB) theory., Comment: Version accepted for publication in PRL. 1 page, 1 figure

We compare the critical behavior of the short-range Ising spin glass with a spin glass with long-range interactions which fall off as a power sigma of the distance. We show that there is a value of sigma of the long-range model for which the critical behavior is very similar to that of the short-range model in four dimensions. We also study a value of sigma for which we find the critical behavior to be compatible with that of the three dimensional model, though we have much less precision than in the four-dimensional case., Comment: 15 pages, 10 figures

The low temperature excitations in the anisotropic antiferromagnetic Fe_{1-x} Zn_x F_2 for x=0.25 and 0.31, at and just above the magnetic percolation threshold concentration x_p=0.25, were measured using inelastic neutron scattering. The excitations were simulated for x=0.31 using a localized, classical excitation model, which accounts well for the energies and relative intensities of the excitations observed in the scattering experiments., Comment: 6 pages, 6 figures

We present a microcanonical Monte Carlo simulation of the site-diluted Potts model in three dimensions with eight internal states, partly carried out in the citizen supercomputer Ibercivis. Upon dilution, the pure model's first-order transition becomes of the second-order at a tricritical point. We compute accurately the critical exponents at the tricritical point. As expected from the Cardy-Jacobsen conjecture, they are compatible with their Random Field Ising Model counterpart. The conclusion is further reinforced by comparison with older data for the Potts model with four states., Comment: Final version. 9 pages, 9 figures

We describe Janus, a massively parallel FPGA-based computer optimized for the simulation of spin glasses, theoretical models for the behavior of glassy materials. FPGAs (as compared to GPUs or many-core processors) provide a complementary approach to massively parallel computing. In particular, our model problem is formulated in terms of binary variables, and floating-point operations can be (almost) completely avoided. The FPGA architecture allows us to run many independent threads with almost no latencies in memory access, thus updating up to 1024 spins per cycle. We describe Janus in detail and we summarize the physics results obtained in four years of operation of this machine; we discuss two types of physics applications: long simulations on very large systems (which try to mimic and provide understanding about the experimental non-equilibrium dynamics), and low-temperature equilibrium simulations using an artificial parallel tempering dynamics. The time scale of our non-equilibrium simulations spans eleven orders of magnitude (from picoseconds to a tenth of a second). On the other hand, our equilibrium simulations are unprecedented both because of the low temperatures reached and for the large systems that we have brought to equilibrium. A finite-time scaling ansatz emerges from the detailed comparison of the two sets of simulations. Janus has made it possible to perform spin-glass simulations that would take several decades on more conventional architectures. The paper ends with an assessment of the potential of possible future versions of the Janus architecture, based on state-of-the-art technology., Comment: 19 pages, 3 figures

Spin glasses are a longstanding model for the sluggish dynamics that appears at the glass transition. However, spin glasses differ from structural glasses for a crucial feature: they enjoy a time reversal symmetry. This symmetry can be broken by applying an external magnetic field, but embarrassingly little is known about the critical behaviour of a spin glass in a field. In this context, the space dimension is crucial. Simulations are easier to interpret in a large number of dimensions, but one must work below the upper critical dimension (i.e., in d<6) in order for results to have relevance for experiments. Here we show conclusive evidence for the presence of a phase transition in a four-dimensional spin glass in a field. Two ingredients were crucial for this achievement: massive numerical simulations were carried out on the Janus special-purpose computer, and a new and powerful finite-size scaling method., Comment: 10 pages, 6 figures

Using the results of large scale numerical simulations we study the probability distribution of the pseudo critical temperature for the three-dimensional Edwards-Anderson Ising spin glass and for the fully connected Sherrington-Kirkpatrick model. We find that the behavior of our data is nicely described by straightforward finite-size scaling relations., Comment: 23 pages, 9 figures. Version accepted for publication in J. Stat. Mech

We study the sample-to-sample fluctuations of the overlap probability densities from large-scale equilibrium simulations of the three-dimensional Edwards-Anderson spin glass below the critical temperature. Ultrametricity, Stochastic Stability and Overlap Equivalence impose constraints on the moments of the overlap probability densities that can be tested against numerical data. We found small deviations from the Ghirlanda-Guerra predictions, which get smaller as system size increases. We also focus on the shape of the overlap distribution, comparing the numerical data to a mean-field-like prediction in which finite-size effects are taken into account by substituting delta functions with broad peaks, Comment: 11 pages, 8 figures; version accepted for publication on PRB and updated Refs

Tethering methods allow us to perform Monte Carlo simulations in ensembles with conserved quantities. Specifically, one couples a reservoir to the physical magnitude of interest, and studies the statistical ensemble where the total magnitude (system+reservoir) is conserved. The reservoir is actually integrated out, which leaves us with a fluctuation-dissipation formalism that allows us to recover the appropriate Helmholtz effective potential with great accuracy. These methods are demonstrating a remarkable flexibility. In fact, we illustrate two very different applications: hard spheres crystallization and the phase transition of the diluted antiferromagnet in a field (the physical realization of the random field Ising model). The tethered approach holds the promise to transform cartoon drawings of corrugated free-energy landscapes into real computations. Besides, it reduces the algorithmic dynamic slowing-down, probably because the conservation law holds non-locally., Comment: 44 pages

We study the critical behavior of the Diluted Antiferromagnet in a Field with the Tethered Monte Carlo formalism. We compute the critical exponents (including the elusive hyperscaling violations exponent $\theta$). Our results provide a comprehensive description of the phase transition and clarify the inconsistencies between previous experimental and theoretical work. To do so, our method addresses the usual problems of numerical work (large tunneling barriers and self-averaging violations)., Comment: 4 pages, 2 figures

We present a tethered Monte Carlo simulation of the crystallization of hard spheres. Our method boosts the traditional umbrella sampling to the point of making practical the study of constrained Gibb's free energies depending on several crystalline order-parameters. We obtain high-accuracy estimates of the fluid-crystal coexistence pressure for up to 2916 particles (enough to accommodate fluid-solid interfaces). We are able to extrapolate to infinite volume the coexistence pressure (p_{co}=11.5727(10) k_B T/\sigma^3) and the interfacial free energy (\gamma_{100}=0.636(11) k_B T/\sigma^2)., Comment: 6 pages, 4 pdf figures. Version to be published in PRL. Appendices contain Supplemental Material

We perform a finite size scaling study of the three-dimensional Heisenberg spin glass in the presence of weak random anisotropic interactions, up to sizes L=32. Anisotropies have a major impact on the phase transition. The chiral-glass susceptibility does not diverge, due to a large anomalous dimension. It follows that the anisotropic spin-glass belongs to a Universality Class different from the isotropic model, which questions the applicability of the chirality scenario., Comment: 5 pages, 4 postscript figures

Diagrammatic techniques to compute perturbatively the spectral properties of Euclidean Random Matrices in the high-density regime are introduced and discussed in detail. Such techniques are developed in two alternative and very different formulations of the mathematical problem and are shown to give identical results up to second order in the perturbative expansion. One method, based on writing the so-called resolvent function as a Taylor series, allows to group the diagrams in a small number of topological classes, providing a simple way to determine the infrared (small momenta) behavior of the theory up to third order, which is of interest for the comparison with experiments. The other method, which reformulates the problem as a field theory, can instead be used to study the infrared behaviour at any perturbative order., Comment: 29 pages

We numerically study the aging properties of the dynamical heterogeneities in the Ising spin glass. We find that a phase transition takes place during the aging process. Statics-dynamics correspondence implies that systems of finite size in equilibrium have static heterogeneities that obey Finite-Size Scaling, thus signaling an analogous phase transition in the thermodynamical limit. We compute the critical exponents and the transition point in the equilibrium setting, and use them to show that aging in dynamic heterogeneities can be described by a Finite-Time Scaling Ansatz, with potential implications for experimental work., Comment: Version accepted for publication in Phys. Rev. Lett. 4 pages, 3 postscript figures

We present a massive equilibrium simulation of the three-dimensional Ising spin glass at low temperatures. The Janus special-purpose computer has allowed us to equilibrate, using parallel tempering, L=32 lattices down to T=0.64 Tc. We demonstrate the relevance of equilibrium finite-size simulations to understand experimental non-equilibrium spin glasses in the thermodynamical limit by establishing a time-length dictionary. We conclude that non-equilibrium experiments performed on a time scale of one hour can be matched with equilibrium results on L=110 lattices. A detailed investigation of the probability distribution functions of the spin and link overlap, as well as of their correlation functions, shows that Replica Symmetry Breaking is the appropriate theoretical framework for the physically relevant length scales. Besides, we improve over existing methodologies to ensure equilibration in parallel tempering simulations., Comment: 48 pages, 19 postscript figures, 9 tables. Version accepted for publication in the Journal of Statistical Mechanics