Papers

We study phonon-mediated pairings in a single-flavor metal with a tunable Berry curvature. In the absence of Berry curvature, we discover an unexpected possibility: $p$-wave superconductivity emerging purely from the retardation effect, while the static BCS approximation fails to predict its existence. The gap function exhibits sign-change behavior in frequency (owing to the dynamical structure of the phonon-mediated interaction in the $p$-wave channel), and $T_c$ obeys a BCS-like scaling. We further show that the Berry curvature stabilizes the chiral $p$-wave superconductivity and can induce transitions to higher-angular-momentum pairings. Our results establish that the phonon-mediated mechanism is a viable pairing candidate in single-flavor systems, such as the quarter-metal superconductivity observed in rhombohedral graphene multilayers.

We identify quantum geometric bounds for observables in non-Hermitian systems. We find unique bounds on non-Hermitian quantum geometric tensors, generalized two-point response correlators, conductivity tensors, and optical weights. We showcase these findings in topological systems with non-Hermitian Chern numbers. We demonstrate that the non-Hermitian geometric constraints on response functions naturally arise in open quantum systems governed by out-of-equilibrium Lindbladian dynamics. Our findings are relevant to experimental observables and responses under the realistic setups that fall beyond the idealized closed-system descriptions.

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 the recursion method for dynamical correlations and transport in fermionic systems on one-, two-, and three-dimensional lattices. The implementation is applicable in the strongly correlated regime, yields results directly in the thermodynamic limit, and covers all time scales, from short times through the thermalization time to the late-time asymptotics. Focusing on two paradigmatic models – interacting spinless fermions and the Hubbard model – we confirm the universal operator growth hypothesis in the fermionic case, compute infinite-temperature current-current autocorrelation functions, and determine the charge diffusion constant and the high-temperature conductivity. The diffusion constant is obtained in both the perturbative and the nonperturbative regime of the interaction strength, and we locate the boundary between them. We compare our approach with the Majorana propagation method, which converges up to intermediate times and agrees with ours where converged. Our results highlight a symbolic computational paradigm in which the most expensive step is performed once, producing reusable symbolic output that yields physical observables for arbitrary model parameters.

We present a diffuse neutron and x-ray scattering study of structural as well as 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 $\Gamma$-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 $\tau_0 \propto \epsilon^{-1.03 \pm 0.02}$ with the dissipation parameter $\epsilon = 1 - e^2$. (ii) The steady-state temperature follows a precise power-law $T_{\mathrm{steady}} \propto \epsilon^{-1.51 \pm 0.03}$, reflecting the non-linear balance between thermostat heating and collisional dissipation. (iii) The velocity autocorrelation function $\bar{A}(\tau_w, \tau)$ demonstrates pronounced aging, with decay rates $\lambda$ following a power-law slowing down $\lambda(\tau_w) \propto \tau_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 mixture's reduced third virial coefficient – which effectively captures three-body correlations and excluded-area constraints – we show that the RCP fraction depends nearly linearly on this parameter, leading to a near-universal collapse of simulation data over a wide range of size ratios and compositions. Comparisons with previous models by Brouwers and Zaccone indicate that the present approach provides more accurate and consistent predictions. The method can be naturally extended to polydisperse mixtures with continuous size distributions and is structurally consistent with the surplus equation-of-state formulation, offering a compact framework for understanding the near 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 $\nu=-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.

How can we learn the laws underlying the dynamics of stochastic systems when their trajectories are sampled sparsely in time? Existing methods either require temporally resolved high-frequency observations, or rely on geometric arguments that apply only to conservative systems, limiting the range of dynamics they can recover. Here, we present a new framework that reconciles these two perspectives by reformulating inference as a stochastic control problem. Our method uses geometry-driven path augmentation, guided by the geometry in the system's invariant density to reconstruct likely trajectories and infer the underlying dynamics without assuming specific parametric models. Applied to overdamped Langevin systems, our approach accurately recovers stochastic dynamics even from extremely undersampled data, outperforming existing methods in synthetic benchmarks. This work demonstrates the effectiveness of incorporating geometric inductive biases into stochastic system identification methods.

Photonic topological insulators (PTIs) offer robust platforms for light manipulation, but reconfigurable control of their topological properties without degrading performance remains a major challenge. While phase-change materials (PCMs) provide large refractive index modulation, widely used materials such as Ge2Sb2Te5 (GST) have been successfully deployed in commercial applications including optical data storage. However, they exhibit significant optical absorption in their crystalline state, which poses a challenge for transmissive photonic devices such as PTIs where high transparency is essential. Here, we overcome this fundamental limitation by integrating the ultra-low-loss PCM antimony triselenide (Sb2Se3) onto a silicon-based 2D PTI. We achieve submicron-scale selective patterning of Sb2Se3 on a photonic crystal for the first time, and demonstrate a topological phase transition induced by the material phase change. Owing to the transparency of Sb2Se3 in both its amorphous and crystalline states, a high Q-factor on the order of 10^3 is preserved-representing nearly an order-of-magnitude improvement over previous GST-based devices. This work resolves the absorption-loss bottleneck in reconfigurable PTIs and paves the way for practical, low-loss, tunable topological photonic devices.

We present a novel semiclassical framework tailored to determine the scaling dimensions of heavy neutral composite operators in conformal field theories (CFTs) which are inaccessible with other current methodologies. It utilizes the state-operator correspondence to map the desired scaling dimensions to the semiclassical energy spectrum of periodic homogeneous field configurations on a cylinder. As concrete applications, we provide detailed analyses for the $\phi^4$ theory near four dimensions and $\phi^6$ near three dimensions, semiclassically determining the full spectrum of neutral operators in the traceless symmetric Lorentz representations. Our methodology is presented pedagogically and is readily applicable to a vast class of CFTs.

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 {\deg}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 80%. The method demonstrated in this work opens up a new way to achieve fast, universal, and experiment-calibrated XANES prediction.

We define a black hole state in a spin chain by studying thermal correlators in holography. Focusing on the Heisenberg model we investigate the thermal and complexity properties of the black hole state by evaluating its entanglement entropy, emptiness formation probability and Krylov complexity. The entanglement entropy grows logarithmically with effective central charge c=5.2. We find evidence for thermalization at infinite temperature.

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 these settings realize "Swap" models exhibiting anomalous CSL phases, 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 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 Streda 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{\deg} 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.