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.

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.

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.

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 Collapse

Key Words Collapse

MESH Headings Collapse

Grants Collapse

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.

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.

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.

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.

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.

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.

Dynamical transition in the D=3 Edwards-Anderson spin glass in an external magnetic field

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 field-programmable gate array (FPGAs), we have been able to reach times equivalent to 0.01 s 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.

Scaling behavior of the Heisenberg model in three dimensions

We report on extensive numerical simulations of the three-dimensional Heisenberg model and its analysis through finite-size scaling of Lee-Yang zeros. Besides the critical regime, we also investigate scaling in the ferromagnetic phase. We show that, in this case of broken symmetry, the corrections to scaling contain information on the Goldstone modes. We present a comprehensive Lee-Yang analysis, including the density of zeros, and confirm recent numerical estimates for critical exponents.

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.

Microcanonical finite-size scaling in second-order phase transitions with diverging specific heat

A microcanonical finite-size ansatz in terms of quantities measurable in a finite lattice allows extending phenomenological renormalization (the so-called quotients method) to the microcanonical ensemble. The ansatz is tested numerically in two models where the canonical specific heat diverges at criticality, thus implying Fisher renormalization of the critical exponents: the three-dimensional ferromagnetic Ising model and the two-dimensional four-state Potts model (where large logarithmic corrections are known to occur in the canonical ensemble). A recently proposed microcanonical cluster method allows simulating systems as large as L=1024 (Potts) or L=128 (Ising). The quotients method provides accurate determinations of the anomalous dimension, eta, and of the (Fisher-renormalized) thermal nu exponent. While in the Ising model the numerical agreement with our theoretical expectations is very good, in the Potts case, we need to carefully incorporate logarithmic corrections to the microcanonical ansatz in order to rationalize our data.

Coagulation reactions in low dimensions: Revisiting subdiffusive A+A reactions in one dimension

We present a theory for the coagulation reaction A+A-->A for particles moving subdiffusively in one dimension. Our theory is tested against numerical simulations of the concentration of A particles as a function of time ("anomalous kinetics") and of the interparticle distribution function as a function of interparticle distance and time. We find that the theory captures the correct behavior asymptotically and also at early times, and that it does so whether the particles are nearly diffusive or very subdiffusive. We find that, as in the normal diffusion problem, an interparticle gap responsible for the anomalous kinetics develops and grows with time. This corrects an earlier claim to the contrary on our part.

Site-diluted Ising model in four dimensions

In the literature, there are five distinct fragmented sets of analytic predictions for the scaling behavior at the phase transition in the random-site Ising model in four dimensions. Here, the scaling relations for logarithmic corrections are used to complete the scaling pictures for each set. A numerical approach is then used to confirm the leading scaling picture coming from these predictions and to discriminate between them at the level of logarithmic corrections.

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.

Nonequilibrium spin-glass dynamics from picoseconds to a tenth of a second

We study numerically the nonequilibrium dynamics of the Ising spin glass, for a time spanning 11 orders of magnitude, thus approaching the experimentally relevant scale (i.e., seconds). We introduce novel analysis techniques to compute the coherence length in a model-independent way. We present strong evidence for a replicon correlator and for overlap equivalence. The emerging picture is compatible with noncoarsening behavior.

Dilute one-dimensional spin glasses with power law decaying interactions

We introduce a diluted version of the one-dimensional spin-glass model with interactions decaying in probability as an inverse power of the distance. In this model, varying the power corresponds to changing the dimension in short-range models. The spin-glass phase is studied in and out of the range of validity of the mean-field approximation in order to discriminate between different theories. Since each variable interacts only with a finite number of others the cost for simulating, the model is drastically reduced with respect to the fully connected version, and larger sizes can be studied. We find both static and dynamic indications in favor of the so-called replica symmetry breaking theory.

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.

First-order transition in a three-dimensional disordered system

We present the first detailed numerical study in three dimensions of a first-order phase transition that remains first order in the presence of quenched disorder (specifically, the ferromagnetic-paramagnetic transition of the site-diluted four states Potts model). A tricritical point, which lies surprisingly near the pure-system limit and is studied by means of finite-size scaling, separates the first-order and second-order parts of the critical line. This investigation has been made possible by a new definition of the disorder average that avoids the diverging-variance probability distributions that plague the standard approach. Entropy, rather than free energy, is the basic object in this approach that exploits a recently introduced microcanonical Monte Carlo method.

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.

Phase transition in tensionless surfaces

We study the critical behavior of the Laplacian roughening model, which describes the growth of tensionless surfaces. This type of growth phenomena is very important, for instance, in biological membranes and in molecular beam epitaxy. We summarize previous analytical and numerical results and point out their contradictions and differences, thus making clear the context of our work. Our contribution, achieved through large-scale numerical simulations, is the confirmation that the model exhibits a single continuous phase transition: the transition takes place between a continuum massless (i.e., with infinite correlation length) bilaplacian behavior and a massive propagator (finite correlation length).

Universality in the off-equilibrium critical dynamics of the three-dimensional diluted Ising model

We study the off-equilibrium critical dynamics of the three-dimensional diluted Ising model. We compute the dynamical critical exponent z and we show that it is independent of the dilution only when we take into account the scaling corrections to the dynamics. Finally, we will compare our results with the experimental data.