arXiv Curator

The recent discovery of high-$T_c$ superconductivity in Ruddlesden-Popper (RP) nickelates has motivated extensive efforts to explore higher $T_c$ superconductors. Here, we systematically investigate Nd-doped La$_3$Ni$_2$O$_7$ using density functional theory (DFT) and renormalized mean-field theory (RMFT). DFT calculations reveal that both the lattice constants and interlayer spacing decrease upon Nd substitution, similar to the effect of physical pressure. However, the in-plane Ni-O-Ni bond angle evolves non-monotonically with doping, increasing to a maximum at 70% (~2/3) Nd doping level and then falling sharply at 80%, which leads to a reduction in orbital overlap. Moreover, Nd doping has a more pronounced effect on the Ni-$d{_{z^2}}$ orbital, demonstrating an orbital-dependent effect of rare-earth substitution. Through the bilayer two-orbital t-J model, RMFT analysis further shows an $s\pm$-wave pairing symmetry, with $T_c$ rising to a maximum at about 70% Nd substitution before declining, in agreement with the transport measurements. The variation in $T_c$ can be traced to the competition between continuously enhanced interlayer superexchange coupling $J_\perp^z$ and a gradual decrease in particle density. These results highlight the delicate interplay among structural tuning, orbital hybridization, and superconductivity, providing important clues to design higher-$T_c$ RP nickelate superconductors.

The recently discovered superconducting infinite-layer nickelates offer a novel platform to explore an exotic pairing mechanism in multi-band systems towards high-temperature superconductivity and associated rich quantum phases, contrasting with cuprates. Here, we show that infinite-layer (La,Sr)NiO$_2$ exhibits strong-coupling superconductivity, resilient to in-plane magnetic fields exceeding 47 T at optimal doping - more than twice the Pauli limit for conventional BCS superconductors. This violation becomes pronounced towards the underdoped regime, implying an intriguing evolution of pairing glue. The unexpected observation of positive Seebeck coefficients in this regime indicates the presence of nontrivial electron correlations. Furthermore, our comprehensive investigation across the superconducting dome reveals that the evolution of (thermo)electric normal-state properties - specifically, the sign changes of the Hall and Seebeck coefficients - coincide with the evolution of superconducting anisotropy and pairing strength. This demonstrates a definitive link between electron correlations and strong-coupling superconductivity in (La,Sr)NiO$_2$, contributing to a unified framework for understanding unconventional superconductivity.

Motivated by the discovery of superconductivity in bilayer La$_3$Ni$_2$O$_7$ at 80 K and the increased superconducting transition temperature, $T_\text{c}$, up to 92 K in single crystals of La$_2$SmNi$_2$O$_7$ under pressure, we systematically study the effect of Sm doping on the superconductivity and structure of La$_{3-x}$Sm$_x$Ni$_2$O$_7$ (0 $\leq$ x $\leq$ 1.5) under pressure. Experimental investigations in polycrystalline samples reveal that Sm doping monotonically decreases the lattice constants $c$ and $a$, thereby enhancing crystal structure distortion and leading to an evolution of the metallic ground state in La$_3$Ni$_2$O$_7$ to an insulating state in La$_{1.5}$Sm$_{1.5}$Ni$_2$O$_7$. The maximum onset $T_\text{c}$ in compounds $x=0.9$ and 1.5 is 89 K, while the pressure that drives the emergence of superconductivity is higher for higher doping levels. The results suggest that the enhancement of $T_\text{c}$ in La$_{3-x}$Sm$_x$Ni$_2$O$_7$ is mainly affected by the compressed $c$ lattice before saturation, and the structure transition is critical for the emergence of superconductivity. Our experimental results provide insight into the influence of elemental substitution on nickelate superconductors, offering a means to increase the transition temperature further.

The discovery of high-temperature superconductivity (HTSC) in La$_3$Ni$_2$O$_7$ has aroused significant interest in exploring the pairing mechanism. Previous studies have proposed an effective d$_{x^2-y^2}$-orbital bilayer $t-J_{\parallel}-J_{\perp}$ model, in which the electrons of the d$_{x^2-y^2}$ orbital are charge carriers, which are subject to the intralayer antiferromagnetic (AFM) superexchange $J_{\parallel}$ and the large interlayer AFM superexchange $J_{\perp}\approx 2J_{\parallel}$, with the latter transferred from the nearly half-filled and hence localized $d_{z^2}$ orbital through the strong Hund's rule coupling. Here we study this model by the variational Monte Carlo (VMC) simulation and find a dominant interlayer s-wave pairing, in which the SC order parameters have a drastic improvement compared with those of the mean field (MF) type of theories. In real materials, the Hund's coupling is finite, leading to reduced $J_{\perp}$, dictating that the MF-type theories have difficulty explaining the HTSC. However, our VMC calculations find that even for effective $J_{\perp}$ as weak as $J_{\perp}=J_{\parallel}$, the interlayer pairing is still considerably large and can be compared with the $T_c$ observed in experiments, which is very weak in MF-type theories. This result indicates the important role of the Gutzwiller projection in improving the $T_c$, which is ignored in the MF-type theories. In addition, our results show that suppressed interlayer hopping can promote interlayer pairing, which is consistent with the fact that the interlayer hopping of the d$_{x^2-y^2}$ orbital in La$_3$Ni$_2$O$_7$ is very weak. Our research offers a new perspective for understanding the pairing mechanism of bilayer nickelates and provides a reference for recent ultra-cold atom experiments in mixed-dimensional systems.

The family of nickelate superconductors have long been explored as analogs of the high temperature cuprates. Nonetheless, the recent discovery that certain stoichiometric nickelates superconduct up to high $T_c$ under pressure came as a surprise. The mechanisms underlying the superconducting state remain experimentally unclear. In addition to the practical challenges posed by working in a high pressure environment, typical samples exhibit anomalously weak diamagnetic responses, which have been conjectured to reflect inhomogeneous `filamentary' superconducting states. We perform wide-field, high-pressure, optically detected magnetic resonance spectroscopy to image the local diamagnetic responses of as grown La$_3$Ni$_2$O$_7$ samples \emph{in situ}, using nitrogen vacancy quantum sensors embedded in the diamond anvil cell. These maps confirm significant inhomogeneity of the functional superconducting responses at the few micron scale. By spatially correlating the diamagnetic Meissner response with both the local tensorial stress environment, also imaged \emph{in situ}, and stoichiometric composition, we unravel the dominant mechanisms suppressing and enhancing superconductivity. Our wide-field technique simultaneously provides a broad view of sample behavior and excellent local sensitivity, enabling the rapid construction of multi-parameter phase diagrams from the local structure-function correlations observed at the sub-micron pixel scale.

This paper presents the temperature dependence of the static spin susceptibility at $\mathbf{q} = (\pi, \pi)$ and $\mathbf{q} = (\pi, 0)$ in a Mott-proximate cuprate model with an antinodal pseudogap -- a model system for high-temperature superconducting (HTSC) cuprates. The results show the susceptibility onset temperature tracks the critical temperature ($T_c$) of HTSCs with a comparable scale across the electron filling factor. Also, as the electron filling decreases and the chemical potential approaches the antinodal van Hove region, the susceptibility at $\mathbf{q}=(\pi,0)$ -- the axial particle-hole response -- grows markedly. The emergence of cuprate superconductivity correlates with a suppression of low-energy antinodal spin response and associated particle-hole excitations, which would otherwise dephase $d$-wave pairing, commonly attributed to spin fluctuations. In this context, the pseudogap partially suppresses antinodal spectral weight near $\omega = 0$, thereby reducing the low-$\omega$ particle-hole phase space.

This paper presents the temperature dependence of the static spin susceptibility at $\mathbf{q} = (\pi, \pi)$ and $\mathbf{q} = (\pi, 0)$ in a Mott-proximate cuprate model with an antinodal pseudogap -- a model system for high-temperature superconducting (HTSC) cuprates. The results show the susceptibility onset temperature tracks the critical temperature ($T_c$) of HTSCs with a comparable scale across the electron filling factor. Also, as the electron filling decreases and the chemical potential approaches the antinodal van Hove region, the susceptibility at $\mathbf{q}=(\pi,0)$ -- the axial particle-hole response -- grows markedly. The emergence of cuprate superconductivity correlates with a suppression of low-energy antinodal spin response and associated particle-hole excitations, which would otherwise dephase $d$-wave pairing, commonly attributed to spin fluctuations. In this context, the pseudogap partially suppresses antinodal spectral weight near $\omega = 0$, thereby reducing the low-$\omega$ particle-hole phase space.

In many unconventional superconductors, the strange-metal regime is thought to emerge from quantum criticality, yet in cuprates this link is obscured by the enigmatic pseudogap. Superconducting infinite-layer nickelates provide a new arena to test this paradigm but are constrained to thin films, precluding calorimetry. We use the Seebeck coefficient as a low-temperature proxy for entropy per carrier and uncover a clear quantum-critical thermodynamic signature: in La$_{1-x}$Sr$_x$NiO$_2$ at the onset of $T$-linear resistivity ($x=0.20$), $S/T$ diverges logarithmically upon cooling, $S/T \propto \log T$. Boltzmann transport based on ARPES-derived band structure reproduces the high-temperature magnitude and sign of $S/T$ and reveals a threefold mass renormalization at the Fermi level. To identify the terminating phase, we analyze Hall data across Nd$_{1-x}$Sr$_x$NiO$_2$ and show that its temperature evolution is quantitatively captured by a minimal two-band model in which a strongly correlated Ni-$d_{x^2-y^2}$ Fermi surface exhibits Planckian $T$-linear scattering while the rare-earth Nd-$s$ pocket remains Fermi-liquid-like. Inverting the zero-temperature Hall response reveals a collapse of the Ni-$d_{x^2-y^2}$ band carrier density from $1+p$ to $p$ holes across the critical doping, without long-range magnetic order -- mirroring the cuprate pseudogap transition in cuprates. These results establish a quantum critical point at the end of a pseudogap-like phase in infinite-layer nickelates and unify the broader paradigm among correlated superconductors that strange metal behaviour is intimately linked to quantum criticality.

The remarkable discovery of high temperature superconductivity in bulk bilayer nickelates under high pressure has prompted the conjecture that epitaxial compressive strain might mimic essential aspects of hydrostatic pressure. The successful realization of superconductivity in films on SrLaAlO4 (001) (SLAO) supports this correspondence, yet it remains unclear whether the rich pressure-temperature phase diagram of bilayer nickelates can be systematically mapped (and studied at ambient pressure) as a function of epitaxial strain. To this end, experimental access near the elusive edge of the superconducting phase boundary would provide invaluable insight into the nature of the superconducting state and the ground state from which it emerges. It would also offer a benchmark for theoretical models. Here we report superconducting bilayer nickelates grown on LaAlO3 (001) (LAO), where the compressive strain required for ambient-pressure superconductivity is nearly halved to -1.2%. These films exhibit a superconducting onset above 10 K and reach zero resistance at 3 K, with normal-state transport properties differing from those of films grown on SLAO. Our results offer a new opportunity to probe emergent phenomena near the superconducting phase boundary in the strain-temperature phase diagram of bilayer nickelates.

The pairing symmetry of the recently discovered bilayer nickelate superconductor La$_3$Ni$_2$O$_7$ is a subject of intense debate in condensed matter physics, with the two leading theoretical candidates being a sign-reversing $s_{\pm}$-wave and a $d$-wave state. To investigate its ground-state properties in the intermediate coupling regime which is critical for real materials, we construct a two-orbital bilayer Hubbard model and employ the constrained-path quantum Monte Carlo method for large-scale simulations. By systematically calculating ground-state pairing correlation functions across parameter spaces, we map its pairing symmetry phase diagram. We find that an increasing Hund's coupling selectively enhances the interlayer $s_{\pm}$-wave pairing while suppressing the intralayer $d$-wave pairing. Similarly, a larger crystal field splitting drives a transition from $d$-wave- to $s_{\pm}$-wave-dominant states. Further analysis reveals that the strength of the intralayer $d$-wave pairing is strongly correlated with the $(\pi, \pi)$ antiferromagnetic spin fluctuations, which are in turn effectively suppressed by a large crystal field splitting, thereby weakening the $d$-wave pairing channel. Additionally, the dominant pairing symmetry transition region roughly overlaps with the inversion of orbital occupancy response to Hubbard $U$, suggesting an intrinsic link between pairing competition and orbital physics. Our results indicate that, within the parameter regime relevant to the actual material, the $s_{\pm}$-wave is the most probable pairing symmetry.

Loop currents, as an orbital magnetism, have been proposed as a possible fluctuation mechanism for superconducting pairing, which always remains elusive. Here, we investigate the role of an interlayer loop current fluctuation in mediating superconductivity using an unbiased bilayer $t-J_{\perp}-V$ model via sign-problem-free projector quantum Monte Carlo simulations. The model spontaneously generates the interlayer loop current by breaking time-reversal and translational symmetries, favored by interlayer Coulomb repusion. With hole doping, the loop current is rapidly suppressed, while its fluctuations give rise to an interlayer $s$-wave superconductivity. Our results establish a phase diagram to demonstrate a superconducting transition due to the interlayer loop current fluctuations. It also provides possible insights into some physics related to bilayer nickelates, with which it shares a similar structure and a large interlayer spin exchange.

Phase fluctuations are a key factor distinguishing nonthermal (ultrafast) and thermal phase transitions. Charge order in cuprates is characterized by short-range coherence while competing with superconductivity, and as such, it provides a representative case to study the role of phase fluctuation in coupled order parameter dynamics. In this work, we investigated the intertwined evolution of charge order and superconductivity in cuprate/manganite heterostructures using time-resolved resonant X-ray scattering. The resulting dynamics are analyzed within a space- and time-dependent nonperturbative model capturing both amplitude and phase dynamics. At low fluence, photo-induced suppression of superconductivity results in a nonthermal enhancement of charge order, underscoring the dynamic competition between charge order and superconductivity. With increasing fluence, the slowing down of melting and recovery dynamics is observed, indicating a critical role of phase fluctuations. At high fluence, both charge order and superconductivity remain suppressed for an extended time window due to decoupling between amplitude and phase dynamics and the delayed recovery of phase coherence. Our work underscores the importance of phase fluctuation for understanding the dynamic competition between order parameters in cuprates.

Simulation of the time-dynamics of fermionic many-body systems has long been predicted to be one of the key applications of quantum computers. Such simulations -- for which classical methods are often inaccurate -- are critical to advancing our knowledge and understanding of quantum chemistry and materials, underpinning a wide range of fields, from biochemistry to clean-energy technologies and chemical synthesis. However, the performance of all previous digital quantum simulations has been matched by classical methods, and it has thus far remained unclear whether near-term, intermediate-scale quantum hardware could offer any computational advantage in this area. Here, we implement an efficient quantum simulation algorithm on Quantinuum's System Model H2 trapped-ion quantum computer for the time dynamics of a 56-qubit system that is too complex for exact classical simulation. We focus on the periodic spinful 2D Fermi-Hubbard model and present evidence of spin-charge separation, where the elementary electron's charge and spin decouple. In the limited cases where ground truth is available through exact classical simulation, we find that it agrees with the results we obtain from the quantum device. Employing long-range Wilson operators to study deconfinement of the effective gauge field between spinons and the effective potential between charge carriers, we find behaviour that differs from predictions made by classical tensor network methods. Our results herald the use of quantum computing for simulating strongly correlated electronic systems beyond the capacity of classical computing.

Recent angle-dependent magnetoresistance measurements in underdoped cuprates have revealed compelling evidence for small hole pockets in the pseudogap regime, including observation of the Yamaji effect in HgBa$_2$CuO$_{4+\delta}$ (Chan et al., Nature Physics 10.1038/s41567-025-03032-2 (2025)). A key distinction between theories is their predicted Fermi volumes, measured as fractions of the square lattice Brillouin zone: $p/4$ per pocket for spin density wave (SDW) versus $p/8$ for fractionalized Fermi liquid (FL*), where $p$ is the hole doping. We calculate the $c$-axis magnetoresistance $\rho_{zz}(\theta, \phi)$ within the semiclassical Boltzmann formalism for both states, and using the ancilla layer model (ALM) for FL* in a single-band Hamiltonian. While the SDW Yamaji azimuthal angle $\theta$ is marginally smaller than the FL* Yamaji $\theta$, the most significant difference is the presence of a second Yamaji peak in the SDW theory near in-plane angle $\phi = 45^\circ$; there is no such second peak in the observations in HgBa$_2$CuO$_{4+\delta}$. Furthermore, we examine the role of anisotropic interlayer hopping present in LSCO, and predict the potential emergence of the Yamaji effect in LSCO at higher magnetic fields. Our results support the FL* interpretation of Fermi arcs in the pseudogap phase, and establish Yamaji angle measurements as a discriminatory tool between theoretical models.

The Anderson impurity model (AIM) is of fundamental importance in condensed matter physics to study strongly correlated effects. However, accurately solving its long-time dynamics still remains a great numerical challenge. An emergent and rapidly developing numerical strategy to solve the AIM is to represent the Feynman-Vernon influence functional (IF), which encodes all the bath effects on the impurity dynamics, as a matrix product state (MPS) in the temporal domain. The computational cost of this strategy is basically determined by the bond dimension $\chi$ of the temporal MPS. In this work, we propose an efficient and accurate method which, when the hybridization function in the IF can be approximated as the summation of $n$ exponential functions, can systematically build the IF as a MPS by multiplying $O(n)$ small MPSs, each with bond dimension $2$. Our method gives a worst case scaling of $\chi$ as $2^{8n}$ and $2^{2n}$ for real- and imaginary-time evolution respectively. We demonstrate the performance of our method for two commonly used bath spectral functions, where we show that the actually required $\chi$s are much smaller than the worst case.

We perform numerically exact determinant quantum Monte Carlo simulations of the Hubbard model and analyze pairing tendencies by evaluating correlation functions at the imaginary-time midpoint ($\tau=\beta/2$), which suppresses high-frequency weight and emphasizes low-energy physics. Using this diagnostic, we identify clear finite-temperature signatures of underlying $d$-wave superconductivity for electron doping, while finding no clear indication upon cooling for hole doping. Our analysis enables direct comparison with ground-state DMRG, revealing consistent real-space pairing patterns. These results provide a practical route to bridge the gap between finite-temperature and ground-state numerically exact simulations of the Hubbard model despite the fermion sign problem.

Artificial Intelligence (AI) has become an exceptionally powerful tool for analyzing scientific data. In particular, attention-based architectures have demonstrated a remarkable capability to capture complex correlations and to furnish interpretable insights into latent, otherwise inconspicuous patterns. This progress motivates the application of AI techniques to the analysis of strongly correlated electrons, which remain notoriously challenging to study using conventional theoretical approaches. Here, we propose novel AI workflows for analyzing snapshot datasets from tensor-network simulations of the two-dimensional (2D) Hubbard model over a broad range of temperature and doping. The 2D Hubbard model is an archetypal strongly correlated system, hosting diverse intriguing phenomena including Mott insulators, anomalous metals, and high-$T_c$ superconductivity. Our AI techniques yield fresh perspectives on the intricate quantum correlations underpinning these phenomena and facilitate universal omnimetry for ultracold-atom simulations of the corresponding strongly correlated systems.

The two-dimensional (2D) Hubbard model has long attracted interest for its rich phase diagram and its relevance to high-$T_c$ superconductivity. However, reliable finite-temperature studies remain challenging due to the exponential complexity of many-body interactions. Here, we introduce an enhanced $1\text{s}^+$ eXponential Tensor Renormalization Group algorithm that enables efficient finite-temperature simulations of the 2D Hubbard model. By exploring an expanded space, our approach achieves two-site update accuracy at the computational cost of a one-site update, and delivers up to 50% acceleration for Hubbard-like systems, which enables simulations down to $T\!\approx\!0.004t$. This advance permits a direct investigation of superconducting order over a wide temperature range and facilitates a comparison with zero-temperature infinite Projected Entangled Pair State simulations. Finally, we compile a comprehensive dataset of snapshots spanning the relevant region of the phase diagram, providing a valuable reference for Artificial Intelligence-driven analyses of the Hubbard model and a comparison with cold-atom experiments.

Competing orders represent a central challenge in understanding strongly correlated systems. In this work, we employ projector quantum Monte Carlo simulations to study a sign-problem-free bilayer extended Hubbard model. In this model, a charge stripe phase, characterized by a peak at momentum $k_x=2\pi\delta$ is induced by highly anisotropic interlayer spin-exchange coupling $J_z$, and strongly suppressed upon introducing the spin-flip term $J_\bot$; in contrast, \(J_\perp\) favors the emergence of interlayer pairing superconductivity. We further demonstrate that the anisotropy of the interlayer spin-exchange directly governs the competition between these two phases, while the on-site interaction \(U\) plays a complex role in tuning both the charge stripe and superconductivity. Our work identifies the key factors driving charge stripe formation, highlights the sensitivity of both the charge stripe and superconducting phases to interaction parameters, and thereby provides valuable insights into competing orders in strongly correlated systems.

This article explores the deep interconnections among three seemingly unrelated concepts in condensed matter physics: electronic liquid crystal phases, multipole expansions, and altermagnetism. At the heart of these phenomena lies a shared foundation: spin-momentum locking in the nonrelativistic regime. Originally proposed in the context of electronic liquid crystal phases, spin-momentum locking was later elegantly incorporated into the formalism of multipole expansions. This framework can be further extended across multiple atomic sites, making it particularly effective for describing altermagnets, which feature localized magnetic moments distributed over at least two sublattices. In the second part of the article, we examine superconducting phenomena that stem from this shared mechanism, focusing on superconductivity in systems with spin-momentum locked Fermi surfaces. We highlight a rich variety of unconventional superconducting states, including finite-momentum pairing, $d$-wave and spin-triplet superconductivity, and topological Bogoliubov Fermi surfaces, among others. Additional related topics are addressed in the concluding section. Overall, this work offers both an accessible introduction to the newly identified magnetic order known as altermagnetism and a conceptual guide for researchers aiming to harness the ensuing unconventional superconductivity in the development of future quantum technologies.

Simulating the time-dynamics of quantum many-body systems was the original use of quantum computers proposed by Feynman, motivated by the critical role of quantum interactions between electrons in the properties of materials and molecules. Accurately simulating such systems remains one of the most promising applications of general-purpose digital quantum computers, in which all the parameters of the model can be programmed and any desired physical quantity output. However, performing such simulations on today's quantum computers at a scale beyond the reach of classical methods requires advances in the efficiency of simulation algorithms and error mitigation techniques. Here we demonstrate programmable digital quantum simulation of the dynamics of the 2D Fermi-Hubbard model -- one of the best-known simplified models of electrons in crystalline solids -- at a scale beyond exact classical simulation. We implement simulations of this model on lattice sizes up to $6\times 6$ using 72 qubits on Google's Willow quantum processor, across a range of physical parameters, including on-site electron-electron interaction strength and magnetic flux, and study phenomena including formation of magnetic polarons, i.e. charge carriers surrounded by local magnetic polarisation, dynamical symmetry breaking in stripe-ordered states, attraction of charge carriers on an entangled state known as a valence bond solid, and the approach to equilibrium through thermalisation. We validate our results against exact calculations in parameter regimes where these are feasible, and compare them to approximate classical simulations performed using tensor network and operator propagation methods. Our results demonstrate that programmable digital quantum simulation of many-body interacting electron models is now competitive on state-of-the-art quantum hardware.

We develop a nonequilibrium multi-orbital extension of the two-particle self-consistent theory and apply it to the bilayer Hubbard model as a minimal platform to investigate correlation effects in the presence of interlayer interactions and tunneling. The method determines vertex corrections in the spin and charge channels self-consistently at the two-particle level, thereby avoiding the spurious finite-temperature phase transitions that limit dynamical mean-field theory in two dimensions. We derive the spectral self-energy and implement the framework directly on the real-frequency axis within the Keldysh nonequilibrium Green's function formalism, enabling the treatment of both equilibrium and non-equilibrium steady states without relying on numerical analytic continuation. As an application, we demonstrate that a pseudogap can emerge in the bilayer Hubbard model when spin, charge, or excitonic fluctuations become sufficiently strong. Instabilities in different channels are also evaluated in an unbiased manner across the parameter space. Remarkably, we find that the excitonic susceptibility grows with increasing interlayer bias, before it gets suppressed at large biases by the charge imbalance between the layers. This work establishes a versatile and computationally efficient framework for investigating correlated multi-orbital systems under nonequilibrium conditions.

We report on the study of the magnetic excitations of Mott insulator La2CuO4 by using resonant inelastic x-ray scattering (RIXS) and cluster calculations within the framework of exact diagonalization. Our results demonstrate experimentally that the bimagnon excitation in Cu L-edge RIXS is enhanced if the incident x-ray energy is slightly above the absorption edge. Through incident-energy-dependent momentum-resolved RIXS, we investigated the excitation of the bimagnon with predominantly A1 symmetry. The bimagnons of La2CuO4 exhibit a nearly flat dispersion with momentum along the Cu-O bond direction. This observation agrees with the bimagnon dispersion from the calculations on a single-band Hubbard model rather than a Heisenberg model with only the nearest neighbor exchange interaction. This means that the effect of the higher-order spin couplings such as the cyclic or ring exchange interactions caused by the coherent motion of electrons beyond nearest-neighbor sites is important for understanding the bimagnon dynamics of cuprates.

Strong evidence of unconventional superconductivity has been very recently reported experimentally in twisted transition metal dichalcogenide bilayer and gathered a significant amount of interest. Here we consider the Hubbard model on a triangular lattice describing the hole-doped moir\'e superlattice emerging in WSe$_{2}$/WSe$_{2}$ twisted homobilayer in the moderately correlated regime. By applying the Density Matrix Renormalization Group, we diagonalize the spin-valley-polarized Hamiltonian and show signatures of coexisting singlet and triplet pairings in the range of hole dopings and displacement fields reported in the experiments. In this view, we show that the superconductivity in the WSe$_{2}$/WSe$_{2}$ twisted homobilayer is likely to be induced by electronic correlations and has a mixed-symmetry character. These predictions can shed light on the nature of the superconducting state observed in the twisted homobilayer of WSe$_{2}$/WSe$_{2}$. We also identify the emerging superconducting orders, which are $d_{xy}(d_{x^2-y^2} \pm id_{xy} )$ and $p_y(p_{x}\mp ip_{y})$ for the singlet and triplet channels in the cylinder of width three(four), respectively.

Doped Mott insulators host intertwined spin-charge phenomena that evolve with temperature and can culminate in stripe order or superconductivity at low temperatures. The two-dimensional $t$-$J$ model captures this interplay yet finite-temperature, infinite-size calculations remain difficult. Using purification represented by a tensor network - an infinite projected entangled-pair state (iPEPS) ansatz - we simulate the $t$-$J$ model at finite temperature directly in the thermodynamic limit, reaching temperatures down to one tenth of the hopping rate and hole concentrations up to one quarter of the lattice sites. Beyond specific heat, uniform susceptibility, and compressibility, we introduce dopant-conditioned multi-point correlators that map how holes reshape local exchange. Nearest-neighbor hole pairs produce a strong cooperative response that reinforces antiferromagnetism on the adjacent parallel bonds, and single holes weaken nearby antiferromagnetic bonds; d-wave pairing correlations remain short-ranged over the same window. These results provide experiment-compatible thermodynamic-limit benchmarks and establish dopant-conditioned correlators as incisive probes of short-range spin-texture reorganization at finite temperature.

Doped Mott insulators host intertwined spin-charge phenomena that evolve with temperature and can culminate in stripe order or superconductivity at low temperatures. The two-dimensional $t$-$J$ model captures this interplay yet finite-temperature, infinite-size calculations remain difficult. Using purification represented by a tensor network - an infinite projected entangled-pair state (iPEPS) ansatz - we simulate the $t$-$J$ model at finite temperature directly in the thermodynamic limit, reaching temperatures down to one tenth of the hopping rate and hole concentrations up to one quarter of the lattice sites. Beyond specific heat, uniform susceptibility, and compressibility, we introduce dopant-conditioned multi-point correlators that map how holes reshape local exchange. Nearest-neighbor hole pairs produce a strong cooperative response that reinforces antiferromagnetism on the adjacent parallel bonds, and single holes weaken nearby antiferromagnetic bonds; d-wave pairing correlations remain short-ranged over the same window. These results provide experiment-compatible thermodynamic-limit benchmarks and establish dopant-conditioned correlators as incisive probes of short-range spin-texture reorganization at finite temperature.

We summarize key experimental studies of the low energy properties of overdoped cuprate high temperature superconductors and conclude that a theoretical understanding of the "essential physics" is achievable in terms of a conventional Fermi-liquid treatment of the normal state, and a BCS mean-field treatment of the (d-wave) superconducting state. For this perspective to be consistent, it is necessary to posit that there is a crossover from a strongly correlated underdoped regime (where a different theoretical perspective is necessary) to the more weakly correlated overdoped regime. It is also necessary to argue that the various observed features of the overdoped materials that are inconsistent with this perspective can be attributed to the expected effects of the intrinsic disorder associated with most of the materials being solid state solutions (alloys). As a test of this idea, we make a series of falsifiable predictions concerning the expected behavior of "ideal" (disorder free) overdoped cuprates.

We report the observation of an emergent charge order mode in the optically-excited cuprate ladder Sr$_{14}$Cu$_{24}$O$_{41}$. Near-infrared light in the ladder plane drives a symmetry-protected electronic metastable state together with a partial melting of the equilibrium charge order. Our time-resolved resonant inelastic x-ray scattering measurements at the upper Hubbard band reveal a gapless collective excitation dispersing from the charge-order wavevector up to 0.8 eV with a slope on the order of the quasiparticle velocity. These findings reveal a regime where correlated carriers acquire itinerant character at finite momentum, and charge order becomes dynamically fluctuating, offering a platform to explore light-induced pairing instabilities.

Controlling nonequilibrium dynamics in quantum materials requires ultrafast probes with spectral selectivity. We report femtosecond reflectivity measurements on the cuprate superconductor Bi$_2$Sr$_2$CaCu$_2$O$_{8+\delta}$ using free-electron laser extreme-ultraviolet (23.5--177~eV) and near-infrared (1.5~eV) pump pulses. EUV pulses access deep electronic states, while NIR light excites valence-band transitions. Despite these distinct channels, both schemes produce nearly identical dynamics: above $T_c$, excitations relax through fast (100--300~fs) and slower (1--5~ps) channels; below $T_c$, a delayed component signals quasiparticle recombination and condensate recovery. We find that when electronic excitations are involved, the ultrafast response is governed mainly by absorbed energy rather than by the microscopic nature of the excitation. In contrast, bosonic driving in the THz or mid-infrared produces qualitatively different dynamics. By demonstrating that EUV excitation of a correlated superconductor yields macroscopic dynamics converging with those from optical pumping, this work defines a new experimental paradigm: FEL pulses at core-level energies provide a powerful means to probe and control nonequilibrium electronic states in quantum materials on their intrinsic femtosecond timescales. This establishes FEL-based EUV pumping as a new capability for ultrafast materials science, opening routes toward soft X-ray and attosecond studies of correlated dynamics.

We developed a practical quantum advantage benchmarking framework that connects the accumulation of entropy in a quantum processing unit and the degradation of the solution to a target optimization problem. The benchmark is based on approximating from below the Gibbs states boundary in the energy-entropy space for the application of interest. We believe the proposed benchmarking technique creates a powerful bridge between hardware benchmarking and application benchmarking, while remaining hardware-agnostic. It can be extended to fault-tolerant scenarios and relies on computationally tractable numerics. We demonstrate its applicability on the problem of finding the ground state of the two-dimensional Fermi-Hubbard.

Engineering van Hove singularities (vHss) near the Fermi level, if feasible, offers a powerful route to control exotic quantum phases in electronic and magnetic behaviors. However, conventional approaches, which rely primarily on chemical and electrical doping, focus mainly on local electrical or optical measurements, limiting their applicability to coupled functionalities. In this study, a vHs-induced insulator-metal transition coupled with a ferromagnetic phase transition was empirically achieved in atomically designed quasi-2D SrRuO3 (SRO) superlattices via epitaxial strain engineering, which has not been observed in conventional 3D SRO systems. Theoretical calculations revealed that epitaxial strain effectively modulates the strength and energy positions of vHs of specific Ru orbitals, driving correlated phase transitions in the electronic and magnetic ground states. X-ray absorption spectroscopy confirmed the anisotropic electronic structure of quasi-2D SRO modulated by epitaxial strain. Magneto-optic Kerr effect and electrical transport measurements demonstrated modulated magnetic and electronic phases. Furthermore, magneto-electrical measurements detected significant anomalous Hall effect signals and ferromagnetic magnetoresistance, indicating the presence of magnetically coupled charge carriers in the 2D metallic regime. This study establishes strain engineering as a promising platform for tuning vHss and resultant itinerant ferromagnetism of low-dimensional correlated quantum systems.

We present a tangent-space Krylov (TaSK) method for efficient computation of zero-temperature real-frequency spectral functions on top of ground state (GS) matrix product states (MPS) obtained from the Density Matrix Renormalization Group. It relies on projecting resolvents to the tangent space of the GS-MPS, where they can be efficiently represented using Krylov space techniques. This allows for a direct computation of spectral weights and their corresponding positions on the real-frequency axis. We demonstrate the accuracy and efficiency of the TaSK approach by showcasing spectral data for various models. These include the 1D Haldane-Shastry and Heisenberg models as benchmarks. As an interesting application, we study the Hubbard model on a cylinder at half-filling, augmented by a density-assisted hopping (DAH) term. We find that DAH leads to particle-hole asymmetric single-particle mobilities and lifetimes in the resulting Mott insulator, and identify the responsible scattering processes. Further, we find that DAH influences the dispersion of Green's function zeros beyond its range, which has a frustrating effect on the Mott insulator studied here.

SrCuO2 has long been considered a near-archetypal realization of a quasi one dimensional (1D) system of interacting electrons with short-range interactions. Within this framework, experimental observations - interpreted through the lens of the 1D Hubbard model-suggest that electron and hole excitations decay into two types of (unphysical) collective bosonic modes: spinons, which carry the spin degree of freedom, and holons, which carry the charge degree of freedom. This model, known as spin-charge separation, is most directly evidenced by angle-resolved photoemission spectroscopy (ARPES), where a photo-induced hole decays into a continuum of these excitations. Here we present an alternative perspective grounded in first-principles, self-consistent, and parameter-free many-body perturbation theory. In this revised quasiparticle framework, ARPES can be understood as a one-body effect arising from mild disorder in a long-range antiferromagnetic ground state. the emergence of the so-called spinon branch arises naturally from spin disorder, the anomalous line widths are accurately captured, and we provide a compelling explanation for the spectral weight observed at the non-magnetic zone boundary. This reinterpretation provides a unified explanation for key experimental signatures previously attributed to spin-charge separation, including features observed in optical conductivity. Additionally, we show that SrCuO2 exhibits a nontrivial interchain coupling that significantly influences both its one-particle and two-particle spectral functions. By comparing the spectral features of SrCuO2 with those of La2CuO4, we argue that SrCuO2 shares notable similarities with the two-dimensional cuprates - both being rooted in a common CuO4 plaquette-based molecular orbital framework.

The fate of spin-charge separation beyond the low energy remains elusive up to now. Here we develop a microscopic theory of the correlation functions using the strong coupling expansion of the Hubbard model and demonstrate its validity down to the experimentally relevant $r_{\rm s}>1$. Evaluating the spectral function, we show the general stability of the nonlinear spin-charge modes in whole energy band and investigate all the nonlinear features systematically. We confirm the general prediction experimentally in semiconductor quantum wires. Furthermore, we observe a signal consistent with a continuum of the nonlinear excitations and with a final spectral density around the $3 k_{\rm F}$ point, indicating the robustness of the Hubbard model predictions for a finite range interaction.

Understanding non-equilibrium spin transport through 2D ferromagnets is a theoretical challenge, as correlations produce a complex electronic structure with coexisting itinerant and localized electrons. We have developed a fully non-equilibrium ab initio method, combining density functional theory, dynamical mean-field theory, and non-equilibrium Green's functions to investigate the transport in Fe$_4$GeTe$_2$, a prototypical high-temperature 2D ferromagnet. We show that, while spin transport remains essentially single-particle under moderate bias, inelastic spin-dependent scattering of carriers with particle-hole excitations drives a distinctive hot-correlated electron regime beyond a critical voltage. This regime is marked by incoherent features in both the electronic spectrum and the conductance, which are experimentally accessible. Our results demonstrates that material-specific many-body non-equilibrium methods are essential for a complete understanding of spin transport in 2D ferromagnets.

We demonstrate that the Mott transition exhibits universal scaling as a consequence of the breaking of a $\mathbb{Z}_2$ symmetry in momentum space. A direct consequence of this discrete symmetry breaking is the charge or Mott gap itself. From extensive numerics, we proffer that it is the charge compressibility that acts as the underlying order parameter as it is zero in the insulator and non-zero in the metallic state. Additionally, the Widom line (temperature of the extremum of the compressibility) obeys a universal scaling of $T_m=0.39U$ deep into the insulating state directly from $Z_2$ universality. Furthermore, the temperature at which the second derivative of the compressibility has a minimum is independent of lattice geometry, exhibiting a universal scaling of $|U-U_c|^\alpha$ where $\alpha\approx 1$. Finally, our computational approach reproduces the key features of the doping dependence of the compressibility demonstrated in recent cold-atom quantum simulators of the Hubbard model, thereby corroborating our conclusions on $\mathbb{Z}_2$ universality.

The Hubbard model on the Kagome lattice is a widely used interacting model for describing the electronic properties of various transition metal-based Kagome materials. We find altermagnetism driven by Coulomb interaction in the Kagome Hubbard model at Dirac filling with no spin-orbit coupling nor explicit spatial symmetry breaking present. We show how this insulating altermagnet is relevant to other lattices with larger unit cells such as the Lieb-Kagome lattice. The ALM found displays a characteristic magnon splitting which can be detected in inelastic neutron scattering experiments on interacting Kagome materials.

X-ray and neutron diffraction are foundational tools for determining crystal structures, but their resolution limits can lead to misassignments, especially in materials with subtle distortions or competing phases. Here, we demonstrate the use of nonlinear transport as a complementary approach to uncover hidden crystal symmetries, using the strongly correlated Ca$_3$Ru$_2$O$_7$ as a case study. Below 48 K (T$_S$), where the magnetic moments of the antiferromagnetic phase reorient from the a- to the b-axis, leading to a pseudogap opening, our measurements, with support of DFT, reveal a previously overlooked lower-symmetry phase. This is manifested by the emergence of longitudinal nonlinear resistance (NLR) along the b-axis below T$_S$, providing direct evidence of combined translational and time-reversal symmetry breaking. This response also suggests a transformation from a conventional antiferromagnet into an altermagnet. The lower-symmetry phase arises from a subtle lattice distortion (~0.1 pm) associated with the magnetic transition at T$_S$, below the detection limit of conventional diffraction. Moreover, this NLR below T$_S$ is accompanied by a nonlinear Hall effect, both of which are enhanced by the large quantum metric associated with Weyl chains near the Fermi surface. Our findings demonstrate nonlinear transport as a sensitive probe of hidden symmetry breaking and altermagnetism, complementing and extending beyond the reach of traditional diffraction and spectroscopic techniques.

Resonant Inelastic X-Ray Scattering (RIXS) is a well-established tool for probing excitations in a wide range of materials. The measured spectra strongly depend on the scattering geometry, via its influence on the polarization of the incoming and outgoing light. By employing a tensor representation of the 4-point response function that governs the RIXS intensity, we disentangle the experimental geometry from the intrinsic material properties. In dipole-dipole RIXS processes and low-symmetry crystals, up to 81 linearly independent fundamental spectra can be measured as a function of light polarization. However, for crystals or molecules with symmetry, the number of independent fundamental spectra that define the RIXS tensor is significantly reduced. This work presents a systematic framework for determining the number of fundamental spectra and expressing the RIXS tensor in terms of these fundamental components. Given a specific experimental geometry, the measured spectrum can be represented as a linear combination of these fundamental spectra. To validate our approach, we performed calculations for different point group symmetries, both with and without an applied magnetic field. Within the same framework, we derived expressions for powder spectra in momentum-independent processes and spectra obtained using Bragg spectrometers. This formalism provides a valuable toolkit for optimizing experiment planning, data interpretation, and RIXS simulation.

Accurate simulations of the Hubbard model are crucial to understanding strongly correlated phenomena, where small energy differences between competing orders demand high numerical precision. In this work, Neural Quantum States are used to probe the strongly coupled and underdoped regime of the square-lattice Hubbard model. We systematically compare the Hidden Fermion Determinant State and the Jastrow-Backflow ansatz, parametrized by a Vision Transformer, finding that in practice, their accuracy is similar. We also test different symmetrization strategies, finding that output averaging yields the lowest energies, though it becomes costly for larger system sizes. On cylindrical systems, we consistently observe filled stripes. On the torus, our calculations display features consistent with a doped Mott insulator, including antiferromagnetic correlations and suppressed density fluctuations. Our results demonstrate both the promise and current challenges of neural quantum states for correlated fermions.

The Stoner criterion for ferromagnetism arises from interaction-driven asymmetric filling of spin bands, requiring that the spin susceptibility: (i) peaks dominantly at $\mathbf{Q}=\bm{0}$; and (ii) diverges at a critical interaction strength. Here, we demonstrate that this Stoner mechanism breaks down due to competition with altermagnetic orders, even when both conditions are met. Altermagnetism in solids is characterized by collinear antiparallel spin alignment that preserves translational symmetry, and inherently fulfills these requirements. As a proof of concept, we study a two-orbital Hubbard model with electron filling near Van Hove singularities at high-symmetry momenta. Our results reveal that orbital-resolved spin fluctuations, amplified by strong inter-orbital hopping, stabilize intrinsic altermagnetic order. A quantum phase transition from altermagnetism to ferromagnetism occurs at critical Hund's coupling $J_H$. We further propose directional spin conductivity anisotropy as a detectable signature of this transition via non-local spin transport. This work establishes the pivotal role of altermagnetism in correlated systems.

The fermionic Hubbard model, when combined with the ingredient of frustration, associated with the breaking of particle-hole symmetry, harbors a rich phase diagram. Aspects of theoretical findings associated with the nature of magnetism and metallicity, in a diverse set of parameter regimes, are now being actively investigated in triangular Hubbard cold atom and solid-state (moir\'e) based emulators. Building on the theoretical work of Haerter and Shastry [Phys. Rev. Lett. 95,087202 (2005)], we explore the impact of kinetically frustrated magnetism, a phenomenon where antiferromagnetic order emerges without any underlying magnetic interactions, at finite hole density. We numerically study the infinite-$U$ triangular Hubbard model using the density matrix renormalization group algorithm and estimate the extent of stability of the kinetically induced $120^{\circ}$ antiferromagnetic state to hole doping. Beyond the Haerter-Shastry regime, we find an intermediate phase with multimer (involving multiple correlated spins) stripes that eventually gives way to a paramagnet. We also find evidence of gapless charge excitations (metallicity) throughout the phase diagram for finite hole density. We discuss the implications at large, but finite and realistic values of $U/t$, and investigate whether kinetic magnetism and superexchange collaborate or compete.

We develop an extension of the fermionic functional renormalization group for systems where strong correlations give rise to projected Hilbert spaces. We use our method to calculate the phase diagram and the electronic spectral function of the Hubbard model at infinite on-site repulsion where many-body states involving doubly occupied lattice sites are eliminated from the physical Hilbert space. For a square lattice with nearest-neighbor hopping we find that the ground state evolves from a paramagnetic Fermi liquid at low densities via a state with antiferromagnetic stripe order at intermediate densities to an extended Nagaoka ferromagnet at high densities. In the strongly correlated magnetic phases, the electrons form an incoherent non-Fermi liquid. Both at high and low densities, the volume of the Fermi surface is not constrained by Luttinger's theorem.

Accurately describing the ground state of strongly correlated systems is essential for understanding their emergent properties. Neural Network Backflow (NNBF) is a powerful variational ansatz that enhances mean-field wave functions by introducing configuration-dependent modifications to single-particle orbitals. Although NNBF is theoretically universal in the limit of large networks, we find that practical gains saturate with increasing network size. Instead, significant improvements can be achieved by using a multi-determinant ansatz. We explore efficient ways to generate these multi-determinant expansions without increasing the number of variational parameters. In particular, we study single-step Lanczos and symmetry projection techniques, benchmarking their performance against diffusion Monte Carlo and NNBF applied to alternative mean fields. Benchmarking on a doped periodic square Hubbard model near optimal doping, we find that a Lanczos step, diffusion Monte Carlo, and projection onto a symmetry sector all give similar improvements achieving state-of-the-art energies at minimal cost. By further optimizing the projected symmetrized states directly, we gain significantly in energy. Using this technique we report the lowest variational energies for this Hamiltonian on $4\times 16$ and $4 \times 8$ lattices as well as accurate variance extrapolated energies. We also show the evolution of spin, charge, and pair correlation functions as the quality of the variational ansatz improves.

Recently, a non-Hermitian Anderson impurity model with one-body loss has been studied in [Phys. Rev. B 111, 125157 (2025)}], and it has been demonstrated that the renormalization effect generated by strong correlations counterintuitively changes the nature of dissipation into an emergent many-body dissipation that causes a Kondo breakdown. In a closely related context, it is also known that two-body loss in a non-Hermitian Kondo model triggers the Kondo breakdown. To elucidate the essence of these phenomena, we study the Anderson impurity model with a non-Hermitian complex hybridization as an effective model that provides a simple understanding of the Kondo breakdown. Using the slave-boson mean-field theory, we show that this model can explain the Kondo breakdown with a single complex parameter. Furthermore, we provide the exact Bethe ansatz solutions that support the results obtained by the slave-boson mean-field theory. Finally, we point out that the Lehmann representation for the non-Hermitian Green function cannot be obtained by the analytic continuation to the complex energy upon the Kondo breakdown, where the analyticity of the non-Hermitian Green function in the half-complex-$\omega$ plane no longer holds.

The Hubbard model has occupied the minds of condensed matter physicists for most part of the last century. This model provides insight into a range of phenomena in correlated electron systems. We wish to examine the paradigm of quantum algorithms for solving such many-body problems. The focus of our current work is on the one-dimensional model which is integrable, meaning that there exist analytical results for determining its ground state. In particular, we demonstrate how to perform a gate-based quantum computer simulation of quantum annealing for the Hubbard Hamiltonian. We perform simulations for systems with up to 40 qubits to study the scaling of required annealing time for obtaining the ground state. We find that for the half-filled cases considered, there is a substantial quantum speed-up over algorithms based on the Bethe-ansatz equations.

The use of magnetic nanographenes as building blocks for artificial spin lattices is enabling the exploration of flagship model Hamiltonians of one-dimensional quantum magnetism with an unprecedented degree of control. The spin-1 Heisenberg model, incorporating both linear and quadratic exchange interactions, was first realized using [3]-triangulenes, where the hallmark Haldane phase with spin fractionalization was observed. Later, the spin-1/2 Heisenberg Hamiltonian with exchange alternation was realized with Clar's goblets, where two additional topological phases were identified. Here we show that spin-1 nanographenes can also be used to explore the topological phase transition between the Haldane phase and a dimerized phase predicted for spin-1 chains with bond-alternation. We first study how the boundary of the phase transition is modified by non-linear exchange, known to be present in spin-1 nanographenes, using density matrix renormalization group (DMRG). Combining multiconfigurational with first-principles calculations, we propose two candidates to realize different topological phases of the model: a recently synthesized extended Clar's goblet, and a passivated [4]-triangulene. Moreover, we show how these two phases can be identified experimentally using inelastic electron tunneling spectroscopy (IETS). This work paves the way for the experimental realization of these topological phases, which can be locally probed with scanning tunneling microscopy.

We experimentally study the two-dimensional Fermi-Hubbard model using a Rydberg-based quantum processing unit in the analog mode. Our approach avoids encoding directly the original fermions into qubits and instead relies on reformulating the original model onto a system of fermions coupled to spins and then decoupling them in a self-consistent manner. We then introduce the auxiliary spin solver: this hybrid quantum-classical algorithm handles a free-fermion problem, which can be solved efficiently with a few classical resources, and an interacting spin problem, which can be naturally encoded in the analog quantum computer. This algorithm can be used to study both the equilibrium Mott transition as well as non-equilibrium properties of the original Fermi-Hubbard model, highlighting the potential of quantum-classical hybrid approaches to study strongly correlated matter.

Recent observations of quantum anomalous Hall effects in moir\'e systems have revealed the emergence of interaction-driven ferromagnetism with significant orbital contributions. To capture this physics, we extend the modern theory of orbital magnetization to Hartree-Fock states and show that the standard expression remains valid with Hartree-Fock orbitals and Hamiltonians. We then benchmark our theory against the Kane-Mele-Hubbard model in a weak field, which yields excellent agreement with direct numerical calculations. Applying our theory to twisted MoTe$_2$ bilayers, we find orbital magnetization of order one Bohr magneton per moir\'e cell with a non-monotonic twist-angle dependence. Our work establishes a general theory of orbital magnetization in interacting moir\'e systems and provides quantitative guidance for interpreting recent experiments.

We introduce TeMFpy, a Python library for converting fermionic mean-field states to finite or infinite matrix product state (MPS) form. TeMFpy includes new, efficient, and easy-to-understand algorithms for both Slater determinants and Pfaffian states. Together with Gutzwiller projection, these also allow the user to build variational wave functions for various strongly correlated electron systems, such as quantum spin liquids. We present all implemented algorithms in detail and describe how they can be accessed through TeMFpy, including full example workflows. TeMFpy is built on top of TeNPy and, therefore, integrates seamlessly with existing MPS-based algorithms.

We study the effect of density-assisted hopping on different dimerized lattice geometries, such as bilayers and ladder structures. We show analytically that the density-assisted hopping induces an attractive interaction in the lower (bonding) band of the dimer structure and a repulsion in the upper (anti-bonding) band. Overcoming the onsite repulsion, this can lead to the appearance of superconductivity. The superconductivity depends strongly on the filling, and present a pairing structure more complex than s-wave pairing. Combining numerical and analytical methods such as the matrix product states ansatz, bosonization and perturbative calculations we map out the phase diagram of the two-leg ladder system and identify its superconducting phase. We characterize the transition from the non-density-assisted repulsive regime to the spin-gapped superconducting regime as a Berezinskii-Kosterlitz-Thouless transition.

Two-dimensional van der Waals (vdW) ferromagnets drive the advancement in spintronic applications and enable the exploration of exotic magnetism in low-dimensional systems. The entanglement of dual $-$ localized and itinerant $-$ nature of electrons lies at the heart of the correlated electron systems giving rise to exotic ground state properties such as complex magnetism, heavy fermionic behavior, Kondo lattice formation, \textit{etc}. Through temperature-dependent electronic structure of vdW ferromagnets, (Co$_{x}$Fe$_{1-x}$)$_{3}$GeTe$_{2}$, probed using high-resolution photoemission spectroscopy and density functional theory combined with dynamical mean field theory (DFT+DMFT), we provide direct evidence of the emergence of Kondo resonance peak driven by complex interplay between localized and itinerant electrons. In overall agreement with experimental electronic structure and magnetic properties, DFT+DMFT also reveals finite spin band splitting well beyond $T_{C}$. Core levels, valence band photoemission spectra together with DFT+DMFT spectral functions reveal insignificant change across $T_{C}$ indicating non-Stoner magnetism in (Co$_{x}$Fe$_{1-x}$)$_{3}$GeTe$_{2}$. Our results provide a way forward to the understanding of complex interplay between electronic structure, exotic magnetism and heavy fermionic behavior leading to Kondo scenerio in 3$d$ vdW ferromagnets.

Altermagnetic (AM) compounds have recently emerged as a promising platform for realizing unconventional quantum phases, enabled by their unique spin-split band structure at zero net magnetization. Here, we present a first-principles investigation of magnetically intercalated transition metal dichalcogenides (TMDs) of the form XY$_4$Z$_8$ (X $=$ Mn, Fe, Co, Ni, Cr, or V; Y $=$ Nb or Ta; and Z $=$ Se or S), identifying a subset of new versatile AM candidates. Our results establish a direct correlation between interatomic geometry, quantified by the ratio of interlayer to intralayer spacing, and the selection of magnetic ground states. Systems with A-type antiferromagnetic order exhibit momentum-dependent spin splitting consistent with AM behavior. Crucially, the combination of the AM spin-splitting and the spin-orbit coupling leads to the emergence of Weyl nodes together with the corresponding topological Fermi arc surface states. Moreover, we identify flat bands near the Fermi level that originate from the intercalant-induced formation of an effective kagome-like sublattice in the TMD layer. These results collectively establish magnetically intercalated TMDs as a promising platform for engineering altermagnetism, flat bands, and Weyl fermions within a single material family, facilitating the development of topological and spintronic applications.

We investigate the Sn-substitution effect on the thermoelectric transport properties of the layered cobaltite [Ca$_2$CoO$_3$]$_{0.62}$[CoO$_2$] single crystals, which exhibit a non-monotonic temperature variation of the electrical resistivity and the Seebeck coefficient owing to the complex electronic and magnetic states. We find that the onset temperature of the short-range spin-density-wave (SDW) formation increases with the substituted Sn content, indicating the impurity-induced stabilization of the SDW order, reminiscent of the disorder/impurity-induced spin order in the cuprate superconductors. The Seebeck coefficient is well related to such impurity effects, as it is slightly enhanced below the onset temperatures, implying a decrease in the carrier concentration due to a pseudo-gap formation associated with the SDW ordering. We discuss the site-dependent substitution effects including earlier studies, and suggest that substitution to the conducting CoO$_2$ layers is essential to increase the onset temperature, consistent with the impurity-induced SDW picture realized in the conducting layers with the cylindrical Fermi surface.

We introduce a Markov Chain Monte Carlo (MCMC) algorithm that dramatically accelerates the simulation of quantum many-body systems, a grand challenge in computational science. State-of-the-art methods for these problems are severely limited by $O(N^3)$ computational complexity. Our method avoids this bottleneck, achieving near-linear $O(N \log N)$ scaling per sweep. Our approach samples a joint probability measure over two coupled variable sets: (1) particle trajectories of the fundamental fermions, and (2) auxiliary variables that decouple fermion interactions. The key innovation is a novel transition kernel for particle trajectories formulated in the Fourier domain, revealing the transition probability as a convolution that enables massive acceleration via the Fast Fourier Transform (FFT). The auxiliary variables admit closed-form, factorized conditional distributions, enabling efficient exact Gibbs sampling update. We validate our algorithm on benchmark quantum physics problems, accurately reproducing known theoretical results and matching traditional $O(N^3)$ algorithms on $32\times 32$ lattice simulations at a fraction of the wall-clock time, empirically demonstrating $N \log N$ scaling. By reformulating a long-standing physics simulation problem in machine learning language, our work provides a powerful tool for large-scale probabilistic inference and opens avenues for physics-inspired generative models.

A multichannel Kondo system, where a single quantum spin couples to multiple channels of an electronic bath, provides one of the simplest examples of a zero-dimensional non-Fermi liquid. It is natural to ask: what happens when an extensive number of such systems are coupled together? A simple renormalization group argument implies that in a chain of SU(N) multichannel quantum systems, where each spin is coupled to its own bath of K channels, the individual spins dynamically decouple at low energy when N>K, resulting in a 'sliding' non-Fermi liquid. Using Quantum Monte Carlo (QMC) simulations, we find evidences of a continuous, 'dimensionality-changing' phase transition out of this non-Fermi liquid into a valence-bond solid phase as the intersite coupling is increased. Remarkably, at the critical point, correlations exhibit a power-law behavior even along the direction in which the spins are coupled, indicating the breakdown of dynamical decoupling at the transition. We also develop an RG scheme to understand the universal aspects of this transition.

We introduce the concept of \emph{orbital altermagnetism}, a symmetry-protected magnetic order of pure orbital degrees of freedom. It is characterized with ordered anti-parallel orbital magnetic moments in real space but momentum-dependent orbital band splittings, analogous to spin altermagnetism. Using a minimal tight-binding model with complex hoppings in a square-kagome lattice, we show that such order inherently arises from staggered loop currents, producing a $d$-wave-like orbital-momentum locking. First-principles calculations show that orbital altermagnetism emerges independent of spin ordering in in-plane ferromagnets of CuBr$_2$ and VS$_2$, so that it can be unambiguously identified experimentally. On the other hand, it may also coexist with spin altermagnetism, such as in monolayer MoO and CrO. The orbital altermagnetism offers an alternative platform for symmetry-driven magnetotransport and orbital-based spintronics, as exemplified by large nonlinear current-induced orbital magnetization.

The presence of a boundary enriches the nature of quantum phase transitions. However, the boundary critical phenomena in topological superconductors remain underexplored so far. Here, we investigate the boundary criticality in a two-dimensional correlated time-reversal-invariant topological superconductor tuned through a quantum phase transition into a trivial time-reversal-breaking superconductor. Using sign-problem-free determinant quantum Monte Carlo simulations, we chart the quantum phase diagram and reveal the boundary criticalities encompassing ordinary, special, and extraordinary transitions. Additionally, using renormalization group analysis, we compute the boundary critical exponent up to two loops. Remarkably, the simulations and two-loop renormalization group calculations consistently demonstrate that the presence of the boundary Majorana fermion at the special transition gives rise to a new type of boundary Gross-Neveu-Yukawa fixed point. We conclude with a discussion of possible experimental realizations in iron-based superconductors.

Altermagnets, a novel class of collinear magnetic materials, exhibit unique spin-split band structures, yet topological insulating states in intrinsic altermagnetic systems are rare. Here, we identify monolayer Fe$_2X_2$O ($X$ = Cl, Br, I) as a new family of 2D altermagnetic real Chern insulators. These materials display robust $d$-wave altermagnetic ordering, semiconducting band gaps, and nontrivial real Chern numbers per spin channel, yielding spin-polarized topological corner modes. They also feature spin-polarized valleys with strong altermagnetism-valley-spin-lattice coupling, enabling valley-selective excitation via linear dichroism and strain-induced valley polarization. In multiferroic Fe$_2$Cl$_2$O, magnetism coexists with ferroelasticity, and an applied strain can switche the N\'eel vector. These findings position 2D iron oxyhalides as a promising platform for exploring altermagnetism and magnetic topological states for spintronics and valleytronics.

Altermagnets, a novel class of collinear magnetic materials, exhibit unique spin-split band structures, yet topological insulating states in intrinsic altermagnetic systems are rare. Here, we identify monolayer Fe$_2X_2$O ($X$ = Cl, Br, I) as a new family of 2D altermagnetic real Chern insulators. These materials display robust $d$-wave altermagnetic ordering, semiconducting band gaps, and nontrivial real Chern numbers per spin channel, yielding spin-polarized topological corner modes. They also feature spin-polarized valleys with strong altermagnetism-valley-spin-lattice coupling, enabling valley-selective excitation via linear dichroism and strain-induced valley polarization. In multiferroic Fe$_2$Cl$_2$O, magnetism coexists with ferroelasticity, and an applied strain can switche the N\'eel vector. These findings position 2D iron oxyhalides as a promising platform for exploring altermagnetism and magnetic topological states for spintronics and valleytronics.

Localized $5d^2$ electrons in a cubic crystal field possess multipoles such as electric quadrupoles and magnetic octupoles. We studied the cubic double perovskite Ba$_2$CaOsO$_6$ containing the Os$^{6+}$ ($5d^2$) ions, which exhibits a phase transition to a `hidden order' below $T^* \sim$ 50 K, by X-ray absorption spectroscopy (XAS) and X-ray magnetic circular dichroism (XMCD) at the Os $L_{2,3}$ edge. The cubic ligand-field splitting between the $t_{2g}$ and $e_g$ levels of Os $5d$ was deduced by XAS to be $\sim$4 eV. The temperature dependence of the XMCD spectra was consistent with a $\sim$18 meV residual cubic splitting of the lowest $J_{\rm eff} =$ 2 multiplet state into the non-Kramers $E_g$ doublet ground state and the $T_{2g}$ triplet excited state. Ligand-field (LF) multiplet calculation under fictitious strong magnetic fields indicated that the exchange interaction between nearest-neighbor octupoles should be as strong as $\sim$1.5 meV if a ferro-octupole order is stabilized in the `hidden-ordered' state, consistent with the exchange interaction of $\sim$1 meV previously predicted theoretically using model and density functional theory calculations.

Skyrmions are emergent many-body excitations that lie at the heart of both multi-component quantum Hall-like systems and deconfined quantum criticality. In a companion article (X. Wang et al., arXiv:2507.22971), we studied a microscopic time-reversal symmetric model of tunable interacting Chern bands using numerically exact determinant quantum Monte Carlo calculations, and presented evidence for the emergence of robust skyrmion excitations. These charged excitations emerge in the vicinity of a many-body insulator at a commensurate filling of the Chern bands, and lead to the onset of superconductivity when doped away from the insulating phase. Here, we present quantum Monte-Carlo results and a complementary field-theoretical analysis for the quantum phase transition(s) that arise between the intertwined phases as a function of two distinct tuning parameters. Our numerical results are consistent with a single continuous quantum phase transition between an insulating Chern antiferromagnet and a fully gapped superconductor, with an emergent SO(5) symmetry at the putative critical point, highly suggestive of deconfined quantum (pseudo-)criticality. We also present a detailed comparison between the momentum-resolved spectral functions associated with the neutral collective modes, single electron and composite spin-polaron excitations obtained using a combination of Monte-Carlo computations and a Bethe-Salpeter analysis built on top of the self-consistent Hartree-Fock calculation. We end with a brief outlook on some of the interesting open problems.

Spin current generated electrically is among the core phenomena of spintronics for driving high-performance spin device applications. Here, on the basis of systematic investigations for the hole doped single-layer Haldane-Hubbard(HH) model, we propose a new bilayer HH model to realize the compensated odd-parity spin splitting and the $T$-even spin Hall conductivity where the two layers are connected by the time reversal transformation. Our results show that the vanishing layer-dependent electric potential $V_{L}$ gives rise to odd-parity ALM protected by the combined symmetry $TM_{xy}$ with $T$ and $M_{xy}$ being the time reversal and mirror reflection perpendicular to $z$ axis, and the $T$-even spin Hall conductivity simultaneously. In addition, though the staggered magnetization within each layer is substantially impacted by the layer-dependent electric potential, small $V_{L}$'s only bring negligible changes to the net magnetization and the spin Hall conductivity, indicating that the alternating spin splitting in momentum space and the spin Hall conductivity are insusceptible to external elements. Most importantly, our work provides a general framework for the simultaneous realization of the compensated odd-parity spin splitting in momentum space and the spin Hall conductivity in collinear magnets, in terms of stacked multi-layer systems.

The development of altermagnets is fundamentally important for advancing spintronic device technology, but remains unpractical for the weak spin splitting in most cases, especially in two-dimensional materials. Based on spin group symmetry analysis and first-principles calculations, a novel hydroxyl rotation strategy in collinear antiferromagnets has been proposed to design altermagnets. This approach achieves a large chirality-reversible spin splitting exceeding $1130$ meV in $\alpha_{60}$-Mn$_2$B$_2$(OH)$_2$ monolayer. The system also exhibits intrinsic features of a node-line semimetal in the absence of spin-orbit coupling. Besides, the angles of hydroxyl groups serve as the primary order parameter, which can switch on/off the altermagnetism coupled with the ferroelastic mechanism. The corresponding magnetocrystalline anisotropy have also been modulated. Moreover, an interesting spin-related transport property with the spin-polarized conductivity of 10$^{19}$ $\Omega^{-1}m^{-1}s^{-1}$ also emerges. These findings uncover the hydroxyl rotation strategy as a versatile tool for designing altermagnetic node-line semimetals and opening new avenues for achieving exotic chemical and physical characteristics associated with large spin splitting.

We propose a new pathway to the quantized anomalous Hall effect (QAHE) by coupling an altermagnet to a topological crystalline insulator (TCI). The former gaps the topological surface states of the TCI, thereby realizing the QAHE in a robust and switchable platform with near- vanishing magnetization. We demonstrate the feasibility of this approach by studying a slab of the TCI SnTe coupled to an altermagnetic RuO2 layer. Our first-principles calculations reveal that the d-wave altermagnetism in RuO2 induces a 7 meV gap to the Dirac surface states on the (110) surface of SnTe, producing a finite anomalous Hall effect. Our approach generalizes to broader classes of altermagnetic materials and TCIs, thereby providing a family of topological altermagnetic heterostructures with small or vanishing magnetization that support nontrivial Chern numbers. Our results highlight a promising new topological platform with great tunability and applications to spintronics.

The interplay between non-trivial band topology and strong electronic correlations is a central challenge in modern condensed matter physics. We investigate this competition on a two-leg ladder model with a p-wave-like hybridisation between the legs. This model hosts a symmetry-protected topological phase in its non-interacting limit. Using the density-matrix renormalisation group algorithm, we compute the comprehensive quantum phase diagram in the presence of a repulsive inter-leg density-density interaction. Our analysis, based on entanglement entropy and the entanglement spectrum, reveals a fascinating dichotomy in the stability of the topological phase. We find a non-trivial change in the value of the edge entanglement entropy as we include interaction. Furthermore, we find that the phase boundary separating a trivial insulator phase and a topological one with winding number two remains robustly pinned at its non-interacting location, irrespective of the interaction strength. Variation of the effective conformal field theory's central charge near the critical line explains the robustness of the gap. In contrast, the transition to an insulating phase with winding number one is heavily renormalised, with the critical line shifting significantly as the interaction increases. By successfully mapping the phase diagram and identifying the distinct behaviours of the phase boundaries, our work clarifies how interactions can selectively preserve or destroy different aspects of a topological phase.

Multiorbital Hubbard models with Hund coupling and crystal-field splitting exhibit an instability toward spin-triplet excitonic order in the parameter regime characterized by strong local spin fluctuations. Upon chemical doping, two distinct types of excitonic ferromagnetism have been reported. Using steady-state nonequilibrium dynamical mean-field theory, we demonstrate that photo-doped half-filled systems can host nonthermal counterparts of these excitonic phases and exhibit a rich phase diagram in the space of photo-doping and crystal field splitting. Photo-doping a spin-triplet excitonic insulator provides a route towards nonequilibrium control of magnetic order.

Altermagnets exhibit momentum-dependent spin splitting without net magnetization, combining characteristics of both ferromagnets and antiferromagnets, making them highly interesting for spintronics applications. CrSb is a prime candidate with a high N\'eel temperature ($\sim700$~K) and a large exchange-driven splitting of $\sim0.6$--1~eV. Using ab-initio calculations, we consider slabs of various orientations in the ultrathin limit. We find that (100) oriented slabs have spin-degenerate bands. In (0001) oriented slabs, the exchange-driven altermagnetic spin splitting collapses, but including spin-orbit coupling restores a residual anisotropic splitting of $\sim70$~meV. In contrast, the (110) oriented slabs show an altermagnetic spin splitting of $\sim400$~meV, and emerges as a robust candidate for realizing large, exchange-driven altermagnetism

Drawing the quantum phase diagram of a many-body system in the parameter space of its Hamiltonian can be seen as a learning problem, which implies labelling the corresponding ground states according to some classification criterium that defines the phases. In this work we adopt unsupervised learning, where the algorithm has no access to any priorly labeled states, as a tool for determining quantum phase diagrams of many-body systems. The algorithm directly works with quantum states: given the ground-state configurations for different values of the Hamiltonian parameters, the process uncovers the most significant way of grouping them based on a similarity criterion that refers to the fidelity between quantum states, that can be easily estimated, even experimentally. We benchmark our method with two specific spin-$\frac{1}{2}$ chains, with states determined via tensor network techniques. We find that unsupervised learning algorithms based on spectral clustering, combined with ``silhouette'' and ``elbow'' methods for determining the optimal number of phases, can accurately reproduce the phase diagrams. Our results show how unsupervised learning can autonomously recognize and possibly unveil novel phases of quantum matter.

We present a machine learning (ML) force-field framework for simulating the non-equilibrium dynamics of charge-density-wave (CDW) order driven by the Peierls instability. Since the Peierls distortion arises from the coupling between lattice displacements and itinerant electrons, evaluating the adiabatic forces during time evolution is computationally intensive, particularly for large systems. To overcome this bottleneck, we develop a generalized Behler-Parrinello neural-network architecture -- originally formulated for ab initio molecular dynamics -- to accurately and efficiently predict forces from local structural environments. Using the locality of electronic responses, the resulting ML force field achieves linear scaling efficiency while maintaining quantitative accuracy. Large-scale dynamical simulations using this framework uncover a two-stage coarsening behavior of CDW domains: an early-time regime characterized by a power-law growth $L \sim t^{\alpha}$ with an effective exponent $\alpha \approx 0.7$, followed by a crossover to the Allen-Cahn scaling $L \sim \sqrt{t}$ at late times. The enhanced early-time coarsening is attributed to anisotropic domain-wall motion arising from electron-mediated directional interactions. This work demonstrates the promise of ML-based force fields for multiscale dynamical modeling of condensed-matter lattice models.

Altermagnets constitute an emerging class of collinear magnets that exhibit zero net magnetization yet host spin-split electronic bands arising from non-relativistic spin-space-group symmetries. Realization of altermagnetism in the two-dimensional (2D) limit remains an outstanding challenge because dimensional reduction suppresses kZ dispersion and destabilizes the symmetry operations essential for spin compensation. Here, we demonstrate genuine 2D altermagnetism in epitaxial unit-cell-thin films of CrSb grown on Bi2Te3. It reveals a thickness-driven transition from a ferrimagnetic state in 1-unit-cell films to an altermagnetic state above a critical thickness of 7/4 unit cell. The transition originates from interfacial symmetry breaking at the Cr-terminated layer that induces local moment imbalance. With increasing thickness the key spin-space-group symmetries [C2||C6Zt] and [C2||MZ] restores, which leads to altermagnetism with zero net magnetization and momentum-dependent spin splitting. Our results provide the first experimental realization of altermagnetism in the 2D regime and establish a route for integrating stray-field-free spin order into nanoscale spintronic architectures.

Surface oxides are associated with two-level systems (TLSs) that degrade the performance of niobium-based superconducting quantum computing devices. To address this, we introduce a predictive framework for selecting metal capping layers that inhibit niobium oxide formation. Using DFT-calculated oxygen interstitial and vacancy energies as thermodynamic descriptors, we train a logistic regression model on a limited set of experimental outcomes to successfully predict the likelihood of oxide formation beneath different capping materials. This approach identifies Zr, Hf, and Ta as effective diffusion barriers. Our analysis further reveals that the oxide formation energy per oxygen atom serves as an excellent standalone descriptor for predicting barrier performance. By combining this new descriptor with lattice mismatch as a secondary criterion to promote structurally coherent interfaces, we identify Zr, Ta, and Sc as especially promising candidates. This closed-loop strategy integrates first-principles theory, machine learning, and limited experimental data to enable rational design of next-generation materials.

The kagome lattice of spin-1/2 Cu atoms in herbertsmithite (ZnCu3(OH)6Cl2) may sustain a quantum spin liquid (QSL) state with spinon quasiparticles. Each kagome plane is separated from its homologues by a layer of spinless Zn atoms. Providentially, however, some spin-1/2 Cu atoms substitute randomly onto these inter-kagome Zn sites. We reconceptualize these 'impurity' atoms as quantum 'witness-spins', an exceptional new interrogative of the conjectured Z2-gauge-symmetric QSL state. Thus we introduce spin-noise spectroscopy to explore herbertsmithite witness-spin dynamics for QSL studies. It reveals the existence, slowing and intensification of spin noise, prefatory to a sharp transition at T* {\approx} 260 mK. Below T* the spin-noise power spectral density S_M({\omega},T) {\propto} {\omega}^{-{\alpha}(T)} stabilizes at {\alpha} {\approx} 1; the spin noise variance {\sigma}_M^2(T) diminishes precipitously; the ultra-low-field magnetic susceptibility {\chi}(T) undergoes a sharp transition into a phase exhibiting an Edwards-Anderson order-parameter and ultra-slow spin-state relaxation. A Z2 QSL theory of spinon-mediated witness-spin interactions corresponds best to all these experimental observations, predicting slowing and intensification of witness-spin fluctuations and noise spectrum S_M({\omega},T) with cooling, with a transition into a unique spinon-mediated phase signified by rapidly diminishing spin noise, with S_M({\omega},T) {\propto} {\omega}^{-1}, a sharp cusp in the DC magnetic susceptibility {\chi}(T), and the appearance of an Edwards-Anderson order-parameter. We rule out numerous other mechanisms for these effects, so that only spinon-mediation by either a Z2 or U(1) QSL is consistent with all present herbertsmithite empirics, with the former model providing a closer match to data.

Surface oxides are associated with two-level systems (TLSs) that degrade the performance of niobium-based superconducting quantum computing devices. To address this, we introduce a predictive framework for selecting metal capping layers that inhibit niobium oxide formation. Using DFT-calculated oxygen interstitial and vacancy energies as thermodynamic descriptors, we train a logistic regression model on a limited set of experimental outcomes to successfully predict the likelihood of oxide formation beneath different capping materials. This approach identifies Zr, Hf, and Ta as effective diffusion barriers. Our analysis further reveals that the oxide formation energy per oxygen atom serves as an excellent standalone descriptor for predicting barrier performance. By combining this new descriptor with lattice mismatch as a secondary criterion to promote structurally coherent interfaces, we identify Zr, Ta, and Sc as especially promising candidates. This closed-loop strategy integrates first-principles theory, machine learning, and limited experimental data to enable rational design of next-generation materials.

We extend the previous study of extracting crystalline symmetry-protected topological invariants to the correlated regime. We construct the interacting Hofstadter model defined on square lattice with the rotation and translation symmetry defects: disclination and dislocation. The model realizes Chern insulator and the charge density wave state as one tunes interactions. Employing the density matrix renormalization group (DMRG) method, we calculate the excess charge around the defects and find that the topological invariants remain quantized in both phases, with the topological quantity extracted to great precision. This study paves the way for utilizing matrix product state, and potentially other quantum many-body computation methods, to efficiently study crystalline symmetry defects on 2D interacting lattice systems.

Excitement about the magnetic and electronic properties of RuO$_2$ is growing, fueled by reports of antiferromagnetism, strain-induced superconductivity, and its recent classification as a member of a newly proposed magnetic class, altermagnets, with RuO$_2$ widely regarded as the paradigmatic example. Nevertheless, the magnetic ground state of RuO$_2$ remains contentious, as several recent experiments report no evidence of magnetic order. To address this discrepancy, we performed resonant elastic scattering measurements on a series of epitaxial RuO$_2$ thin films grown on the (100)-plane of TiO$_2$ substrates across a range of strain states. Leveraging full polarization control and azimuthal scans of the structurally forbidden 100 Bragg reflection, we systematically tested for signatures of colinear antiferromagnetic order. We found that the resonant elastic scattering signal in RuO$_2$ thin films likely originates from anisotropic charge scattering, not long-range antiferromagnetic order. Using angle-resolved photoemission spectroscopy we uncover a band structure without altermagnetic band splitting that is consistent with a nonmagnetic phase. Similarly, anisotropic magnetoresistance results show no evidence of magnetism. The combination of three independent measurements suggests the absence of altermagnetism in RuO$_2$.

Altermagnets are fully compensated collinear antiferromagnets that lack the combined time-reversal and translation symmetry. Here we show that their symmetry allows for a switchable ferro-spinetic polarization - the spin analogue of ferroelectricity - in a direction dictated by the lattice symmetry. We demonstrate this effect first in its purest form in an interacting altermagnetic fermion model, in which a many-body chiral symmetry forbids any charge polarization. Our quantum Monte Carlo simulations reveal edge-localized, reversible spin accumulations fully consistent with this symmetry locking. Breaking the chiral symmetry releases the charge sector: a ferroelectric polarization emerges orthogonal to the ferro-spinetic one, yielding mutually perpendicular switchable spin- and charge-polarized responses. We identify Mn-based metal-organic frameworks as realistic hosts for this effect, offering a practical route for experimental verification.

We investigate non-stabilizerness, also known as ``magic,'' to understand criticality and exceptional points in non-Hermitian quantum many-body systems. Our focus is on parity-time ($\mathcal{PT}$) symmetric spin chains, specifically the non-Hermitian transverse-field Ising and XX models. We calculate stabilizer R\'enyi entropies in their ground states using non-Hermitian matrix product state methods. Our findings show that magic exhibits unique and model-specific signs of phase transitions. In the Ising chain, it peaks along the regular Hermitian-like critical line but disappears across exceptional points. In contrast, in the XX chain, it reaches its maximum at the exceptional line where $\mathcal{PT}$ symmetry is broken. Finite-size scaling reveals that these effects become more pronounced with larger systems, highlighting non-stabilizerness as a sensitive marker for both quantum criticality and non-Hermitian spectral degeneracies. We also investigate magic in momentum space for the XX model analytically and find that is reaches a minimum around exceptional points. Our results indicate that magic takes extremal values at the exceptional points and serves as a valuable tool for examining complexity, criticality, and symmetry breaking in non-Hermitian quantum matter.

A field-induced quantum spin liquid (QSL) state is an extraordinary phenomenon, hitherto unobserved in metallic frustrated compounds. Recent bulk measurements have revealed intriguing field-induced magnetic states in metallic frustrated CePdAl. However, the nature of these field-induced states, potentially including a QSL state, remains unclear due to the lack of detailed microscopic investigation. To elucidate these field-induced states, we employed the transverse-field muon spin relaxation/rotation (TF-$\mu$SR) technique, applying various magnetic fields parallel to the c-axis in single-crystalline CePdAl over a broad temperature range (100~K-100~mK). Our $\mu$SR data indicate that field-induced low-temperature states for fields B$\leq B_{c2}(=3.4~T)$ exhibit long-range magnetic order, whereas for B>$B_{c2}$ they yield contrasting behavior. Notably, at 3.75 T, the transverse relaxation rate ($\lambda_T$) diverges following a power-law dependence below 800~mK along with an indication of finite frustration, whereas the Knight shift is temperature independent. These observations corroborate the signature of a critical spin-liquid (CSL) with antiferromagnetic spin fluctuations. Furthermore, at 4.3 T, a non-Fermi liquid state is observed where frustration is absent. This comprehensive microscopic study strongly suggests the existence of a CSL state in a metallic frustrated system.

Resonant Inelastic X-ray scattering (RIXS) is a synchrotron-based spectroscopy that has seen growing interest across a range of scientific disciplines beyond fundamental physics. The interpretation of experimental RIXS data requires theoretical calculations based on the Kramers-Heisenberg formula. However, due to the dependence of RIXS on both the incident and scattered photon properties, a tractable treatment of the angular dependence in this formula has been lacking. In this work, within the electric dipole approximation, we determine the number of fundamental spectra contributing to the RIXS cross-section for all crystallographic point groups. We then derive a general expression for the RIXS cross-section of isotropic samples such as un-textured powders, homogeneous glasses or liquids, explicitly accounting for the polarization and propagation directions of both the incident and scattered photons. Simplified forms of the RIXS expressions are subsequently obtained for most common point groups. Finally, we demonstrate the applicability of our formalism through a case study of uranium 3d4f RIXS.

The magnetic ground state of RuO$_2$ has been under intense debate. Using first-principles calculations, we show that compressive strain along [001] direction stabilizes an altermagnetic phase in RuO$_2$ thin films grown on (100) and (110) TiO$_2$ substrates. We further identify that compressive strain enhances the density of states near the Fermi level, resulting in a Fermi surface instability and the emergence of altermagnetism. The magnitude of strain and the associated increase in the density of states can be tuned by varying the film thickness, as systematically confirmed by x-ray diffraction and photoemission spectroscopy measurements. Symmetry analysis further reveals that (100) RuO$_2$ hosts an ideal altermagnetic order, whereas broken symmetry in (110) films leads to an uncompensated ferrimagnetic state. Finally, we discuss the effects of Hubbard $U$ parameters and evaluate the realistic tunneling magnetoresistance of (100) RuO$_2$.

We investigate the nonequilibrium dynamics of charge density waves in a pumped one-dimensional Hubbard superlattice with staggered onsite Coulomb interactions at half-filling, using time-dependent exact diagonalization. In equilibrium, the system exhibits commensurate charge correlations consistent with the superlattice periodicity. Under laser excitation, the charge correlation function exhibits distinct behaviors across four representative frequencies, spanning both linear and nonlinear optical regimes. Notably, we observe a laser-induced commensurate-to-incommensurate transition in the charge order, manifested by a shift in the peak wavevector of the charge structure factor. This transition is driven by sublattice-selective doublon-holon dynamics, where the laser frequency and intensity determine whether excitations predominantly destabilize the charge order on the weakly or strongly interacting sublattice. Our analysis of the excitation spectrum and site-resolved correlation dynamics reveals the underlying mechanisms of this transition. These results suggest a promising optical strategy for controlling charge order in superlattice-based quantum materials.

Altermagnets have recently garnered significant interest due to their vanishing net magnetic moment and non-relativistic momentum-dependent spin splitting. However, altermagnetic (AM) multiferroics especially triferroics remain scarce. We investigate the experimentally synthesized non-van der Waals CrSb as a model system to explore the effects of dimensionality and facet orientation on its ferroic properties. NiAs, MnP, wurtzite (WZ), zincblende (ZB), and rocksalt (RS) phases are considered. Using first-principles calculations, we predict the altermagnetism of CrSb in MnP phase which has comparable stability with experimental NiAs phase. Both NiAs- and MnP-phase (110) facets exhibit AM-ferroelastic (FC) biferroics, while the WZ-phase bulk and (001) facets host ferromagnetic (FM) or AM-ferroelectric (FE) biferroics. Notably, the WZ-phase (110) facets are identified as FM/AM-FE-FC triferroics, with moderate energy barriers of 0.129 and 0.363 eV atom-1 for FE and FC switching, respectively. Both FE and FC switching can reverse the AM spin splitting in antiferromagnetic (AFM) configurations while preserving the high spin polarization in FM states. The magnetic anisotropy is highly tunable, exhibiting either uniaxial or in-plane behavior depending on the phase, dimension, and facet. This work establishes a framework for designing AM multiferroics through polymorphic, dimensional, and facet engineering, offering promising avenues for multifunctional spintronic applications.

The last two decades have witnessed a tremendous number of computational predictions of hydride-based (phonon-mediated) superconductors, mostly at extremely high pressures, i.e., hundreds of GPa. These discoveries were heavily driven by Migdal-\'{E}liashberg theory (and its first-principles computational implementations) for electron-phonon interactions, the key concept of phonon-mediated superconductivity. Dozens of predictions were experimentally synthesized and characterized, triggering not only enormous excitement in the community but also some debates. In this Article, we review the computational-driven discoveries and the recent developments in the field from various essential aspects, including the theoretical, computational, and, specifically, artificial intelligence (AI)/machine learning (ML) based approaches emerging within the paradigm of materials informatics. While challenges and critical gaps can be found in all of these approaches, AI/ML efforts specifically remain in its infant stage for good reasons. However, opportunities exist when these approaches can be further developed and integrated in concerted efforts, in which AI/ML approaches could play more important roles.

We present a comprehensive theoretical framework for quantum criticality in the non-Hermitian detuned PXP model, and establish the complete phase diagram, which had remained elusive in previous studies. Starting from a numerically identified phase transition point, we construct an exact second-order phase transition boundary through a similarity transformation in the real-energy regime. By introducing the biorthogonal entanglement entropy and biorthogonal Loschmidt echo, we demonstrate from both equilibrium and nonequilibrium perspectives that this transition belongs to the Ising universality class. Using the correlation function, we further distinguish between confined and deconfined phases within the $\mathcal{PT}$-symmetric region. In the complex-energy regime, we identify both a full $\mathcal{PT}$ transition and a first-excited-state $\mathcal{PT}$ transition, respectively. Moreover, we identify the location of the Yang-Lee edge singularity (YLES) using both the associated-biorthogonal and self-normal Loschmidt echoes, and extract the corresponding critical exponent, which agrees with the predictions of non-unitary conformal field theory. Finally, we propose an experimental scheme to observe the YLES in Rydberg atomic arrays, which offers a promising route to exploring non-Hermitian critical phenomena and singularities in future experimental settings.

We introduce UEIPNet, an equivariant graph neural network designed to predict both interatomic potentials and tight-binding (TB) Hamiltonians for an atomic structure. The UEIPNet is trained using density functional theory calculations followed by Wannier projection to predict energies and forces as node-level targets and Wannier-projected TB matrices as edge-level targets. This enables physically consistent modeling of coupled mechanical electronic responses with near-DFT accuracy. Trained on bilayer graphene and monolayer MoS2 DFT data, UEIPNet captures key deformation-electronic effects: in twisted bilayer graphene, it reveals how interlayer spacing, in-plane strain, and out-of-plane corrugation drive isolated flat-band formation, and further shows that modulating substrate interaction strength can generate flat bands even away from the magic angle. For monolayer MoS2, the UEIPNet accurately reproduces phonon dispersions, strain-dependent band-gap evolution, and local density of states modulations under non-uniform strain. The UEIPNet offers a generalized, efficient, and scalable framework for studying deformation-electronic coupling in large-scale atomistic systems, bridging classical atomistic simulations and electronic-structure calculations.

Materials with bespoke properties have long been identified by computational searches, and their experimental realisation is now coming within reach through autonomous laboratories. Scattering experiments are central to verifying the atomic structures of autonomously synthesised materials. Yet, interpreting these measurements typically requires user expertise and manual processing, or machine learning (ML) models trained on predefined datasets, limiting fully autonomous materials discovery. Here, we introduce a differentiable optimisation framework that treats scattering calculations, energetics, and chemical constraints as a unified refinement problem. Capability demonstrations across molecules, crystal structures, nanoparticles, and amorphous matter show that this data-driven approach resolves structural degeneracies with multi-modal inputs - suggesting its usefulness for informing, and ultimately guiding, the operation of autonomous laboratories.

Multiferroics, which combine ferroelectric and magnetic order, offer a transformative platform for next-generation electronic devices. However, the intrinsic competition between the mechanisms driving ferroelectricity and magnetism in single-phase materials severely limits their performance, typically resulting in weak magnetoelectric coupling at room temperature. Here, we propose a solution to this long-standing challenge through the novel concept of fractional quantum multiferroics (FQMF), where strong magnetoelectric coupling is naturally realized by coupling fractional quantum ferroelectricity (FQFE) with altermagnetism (AM). Symmetry analysis shows that reversing the FQFE polarization necessarily inverts the AM spin splitting under parity-time ($\mathcal{PT}$) or time-reversal ($\mathcal{T}\tau$) operations. A minimal tight-binding model reproduces this effect, demonstrating electrically driven spin control without rotating the N\'eel vector. First-principles calculations further identify a broad family of candidate materials in two and three dimensions including bulk MnTe, Cr$_2$S$_3$, Mn$_4$Bi$_3$NO$_{15}$ and two-dimensional AB$_2$ bilayers such as MnX$_2$ (X=Cl, Br, I), CoCl$_2$, CoBr$_2$, and FeI$_2$. Notably, MnTe exhibits a high N\'eel temperature ($\sim$300 K) and a large electrically switchable spin splitting ($\sim$0.8 eV), demonstrating room-temperature magnetoelectric performance that surpasses that of conventional multiferroics. To further showcase the technological potential, we propose an electric-field-controlled FQMF tunnel junction based on MnTe that achieves tunneling magnetoresistance exceeding 300\%. This work establishes FQMF as a distinct and promising route to achieving room-temperature strong magnetoelectric coupling, opening a new avenue for voltage-controlled spintronics.