2025-12-30 | | Total: 139
Gapless quantum phases can become distinct when internal symmetries are enforced, in analogy with gapped symmetry-protected topological (SPT) phases. However, this distinction does not always lead to protected edge modes, raising the question of how the bulk-boundary correspondence is generalized to gapless cases. We propose that the spatial interface between gapless phases -- rather than their boundaries -- provides a more robust fingerprint. We show that whenever two 1+1d conformal field theories (CFTs) differ in symmetry charge assignments of local operators or twisted sectors, any symmetry-preserving spatial interface between the theories must flow to a non-invertible defect. We illustrate this general result for different versions of the Ising CFT with $\mathbb{Z}_2 \times \mathbb{Z}_2^T$ symmetry, obtaining a complete classification of allowed conformal interfaces. When the Ising CFTs differ by nonlocal operator charges, the interface hosts 0+1d symmetry-breaking phases with finite-size splittings scaling as $1/L^3$, as well as continuous phase transitions between them. For general gapless phases differing by an SPT entangler, the interfaces between them can be mapped to conformal defects with a certain defect 't Hooft anomaly. This classification also gives implications for higher-dimensional examples, including symmetry-enriched variants of the 2+1d Ising CFT. Our results establish a physical indicator for symmetry-enriched criticality through symmetry-protected interfaces, giving a new handle on the interplay between topology and gapless phases.
Ba$_2$IrO$_4$ has been refined in the tetragonal $I4/mmm$ phase without octahedral rotations, and its physical properties have been interpreted in this high-symmetry structure. However, the dynamical stability of this undistorted phase has not previously been questioned. It is important to establish whether other lower-symmetry structures are energetically more favorable because octahedral rotations control electronic bandwidths and constrain which magnetic interactions are allowed by symmetry. Here I compute first-principles phonon dispersions of $I4/mmm$ Ba$_2$IrO$_4$ including spin-orbit interaction. I find a nearly-flat nondegenerate unstable branch along the Brillouin-zone boundary segment $XP$ associated with inplane rotations of the IrO$_6$ octahedra. Using group-theoretical analysis, I enumerate the symmetry-allowed distortions associated with the $X_2^+$ and $P_4$ instabilities and fully relax the resulting structures. Only five of the twelve possible distortions can be stabilized, and the energy gain scales with the number of layers that exhibit octahedral rotations: phases with rotations in every IrO$_6$ layer are lower by $-5.8$ meV/atom and are nearly degenerate with respect to the stacking phase. Electronic structure calculations show that these rotated phases host a narrow and well-separated half-filled $J_{\textrm{eff}} = 1/2$ manifold, whereas structures with rotations only in alternate layers have broader and more entangled bands. This motivates a reinvestigation of the crystal structure of Ba$_2$IrO$_4$ and indicates that octahedral rotations should be considered in modeling its correlated electronic and magnetic properties.
Two-band Hamiltonians provide a typical description of topological band structures, in which the eigenfunctions can be characterized by a %Bloch vector field whose winding number that defines an integer topological invariant. This winding number is quantized and protected against continuous deformations of the Hamiltonian. Here we show that the Bloch vector and its winding number can be directly related to the gradient of the energy dispersion. Since the energy gradient is proportional to the group velocity, our result establishes an experimentally accessible correspondence between the Bloch vector field and angle-resolved photoemission spectroscopy measurements. We discuss a mapping between the gradient of the energy dispersion and the Bloch vector. This implies a direct and measurable relation between two-band Hamiltonians and their underlying topological structures.
We present a symbolic implementation of recursion method for the dynamics of strongly correlated fermions on one-, two- and three-dimensional lattices. Focusing on two paradigmatic models, interacting spinless fermions and the Hubbard model, we first directly confirm that the universal operator growth hypothesis holds for interacting fermionic systems, manifested by the linear growth of Lanczos coefficients. Equipped with symbolically computed Lanczos coefficients and knowledge of their asymptotics, we are able to compute infinite-temperature autocorrelation functions up to times long enough for thermalization to occur. In turn, the knowledge of autocorrelation functions unlocks transport properties. We compute with high precision the charge diffusion constant over a broad range of interaction strengths, $V$. Surprisingly, we observe that these results are well described by a simple functional dependence $\sim 1/V^2$ universally valid both for small and large $V$. All results are obtained directly in the thermodynamic limit. Our results highlight the promise of symbolic computational paradigm where the most costly step is performed once and outputs symbolic results that can further be used multiple times to easily compute physical quantities for specific values of model parameters.
We present a diffuse neutron and x-ray scattering study of structural, spin- and charge-density-wave fluctuations in the electrical insulator La$_{2-x}$Sr$_x$NiO$_4$. This lamellar nickelate is an isostructural analogue of the high-temperature cuprate superconductor La$_{2-x}$Sr$_x$CuO$_4$, for which recent experiments uncovered evidence for unusual structural and superconducting fluctuations indicative of rare-region physics due to inherent inhomogeneity unrelated to common point disorder effects. We find closely analogous nanoscale orthorhombic fluctuation behavior in La$_{2-x}$Sr$_x$NiO$_4$, including exponential scaling of the diffuse scattering intensity and power-law scaling of the characteristic length with relative temperature. Moreover, our neutron and x-ray scattering data reveal similar behavior for short-range magnetic and charge fluctuations above the respective ordering temperatures. These observations indicate that rare-region effects are a generic feature of perovskite-related structures and lead to universal fluctuations of both structural and electronic degrees of freedom over extended temperature ranges.
We investigate the dynamics of a run-and-tumble particle in a double-well potential and demonstrate that, in stark contrast to Brownian particles, active dynamics can lead to strong ergodicity breaking. When the barrier height exceeds a critical threshold, the long-time position distribution depends crucially on the initial condition: if the particle starts within the basin of attraction of one well, it remains trapped there, while if it begins between the two basins, it can reach either well with a finite probability, which we compute exactly via hitting probabilities. Below the critical barrier height, ergodicity is restored and the system converges to a unique stationary distribution, which we derive analytically. Using this result, we also estimate the characteristic barrier crossing time and show that it violates Kramer's-Arrhenius law, and displays a divergence near the critical height following a Vogel-Fulcher-Tammann-like form with an anomalous exponent $1/2$.
Recent work that analyzed the effect of vacancy disorder on a short-range resonating valence bond spin liquid state of kagome-lattice antiferromagnets argued that such spin liquids are stable to vacancy disorder. The argument relied crucially on a numerical study that identified the following property of the site-diluted kagome lattice: maximum-density dimer packings (maximum matchings) of any connected component of such site-diluted kagome lattices have at most one unmatched vertex that hosts a monomer. Here, we provide an inductive proof of a stronger result that implies this property: If a connected cluster of such a lattice has an odd number of vertices, its Gallai-Edmonds decomposition~\cite{Lovas_Plummer_1986} has exactly one ${\mathcal R}$-type region that spans the entire connected cluster and hosts a single monomer of any maximum-density dimer packing. If on the other hand it has an even number of sites, it admits perfect matchings (fully-packed dimer coverings with no monomers) and its Gallai-Edmonds decomposition consists of a single ${\mathcal P}$-type region that spans the entire cluster. Our proof also applies to the site-diluted Archimedean star lattice, the site-diluted pyrochlore lattice (corner-sharing tetrahedra), the site-diluted hyperkagome lattice, and, more generally, to any lattice satisfying a certain local connectivity property. It does not apply to bond-diluted versions of such lattices.
The recent discovery of high-temperature superconductivity in La$_3$Ni$_2$O$_7$ under ambient-pressure in strained thin films raises the question of how superconductivity can be optimized through strain. In this work, we investigate the strain-dependent electronic structure and superconducting transition temperature ($T_c$) of La$_3$Ni$_2$O$_7$ using density functional theory combined with random phase approximation spin-fluctuation calculations. We find that biaxial strain acts as a tuning parameter for Fermi surface topology and magnetic correlations. Large tensile strain drives a Lifshitz transition characterized by a $d_{z^2}$ band crossing, leading to a sharp increase in the density of states and theoretical pairing strength. However, this is accompanied by a large increase in magnetic proximity, suggesting strong competition with spin-density-wave order. Conversely, under compressive strain, we identify a structurally selective $T_c$ enhancement restricted to the high-symmetry $I4/mmm$ phase. This effect is driven by the straightening of Ni--O--Ni bonds and the emergence of a $Γ$-centered hole pocket, yielding $T_c$ values consistent with recent thin-film experiments. Our results highlight the balance between structural symmetry, electronic topology, and magnetic instability in nickelates, and provides a theoretical framework for optimizing superconductivity via strain engineering.
We establish universal scaling laws and quantify aging in three-dimensional uniformly heated hard sphere granular gases through large-scale event-driven molecular dynamics ($N=500{,}000$). We report three primary quantitative discoveries: (i) The characteristic energy decay time exhibits a universal inverse scaling $τ_0 \propto ε^{-1.03 \pm 0.02}$ with the dissipation parameter $ε= 1 - e^2$. (ii) The steady-state temperature follows a precise power-law $T_{\mathrm{steady}} \propto ε^{-1.51 \pm 0.03}$, reflecting the non-linear balance between thermostat heating and collisional dissipation. (iii) The velocity autocorrelation function $\bar{A}(τ_w, τ)$ demonstrates pronounced aging, with decay rates $λ$ following a power-law slowing down $λ(τ_w) \propto τ_w^{-0.82 \pm 0.05}$. These results establish the first 3D quantitative benchmarks for aging in driven dissipative gases, where near-Gaussian statistics persist despite extreme structural clustering.
We propose a simple and accurate approach to estimate the random close packing (RCP) fraction of binary hard-disk mixtures. By introducing a parameter based on the reduced third virial coefficient of the mixture, we show that the RCP fraction depends nearly linearly on this parameter, leading to a universal collapse of simulation data across a wide range of size ratios and compositions. Comparisons with previous models by Brouwers and Zaccone demonstrate that our approach provides the most consistent and accurate predictions. The method can be naturally extended to polydisperse mixtures with continuous size distributions, offering a robust framework for understanding the universality of RCP in hard-disk systems.
The experimental discovery of fractional quantum anomalous Hall (FQAH) states in tunable moiré superlattices has sparked intense interest in exploring the interplay between topological order and symmetry breaking phases. In this paper, we present a comprehensive numerical study of this interplay through large-scale density matrix renormalization group (DMRG) simulations on a minimal two-band lattice model of twisted bilayer MoTe$_2$ at filling $ν=-2/3$. We find robust FQAH ground states and provide clear numerical evidences for anyon excitations with fractional charge and pronounced real-space density modulations, directly supporting the recently proposed anyon density-wave halo picture. We also map out the displacement field dependent phase diagram, uncovering a rich landscape of charge ordered states emerging from the FQAH, including a quantum anomalous Hall crystal (QAHC) with an integer quantized Hall conductance. We expect our work to inspire further research interest of intertwined correlated topological phases in moiré systems.
We report an emergent ac Josephson-like effect arising without external bias, driven by the interplay between nonreciprocity and nonlinearity in coupled condensates. Using a minimal model of three mutually nonreciprocally coupled condensates, we uncover a rich landscape of dynamical phases governed by generalized Josephson equations. This goes beyond the Kuramoto framework owing to inherent nonreciprocity and dynamically evolving effective couplings, leading to static and dynamical ferromagnetic and (anti)vortex states with nontrivial phase winding. Most strikingly, we identify an ac phase characterized by the emergence of two distinct frequencies, which spontaneously break the time-translation symmetry: one associated with the precession of the global U(1) Goldstone mode and the other with a stabilized limit cycle in a five-dimensional phase space. This phase features bias-free autonomous oscillatory currents beyond conventional Josephson dynamics. We further examine how instabilities develop in the ferromagnetic and vortex states, and how they drive transitions into the ac regime. Interestingly, the transition is hysteretic: phases with different winding numbers destabilize under distinct conditions, reflecting their inherently different nonlinear structures. Our work lays the foundation for exploring nonreciprocity-driven novel dynamical phases in a broad class of condensate platforms.
Recent years have seen a proliferation in investigations on Altermagnetism due to its exciting prospects both from an applications perspective and theoretical standpoint. Traditionally, altermagnets are distinguished from collinear antiferromagnets using the central concept of halving subgroups within the spin space group formalism. In this work, we propose the Fundamental Lemma of Altermagnetism (FLAM) deriving the exact conditions required for the existence of altermagnetic phase in a magnetic material on the basis of site-symmetry groups and halving subgroups for a given crystallographic space group. The spin group formalism further clubs ferrimagnetism with ferromagnetism since the same-spin and opposite-spin sublattices lose their meaning in the presence of multiple magnetic species. As a consequence of FLAM, we further propose a class of fully compensated ferrimagnets, termed as Alterferrimagnets (AFiMs), which can show alternating momentum-dependent spin-polarized non-relativistic electronic bands within the first Brillouin zone. We show that alterferrimagnetism is a generalization of traditional collinear altermagnetism where multiple magnetic species are allowed to coexist forming fully compensated magnetic-sublattices, each with individual up-spin and down-spin sublattices.
Nanoconfined liquid crystals (LCs) and their nanocomposites are driving the next generation of photonic applications. Consequently, deepening our understanding of mesophase stability, defect topology, and the dynamic response of LCs at the nanoscale requires the development of novel characterization approaches. This motivates us to perform in situ observations on model 4'-octyl-4-cyanobiphenyl (8CB) LC using liquid-phase scanning transmission electron microscopy (LP-STEM). We find that the electron beam induced consecutive phase changes from smectic A to nematic (SmA-N) and from nematic to isotropic (N-I). The kinetic dependence of the phase transition on dose rate shows that the time between SmA-N and N-I shortens with increasing rate, revealing the hypothesis that a higher electron dose rate increases the energy dissipation rate, leading to substantial heat generation in the sample. We report on the spontaneous formation of disclinations, ordering effects, and complete process reversibility. Radiolytic effects of the electron beam are discussed in detail, and additional experiments with external heating indicate that the observed phenomena are mainly thermal in nature. The results are supported by calculations of heat diffusion, suggesting the nanoconfined 8CB differs significantly in thermal properties compared to the bulk one. This is the first detailed study of LC phase transitions using LP-STEM, which paves the way for further studies of nanoconfined LCs and for the development of the technique for advanced LC materials research.
Phosphorus-based lubricant additives are used for protecting metallic contacts under boundary lubrication by forming surface films that reduce wear and friction. Despite their importance, the molecular mechanisms driving their friction-reducing effects remain unclear, especially for phosphate esters, whose molecular structure critically impact tribological behavior. In this study, we use machine learning-based molecular dynamics simulations to investigate the tribological performance of three representative phosphorus-based additives, Dibutyl Hydrogen Phosphite (DBHP), Octyl Acid Phosphate (OAP), and Methyl Polyethylene Glycol Phosphate (mPEG-P), on iron surfaces. The mPEG-P family is further analyzed by varying esterification degree and chain length. DBHP exhibits the lowest friction and largest interfacial separation, resulting from steric hindrance and tribochemical reactivity, as indicated by P-O bond cleavage and enhanced O-Fe interactions. In contrast, OAP and mPEG-P monoesters produce higher friction due to limited steric protection and reduced resistance to shear, leading to partial loss of surface coverage under extreme conditions. Within the mPEG-P family, multi-ester and longer-chain molecules significantly lower friction by maintaining larger separations, demonstrating that steric effects can outweigh surface reactivity under severe confinement. Overall, these results provide atomistic insights into how molecular architecture controls additive performance and support the design of phosphorus-based lubricants combining reactive anchoring with optimized steric structures for durable, low-friction interfaces.
We show that nonequilibrium spin injection into a superconductor can generate an anomalous supercurrent or induce a phase gradient, even for spin voltages below the superconducting gap. Our mechanism does not require breaking time-reversal symmetry in the effective superconducting Hamiltonian, but instead relies on nonequilibrium spin injection. We further demonstrate that superconductivity enhances spin injection due to the large quasiparticle density of states near the pairing gap, an effect that persists well below the gap. This contrasts with earlier works predicting the absence of spin injection at zero temperature and small spin voltages. Our results provide a natural explanation for long-standing experimental observations of spin injection in superconductors and predict novel effects arising from spin-charge coupling, including the electrical control of anomalous phase gradients in superconducting systems with spin-orbit coupling. These effects are broadly testable in a variety of materials and hybrid superconducting structures.
Germanium hole spin qubits based on strained Ge/SiGe quantum well have attracted much research attention due to the strong spin-orbit coupling. In particular, the strain dependence of the heavy-hole--light-hole mixing and thus the $g$-tensor anisotropy offer unique opportunities for acoustic driving and spin-phonon coupling. In this work we numerically simulate the coherent control of a Ge hole spin using surface acoustic waves. The periodic strain dynamically modulates the $g$-tensor matrix and causes fast spin rotation under a small acoustic amplitude. Moreover, we show a strong anisotropy and confinement dependence of the Rabi frequency coming from the phase-shifted longitudinal and shear strain components. Our work lays the foundations for acoustic-driven spin control and spin-phonon coupling using Ge hole spin qubits.
In this paper we review the recent progress in theoretical understanding of the peculiarities of photogalvanic phenomena, photon drag and inverse Faraday effects in superconductors and hybrid superconducting structures. Our study is based on the time-dependent Ginzburg-Landau (TDGL) theory with a complex relaxation constant which provides the simplest description of the mechanisms of the second-order nonlinear effects in the electrodynamic response and related mechanisms of generation of dc photocurrents, magnetic moment and switching between different current states under the influence of electromagnetic radiation of various polarization.
Using quasi-simultaneous synchrotron X-ray diffraction and tomography techniques, we have studied in-situ and in real-time the nucleation and co-growth dynamics of the peritectic structures in an Al-Mn alloy during solidification. We collected ~30 TB 4D datasets which allow us to elucidate the phases' co-growth dynamics and their spatial, crystallographic and compositional relationship. The primary Al4Mn hexagonal prisms nucleate and grow with high kinetic anisotropy -70 times faster in the axial direction than the radial direction. In all cases, a ~5 um Mn-rich diffusion layer forms at the liquid-solid interface, creating a sharp local solute gradient that governs subsequent phase transformation. The peritectic Al6Mn phases nucleate epitaxially within this diffusion zone, initially forming a thin shell surrounding the Al4Mn with an orientation relationship of {10-10}HCP // {110}O, [0001]HCP // [001]O. Such ~5 um Mn-rich diffusion layers also cause solute depletion at the liquid side of the liquid-solid interface, limiting further epitaxial phase growth, but prompting phase re-nucleation and branching at crystal edges, resulting tetragonal prism structures that no longer follow the initial orientation relationship. The anisotropic diffusion also led to the formation of core defects at the centre of both phases. Furthermore, increasing cooling rate from 0.17 to 20 °C/s can disrupt the stability of the solute diffusion zone, effectively suppressing the formation of the core defects and forcing a transition from faceted to non-faceted morphologies. Our work establishes a new theoretical framework for how to tailor and control the peritectic structures in metallic alloys through solidification processes.
Controllable quasiparticle radiation in two-dimensional (2D) semiconductors is essential for efficient carrier recombination, tunable emission, and modulation of valley polarization which are strongly determined by both the density and nature of underlying excitonic species. Conventional chemical doping techniques, however, often hinder the reversibility and density of excitonic charge states (exciton and trion) due to unfavorable interactions between dopant and 2D materials. In this work, efficient excitonic charge state conversion is achieved by doping monolayer WS$_2$ using water rinsed PVA and the quasiparticle densities are further enhanced by applying high periodic biaxial strain (up to 2.3%) through a 2D silica microsphere array. The method presented here enables nearly 100% reversible trion-to-exciton conversion without the need of electrostatic gating, while delivering thermally stable trions with a large binding energy of ~56 meV and a high free electron density of ~3$\times$10$^{13}$ cm$^{-2}$ at room temperature. Strain-induced funneling of the PVA-injected free electrons substantially increases the excitonic quasiparticle densities and boosts the trion emission by 41%. Overall, this approach establishes a versatile platform for excitonic charge state conversion and enhanced quasiparticle density in 2D materials, offering promising opportunities for future optical data storage, quantum-light and display technologies.
Theoretical simulation is helpful for accurate interpretation of experimental X-ray absorption near-edge structure (XANES) spectra that contain rich atomic and electronic structure information of materials. However, current simulation methods are usually too complex to give the needed accuracy and timeliness when a large amount of data need to be analyzed, such as for in-situ characterization of battery materials. To address these problems, artificial intelligence (AI) models have been developed for XANES prediction. However, instead of using experimental XANES data, the existing models are trained using simulated data, resulting in significant discrepancies between the predicted and experimental spectra. Also, the universality across different elements has not been well studied for such models. In this work, we firstly establish a crystal graph neural network, pre-trained on simulated XANES data covering 48 elements, to achieve universal XANES prediction with a low average relative square error of 0.020223; and then utilize transfer learning to calibrate the model using a small experimental XANES dataset. After calibration, the edge energy misalignment error of the predicted S, Ti and Fe K edge XANES is significantly reduced by about 55%. The method demonstrated in this work opens up a new way to achieve fast, universal, and experiment-calibrated XANES prediction.
At high-frequency a periodically-driven quantum spin-1/2 system can emulate a chiral spin liquid (CSL) described by an effective static local chiral hamiltonian. In contrast, at low-frequency {\it anomalous} CSL can be realized in Swap Models, in which one-way spin transport occurs at the edge although the bulk time-evolution operator over one period is trivial. In this work we explicitly construct a family of Floquet quantum spin-1/2 models on the square lattice implementing Swap Models to investigate the stability of the anomalous CSL under frequency detuning and the transition to the high-frequency regime. We have used the average-energy spectrum on finite-size torus and cylinders to unfold the Floquet quasi-energy spectrum over the whole frequency range and obtain the geometrical Berry phases. This enabled us to identify three regimes upon increasing detuning: i) a finite-size regime (with no folding of the Floquet spectrum), ii) an intermediate (narrow) regime with folding and very few resonances and iii) a regime with an increased density of resonances suggesting heating. At small detuning, edge modes are revealed by spectroscopic tools and from the diamagnetic response of the system giving access to the anomalous winding number. The analysis of all the data suggests that the anomalous CSL is not continuously connected to the high-frequency CSL. We also discuss the possible occurrence of a long-lived prethermal anomalous CSL.
Ultrafast electric manipulation of magnetic order in solids is critical for the development of future terahertz data processing. A fascinating concept for such high-speed operation is offered in metallic antiferromagnets by Néel spin-orbit torque. It should allow one to coherently rotate the ordered spins by simply applying an electric current of suitable amplitude and polarity. However, such switching has been severely hampered by competing heat-induced effects, and it has not yet been achieved on the intrinsically ultrafast time scales of antiferromagnets. Here, we report robust, direction-controlled and non-thermal rotation of the Néel vector $\mathbf{L}$ by $\pm$90° at room temperature in the antiferromagnet Mn$_2$Au driven by phase-locked terahertz current pulses. All observed features are consistent with ultrafast Néel spin-orbit torque: First, nonlinear optical imaging reveals that the terahertz current direction sets the final orientation of $\mathbf{L}$ in the absence of any bias field for at least two months. Second, transient optical birefringence shows that the switching proceeds ultrafast in less than 15 picoseconds. Finally, atomistic spin-dynamics simulations reproduce the observed dynamics and confirm the minor role of thermal effects. While the switching is already one order of magnitude faster than in ferromagnets at comparable dissipated energy, our simulations predict routes toward switching times and energies which are another order of magnitude lower. Our approach can be transferred to electric-field-driven switching in many more antiferromagnets, including magnetoelectric insulators. The engineering of spin torques, resonance frequencies and read-out mechanisms provides an exciting pathway toward on-chip applications of terahertz antiferromagnetic spin-orbitronics.
The quantum anomalous Hall (QAH) effect holds fundamental importance in topological physics and technological promise for electronics. It is generally believed that the QAH effect can only be realized in insulators. In this Letter, we theoretically demonstrate that the QAH effect can also be realized in metallic systems, representing a phase distinct from the conventional QAH phase in insulators. This phase is characterized by the coexistence of chiral edge channels and isotropic bulk conduction channels without a bulk energy gap. Notably, in a six-terminal Hall bar, our calculations show that, the quantized Hall conductivity and nonzero longitudinal conductivity can emerge due to dephasing, despite the Hall resistivity itself never becoming quantized. Furthermore, the quantized Hall conductivity exhibits remarkable robustness against disorder. Our findings not only extend the range of materials capable of hosting the QAH effect from insulators to metals, but also provide insights that may pave the way for the experimental realization of the QAH effect at elevated temperatures.
Order versus disorder in the structure of materials plays a key role in the theoretical prediction of their properties. However, this structural description appears to be ineffective for new families of materials such as high entropy alloys (HEAs), which combine crystallographic order with chemical disorder. Here, we demonstrate for five-element HEAs as pure solid solutions that the chemical disorder of the elements decorating their cubic structure underlies the generation of second optical harmonics, overcoming the theoretical limit imposed on centrosymmetric crystals. Moreover, we discover that this disorder, inherent to HEAs, sets a threshold for non-linear light emission from the 4th to the 26th order. As a consequence of the 0.5 eV broadening of the energy levels of the five elements of the HEA, the emission spectrum covers broad visible (400-650 nm) and infrared (800-1600 nm) ranges. In addition to the challenge of theoretically predicting non-linear effects in unconventional materials, the duality of structural order and chemical disorder in HEAs offers the opportunity to design sustainable alternatives to urgently needed optical materials.