2026-07-28 | | Total: 137
Although anyon condensation is a standard mechanism for relating topological orders, anyon condensation in symmetry-enriched topological (SET) phases is more intricate because the condensate must also be compatible with the global symmetry. We study symmetry-preserving anyon condensation in SET phases described by the enlarged Hu-Geer-Wu (HGW) string-net model with multifusion-category input data. We show that a $G$-preserving condensation is characterized by a compatible grading of the input multifusion category, and that this grading constructs the multifusion-category input of the child SET phase. To make this construction concrete, we consider the case where the relevant data come from a finite group extension $E$ of $G$ by $N$: an $E$-graded fusion category induces a $G$-SET input by passing to the quotient symmetry $G$, and the resulting input naturally carries a compatible $N$-grading that implements the further condensation inside the SET phase while preserving $G$. We illustrate the construction using three quantum-double examples: the trivial extension $\mathbb{Z}_2\times\mathbb{Z}_2$, the non-Abelian semidirect product $S_3$, and the nontrivial central extension $\mathbb{Z}_4$. The $\mathbb{Z}_4$ example further shows that symmetry fractionalized anyons do not obstruct symmetry-preserving condensation once the condensate is treated as a physical coherent state.
Fluctuation theorems for work provide universal constraints on nonequilibrium fluctuations, yet their quantum generalizations often rely on externally prescribed classical driving protocols. While for classical systems, fluctuation theorems have been extended to autonomous work, where the dynamics of the work source is subject to the backaction of the system, their generalization to the quantum regime is constrained by the uncertainty principle. Here, we extend fluctuation theorems for autonomous work from the classical regime to the quantum regime. By performing successive projective measurements over the work source and the system, we derive Jarzynski-type and Crooks-type fluctuation theorems for autonomous inclusive work from initial mixed thermal states. These relations are analogous to fluctuation theorems for autonomous work in the classical regime and explicitly incorporate the fluctuations of the work source. However, quantum noncommutativity prevents a consistent reduction to the nonautonomous counterparts, even in the limit of a large work source and correspondingly negligible backaction. By contrast, under the exclusive work definition, the nonautonomous limit is recovered when the measured observable of the work source commutes with its bare Hamiltonian and the backaction of the system on the work source is negligible. Our results are illustrated with the Dicke model, where a single-mode radiation field and an ensemble of two-level atoms act as the system of interest and the work source, respectively.
Hexagonal group-IV semiconductors have attracted increasing interest owing to their unconventional electronic and optical properties compared to the cubic diamond phase. However, the stabilization of these metastable allotropes in planar heterostructures remains a major challenge. In this work, we demonstrate the epitaxial growth of planar hexagonal germanium on non-basal $m$-plane CdS substrates by low-energy plasma-enhanced chemical vapor deposition. The role of growth temperature in the formation and stabilization of the hexagonal phase is investigated. X-ray diffraction, scanning transmission electron microscopy, and polarization-resolved Raman spectroscopy reveal the formation of epitaxial hexagonal germanium with the expected crystal symmetry. In particular, the Raman response exhibits the characteristic polarization selection rules of the E$_\text{2g}$ phonon mode of hexagonal Ge. Conversely, photoluminescence spectroscopy does not reveal any Ge-related emission feature. Combined transmission electron microscopy observations and atomistic modeling show that strain relaxation is governed by a limited set of dislocation mechanisms, which efficiently relieve most of the mismatch strain within a few nanometers from the interface and can involve localized cubic stacking insertions. At greater distances from the interface, the progressive loss of hexagonal order is increasingly dominated by stacking-fault disorder, particularly I3-type defects. These results establish CdS as a promising template for the planar stabilization of hexagonal Ge and provide insight into the defect mechanisms governing strain relaxation in metastable group-IV heterostructures.
The sign problem breaks the polynomial scaling of Quantum Monte Carlo methods. We alleviate it by rotation of the local basis of the Hilbert space, such that the phase of off-diagonal matrix elements of the rotated Hamiltonian is minimized. These minima coincide with optima of the average sign. We benchmark our method for frustrated Heisenberg antiferromagnets in one dimension and on the two-dimensional maple-leaf lattice. Our approach reveals more efficient bases with improved sign in a large part of the parameter space for all models considered, enabling us to reach lower temperatures, on par with state-of-the-art numerical linked cluster expansions which we directly compare to.
Our study provides the disorder-induced contribution to the orbital Hall conductivity in three-dimensional Weyl semimetals with broken time-reversal symmetry. Using the quantum kinetic approach, we analyse the impact of side-jump and skew scattering contributions to the system. The dependence of the orbital Hall conductivity on both disorder potential and the Fermi energy is explicitly demonstrated. Furthermore, we demonstrate that the higher-order disorder contribution, especially from the third power of disorder potential, dominates the orbital Hall conductivity under an oscillating electric field in a time-reversal symmetry broken Weyl semimetal, suppressing other scattering mechanisms, including the side jump contributions. We can enhance the extrinsic orbital Hall conductivity by tuning the strength of the disorder potential, applied energy, and choosing the system with appropriate Weyl node separation. Finally, our results are supported by numerical estimations and highlight potential experimental relevance for advancing orbitronics device technologies.
We report a new approach to interfacial equilibria among multiple solution phases in the phase-field method. It employs auxiliary non-conserved variables to describe the composition differences among different phases, and their temporal evolution is driven by the differences among different chemical diffusional potentials. It is reduced to the Wheeler-Boettinger-McFadden (WBM) model of equal interfacial chemical compositions and to the Kim-Kim-Suzuki (KKS) model of equal interfacial chemical diffusion potentials as two limiting cases. It can directly incorporate thermodynamic databases without any further approximations and simplifications. It is generally applicable to a wide range of problems involving composition evolution and chemical equilibria, including processes such as interdiffusion between two ordered phases not sharing any common composition range, which would pose difficulty to the existing treatments using WBM and KKS models.
We present a lattice model for morphogen-controlled branching morphogenesis which combines ideas and concepts from non-equilibrium physics and developmental biology. In this model, the stochastic occupation dynamics of cells is coupled with signaling molecules (morphogens) produced by the cells. We investigate growth patterns governed by morphogen concentration gradients, spanning regimes ranging from the diffusion-limited aggregation limit to the stochastic surface growth (Eden model) limit. Moreover, we introduce control over topology by a local operator and study growth, degrowth, and steady-state dynamics of branched patterns. The topology-preserving steady-state clusters exhibit a power-law scaling of radius of gyration with cluster size, yielding the exponents $0.68\pm0.01$ in square lattice and $0.67\pm 0.01$ in hexagonal lattice.
Altermagnetism, characterized by its ferromagnetism-like spin-splitting band structure and antiferromagnetism-like magnetic order, has garnered considerable attention recently. Although hole doping may promote magnetism in the debated altermagnet candidate RuO2, the evolution of its electronic and magnetic properties under hole doping remains poorly understood. Based on first-principles calculations, we employ quantum oscillations to study hole-doped RuO2. We find that hole doping can enhance spin splitting and reconstruct the Fermi surface in RuO2, which is revealed by the angle-dependent quantum oscillation frequency. By tracking a pair of closed Fermi-surface pockets, we identify a meaningful correlation between the magnetic moment of Ru and a quantum-oscillation-based signature of spin splitting. This correlation follows a quasi-linear trend over a broad doping range, which can be captured by a minimal two-dimensional d-wave altermagnetic model. In addition, the hole-doped RuO2 exhibits a transition from a nonmagnetic to an altermagnetic state via an intermediate state. The quasi-linear correlation through quantum oscillation and the distinct quantum oscillation frequency of the stable altermagnetic state can serve as useful signatures for identifying the altermagnetic state in RuO2. Our results provide a comprehensive framework for understanding hole-doped RuO2, offering new insights into altermagnetic transitions and their identification.
We study the real-frequency structure of the self-energy in the bilayer Hubbard model, using the dynamical cluster approximation. At half filling, the Mott insulator-band insulator (MI-BI) crossover involves a rearrangement of self-energy poles between the bonding and antibonding bands; these poles cross as the interlayer hopping $t_{\perp}$ is varied. Upon doping, this pole structure produces a band-selective pseudogap and enhances $s^{\pm}$-wave superconductivity. The order parameter is maximized near the MI-BI boundary, where low-energy anomalous self-energy poles develop simultaneously in both bands and cooperatively enhance the pairing. We further show that these self-energy poles can be interpreted as emergent fermionic excitations, offering an enhanced-pairing mechanism in common with the single-layer Hubbard model. The controllability of these poles through $t_{\perp}$ makes the bilayer system an unconventional platform for optimizing strongly correlated superconductivity.
Quasi-two-dimensional (2D) van der Waals (vdW) ferromagnets such as the series Fe$_{3-x}$GeTe$_2$, with a relatively high Curie temperature and robust metallicity, offer an ideal platform for investigating itinerant magnetism in reduced dimensions. Here, we present a comprehensive electron spin resonance (ESR) investigation of single-crystalline almost-stoichiometric Fe$_{3.03 \pm 0.03}$GeTe$_{2}$ across wide ranges of frequencies, temperatures, and magnetic fields to gain quantitative insights into its magnetic anisotropy and spin dynamics. Frequency-dependent ESR measurements establish Fe$_{3}$GeTe$_{2}$ as an easy-axis ferromagnet. Temperature-dependent high-field ESR reveals a large internal field that gradually decreases at higher temperatures. Remarkably, this internal field persists even above $T_\mathrm{C}$, evidencing short-range spin correlations in Fe$_{3}$GeTe$_{2}$. Analysis of spin-wave modes yields a strong uniaxial magnetocrystalline anisotropy $K_{\text{int}} \approx - 5 \times 10^6$ erg cm$^{-3}$ at 3 K and a large magnon gap $Δ(3 \mathrm{K})$ $\approx$ 87.8 $\pm$ 13.7 GHz ($\approx$ 0.363 $\pm$ 0.057 meV). Our study highlights that small variations in Fe content in Fe$_{3}$GeTe$_{2}$ lead to substantial changes in the magnon gap at low temperatures, indicating the extreme sensitivity of spin dynamics to the chemical composition. These results establish Fe$_{3}$GeTe$_{2}$ as a model vdW ferromagnet for exploring tunable anisotropies and magnon excitations in metallic 2D magnets.
Bayesian calibration of material constitutive parameters from multimodal mechanical test data is often limited by the need to specify a joint likelihood across measurement modalities that differ in dimensionality, noise structure, and physical units. The resulting posteriors are often broad or strongly correlated, causing standard Markov Chain Monte Carlo (MCMC) samplers to mix poorly. Here, we present an amortized, likelihood-free framework that combines Gaussian process (GP) surrogates with Conditional Flow Matching (CFM) to learn conditional posteriors over constitutive parameters directly from synthetic multimodal parameter--observation pairs, avoiding hand-crafted likelihoods and repeated MCMC sampling. Once trained, the GP--CFM model generates posterior samples for each new specimen at negligible cost. The utility of this novel approach is demonstrated in this paper by estimating the values of Young's modulus and yield strength from the early portion of the force--displacement ($F$--$D$) curve and a Digital Image Correlation (DIC)-based displacement field measured in a Small Punch Test (SPT). It is observed that the $F$--$D$ data alone produce broad posteriors, consistent with limited parameter discrimination in the global response. Adding the DIC-measured displacement field was seen to contract the posteriors and shift them towards the independently measured tensile reference values. This work establishes a robust likelihood-free framework for the inference of material constitutive parameters from multimodal data, demonstrated through SPT--DIC integration.
We investigate high-dimensional, non-linear dynamical systems when exposed to incoherent periodic inputs and Hebbian-like synaptic plasticity. Our findings reveal a striking phenomenon: depending on the interplay between the strength of the periodic drive and synaptic plasticity, the system's phase diagram can give rise to a region where, once both inputs are removed, the collective dynamics spontaneously settles into a limit cycle. This suggests that periodic drives can imprint lasting rhythmic patterns into the network through plasticity, effectively teaching it to oscillate on its own. Numerical simulations on finite size systems show that the limit cycle phase can be easily detected on single-sample trajectories, while averaged curves are affected by strong finite size effects due to sample-to-sample fluctuations of the period of the limit cycles.
We introduce a data-efficient neural-vision assisted method to automate crystallographic alignment during molecular beam epitaxy (MBE) growth. Trained on reflection high-energy electron diffraction (RHEED) patterns from only 15 CdTe structures, our model - enabled by physics-aware postprocessing - reliably infers crystallographic directions, replacing manual frame-by-frame inspection. To this end, we design, test, and critically compare neural-network architectures based on 2D and 3D ResNet configurations, both with and without postprocessing that leverages the physical constraints of RHEED image acquisition. Our work delivers (i) a fully trained neural system ready for closed-loop deployment in future CdTe growth experiments and (ii) a generalizable pipeline for new materials where access to diverse RHEED datasets is limited. More broadly, this study represents a step toward AI-driven MBE growth and demonstrates the potential of machine-learning-assisted automation in thin-film synthesis.
Long-wavelength Bogoliubov phonons in an ordinary condensate realize the standard acoustic Lorentz metric. We show that a condensate formed in a hyperbolic polariton band realizes a different collective geometry with an indefinite Bogoliubov metric whose spatial signature is inherited from the opposite signs of the band curvatures. This metric converts the acoustic light cone into a hyperbolic stability wedge, separating propagating quasiparticles from dynamically unstable ones. In a driven-dissipative condensate, gain saturation turns this metric relation into a non-Hermitian Bogoliubov dispersion. The zero-discriminant surface becomes an exceptional cone in the parametric space $(q_x,q_y,Δ_{\rm NH})$, appearing as an exceptional hyperbola at fixed gain saturation. Across this surface the Bogoliubov branches coalesce, the biorthogonal phase rigidity collapses, and the spectrum changes from propagation to amplified or overdamped dynamics. Our results identify hyperbolic polariton condensates as a controllable setting where non-Hermitian exceptional degeneracies are organized by an effective Bogoliubov metric.
We introduce a general framework to solve a class of combinatorial opti- mization problems, including the matching problem, the Traveling Salesman Problem, and also the minimum weight k-factor problem. By reformulating these problems as an arrangement model, we recast the optimization task into a grand-canonical ensemble, where chemical potentials are used to re- lax strict topological constraints. The analytical solution found can serve as a polynomial-time algorithm to compute an approximate minimum cost for arbitrary k and link-weight distributions. Our framework is complementary to existing approaches and reveals new connections between combinatorial optimization and the statistical physics of disordered systems.
We propose an average-current interferometer for the directional causal response of fractional quantum Hall edges. Three coherent quantum point contacts (QPCs) form a flux-enclosing loop, so the leading Aharonov-Bohm harmonic is cubic in tunneling. Interference between a direct transfer and a coherent two-step path resolves downstream and upstream propagation. For a local Laughlin edge at $ν=1/m$, the upstream coefficient vanishes exactly, while the downstream amplitude scales as $E^{3ν-2}$. A weak finite-range nonlocal density interaction spanning the tunneling points activates the upstream coefficient without creating an upstream mode. At low temperature and unresolved delay, its amplitude scales as $E^{2ν-1}$ and its aligned phase relative to the downstream reference is $-π(1-ν)/2+χ_a$ modulo $2π$, with $χ_a=0$ or $π$. Opposite cyclic voltage orderings isolate the two directions, while a folded same-filling Laughlin edge with a neutral weak link realizes a sign-tunable bridge. For a general Abelian edge, the device probes the directional content of the selected tunneling vertex. If $δ_\pm$ are its downstream and upstream weights, its nonzero directional amplitudes scale as $E^{3Δ_\ell-2}$, with $Δ_\ell=δ_++δ_-$, and obey $A^u_\ell/A^d_\ell=|\sin(πδ_-)/\sin(πδ_+)|$. When both are nonzero, positive-flux alignment gives, after removal of a known fixed sign, the relative phase $πΔ_\ell=2πh_\ell$. In the same-vertex coherent unresolved-flight regime, the scaling dimension and directional ratio determine the exchange angle $θ_\ell=π(δ_+-δ_-)$ modulo $2π$. The Aharonov-Bohm flux frequency additionally gives the charge, enabling separate extraction of quasiparticle charge, scaling dimension, and exchange angle.
The study of electron hydrodynamics provides a powerful framework for understanding transport in ultraclean conductors, yet experimental evidence has thus far been largely restricted to the linear-response regime. Here, we report the direct observation of a strongly nonlinear transport regime in a high-mobility two-dimensional electron system. By engineering devices with a Venturi-shaped wedge geometry specifically designed to enhance convective nonlinearities, we uncover a pronounced nonlinear voltage response and large diodicity in the current-voltage characteristics. Our experimental findings show quantitative agreement with a theoretical model that attributes the observed nonlinearity to the convective acceleration of the electron fluid, analogous to the Bernoulli effect. These results provide compelling evidence for the applicability of the hydrodynamic framework to two-dimensional electron transport and open new avenues for exploring nonlinear and preturbulent phenomena in solid-state systems.
As an homage to the memory of Hans Willem Capel, I present some personal reminiscences and thoughts that come to my mind in this special occasion.
Two-dimensional topological insulators (2DTIs) harbor spin-polarized edge states that are topologically protected by time-reversal symmetry against non-magnetic structural disorder. However, coupling to magnetic impurities breaks this symmetry, inducing backscattering and destroying perfect quantization. While the impact of isolated dilute magnetic impurities is well understood, the transport properties in the presence of dense, disordered ensembles of magnetic moments remain poorly understood. In this work, we develop an analytical framework, supported by extensive numerical simulations, that captures the behavior of edge transport in two-dimensional topological insulators (2DTIs) with a finite concentration of magnetic impurities. We predict the onset of Anderson localization and uncover an anomalous localization regime characterized by a sub-exponential decay of the conductance, scaling as $\ln {\cal G} \propto -\sqrt{L}$, where $L$ is the system length. Furthermore, we demonstrate that the transport exhibits a universal scaling behavior governed solely by the effective impurity concentration. Applying our model to Mn-doped HgTe quantum wells, we find excellent agreement with experimental data. These findings provide a theoretical foundation for understanding anomalous localization phenomena in magnetically doped topological phases.
Our work explores large-area CVD graphene as a scalable platform for photodetectors. We show that substrate-induced photogating in graphene/SiO$_2$/Si structures enables a photoconductivity signal in mono-, bi-, and trilayer CVD graphene. The devices respond to illumination with photon energies above the Si band gap and do not respond to telecom frequency, suggesting the photoconductivity mechanism through the photogating. Large sample area and gate voltage-dependent resistivity measurements allow us to prove the dominant role of the photogating directly, demonstrating similarity of the capacitive recharging current and the photoresponse. The sensitivity is the highest (995 A/W) for the monolayer graphene. Increasing the number of graphene layers reduces the sensitivity due to screening and parallel conduction, but can improve detectivity by lowering the noise level. The sensitivity of CVD graphene device depends on fabrication route of photolitography: device with metal contacts fabricated before graphene transfer show higher photoresponce. Thus, our work demonstrate the potential of industrially relevant CVD graphene for sensitive photogating-based photodetectors.
Sliding ferroelectrics provide a nonvolatile platform for the electrical control of unconventional magnetism through reversible interlayer sliding. However, the coupling between sliding ferroelectricity and hybrid-parity nonrelativistic spin splitting (NSS) remains largely unexplored. Here, we introduce a class of hybrid-parity sliding multiferroics in which the spontaneous ferroelectric polarization is coupled to certain NSS components through interlayer sliding, allowing these components to be reversibly switched in an electrical way. Symmetry analysis identifies coplanar magnets as natural platforms for realizing this form of sliding multiferroicity. First-principles calculations establish bilayer VBr$_2$ as a representative example, demonstrating the coupled reversal of the out-of-plane ferroelectric polarization ($\pm$0.12 pC/m) and the signs of both even- and odd-parity NSS components via an interlayer-sliding pathway with an ultralow barrier of 6 meV/f.u. The signs of these NSS components are locked to the sliding-switchable ferroelectric polarization and encoded in the spin-current responses, providing a signature of the coupled ferroic switching. Our findings expand the scope of sliding multiferroics and the functionality of sliding ferroelectrics for low-energy, nonvolatile logic devices.
Heterogeneous density functional theory (DFT) calculations, particularly plane-wave implementations, introduce systematic formation energy errors ranging from tens to hundreds of meV/atom, depending on the selection of exchange-correlation functionals, kinetic energy cutoffs, pseudopotentials, and dispersion corrections. As demonstrated by the MatPES dataset, identical structures can exhibit an average energy discrepancy of 107 meV/atom between PBE and r2SCAN calculations. Such method-dependent discrepancies hinder the integration of multi-source DFT data, greatly limiting the scale and quality of datasets for training robust materials AI models. Here, we resolve this fundamental data silo barrier via graph-based transfer learning. Leveraging 380,190 structurally paired PBE-r2SCAN entries from the MatPES database, we train a structure-aware graph neural network to predict cross-functional energy residuals and align inconsistent DFT energy scales. By adopting GPTFF model architecture, the model converts conventional PBE energies to r2SCAN-level accuracy with a mean absolute error of 14.3 meV/atom, compared with 18.2 meV/atom achieved by CHGNet. This versatile approach effectively upgrades massive legacy PBE datasets to high-precision r2SCAN standards. It enables reliable predictions of phase stability, battery voltage profiles, and reaction thermodynamics, while allowing the integration of multi-source DFT data to advance the development of high-performance materials foundation models.
We develop from first principles the hydrodynamics of a two-dimensional chiral fluid, i.e. one carrying a net microscopic angular-momentum (spin) field. Enforcing angular-momentum conservation without imposing stress-tensor symmetry, we derive the full form of the stress tensor and of the spin flux, and we show that the entire chiral response is generated from the classical Newtonian one by a single operation of direct physical origin --- the $90^{\circ}$ rotation through which chirality acts, mirrored at the particle level by transverse forces (Caprini & Marini Bettolo Marconi 2025). Applied to the irreducible (deviatoric) decomposition, the rotation assigns to each classical channel a chiral partner --- pressure to chiral pressure, bulk and shear viscosities to their odd counterparts, the spin-flux gradient to its rotated image --- one coefficient and one mechanical action per channel, with no further cross-couplings. In this representation the steady base states become elementary. Quiescent states are organised by a holomorphic chiral complex potential, the mechanical and chiral pressures forming a conjugate harmonic pair subject to a topological existence condition; inhomogeneous activity forces azimuthal flows; and a boundary-driven confined flow, the chiral Stokes cavity, obeys a modified Helmholtz--Poisson system, solved in closed form in a circular domain and numerically in a square one. A single dimensionless group controls both geometries and sets the crossover from a screened, single-vortex regime to a sequence of sign-reversing vortical structures. The theory yields quantitative predictions, amenable to direct comparison with experiments on air-fluidised chiral disks (López-Castaño et al. 2022), and recovers the phenomenological frameworks of the chiral-fluid literature as particular cases.
The vast chemical design space and complex, interdependent design variables make catalyst discovery for targeted properties highly labor- and resource-intensive. Although generative models have emerged as a promising solution, existing approaches are generally limited to single-property conditioning or narrow chemical spaces. Here, we present Catalyst Diffusion Transformer (CatDiT), a unified framework for inverse catalyst design that generates valid and novel structures ranging from intermetallic alloys to oxide surfaces. By learning compressed latent representations, CatDiT enables efficient training and rapid sampling while supporting simultaneous conditioning on adsorbate type, binding energy, and catalyst class. The model provides reliable control of discrete properties and directional control of continuous properties, enriching candidate pools for reaction-specific catalyst discovery. As a representative application, multi-conditional generation for the nitrogen reduction reaction (NRR) yields 28 density functional theory (DFT)-relaxed alloy candidates that satisfy the target activity window and lie above the pure-metal *N-*H scaling line, corresponding to a ~1.5-fold enrichment over the source distribution. These results establish CatDiT as a practical and scalable approach for property-directed catalyst inverse design and targeted catalyst generation.
Non-reciprocal interactions find broad applicability in non-equilibrium and living systems. Their canonical implementation involves asymmetric couplings between two entities, which generally induce spatio-temporal patterns and time-dependent steady states that break time-translational invariance, representing a clear deviation from equilibrium physics. Although their phenomenology is well understood, whether non-reciprocal interactions induce new universality classes, and if so, under what conditions, remains an open question. In the present work, we perform a field-theoretic renormalization group (RG) analysis of the dynamics of two non-reciprocally coupled $n$-vector order parameters possessing a $U(n)$ symmetry, generalizing previous results to order parameters with multiple components, a feature that has been shown to generate novel non-equilibrium critical behavior in certain non-reciprocal systems. To lowest order in $ε= 4-d$, we find that the non-reciprocal coupling is RG-irrelevant, and the critical behavior is governed by the equilibrium fixed point of the Model A universality class of Hohenberg and Halperin with $2n$ vector components, even though the transition is into a non-equilibrium state. Our results demonstrate that non-reciprocity alone may not be sufficient to induce novel universality classes.