Scaling Law Describes the Spin-Glass Response in Theory, Experiments, and Simulations.

The correlation length ξ, 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.

Aging Rate of Spin Glasses from Simulations Matches Experiments.

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.

Numerical Construction of the Aizenman-Wehr Metastate.

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.

Matching Microscopic and Macroscopic Responses in Glasses.

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)].PRLTAO0031-900710.1103/PhysRevLett.118.157203 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, nonlinear susceptibilities, recently measured in glass-forming liquids, scale as powers of the same microscopic length.

Erratum: Large-scale kinetic roughening behavior of coffee-ring fronts [Phys. Rev. E 106, 044801 (2022)].

This corrects the article DOI: 10.1103/PhysRevE.106.044801.

Evidence of a second-order phase transition in the six-dimensional Ising spin glass in a field.

The very existence of a phase transition for spin glasses in an external magnetic field is controversial, even in high dimensions. We carry out massive simulations of the Ising spin-glass in a field, in six dimensions (which, according to classical-but not generally accepted-field-theoretical studies, is the upper critical dimension). We obtain results compatible with a second-order phase transition and estimate its critical exponents for the simulated lattice sizes. The detailed analysis performed by other authors of the replica symmetric Hamiltonian, under the hypothesis of critical behavior, predicts that the ratio of the renormalized coupling constants remain bounded as the correlation length grows. Our numerical results are in agreement with this expectation.

Universal interface fluctuations in the contact process.

We study the interface representation of the contact process at its directed-percolation critical point, where the scaling properties of the interface can be related to those of the original particle model. Interestingly, such a behavior happens to be intrinsically anomalous and more complex than that described by the standard Family-Vicsek dynamic scaling Ansatz of surface kinetic roughening. We expand on a previous numerical study by Dickman and Muñoz [Phys. Rev. E 62, 7632 (2000)10.1103/PhysRevE.62.7632] to fully characterize the kinetic roughening universality class for interface dimensions d=1,2, and 3. Beyond obtaining scaling exponent values, we characterize the interface fluctuations via their probability density function (PDF) and covariance, seen to display universal properties which are qualitatively similar to those recently assessed for the Kardar-Parisi-Zhang (KPZ) and other important universality classes of kinetic roughening. Quantitatively, while for d=1 the interface covariance seems to be well described by the KPZ, Airy_{1} covariance, no such agreement occurs in terms of the fluctuation PDF or the scaling exponents.

Erratum: Numerical test of the replica-symmetric Hamiltonian for correlations of the critical state of spin glasses in a field [Phys. Rev. E 105, 054106 (2022)].

This corrects the article DOI: 10.1103/PhysRevE.105.054106.

Critical behavior of the three-dimensional random-anisotropy Heisenberg model.

We have studied the critical properties of the three-dimensional random anisotropy Heisenberg model by means of numerical simulations using the Parallel Tempering method. We have simulated the model with two different disorder distributions, cubic and isotropic ones, with two different anisotropy strengths for each disorder class. For the case of the anisotropic disorder, we have found evidence of universality by finding critical exponents and universal dimensionless ratios independent of the strength of the disorder. In the case of isotropic disorder distribution the situation is very involved: we have found two phase transitions in the magnetization channel which are merging for larger lattices remaining a zero magnetization low-temperature phase. Studying this region using a spin-glass order parameter we have found evidence for a spin-glass phase transition. We have estimated effective critical exponents for the spin-glass phase transition for the different values of the strength of the isotropic disorder, discussing the crossover regime.

Spreading fronts of wetting liquid droplets: Microscopic simulations and universal fluctuations.

We have used kinetic Monte Carlo (kMC) simulations of a lattice gas to study front fluctuations in the spreading of a nonvolatile liquid droplet onto a solid substrate. Our results are consistent with a diffusive growth law for the radius of the precursor layer, R∼t^{δ}, with δ≈1/2 in all the conditions considered for temperature and substrate wettability, in good agreement with previous studies. The fluctuations of the front exhibit kinetic roughening properties with exponent values which depend on temperature T, but become T independent for sufficiently high T. Moreover, strong evidence of intrinsic anomalous scaling has been found, characterized by different values of the roughness exponent at short and large length scales. Although such a behavior differs from the scaling properties of the one-dimensional Kardar-Parisi-Zhang (KPZ) universality class, the front covariance and the probability distribution function of front fluctuations found in our kMC simulations do display KPZ behavior, agreeing with simulations of a continuum height equation proposed in this context. However, this equation does not feature intrinsic anomalous scaling, at variance with the discrete model.

Large-scale kinetic roughening behavior of coffee-ring fronts.

We have studied the kinetic roughening behavior of the fronts of coffee-ring aggregates via extensive numerical simulations of the off-lattice model considered for this context [Dias et al., Soft Matter 14, 1903 (2018)1744-683X10.1039/C7SM02136D]. This model describes ballistic aggregation of patchy colloids and depends on a parameter r_{AB} which controls the affinity of the two patches, A and B. Suitable boundary conditions allow us to elucidate a discontinuous pinning-depinning transition at r_{AB}=0, with the front displaying intrinsic anomalous scaling, but with unusual exponent values α≃1.2, α_{loc}≃0.5, ß≃1, and z≃1.2. For 00.01 and the system suffers a strong crossover dominated by the r_{AB}=0 behavior for r_{AB}≤0.01. A detailed analysis of correlation functions shows that the aggregate fronts are always in the moving phase for 0

Numerical test of the replica-symmetric Hamiltonian for correlations of the critical state of spin glasses in a field.

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 40 years later, we confirm the tenets of Bray and Roberts's 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.

Static versus dynamic heterogeneities in the D = 3 Edwards-Anderson-Ising spin glass.

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.

Ising spin-glass transition in a magnetic field outside the limit of validity of mean-field theory.

The spin-glass transition in a magnetic field is studied both in and out of the limit of validity of mean-field theory on a diluted one dimensional chain of Ising spins where exchange bonds occur with a probability decaying as the inverse power of the distance. Varying the power in this long-range model corresponds, in a one-to-one relationship, to changing the dimension in spin-glass short-range models. Evidence for a spin-glass transition in a magnetic field is found also for systems whose equivalent dimension is below the upper critical dimension in a zero magnetic field.

Mean field model of coagulation and annihilation reactions in a medium of quenched traps: subdiffusion.

We present a mean field model for coagulation (A+A-->A) and annihilation (A+A-->0) reactions on lattices of traps with a distribution of depths reflected in a distribution of mean escape times. The escape time from each trap is exponentially distributed about the mean for that trap, and the distribution of mean escape times is a power law. Even in the absence of reactions, the distribution of particles over sites changes with time as particles are caught in ever deeper traps, that is, the distribution exhibits aging. Our main goal is to explore whether the reactions lead to further (time dependent) changes in this distribution.

Computation of the dynamic critical exponent of the three-dimensional Heisenberg model.

Working in and out of equilibrium and using state-of-the-art techniques we have computed the dynamic critical exponent of the three-dimensional Heisenberg model. By computing the integrated autocorrelation time at equilibrium, for lattice sizes L≤64, we have obtained z=2.033(5). In the out-of-equilibrium regime we have run very large lattices (L≤250) obtaining z=2.04(2) from the growth of the correlation length. We compare our values with that previously computed at equilibrium with relatively small lattices (L≤24), with that provided by means a three-loops calculation using perturbation theory and with experiments. Finally we have checked previous estimates of the static critical exponents, η and ν, in the out-of-equilibrium regime.

Scaling analysis of the site-diluted Ising model in two dimensions.

A combination of recent numerical and theoretical advances are applied to analyze the scaling behavior of the site-diluted Ising model in two dimensions, paying special attention to the implications for multiplicative logarithmic corrections. The analysis focuses primarily on the odd sector of the model (i.e., that associated with magnetic exponents), and in particular on its Lee-Yang zeros, which are determined to high accuracy. Scaling relations are used to connect to the even (thermal) sector, and a first analysis of the density of zeros yields information on the specific heat and its corrections. The analysis is fully supportive of the strong scaling hypothesis and of the scaling relations for logarithmic corrections.

Target problem with evanescent subdiffusive traps.

We calculate the survival probability of a stationary target in one dimension surrounded by diffusive or subdiffusive traps of time-dependent density. The survival probability of a target in the presence of traps of constant density is known to go to zero as a stretched exponential whose specific power is determined by the exponent that characterizes the motion of the traps. A density of traps that grows in time always leads to an asymptotically vanishing survival probability. Trap evanescence leads to a survival probability of the target that may go to zero or to a finite value indicating a probability of eternal survival, depending on the way in which the traps disappear with time.

Universal behavior of crystalline membranes: Crumpling transition and Poisson ratio of the flat phase.

We revisit the universal behavior of crystalline membranes at and below the crumpling transition, which pertains to the mechanical properties of important soft and hard matter materials, such as the cytoskeleton of red blood cells or graphene. Specifically, we perform large-scale Monte Carlo simulations of a triangulated two-dimensional phantom network which is freely fluctuating in three-dimensional space. We obtain a continuous crumpling transition characterized by critical exponents which we estimate accurately through the use of finite-size techniques. By controlling the scaling corrections, we additionally compute with high accuracy the asymptotic value of the Poisson ratio in the flat phase, thus characterizing the auxetic properties of this class of systems. We obtain agreement with the value which is universally expected for polymerized membranes with a fixed connectivity.