Renormalization Group Predictions Research Articles

Experimental Test of Sinai's Model in DNA Unzipping.

The experimental measurement of correlation functions and critical exponents in disordered systems is key to testing renormalization group (RG) predictions. We mechanically unzip single DNA hairpins with optical tweezers, an experimental realization of the diffusive motion of a particle in a one-dimensional random force field, known as the Sinai model. We measure the unzipping forces F_{w} as a function of the trap position w in equilibrium and calculate the force-force correlator Δ_{m}(w), its amplitude, and correlation length, finding agreement with theoretical predictions. We study the universal scaling properties since the effective trap stiffness m^{2} decreases upon unzipping. Fluctuations of the position of the base pair at the unzipping junction u scales as u∼m^{-ζ}, with a roughness exponent ζ=1.34±0.06, in agreement with the analytical prediction ζ=4/3. Our study provides a single-molecule test of the functional RG approach for disordered elastic systems in equilibrium.

Puzzle of bicriticality in the XXZ antiferromagnet

Renormalization-group theory predicts that the XXZ antiferromagnet in a magnetic field along the easy Z-axis has asymptotically either a tetracritical phase-diagram or a triple point in the field-temperature plane. Neither experiments nor Monte Carlo simulations procure such phase diagrams. Instead, they find a bicritical phase-diagram. Here this discrepancy is resolved: after generalizing a ubiquitous condition identifying the tetracritical point, we employ new renormalization-group recursion relations near the isotropic fixed point, exploiting group-theoretical considerations and using accurate exponents at three dimensions. These show that the experiments and simulations results can only be understood if their trajectories flow towards the fluctuation-driven first order transition (and the associated triple point), but reach this limit only for prohibitively large system sizes or correlation lengths. In the crossover region one expects a bicritical phase diagram, as indeed is observed. A similar scenario may explain puzzling discrepancies between simulations and renormalization-group predictions for a variety of other phase diagrams with competing order parameters.

Tensor Network Calculation of the Logarithmic Correction Exponent in the XY Model

We study the logarithmic correction to the scaling of the first Lee-Yang (LY) zero in the classical $XY$ model on square lattices by using tensor renormalization group methods. In comparing the higher-order tensor renormalization group (HOTRG) and the loop-optimized tensor network renormalization (LoopTNR), we find that the entanglement filtering in LoopTNR is crucial to gaining high accuracy for the characterization of the logarithmic correction, while HOTRG still proposes approximate bounds for the zero location associated with two different bond-merging algorithms of the higher-order singular value decomposition and the oblique projectors. Using the LoopTNR data computed up to the system size of $L=1024$ in the $L \times L$ lattices, we estimate the logarithmic correction exponent $r = -0.0643(9)$ from the extrapolation of the finite-size effective exponent, which is comparable to the renormalization group prediction of $r = -1/16$.

Equilibrium to off-equilibrium crossover in homogeneous active matter

We study the crossover between equilibrium and off-equilibrium dynamical universality classes in the Vicsek model near its ordering transition. Starting from the incompressible hydrodynamic theory of Chen et al \cite{chen2015critical}, we show that increasing the activity leads to a renormalization group (RG) crossover between the equilibrium ferromagnetic fixed point, with dynamical critical exponent $z = 2$, and the off-equilibrium active fixed point, with $z = 1.7$ (in $d=3$). We run simulations of the classic Vicsek model in the near-ordering regime and find that critical slowing down indeed changes with activity, displaying two exponents that are in remarkable agreement with the RG prediction. The equilibrium-to-off-equilibrium crossover is ruled by a characteristic length scale beyond which active dynamics takes over. Such length scale is smaller the larger the activity, suggesting the existence of a general trade-off between activity and system's size in determining the dynamical universality class of active matter.

Lattice Abelian-Higgs model with noncompact gauge fields

We consider a noncompact lattice formulation of the three-dimensional electrodynamics with $N$-component complex scalar fields, i.e., the lattice Abelian-Higgs model with noncompact gauge fields. For any $N\ensuremath{\ge}2$, the phase diagram shows three phases differing for the behavior of the scalar-field and gauge-field correlations: The Coulomb phase (short-ranged scalar and long-ranged gauge correlations), the Higgs phase (condensed scalar-field and gapped gauge correlations), and the molecular phase (condensed scalar-field and long-ranged gauge correlations). They are separated by three transition lines meeting at a multicritical point. Their nature depends on the phases they separate, and on the number $N$ of components of the scalar field. In particular, the Coulomb-to-molecular transition line (where gauge correlations are irrelevant) is associated with the Landau-Ginzburg-Wilson ${\mathrm{\ensuremath{\Phi}}}^{4}$ theory sharing the same SU($N$) global symmetry but without explicit gauge fields. On the other hand, the Coulomb-to-Higgs transition line (where gauge correlations are relevant) turns out to be described by the continuum Abelian-Higgs field theory with explicit gauge fields. Our numerical study is based on finite-size scaling analyses of Monte Carlo simulations with ${C}^{*}$ boundary conditions (appropriate for lattice systems with noncompact gauge variables, unlike periodic boundary conditions), for several values of $N$, i.e., $N=2$, 4, 10, 15, and 25. The numerical results agree with the renormalization-group predictions of the continuum field theories. In particular, the Coulomb-to-Higgs transitions are continuous for $N\ensuremath{\gtrsim}10$, in agreement with the predictions of the Abelian-Higgs field theory.

Renormalization group flows for Wilson-Hubbard matter and the topological Hamiltonian

Understanding the robustness of topological phases of matter in the presence of interactions poses a difficult challenge in modern condensed matter, showing interesting connections to high energy physics. In this work, we leverage these connections to present a complete analysis of the continuum long-wavelength description of a generic class of correlated topological insulators: Wilson-Hubbard topological matter. We show that a Wilsonian renormalization group (RG) approach, combined with the so-called topological Hamiltonian, provide a quantitative route to understand interaction-induced topological phase transitions that occur in Wilson-Hubbard matter. We benchmark two-loop RG predictions for a quasi-1D Wilson-Hubbard model by means of exhaustive numerical simulations based on matrix product states (MPS). The agreement of the RG predictions with MPS simulations motivates the extension of the RG calculations to higher-dimensional topological insulators.

Nonuniversal critical aging scaling in three-dimensional Heisenberg antiferromagnets

We numerically investigate the stationary and non-equilibrium critical dynamics in three-dimensional isotropic Heisenberg antiferromagnets. Since the non-conserved staggered magnetization couples dynamically to the conserved magnetization density, we employ a hybrid simulation algorithm that combines reversible spin precession with relaxational Kawasaki spin exchange processes. We measure the dynamic critical exponent and identify a suitable intermediate time window to obtain the aging scaling exponents. Our results support an earlier renormalization group prediction: While the critical aging collapse exponent assumes a universal value, the associated temporal decay exponent in the two-time spin autocorrelation function depends on the initial distribution of the conserved fields; here, specifically on the width of the initial spin orientation distribution.

Local roughness exponent in the nonlinear molecular-beam-epitaxy universality class in one dimension.

We report local roughness exponents, α_{loc}, for three interface growth models in one dimension which are believed to belong to the nonlinear molecular-beam-epitaxy (nMBE) universality class represented by the Villain-Lais-Das Sarma (VLDS) stochastic equation. We applied an optimum detrended fluctuation analysis (ODFA) [Luis et al., Phys. Rev. E 95, 042801 (2017)2470-004510.1103/PhysRevE.95.042801] and compared the outcomes with standard detrending methods. We observe in all investigated models that ODFA outperforms the standard methods providing exponents in the narrow interval α_{loc}^{}∈[0.96,0.98] quantitatively consistent with two-loop renormalization group predictions for the VLDS equation. In particular, these exponent values are calculated for the Clarke-Vvdensky and Das Sarma-Tamborenea models characterized by very strong corrections to the scaling, for which large deviations of these values had been reported. Our results strongly support the absence of anomalous scaling in the nMBE universality class and the existence of corrections in the form α_{loc}^{}=1-ε of the one-loop renormalization group analysis of the VLDS equation.

Disorder and critical phenomena: the $$\alpha =0$$ α = 0 copolymer model

The copolymer model is a disordered system built on a discrete renewal process with inter-arrival distribution that decays in a regularly varying fashion with exponent $$1+ \alpha \;\geqslant \;1$$ . It exhibits a localization transition which can be characterized in terms of the free energy of the model: the free energy is zero in the delocalized phase and it is positive in the localized phase. This transition, which is observed when tuning the mean h of the disorder variable, has been tackled in the physics literature notably via a renormalization group procedure that goes under the name of strong disorder renormalization. We focus on the case $$\alpha =0$$ —the critical value $$h_c(\beta )$$ of the parameter h is exactly known (for every strength $$\beta $$ of the disorder) in this case—and we provide precise estimates on the critical behavior. Our results confirm the strong disorder renormalization group prediction that the transition is of infinite order, namely that when $$h\searrow h_c(\beta )$$ the free energy vanishes faster than any power of $$h-h_c(\beta )$$ . But we show that the free energy vanishes much faster than the physicists’ prediction.

Review of Recent Developments in the Random-Field Ising Model

A lot of progress has been made recently in our understanding of the random-field Ising model thanks to large-scale numerical simulations. In particular, it has been shown that, contrary to previous statements: the critical exponents for different probability distributions of the random fields and for diluted antiferromagnets in a field are the same. Therefore, critical universality, which is a perturbative renormalization-group prediction, holds beyond the validity regime of perturbation theory. Most notably, dimensional reduction is restored at five dimensions, i.e., the exponents of the random-field Ising model at five dimensions and those of the pure Ising ferromagnet at three dimensions are the same.

Restoration of dimensional reduction in the random-field Ising model at five dimensions.

The random-field Ising model is one of the few disordered systems where the perturbative renormalization group can be carried out to all orders of perturbation theory. This analysis predicts dimensional reduction, i.e., that the critical properties of the random-field Ising model in D dimensions are identical to those of the pure Ising ferromagnet in D-2 dimensions. It is well known that dimensional reduction is not true in three dimensions, thus invalidating the perturbative renormalization group prediction. Here, we report high-precision numerical simulations of the 5D random-field Ising model at zero temperature. We illustrate universality by comparing different probability distributions for the random fields. We compute all the relevant critical exponents (including the critical slowing down exponent for the ground-state finding algorithm), as well as several other renormalization-group invariants. The estimated values of the critical exponents of the 5D random-field Ising model are statistically compatible to those of the pure 3D Ising ferromagnet. These results support the restoration of dimensional reduction at D=5. We thus conclude that the failure of the perturbative renormalization group is a low-dimensional phenomenon. We close our contribution by comparing universal quantities for the random-field problem at dimensions 3≤D<6 to their values in the pure Ising model at D-2 dimensions, and we provide a clear verification of the Rushbrooke equality at all studied dimensions.

Finite-size scaling study of dynamic critical phenomena in a vapor-liquid transition.

Via a combination of molecular dynamics (MD) simulations and finite-size scaling (FSS) analysis, we study dynamic critical phenomena for the vapor-liquid transition in a three dimensional Lennard-Jones system. The phase behavior of the model has been obtained via the Monte Carlo simulations. The transport properties, viz., the bulk viscosity and the thermal conductivity, are calculated via the Green-Kubo relations, by taking inputs from the MD simulations in the microcanonical ensemble. The critical singularities of these quantities are estimated via the FSS method. The results thus obtained are in nice agreement with the predictions of the dynamic renormalization group and mode-coupling theories.

Long-range temporal correlations in the Kardar–Parisi–Zhang growth: numerical simulations

To analyze long-range temporal correlations in surface growth, we study numerically the (1 + 1)-dimensional Kardar–Parisi–Zhang (KPZ) equation driven by temporally correlated noise, and obtain the scaling exponents based on two different numerical methods. Our simulations show that the numerical results are in good agreement with the dynamic renormalization group (DRG) predictions, and are also consistent with the simulation results of the ballistic deposition (BD) model.

Monte Carlo investigation of critical properties of the Landau point of a biaxial liquid-crystal system.

Extensive Monte Carlo simulations are performed to investigate the critical properties of a special singular point usually known as the Landau point. The singular behavior is studied in the case when the order parameter is a tensor of rank 2. Such an order parameter is associated with a nematic-liquid-crystal phase. A three-dimensional lattice dispersion model that exhibits a direct biaxial nematic-to-isotropic phase transition at the Landau point is thus chosen for the present study. Finite-size scaling and cumulant methods are used to obtain precise values of the critical exponent ν=0.713(4), the ratio γ/ν=1.85(1), and the fourth-order critical Binder cumulant U^{*}=0.6360(1). Estimated values of the exponents are in good agreement with renormalization-group predictions.

Criterion for Many-Body Localization-Delocalization Phase Transition

We propose a new approach to probing ergodicity and its breakdown in quantum many-body systems based on their response to a local perturbation. We study the distribution of matrix elements of a local operator between the system's eigenstates, finding a qualitatively different behaviour in the many-body localized (MBL) and ergodic phases. To characterize how strongly a local perturbation modifies the eigenstates, we introduce the parameter ${\cal G}(L)=\langle \ln (V_{nm}/\delta) \rangle$, which represents a disorder-averaged ratio of a typical matrix element of a local operator $V$ to the energy level spacing, $\delta$; this parameter is reminiscent of the Thouless conductance in the single-particle localization. We show that the parameter ${\cal G}(L)$ decreases with system size $L$ in the MBL phase, and grows in the ergodic phase. We surmise that the delocalization transition occurs when ${\cal G}(L)$ is independent of system size, ${\cal G}(L)={\cal G}_c\sim 1$. We illustrate our approach by studying the many-body localization transition and resolving the many-body mobility edge in a disordered 1D XXZ spin-1/2 chain using exact diagonalization and time-evolving block decimation methods. Our criterion for the MBL transition gives insights into microscopic details of transition. Its direct physical consequences, in particular logarithmically slow transport at the transition, and extensive entanglement entropy of the eigenstates, are consistent with recent renormalization group predictions.

Quantum criticality of a spin-1 XY model with easy-plane single-ion anisotropy via a two-time Green function approach avoiding the Anderson–Callen decoupling

In this work we study the quantum phase transition, the phase diagram and the quantum criticality induced by the easy-plane single-ion anisotropy in a d-dimensional quantum spin-1 XY model in absence of an external longitudinal magnetic field. We employ the two-time Green function method by avoiding the Anderson–Callen decoupling of spin operators at the same sites which is of doubtful accuracy. Following the original Devlin procedure we treat exactly the higher order single-site anisotropy Green functions and use Tyablikov-like decouplings for the exchange higher order ones. The related self-consistent equations appear suitable for an analysis of the thermodynamic properties at and around second order phase transition points. Remarkably, the equivalence between the microscopic spin model and the continuous O(2)-vector model with transverse-Ising model (TIM)-like dynamics, characterized by a dynamic critical exponent z=1, emerges at low temperatures close to the quantum critical point with the single-ion anisotropy parameter D as the non-thermal control parameter. The zero-temperature critic anisotropy parameter Dc is obtained for dimensionalities d>1 as a function of the microscopic exchange coupling parameter and the related numerical data for different lattices are found to be in reasonable agreement with those obtained by means of alternative analytical and numerical methods. For d>2, and in particular for d=3, we determine the finite-temperature critical line ending in the quantum critical point and the related TIM-like shift exponent, consistently with recent renormalization group predictions. The main crossover lines between different asymptotic regimes around the quantum critical point are also estimated providing a global phase diagram and a quantum criticality very similar to the conventional ones.

Bose-Einstein condensation and critical behavior of two-component bosonic gases

We study Bose-Einstein condensation (BEC) in three-dimensional two-component bosonic gases, characterizing the universal behaviors of the critical modes arising at the BEC transitions. For this purpose, we use field-theoretical (FT) renormalization-group (RG) methods and perform mean-field and numerical calculations. The FT RG analysis is based on the Landau-Ginzburg-Wilson Phi4 theory with two complex scalar fields which has the same symmetry as the bosonic system. In particular, for identical bosons with exchange Z_2,e symmetry, coupled by effective density-density interactions, the global symmetry is Z_2e X U(1) X U(1). At the BEC transition it may break into Z_2,e X Z_2 X Z_2 when both components condense simultaneously, or to U(1) X Z_2 when only one component condenses. This implies different universality classes for the corresponding critical behaviors. Numerical simulations of the two-component Bose-Hubbard model in the hard-core limit support the RG prediction: when both components condense simultaneously, the critical behavior is controlled by a decoupled XY fixed point, with unusual slowly-decaying scaling corrections arising from the on-site inter-species interaction.

Monte Carlo study of anisotropic scaling generated by disorder.

We analyze the critical properties of the three-dimensional Ising model with linear parallel extended defects. Such a form of disorder produces two distinct correlation lengths, a parallel correlation length ξ(∥) in the direction along defects and a perpendicular correlation length ξ(⊥) in the direction perpendicular to the lines. Both ξ(∥) and ξ(⊥) diverge algebraically in the vicinity of the critical point, but the corresponding critical exponents ν(∥) and ν(⊥) take different values. This property is specific for anisotropic scaling and the ratio ν(∥)/ν(⊥) defines the anisotropy exponent θ. Until now, estimates of quantitative characteristics of the critical behavior for such systems have been obtained only within the renormalization group approach. We report a study of the anisotropic scaling in this system via Monte Carlo simulation of the three-dimensional system with Ising spins and nonmagnetic impurities arranged into randomly distributed parallel lines. Several independent estimates for the anisotropy exponent θ of the system are obtained, as well as an estimate of the susceptibility exponent γ. Our results corroborate the renormalization group predictions obtained earlier.

Emerging criticality in the disordered three-color Ashkin-Teller model

We study the effects of quenched disorder on the first-order phase transition in the two-dimensional three-color Ashkin-Teller model by means of large-scale Monte Carlo simulations. We demonstrate that the first-order phase transition is rounded by the disorder and turns into a continuous one. Using a careful finite-size-scaling analysis, we provide strong evidence for the emerging critical behavior of the disordered Ashkin-Teller model to be in the clean two-dimensional Ising universality class, accompanied by universal logarithmic corrections. This agrees with perturbative renormalization-group predictions by Cardy. As a byproduct, we also provide support for the strong-universality scenario for the critical behavior of the two-dimensional disordered Ising model. We discuss consequences of our results for the classification of disordered phase transitions as well as generalizations to other systems.

Diluted antiferromagnetic 3D Ising model in a field

We present numerical simulations for the diluted antiferromagnetic 3D Ising model in an external magnetic field (DAFF) at zero temperature. Our results are compatible with the DAFF being in the same universality class as the random field Ising model, in agreement with the renormalization group prediction.