Universal Scaling Function Research Articles

Dynamic fluctuations of current and mass in nonequilibrium mass transport processes

We study steady-state dynamic fluctuations of current and mass in several variants of random average processes on a ring of L sites. These processes violate detailed balance in the bulk and have nontrivial spatial structures: their steady states are not described by the Boltzmann–Gibbs distribution and can have nonzero spatial correlations. Using a microscopic approach, we exactly calculate the second cumulants, or the variance, ⟨Qi2(T)⟩c and ⟨Qsub2(l,T)⟩c , of cumulative (time-integrated) currents up to time T across the ith bond and across a subsystem of size l (summed over bonds in the subsystem), respectively. We also calculate the (two-point) dynamic correlation function of the subsystem mass. In particular, we show that, for large L≫1 , the second cumulant ⟨Qi2(T)⟩c of the cumulative current up to time T across the ith bond grows linearly as ⟨Qi2⟩c∼T for initial times T∼O(1) , subdiffusively as ⟨Qi2⟩c∼T1/2 for intermediate times 1≪T≪L2 , and then again linearly as ⟨Qi2⟩c∼T for long times T≫L2 . The scaled cumulant liml→∞,T→∞⟨Qsub2(l,T)⟩c/2lT of current across the subsystem of size l and up to time T converges to the density-dependent particle mobility χ(ρ) when the large subsystem-size limit is taken first, followed by the large-time limit; when the limits are reversed, it simply vanishes. Remarkably, regardless of the dynamical rules, the scaled current cumulant D⟨Qi2(T)⟩c/2χL≡W(y) as a function of scaled time y=DT/L2 can be expressed in terms of a universal scaling function W(y) , where D is the bulk-diffusion coefficient; interestingly, the intermediate-time subdiffusive and long-time diffusive growths can be connected through a single scaling function W(y) . The power spectra for current and mass are also exactly characterized by the respective scaling functions. Furthermore, we provide a microscopic derivation of equilibrium-like Green–Kubo and Einstein relations that connect the steady-state current fluctuations to an ‘operational’ mobility (i.e. the response to an external force field) and mass fluctuation, respectively.

Corrections to Diffusion in Interacting Quantum Systems

The approach to equilibrium in interacting classical and quantum systems is a challenging problem of both theoretical and experimental interest. One useful organizing principle characterizing equilibration is the dissipative universality class, the most prevalent one being diffusion. In this paper, we use the effective field theory (EFT) of diffusion to systematically obtain universal power-law corrections to diffusion. We then employ large-scale simulations of classical and quantum systems to explore their validity. In particular, we find universal scaling functions for the corrections to the dynamical structure factor ⟨n(x,t)n⟩, in the presence of a single U(1) or SU(2) charge in systems with and without particle-hole symmetry, and present the framework to generalize the calculation to multiple charges. Classical simulations show remarkable agreement with EFT predictions for subleading corrections, pushing precision tests of effective theories for thermalizing systems to an unprecedented level. Moving to quantum systems, we perform large-scale tensor-network simulations in unitary and noisy 1D Floquet systems with conserved magnetization. We find a qualitative agreement with EFT, which becomes quantitative in the case of noisy systems. Additionally, we show how the knowledge of EFT corrections allows for fitting methods, which can improve the estimation of transport parameters at the intermediate times accessible by simulations and experiments. Finally, we explore nonlinear response in quantum systems and find that EFT provides an accurate prediction for its behavior. Our results provide a basis for a better understanding of the nonlinear phenomena present in thermalizing systems. Published by the American Physical Society 2024

Effect of constraint relaxation on dynamic critical phenomena in minimum vertex cover problem

The effects of constraint relaxation on dynamic critical phenomena in the Minimum Vertex Cover (MVC) problem on Erdős-Rényi random graphs are investigated using Markov chain Monte Carlo simulations. Following our previous work that revealed the reduction of the critical temperature by constraint relaxation based on the penalty function method, this study focuses on investigating the critical properties of the relaxation time along its phase boundary. It is found that the dynamical correlation function of MVC with respect to the problem size and the constraint strength follows a universal scaling function. The analysis shows that the relaxation time decreases as the constraints are relaxed. This decrease is more pronounced for the critical amplitude than for the critical exponent, and this result is interpreted in terms of the system's microscopic energy barriers due to the constraint relaxation.

Observation of Subdiffusive Dynamic Scaling in a Driven and Disordered Bose Gas.

We explore the dynamics of a tuneable box-trapped Bose gas under strong periodic forcing in the presence of weak disorder. In absence of interparticle interactions, the interplay of the drive and disorder results in an isotropic nonthermal momentum distribution that shows subdiffusive dynamic scaling, with sublinear energy growth and the universal scaling function captured well by a compressed exponential. We explain that this subdiffusion in momentum space can naturally be understood as a random walk in energy space. We also experimentally show that for increasing interaction strength, the gas behavior smoothly crosses over to wave turbulence characterized by a power-law momentum distribution, which opens new possibilities for systematic studies of the interplay of disorder and interactions in driven quantum systems.

Nonlinear response in diffusive systems

Nonintegrable systems thermalize, leading to the emergence of fluctuating hydrodynamics. Typically, this hydrodynamics is diffusive. We use the effective field theory (EFT) of diffusion to compute higher-point functions of conserved densities. We uncover a simple scaling behavior of correlators at late times, and, focusing on three and four-point functions, derive the asymptotically exact universal scaling functions that characterize nonlinear response in diffusive systems. This allows for precision tests of thermalization beyond linear response in quantum and classical many-body systems. We confirm our predictions in a classical lattice gas.

Evidence of second-order transition and critical scaling for the dynamical ordering transition in current-driven vortices

Dynamical ordering from a disordered plastic flow to an anisotropically ordered smectic flow induced by a dc force has been studied in various many-particle systems, including vortices in type-II superconductors. However, it remains unclear whether the dynamical ordering is a true phase transition because of lack of suitable experimental methods. Here, we study the response of vortex flow to the transverse force using a cross-shaped amorphous Mox\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$_{x}$$\\end{document}Ge1-x\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$_{1-x}$$\\end{document} film. From transverse current-voltage (force-velocity) characteristics under various longitudinal currents, we find a change of the transverse response in low voltage (velocity) regions from a nonlinear to linear behavior at a well-defined longitudinal current that marks the dynamical ordering transition. We also find the scaling collapse of the transverse current-voltage curves to a universal scaling function, providing evidence of the second-order transition for the dynamical ordering transition.

Pre-equilibrium Photon and Dilepton Production

We use QCD kinetic theory to compute photon and dilepton production in the chemically equilibrating out-of-equilibrium quark-gluon plasma created in the early stages of high-energy heavy-ion collisions. We derive universal scaling functions for the pre-equilibrium spectra of photons and dileptons. These scaling functions can be used to make realistic predictions for the pre-equilibrium emission and consequently establish the significance of the preequilibrium phase for the production of electromagnetic probes in heavy-ion collisions.

Phase diagram near the quantum critical point in Schwinger model at θ = π: analogy with quantum Ising chain

Abstract The Schwinger model, 1D quantum electrodynamics, has CP symmetry at θ = π due to the topological nature of the θ term. At zero temperature, it is known that as the fermion mass increases, the system undergoes a second-order phase transition to the CP broken phase, which belongs to the same universality class as the quantum Ising chain. In this paper, we obtain the phase diagram near the quantum critical point (QCP) in the temperature and fermion mass plane using first-principle Monte Carlo simulations, while avoiding the sign problem by using the lattice formulation of the bosonized Schwinger model. Specifically, we perform a detailed investigation of the correlation function of the electric field near the QCP and find that its asymptotic behavior can be described by the universal scaling function of the quantum Ising chain. This finding indicates the existence of three regions near the QCP, each characterized by a specific asymptotic form of the correlation length, and demonstrates that the CP symmetry is restored at any nonzero temperature, entirely analogous to the quantum Ising chain. The range of the scaling behavior is also examined and found to be particularly wide.

Flow Rate Dependency of Steady-State Two-Phase Flows in Pore Networks: Universal, Relative Permeability Scaling Function and System-Characteristic Invariants

The phenomenology of steady-state two-phase flow in porous media is conventionally recorded by the relative permeability diagrams in terms of saturation. Yet, theoretical, numerical and laboratory studies of flow in artificial pore network models and natural porous media have revealed a significant dependency on the flow rates—especially when the flow regime is capillary to capillary/viscous and part of the disconnected non-wetting phase remains mobile. These studies suggest that relative permeability models should incorporate the functional dependence on flow intensities. In the present work, a systematic dependence of the pressure gradient and the relative permeabilities on flow rate intensity is revealed. It is based on extensive simulations of steady-state, fully developed, two-phase flows within a typical 3D model pore network, implementing the DeProF mechanistic–stochastic model algorithm. Simulations were performed across flow conditions spanning 5 orders of magnitude, both in the capillary number, Ca, and the flow rate ratio, r, and for different favorable /unfavorable viscosity ratio fluid systems. The systematic, flow rate dependency of the relative permeabilities can be described analytically by a universal scaling function along the entire domain of the independent variables of the process, Ca and r. This universal scaling comprises a kernel function of the capillary number, Ca, that describes the asymmetric effects of capillarity across the entire flow regime—from capillarity-dominated to mixed capillarity/viscosity- to viscosity-dominated flows. It is shown that the kernel function, as well as the locus of the cross-over relative permeability values, are single-variable functions of the capillary number; they are both identified as viscosity ratio invariants of the system. Both invariants can be correlated with the structure of the pore network, through a function of Ca. Consequently, the correlation is associated with the wettability characteristics of the system. Among the potential applications, the proposed, universal, flow rate dependency scaling laws are the improvement of core analysis and dynamic rock-typing protocols, as well as integration into field-scale simulators or associated machine learning interventions for improved specificity/accuracy.

Universal scaling of shear thickening transitions

Nearly, all dense suspensions undergo dramatic and abrupt thickening transitions in their flow behavior when sheared at high stresses. Such transitions occur when the dominant interactions between the suspended particles shift from hydrodynamic to frictional. Here, we interpret abrupt shear thickening as a precursor to a rigidity transition and give a complete theory of the viscosity in terms of a universal crossover scaling function from the frictionless jamming point to a rigidity transition associated with friction, anisotropy, and shear. Strikingly, we find experimentally that for two different systems—cornstarch in glycerol and silica spheres in glycerol—the viscosity can be collapsed onto a single universal curve over a wide range of stresses and volume fractions. The collapse reveals two separate scaling regimes due to a crossover between frictionless isotropic jamming and frictional shear jamming, with different critical exponents. The material-specific behavior due to the microscale particle interactions is incorporated into a scaling variable governing the proximity to shear jamming, that depends on both stress and volume fraction. This reformulation opens the door to importing the vast theoretical machinery developed to understand equilibrium critical phenomena to elucidate fundamental physical aspects of the shear thickening transition.

Role of finite probe size in measuring growth exponent in film deposition

We use computer simulations to investigate the effects of the tip diameter of an electrostatic force microscope (EFM) operating at a constant force on the extraction of the growth exponent β during film growing in a one-dimensional substrate. Laplace’s equation is solved in the EFM simulation using the finite element method to determine the electrostatic force between the tip and the film interface. Importantly, for EFM tips with sufficiently large apex diameters, the topographies calculated with EFM and those computed with the transformed mean height profile (TMHP) method, where the interface is divided into bins of the same tip diameter size and the average height within each bin is used to transform the original interface, are almost identical. This was shown in the context of lattice models of the Kardar–Parisi–Zhang (KPZ) and Villain–Lai–Das–Sarma (VLDS) classes. The global roughness of the film surface, W, scales with the diameter of the EFM tip, ε, as W/a=(ε/a)αg[Ψ], where a is the lattice parameter, α is the KPZ/VLDS roughness exponent, and g is a universal scaling function of the argument Ψ≡t/(ε/a)z, where t and z are the reduced time of deposition and the KPZ/VLDS dynamic exponent, respectively. These results provide a limit for ε from which a KPZ/VLDS growth exponent can be reliably determined with EFM at a constant force. When the EFM tip diameter is larger than the surface correlation length, a misleading effective growth exponent consistent with uncorrelated growth is found.

Sensitivity of the thermodynamics of two-dimensional systems towards the topological classes of their surfaces

Using Monte Carlo simulations we study the two-dimensional Ising model on closed manifolds of various topologies near the critical point of the planar bulk system. To this end we consider triangular, square, and hexagonal lattices. We find that the corresponding universal scaling functions of the magnetic susceptibility χ, as well as those of the specific heat C, differ for distinct Euler characteristics K of the manifold. In particular, for each lattice decoration the maxima of the scaling functions of χ grow as K increases. For the specific heat this relationship is inverted: e.g., for each lattice decoration the maxima of the scaling functions of C on the spherical surface (Euler characteristic K=2) is smaller than on the projective plane (K=1) which, in turn, is smaller than on the torus and on the Klein bottle (both with K=0). We find that if the aspect ratio of the lattices is conform to a square geometrical shape, the magnetic susceptibility scaling functions for different lattice types differ only by a non-universal amplitude.

Many-body quantum chaos in stroboscopically-driven cold atoms

In quantum chaotic systems, the spectral form factor (SFF), defined as the Fourier transform of two-level spectral correlation function, is known to follow random matrix theory (RMT), namely a ‘ramp’ followed by a ‘plateau’ in late times. Recently, a generic early-time deviation from RMT, so-called the ‘bump’, was shown to exist in random quantum circuits as toy models for many-body quantum systems. We demonstrate the existence of ‘bump-ramp-plateau’ behavior in the SFF for a number of paradigmatic and stroboscopically-driven 1D cold-atom models: spinless and spin-1/2 Bose-Hubbard models, and nonintegrable spin-1 condensate with contact or dipolar interactions. We find that the scaling of the many-body Thouless time tTh —the onset of RMT—, and the bump amplitude are more sensitive to variations in atom number than the lattice size regardless of the hyperfine structure, the symmetry classes, or the choice of driving protocol. Moreover, tTh scaling and the increase of the bump amplitude in atom number are significantly slower in spinor gases than interacting bosons in 1D optical lattices, demonstrating the role of locality. We obtain universal scaling functions of SFF which suggest power-law behavior for the bump regime in quantum chaotic cold-atom systems, and propose an interference measurement protocol.

Replica approach to the generalized Rosenzweig-Porter model

The generalized Rosenzweig-Porter model with real (GOE) off-diagonal entries arguably constitutes the simplest random matrix ensemble displaying a phase with fractal eigenstates, which we characterize here by using replica methods. We first derive analytical expressions for the average spectral density in the limit in which the size NN of the matrix is large but finite. We then focus on the number of eigenvalues in a finite interval and compute its cumulant generating function as well as the level compressibility, i.e., the ratio of the first two cumulants: these are useful tools to describe the local level statistics. In particular, the level compressibility is shown to be described by a universal scaling function, which we compute explicitly, when the system is probed over scales of the order of the Thouless energy. Interestingly, the same scaling function is found to describe the level compressibility of the complex (GUE) Rosenzweig-Porter model in this regime. We confirm our results with numerical tests.

Scaling of Average Avalanche Shapes for Acoustic Emission during Jerky Motion of Single Twin Boundary in Single-Crystalline Ni2MnGa

Temporal average shapes of crackling noise avalanches, U(t) (U is the detected parameter proportional to the interface velocity), have self-similar behavior, and it is expected that by appropriate normalization, they can be scaled together according to a universal scaling function. There are also universal scaling relations between the avalanche parameters (amplitude, A, energy, E, size (area), S, and duration, T), which in the mean field theory (MFT) have the form E∝A3, S∝A2, S∝T2. Recently, it turned out that normalizing the theoretically predicted average U(t) function at a fixed size, U(t)=atexp-bt2 (a and b are non-universal, material-dependent constants) by A and the rising time, R, a universal function can be obtained for acoustic emission (AE) avalanches emitted during interface motions in martensitic transformations, using the relation R~A1-φ too, where φ is a mechanism-dependent constant. It was shown that φ also appears in the scaling relations E~A3-φ and S~A2-φ, in accordance with the enigma for AE, that the above exponents are close to 2 and 1, respectively (in the MFT limit, i.e., with φ= 0, they are 3 and 2, respectively). In this paper, we analyze these properties for acoustic emission measurements carried out during the jerky motion of a single twin boundary in a Ni50Mn28.5Ga21.5 single crystal during slow compression. We show that calculating from the above-mentioned relations and normalizing the time axis of the average avalanche shapes with A1-φ, and the voltage axis with A, the averaged avalanche shapes for the fixed area are well scaled together for different size ranges. These have similar universal shapes as those obtained for the intermittent motion of austenite/martensite interfaces in two different shape memory alloys. The averaged shapes for a fixed duration, although they could be acceptably scaled together, showed a strong positive asymmetry (the avalanches decelerate much slower than they accelerate) and thus did not show a shape reminiscent of an inverted parabola, predicted by the MFT. For comparison, the above scaling exponents were also calculated from simultaneously measured magnetic emission data. It was obtained that the φ values are in accordance with theoretical predictions going beyond the MFT, but the AE results for φ are characteristically different from these, supporting that the well-known enigma for AE is related to this deviation.

Numerical investigation of localization in two-dimensional quasiperiodic mosaic lattice

A one-dimensional lattice model with mosaic quasiperiodic potential is found to exhibit interesting localization properties, e.g. clear mobility edges (Wang et al 2020 Phys. Rev. Lett. 125 196604). We generalize this mosaic quasiperiodic model to a two-dimensional version, and numerically investigate its localization properties: the phase diagram from the fractal dimension of the wavefunction, the statistical and scaling properties of the conductance. Compared with disordered systems, our model shares many common features but also exhibits some different characteristics in the same dimensionality and the same universality class. For example, the sharp peak at of the critical distribution and the large g limit of the universal scaling function β resemble those behaviors of three-dimensional disordered systems.

Critical Probability Distributions of the Order Parameter from the Functional Renormalization Group.

We show that the functional renormalization group (FRG) allows for the calculation of the probability distribution function of the sum of strongly correlated random variables. On the example of the three-dimensional Ising model at criticality and using the simplest implementation of the FRG, we compute the probability distribution functions of the order parameter or, equivalently, its logarithm, called the rate functions in large deviation theory. We compute the entire family of universal scaling functions, obtained in the limit where the system size L and the correlation length of the infinite system ξ_{∞} diverge, with the ratio ζ=L/ξ_{∞} held fixed. It compares very accurately with numerical simulations.

Universal scaling for disordered viscoelastic matter near the onset of rigidity.

The onset of rigidity in interacting liquids, as they undergo a transition to a disordered solid, is associated with a rearrangement of the low-frequency vibrational spectrum. In this Letter, we derive scaling forms for the singular dynamical response of disordered viscoelastic networks near both jamming and rigidity percolation. Using effective-medium theory, we extract critical exponents, invariant scaling combinations, and analytical formulas for universal scaling functions near these transitions. Our scaling forms describe the behavior in space and time near the various onsets of rigidity, for rigid and floppy phases and the crossover region, including diverging length scales and timescales at the transitions.

Unified modeling and experimental realization of electrical and thermal percolation in polymer composites

Correlations between electrical and thermal conduction in polymer composites are blurred due to the complex contribution of charge and heat carriers at the nanoscale junctions of filler particles. Conflicting reports on the lack or existence of thermal percolation in polymer composites have made it the subject of great controversy for decades. Here, we develop a generalized percolation framework that describes both electrical and thermal conductivity within a remarkably wide range of filler-to-matrix conductivity ratios (Yf/Ym), covering 20 orders of magnitude. Our unified theory provides a genuine classification of electrical conductivity with typical Yf/Ym≥1010 as insulator–conductor percolation with the standard power-law behavior and of thermal conductivity with 102≤Yf/Ym≤104 as poor–good conductor percolation characterized by two universal critical exponents. Experimental verification of the universal and unified features of our theoretical framework is conducted by constructing a 3D segregated and well-extended network of multiwalled carbon nanotubes in polypropylene as a model polymer matrix under a carefully designed fabrication method. We study the evolution of the electrical and thermal conductivity in our fabricated composites at different loading levels up to 5 vol. %. Significantly, we find an ultralow electrical percolation threshold at 0.02 vol. % and a record-low thermal percolation threshold at 1.5 vol. %. We also apply our theoretical model to a number of 23 independent experimental and numerical datasets reported in the literature, including more than 350 data points, for systems with different microscopic details, and show that all collapse onto our proposed universal scaling function, which depends only on dimensionality.

Creep failure of amorphous solids under tensile stress.

Applying constant tensile stress to a piece of amorphous solid results in a slow extension, followed by an eventual rapid mechanical collapse. This "creep" process is of paramount engineering concern, and as such was the subject of study in a variety of materials, for more than a century. Predictive theories for τ_{w}, the expected time of collapse, are incomplete, mainly due to its dependence on a bewildering variety of parameters, including temperature, system size, tensile force, but also the detailed microscopic interactions between constituents. The complex dependence of the collapse time on all the parameters is discussed below, using simulations of strip of amorphous material. Different scenarios are observed for ductile and brittle materials, resulting in serious difficulties in creating an all-encompassing theory that could offer safety measures for given conditions. A central aim of this paper is to employ scaling concepts, to achieve data collapse for the probability distribution function (pdf) of lnτ_{w}. The scaling ideas result in a universal function which provides a prediction of the pdf of lnτ_{w} for out-of-sample systems, from measurements at other values of these parameters. The predictive power of the scaling theory is demonstrated for both ductile and brittle systems. Finally, we present a derivation of universal scaling function for brittle materials. The ductile case appears to be due to a plastic necking instability and is left for future research.