arXiv Curator

By combining density functional theory (DFT) and density matrix renormalization group calculations, we investigate the unusual pressure dependence of superconducting transition temperature ($T_c$) in the nickelate superconductor La$_{3}$Ni$_{2}$O$_{7}$. Using the hopping integrals and on-site potentials obtained by fitting the DFT band structures, we map a quantum phase diagram of a bilayer two-orbital Hubbard model with increasing pressure in a ladder geometry, which has an intermediate Hubbard repulsion and a Hund's coupling. Near $3/8$ filling, we find a strong spin density wave order, which at $3/8$ filling shows a real-space spin pattern similar to the spin-charge stripe order along a lattice direction. At $21/64$ filling, we find a superconducting phase with interlayer superconductivity (SC) in both the $d_{z^2}$ and $d_{x^2-y^2}$ orbitals, as well as in-plane SC in the $d_{z^2}$ orbital. Intriguingly, the SC is weakened with increasing pressure and transits to a Luttinger liquid above $80$ GPa, which qualitatively agrees with the experimental observations of decreasing $T_c$ with increasing pressure and a transition to Fermi liquid above $80$ GPa in La$_{3}$Ni$_{2}$O$_{7}$. Through a comparative study, we further show that the ratio of interaction to hopping integral, which reduces moderately with increasing pressure, may play a dominant role in the weakening of SC. Our results of this experimentally relevant model not only find a robust SC through suppressing the competing spin density wave order, but also give new insight into the unusual pressure dependence of SC in La$_{3}$Ni$_{2}$O$_{7}$.

The precise crystal structure of La3Ni2O7 in its high-pressure superconducting state has been a subject of intense debate, with proposed models including both orthorhombic and tetragonal symmetries. Using high-pressure Raman spectroscopy combined with frst-principles calculations, we unravel the structural evolution of La3Ni2O7 under pressure up to 32.7 GPa. We identify a clear structural transition sequence: from the orthorhombic Amam phase to a mixed Amam+Fmmm phase at 4 GPa, followed by a complete transition to the tetragonal I4/mmm phase at 14.5 GPa, which is signaled by a pronounced phonon renormalization. The emergence of bulk superconductivity is found to coincide precisely with this transition to the I4/mmm phase. Our results de nitively establish the tetragonal I4/mmm structure as the host of superconductivity in La3Ni2O7, resolving a central controversy and providing a critical foundation for understanding the superconducting mechanism in nickelates.

We uncover a dual mechanism by which oxygen defects suppress superconductivity in the bilayer nickelate La$_3$Ni$_2$O$_{7+δ}$ using density functional theory, dynamical mean-field theory, and functional renormalization group analysis. Apical vacancies and interbilayer interstitials emerge as the dominant low-energy defect species and are further stabilized by orthorhombic domain walls. These two defect classes drive the electronic structure in opposing directions. Vacancy-induced disorder generates local magnetic moments and promotes Anderson localization at moderate concentrations, whereas periodic interstitial ordering yields a coherent but weakly correlated metallic background that fails to support superconductivity. These findings highlight the decisive role of oxygen defects in shaping the superconducting and provide microscopic guidance for improving superconductivity through controlled defect engineering.

Recently, a preliminary magnetic field-reentrant superconductivity manifested in high critical-temperature (Tc) Eu-doped infinite-layer (IL) nickelates, beyond analogous discoveries exclusively in low-Tc systems. This arises more intriguing fundamental issues about potential quantum-phase boundary and criticality between unconventional superconductivity and field-reentrant-one, which are inexplicable owing to formidable challenges in growing IL-nickelates towards later-series rare-earths. Herein, we demonstrate the 4f-orbital related quantum rivalries between high-Tc and reentrant superconductivity in (Nd1-yRE'y)1-xEuxNiO2 (RE': Pr, Nd, Sm, Gd and Dy) system, via opening up the chemical avenue for a quantum leap in their growths. Robust magnetic field-reentrant superconductivity with uniaxial anisotropy is validated to fringe at boundary of the superconducting dome with optimal-Tc near 40 K. Reinforced reentry is realized via introducing Gd3+ with half-filled deeper 4f-orbital energy-states, compared to Eu2+, that strengthens on-site magnetic-moment. Our findings largely enrich the superconducting phase-diagram for nickelates, establishing an ideal platform for studying 4f-related unconventional superconductivity and quantum criticality.

We introduce a neural network impurity solver for real-frequency DMFT that employs a multihead cross-attention mechanism to map hybridization functions to spectral functions, conditioned on impurity parameters. Trained on high-quality MPS data from complex contour time evolution and incorporating derivative constraints with respect to the complex-time angle, our model achieves smooth generalization to the real-frequency axis. Benchmarking on the single-band Hubbard model for the Bethe lattice demonstrates quantitative accuracy across metallic, strongly correlated, and insulating regimes.

Recent studies suggest that the tetragonal phase of the Ruddlesden-Popper (RP) bilayer nickelate, La$_3$Ni$_2$O$_7$ or La$_2$PrNi$_2$O$_7$, which is stabilized under high pressures, is responsible for high-temperature superconductivity (HTSC). In this context, realization of the tetragonal phase at ambient pressure could be a rational step to achieve the goal of ambient-pressure HTSC in the nickelate system. By employing the concept of Goldschmidt tolerance factor, we succeed in stabilizing the tetragonal phase by aluminum doping together with post annealing under moderately high oxygen pressure. X-ray and neutron diffractions verify the tetragonal $I4/mmm$ structure for the post-annealed samples La$_3$Ni$_{2-x}$Al$_x$O$_{7-δ}$ (0.3 $\leq x \leq$ 0.5). The Al-doped samples, including the tetragonal ones, show semiconducting properties, carry localized magnetic moments, and exhibit spin-glass-like behaviors at low temperatures, all of which can be explained in terms of charge carrier localization. Furthermore, high-pressure resistance measurements on post-annealed samples reveal that even a low Al doping ($x$ = 0.05) suppresses superconductivity almost completely. This work gives information about the effect of nonmagnetic impurity on metallicity as well as superconductivity in bilayer nickelates, which would contribute to understanding the superconducting mechanism in RP nickelates.

To this day, high-temperature cuprate superconductors remain an unparalleled platform for studying the competition and coexistence of emergent, static and dynamic, quantum phases of matter exhibiting high transition temperature non-s-wave superconductivity, non-Fermi liquid transport and a still enigmatic pseudogap regime. However, how superconductivity emerges alongside and competes with the pseudogap regime remains an open question. Here, we present a high-resolution, time- and angle-resolved photoemission study of the near-antinodal region of optimally-doped Bi$_2$Sr$_2$CaCu$_2$O$_{8+δ}$. For a sufficiently high excitation fluence, we disrupt superconductivity and drive a transient change from a symmetric superconducting-like to an asymmetric pseudogap-like density of states, for electronic temperatures well below the equilibrium superconducting critical temperature. Conversely, when the superconductivity is fully restored, the pseudogap is suppressed, as signaled by a fully particle-hole symmetric density of states. A unique aspect of our experiments is that the pseudogap coexists with superconducting features at intermediate times or at intermediate fluence. Our findings challenge the paradigm that superconductivity emerges by establishing phase coherence in the pseudogap. Instead, our experimental results, supported by phenomenological theory, demonstrate that the two states compete, and that the low-temperature ground state of the cuprates originates from a competition between superconducting and pseudogap states.

The Fermi-Hubbard model is the starting point for the simulation of many strongly correlated materials, including high-temperature superconductors, whose modelling is a key motivation for the construction of quantum simulation and computing devices. However, the detection of superconducting pairing correlations has so far remained out of reach, both because of their off-diagonal character-which makes them inaccessible to local density measurements-and because of the difficulty of preparing superconducting states. Here, we report measurement of significant pairing correlations in three different regimes of Fermi-Hubbard models simulated on Quantinuum\'s Helios trapped-ion quantum computer. Specifically, we measure non-equilibrium pairing induced by an electromagnetic field in the half-filled square lattice model, d-wave pairing in an approximate ground state of the checkerboard Hubbard model at $1/6$-doping, and s-wave pairing in a bilayer model relevant to nickelate superconductors. These results show that a quantum computer can reliably create and probe physically relevant states with superconducting pairing correlations, opening a path to the exploration of superconductivity with quantum computers.

Inspired by nature, this study employs the Materials Genome Initiative to identify key components of HTSC superconductors. Integrating AI with high-throughput screening, we uncover crucial superconducting "genes". Through HTS techniques and advanced machine learning models, we demonstrate that Functional Convolutional Neural Networks (CNNs) ensure accurate extrapolation of potential compounds. Leveraging extensive datasets from the ICSD, the Materials Project and COD, our implemented HTS pipeline classifies superconductors, with CNN and long short-term memory (LSTM) models predicting Tc and their foundational elements. We address the scarcity of non-superconducting material data by compiling a dataset of 53,196 non-superconducting materials (DataG Non-Sc) and introduce a novel neural network architecture using Functional API for improved prediction, offering a powerful tool for future superconductor discovery. Our findings underscore the transformative potential of combining HTS with AI-driven models in advancing HTSC materials, highlighting Pu and H elements (with Tc nearly 100 K) as significant predictors of high-temperature superconductivity, suggesting their role as a crucial gene in these materials.

Large Language Models (LLMs) show great promise as a powerful tool for scientific literature exploration. However, their effectiveness in providing scientifically accurate and comprehensive answers to complex questions within specialized domains remains an active area of research. Using the field of high-temperature cuprates as an exemplar, we evaluate the ability of LLM systems to understand the literature at the level of an expert. We construct an expert-curated database of 1,726 scientific papers that covers the history of the field, and a set of 67 expert-formulated questions that probe deep understanding of the literature. We then evaluate six different LLM-based systems for answering these questions, including both commercially available closed models and a custom retrieval-augmented generation (RAG) system capable of retrieving images alongside text. Experts then evaluate the answers of these systems against a rubric that assesses balanced perspectives, factual comprehensiveness, succinctness, and evidentiary support. Among the six systems two using RAG on curated literature outperformed existing closed models across key metrics, particularly in providing comprehensive and well-supported answers. We discuss promising aspects of LLM performances as well as critical short-comings of all the models. The set of expert-formulated questions and the rubric will be valuable for assessing expert level performance of LLM based reasoning systems.

Infinite-layer nickelates are among the most promising cuprate-akin superconductors, although relevant differences from copper oxides have been reported. Here, we present momentum- and polarization-resolved RIXS measurements on chemically undoped, superconducting PrNiO2, and compare its magnetic and orbital excitations with those of the reference infinite layer cuprate CaCuO2. In PrNiO2, the in-plane magnetic exchange integrals are smaller than in CaCuO2, whereas the out-of-plane values are similar, indicating that both materials support a three-dimensional antiferromagnetic order. Orbital excitations, associated to the transitions within 3d states of the metal, are well reproduced within a single-ion model and display similar characteristics, except for the Ni-dxy peak which, besides lying at significantly lower energy, shows an opposite dispersion to that of Cu-dxy. This is interpreted as a consequence of orbital superexchange coupling between nearest neighbor sites, which drives the orbiton propagation. Our observations demonstrate that infinitelayer cuprates and nickelates share most of the spin and orbital properties, despite their markedly different charge-transfer energy Delta.

Infinite-layer nickelates are among the most promising cuprate-akin superconductors, although relevant differences from copper oxides have been reported. Here, we present momentum- and polarization-resolved RIXS measurements on chemically undoped, superconducting PrNiO2, and compare its magnetic and orbital excitations with those of the reference infinite layer cuprate CaCuO2. In PrNiO2, the in-plane magnetic exchange integrals are smaller than in CaCuO2, whereas the out-of-plane values are similar, indicating that both materials support a three-dimensional antiferromagnetic order. Orbital excitations, associated to the transitions within 3d states of the metal, are well reproduced within a single-ion model and display similar characteristics, except for the Ni-dxy peak which, besides lying at significantly lower energy, shows an opposite dispersion to that of Cu-dxy. This is interpreted as a consequence of orbital superexchange coupling between nearest neighbor sites, which drives the orbiton propagation. Our observations demonstrate that infinitelayer cuprates and nickelates share most of the spin and orbital properties, despite their markedly different charge-transfer energy Delta.

Infinite-layer nickelates are among the most promising cuprate-akin superconductors, although relevant differences from copper oxides have been reported. Here, we present momentum- and polarization-resolved RIXS measurements on chemically undoped, superconducting PrNiO2, and compare its magnetic and orbital excitations with those of the reference infinite layer cuprate CaCuO2. In PrNiO2, the in-plane magnetic exchange integrals are smaller than in CaCuO2, whereas the out-of-plane values are similar, indicating that both materials support a three-dimensional antiferromagnetic order. Orbital excitations, associated to the transitions within 3d states of the metal, are well reproduced within a single-ion model and display similar characteristics, except for the Ni-dxy peak which, besides lying at significantly lower energy, shows an opposite dispersion to that of Cu-dxy. This is interpreted as a consequence of orbital superexchange coupling between nearest neighbor sites, which drives the orbiton propagation. Our observations demonstrate that infinite layer cuprates and nickelates share most of the spin and orbital properties, despite their markedly different charge-transfer energy Delta.

Optimization of unconventional superconductivity involves a balance of interaction strengths. Precise determination of correlation strength across different material families is therefore important. Here, we present a combined X-ray absorption spectroscopy (XAS) and resonant inelastic X-ray scattering (RIXS) study of infinite-layer PrNiO$_2$ and SrCuO$_2$ that enables fair comparison of their interaction strengths. For both compounds, we study the orbital and magnetic excitations and extract their dispersions along high-symmetry directions. Using a single-band Hubbard model and including higher-order exchange interactions, we derive the correlation factor $U/t$ for both compounds. A key finding is that despite a smaller Coulomb repulsion $U$, PrNiO$_2$ exhibits a correlation strength that is 20% stronger than that of its isostructural cuprate counterpart SrCuO$_2$. This indicates that a moderation of the correlation strength may further optimize superconductivity in nickelates.

We propose a structure-property relation that could be key to the pseudogap phenomenology in cuprates. The underlying nonsymmorphic crystal structure in the low-temperature orthogonal phase endows the lattice with a sublattice structure that gives rise to two electronic bands near the Fermi surface. In the presence of spin-orbit coupling, the hybridization of these two bands generates small Fermi pockets and, correspondingly, a change in the number of free carriers. A sublattice structure also leads to angle-resolved photoemission spectroscopy (ARPES) matrix-element interference, naturally explaining the emergence of Fermi arcs and their consistency with closed Fermi pockets. We employ a symmetry analysis to highlight the expected Fermi surface properties, and complement it with density functional theory (DFT) for a quantitative discussion and comparison with recent experiments in doped La$_2$CuO$_4$. The proposed mechanism can consistently account for the most salient features of the pseudogap in the cuprates, namely, the Fermi surface reconstruction with the formation of small Fermi pockets and the corresponding change in carrier density, and the observation of Fermi arcs by ARPES.

Whether the ground state of the square lattice Hubbard model exhibits superconductivity remains a major open question, central to understanding high temperature cuprate superconductors and ultra-cold fermions in optical lattices. Numerical studies have found evidence for stripe-ordered states and superconductivity at strong coupling but the phase diagram remains controversial. Here, we show that one can resolve the subtle energetics of metallic, superconducting, and stripe phases using a new class of neural quantum state (NQS) wavefunctions that extend hidden fermion determinant states to Pfaffians. We simulate several hundred electrons using fast Pfaffian algorithms allowing us to measure off-diagonal long range order. At strong coupling and low hole-doping, we find that a non-superconducting filled stripe phase prevails, while superconductivity coexisting with partially-filled stripes is stabilized by a negative next neighbor hopping t-prime, with |t-prime| > 0.1. At larger doping levels, we introduce momentum-space correlation functions to mitigate finite size effects that arise from weakly-bound pairs. These provide evidence for uniform d-wave superconductivity at U = 4, even when t-prime = 0. Our results highlight the potential of NQS approaches, and provide a fresh perspective on superconductivity in the square lattice Hubbard model.

The frequency structure of the superconducting correlations in cuprates gives insights on the pairing mechanism. Here we present an exhaustive study of this problem in the two-dimensional Hubbard model with cellular dynamical mean-field theory. To this end, we systematically quantify the dependence on doping and interaction strength of the superconducting gap, of the frequency scales where pairing occurs, and of their relative contribution to pairing. We find pair-forming alternating with pair-breaking processes as a function of frequency, providing new evidence that pairing can arise in principle from both low- and high-frequency processes. However, we find that the high-frequency pair-forming processes are outweighed by pair-breaking ones. Hence, the net contribution to pairing comes only from the low-frequency pair-forming processes, a result that underscores their key role in the pairing mechanism.

We elucidate the electronic structure and quantum many-body instabilities of the monolayer nickelate La2NiO4 under hydrostatic pressure using a combination of density functional theory, dynamical mean-field theory (DFT+DMFT), and random phase approximation (RPA). Our DFT+DMFT calculations reveal non-Fermi-liquid behavior and coherence loss near the Fermi level at low pressures, driven by strong electron correlations within the Ni-e_g orbital manifold, which is analogous to the low-energy electronic properties observed in La3Ni2O7. However, multi-orbital spin susceptibility analysis demonstrates an exceptionally suppressed critical Stoner parameter U_c (about 0.4~0.7 eV), indicating robust magnetic order that dominates the ground state and precludes superconductivity in the pristine system. Below U_c, superconducting instabilities exhibit a pressure-driven symmetry transition: the d_(x2-y2)-wave pairing prevails at ambient and low pressure, while the s+g-wave symmetry occurs above 75 GPa. This transition is attributed to pressure-induced self-doping effect. The high-angular-momentum g-wave component incurs significant energetic penalties, rendering high-Tc superconductivity unlikely. We conclude that the absence of superconductivity in La2NiO4 arises from its robust intrinsic magnetism and the unfavorable pairing symmetry under pressure, suggesting that alternative routes-such as chemical doping or epitaxial strain-are necessary to suppress magnetism and unlock superconducting states.

The intricate interplay of structural, charge and spin orders in layered cuprates leads to emergent phenomena, most notably high-temperature superconductivity. However, there is growing awareness that both the structure and electronic ordering that underpin them are not fully homogeneous. Here, we employ scanning three-dimensional X-ray diffraction to spatially resolve structural distortions in La$_{1.675}$Eu$_{0.2}$Sr$_{0.125}$CuO$_{4}$, a prototypical system that exhibits a strong competition between superconductivity and charge density wave order, across its low-temperature orthorhombic to tetragonal phase transition. We uncover two forms of intrinsic microstructural heterogeneity: at 300 K, we reveal remarkably wide tetragonal-like domain wall regions within the nominally orthorhombic crystal structure, and, upon cooling to 100 K, a fine microstructure of orthorhombic-like stripes embedded within the tetragonal matrix emerges. This spatially resolved view directly defines the microstructural architecture of phase coexistence in this system, demonstrating how structural distortions generate intrinsic heterogeneity that shapes the balance between superconductivity and charge density wave order as constrained by the restrictive orthorhombic and tetragonal structures, offering a pathway to control correlated phenomena.

Many emergent phenomena appear in doped Mott insulators near the insulator-to-metal transition. In high-temperature cuprate superconductors, superconductivity arises when antiferromagnetic (AFM) order is gradually suppressed by carrier doping, and a $\textit{d}$-wave superconducting gap forms when an enigmatic nodal gap evolves into a point node. Here, we examine electron-doped Sr$_{2}$IrO$_{4}$, the 5$\textit{d}$-electron counterpart of cuprates, using angle-resolved photoemission spectroscopy. At low doping levels, we observe the formation of electronic states near the Fermi level, accompanied by a gap at the AFM zone boundary, mimicking the AFM gap in electron-doped cuprates. With increasing doping, a distinct gap emerges along the (0,0)-($\pi$,$\pi$) nodal direction, paralleling that observed in hole-doped cuprates. This anomalous nodal gap persists after the collapse of the AFM gap and gradually decreases with further doping. It eventually vanishes into a point node of the reported $\textit{d}$-wave gap. These observations replicate the characteristic features in both electron- and hole-doped cuprates, indicating a unified route toward nodal metallicity in doped spin-1/2 AFM Mott insulators.

Angle-resolved photoemission spectroscopy in combination with band structure calculations have shown that the layered transition metal dichalcogenide NiTe$_2$ is a type-II Dirac semimetal. However, there are conflicting conclusions in the literature regarding the role of electron correlations in NiTe$_2$. We study the core-level and valence band electronic structure of single crystal NiTe$_2$ using soft and hard X-ray photoemission spectroscopy (SXPES, HAXPES), X-ray absorption spectroscopy (XAS) and Ni $2p-3d$ resonant photoemission spectroscopy(resonant-PES) to quantify electronic parameters in NiTe$_2$. The on-site Coulomb energy ($U_{dd}$) in the Ni $3d$ states is quantified from measurements of the Ni $3d$ single particle density of states and the two-hole correlation satellite. The Ni $2p$ core level and $L$-edge XAS spectra are analyzed by charge transfer cluster model calculations using the experimentally estimated $U_{dd}$ (= 3.7 eV), and the results show that NiTe$_2$ exhibits a negative charge-transfer energy ($\Delta$ = -2.8 eV). The same type of cluster model analysis of NiO $L$-edge XAS confirms its well-known strongly correlated charge-transfer insulator character, with $U_{dd}$ = 7.0 eV and $\Delta$ = 6.0 eV. The $d$-$p$ hybridization strength $T_{eg}$ for NiTe$_2$$<$NiO, and indicates that the reduced $U_{dd}$ in NiTe\textsubscript{2} compared to NiO is not due to an increase in $T_{eg}$. The increase in $d^n$ count on the Ni site in NiTe$_{2}$ by nearly one electron is attributed to negative-$\Delta$ and a reduced $U_{dd}$. However, since $U_{dd}$$>$$|\Delta|$, the results indicate the important role of a finite repulsive $U_{dd}$ in making NiTe$_{2}$ a moderately correlated $p$-type Dirac semi-metal.

Rare-earth nickelates RNiO3 (R=rare-earth element) exhibit three kinds of phase transitions with decreasing temperature: a structural transition from a pseudo-cubic to a monoclinic phase, a metal- insulator transition (MIT), and a magnetic transition from a paramagnetic state to an ordered one. The first two occur at the same temperature, which has led to a consensus that the MIT is driven by lattice distortions. We show here that the primary driving force for the MIT is magnetic; however because of the unusual d7 configuration of Ni, additional flexibility in spin configurations are also needed which symmetry-lowing structural deformations make possible. The latter enable Ni to disproportionate into two kinds: a high-spin and a low-spin configuration, which allow the system to reduce its unfavorable orbital moment and also open a gap.

The Hubbard model has attracted considerable interest due to its prototypical role in describing strongly interacting electronic systems, such as high-critical-temperature superconductors as well as many novel quantum materials. By introducing next-nearest-neighbor (NNN) hoppings to the Hubbard model, the phase diagram becomes richer, and fascinating phenomena arise in both, one-dimensional chains and square lattices, such as: antiferromagnetism (AFM), ferromagnetism (FM), superconductivity (SC), as well as charge orders, among others. Moreover, NNN hoppings play a fundamental role in understanding effects of doping on magnetism and pairing orders in strongly interacting regimes. In this article, we review the recent progress in understanding the different competing phases of this model in one and two dimensions from a computational perspective. We comment on the pressing technical challenges, illustrate the controversial results concerning the emergence of the SC phase, and conclude with our perspectives on future explorations.

Characterizing the spin and charge degrees of freedom in high-temperature superconducting cuprates under non-equilibrium conditions provides new insights into their electronic correlations. However, their collective dynamics have been largely unexplored due to experimental challenges. Here, we use time-resolved resonant inelastic X-ray scattering (trRIXS) at the Cu $L_3$-edge to simultaneously track the collective spin (paramagnon) and charge (acoustic plasmon) dynamics in an optimally electron-doped cuprate driven out-of-equilibrium by a femtosecond pump laser pulse. Upon pumping, we observed an anti-Stokes signal associated with paramagnon generation, which modifies the paramagnon dispersion near the zone center, though the bandwidth remained unchanged, suggesting no significant alteration to spin exchange interactions. Simultaneously, in the charge sector, the acoustic plasmon's energy and spectral weight decreased, suggesting a light-induced redistribution of charge carriers. The variations of both the paramagnon and the plasmon were locked in time, demonstrating a robust intertwining between the spin and charge degrees of freedom on a femtosecond timescale, even in this non-equilibrium state.

We investigated the Cu-spin correlation in the overdoped regime of the electron-doped high-Tc cuprate thin films of La_2-x_Ce_x_CuO_4_, changing the reduction condition from muon spin relaxation using low-energy muons. The Cu-spin correlation developed at low temperatures for optimally reduced films with x=0.13 as well as x=0.17 where the superconductivity was almost suppressed. These results are contrary to those observed in the hole-doped high-Tc cuprates where the development of the antiferromagnetic Cu-spin correlation disappears together with the suppression of superconductivity. The Cu-spin correlation developed at low temperatures in x=0.17 may be understood in terms of antiferromagnetism, but it may be related to a ferromagnetic order recently suggested in the nonsuperconducting heavily overdoped La$_2-x_Ce_x_CuO_4_ with x~0.18.

The subject of the present study is the double-layer square-lattice system with the intralayer phonon modulations. We investigate the superconducting and excitonic pairings, as well as their coexistence, as functions of various physical parameters in the system. These parameters include temperature, intralayer and interlayer Coulomb interactions, the electron-phonon coupling parameter, doping and the applied electric field. The existence of superconductivity is demonstrated by considering the influence of intralayer phonons on the total charge density leading to the modification of the total energy of the electrons. Our results provides insights into the long-standing problem of the mechanism of superconductivity in high-$T_c$ cuprate superconductors.

To clarify the origin of the pseudogap and strange metal states as well as their mutual relationship in cuprate superconductors, a comprehensive study on the spectral function, Fermi surface, resistivity and dynamical spin susceptivity of the Hubbard model on the square lattice has been conducted by means of the ladder dual-fermion approximation with an electron self-energy correction. It is found that the appearance of these two states requires that the characteristic hole concentration below which the Mott-Heisenberg and Slater mechanisms of electron localization occurs $p_{\rm MS}$ nearly coincides with the hole concentration where the Van Hove singularity (VHS) point is placed at the vicinity of the Fermi level. When this condition is met the VHS point is pined at which the nesting condition of the antiferromagnetic (AFM) fluctuation is fulfilled almost everywhere on the Fermi surface in wide range of the hole concentration in a metallic state, i.e., the strange metal state. The spin fluctuation of the strange metal state is nearly quantum critical and the dynamical spin susceptivity is well described by overdamped spin wave having the $ω/T$ scaling with the relaxation rate at the Planckian limit. Because of these distinctive features of the strange metal state, the $k$ dependence of scattering rate of electrons is small and electrons behave as the marginal Fermi liquid, resulting in $T$-linear resistivity. In contrast, the pseudogap state is magnetically in the renormalized classical regime and the pseudogap is formed near the X point where the nesting condition of the short-range AFM order is fulfilled.

The tension between fermion pairing and magnetism affects numerous strongly correlated electron systems, from high-temperature cuprates to twisted bilayer graphene. Exotic forms of fermion pairing and superfluidity are predicted when attraction between fermions competes with spin doping. Here, we follow the evolution of fermion pairing and charge and spin order in a spin-imbalanced attractive Hubbard gas of fermionic $^{40}$K atoms, covering a wide range of densities, magnetizations, and interactions with single-atom resolution. At low spin imbalance and weak interactions, we find a mixture of nonlocal fermion pairs coexisting with itinerant excess fermions. For stronger interactions an effective hard-core Bose-Fermi mixture emerges. Spin doping drives a crossover from charge-density wave correlations to a Fermi liquid of polarons. Beyond a certain spin imbalance and interaction strength, we find evidence for the onset of combined spin- and pair-density wave order, a possible precursor for the existence of magnetized superfluidity in the attractive Hubbard system.

We demonstrate that a machine learning framework based on kernel ridge regression can encode and predict the self-energy of one-dimensional Hubbard models using only mean-field features such as static and dynamic Hartree-Fock quantities and first-order GW calculations. This approach is applicable across a wide range of on-site Coulomb interaction strengths $U/t$, ranging from weakly interacting systems ($U/t \ll 1$) to strong correlations ($U/t > 8$). The predicted self-energy is transformed via Dyson's equation and analytic continuation to obtain the real-frequency Green's function, which allows access to the spectral function and density of states. This method can be used for nearest-neighbor interactions $t$ and long-range hopping terms $t'$, $t''$, and $t'''$.

The issue of whether $d$-wave superconductivity (SC) occurs in the square-lattice Hubbard model with $U$ of order of the bandwidth has been one of the most debated issues to emerge from the study of high temperature SC. Here, we report variational results on eight-leg cylinders with next-nearest-neighbor hopping in the range $-0.5 t \leq t'\leq 0.25 t$ with $U = 8t$ and $12t$ and doped hole concentrations $δ=1/12$ and $1/8$. For $t'\leq 0$, the ground-state appears to be a charge-density wave (CDW) of one sort or another with SC correlations that are extremely short-ranged. In contrast, in some cases, the local magnetic order has a correlation length greater than half the cylinder width - suggestive that magnetic order might also arise in the 2D limit. For $t'>0$, our results depend more strongly on boundary conditions (periodic vs antiperiodic), making it still harder to correctly guess whether SC or CDW correlations dominate in the 2D limit. These results were obtained employing matrix-product states with bond dimensions large enough that energy differences as small as $10^{-3}t$ per site can be resolved.

Numerical simulations of strongly correlated fermions at finite temperature are essential for studying high-temperature superconductivity and other quantum many-body phenomena. The recently developed tangent-space tensor renormalization group (tanTRG) provides an efficient and accurate framework by representing thermal density operators as matrix product operators. However, the particle number generally varies during the cooling process. The conventional strategy of fine-tuning chemical potentials to reach a target filling is computationally demanding. Here we propose a fixed-$N$ tanTRG algorithm that stabilizes the average particle number by adaptively tuning the chemical potential within the imaginary-time evolution. We benchmark its accuracy on exactly solvable free fermions, and further apply it to the square-lattice Hubbard model. For hole-doped cases, we study the temperature evolution of charge and spin correlations, identifying several characteristic temperature scales for stripe formation. Our results establish fixed-$N$ tanTRG as an efficient and reliable tool for finite-temperature studies of correlated fermion systems.

We consider the SU(2) gauge theory for spin fluctuations in the two-dimensional Hubbard model, where the electron field is fractionalized in terms of spinons and chargons. In this theory, spinons are described by a non-linear sigma model, while chargons are treated as fermions at a mean-field level. We investigate the instability to a superconducting state SC*, arising from a fractionalized Fermi liquid (FL*) where pairing between chargons occurs. Consistent with previous studies, our analysis reveals a coexisting phase characterized by both magnetic and superconducting order for the chargons. The central contribution of this work is the calculation of the feedback of superconductivity on spatial and temporal spin stiffnesses, thereby quantifying its impact on spin fluctuations. Our key finding is that superconductivity significantly suppresses these spin stiffnesses, enhancing quantum spin fluctuations. This enhancement suggests that superconductivity can play a crucial role in stabilizing quantum disorder against long-range magnetic ordering.

We present a theory for the spin fluctuation spectrum of the hole-doped cuprates in a ground state with period 4 unidirectional charge density wave (`stripe') order. Motivated by recent experimental evidence for a fractionalized Fermi liquid (FL*) description of the intermediate temperature pseudogap metal, we employ a theory of fermionic spinons which are confined with the onset of stripe order at low temperatures. The theory produces the `hourglass' spectrum near stripe-ordering wavevector observed by neutron scattering. Additional scattering from spinon continua and bound states appears at higher energies and elsewhere in the Brillouin zone, and could be observed by neutron or X-ray scattering.

We study interaction-driven crossover from a correlated metallic state to Mott-insulating state on the half-filled 3x3 square cluster using exact diagonalization in the S=1/2 and S=3/2 sectors. Using complementary diagnostics such as order parameters, principal component analysis (PCA) of correlation matrices and quantum geometry, we obtain a unified, data-driven characterization of the crossover. All three diagnostics consistently identify a broad crossover centered in the weak-to-intermediate coupling regime U~2-6; charge gaps open rapidly, double occupancy suppresses, local moments form, entanglement is significantly reduced, PCA concentrates variance into a few dominant modes and the distance matrices reveals rapid wavefunction reorganization. Our multi-pronged approach establishes the Mott crossover as a smooth but well defined reorganization of electronic correlations with pronounced spin-sector dependence in finite systems.

We present a framework to characterize Mott insulating phases within the interacting one-body picture, focusing on the Hubbard diamond chain featuring both Hubbard interactions and spin-orbit coupling simulated within cellular dynamical mean field theory. Using symmetry analysis of the single-particle Green's function, we classify spectral functions by irreducible representations at high-symmetry points of the Brillouin zone. Complementarily, we calculate the one-body reduced density matrix which allows us to reach both a qualitative description of charge distribution and an analysis of the state purity. Moreover, within the Tensor Network framework, we employ a Density Matrix Renormalization Group approach to confirm the presence of three distinct phases and their corresponding phase transitions. Our results highlight how symmetry-labelled spectral functions and effective orbital analysis provide accessible single-particle tools for probing correlation-driven insulating phases.

Altermagnetism is a recently identified phase with a d, g or i-wave spin symmetry of magnetic ordering. Its discovery opens new research fronts at intersections of magnetism and spintronics with fields ranging from superconductivity to topological and relativistic quantum physics. Here we demonstrate an unconventional relativistic spin polarization in an altermagnet by spin and angle resolved photoemission spectroscopy of electronic bands in single-domain MnTe. The relativistic spin-orbit coupling origin is revealed by observing that the alternating momentum-dependent spin polarization is orthogonal to the magnetic-ordering vector. The collinearity, even-parity and time-reversal-odd nature of the demonstrated relativistic spin polarization in the altermagnet is unparalleled in conventional forms of the relativistic spin polarization. Our experimental results and methodology are supported by non-relativistic spin-symmetry and relativistic magnetic-symmetry analyses, and microscopic ab initio ground-state and photoemission theory.

Altermagnets are characterized by anisotropic band/spin splittings in momentum space, dictated by their spin-space group symmetries. However, the real-space modulations of altermagnetism are often neglected and have not been explored experimentally. Here we combine neutron diffraction, angle-resolved photoemission spectroscopy (ARPES), spin-resolved ARPES and density functional theory to demonstrate that Cs$_{1-δ}$V$_2$Te$_2$O realizes a spatially modulated form of altermagnetism, i.e., hidden altermagnetism. Such a state in Cs$_{1-δ}$V$_2$Te$_2$O results from its G-type antiferromagnetism and two-dimensional electronic states, allowing for the development of spatially alternating altermagnetic layers, whose local spin polarizations are directly verified by spin-resolved ARPES measurements. Our experimental discovery of hidden altermagnetism broadens the scope of unconventional magnetism and opens routes to exploring emergent phenomena from real-space modulations of altermagnetic order.

In this work, we revisit the electron-hole asymmetry of antiferromagnetism in cuprates by studying the three-band Emery model. Using parameters relevant to La$_2$CuO$_4$, we benchmark the anti-ferromagnetic response for a large range of dopings with variational Monte Carlo, determinant quantum Monte Carlo, constrained-path auxiliary-field quantum Monte Carlo, density-matrix embedding theory, and the Gutzwiller approximation. Across methods and accessible sizes/temperatures, we find no significant electron-hole asymmetry if we consider only Neel anti-ferronagnetic response and ignore other possible orders such as stripe state. This result is robust to a moderate oxygen-site repulsion $U_p$ and to parameter sets of Nd$_2$CuO$_4$. Incorporating dopant-induced local potentials reveals an extrinsic route to asymmetry: Cu-site defects enhance AFM on the electron-doped side, whereas O-site defects suppress it on the hole-doped side. These results indicate that dopant-driven effects make a non-negligible contribution to apparent electron-hole asymmetry in the general phase diagram of cuprates and should be included when analyzing competing orders in cuprates.

Quantum spin liquids (QSL) are phases of matter which are distinguished not by the symmetries they break, but rather by the patterns of entanglement within them. Although these entanglement properties have been widely discussed for ground states, the way in which QSL form at finite temperature remains an open question. Here we introduce a method of characterizing both the depth and spatial structure of entanglement, and use this to explore how patterns of entanglement form as temperature is reduced in two widely studied models of QSL, the Kitaev honeycomb model, and the spin-1/2 Heisenberg antiferromagnet on a Kagome lattice. These results enable us to evaluate both the temperature at which spins within the high-temperature paramagnet first become entangled, and the temperature at which the system first develops the structured, multipartite entanglement characteristic of its QSL ground state.

Emergence of charge-2$e$ bosonic carriers as tightly bound electrons presents perhaps the simplest route to understand the non-Fermi liquid behaviors widely observed in functional materials. However, such scenario is typically discarded when unbound electrons are observed near the chemical potential. Here, using attractive Hubbard model as a representative example, we demonstrate the emergence of such ''bosonic'' bound pairs in the presence of residual low-energy electrons through determinant quantum Monte Carlo computation of their propagators. Furthermore, above the superfluid temperature the electronic spectral function is found to display a typical non-Fermi liquid behavior, namely a pseudogap near the chemical potential. This study provides the microscopic foundation for scenarios of boson-fermion mixed liquid as effective descriptions for some of the strongly correlated functional materials.

We investigate the periodic Anderson model composed of an itinerant $c$-band and a strongly localized $f$-band, featuring on-site electron-electron interactions in the $f$-orbitals. The two bands interact via a hybridization term with spin-orbit coupling, which enables spin-flip processes. In the non-interacting limit, these profoundly alter the electronic structure, leading to the emergence of flat bands, van Hove singularities, and, most notably, Dirac cones within a single Kondo-Dirac semimetal order. The strongly interacting regime is explored via the determinant quantum Monte Carlo method, in the absence of the sign problem, where we unveil a complete ground-state phase diagram revealing two distinct phases, the Kondo-Dirac semimetal phase and a novel antiferromagnetic semimetal phase. Their characterization by the spectral functions establishes an orbital-selective Mott transition in the antiferromagnetic semimetal phase, marked by the opening of a gap exclusively in the $f$-orbital while Dirac cones persist in the $c$-orbital. Conversely, in the Kondo-Dirac semimetal phase, both $c$- and $f$-orbitals sustain robust Dirac cones. We establish that spin-orbit coupling in the hybridization term gives rise to Dirac cones, which, combined with additional symmetry-breaking conditions, can generate novel topological states.

Light-matter interaction in frustrated Kagome metals enables access to hidden quantum states, yet the microscopic origin of symmetry breaking under ultrafast excitation remains elusive. Here, we uncover a microscopic mechanism for laser-induced symmetry breaking in RbV3Sb5 through first-principles real-time simulations. Selective excitation of a single-QM phonon mode dynamically breaks both rotational and time-reversal symmetries within the 2X2X1 charge density wave (CDW) superlattice. The resulting anisotropic lattice distortion lifts geometric frustration and stabilizes a nonequilibrium ferrimagnetic phase, accompanied by a sizable intrinsic anomalous Hall effect. Distinct from prior interpretations based on orbital antiferromagnetism or extrinsic perturbations, our findings reveal a spin-driven pathway for symmetry breaking under strong optical fields. These results provide a microscopic foundation for exploring how spin, lattice and charge degrees of freedom are intertwined in nonequilibrium correlated states.

The single-layered perovskite Ca$_2$RuO$_4$, characterized by a 4$d^4$ electron configuration, has been studied from first principles using density functional theory (DFT) using the generalized gradient approximation, with inclusion of strong on-site Coulomb interactions and spin-orbit coupling (GGA+SO+$U$), in the framework of the fully relativistic, spin-polarized Dirac linear muffin-tin orbital (LMTO) band-structure method. This approach enabled a comprehensive investigation of the electronic structure of Ca$_2$RuO$_4$ through the modeling of relevant spectra obtained from synchrotron-based techniques widely used to probe electronic properties, with a primary focus on resonant inelastic X-ray scattering (RIXS) at the Ru $L_3$ and O $K$ edges. The calculated spectra were thoroughly analyzed with available experimental data reported in the literature. The good agreement between our results and experimental observations for Ca$_2$RuO$_4$ enables a conclusive interpretation of key features in the spectra obtained from the aforementioned techniques. Consequently, this enables us to describe its electronic properties and to establish a solid theoretical approach suitable for routine modeling of spectra, particularly from RIXS, aimed at characterizing the electronic structure and properties of similar or more complex strongly correlated, technologically relevant materials.

The single-layered perovskite Ca$_2$RuO$_4$, characterized by a 4$d^4$ electron configuration, has been studied from first principles using density functional theory (DFT) using the generalized gradient approximation, with inclusion of strong on-site Coulomb interactions and spin-orbit coupling (GGA+SO+$U$), in the framework of the fully relativistic, spin-polarized Dirac linear muffin-tin orbital (LMTO) band-structure method. This approach enabled a comprehensive investigation of the electronic structure of Ca$_2$RuO$_4$ through the modeling of relevant spectra obtained from synchrotron-based techniques widely used to probe electronic properties, with a primary focus on resonant inelastic X-ray scattering (RIXS) at the Ru $L_3$ and O $K$ edges. The calculated spectra were thoroughly analyzed with available experimental data reported in the literature. The good agreement between our results and experimental observations for Ca$_2$RuO$_4$ enables a conclusive interpretation of key features in the spectra obtained from the aforementioned techniques. Consequently, this enables us to describe its electronic properties and to establish a solid theoretical approach suitable for routine modeling of spectra, particularly from RIXS, aimed at characterizing the electronic structure and properties of similar or more complex strongly correlated, technologically relevant materials.

Recently plasmon excitations in bilayer lattice systems were studied extensively in the weak-coupling regime. Unlike single-layer systems, these bilayers exhibit two distinct modes, $ω_{\pm}$, which show characteristic dependences upon the momentum and hopping integrals along the $z$ direction. To apply them to cuprates, strong correlation effects should be considered, but a comprehensive analysis has not yet been investigated. In this work, we present a strong-coupling theory to analyze the charge dynamics of a bilayer system, utilizing the $t$-$J$-$V$ model, which includes the long-range Coulomb interaction, $V$, on a lattice. Although our theoretical framework is fundamentally different from the weak-coupling approach, we find that resulting plasmon excitations are similar to those of a weak-coupling theory. A key distinction is that our strong-coupling framework reveals a noticeable suppression of particle-hole excitations, which allows the plasmon modes to remain well-defined over a wider region of momentum. We suggest that the experimentally reported plasmon excitations in Y-based cuprates can be described by the $ω_{-}$ mode, although we call for more systematic experiments to verify this.

We construct parametrized isometric tensor network states -- referred to as skeletons -- that allow us to explore phases of abelian topological order and can be efficiently implemented on quantum processors. We obtain stable finite correlation length deformations of string-net fixed points, which are constructed both by conserving virtual symmetries of the tensor and by imposing local isometry constraints. They connect distinct topological phases via a shared critical point, thereby providing analytically tractable examples of phase transitions beyond anyon condensation. By mapping such classes of 2D tensor networks to 1D stochastic automata with local update rules, we show that expectation values of generalized Pauli strings of arbitrary weight can be efficiently computed using classical methods. Therefore these skeletons not only serve as an organizing principle for abelian topological order but also provide a non-trivial testbed for quantum processors.

Topological materials--including insulators (TIs) and semimetals (TSMs)--hold immense promise for quantum technologies, yet their discovery remains constrained by the high computational cost of first-principles calculations and the slow, resource-intensive nature of experimental synthesis. Here, we introduce TXL Fusion, a hybrid machine learning framework that integrates chemical heuristics, engineered physical descriptors, and large language model (LLM) embeddings to accelerate the discovery of topological materials. By incorporating features such as space group symmetry, valence electron configurations, and composition-derived metrics, TXL Fusion classifies materials across trivial, TSM, and TI categories with improved accuracy and generalization compared to conventional approaches. The framework successfully identified new candidates, with representative cases further validated through density functional theory (DFT), confirming its predictive robustness. By uniting data-driven learning with chemical intuition, TXL Fusion enables rapid and interpretable exploration of complex materials spaces, establishing a scalable paradigm for the intelligent discovery of next-generation topological and quantum materials.

A two-dimensional resonant inelastic x-ray scattering (2D-RIXS) microscopy system has been developed at the beamline BL02U of NanoTerasu. The instrument combines a Wolter type-I mirror for spatial imaging with a varied-line-spacing grating spectrometer, simultaneously achieving micrometer-scale spatial resolution and ultrahigh energy resolution in the soft x-ray regime. Test chart measurements confirm a vertical spatial resolution of 1.0 um near the field-of-view center, and the horizontal resolution determined by the incident beam footprint is 0.8 um. RIXS imaging capabilities have been demonstrated by the measurements of a patterned NanoTerasu logo and exfoliated NiPS${}_3$ nanoflakes, highlighting its efficiency in locating specific microscale regions within inhomogeneous samples. These results establish 2D-RIXS microscopy as a position-sensitive probe of elementary excitations in quantum materials and functional devices.

We present a novel framework for controlling Higgs mode and vortex dynamics in superconductors using structured light. We propose a phenomenon analog of the Kapitza-Dirac effect in superconductors, where Higgs waves scatter off light-induced vortex lattices, generating interference patterns akin to matter wave diffraction. We also find that the vortices enable the linear coupling of Higgs mode to the electromagnetic field. This interplay between light-engineered Higgs excitations and emergent vortex textures opens a pathway to probe nonequilibrium superconductivity with unprecedented spatial and temporal resolution. Our results bridge quantum optics and condensed matter physics, offering new examples of quantum printing where one uses structured light to manipulate the collective modes in correlated quantum fluids.

We report two types of singularities that arise from fluctuations during the formation of charge- or spin-density waves. The first is the exceptional point (EP), corresponding to a higher-order pole of the retarded Green's function. Such EPs lead to algebraic corrections in the decay of quasiparticle occupations and are observable through time-resolved angle-resolved photoemission spectroscopy (Tr-ARPES). The second is a spectral cusp, defined by the coalescence of three extrema in the real-frequency spectral function $A(\mathbf{k}, ω)$. This cusp enforces the formation of Fermi arcs and induces a "threading" structure in the nearby band structure, both of which are directly observable in ARPES.

The study of real-time dynamics of fermions remains one of the last frontiers beyond the reach of classical simulations and is key to our understanding of quantum behavior in chemistry and materials, with implications for quantum technology. Here we introduce a Heisenberg-picture algorithm that propagates observables expressed in a Majorana-string basis using a truncation scheme that preserves Trotter accuracy and aims at maintaining computational efficiency. The framework is exact for quadratic Hamiltonians--remaining restricted to a fixed low-weight sector determined by the physical observable--admits variational initial states, and can be extended to interacting regimes via systematically controlled truncations. We benchmark our approach on one- and two-dimensional Fermi-Hubbard quenches, comparing against tensor network methods (MPS and fPEPS) and recent experimental data. The method achieves high accuracy on timescales comparable to state-of-the-art variational techniques and experiments, demonstrating that controlled Majorana-string truncation is a practical tool for simulating two-dimensional fermionic dynamics.

We study the interplay of magnetic order and superconductivity in the square-lattice Hubbard model under periodic driving with circularly polarized light. Formulating diagrammatic techniques based on the random-phase approximation in terms of Floquet Green's functions, allows us to analyze fluctuation-driven unconventional pairing for weak-to-moderate interactions. The interplay of repulsive interactions and photo-assisted hopping of electrons gives rise to a rich magnetic phase diagram featuring an antiferromagnetic-to-ferromagnetic crossover prior to a Floquet Lifshitz transition. Close to the antiferromagnetic transition, topological $d+id$-wave superconductivity prevails in the phase diagram for a wide range of drive parameters. At intermediate-to-high frequency driving near the Floquet Lifshitz transition, superconducting orders are tuned from spin-singlet $d$-wave to spin-triplet $p$-wave character, providing an effective protocol for Floquet engineering topological superconductivity.

We study the spin-1/2 Kitaev honeycomb gapless spin liquid in the presence of Stone-Wales-type local lattice defects with odd-sided plaquettes. While the clean Kitaev model has no finite-temperature phase transitions, we find that introducing a finite defect density $n_d\approx 10^{-4}$--$10^{-2}$ produces a true phase transition with a sizeable $T_c \approx 2 n_d$ in units of the Kitaev exchange. The resulting non-Abelian chiral quantum spin liquid exhibits scalar spin chirality and electron orbital magnetization which peak near lattice defects. This disorder-driven instability relies on an emergent long range ferromagnetic interaction $r^{-γ}$ ($γ\approx 2.7$) between defect chiralities, mediated by the nearly-gapless fermions, with implications for topology generation in Dirac cones with fluctuating mass terms.

Quasiparticle hybridization remains a major challenge to realizing and controlling exotic states of matter in existing quantum simulation platforms. We report the absence of hybridization for compact localized states (CLS) emerging in the chiral spin liquid described by the Yao-Kivelson model. The CLS form due to destructive quantum interference at fine-tuned coupling constants and populate perfectly flat quasiparticle bands on an effective kagome lattice. Using a formalism for general Majorana-hopping Hamiltonians, we derive exact expressions for CLS for various flux configurations and both for the topological and trivial phases of the model. In addition to finite-energy matter fermions with characteristic spin-spin correlations, we construct compact localized Majorana zero modes attached to $\pi$-flux excitations, which enable non-Abelian braiding of Ising anyons with minimal separation. Our results inform the quantum simulation of topologically ordered states of matter and open avenues for exploring flat-band physics in quantum spin liquids.

We revisit gradient-based optimization for infinite projected entangled pair states (iPEPS), a tensor network ansatz for simulating many-body quantum systems. This approach is hindered by two major challenges: the high computational cost of evaluating energies and gradients, and an ill-conditioned optimization landscape that slows convergence. To reduce the number of optimization steps, we introduce an efficient preconditioner derived from the leading term of the metric tensor. We benchmark our method against standard optimization techniques on the Heisenberg and Kitaev models, demonstrating substantial improvements in overall computational efficiency. Our approach is broadly applicable across various contraction schemes, unit cell sizes, and Hamiltonians, highlighting the potential of preconditioned optimization to advance tensor network algorithms for strongly correlated systems.

We study the extended Hubbard model (EHM) with both onsite Hubbard interaction and the intersite density-density interaction between nearest-neighbors using the standard Hartree mean-field approximation (MFA) on the Bethe lattice. We found that, at the ground state, the system can be in a charge-ordered insulating (COI), a charge-order metallic (COM) or a non-charge-ordered (NO) state. Moreover, the finite-temperature phase diagrams are presented. Several observables like a charge-order parameter, a spectral function, and particularly at finite temperatures, a charge carrier concentration (to visualize the degree of metallicity) are analyzed. The results show that increasing onsite repulsion suppresses charge order and change the properties of the system from insulating to metallic. Worth noting, that a number of phenomena can be found within the MFA, where their analysis is much simpler than in more advanced approaches. The method used for the EHM on the Bethe lattice also allows for a series of analytical derivations and simplification to see general geometry-independent features and analytical results, avoiding the numerical inaccuracies and other issues that appear with a purely numerical solution.

Understanding charge transfer at oxide interfaces is crucial for designing materials with emergent electronic and magnetic properties, especially in systems where strong electron correlations and spin-orbit coupling coexist. SrIrO$_3$/SrCoO$_3$ (SIO/SCO) superlattices offer a unique platform to explore these effects due to their contrasting electronic structures and magnetic behaviors. Building on past theory based on continuity of O 2p band alignment, we employ density functional theory (DFT) to model electron transfer from Ir to Co across the SIO/SCO interface. To characterize these effects, we synthesized epitaxial SIO/SCO superlattices via molecular beam epitaxy. Structural and transport measurements confirmed high crystallinity, metallic behavior, and suppression of Kondo scattering that has been reported in uniform SIO films. Further characterization via X-ray absorption spectroscopy (XAS) revealed orbital anisotropy and valence changes consistent with interfacial charge transfer. Co K- and L$_{2,3}$-edge and Ir L$_2$-edge spectra verified electron donation from Ir to Co, stabilizing the perovskite SCO phase and tuning the electronic structure of SIO via hole-doping. O K-edge XAS showed band alignment shifts in the SIO layer consistent with DFT predictions. Our work here provides a pathway for engineering oxide heterostructures with tailored magnetic and electronic properties.

Altermagnets (AMs) are a recently identified class of unconventional collinear compensated antiferromagnets that exhibit momentum-dependent spin splitting despite having zero net magnetization. This unconventional magnetic order gives rise to a range of phenomena, including the anomalous Hall effect, chiral magnons, and nonlinear photocurrents. Here, using spin space group (SSG) symmetry analysis and first-principles calculations, we demonstrate an efficient strategy to control altermagnetism in two-dimensional multiferroics through ferroelectric polarization and interlayer sliding. For material realization, we find that monolayer and bilayer FeCuP2S6 exhibit finite spin splitting when ferroelectric sublattices are connected by nonsymmorphic screw-axis operations rather than pure translation or inversion symmetry. Interlayer sliding further enables reversible switching or suppression of spin splitting through modifications of the SSG. Our calculations further reveal that the anomalous Hall response serves as a direct probe of these spin-split states. These findings establish two-dimensional van der Waals multiferroics as promising platforms for realizing electrically controllable altermagnetism and advancing next-generation spintronic and magnetoelectric technologies.

Symmetry breaking (SB) has proven to be a powerful approach for describing quantum materials: strong correlation, mass renormalization, and complex phase transitions are among the phenomena that SB can capture, even when coupled to a mean-field-like theory. Traditionally, corrective schemes were required to account for these effects; however, SB has emerged as an alternative that can also successfully describe the intricate physics of quantum materials. Here, we explore spin SB on EuB6 and SmB6 and how its relation to the exchange field can determine onsite properties, depending on the type of symmetry breaking. Using spin-polarized Density Functional Theory (DFT) calculations with the r2SCAN functional, we systematically compare four magnetic configurations, one totally symmetric - non-magnetic (NM) configuration - and three with different types of symmetry breaking: ferromagnetic (FM), antiferromagnetic (AFM) and a paramagnetic (PM) configuration - modeled through a Special Quasirandom Structure (SQS) method - to capture local symmetry-breaking effects. Our results show that the PM configuration produces distinct magnetic environments for the rare-earth atoms, leading to different exchange fields. These, in turn, induce symmetry breaking in the electronic and magnetic properties of Eu and Sm. Those results provide an alternative explanation for the experimental results on both materials, EuB6 and SmB6, where X-ray Absorption Spectroscopy (XAS) and X-ray Absorption Near Edge Structure (XANES) measurements suggest the presence of multiple atomic environments, previously attributed to a mixed-valence configuration.

The accurate simulation of strongly correlated electron systems remains a central challenge in condensed matter physics, motivating the development of various non-perturbative many-body methods. Such methods are typically benchmarked against the numerical exact determinant quantum Monte Carlo (DQMC) in the Hubbard model; however, DQMC is limited by the fermionic sign problem and the uncertainties of numerical analytic continuation. To address these issues, we use the exactly solvable Hatsugai-Kohmoto (HK) model as a benchmarking platform to evaluate three many-body approximations: $GW$, $HGW$, and $SGW$. We compare the Green's functions, spectral functions, and response functions obtained from these approximations with the exact solutions. Our analysis shows that the $GW$ approximation, often considered insufficient for describing strong correlation, exhibits a previously unreported solution branch that accurately reproduces Mott physics in the HK model. In addition, using a covariant formalism, we find that $HGW$ provides an accurate description of charge response, while $SGW$ performs well for spin correlations. Overall, our work demonstrates that the HK model can effectively benchmark many-body approximations and helps refine the understanding of $GW$ methods in strongly correlated regimes.

The accurate treatment of electron correlation in extended molecular systems remains computationally challenging using classical electronic structure methods. Hybrid quantum-classical algorithms offer a potential route to overcome these limitations; however, their practical deployment on existing quantum computers requires strategies that both reduce problem size and mitigate hardware noise. In this work, we combine Density Matrix Embedding Theory (DMET) with Sample-based Quantum Diagonalization (SQD) to compute ground-state energies of a set of natural ligand-like molecules in the minimal Slater Type Orbital (STO-3G) basis set. DMET provides a systematic fragmentation of a molecule into embedded impurity subproblems, while SQD enables construction and classical diagonalization of reduced configuration spaces through quantum sampling enhanced by iterative configuration recovery. The resulting embedded Hamiltonians are solved on IBM's Eagle R3 superconducting quantum hardware (IBM Sherbrooke). The DMET-SQD energies obtained for all systems considered exhibit strong agreement with DMET-FCI benchmark values within chemical accuracy (1 kcal/mol). These results demonstrate that sample-based quantum methods, when integrated with a robust embedding framework, can reliably extend quantum computation towards simulation of chemically relevant molecular systems, showcasing potential applications in the field of drug discovery.

$n$-particle reduced density matrices ($n$-RDMs) play a central role in understanding correlated phases of matter. Yet the calculation of $n$-RDMs is often computationally inefficient for strongly-correlated states, particularly when the system sizes are large. In this work, we propose to use neural network (NN) architectures to accelerate the calculation of, and even predict, the $n$-RDMs for large-size systems. The underlying intuition is that $n$-RDMs are often smooth functions over the Brillouin zone (BZ) (certainly true for gapped states) and are thus interpolable, allowing NNs trained on small-size $n$-RDMs to predict large-size ones. Building on this intuition, we devise two NNs: (i) a self-attention NN that maps random RDMs to physical ones, and (ii) a Sinusoidal Representation Network (SIREN) that directly maps momentum-space coordinates to RDM values. We test the NNs in three 2D models: the pair-pair correlation functions of the Richardson model of superconductivity, the translationally-invariant 1-RDM in a four-band model with short-range repulsion, and the translation-breaking 1-RDM in the half-filled Hubbard model. We find that a SIREN trained on a $6\times 6$ momentum mesh can predict the $18\times 18$ pair-pair correlation function with a relative accuracy of $0.839$. The NNs trained on $6\times 6 \sim 8\times 8$ meshes can provide high-quality initial guesses for $50\times 50$ translation-invariant Hartree-Fock (HF) and $30\times 30$ fully translation-breaking-allowed HF, reducing the number of iterations required for convergence by up to $91.63\%$ and $92.78\%$, respectively, compared to random initializations. Our results illustrate the potential of using NN-based methods for interpolable $n$-RDMs, which might open a new avenue for future research on strongly correlated phases.

Characterizing topological phases for strongly interacting fermions in the mixed-state regime remains a major challenge. Here we introduce a general and numerically efficient framework to diagnose mixed-state topological phases in strongly interacting systems via the disorder parameter (DP) of the U(1) charge operator. Specifically, from the finite-size scaling of the second derivative of the DP generating function, we introduce the topological scaling indicator, which exhibits a characteristic linear scaling with the system's linear dimension for topological phases, a signature that vanishes upon transition into a topologically trivial phase. Crucially, we develop an efficient determinant Quantum Monte Carlo algorithm that facilitates the evaluation of this indicator in interacting systems. We apply our approach to two paradigmatic models: for the Kane-Mele-Hubbard model, we successfully map the interaction-driven transition from a quantum spin Hall insulator to a trivial Mott insulator. Furthermore, our method circumvents the limitations imposed by the severe sign problem in the Haldane-Hubbard model, enabling robust identification of the quantum anomalous Hall phase at accessible temperatures. This work provides a powerful and accessible tool for the numerical exploration of topological phenomena in interacting mixed states, opening a pathway to study systems previously inaccessible due to computational obstacles.

We present a data-driven analysis of correlation hierarchies in the quarter-filled simple and extended Hubbard models by applying principal component analysis (PCA) to exact-diagonalization (ED) data on 3x4 and 4x4 cylindrical clusters. While the non-interacting limit (U=0) provides a finite-size reference, increasing on-site repulsion U induces localization and reorganizes the low-energy spectrum. For the extended model, we examine moderate (U=4) and strong (U=10) coupling regimes, where conventional structure factors reveal familiar crossovers among charge, spin and local-pairing correlations. PCA of the corresponding correlation matrices captures these crossovers directly from the data, without assuming predefined order parameters by identifying charge-dominated, spin-dominated and pairing-dominated regimes through variance condensation into leading components. This establishes PCA as a transparent, model-agnostic framework for uncovering the hierarchy and competition of correlation channels in finite Hubbard lattice clusters, providing a bridge between exact diagonalization and modern machine-learning diagnostics in strongly correlated systems.

We use the Gaussian Phase-Space Representation to solve the real-time dynamic of interacting fermions in 1D, 2D, and 3D systems. The method is exact up to a spiking point, which represents a limit on the practical simulation time. The spiking can be delayed, and the practical simulation time extended, by adjusting the gauges of the representation, resulting in different equivalent stochastic differential equations. Here, we work on the so-called diffusion gauge and propose an algorithm to find efficiently new implementations of the noise terms. Compared with our initial results [F. Rousse \textit{et al.} 2024, J. Phys. A: Math. Theor. \textbf{57}, 015303], the new method achieves a significantly longer practical simulation time and can be applied to significantly larger systems.

The interplay between unconventional density waves and exotic superconductivity has attracted growing interest. Kagome superconductors $A\rm{V}_3\rm{Sb}_5$ ($A = \rm{K}, \rm{Rb}, \rm{Cs}$) offer a platform for studying quantum phase transitions and the resulting symmetry breaking. Among these quantum phases, the $4a_0$ stripe charge-density-wave (CDW) has been widely observed for $A=\rm{Rb}$ and $\rm{Cs}$ by scanning tunneling microscopy (STM) and nuclear magnetic resonance (NMR) measurements. However, the microscopic origin of the $4a_0$ stripe CDW remains elusive, and no theoretical studies addressing this phenomenon have been reported so far. In this paper, we propose a microscopic mechanism for the emergence of the $4a_0$ stripe CDW. We analyze the CDW instability in the 12-site kagome lattice Hubbard model with the $2\times2$ bond order driven by the paramagnon-interference mechanism by focusing on the short-range magnetic fluctuations due to the geometrical frustration of kagome lattice. We reveal that the nesting vector of the reconstructed Fermi surface, formed by the $2\times 2$ bond order, gives rise to a $4a_0$-period CDW. Remarkably, the obtained stripe CDW is composed of both the off-site hopping integral modulations and on-site potentials. The real-space structure of the stripe CDW obtained here is in good qualitative agreement with the experimentally observed stripe pattern.

We consider the electronic and magnetic properties of the Hubbard model on a square lattice with the Fermi level near van Hove singularity and the ratio of the next-nearest-neighbor and nearest-neighbor hoppings $t'/t=-0.45$, which favours the ferromagnetic instability. We find, that a self-consistent consideration of the ferromagnetic fluctuations within the D-TRILEX approach results in the splitting of the electronic spectral function at low temperatures. This splitting exhibits only a weak momentum dependence, and only one of the split bands crosses the Fermi level. As a result, the Fermi surface itself remains unsplit, but its area increases, reflecting the presence of non-Fermi-liquid electronic excitations. We show that both the self-consistent account of the non-local contributions to the electronic self-energy and the proper treatment of electron interaction vertices in D-TRILEX are important to obtain this behaviour.

An altermagnet is a recently identified class of magnets that exhibit a zero net magnetic moment but break symmetry under the combined operations of parity and time reversal. It typically consists of two magnetic sites of opposite spins related by rotation within the unit cell. Here, we use circular dichroism (CD) in resonant inelastic X-ray scattering (RIXS) to identify a new form of altermagnetism, namely bi-altermagnetism, in the correlated insulator Fe2Mo3O8, which comprises two altermagnetic sublattices: one with alternating quasi-octahedral Fe environments and the other with alternating tetrahedral Fe environments. We experimentally revealed the emergence of CD in an achiral, zero-magnetization system, thereby probing mirror-symmetry breaking associated with altermagnetic order. Notably, the CD appeared at sublattice-specific excitations of the octahedral and tetrahedral sites, indicating symmetry breaking in both altermagnetic sublattices. Calculations based on a model with the bi-altermagnetic order along the c axis successfully reproduce the observed CD. Our findings provide compelling evidence for bi-altermagnetism in Fe2Mo3O8, and showcase the use of RIXS-CD as a probe of magnetic sublattices in systems with zero net magnetization.

Bond-dependent magnetic interactions, particularly those described by the Kitaev model, have emerged as a key pathway toward realizing unconventional magnetic states such as quantum spin liquids and topologically nontrivial excitations, including skyrmions. These interactions frustrate conventional magnetic order and give rise to rich collective behavior that continues to challenge both theory and experiment. While Kitaev physics has been extensively explored in the context of honeycomb magnets, direct evidence for its role in real materials remains scarce. Magnetic van der Waals (vdW) materials have emerged as a versatile platform for exploring low-dimensional electrical, magnetic, and correlated electronic phenomena, and provide a fertile ground for potential applications ranging from spintronics to multiferroic devices and quantum information technologies. Here, we demonstrate, through magneto-transmission, Faraday angle rotation, and magnetic circular dichroism measurements, that the magnetic excitation spectrum of NiI$_2$, a van der Waals multiferroic material, is more accurately captured by a Kitaev-based spin model than by the previously invoked helical spin framework.

Fractionalized excitations such as spinons and anyons have emerged as a central theme in condensed matter physics with broad implications for superconductivity, quantum statistics, and quantum computation. The nearly ideal one-dimensional $S=1/2$ system $β$-VOSO$_4$ without long-range order down to 85 mK provides a promising platform to experimentally explore such fractionalized excitations. Here, we employ Raman spectroscopy to probe magnetic excitations and the evolution of spin correlations in $β$-VOSO$_4$. Spinon signatures are found along the chain direction, evidenced by a broad, gapless scattering continuum at low temperatures. The temperature dependence of the spinon spectral weight aligns considerably with numerical density matrix renormalization group calculations. By comparing the experimental spinon spectral weight with calculated results and evaluating the associated quantum Fisher information (QFI) therefrom, we observe a steep increase in QFI upon cooling, indicating rapidly growing correlation lengths. Our study showcases QFI as a probe of spin correlations in quantum magnets.

We establish a rigorous density functional theory (DFT) framework for core-level X-ray absorption spectroscopy (XAS) by formulating a constrained search for core-excited states based on the Gunnarsson-Lundqvist theorem. Within this framework, the explicit-core Delta SCF scheme enables shift-free absolute edge alignment and a consistent treatment of L/M edges with spin-orbit-resolved projectors. In addition, by exploiting dipole selection rules, we recast the evaluation of the dipole matrix elements, which otherwise requires many independent Slater determinant calculations, into a compact single determinant form. This reduces the computational scaling from $\mathcal{O}(N^4)$ to $\mathcal{O}(N^3)$, where $N$ is the number of electrons, without introducing additional approximations. Across representative C, B, O, and Li K-edge benchmarks in molecules and solids, the method reproduces line shapes, polarization anisotropies, and absolute onsets without empirical shifts, providing a robust and scalable route to quantitatively reliable XAS simulations within DFT.

Light-matter interaction in frustrated Kagome metals enables access to hidden quantum states, yet the microscopic origin of symmetry breaking under ultrafast excitation remains elusive. Here, we uncover a microscopic mechanism for laser-induced symmetry breaking in RbV3Sb5 through first-principles real-time simulations. Selective excitation of a single-QM phonon mode dynamically breaks both rotational and time-reversal symmetries within the 2X2X1 charge density wave (CDW) superlattice. The resulting anisotropic lattice distortion lifts geometric frustration and stabilizes a nonequilibrium ferrimagnetic phase, accompanied by a sizable intrinsic anomalous Hall effect. Distinct from prior interpretations based on orbital antiferromagnetism or extrinsic perturbations, our findings reveal a spin-driven pathway for symmetry breaking under strong optical fields. These results provide a microscopic foundation for exploring how spin, lattice and charge degrees of freedom are intertwined in nonequilibrium correlated states.

The functional renormalization group (fRG) is acknowledged as a powerful tool in quantum many-body physics and beyond. On the technical side, conventional implementations of the fRG rely on regulators for bare propagators only. Starting from Schwinger--Dyson and Bethe--Salpeter equations, we develop here an fRG formulation where both bare propagators and bare interactions can be dressed with regulators. The approach thus obtained is a generalization of the multiloop fRG recently introduced for many-fermion systems. Using the single-boson exchange decomposition, we show that the underlying flow equations are simply interpreted as adding a regulator to the bosonic propagator and that such an extension scarcely changes the original structure of the flow equations. Overall, we provide a framework for implementing approaches that cannot be realized with conventional fRG methods, such as temperature flows for models with retarded interactions. For concrete applications, we analyze the loop convergence of our scheme against conventional cutoff schemes for the Hubbard atom and the Anderson impurity model. Finally, we devise a new temperature-flow scheme that implements a cutoff in both the propagator and the bare interaction, and demonstrate its validity on a model of an Anderson impurity coupled to a phonon.

4Hb-TaS_2, a van der Waals superlattice comprising alternate stacked Ising superconducting 1H-TaS_2 and cluster Mott insulating 1T-TaS_2, exhibits emergent properties beyond those of its constituent layers. Notable phenomena include time-reversal-symmetry-breaking superconductivity and spontaneous vortex phases, which are driven by nontrivial interlayer interactions that remain debated. Using area-selective angle-resolved photoemission spectroscopy, we provide direct spectroscopic evidence of such interaction by systematically probing the electronic structures of 1T- and 1H-terminted surfaces of 4Hb-TaS_2. The metallic states of subsurface 1H-layers are folded to the Brillouin zone center by the sqrt(13) by sqrt(13) modulation of the surface 1T-layer, forming chiral "windmill" Fermi surfaces via Umklapp scattering. These conducting states further hybridize with the incipient flat band of the surface 1T-layer, producing a Kondo-like peak at the Fermi level. Interlayer charge transfer induces distinct 3 by 3 and 2 by 2 charge orders on the surface and subsurface 1H-layers, respectively, which result in characteristic segmented Fermi surfaces and dichotomously shift the van Hove singularities. These findings reconcile the competing Kondo and Mott-Hubbard models in this material and emphasize the interplay of flat bands, van hove singularities, charge orders, and unconventional superconductivity in correlated superlattices.

The doped quantum spin liquid on the kagome lattice provides a fascinating platform to explore exotic quantum states, such as the reported holon Wigner crystal at low doping. By extending the doping range to $δ= 0.027$ - $0.36$, we study the kagome-lattice $t$-$J$ model using the state-of-the-art density matrix renormalization group calculation. On the $L_y=3$ cylinder ($L_y$ is the number of unit cells along the circumference direction), we establish a quantum phase diagram with increasing doping level. In addition to the charge density wave (CDW) states at lower doping, we find an emergent Fermi-liquid-like phase by melting the holon Wigner crystal at $δ\approx 0.15$, which is characterized by suppression of charge density oscillation and power-law decay of various correlation functions. On the wider $L_y = 4$ cylinder, the bond-dimension extrapolated correlation functions also support such a Fermi-liquid-like state, suggesting its stability with increasing system size. In a narrow doping range near $δ= 1/3$ on the $L_y = 3$ cylinder, we find a state with an exponential decay of single-particle correlation but the other correlation functions preserving the features in the Fermi-liquid-like phase, which may be a precursor of a superconducting state. Nevertheless, this peculiar state near $δ= 1/3$ disappears on the $L_y = 4$ cylinder, implying a possible lattice size dependence. Our results reveal a quantum melting from a holon Wigner crystal to a Fermi-liquid-like state with increasing hole density, and suggest a doping regime to explore superconductivity for future study.

The precise role of e-ph coupling in graphene and related materials on a honeycomb lattice is not yet fully understood, despite extensive research on these systems. Here, we perform sign-problem-free determinant quantum Monte Carlo (DQMC) simulations of the optical Su-Schrieffer-Heeger (oSSH)-Hubbard model on the honeycomb lattice, focusing on the parameters relevant to graphene. Performing finite-size scaling analyzes, we obtain the model's ground state phase diagram, which includes the semi-metal (SM), Kekulé Valence Bond Solid (KVBS), and anti-ferromagnetic (AFM) phases, as well as indications of a small KVBS/AFM coexistence region. We find that a weak to moderate Hubbard repulsion, tuned toward the SM-AFM critical value in the pure honeycomb Hubbard model, enhances KVBS correlations and can even stabilize the KVBS phase. Estimating the effective parameters for graphene places it in the SM region of the phase diagram, but near the SM-KVBS phase boundary. Notably, we predict that increasing either the on-site Hubbard repulsion or the e-ph coupling strength drives graphene toward the KVBS phase rather than the AFM phase, highlighting a synergistic effect that can be exploited to further control the remarkable properties of graphene and related materials.

Altermagnets are a class of collinear magnets that exhibit non-relativistic spin splitting (NRSS) of electronic bands in the absence of net magnetization. Their potential to generate large spin polarization without spin-orbit coupling has created strong interest in probes that access the underlying order parameter directly. In this Perspective, we show that linear magneto-birefringence (LMB) provides a natural and broadly applicable route to detecting altermagnetic order. Building on the correspondence between the momentum-space structure of NRSS and the ferroic ordering of magnetic multipoles in real space, we demonstrate how $d$-wave and $g$-wave NRSS textures yield distinct LMB responses. We present a symmetry-based framework that identifies the optical geometries and field configurations required to isolate specific multipole components, enabling domain imaging and providing benchmarks for theoretical models of LMB.

We investigate the interplay between altermagnetic band structures and electronic correlations by focusing on the $d_{x^2-y^2}$ altermagnetic generalization of the Hatsugai-Kohmoto model. We find that with increasing interaction, a many-body Lifshitz transition takes place when doubly occupied regions disappear from the Fermi surface and each momentum state becomes fully spin polarized. The spin susceptibility is directly evaluated from the Kubo formula in terms of many-body occupation probabilities. We find that the dynamical susceptibility, which possesses only transverse non-zero components for small wavevectors, develops a gap proportional to the interaction strength, and displays a sharp peak at a frequency increasing with the interaction. %with increasing frequency. Above the Lifshitz transition, this peak moves to the lower gap edge and becomes log-divergent. The signal intensity increases with the interaction up until the Lifshitz transition and saturates afterwards. The static susceptibility remains unaffected by the correlations and altermagnetism reduces the static transverse response.

Quantum spin liquids can arise from Kitaev magnetic interactions, and exhibit fractionalized excitations with the potential for a topological form of quantum computation. This review surveys recent experimental and theoretical progress on the pursuit of phenomena related to Kitaev magnetism in layered and exfoliatable materials, which offer numerous opportunities to apply powerful techniques from the field of atomically thin materials. We primarily focus on the antiferromagnetic Mott insulator $α$-RuCl3, which exhibits Kitaev couplings and is readily exfoliated to single- or few-layer sheets, and thus serves as a test bed for developing probes of Kitaev phenomena in atomically thin materials and devices. We introduce the Kitaev model and how it is realized in $α$-RuCl3 and other material candidates; and cover $α$-RuCl3 synthesis and fabrication into van der Waals heterostructure devices. A key discovery is a work-function-mediated charge transfer that heavily dopes both the $α$-RuCl3 and proximate materials, and can enhance Kitaev interactions by up to 50%. We further discuss a wide range of recent results in electronic transport and optical and tunneling spectroscopies of $α$-RuCl3 devices. The experimental techniques and theoretical insights developed for $α$-RuCl3 establish a framework for discovering and engineering superior two-dimensional Kitaev materials that may ultimately realize elusive quantum spin liquid phases.

Altermagnetism, characterized by zero net magnetization and symmetry-protected spin-split band structures, has recently emerged as a promising platform for spintronics. In altermagnetic tunnel junctions (AMTJs), the suppression of tunneling in the antiparallel configuration relies on the mismatch between spin-polarized conduction channels in momentum space. However, ideal nonoverlapping spin-polarized Fermi surfaces are rarely found in bulk altermagnets. Motivated by the critical influence of Fermi surface geometry on tunneling magnetoresistance (TMR), we investigate three experimentally synthesized altermagnets -- bulk $\mathrm{V_2Te_2O}$, $\mathrm{RbV_2Te_2O}$, and $\mathrm{KV_2Se_2O}$ -- to elucidate how flatband-driven Fermi surfaces minimize spin-channel overlap and boost AMTJ performance. Notably, $\mathrm{RbV_2Te_2O}$ and $\mathrm{KV_2Se_2O}$ host flat altermagnetic Fermi sheets, which confine spin degeneracy to minimal arc-like or nodal-like regions. Such Fermi surface geometry drastically reduces spin overlap, resulting in an unprecedented intrinsic TMR well over $10^3\%$ in the $\mathrm{KV_2Se_2O}$-based AMTJ. Incorporating an insulating barrier further enhances the TMR to $\sim10^5\%$, surpassing most conventional MTJs. These results not only establish $\mathrm{KV_2Se_2O}$ as a compelling candidate AMTJ material, but also highlight the critical role of flatband Fermi surface geometry in achieving high-performance altermagnetic-spintronic device technology.

Materials which show a strong time-reversal symmetry-breaking response leading to spin-polarization phenomena, in conjunction with antiparallel magnetic alignments producing zero net magnetization, have recently been identified, classified, and been given the name 'altermagnets'. However, measuring and diagnosing possible candidates as altermagnetics still remains a challenge. From the uncertainty of the material being an altermagnet, additional experimental probes are essential to resolve this. Here, we propose using spin-dependent and magnetic momentum density probes such as spin-polarised positron annihilation and revisiting magnetic Compton scattering. By looking at the previously claimed altermagnetic candidates RuO2, CrSb and MnTe, we present theoretical altermagnetic calculations of the experimental quantities measured by these probes. We show that these quantities should produce a measurable signal and unequivocally confirm the altermagnetic state. We also highlight the additional benefits from these probes such as extracting spin-resolved Fermi surfaces which are key for further understanding the nature of the altermagnetic state.

Even in the quantum limit, non-equilibrium quasiparticle (QP) populations induce QP poisoning that irreversibly relaxes the quantum state and significantly degrades the coherence of transmon qubits. A particularly detrimental yet previously unexplored mechanism arises from QP multiplication facilitated by vortex trapping in superconducting quantum circuits, where a high-energy QP relaxes by breaking additional Cooper pairs and amplifying the QP population due to the locally reduced excitation gap and enhanced quantum confinement within the vortex core. Here we directly resolve this elusive QP multiplication process by revealing vortex-controlled QP self-generation in a highly nonequilibrium regime preceding the phonon bottleneck of QP relaxation. At sufficiently low fluence, femtosecond-resolved magneto-reflection spectroscopy directly reveals a continuously increasing QP population that is strongly dependent on magnetic-field-tuned vortex density and absent at higher excitation fluences. Quantitative analysis of the emergent QP pre-bottleneck dynamics further reveals that, although the phonon population saturates within $\simeq$10~ps, both free and trapped QPs continue to grow in a self-sustained manner--hallmarks of the long-anticipated QP-vortex interactions in nonequilibrium superconductivity. We estimate a substantial increase of $\sim$34\% in QP density at vortex densities of $\sim$ 100 magnetic flux quanta per $\mathrm{\mu m^{2}}$. Our findings establish a powerful spectroscopic tool for uncovering QP multiplication and reveal vortex-assisted QP relaxation as a critical materials bottleneck whose mitigation will be essential for resolving QP poisoning and enhancing coherence in superconducting qubits.

We apply the methodology of our recent paper 'The Dynamics of the Hubbard Model through Stochastic Calculus and Girsanov Transformation' [1] to thermodynamic correlation functions in the Fermi-Hubbard model. They can be obtained from a stochastic differential equation (SDE) system. To this SDE system, a Girsanov transformation can be applied. This has the effect that the usual determinant or pfaffian which shows up in a pfaffian quantum Monte Carlo (PfQMC) representation [2] basically gets absorbed into the new integration variables and information from that pfaffian moves into the drift part of the transformed SDE system. While the PfQMC representation depends heavily on the choice of how the quartic interaction has been factorized into quadratic quantities in the beginning, the Girsanov transformed formula has the very remarkable property that it is nearly independent of that choice, the drift part of the transformed SDE system as well as a remaining exponential which has the obvious meaning of energy are always the same and do not depend on the Hubbard-Stratonovich details. The resulting formula may serve as a starting point for further theoretical or numerical investigations. Here we consider the spin-spin correlation at half-filling on a bipartite lattice and obtain an analytical proof that the signs of these correlations have to be of antiferromagnetic type, at arbitrary temperatures. Also, by checking against available benchmark data [3], we find that approximate ground state energies can be obtained from an ODE system. This may even hold for exact ground state energies, but future work would be required to prove or disprove this. As in [1], the methodology is generic and can be applied to arbitrary quantum many body or quantum field theoretical models.

The momentum-space signatures of excitons can be experimentally accessed through time-resolved (pump-probe) photoelectron spectroscopy. In this work, we develop a computational framework for exciton photoemission orbital tomography (exPOT) in periodic systems, enabling the simulation and interpretation of experimental observables within many-body perturbation theory. By connecting the $GW$+Bethe-Salpeter Equation (BSE) approach to photoemission tomography, our formalism captures exciton photoemission in periodic systems, explicitly incorporating photoemission matrix element effects induced by the probe pulse. The correlated nature of electrons and holes introduces distinct consequences for excitonic photoemission, including a dependence on pump pulse polarization. Using the prototypical two-dimensional material hexagonal boron nitride, we demonstrate these effects and show how our framework extends to excitons with finite center-of-mass momentum, making it well-suited to studying momentum-dark excitons. This provides valuable insights into the microscopic nature of excitonic phenomena in quantum materials.