2026-05-07 | | Total: 82
Altermagnetism has so far mainly been understood in its even- and odd-parity forms. We show that collinear antiferromagnets with zero net magnetization can also host mixed-parity spin splitting, namely neither purely even nor purely odd in momentum. We identify the symmetry conditions for such mixed-parity altermagnetism and show that, in two dimensions, it can arise in spin-orbital magnets when the two antiparallel spin sectors are related by a single mirror symmetry. Using a two-sublattice two-orbital model, we demonstrate that circularly polarized light induces mixed-parity altermagnetism at finite staggered potential and odd-parity spin-orbital altermagnetism at zero staggered potential. Mixed-parity altermagnetism thereby emerges as the intermediate spin-split regime between even- and odd-parity altermagnetism when spin splitting and zero net magnetization are maintained. Spin-resolved orbital Edelstein effects provide a complementary electrical probe of the underlying spin-orbital order.
Rhombohedral graphene exhibits an exceptionally diverse array of correlated phases that depend sensitively on the displacement field. Compiling reported phases into a unified phase diagram reveals a pronounced field-dependent electron-hole asymmetry: correlated states on the hole-doped side emerge at small displacement fields, whereas the fractional quantum anomalous Hall effect (FQAHE) is observed exclusively on the electron-doped side under large displacement fields. This stark asymmetry highlights the need to understand how flat bands evolve with displacement fields. Here, we directly visualize the field-induced electron-hole asymmetric band flattening in rhombohedral pentalayer graphene (R5G) using nanospot angle-resolved photoemission spectroscopy with electrostatic gating. Beyond gap opening and spectral weight redistribution indicative of layer polarization, the gating field drives a strongly asymmetric modification of the flat bands: the flat valence band (FVB) evolves into an M-shaped dispersion at high field, whereas the flat conduction band (FCB) progressively flattens with increasing field. Comparison with calculations identifies critical parameters governing the band curvature of R5G, from which the resulting finite Berry curvature and near-ideal quantum geometry support the emergence of topological phases under electron doping at large fields. These results establish a direct link between the asymmetric phase diagram, band structure evolution, and quantum geometry, providing a microscopic framework for understanding correlated and topological phases in rhombohedral graphene.
Spin-orbital generalization of Kitaev model provides a robust extension to the original Kitaev model. However, real materials often exhibit competing interactions that break exact solvability which can give rise to new phases. Motivated by recent microscopic proposals of coexisting Yao-Lee and Kitaev couplings, we investigate the fate of the ground state when two independent exactly solvable spin liquid Hamiltonians each originally formulated on different lattice geometries are combined on a common lattice environment. We first focus on the hybrid Kitaev's honeycomb and square-lattice model. Using self-consistent mean-field analysis and perturbative calculation, we show that the strong-Kitaev regime yields magnetic order in the spin sector, while the orbital sector retains its topological order. We further analyze the hybridization of the Yao-Lee and square-lattice models and find that the model exhibits a rich evolution of Majorana Dirac bands and Lifshitz transitions. Remarkably, when the Yao-Lee and square-lattice couplings are equal and opposite, the model restores its exact solvability with a single itinerant Majorana flavor. These results demonstrate that hybrid spin liquid platforms may host various emergent phases beyond conventional exactly solvable limits.
Hopf bifurcations are a universal route to self-sustained oscillations in driven systems. Despite the absence of any singular stationary state, we show that time-averaged observables generically exhibit singularities at the onset of oscillations. The origin of this behavior is geometric: phase averaging over the emergent periodic attractor eliminates odd powers of the oscillation amplitude, while the squared amplitude varies smoothly with the distance from the bifurcation. Consequently, the excess of any smooth time-averaged observable admits an integer-power expansion; observables remain finite but display discontinuities in finite-order derivatives. This yields an Ehrenfest-like hierarchy of Hopf singularities, in which the first nonanalytic derivative is determined by the lowest-order coupling between the observable and the limit-cycle waveform that survives phase averaging. Generic observables therefore exhibit kink singularities, while symmetry or geometric cancellations can suppress lower-order couplings and shift nonanalyticity to higher derivatives. We demonstrate this mechanism in chemical, electronic, and climate oscillators. Our results identify supercritical Hopf bifurcations as a universal mechanism for nonanalytic observable behavior, where singular features arise without any underlying singular stationary state. They thus provide a generic setting for singular behavior without divergence.
We study dynamical spin correlations in a dissipative XXZ spin chain subject to uniform local spin-loss and pumping. Starting from a mixed steady state that is featureless albeit possessing finite magnetization, rich dynamics emerges in time-dependent two-point correlators evaluated on top of it. For unitary evolution in which the reservoir is absent, the longitudinal correlators reproduce the established hierarchy of spin-transport universality classes - ballistic, Kardar-Parisi-Zhang (KPZ) superdiffusive, and diffusive - across the phase diagram. However, for finite magnetization, additional ballistic light cone propagation gets superimposed on the previous universality classes, arising from magnon propagation.The transverse correlator displays very fast, exponential decay of correlations without wavefront propagation in the easy-plane case. At the isotropic point, it follows KPZ scaling due to $SU(2)$ symmetry, while in the easy-axis regime, it is characterized by ballistic spreading of correlations. Under full Lindbladian dynamics, the universality classes are preserved at early times, while the correlations acquire an overall exponential damping in the long-time limit. In terms of methods, we have used vectorized TEBD for numerical simulations and exact analytical results obtained via a Pfaffian representation and the third-quantization framework for the noninteracting XX case.
Engineering the dispersion relation is one of the key ingredients enabling the application of spin waves in computational elements. One way to engineer the spin-wave band structure is to create an artificial magnonic crystal, which can be used to design specific band gaps or dispersion branches. However, creating a two-dimensional magnonic crystal usually requires removing material, which dramatically decreases the decay lengths of spin waves. Here, we present a method to manipulate the demagnetizing field landscape by utilizing large-area curvilinear nanotemplates consisting of three-dimensional nanopyramids arranged in a square lattice with a period of 400 nm. In a 50-nm-thick Permalloy film grown on these curvilinear templates, we experimentally observe a complete in-plane band gap together with flat-band modes that exhibit strong real-space localization of the spin waves in the pyramid valleys. Micro-focused Brillouin light scattering measurements corroborate the numerically predicted dispersion and reveal the possibility of opening and closing this gap by varying the external magnetic field. Our results establish three-dimensional-templated continuous films as a versatile platform for two-dimensional signal processing and magnonic computing elements.
We prove several rigorous results about the asymptotic behaviour of the numbers of tadpoles (or lassos) embedded in a lattice, including cases where the head is subject to a constraint like being unknotted, or where the tail pierces the surface spanned by the head.
When a quantum system exhibiting a second order phase transition is quenched across the critical point in large but finite time, the dynamics are not adiabatic in the critical region and the Kibble-Zurek (KZ) mechanism provides a framework to determine local observables such as the mean defect density. However, to find higher-point functions, one has to go beyond the KZ paradigm asshown in recent works on one-dimensional transverse field Ising model (TFIM) following a linear quench. It has been found that (i) besides the KZ scale, the quench dynamics depend on another length scale that arises due to the finite phase difference between the low energy modes, and (ii) contrary to the expectations based on the KZ mechanism, in general, the correlation functions do not decay exponentially with distance. Motivated by these results for the linear quench, we are interested in understanding if these properties are universal, and consider the 1D TFIM when the transverse field varies algebraically in the vicinity of the critical field. We focus on the equal-time,longitudinal kink-kink correlation function at the end of the quench from the paramagnetic to the ferromagnetic phase, and find that (i) the correlator depends only on the KZ length for superlinear quenches, otherwise an additional dephasing length is required to describe it, and (ii) the dephased correlator decays as a compressed exponential with an exponent that changes continuously with the quench exponent. Our results are obtained using an adiabatic perturbation theory, analytical arguments and exact numerical integration of the relevant equations.
Materials science data collection can be expensive, making the reuse and long-term utility of datasets critical important for future discovery campaigns. In practice, researchers prioritize a subset of properties due to research interests. However, ignoring a subset of outcomes in data collection campaigns potentially generate datasets poorly suited for future learning tasks. Here, we present a framework for dataset construction that maximizes informativeness for target properties of interest while preserving performance on untargeted ones. Our approach uses diversity-aware selection to ensure broad coverage of the materials space. In noisy experimental dataset construction, we find that without our diversity-aware framework, prediction performance on untargeted properties can degrade by up to 40% relative to random sampling, whereas applying our framework yields improvements of up to 10% . For targeted properties, performance can degrade with respect to random sampling by up to 12.5% without diversity, while our framework achieves gains of up to 25%. Incorporating diversity into dataset construction not only preserves informativeness for the targeted properties, but also improves materials coverage for potential future objectives. As a result, the constructed datasets remain broadly informative across considered and unconsidered outcomes, ensuring unbiased quality entries and mitigating cold-start limitations in subsequent modeling and discovery campaigns.
Dissipationless charge transport without any energy loss is one of the most fascinating phenomena in condensed matter physics. This extraordinary state manifests in two well-established systems: superconductors and quantum Hall systems. A proposed third category is associated with chiral loop current order, characterized by the spontaneous formation of microscopic electric current loops. The microscopic origin of these currents stems from imaginary hopping terms, conceptualized as an imaginary charge density wave (iCDW). Despite extensive investigations, its existence remains highly controversial. Here we report site-selective spectroscopic evidence for a pure iCDW in the kagome nonmagnetic metal CsV$_3$Sb$_5$. Nuclear quadrupole resonance spectra at out-of-plane $^{121}$Sb site sensitive to in-plane currents reveal anomalous broadening below $T^*\approx$120 K, coinciding with the nematic transition well above the real charge density wave (CDW). Under magnetic fields, the spectra exhibit asymmetric lineshapes, demonstrating that this broadening purely originates from magnetic effects rather than from electric quadrupolar effects associated with CDW fluctuations. The observed lineshapes are quantitatively consistent with ~1 mT local fields induced by chiral loop currents, indicating spontaneous time-reversal symmetry breaking. This microscopic identification of the long-sought pure iCDW establishes a novel form of quantum order, potentially revolutionizing our understanding of exotic electronic states in quantum materials.
Using Hartree-Fock mean-field approach, we study the role of on-site ($U$) and extended ($V$) Hubbard interactions on the existence and evolution of edge modes in a half-filled Su-Schrieffer-Heeger (SSH) chain. We analyze the energy spectrum, local probability amplitudes, and site-resolved charge and spin density profiles across topological, critical, and trivial hopping regimes. We find that the localization of bound states is controlled by the ratio $2V/U$, with edge spin-density-wave modes for $U>2V$ and mid-chain charge-density-wave domain walls for $U<2V$, independent of band topology. These results establish the correlation-driven origin of localized states in finite one-dimensional chains.
We develop a unified fluctuation-response theory in the frequency domain for nonequilibrium steady states governed by overdamped Langevin dynamics and Markov jump processes. The relation expresses the power spectrum of general observables exactly as a quadratic form of local responses measured at the same frequency, thereby extending static nonequilibrium fluctuation-response relations to finite frequencies. The decomposition is spatial for Langevin systems and edge-resolved for Markov jump processes, and applies uniformly to state-dependent observables, current-like observables, and their combinations. As consequences of the same identity, we derive frequency-domain response uncertainty relations, kinetic and thermodynamic uncertainty relations, the equilibrium fluctuation-dissipation theorem, and Harada-Sasa-type relations. Applications to stochastic networks and driven diffusive systems illustrate how the theory resolves fluctuation spectra into edge-wise contributions and reveals frequency-dependent tradeoffs between fluctuations, response, and dissipation.
When calculating residual strain via neutron or X-ray diffraction, uncertainties propagated from the peak fit are often inadequate to describe the true scatter of measurements about a singular strain state, such as one that should describe a macroscopic continuum. Because diffraction is inherently a selective process, orientation dependent scatter arises from the sub-sampling of strong microstructure and strain gradients. This paper investigates the appropriateness of propagated uncertainties with reference to their original intention, i.e., noise about a mean value. Thirty-six unique orientations of strain measurements are taken at multiple locations within an additive friction-stir deposition component with fine-scale gradients (~200 um) of plastic strain, texture, and residual elastic strain. Multiple strain and stress calculation pathways are compared: direct substitution of three measurements into Hooke's law, direct inversion of any six unique orientations into the strain state tensor, and thirty-six measurement least-squares estimation. For the latter two cases, the appropriateness of the uncertainty interval is statistically evaluated based on a physical constraint: common agreement under the strain transformation law. For this sample, the direct inversion of six measurements retains a conservative estimate of the uncertainty. However, propagated uncertainties in the least-squares solution greatly underestimate the true experimental scatter. A simple pathway to estimate appropriate uncertainty intervals is suggested. These results demonstrate that interpretation of uncertainty in residual strain is strongly dependent on intrinsic, sample-dependent effects, and that oversampling orientations and statistical analysis can give more accurate results with realistic uncertainties.
We study how the topological properties of a one-dimensional staggered lattice, loaded into states with orbital angular momentum $l=1$, can be controlled simply by tuning the relative angle between sites. The original system under consideration can be depicted as a Creutz ladder model when unwrapping the different state circulations in a synthetic dimension. Depending on the hopping strengths of the chain, different topological regimes may be accessed by changing the ladder angle, as determined by the value of the winding number of the chain. We analytically and numerically explore the different available regimes, and determine the number of topologically protected edge states that exist in each case. We also study the emergence of band inversion across topological transitions and show that it agrees with the winding number calculations, thus serving as an additional topological marker. Then, we propose a realistic experimental implementation in a photonic waveguide system, where the topological transition manifests as a sudden change of the behavior of the propagation of light in the system.
All solids, whether crystalline or disordered, support elastic wave propagation with a linear dispersion relation in the long-wavelength limit. These waves, corresponding to low-frequency phonons, feature a vibrational density of states that follows Debye's classical model. Deviations from Debye's predictions with increasing frequency can emerge from phonon dispersion nonlinearity and from non-phononic vibrational modes, which exist in non-crystalline solids due to structural disorder. Both nonlinear phonon dispersion in disordered solids and its relative contribution to non-Debye anomalies, most notably manifested by the controversial boson peak, remain poorly understood. Here we show that nonlinear phonon dispersion in a broad range of disordered solids, including elastic networks and various glasses, emerge from a mesoscopic, disorder-induced lengthscale, which also controls wave attenuation. We subsequently use analysis and large-scale computer simulations to quantitatively determine the relative contributions of nonlinear phonon softening and non-phononic vibrations to the onset of non-Debye anomalies and to the boson peak. We show that the relative magnitude of the two contributions strongly depends on the strength of disorder of the solid, e.g., controlled by the thermal history upon glass formation, and that for realistic laboratory glasses both pieces of physics significantly contribute to the boson peak. These findings constitute basic progress in understanding disordered solids.
We derive a discrete spectral representation of the single-particle self-energy using a discrete evaluation of Kugler's symmetric improved estimator. Our construction can be used on both the real and the complex (Matsubara) frequency axis. It is guaranteed to remain causal at the numerical level, in contrast to standard approaches that may generate unphysical negative spectral weight or require additional broadening. Our representation can be used for any Hamiltonian; here we apply it to quantum impurity models and in dynamical mean-field theory. The latter is formulated with a discrete hybridization function throughout its self-consistency loop. In both cases and across various numerical methods, we obtain significantly improved accuracy for a range of impurity properties.
Using B2 MoTa, the strongest B2 former among group V/VI refractory binaries, as a model system, we compute temperature-dependent Gibbs energies of formation of vacancies and antisites from ab initio calculations up to 3000 K at the stoichiometric composition. We explicitly account for thermal electronic excitations, vibrational anharmonicity, and electron-vibration coupling. The key finding is that the temperature dependence of the Gibbs energies of vacancy formation exhibits a pronounced sublattice asymmetry. Specifically, the Gibbs energy of formation of a Mo-site vacancy decreases by 1.1 eV from 0 to 3000 K, whereas the decrease for a Ta-site vacancy amounts to 2.1 eV, almost a factor of two larger. Two contributions of distinct origin govern the temperature dependence of this asymmetry: a quasiharmonic contribution associated with the chemical-potential imbalance set by the two antisite structures and an anharmonic contribution associated with the local vibrational response of the vacancy structures. The asymmetry in anharmonic vibrations is traced back to an enlarged local vibrational distribution of the first-nearest neighbors around the Ta-site vacancy. In contrast to the vacancies, the Gibbs energies of antisite formation vary only weakly with temperature.
This work establishes key technological guidelines for designing diamond-based optoelectronic devices, derived from a first-principles investigation of two architectures: a PIN junction with a boron-vacancy-boron (BVB) intermediate-band absorber, and a PN junction based on phosphorus-vacancy (PV) defects. For the PIN solar cell, practical design principles include: i) aligning incident light in the xz-plane to exploit anisotropic absorption; ii) using graded junctions to mitigate tunnelling losses at abrupt interfaces; iii) targeting an absorber thickness of ~500 nm to balance absorption and carrier extraction; and iv) leveraging the high transparency of both contact layers for bifacial device configurations. For the PN diode, the PV-doped diamond operates via impurity-band conduction, making it suitable for degenerate p-type applications such as tunnel diodes or asymmetric junctions, while its temperature-dependent Seebeck anisotropy and sign-reversal offer opportunities for thermal management applications. When paired with phosphorus-doped n-type regions, these defects enable single-dopant junctions that significantly simplify device manufacturing. Using density functional theory with GW corrections, Bethe-Salpeter equation calculations and carrier transport modelling coupled to device electrostatics via a Poisson solver, we show that the BVB defect introduces intermediate bands without degrading diamond's high carrier mobility or thermal conductivity, while PV-doping provides high conductivity at room temperature through impurity-band transport. Overall, both defect-engineered systems preserve diamond's superior transport and thermal properties even after doping, offering viable pathways for high-performance diamond optoelectronics. These guidelines provide a practical foundation for fabricating efficient diamond-based photovoltaic and diode devices.
Isolated microscopic magnetic particles are used to induce local perturbations in dense colloidal suspensions by rotating an external magnet. Confocal microscopy enables tracking of both the magnetic probe particle and adjacent colloidal particles. A probe particle moves with a circular trajectory. Knowing the external force and measuring the amplitude and phase of the probe motion allows us to infer the storage and loss moduli of colloidal suspensions at various volume fractions. These measurements are in qualitative agreement with previous results from conventional rheology. To further analyze the system's response, the oscillatory amplitude of colloidal particles is evaluated as a function of distance from the probe, revealing a 1/r decay in amplitude, consistent with a homogeneous viscoelastic material. These observations confirm that continuum descriptions of the colloidal samples are effective down to length scales comparable to the particle diameter.
We propose a response tensor $\mathbf{\hat χ}$ to characterize the non-reciprocal critical current response of the superconducting (Josephson) diode effect. It describes the coupling between the dipole component of the angular distribution of the critical current and the applied magnetic field -- an analogue to the Hall response in the normal state. In quasi-2D systems with Rashba spin-orbit coupling and point group symmetries $C_{3v}$, $C_{4v}$ or $C_{6v}$, this tensor takes a fully antisymmetric form. When nematicity is present, a symmetric contribution emerges, providing an indicator of the nematic order in the superconducting state. In contrast, for systems exhibiting Dresselhaus spin-orbit coupling with the $D_{2d}$ symmetry, the tensor becomes diagonal traceless, and nematicity brings in a trace part. Our analysis not only accounts for the superconducting diode effect under external applied or intrinsic effective magnetic fields, but also predicts the symmetry conditions for realizing the diode effect when the magnetic field is aligned with the current. Beyond this, the proposed tensor provides a promising tool for detecting nematicity and potential nematic transitions deep within the superconducting phase. It may also encode additional information about the underlying electronic structure and symmetry-breaking orders, warranting further experimental investigation.
Defect centers in GaN emerge as bright sources of single-photons which recently have been demonstrated to optically interface a localized spin. However, the structure and composition of these defects as well as their efficient excitation techniques were not a subject of thorough studies. This work presents evidence that by tuning the excitation laser energy to specific resonance values the excitation efficiency can be enhanced, resulting in relative increase of photoluminescence intensity by up to an order of magnitude. The resonances can be selectively addressed with linearly polarized light, while the emission dipole remains unchanged, enabling polarization-controlled enhancement. These results establish an efficient way of excitation for GaN-based emitters, thereby increasing the generation rate of photons. The data is consistent with excitation via localized vibrational modes associated with point-defect complexes, establishing a practical quasi-resonant route to brighter, polarization-addressable operation of GaN defect emitters and clarifying their energy-level structure.
Periodic $π$-vacancies in graphene superlattices (GSLs) provide a symmetry-based route to band-gap opening in graphene by modifying the $π$-band dispersion. However, the symmetry conditions that determine whether a vacancy motif can open a band gap remain unclear. Here, we investigate periodic $π$-vacancy GSLs using a nearest-neighbor tight-binding model with one $p_z$ orbital per carbon site to identify the symmetry requirements for gap opening. $π$-vacancies, representing functionalized, substituted, or missing carbon sites, are modeled as site deletions in the $π$ basis, with all hopping matrix elements to and from the deleted sites set to zero. We focus on $π$-vacancy motifs with $C_2$ and $C_3$ point-group symmetry. A $3n \times 3n$ GSL, where $n=1,2,3,\ldots$ is the integer scaling factor multiplying the honeycomb primitive-cell vectors, folds $K$ and $K'$ to $Γ$ and can therefore open a band gap. For $C_3$-type vacancies, the Dirac cones remain pinned at high-symmetry points and thus stay at $Γ$ in folded $3n$ GSLs. In contrast, $C_2$-type vacancies that reduce the global point group of the GSL to $D_{2h}$ by preserving a pair of perpendicular mirror symmetries, $σ_v \perp σ_d$, can also constrain the Dirac cones to $Γ$. When the $σ_v$ and $σ_d$ mirror planes are absent, the cones are allowed to shift away from $Γ$ to $(\pm Δq,\pm Δq)$ in the $3n$ superlattice.
Efficient sub-10 nm electric transport remains a major challenge for nanoelectronics due to high losses and impedance mismatches in conventional Drude metals. Despite their promise of dissipationless, reflection-free conduction, topologically protected chiral edge modes remain little explored in their nanoscale spin polarized transport-particularly regarding real-space visualization, magnetic field tunability, and high-frequency edge conductivity. Here, we report magnetic brightening and nanoscale visualization of highly spin-polarizable infrared helical edge states using cryogenic magneto-infrared scattering-type scanning near-field optical microscopy (cm-IR-sSNOM). Our measurements reveal magnetic field-induced near-field conductivity at step edges, uncovering quantum spin Hall spin-splitting modes with enhanced infrared polarizability and slightly narrowed near-field profiles. In addition, the infrared edge electrodynamic response scales nearly linearly with atomic layer number, providing compelling evidence that magnetic-field-induced gaps do not disrupt individual-layer edge states at energies of around 100 meV. These results sharply contrast with microwave and DC transport, where even small magnetically induced gaps decrease edge conduction. Magnetically tunable, topologically robust high-frequency edge modes open a pathway toward ultralow-loss nanoscale interconnects and quantum logic architectures for next-generation microelectronics, spintronics and quantum information science.
Complex colloidal cluster morphologies are desirable for the fabrication of advanced materials, such as photonic crystals and meta-materials, and can be formed through evaporation-driven packing. By coupling lattice Boltzmann and discrete element methods, here we elucidate the rich interplay between fluid and particle dynamics during evaporation-driven self-assembly of spherical colloidal particles. We construct a regime diagram for a wide range of evaporation rates, interparticle friction coefficients, and particle numbers, identifying parameter regimes for open, closed, and minimal moment of inertia cluster configurations. Analyzing the competition between capillary, hydrodynamic, normal, and friction forces, we show that interparticle friction can exert a disproportionately strong influence on the final packing outcome despite being considerably smaller in magnitude than other forces at play. Our simulation results further highlight the potential for tuning colloidal cluster configurations via their dynamic trajectories.
We report on a two-step ultrahigh vacuum chemical vapor deposition synthesis of a vertical Ir(111)/borophene/hexagonal boron nitride heterostructure, using borazine as a single-source precursor. The process takes advantage of the finite solubility of boron in Ir: low precursor pressure at high temperature first establishes a boron reservoir in the near-surface region of the substrate, whereas subsequent growth at higher precursor pressure promotes the formation of a closed hexagonal boron nitride monolayer. During cooldown, the reduced boron solubility drives segregation to the surface, resulting in the formation of a borophene monolayer beneath the hexagonal boron nitride overlayer. The heterostructure, with micron sized grains, homogeneously covers the entire Ir substrate. The study is performed by complementary spot profile analysis low-energy electron diffraction, low-energy electron microscopy, and scanning tunneling microscopy measurements. This intrinsic segregation-assisted growth concept provides a promising route toward scalable synthesis of high-quality, vertical heterostructures of two-dimensional materials.