arXiv Curator

The discovery of high-temperature superconductivity in pressurized bulk $\mathrm{La}_{3}\mathrm{Ni}_{2}\mathrm{O}_{7}$ has ignited significant interest in nickelate superconductors. Unlike cuprates, where superconductivity predominantly originates from the $\mathrm{3}d_{x^2-y^2}$ orbital, nickelates exhibit additional complexities involving contributions from the $\mathrm{3}d_{z^2}$ orbital, prompting fundamental questions about their pairing mechanisms. Despite recent progress in stabilizing superconductivity in $\mathrm{La}_{3}\mathrm{Ni}_{2}\mathrm{O}_{7}$ thin films at ambient pressure, direct spectroscopic evidence of the superconducting gap opening remains elusive. Here, we present an in-situ angle-resolved photoemission spectroscopy study of Sr-doped superconducting $\mathrm{La}_{3}\mathrm{Ni}_{2}\mathrm{O}_{7}$ thin films. Fermi surface mapping reveals Ni-$\mathrm{3}d_{x^2-y^2}$-derived $\alpha$ and $\beta$ pockets, with orbital fillings of 0.11$\pm$0.02 electrons and 0.66$\pm$0.03 holes per Ni, respectively, resulting in a total of 0.45$\pm$0.04 electrons for each Ni. These bands exhibit moderate electron correlations, characterized by a band renormalization factor of 3-4. Notably, both $\alpha$ and $\beta$ bands exhibit leading-edge shifts across the superconducting transition, with gap magnitude of ~1-2 meV at Fermi momenta along the Brillouin zone diagonal and slightly away from the zone diagonal, deviating from the conventional $d_{x^2-y^2}$-wave gap structure. Additionally, the Ni-$\mathrm{3}d_{z^2}$-derived $\gamma$ band lies ~75 meV below the Fermi level, indicating a $\mathrm{3}d_{x^2-y^2}$-dominated fermiology in this compound.

An accurate understanding of the ground state electronic properties of La$_3$Ni$_2$O$_7$, a high-temperature superconductor under pressure, is key for unveiling the origin of its superconductivity. In this paper, we conduct a theoretical study of the electronic structure of the bilayer polymorph of La$_3$Ni$_2$O$_7$ using the hybrid functional approach, which is well suited to tackle the non-local correlation effects arising in this system from the molecular orbital splitting of the Ni $3d_{3z^2-r^2}$ states inside Ni-Ni dimers. Our calculations reveal that bilayer La$_3$Ni$_2$O$_7$ is a strongly correlated magnetic system with robust Ni spin moments. Spin moments on individual Ni sites take on unusually small values because of the electron delocalization over molecular orbitals involving multiple Ni and O sites. We further find that the magnetism of bilayer La$_3$Ni$_2$O$_7$ is intimately linked with charge distribution between different Ni and O orbitals. Two distinct regimes are identified in this regard. In one, molecular orbital physics drives the Ni $3d_{x^2-y^2}$ band towards half-filling, which is a well-established condition for unconventional high-temperature superconductivity upon hole doping in cuprates. In the other, the Ni $3d_{x^2-y^2}$ band is quarter-filled favouring spin- and charge-density wave states and Ni-O bond-disproportionation, which is consistent with several recent experimental claims. It is possible that superconductivity in La$_3$Ni$_2$O$_7$ occurs as a result of a pressure-induced transition between these two competing regimes. Since none of the low energy phases discovered in this study are metallic, non-stoichiometry would be required for superconductivity to occur.

We present a real-frequency third-order strong-coupling impurity solver which employs quantics tensor cross interpolation (QTCI) for an efficient evaluation of the diagram weights. Applying the method to dynamical mean-field theory (DMFT) calculations of the single-band Hubbard model on the Bethe lattice, we clarify the interaction and temperature range in which the third-order approach yields accurate results. Since the calculations are implemented on the real-time/frequency axis, the detailed structure of spectral functions can be obtained without analytical continuation, as we demonstrate with examples for paramagnetic, antiferromagnetic and photo-doped states. Our work establishes a viable path toward high-order, real-frequency impurity solvers for both equilibrium and non-equilibrium DMFT studies.

Cuprate superconductors have attracted extensive attention due to high critical temperatures. Conventional cuprates typically contain perfect CuO$_2$ planes which are considered as a key factor to superconductivity since the superconductivity takes place in them. However, in Ba$_2$CuO$_{3+\delta}$ with $\delta=0.2$ and O-depleted CuO$_2$ planes, superconductivity still arises even with a transition temperature as high as 73 K. Using combined density functional theory and dynamical mean-field theory (DFT+DMFT) calculations, we investigated the electronic correlation and electronic structure of Ba$_2$CuO$_{3.25}$ with alternating quasi-one-dimensional (1D) CuO planes and O-depleted CuO$_2$ planes. We find that although different from the usual cuprates, the Cu atoms are still dominated by a 3$d^9$ configuration and the system is of a new kind of correlated single-orbital physics. The quasi-1D CuO planes, composed of parallel Cu-O chains, are slightly hole-doped quasi-1D Mott insulator, while the O-depleted CuO$_2$ planes are more hole doped, with a 2D correlated electronic structure, and may host superconductivity.

We study a square lattice SU(2) gauge theory with link variables $U$, formulated to describe the intermediate temperature pseudogap metal phase of underdoped cuprates of hole density $p$, and its evolution into nodal $d$-wave superconductivity upon cooling. The theory features two fermionic matter sectors: one consists of SU(2)-gauge-neutral, electron-like quasiparticles forming hole pockets; the other is a critical quantum spin liquid composed of electrically neutral spinons, carrying spin $S = 1/2$, Dirac nodes, and transforming in the fundamental representation of the SU(2) gauge group. These sectors interact via a Yukawa coupling to a spinless Higgs boson $B$, which carries electric charge $+e$ and also transforms under the fundamental representation of SU(2). The theory arises from a canonical transformation of the square lattice Hubbard model and a pair of ancilla qubits on each site, along with a specific choice of the spin liquid. We employ a Born-Oppenheimer-type approximation: thermal fluctuations of the bosons $B$ and $U$ are treated classically, while the fermions are fully quantum. Monte Carlo simulations of an effective boson energy functional are consistent with a Kosterlitz-Thouless transition involving $h/(2e)$ vortices, driving the onset of $d$-wave superconductivity. Each vortex is surrounded by a halo of charge order, consistent with scanning tunneling microscopy observations by Hoffman et al. (arXiv:cond-mat/0201348). Above the transition, the electron spectral function exhibits "Fermi arcs", in agreement with ARPES data by Norman et al. (arXiv:cond-mat/9710163) and Shen et al. (Science 307, 901 (2005)). Our model also displays quantum oscillations with area $p/8$ under a magnetic field, aligning with recent evidence via the Yamaji effect by Chan et al. (arXiv:2411.10631).

We study a square lattice SU(2) gauge theory with link variables $U$, formulated to describe the intermediate temperature pseudogap metal phase of underdoped cuprates of hole density $p$, and its evolution into nodal $d$-wave superconductivity upon cooling. The theory features two fermionic matter sectors: one consists of SU(2)-gauge-neutral, electron-like quasiparticles forming hole pockets; the other is a critical quantum spin liquid composed of electrically neutral spinons, carrying spin $S = 1/2$, Dirac nodes, and transforming in the fundamental representation of the SU(2) gauge group. These sectors interact via a Yukawa coupling to a spinless Higgs boson $B$, which carries electric charge $+e$ and also transforms under the fundamental representation of SU(2). The theory arises from a canonical transformation of the square lattice Hubbard model and a pair of ancilla qubits on each site, along with a specific choice of the spin liquid. We employ a Born-Oppenheimer-type approximation: thermal fluctuations of the bosons $B$ and $U$ are treated classically, while the fermions are fully quantum. Monte Carlo simulations of an effective boson energy functional are consistent with a Kosterlitz-Thouless transition involving $h/(2e)$ vortices, driving the onset of $d$-wave superconductivity. Each vortex is surrounded by a halo of charge order, consistent with scanning tunneling microscopy observations by Hoffman et al. (arXiv:cond-mat/0201348). Above the transition, the electron spectral function exhibits "Fermi arcs", in agreement with ARPES data by Norman et al. (arXiv:cond-mat/9710163) and Shen et al. (Science 307, 901 (2005)). Our model also displays quantum oscillations with area $p/8$ under a magnetic field, aligning with recent evidence via the Yamaji effect by Chan et al. (arXiv:2411.10631).

Despite exhibiting a similar $d_{x^2-y^2}$ band character to cuprates, infinite-layer nickelates host additional electron pockets that distinguish them from single-band cuprates. The elusive orbital origin of these electron pockets has led to competing theoretical scenarios. Here, using polarization-dependent and resonant angle-resolved photoemission spectroscopy (ARPES), we determine the orbital character of the Fermi surfaces in NdNiO$_2$. Our data reveal that the electron-like pocket arises predominantly from interstitial $s$ states, with negligible contributions from rare-earth 5$d$ and 4$f$ orbitals near the Fermi level. The observation of well-defined quantum well states indicates a uniform distribution of these interstitial electrons throughout the film thickness. By comparing with electronic structure of LaNiO$_2$, we find that the rare-earth element modulates the Ni-derived bands and hopping integrals through a chemical pressure effect. These findings clarify the role of rare-earth elements in shaping the low-energy electronic structure and establish the presence of electride-like interstitial $s$ states in a correlated oxide system, where electrons occupy lattice voids rather than atomic orbitals. The electride-like character offer new insight into the self-doping and superconductivity in infinite-layer nickelates.

Recently, the films of the Ruddlesden-Popper (RP) nickelate superconductors, in which the (La,Pr)$_3$Ni$_2$O$_7$ system exhibits a remarkable transition temperature $T_c$ exceeding 40 K, were synthesized at ambient pressure. We systematically investigate the band structures and electronic correlation effect to identify the key factors controlling superconductivity and pathways to enhance $T_c$. Based on density functional theory (DFT) calculations, we construct a bilayer two-orbital ($3d_{3z^2-r^2}$ and $3d_{x^2-y^2}$) tight-binding model for a series of in-plane compression mimicking the substrate effect. We find the band energy at the $M$ point drops with the compression, leading to increase of the density of states at the Fermi level, in stark contrast to the behavior of the bulk under pressure. We then apply functional renormalization group (FRG) method to study the electronic correlation effect on the superconductivity. We find the $s_\pm$-wave pairing symmetry remains robust in the films, the same as the bulk. But somewhat surprisingly, for the films, we find $T_c$ can be enhanced by reducing the in-plane lattice constant, increasing the out-of-plane lattice constant, or further electron-doping. These findings are consistent with the itinerant picture of the superconductivity induced by spin-fluctuations and provide theoretical support for further boosting $T_c$ in future experiments.

The recent discovery of high-temperature superconductivity in hole-doped SmNiO$_2$, exhibiting the record-high transition temperature $T_c$ among infinite-layer (IL) nickelates, has opened a new avenue for exploring design principles of superconductivity. Experimentally determining the electronic structure and magnetic interactions in this new system is crucial to elucidating the mechanism behind the enhanced superconductivity. Here, we report a Ni $L$-edge resonant inelastic x-ray scattering (RIXS) study of superconducting Sm-based IL nickelate thin films Sm$_{1-x-y-z}$Eu$_x$Ca$_y$Sr$_z$NiO$_2$ (SECS). Dispersive paramagnonic excitations are observed in both optimally and overdoped SECS samples, supporting a spin-fluctuation-mediated pairing scenario. However, despite the two-fold enhancement of $T_c$ in the Sm-based nickelates compared to their Pr-based counterparts, the effective exchange coupling strength is reduced by approximately $20\%$. This behavior contrasts with hole-doped cuprates, where magnetic interactions correlate positively with $T_c$, highlighting essential differences in their superconducting mechanisms.

The minimal model for infinite-layer nickelates remains under debate, particularly regarding the hybridization between itinerant interstitial-$s$ and the correlated Ni-3$d_{x^2-y^2}$ orbitals, as well as the interaction between $d_{x^2-y^2}$ and other $3d$ orbitals. Additionally, how the doped holes in La$_{1-x}$Sr$_x$NiO$_2$ are distributed among different orbitals remain unresolved. Motivated by recent angle resolved photoemission spectroscopy (ARPES) experiments, we theoretically study the electronic structure of infinite-layer La$_{1-x}$Sr$_x$NiO$_2$ at various doping levels. We find that, unlike the expectation from a rigid band shift, holes are equally distributed to Ni-3$d_{x^2-y^2}$ and interstitial-$s$ orbitals. The role of interstitial-$s$ orbital is further confirmed from the renormalization of Ni-3$d_{x^2-y^2}$ band, for which the coupling between interstitial-$s$ and Ni-3$d_{x^2-y^2}$ exerts a non-negligible impact on the orbital-selective renormalization observed in ARPES. We also discuss the implication of our results to the single-band model, where the interstitial-$s$ orbital in the normal state of La$_{1-x}$Sr$_x$NiO$_2$ acts as charge donator enhancing the correlation of Ni-3$d_{x^2-y^2}$ by increasing its concentration close to half-filling.

Bilayer nickelates exhibit high-temperature superconductivity under proper hydrostatic pressure or epitaxial strain, signifying the emergence of quantum phase transitions whose physical mechanisms remain unclear. Using a minimal bilayer Hubbard model incorporating only the essential Ni-$d_{3z^2-r^2}$ orbitals, we demonstrate that a phase transition naturally arises from tuning the ratio of intralayer to interlayer hopping amplitudes. The transition point separates regimes with a rich interplay between superconducting and density-wave orders. In the regime with weaker intralayer hopping, quasi-long-range spin-density-wave and pair-density-wave orders coexist, with the former being dominant. Across the transition, the spin-density-wave order becomes short-ranged, accompanied by the emergence of the quasi-long-range charge-density-wave order. Most significantly, superconductivity is dramatically enhanced in this regime, though it no longer exhibits the pair-density-wave signature. This quantum phase transition, driven by the competition between intralayer and interlayer hopping, provides a plausible microscopic explanation for the experimentally observed correlation between the superconducting transition temperature and ratio of out-of-plane to in-plane lattice constants. Our findings may assist future efforts in optimizing experimental conditions to further enhance superconductivity in bilayer nickelates.

Near-term quantum processors are limited in terms of the number of qubits and gates they can afford. They nevertheless give unprecedented access to programmable quantum systems that can efficiently, although imperfectly, simulate quantum time evolutions. Dynamical mean field theory, on the other hand, maps strongly-correlated lattice models like the Hubbard model onto simpler, yet still many-body models called impurity models. Its computational bottleneck boils down to investigating the dynamics of the impurity upon addition or removal of one particle. This task is notoriously difficult for classical algorithms, which has warranted the development of specific classical algorithms called "impurity solvers" that work well in some regimes, but still struggle to reach some parameter regimes. In these lecture notes, we introduce the tools and methods of quantum computing that could be used to overcome the limitations of these classical impurity solvers, either in the long term -- with fully quantum algorithms, or in the short term -- with hybrid quantum-classical algorithms.

Motivated by the bilayer cuprate superconductors and nickelate superconductor La$_3$Ni$_2$O$_7$, we investigate the evolution from intralayer to interlayer superconductivity based on a bilayer two-leg $t$-$J$-$J_{\bot}$ model, where $t$ is the in-plane electron hopping, $J$ is the in-plane spin interaction, and $J_{\bot}$ is the inter-plane spin interaction. By means of the density matrix renormalization group calculations, we obtain the quantum phase diagram of the system by tuning $J_{\bot}$ in a large doping range $\delta = 1/8 - 1/2$. We find that a large $J_{\bot}$ can always drive an interlayer superconductivity by coupling the two layers in both the Luther-Emery liquid and Luttinger liquid states. By coupling two Luther-Emery liquid states, the in-plane superconductivity evolves to inter-plane superconductivity either through an intermediate charge density wave (CDW) phase or directly, depending on doping ratio. This emergent CDW phase, which exists over a finite doping range, appears to develop from the CDW state of the two-leg ladder at $\delta = 1/4$. By coupling two Luttinger liquids, the in-plane Luttinger liquids show a transition to the inter-plane superconducting phase at large $J_{\bot}$, as reported in previous literature. Interestingly, in the intermediate $J_{\bot}$ regime we find that while the in-plane Luttinger-liquid features remain stable, the inter-plane superconductivity can develop an enhanced quasi-long-range order with the power exponent $K^{zz}_{\rm SC} \sim 1$. At last, we show that the interlayer superconductivity is also stable by coupling the bilayer three-leg $t$-$J$ ladders by a strong $J_{\bot}$ interaction, from both the Luther-Emery liquid and Luttinger-liquid states.

By partitioning the electronic system of the optimally doped cuprates in two electronic components: (1) mobile electrons on oxygen sub-lattice; and (2) localized spins on copper sub-lattice, and considering the scattering of mobile electrons (on oxygen sub-lattice) via generation of paramagnons in the localized sub-system (copper spins), we ask what should be the electron-paramagnon coupling matrix element $M_q$ so that T-linear resistivity results. This 'reverse engineering approach' leads to $|M_q|^2 \sim \frac{1}{q^2+\xi(T)^{-2}}$. We comment how can such exotic coupling emerge in 2D systems where short range magnetic fluctuations resides. In other words, the role of quantum criticality is found to be crucial. And the T-linear behaviour of resistivity demands that the magnetic correlation length scales as $\xi(T)\propto\frac{1}{T}$, which seems to be a reasonable assumption in the quantum critical regime of cuprates (that is, near optimal doping where T-linear resistivity is observed).

Recent discovery of unconventional superconductivity on multilayer nickelates under high pressure have stimulated great interest. The magnetic and charge configurations in the normal phase of multilayer nickelates which closely relate to the spin-charge fluctuations remain under strongly debated. In this work, we focus the normal-state magnetic, charge and orbital configurations and its evolutions with hydrostatic pressure and Sr-doping concentration of monolayer nickelate La$_2$NiO$_4$. We reveal that in the ambient pressure tetragonal La2NiO4 displays Neel-type antiferromagnetic order in planer Ni spins with negligible interlayer magnetic coupling; with the increase of pressure, La2NiO4 evolves from insulate into metallic state with the critical pressure about P = 50 GPa; with the substitution of La by Sr, the magnetic ground state evolves from G-type, A-type, to double spin striped antiferromagnetic states, as well as weak charge and orbital orders, at x=1; These results can greatly promote our understanding of the magnetic properties of Ruddlesden-Popper nickelates, and shed the light of the pairing mechanism of unconventional nickelate superconductors.

Strong electron correlations drive Mott insulator transitions. Yet, there exists no framework to classify Mott insulators by their degree of correlation. Cuprate superconductors, with their tunable doping and rich phase diagrams, offer a unique platform to investigate the evolution of those interactions. However, spectroscopic access to a clean half-filled Mott-insulating state is lacking in compounds with the highest superconducting onset temperature. To fill this gap, we introduce a pristine, half-filled thallium-based cuprate system, Tl$_2$Ba$_5$Cu$_4$O$_{10+x}$ (Tl2504). Using high-resolution resonant inelastic x-ray scattering (RIXS), we probe long-lived magnon excitations and uncover a pronounced kink in the magnon dispersion, marked by a simultaneous change in group velocity and lifetime broadening. Modeling the dispersion within a Hubbard-Heisenberg approach, we extract the interaction strength and compare it with other cuprate systems. Our results establish a cuprate universal relation between electron-electron interaction and magnon zone-boundary dispersion. Superconductivity seems to be optimal at intermediate correlation strength, suggesting an optimal balance between localization and itinerancy.

Intertwined charge and spin correlations are ubiquitous in a wide range of transition metal oxides and are often perceived as intimately related to unconventional superconductivity. Theoretically envisioned as driven by strong electronic correlations, the intertwined order is usually found to be strongly coupled to the lattice as signaled by pronounced phonon softening. Recently, both charge/spin density waves (CDW/SDW) and superconductivity have been discovered in several Ruddlesden-Popper (RP) nickelates, in particular trilayer nickelates RE4Ni3O10 (RE=Pr, La). The nature of the intertwined order and the role of lattice-charge coupling are at the heart of the debate about these materials. Using inelastic X-ray scattering, we mapped the phonon dispersions in RE4Ni3O10 and found no evidence of phonon softening near the CDW wavevector over a wide temperature range. Calculations of the electronic susceptibility revealed a peak at the observed SDW ordering vector but not at the CDW wavevector. The absence of phonon softening is in sharp contrast to that in canonical oxide materials, notably cuprates. Our experimental and theoretical findings highlight the crucial role of the spin degree of freedom and establish a foundation for understanding the interplay between superconductivity and density-wave transitions in RP nickelate superconductors and beyond.

In copper oxides (cuprates) with single CuO$_2$ layer such as La$_{2-x}$Ba(Sr)$_x$CuO$_4$, antiferromagnetism coexists with superconductivity at small doping levels $x$, where chemical disorders are significant. Here, we report that superconductivity occurs in a uniform and fully ordered N${\'e}$el state in a single-layer cuprate La$_2$CuO$_{4+\delta}$ with a small amount of excess oxygen $(\delta = 0.015)$ as demonstrated by the $^{139}$La nuclear quadrupole resonance measurement. A uniform oxygen distribution in the crystal is crucial for achieving microscopic phase coexistence and overcoming the miscibility gap associated with the staging instability; self-organized periodic oxygen arrangement driven by mobile oxygen atoms. This finding prompts the reconsideration of superconductivity in cuprates, highlighting that it can emerge in a robust N${\'e}$el state that retains sizable magnetic moments and hosts only a small carrier density.

The rapid development of neural quantum states (NQS) has established it as a promising framework for studying quantum many-body systems. In this work, by leveraging the cutting-edge transformer-based architectures and developing highly efficient optimization algorithms, we achieve the state-of-the-art results for the doped two-dimensional (2D) Hubbard model, arguably the minimum model for high-Tc superconductivity. Interestingly, we find different attention heads in the NQS ansatz can directly encode correlations at different scales, making it capable of capturing long-range correlations and entanglements in strongly correlated systems. With these advances, we establish the half-filled stripe in the ground state of 2D Hubbard model with the next nearest neighboring hoppings, consistent with experimental observations in cuprates. Our work establishes NQS as a powerful tool for solving challenging many-fermions systems.

Significant research advances have led to a consensus that the Fermi-Hubbard model and its extended variants are archetypal frameworks for elucidating the intertwined relationship between stripe orders and superconductivity in hole-doped high-$T_c$ materials. Notably, the Hubbard quantum simulator has recently achieved several remarkable breakthroughs, e.g., being successfully cooled down to the cryogenic regime and enabling the observation of stable fluctuating stripes. However, the microscopic mechanism behind $d$-wave pairing of electrons in the presence of stripes at low temperatures remains poorly understood due to the intricate interplay among the strongly correlated effects and non-negligible thermal fluctuations. Here, we conduct a close investigation of a partially-filled stripe in the representative $t$-$J$ and Fermi-Hubbard models with both numerical and analytical methods. Analogous to quantum gas microscopy, the perfect sampling technique allows us to obtain the high-confidence statistics of the Fock basis states appearing in the ground-state wavefunction. In a refreshing physical paradigm, these data demonstrate that two spinons with opposite chiralities tend to spontaneously pair into a singlet state, which naturally gives rise to the $d$-wave pairing pattern. Then, using the effective theory of quantum colored string, we reconstruct the wavefunction and determine the nature of spinon pairing and its connection to the $d$-wave pairing pattern. Furthermore, spinon singlet pairs enable the establishment of a long-range pair-pair correlation between double stripes. Our work offers new insights into the role of stripe orders in mediating $d$-wave superconductivity and paves the way for further exploration of multi-stripe-mediated pairing mechanisms in the Fermi-Hubbard model.

A central challenge for understanding unconventional superconductivity in most strongly correlated electronic materials is their complicated band structure and presence of competing orders. In contrast, quasi-two-dimensional organic spin liquids are single-band systems with superconductivity arising near the bandwidth-tuned Mott metal-insulator transition in the absence of other orders. Here, we study chemically substituted $\kappa$-organics in which superconducting fluctuations emerge in the phase coexistence region between the Mott insulator and the Fermi liquid. Using magnetotransport and ac susceptibility measurements, we find that global superconductivity fails to set in as temperature $T\rightarrow 0$. Our results indicate instead the presence of superconducting domains embedded in the metallic percolating cluster that undergo a magnetic field-tuned quantum superconductor-to-metal phase transition. Surprisingly, albeit consistent with the percolation picture, universal conductance fluctuations are seen at high fields in macroscopic samples. The observed interplay of the intrinsic inhomogeneity and quantum phase fluctuations provides a new insight into failed superconductivity, a phenomenon seen in various conventional and unconventional superconductors, including cuprates.

Non-Hermitian quantum systems, characterized by their ability to model open systems with gain and loss, have unveiled striking phenomena such as the non-Hermitian skin effect (NHSE), where eigenstates localize at boundaries under open boundary conditions. While extensively studied in non-interacting systems, the interplay between NHSE and strong electron correlations remains largely unexplored. Here, we investigate the non-Hermitian Hubbard model with asymmetric hopping, employing real-space dynamical mean-field theory (R-DMFT), a novel extension to such non-Hermitian correlated models-to capture both local correlations and spatial inhomogeneities. By analyzing the end-to-end Green's function as probes of directional amplification, we demonstrate that strong correlations can suppress the skin effect, leading to a crossover from boundary-dominated to correlation-driven dynamics. Our systematic study reveals that correlations dominate at small to intermediate strength of asymmetric hopping, inducing exponential decay in the end-to-end Green's function, but higher strength of asymmetric hopping can restore amplification even at strong interaction. These results illuminate the tunable interplay between correlations and non-Hermitian physics, suggesting avenues for engineering non-reciprocal transport in correlated open quantum systems.

Electronic systems with quadratic band touchings, commonly found in two- and three-dimensional materials such as Bernal-stacked bilayer graphene, kagome metals, HgTe, and pyrochlore iridates, have attracted significant interest concerning the role of interactions in shaping their electronic properties. However, even in the simplest model of spinless fermions on a two-dimensional checkerboard lattice, the quantum phase diagram as a function of nearest-neighbor interaction remains under debate. We employ constrained-path quantum Monte Carlo simulations (CP-QMC) simulations to investigate the problem using a two-dimensional torus geometry. We cross-validate our results on small lattices by comparing them with density-matrix renormalization group calculations, finding quantitative agreement. In particular, we implement an improved optimization scheme within the CP-QMC simulations, enabling the identification of a bond-nematic Dirac semimetal phase that was found in tensor-network studies on cylindrical geometries, but remains inaccessible to Hartree-Fock mean-field methods. The CP-QMC approach makes it possible to establish the emergence of this phase in a geometry that preserves lattice rotational symmetry and permits extrapolation to the thermodynamic limit. Our results show that the quantum phase diagram of spinless fermions on the checkerboard lattice with nearest-neighbor repulsion features three interaction-induced phases at half filling: a quantum anomalous Hall insulator at weak coupling, a bond-nematic Dirac semimetal at intermediate coupling, and a site-nematic insulator at strong coupling.

Recently, a variational wave-function based on the tensor representation of backflow corrections has achieved state-of-the-art energy precision in solving Fermi-Hubbard-type models. However, the Fermi-Hubbard model is very challenging to solve, and the validity of a method relies on investigating the ground state's physical property. For simplicity, we name the tensor representation of backflow corrections as the Tensor-Backflow in this work. We apply the Tensor-Backflow method to investigate the Fermi-Hubbard model on two-dimensional lattices as large as 256 sites, under various interaction strengths $U$, electron fillings $n$ and boundary conditions. Energy precision can be further improved by considering more backflow terms, such as considering backflow terms from next-nearest-neighbours or from all sites. Energy extrapolations on 64-site lattices give competitive results to the gradient optimized fPEPS with the bond dimension as large as $D$=20. For cases of $n$=0.875 and $U$=8 on the $16\times 16$ lattice under open boundary condition, by considering nearest-neighbour backflow terms, obtained energy is only $4.5\times 10^{-4}$ higher than the state-of-the-art method fPEPS with the bond dimension $D$=20. For periodic boundary condition, the variational wave-function is not enforced on any prior symmetry, meanwhile linear stripe order is successfully obtained. Under the same filling and $U$=10,12, energies obtained from initializations with the obtained wave-function for $U$=8 are lower than that from direct optimizations, meanwhile energies are competitive to or even better than state-of-the-art tensor network methods. For cases of $n$=0.8 and 0.9375, results consistent with the phase diagram from AFQMC are obtained by direct optimizations.

Solid matter is classified through symmetry of ordering phenomena. Experimentally, this approach is straightforward, except when distinct orderings occur with identical or almost identical symmetry breaking. Here we show that the cuprate system Y$_{1-x}$Pr$_x$Ba$_2$Cu$_3$O$_{6+y}$ hosts two distinct orderings with almost identical translational symmetry breaking. Only when applying site-sensitive resonant elastic x-ray scattering (REXS), charge ordering can be conclusively differentiated from a super-lattice structure. These two orderings occur with almost the same in-plane symmetry but manifest at different atomic sites and display different temperature dependence. Differentiating these orders provides an important clue to the anomalous behavior of PrBa$_2$Cu$_3$O$_7$ within the 123-series of high-temperature superconductors. We conclude that the symmetry breaking at the Pr-site is unfavorable for superconducting pairing.

We investigate the microscopic mechanisms behind the stabilization of three-dimensional (3D) charge order by Pr doping in YBa$_2$Cu$_3$O$_7$ (YBCO7). Density-functional-theory calculations locate the lowest-energy Pr superlattices for both Ba- and Y-site substitution. In the Ba-site case, the smaller Pr ion pulls the surrounding atoms inward. This breathing-mode distortion pins charge-stripe walls to the Pr columns and forces them to align along the $c$ axis. Y-site Pr is larger than the host ion, produces an outward distortion, and fails to pin the stripes. Coarse-grained Monte-Carlo simulations show that the stripe correlation length rises in step with the structural correlation length of the Pr dopant as observed in prior experiments. We thus identify dopant-induced lattice pinning as the key mechanism behind 3D charge order in Pr-doped YBCO7. This approach provides quantitative guidelines for engineering electronic orders through targeted ionic substitution.

We report a Cu $K$-edge resonant inelastic x-ray scattering (RIXS) study of spin excitations in the $S$=1/2 antiferromagnetic chain system Sr$_2$CuO$_3$. The spectral weight observed below the charge-transfer gap appears in two-spinon continuum, indicating the fractionalization of a spin-conserving ($\Delta S = 0$) magnetic excitation into two-spinon states. The intensity of these excitations reaches a maximum near the midpoint between the center and the boundary of the Brillouin zone, and decreases toward the zone boundary; this behavior contrasts with that of the spin-flip ($\Delta S = 1$) excitations typically observed via inelastic neutron scattering or Cu $L_3$-edge RIXS. The momentum dependence of the intensity is described by the spin-exchange dynamical structure factor. A phenomenological analysis of the symmetry between the polarization and the $d$ orbital explains the resonance condition for the two-spinon excitations.

Motivated by the highest superconducting transition temperature ($T_c$) in multilayer cuprates,we investigated the trilayer Hubbard model by adopting the large-scale dynamical cluster quantum Monte Carlo simulations. Focusing on the systems with hole dopings within the two outer layers (OL) higher than the inner layer (IL), which is believed to be relevant to the realistic multilayer cuprates, our exploration discovered that the IL and OL manifest strong differentiation in a wide range of hole doping combinations. Specifically, the OLs remain metallic while the IL shows a distinct transition from the pseudogap to superconducting state. More importantly, the highest $T_c$ of the composite trilayer system can be largely enhanced compared to the single layer model and the imbalanced hole dopings between IL and OL are generically beneficial for global SC. We further provide strong numerical evidence on the possibility of $d$-wave superconductivity solely hosted in the IL. Our investigation provides new insight into the origin of highest $T_c$ in multilayer cuprates.

We use symmetry arguments to show that the matrix elements of electron-electron interaction on a lattice reach extrema in states composed of wavevectors near high-symmetry points of the Brillouin zone. The mechanism is illustrated by minimal models of cuprates and Fe-based superconductors, where this dependence originates from the wavevector-dependent orbital composition of wavefunctions. We discuss how these dependences can facilitate finite-momentum pairing. Our results provide symmetry-based guidance for the search for new high-temperature superconductors.

The rapid development of neural quantum states (NQS) has established it as a promising framework for studying quantum many-body systems. In this work, by leveraging the cutting-edge transformer-based architectures and developing highly efficient optimization algorithms, we achieve the state-of-the-art results for the doped two-dimensional (2D) Hubbard model, arguably the minimum model for high-Tc superconductivity. Interestingly, we find different attention heads in the NQS ansatz can directly encode correlations at different scales, making it capable of capturing long-range correlations and entanglements in strongly correlated systems. With these advances, we establish the half-filled stripe in the ground state of 2D Hubbard model with the next nearest neighboring hoppings, consistent with experimental observations in cuprates. Our work establishes NQS as a powerful tool for solving challenging many-fermions systems.

Motivated by the highest superconducting transition temperature ($T_c$) in multilayer cuprates,we investigated the trilayer Hubbard model by adopting the large-scale dynamical cluster quantum Monte Carlo simulations. Focusing on the systems with higher hole dopings within the two outer layers (OL) than the inner layer (IL), which is believed to be relevant to the realistic multilayer cuprates, our exploration discovered that the IL and OL manifest strong differentiation in a wide range of hole doping combinations. Specifically, the OLs remain metallic while the IL shows a distinct transition from the pseudogap to superconducting state. More importantly, the highest $T_c$ of the composite trilayer system can be largely enhanced compared to the single layer model and the imbalanced hole dopings between IL and OL is beneficial to the global SC. We further provide strong numerical evidence on the picture that the IL itself can drive the $d$-wave superconductivity while the OLs only serve as the charge reservoir. Our investigation provides new insight into the origin of highest $T_c$ in multilayer cuprates.

Quantum phase transitions are typically marked by changes in quantum correlations across various spatial scales within the system. A key challenge lies in the fact that experimental probes are generally restricted to diagonal quantities at the single-site scale, which are widely believed to be insufficient for detecting phases with off-diagonal long-range order, such as superconducting states. In a striking departure from conventional expectations, we show that single-site diagonal descriptors -- charge and spin fluctuations, occupation probabilities, and entanglement -- can capture the emergence of off-diagonal long-range order in the one-dimensional extended Hubbard model at half-filling. These single-site quantities display clear critical signatures of the superconducting transition, preceded by a continuous breaking of particle-hole symmetry, consistent with a second-order phase transition. While this symmetry breaking has a negligible effect on single-site descriptors, it allows a direct connection between local fluctuations and nonlocal correlations.

The interference patterns of ultracold atoms, observed after ballistic expansion from optical lattices, encode essential information about strongly correlated lattice systems, including phase coherence and non-local correlations. While the interference of lattice bosons has been extensively investigated, quantitative studies of the lattice fermion interference remain challenging. Here, we report the observation and quantitative characterization of interference patterns in low-temperature, homogeneous fermionic Hubbard gases. We develop a novel method to extract first-order correlations from interference patterns, which directly reflect the short-range phase coherence of lattice fermions. Mapping the nearest-neighbor correlations as a function of lattice filling and interaction strength, we observe a crossover from a metal to a Mott insulator. Moreover, at half filling, the measured correlations agree well with quantum Monte Carlo calculations and remain finite in the regime of strong repulsion, revealing virtual tunneling processes driven by quantum fluctuations.

Understanding thermal properties of materials is fundamental to technological applications and to discover new phenomena. In particular, advances of experimental techniques such as cold-atoms measurements allow the simulation of paradigmatic Hamiltonians with great control over model parameters, such as the Fermi-Hubbard model. One aspect of this model which is not much explored is the behavior of the specific heat as function of density. In this work, we perform Determinant Quantum Monte Carlo simulations of the Hubbard model to analyze the specific heat as the filling, interaction and temperature of the system is changed. We found that, with strong correlations, the specific heat presents a three-maxima structure as a function of filling, with local minima between them. This effect can be explained by a decomposition of kinetic and potential contributions to the specific heat. Moreover, by analyzing the kinetic contribution in momentum space we show that, connected to this specific heat behavior, there is a density anomaly detected through the thermal expansion coefficient. These momentum-space quantities are accessible using cold-atoms experiments. Finally, we map the location of these phenomena and connect the thermal expansion anomaly with the well-known Seebeck coefficient change of sign. Our results provide a new perspective to analyze this change of sign. Moreover, our specific heat results demonstrate interesting phenomena away from the commonly studied half-filling regime.

Conductivity is one of the most direct probes of electronic systems, yet its theoretical description remains challenging in the presence of strong non-local correlations. In this Letter, we analyze the conductivity of the half-filled single-band Hubbard model and identify the role of spatial correlations across the Mott transition. We show that in the correlated metallic regime, an accurate description of the conductivity requires not only the correct spectral function but also the inclusion of complex multi-electron processes encoded in vertex corrections. The crossover to the Mott insulating regime is marked by a vanishing contribution of vertex corrections to the DC conductivity. However, in the Mott insulating case, vertex corrections remain significant for the optical conductivity.

Density functional theory is used to investigate the effect of biaxial strain on the structural, electronic and magnetic properties of [111]-oriented NdNiO$_3$, as a representative of the rare-earth perovskites that undergo metal-to-insulator transitions. We find that this constraint on the system induces unique structural phase transitions not previously observed under the well-studied bulk or [001]-oriented strained systems. We also report unique electronic behaviour, including amplification of the electronic band-gap with tensile strain, and insulating, charge-ordered phases with non-orthorhombic tilt patterns. To provide clarity to the trends we observe, we also investigate the coupling between the breathing mode and strain, where we observe certain strains to directly favour and disfavour the creation of the breathing mode (and thus the associated charge-ordering). The amplification of the band gap with strain is understood in terms of a cooperative coupling between the elastic constraint and octahedral breathing, which expands on the previously reported triggered mechanism mediated by octahedral tilting.

In this paper, we investigate the impact of nonlocal correlations on charge fluctuations in the two-dimensional single-band Hubbard model close to the Mott metal-to-insulator transition. At half-filling and for interaction strengths and temperatures where the system is in the Mott insulating phase, charge fluctuations are strongly suppressed. Under these conditions, dynamical mean-field theory (DMFT) calculations predict a strong enhancement of the charge susceptibility at small (electron or hole) doping. However, these DMFT results include only the effects of purely local correlations despite the importance of nonlocal correlations in two-dimensional systems. We have, hence, carried out dynamical vertex approximation (D$\Gamma$A) simulations which allow for the inclusion of such nonlocal correlation effects while retaining the local ones of DMFT. Our numerical data show that close to half-filling the large uniform charge susceptibility of DMFT is strongly suppressed by nonlocal fluctuations but gradually increases with (electron) doping. At a certain doping value charge fluctuations become eventually larger in D$\Gamma$A with respect to DMFT indicating that the absence of nonlocal correlations underestimates the mobility of the charge carriers in this parameter regime. This metallization effect is also reflected in an enhancement of the D$\Gamma$A kinetic and potential energies and a corresponding reduction of the (absolute value of the) D$\Gamma$A Matsubara self-energy with respect to DMFT.

In this paper, we investigate the impact of nonlocal correlations on charge fluctuations in the two-dimensional single-band Hubbard model close to the Mott metal-to-insulator transition. At half-filling and for interaction strengths and temperatures where the system is in the Mott insulating phase, charge fluctuations are strongly suppressed. Under these conditions, dynamical mean-field theory (DMFT) calculations predict a strong enhancement of the charge susceptibility at small (electron or hole) doping. However, these DMFT results include only the effects of purely local correlations despite the importance of nonlocal correlations in two-dimensional systems. We have, hence, carried out dynamical vertex approximation (D$\Gamma$A) simulations which allow for the inclusion of such nonlocal correlation effects while retaining the local ones of DMFT. Our numerical data show that close to half-filling the large uniform charge susceptibility of DMFT is strongly suppressed by nonlocal fluctuations but gradually increases with (electron) doping. At a certain doping value, charge fluctuations eventually become larger in D$\Gamma$A with respect to DMFT indicating that the absence of nonlocal correlations underestimates the mobility of the charge carriers in this parameter regime. This metallization effect is also reflected in an enhancement of the D$\Gamma$A kinetic and potential energies and a corresponding reduction of the (absolute value of the) D$\Gamma$A Matsubara self-energy with respect to DMFT.

We propose a mixed-configuration approximation based on single-band impurity solvers to efficiently study nonequilibrium multi-orbital systems at moderate computational cost. In this work, we merge the approach with the so-called auxiliary master equation approach. As benchmark, we first show that our approach reproduces the results of quantum Monte Carlo (QMC) for two-orbital impurity models at equilibrium with overall good accuracy, especially for non-degenerate orbitals. We then use our approach as impurity solver for dynamical mean-field theory (DMFT) to address the case of a two-orbital, realistic layered structure, recovering the strong crystal-field-driven charge polarization observed by solving the DMFT self-consistent cycle with QMC, albeit slightly reduced. Finally, we address a prototype nonequilibrium setup by sandwiching this layer between metallic contacts subject to a bias voltage described by different chemical potentials. This simplified model demonstrates out method's potential to access nonequilibrium steady-state behavior of multi-orbital, realistic materials. These findings provide a first-step basis for theoretical studies of nonequilibrium properties of multi-orbital compounds directly in the real frequency domain.

Recently, a variational wave-function based on the tensor representation of backflow corrections has achieved state-of-the-art energy precision in solving Fermi-Hubbard-type models. However, the Fermi-Hubbard model is very challenging to solve, and the validity of a method relies on investigating the ground state's physical property. For simplicity, we name the tensor representation of backflow corrections as the Tensor-Backflow in this work. We apply the Tensor-Backflow method to investigate the Fermi-Hubbard model on two-dimensional lattices as large as 256 sites, under various interaction strengths $U$, electron fillings $n$ and boundary conditions. Energy precision can be further improved by considering more backflow terms, such as considering backflow terms from next-nearest-neighbours or from all sites. Energy extrapolations on 64-site lattices give competitive results to the gradient optimized fPEPS with the bond dimension as large as $D$=20. For cases of $n$=0.875 and $U$=8 on the $16\times 16$ lattice under open boundary condition, by considering nearest-neighbour backflow terms, obtained energy is only $4.5\times 10^{-4}$ higher than the state-of-the-art method fPEPS with the bond dimension $D$=20. For periodic boundary condition, the variational wave-function is not enforced on any prior symmetry, meanwhile linear stripe order is successfully obtained. Under the same filling and $U$=10,12, energies obtained from initializations with the obtained wave-function for $U$=8 are lower than that from direct optimizations, meanwhile energies are competitive to or even better than state-of-the-art tensor network methods. For cases of $n$=0.8 and 0.9375, results consistent with the phase diagram from AFQMC are obtained by direct optimizations.

In this work, we extend the theoretical approach known as "D-TRILEX", developed for solving correlated electronic systems, to a cluster reference system for the diagrammatic expansion. This framework allows us to consistently combine the exact treatment of short-range correlation effects within the cluster, with an efficient diagrammatic description of the long-range charge and spin collective fluctuations beyond the cluster. We demonstrate the effectiveness of our approach by applying it to the one-dimensional nano-ring Hubbard model, where the low dimensionality enhances non-local correlations. Our results show that the cluster extension of D-TRILEX accurately reproduces the electronic self-energy at momenta corresponding to the Fermi energy, in good agreement with the numerically exact quantum Monte Carlo solution of the problem, and outperforms significantly more computationally demanding approach based on the parquet approximation. We show that the D-TRILEX diagrammatic extension drastically reduces the periodization ambiguity of cluster quantities when mapping back to the original lattice, compared to cluster dynamical mean-field theory (CDMFT). Furthermore, we identify the CDMFT impurity problem as the main source of the translational-symmetry breaking and propose the computational scheme for improving the starting point for the cluster-diagrammatic expansion.

In multilayered high-Tc cuprates with three or more CuO2 layers per unit cell, the inner CuO2 planes (IPs) are spatially separated from the dopant layers and thus remain cleaner than the outer planes (OPs). While both interlayer coupling and the presence of clean IPs have been proposed as key factors enhancing superconductivity, their individual roles have been difficult to disentangle, as IPs and OPs typically become superconducting simultaneously. Here we investigate five-layer (Cu,C)Ba2Ca4Cu5Oy (Cu1245) with Tc = 78 K and three-layer Ba2Ca2Cu3O6(F,O)2 (F0223) with Tc = 100 K using ARPES, and uncover an unprecedented situation, in which only the IPs become superconducting while the OPs remain metallic at low temperatures. Model calculations indicate that more than 95% of the OP wavefunction remains confined to OP itself, with minimal hybridization from the superconducting IPs. In particular, we experimentally realize an ideal configuration: a single superconducting CuO2 layer sandwiched between heavily overdoped metallic outer layers, which screen disorder originating from the dopant layers. Strikingly, this clean CuO2 layer exhibits the largest superconducting gap among all known cuprates and coherent Bogoliubov peaks extending beyond the antiferromagnetic zone boundary -- long regarded as the boundary beyond which coherence vanishes in heavily underdoped cuprates. Furthermore, a widely extended coherent flat band emerges at the Brillouin zone edge, overcoming the pseudogap damping effect. Our results introduce a new physical parameter, the degree of screening, to investigate the competition between superconductivity and the pseudogap, potentially shedding new light on its origin. The nearly disorder-free superconducting CuO2 layers offer a model platform for bridging the gap between disordered real materials and idealized theoretical models, which generally neglect disorder effects.

LiNiO2 is a promising cathode material for Li-ion battery but its atomic and electronic structure is under debate. Indeed, two sets of Ni-O distances are identified from local structural probes that are related with either Jahn-Teller distortion or bond disproportionation of NiO6 octahedra. Moreover, LiNiO2 undergoes a monoclinic to rhombohedral transition at 200 K which origin is still unclear. On the other hand, isostructural NaNiO2 shows differences from LiNiO2, as it is a well-known Jahn-Teller distorted system, and it undergoes monoclinic to rhombohedral transition at 500 K associated to the loss of the Jahn-Teller distortion. To understand better these differences, we report here Ni L3-edge Resonant inelastic X-ray scattering experiments on LiNiO2 and NaNiO2 at different temperatures (25 to 520 K) and follow the spectral changes below and above the phase transition temperatures. The observed RIXS spectra and the mapping indicate strong spectral changes for NaNiO2 confirming the disappearance of Jahn-Teller distortion during phase transition while the changes are minor for LiNiO2 suggesting very few modifications in the local structure. Theoretical simulations of RIXS spectra are required for further understanding, however, we believe that the reported dataset can be a crucial resource for developing advanced simulations that are essential to deepening our understanding of the atomic and electronic structure of these nickelates.

Recent experimental study on Bi$_2$Sr$_2$CaCu$_2$O$_{8+x}$ superconductors has revealed an unexpected quantum phase transition from superconducting state to insulatinglike state under pressure [Zhou et al., Nat. Phys. 18, 406 (2022)]. To better understand the physical origin of this pressure-induced phenomenon, here we have studied the structural, electronic, and magnetic properties of undoped and O-doped Bi$_2$Sr$_2$CaCu$_2$O$_{8+x}$ (Bi2212) under pressures based on density-functional theory calculations. We first identified the crystal structure of undoped Bi2212 with the armchair distortions in the BiO layers and reproduced the insulating feature of the parent compound. Then we added an extra O atom to the parent compound to simulate the hole-doping effect and found that the structure with O dopant located in the van der Waals (vdW) gap is energetically the most stable. Further calculations on O-doped (0.125 holes/Cu) Bi2212 revealed that the pressure can induce charge redistributions between CuO$_2$ planes and BiO layers; specifically, Cu-$d_{x^2-y^2}$ orbitals gain electrons and Cu atoms rather than O atoms dominate around the Fermi level under high pressure. Along with the increasing pressure, the density of states at the Fermi level first reaches the maximum at $\sim$ 10 GPa and then shows a valley near the Fermi level above 20 GPa, which may be responsible for the insulatinglike state observed in recent experiments. We suggest that the competition among several factors, such as the increase of electrons in the CuO$_2$ plane, the variation of in-plane hopping due to the shortened Cu-O distance, and the enhanced Coulomb repulsion among the Cu-3$d$ electrons, could lead to the exotic transition under pressure. Our work provides an explanation of the high-pressure behaviors of Bi2212, which may facilitate a comprehensive understanding of cuprate superconductors.

We report results of a study of strongly correlated $d$-wave superconducting slabs within a t-J model solved using the statistically consistent Gutzwiller approach. For different dopings, we model slabs cut at 45$^\circ$ relative to the main crystallographic $ab$-axes, i.e., [110] orientation, and conclude that quasiparticle charge is drawn to the edges. Thus, the correlations are locally strengthened, and a region near the edge is closer to the Mott insulating state with reduced hopping amplitudes on the links. Superconductivity is locally weakened near the edges in the correlated state, and the spectral weight of zero-energy Andreev bound states is substantially reduced compared to weak coupling theory. The renormalization of the edge states leaves no room for an extended $s$-wave component to form for hole dopings ranging from strongly underdoped to slightly overdoped.

The Worldvolume Hybrid Monte Carlo (WV-HMC) method [arXiv:2012.08468] is an efficient and low-cost algorithm for addressing the sign problem. It mitigates the sign problem while avoiding the ergodicity issues that are intrinsic to algorithms based on Lefschetz thimbles. In this study, we apply the WV-HMC method to the Hubbard model away from half filling, which is known to suffer from a severe sign problem. We compute the number density on lattices of spatial size $6 \times 6$ and $8 \times 8$ at inverse temperature $\beta = 6.4$ using $N_t = 20$ Trotter steps. Our results show that the WV-HMC method remains effective even in parameter regions where non-thimble Monte Carlo methods fail due to severe sign problems. In this work, we employ direct solvers for fermion matrix inversion, with a computational cost of $O(N^3)$, where $N$ is the number of degrees of freedom and proportional to the spacetime lattice volume. An alternative algorithm employing pseudofermions and iterative solvers, which reduces the cost to $O(N^2)$ at the expense of careful parameter tuning, will be discussed in a separate publication.

Correlated electron physics is intrinsically a multiscale problem, since high-energy electronic states screen the interactions between the correlated electrons close to the Fermi level, thereby reducing the magnitude of the interaction strength and dramatically shortening its range. Thus, the handling of screening is an essential ingredient in the first-principles modelling of correlated electron systems. Screening is an intrinsically dynamic process and the corresponding downfolding methods such as the constrained Random Phase Approximation indeed produce a dynamic interaction. However, many low-energy methods require an instantaneous interaction as input, which makes it necessary to map the fully dynamic interaction to an effective instantaneous interaction strength. It is a priori not clear if and when such an effective model can capture the physics of the one with dynamic interaction and how to best perform the mapping. Here, we provide a systematic benchmark relevant to correlated materials, in the form of the Anderson impurity model. Overall, we find that a moment-based approach recently proposed by Scott and Booth performs well. We also identify physical regimes, especially under doping, where an instantaneous interaction cannot capture all of the relevant physics.

Developing accurate numerical methods for strongly interacting fermions is crucial for improving our understanding of various quantum many-body phenomena, especially unconventional superconductivity. Recently, neural quantum states have emerged as a promising approach for studying correlated fermions, highlighted by the hidden fermion and backflow methods, which use neural networks to model corrections to fermionic quasiparticle orbitals. In this work, we expand these ideas to the space of Pfaffians, a wave-function that naturally expresses superconducting pairings, and propose the hidden fermion Pfaffian state (HFPS), which flexibly represents both unpaired and superconducting phases and scales to large systems with favorable asymptotic complexity. In our numerical experiments, HFPS provides state-of-the-art variational accuracy in different regimes of both the attractive and repulsive Hubbard models. We show that the HFPS is able to capture both s-wave and d-wave pairing, and therefore may be a useful tool for modeling phases with unconventional superconductivity.

Discrete fermionic and bosonic models for hyperbolic lattices have attracted significant attention across a range of fields since the experimental realization of hyperbolic lattices in metamaterial platforms, sparking the development of hyperbolic crystallography. However, a fundamental and experimentally consequential aspect remains unaddressed: fermions propagating in curved space inherently couple to the underlying geometry via the spin connection, as required by general covariance - a feature not yet incorporated in studies of hyperbolic crystals. Here, we introduce a symmetry-based framework for Dirac fermions on two-dimensional hyperbolic lattices, explicitly incorporating spin-curvature coupling via a discrete spin connection. Starting from the continuous symmetries of the Poincar\'e disk, we classify the irreducible representations and construct a symmetry-adapted basis, establishing a direct correspondence to the continuum Dirac theory. We show that this continuum theory predicts a finite density of states at zero energy for any finite curvature in $D-$dimensional hyperbolic space with $2\leq D \leq 4$, suggesting enhanced susceptibility of Dirac fermions to interaction-driven instabilities at weak coupling. We then derive explicit forms of discrete translational and rotational symmetries for lattices characterized by Schl\"afli symbols $\{p,q\}$, and explicitly construct the discrete spin connection, represented as hopping phases, via parallel transport. Our results pave the way for experimental realization of spin-curvature effects in metamaterial platforms and systematic numerical studies of correlated Dirac phases in hyperbolic geometries.

Altermagnetism, a newly discovered magnetic phase, has spurred growing research activity. Studies from a perspective of dynamical electronic correlation still remain scarce. Employing density functional theory plus dynamical mean-field theory (DFT+DMFT) that incorporates dynamical electronic correlation, we demonstrate that CaCrO$_3$ is a strongly correlated altermagnet. Our DFT+DMFT calculations successfully reproduce the correlated metallic behavior of CaCrO$_3$ and quantitatively capture the incoherence state observed experimentally. We also identify that the altermagnetic CaCrO$_3$ is a Hund's metal with a non-Fermi liquid behavior. Moreover, we find that altermagnetism can induce flat bands, and these incipient flat bands are further promoted by the strong renormalization from Hundness, which drives a heavy-fermion behavior. Hence, we propose that Hund's metals provide an ideal platform for hosting novel strongly correlated altermagnetism. Our work will promote the study of strongly correlated altermagnetism physics.

Charge ordering induced by strong short-range repulsion in itinerant fermion systems typically follows a two-sites alternation pattern. However, the covariance matrix spectrum of the one-dimensional, half-filled, spinless $t$-$V$ model reveals a post-Hartree-Fock picture at strong repulsion, with emergent four-site disruptions of the underlying staggered mean-field state. These disruptions are captured in a thermodynamically extensive manner by a compact four-fermion Coupled Cluster (doubles) state (CCS). Remarkably, all properties of this state may be computed analytically by combinatorial means, and also derived from an exactly solvable correlated hopping Hamiltonian. Furthermore, this Coupled Cluster state can be re-expressed as a low-rank Matrix Product State (MPS) with bond dimension exactly four. In addition, we unveil a hidden connection between this Coupled Cluster ansatz and a Rokhsar-Kivelson state (RKS), which is the ground state of a solvable parent quantum tetramer model. The broad picture that we uncover here thus provides deep connections between several core concepts of correlated fermions and quantum chemistry that have previously enjoyed limited synergy. In contrast to a recent perturbative treatment on top of Hartree-Fock theory, our approach asymptotically captures the correct correlations in the $t$-$V$ model at small $t/V$, and remains a qualitatively accurate approximation even outside the perturbative regime. Our results make the case for further studies of the covariance matrix for correlated electron systems in which ground states have non-trivial unit-cell structure.

Identifying model Hamiltonians is a vital step toward creating predictive models of materials. Here, we combine Bayesian optimization with the EDRIXS numerical package to infer Hamiltonian parameters from resonant inelastic X-ray scattering (RIXS) spectra within the single atom approximation. To evaluate the efficacy of our method, we test it on experimental RIXS spectra of NiPS3, NiCl2, Ca3LiOsO6, and Fe2O3, and demonstrate that it can reproduce results obtained from hand-fitted parameters to a precision similar to expert human analysis while providing a more systematic mapping of parameter space. Our work provides a key first step toward solving the inverse scattering problem to extract effective multi-orbital models from information-dense RIXS measurements, which can be applied to a host of quantum materials. We also propose atomic model parameter sets for two materials, Ca3LiOsO6 and Fe2O3, that were previously missing from the literature.

The SU(N) Hubbard model is a natural extension of the SU(2) model. However, even the N=3 case remains poorly understood. We report a substantially new ground-state phase diagram of the square lattice SU(3) Fermi-Hubbard model. Using a backflow ansatz, we identify strong signatures of a Mott transition and a subsequent magnetic transition, and the suppression of a previously predicted magnetic phase. We study the hole-doped model, identifying a transition from magnetic to paramagnetic behavior in the strong-coupling regime. Our findings offer a qualitatively new ground state picture. More broadly, our work suggests a path to study general SU(N) Hubbard models with arbitrary filling and geometry.

The electronic structure of the lightly hole-doped triangular-lattice moir\'e Hubbard model is studied within cluster perturbation theory (CPT) using 13-site clusters for a fixed doping concentration $p=1/13$ varying the Coulomb parameter $U$ and the hopping phase parameter $\phi$ related to the spin-orbital interaction. We have also developed a rather simple generalized mean-field approximation (GMFA) containing the amplitude of the spin correlations as a free parameter to fit the CPT data.The evolution of the Fermi surface and the pseudogap with the parameters $\phi$ and $U$ is explained from the viewpoint of the short-range magnetic order. The geometric frustration and the additional model parameter related to the spin-orbital interaction result in a more rich physics of the pseudogap state compared to the case of a more conventional square lattice.

Neural-network quantum states (NQS) are powerful neural-network ans\"atzes that have emerged as promising tools for studying quantum many-body physics through the lens of the variational principle. These architectures are known to be systematically improvable by increasing the number of parameters. Here we demonstrate an Adaptive scheme to optimize NQSs, through the example of recurrent neural networks (RNN), using a fraction of the computation cost while reducing training fluctuations and improving the quality of variational calculations targeting ground states of prototypical models in one- and two-spatial dimensions. This Adaptive technique reduces the computational cost through training small RNNs and reusing them to initialize larger RNNs. This work opens up the possibility for optimizing graphical processing unit (GPU) resources deployed in large-scale NQS simulations.

The electron in a solid can be considered a bound state of the three independent, fundamental degrees of freedom creating quasi-particles: spinons, carrying the electron spin; plasmons carrying the collective charge mode and orbitons carrying its orbital degree of freedom. These fundamental degrees of freedom could form ordering states in which dynamics or collective motions could occurr and manifest as low-energy excitations. The exotic properties that appear in the materials exhibiting these electronic orderings are associated with these low-energy excitations. Although the orbital order (OO) and its coupling to the spin system creates very interesting phenomena, the microscopic origin of OO has been much less explored than other electronic properties as it is very difficult to directly access experimental information from OO. Due to the recent improvement in energy resolution and flux, soft x-ray resonant inelastic scattering (RIXS) allows for a re-examination of orbital excitations in manganites. Here, we present a study of low energy excitations in half doped A-site ordered SmBaMn_{2}O_{6} through a combination of RIXS and soft x-ray resonant elastic scattering (REXS) measurements. The obtained experimental data confirm the OO at \mathbf{q} = (0.25, 0.25, 0) and find various low energy excitations below 200 meV. while several excitations can be assigned to be of magnetic and phononic origin, a group of excitations between 80 and 200 meV show a temperature dependence distinctively following that of the OO making them possible candidates for orbital excitations.

We explore topological electronic structure of magnetic phases in a multi-orbital Hubbard model with spin-orbit interactions. To account for more general antiferromagnetic orders that go beyond the collinear N\'eel order, two different spin-orbit interactions, Dresselhaus and Rashba spin-orbit interactions, are considered. By performing the canonical transformation, we derive the corresponding generalized t-J model. At half filling, employing self-consistent magnetic order calculations, we find distinctive spin arrangements under Dresselhaus or Rashba spin-orbit interactions. For the Dresselhaus spin-orbit interaction, the spin configuration exhibits collinear antiferromagnetic order. On the other hand, Rashba interaction results in spins antiferromagnetically aligning in xy-plane and a small interaction controlled by hopping parameter induces spin tilting, causing antiferromagnetic alignment in xy-plane but ferromagnetic alignment in z-direction. We categorize topological properties of these phases for low doping in the generalized t-J model.: for 3D collinear antiferromagnetic order, the system possesses a modified time-reversal symmetry, characterized by the Z2 index. In contrast, for systems with tilted antiferromagnetic orders, it is protected by inversion symmetry and characterized by the Z4 index. We further examine the bulk-edge correspondence for non-collinear magnetic phases, revealing that the surface state becomes gapless when the surface is parallel to the ferromagnetic component of tilted antiferromagnetic order; otherwise, the surface state exhibits a gap. Our findings offer a comprehensive topological characterization for doped and canted antiferromagnetic insulators with spin-orbit interactions, providing valuable insights into the interplay between spin arrangements, symmetries, and topological properties in systems governed by the multi-orbital Hubbard model.

Altermagnetism (AM), the recently discovered third class of collinear magnetic order, is characterized by non-relativistic momentum-dependent spin-split electronic structure with compensated zero net magnetization. It can arise from the conventional antiferromagnetism by introducing local anisotropy on the two opposite-spin sublattices, either through structural changes in local crystallographic symmetry or spontaneous emergence of local staggered orbital order from electron correlations in multi-orbital systems. Here, we demonstrate on the two-dimensional square lattice that a $d$-wave AM can emerge spontaneously in the single-orbital extended Hubbard model, without invoking the spin-orbital coupling and multi-orbital physics. We carry out mean-field studies on the concrete single-orbital $t$-$U$-$V$ model with $U$ and $V$ the onsite and nearest-neighbor Coulomb interactions, obtaining the ground states, analyzing their properties, and determining the phase diagram in the $U$-$V$ plane. The $d$-wave AM with novel spin-transport behavior is found to be stabilized in a wide region of the phase diagram when the system is doped away from half-filling, actualized by the coexistence of onsite antiferromagnetic order and complex $d$-wave nearest-neighbor spin bond orders. Our findings provide an alternative route to achieve AM and substantially expand the range of candidate AM materials.

Motivated by recent developments in the study of quasi-one-dimensional iron systems with Fe$^{2+}$, we comprehensively study the Fe$^{3+}$ chalcogenide chain system. Based on first-principles calculations, the Fe$^{3+}$ chain has a similar electronic structure as discussed before in the iron 2+ chain, due to similar Fe$X_4$ ($X$ = S or Se) tetrahedron chain geometry. Furthermore, a three-orbital electronic Hubbard model for this chain was constructed by using the density matrix renormalization group method. A robust antiferromagnetic coupling was unveiled in the chain direction. In addition, in the intermediate electronic correlation $U/W$ region, we found an interesting orbital-selective Mott phase with the coexistence of localized and itinerant electrons ($U$ is the on-site Hubbard repulsion, while $W$ is the electronic bandwidth). Furthermore, we do not observe any obvious pairing tendency in the Fe$^{3+}$ chain in the electronic correlation $U/W$ region, where superconducting pairing tendencies were reported before in iron ladders. This suggests that superconductivity is unlikely to emerge in the Fe$^{3+}$ systems. Our results establish with clarity the similarities and differences between Fe$^{2+}$and Fe$^{3+}$ iron chains, as well as iron ladders.

The Grassmann time-evolving matrix product operator method has shown great potential as a general-purpose quantum impurity solver, as its numerical errors can be well-controlled and it is flexible to be applied on both the imaginary- and real-time axis. However, a major limitation of it is that its computational cost grows exponentially with the number of impurity flavors. In this work, we propose a multi-flavor extension of it to overcome this limitation. The key insight is that to calculate multi-time correlation functions on one or a few impurity flavors, one could integrate out the degrees of freedom of the rest flavors before hand, which could greatly simplify the calculation. The idea is particularly effective for quantum impurity problems with diagonal hybridization function, i.e., each impurity flavor is coupled to an independent bath, a setting which is commonly used in the field. We demonstrate the accuracy and scalability of our method for the imaginary time evolution of impurity problems with up to three impurity orbitals, i.e., 6 flavors, and benchmark our results against continuous-time quantum Monte Carlo calculations. Our method paves the way of scaling up tensor network algorithms to solve large-scale quantum impurity problems.

We systematically investigate the stability of the symmetry-protected topological (SPT) Haldane phase in spin-1/2 Heisenberg and half-filled Hubbard ladders coupled by ferromagnetic Hund's interactions. Using density-matrix renormalization group (DMRG) method, we analyze key signatures of the Haldane phase: long-range string order, finite spin gap, and characteristic entanglement spectrum degeneracies. In spin-only Heisenberg ladders, we find immediate onset and continuous strengthening of the Haldane phase with increasing Hund's coupling. In contrast, the inclusion of charge fluctuations in Hubbard ladders leads to a nontrivial stability regime, revealing a robust yet bounded region where SPT order persists despite significant charge fluctuations. We identify distinct boundaries separating a trivial insulating phase from the Haldane SPT phase, governed by both Coulomb repulsion and Hund's coupling. Our results highlight the subtle interplay of spin and charge degrees of freedom in correlated itinerant systems and establish essential criteria for observing Haldane physics experimentally in fermionic ladder materials.

The Self-Learning Monte Carlo (SLMC) method is a Monte Carlo approach that has emerged in recent years by integrating concepts from machine learning with conventional Monte Carlo techniques. Designed to accelerate the numerical study of interacting many-body systems, SLMC significantly improves sampling efficiency by constructing an effective model -- via machine learning methods -- based on configurations generated by conventional Monte Carlo methods and then proposes global updates based on the effective model. This enhancement leads to a substantial reduction in autocorrelation time, especially near the critical region, where traditional methods typically suffer from critical slowing down and increased computation complexity. Moreover, SLMC maintains statistical accuracy by implementing a cumulative update scheme that rigorously satisfies the detailed balance condition. And more recent applications have extended the SLMC to convolutional neural networks with applications not only in condensed matter physics but also high-energy physics, quantum chemistry, and quantum simulations. The generic applicability and high computational efficiency make SLMC a powerful and scalable framework for quantum Monte Carlo simulations of strongly correlated electron systems, extending the reach of numerical investigations beyond the limitations of conventional techniques.

Superconductivity arising from fully spin-polarized, repulsively interacting electrons can host intrinsically chiral Cooper pairs and Majorana zero modes, yet no concrete microscopic route to such a state has been established. Motivated by recent observations in twisted MoTe$_2$ and rhombohedral pentalayer graphene, where fractional Chern insulators (FCIs) appear adjacent to spin-valley-polarized superconductors, we investigate a minimal model: spinless electrons in the lowest Landau level subject to a tunable moire potential. Large-scale density-matrix renormalization group (DMRG) calculations show that, as the FCI gap closes, two nearly degenerate phases emerge before a metallic state forms. A chiral $f+if$ superconductor and a commensurate $\sqrt{3} \times \sqrt{3}$ charge-density wave (CDW) differ in energy by less than $1\%$, reproducing the distinct superconducting and re-entrant integer quantum Hall (RIQH) phases experimentally observed near the FCI regime. The superconducting dome remains robust against realistic Coulomb screening, light doping, and variations in lattice geometry. Melting the FCI therefore provides a new mechanism for realizing spin-polarized chiral superconductivity and RIQH order. We predict that twisted MoTe$_2$ at larger twist angles will develop a superconducting dome even at filling $\nu = 2/3$, and suppressing this superconductivity with a magnetic field should drive the system into an RIQH state.

The motion of dopants in magnetic spin lattices has received tremendous attention for at least four decades due to its connection to high-temperature superconductivity. Despite these efforts, we lack a complete understanding of their behavior, especially out-of-equilibrium and at nonzero temperatures. In this Article, we take a significant step towards a much deeper understanding based on state-of-the-art matrix-product-state calculations. In particular, we investigate the non-equilibrium dynamics of a dopant in two-leg $t$--$J$ ladders with antiferromagnetic XXZ spin interactions. In the Ising limit, we find that the dopant is \emph{localized} for all investigated \emph{nonzero} temperatures due to an emergent disordered potential, with a localization length controlled by the underlying correlation length of the spin lattice, whereby it only delocalizes asymptotically in the zero temperature limit. This greatly generalizes the localization effect discovered recently in Hilbert space fragmented models. In the presence of spin-exchange processes, the dopant delocalizes according to a power-law behavior, which is strongly sub-diffusive for weak spin-exchange but which eventually becomes diffusive for strong enough exchange. Moreover, we show that the underlying spin dynamics at infinite temperature behaves qualitatively the same, albeit with important quantitative differences. We substantiate these findings by showing that the dynamics shows self-similar scaling behavior, which strongly deviates from the Gaussian behavior of regular diffusion. Finally, we show that the diffusion coefficient follows an Arrhenius relation at high temperatures, whereby it exponentially decreases for decreasing temperatures.

We construct a dual symmetry-protected topological (SPT) Hamiltonian for the $U(1)$ symmetric anisotropic spin-$\frac{1}{2}$ Heisenberg chain-a model that has traditionally been used to study spontaneous symmetry breaking (SSB) in both ferromagnetic and antiferromagnetic phases, with an intervening extended Luttinger liquid phase. By performing a non-local unitary transformation, we explicitly construct a local fermionic Hamiltonian that exhibits two nontrivial topological phases separated by an extended Luttinger liquid regime. We demonstrate the topological nature of these phases by analyzing the entanglement structure, deriving a non-local string order parameter, and constructing an exact zero mode operator that connects states in different fermionic parity sectors.

Low-symmetry two-dimensional (2D) topological materials such as MoTe$_2$ host efficient charge-to-spin conversion (CSC) mechanisms that can be harnessed for novel electronic and spintronic devices. However, the nature of the various CSC mechanisms and their correlation with underlying crystal symmetries remain unsettled. In this work, we use local spin-sensitive electrochemical potential measurements to directly probe the spatially dependent nonequilibrium spin accumulation in MoTe$_2$ flakes down to four atomic layers. We are able to clearly disentangle contributions originating from the spin Hall and Rashba-Edelstein effects and uncover an abundance of unconventional spin polarizations that develop uniquely in the sample bulk and edges with decreasing thickness. Using ab-initio calculations, we construct a unified understanding of all the observed CSC components in relation to the material dimensionality and stacking arrangement. Our findings not only illuminate previous CSC results on MoTe$_2$ but also have important ramifications for future devices that can exploit the local and layer-dependent spin properties of this 2D topological material.

Recently, radical progress in machine learning (ML) has revolutionized computational materials science, enabling unprecedentedly rapid materials discovery and property prediction, but the quantum many-body problem -- which is the key to understanding excited-state properties, ranging from transport to optics -- remains challenging due to the complexity of the nonlocal and energy-dependent interactions. Here, we propose a symmetry-aware, grid-free, transformer-based model, MBFormer, that is designed to learn the entire many-body hierarchy directly from mean-field inputs, exploiting the attention mechanism to accurately capture many-body correlations between mean-field states. As proof of principle, we demonstrate the capability of MBFormer in predicting results based on the GW plus Bethe Salpeter equation (GW-BSE) formalism, including quasiparticle energies, exciton energies, exciton oscillator strengths, and exciton wavefunction distribution. Our model is trained on a dataset of 721 two-dimensional materials from the C2DB database, achieving state-of-the-art performance with a low prediction mean absolute error (MAE) on the order of 0.1-0.2 eV for state-level quasiparticle and exciton energies across different materials. Moreover, we show explicitly that the attention mechanism plays a crucial role in capturing many-body correlations. Our framework provides an end-to-end platform from ground states to general many-body prediction in real materials, which could serve as a foundation model for computational materials science.

We perform the variational Monte Carlo calculation for recently proposed chiral superconducting states driven by strong Coulomb interactions. We compare the resulting energetics of these electronic phases for the electron dispersion relation $E_k = c_2 k^2+c_4 k^4$. Motivated by the recent discovery of chiral superconductivity in rhombohedral graphene systems, we apply our analysis to relevant parameter regimes. We demonstrate that topological chiral superconducting phases (including a spin-unpolarized state) can be energetically favored over the spin-valley polarized Fermi liquid above the density of Wigner crystal phase. Our results show that the preference for chiral superconductivity is strongest when $c_2$ lies between zero and a negative value, corresponding to a system on the verge of forming a hole pocket around $k=0$. This finding suggests that superconductivity can arise from pure repulsive Coulomb interactions in systems with an almost flat band bottom, without relying on the pairing instability of a Fermi surface. This mechanism opens a new pathway to superconductivity beyond the conventional BCS mechanism.

Altermagnets present a new class of fully compensated collinear magnetic order, where the two sublattices are not related merely by time-reversal combined with lattice translation or inversion, but require an additional lattice rotation. This distinctive symmetry leads to a characteristic splitting of the magnon bands; however the splitting is only partial -- residual degeneracies persist along certain lines in the Brillouin zone as a consequence of the underlying altermagnetic rotation. We consider a two-dimensional $d$-wave altermagnetic spin model on the checkerboard lattice and introduce additional interactions such as an external magnetic field and Dzyaloshinskii-Moriya interactions, that lift these degeneracies. The resulting magnon bands become fully gapped and acquire non-trivial topology, characterized by nonzero Chern numbers. We demonstrate the crucial role of altermagnetism for the generation of the Berry curvature. As a direct consequence of the topological magnons, we find finite thermal Hall conductivity $\kappa_{xy}$, which exhibits a characteristic low-temperature scaling, $\kappa_{xy}\propto T^4$. Moreover, $\kappa_{xy}$ changes sign under reversal of the magnetic field, exhibiting a sharp jump across zero field at low temperatures. We also demonstrate topologically protected chiral edge modes in a finite strip geometry.

Non-equilibrium control of electronic properties in condensed matter systems can result in novel phenomena. In this work, we provide a novel non-equilibrium route to realize half-metallic phases. We explore the periodically driven Hubbard model on a bipartite lattice and demonstrate that a periodic drive can transform a weakly interacting metal into a ferrimagnetic half-metal. We consider a Fermi-Hubbard model with only nearest-neighbour hopping and stabilize the elusive phase simply by driving the site potentials periodically. The drive induces staggered second and third-neighbor hopping and a staggered potential between two sublattices in the Floquet Hamiltonian, whose ground state is explored in this work. Close to the dynamical freezing point, due to the suppression of nearest neighbor hopping in the driven system, an effective enhancement of various terms in the Floquet Hamiltonian, including the e-e interactions, occurs. This helps in stabilizing a broad ferrimagnetic half-metallic phase for a wide range of system parameters. The half-metallic phase achieved in the presence of high drive frequency should be stable for exponentially large time scales in drive frequency and could be perpetually stable beyond a strong enough drive amplitude owing to dynamical freezing. It can hence have potential applications in stable spintronics and other upcoming quantum technologies.

Motivated by experiments in which the quantum anomalous Hall (QAH) charge gap greatly exceeds the Curie temperature, we apply determinantal quantum Monte Carlo to two complementary lattice models with different geometries: the checkerboard quantum-spin-Hall-Hubbard model and the generalized Kane-Mele-Hubbard model. In both cases an incompressible QAH phase with total Chern number $C=\pm 1$ emerges at quarter filling well above the Curie temperature. Each spin channel carries its own nonquantized Chern number while remaining only partly filled, revealing a topological Mott insulator where Mottness opens the charge gap before magnetic order appears. Charge excitations bind triplet excitons that suppress net magnetization, a many-body dynamical effect absent in mean-field theory. The concurrence of these results on two very different models shows that coupling Mottness with band topology generically yields a high-temperature QAH insulator whose charge excitations are dressed by triplet excitons.

Motivated by experiments in which the quantum anomalous Hall (QAH) charge gap greatly exceeds the Curie temperature, we apply determinantal quantum Monte Carlo to two complementary lattice models with different geometries: the checkerboard quantum-spin-Hall-Hubbard model and the generalized Kane-Mele-Hubbard model. In both cases an incompressible QAH phase with total Chern number $C=\pm 1$ emerges at quarter filling well above the Curie temperature. Each spin channel carries its own nonquantized Chern number while remaining only partly filled, revealing a topological Mott insulator where Mottness opens the charge gap before magnetic order appears. Charge excitations bind triplet excitons that suppress net magnetization, a many-body dynamical effect absent in mean-field theory. The concurrence of these results on two very different models shows that coupling Mottness with band topology generically yields a high-temperature QAH insulator whose charge excitations are dressed by triplet excitons.

We present models of topological insulating Hamiltonians exhibiting intrinsic altermagnetic features, protected by combined three-fold or four-fold rotational symmetries with time-reversal. We demonstrate that the spin Chern number serves as a robust topological invariant in two-dimensional systems, while for three-dimensional structures, the topological nature is characterized by the spin Chern numbers computed on the $k_z$=$0$ and $k_z$=$\pi$ planes. The resulting phases support symmetry-protected boundary modes, including corner, hinges and surface states, whose structure is determined by the magnetic symmetry and the local magnetic moments. Our findings bridge the fields of altermagnetism and topological quantum matter, and establish a theoretical framework for engineering spintronic topological systems without net magnetization.

The combination of a centrosymmetric crystallographic structure with local structural alternations and collinear antiferromagnetism can lead to broken PT (Parity $\times$ Time-reversal) symmetry, resulting in altermagnets with non-relativistic spin-split bands. The 6H perovskites with composition A$_3$BB'$_2$O$_9$ exhibit unique layered structural alternations and typically adopt an antiferromagnetic ground state. Here, we report the discovery that several 6H perovskites are indeed altermagnets exhibiting non-relativistic spin-split bands. We also explore the possible presence of net magnetization due to spin-orbit coupling in these materials, as well as the manifestation of giant piezomagnetism. Since the single crystals of 6H perovskites can be readily grown and cleavable, our findings provide a new avenue to study the cleaved atomically-flat surfaces of altermagnets with advanced experimental techniques such as spin-resolved scanning tunneling microscopy (STM) or spin-resolved angle-resolved photoemission spectroscopy (ARPES) to explore their spin splitting nature.

We introduce an efficient approach to implement neural network quantum states (NNQS) as trial wavefunctions in auxiliary-field quantum Monte Carlo (AFQMC). NNQS are a recently developed class of variational ans\"atze capable of flexibly representing many-body wavefunctions, though they often incur a high computational cost during optimization. AFQMC, on the other hand, is a powerful stochastic projector approach for ground-state calculations, but it normally requires an approximate constraint via a trial wavefunction or trial density matrix, whose quality affects the accuracy. Recent progress has shown that a broad class of highly correlated wavefunctions can be integrated into AFQMC through stochastic sampling techniques. In this work, we present a direct integration of NNQS with AFQMC, allowing NNQS to serve as high-quality trial wavefunctions for AFQMC with manageable computational cost. We test the NNQS-AFQMC method on the challenging nitrogen molecule (N$_2$) at stretched geometries. Our results demonstrate that AFQMC with an NNQS trial wavefunction can attain near-exact total energies, highlighting the potential of AFQMC with NNQS to overcome longstanding challenges in strongly correlated electronic structure calculations. We also outline future research directions for improving this promising methodology.

Recently, the study of novel physical properties arising from the combination of altermagnetism and other matter phases has attracted widespread attention, such as the integration of altermagnetism and topology. However, research on the combination of altermagnetism and charge density waves remains relatively sparse. In this letter, based on symmetry analysis and first-principles calculations, we demonstrate for the first time that altermagnetism and charge density waves can coexist in a two-dimensional material and predict monolayer WO to be such a material. Moreover, our calculations reveal that the altermagnetic order in monolayer WO stabilizes the $\sqrt{2}\times\sqrt{2}$ charge density wave. Further, the $\sqrt{2}\times\sqrt{2}$ charge density wave is not driven by Fermi-surface nesting but rather by strong electron-phonon coupling. More importantly, the $\sqrt{2}\times\sqrt{2}$ charge density wave in monolayer WO leads to an anomalous transition from semimetal to metal. Therefore, we realize an anomalous charge density wave phase in altermagnetic WO. Considering the strong electron-phonon coupling and good metallic properties in the altermagnetic charge density wave state, our work may provide new insights into the realization of nontrivial altermagnetic superconductivity.

We consider BCS-like states in the 2D Hubbard model which, starting from some arbitrary point in field space in the neighborhood of a Hartree-Fock ground state, are relaxed within that BCS ansatz to local minima of the energy. As in the Hartree-Fock approximation there are a vast number of local minima. What is new, and unlike the conventional Hartree-Fock states, is that there is a region in parameter (coupling-density) space, where these local minima are associated with d-wave condensates.

Solving the many-electron problem, even approximately, is one of the most challenging and simultaneously most important problems in contemporary condensed matter physics with various connections to other fields. The standard approach is to follow a divide and conquer strategy that combines various numerical and analytical techniques. A crucial step in this strategy is the derivation of an effective model for a subset of degrees of freedom by a procedure called downfolding, which often corresponds to integrating out energy scales far away from the Fermi level. In this work we present a rigorous formulation of this downfolding procedure, which complements the renormalization group picture put forward by Honerkamp [PRB 85, 195129 (2012)}]. We derive an exact effective model in an arbitrarily chosen target space (e.g. low-energy degrees of freedom) by explicitly integrating out the the rest space (e.g. high-energy degrees of freedom). Within this formalism we state conditions that justify a perturbative truncation of the downfolded effective interactions to just a few low-order terms. Furthermore, we utilize the exact formalism to formally derive the widely used constrained random phase approximation (cRPA), uncovering underlying approximations and highlighting relevant corrections in the process. Lastly, we detail different contributions in the material examples of fcc Nickel and the infinite-layer cuprate SrCuO$_2$. Our results open up a new pathway to obtain effective models in a controlled fashion and to judge whether a chosen target space is suitable.

Quantum phase transitions between topologically ordered and symmetry-broken phases lie beyond Landau theory. A prime example is the conjectured continuous transition from the bosonic $\nu = 1/2$ Laughlin state to a superfluid, proposed to be governed by a QED$_3$--Chern--Simons (CS) critical point whose stability remains uncertain. We study half-filled bosons in the lowest Landau level subject to a lattice potential. Infinite-cylinder DMRG reveals a single continuous Laughlin--to--superfluid transition. Adiabatic flux insertion collapses the many-body gap and exposes massless Dirac quasiparticles, while momentum-resolved correlation lengths show that three lattice-related density modes share the same critical exponent, evidencing an emergent $SO(3)$ symmetry. The joint appearance of Dirac dispersion and symmetry enlargement provides microscopic support for a stable QED$_3$--CS fixed point. Our numerical strategy also offers a blueprint for exploring Landau-forbidden transitions in fractional Chern insulators and composite Fermi liquids realised in moire and cold-atom systems.

Superlight pairing of bipolarons driven by electron-phonon interactions (EPIs) in face-center-cubic (FCC) and body-center-cubic (BCC) lattices is investigated using a continuous-time path-integral quantum Monte Carlo (QMC) algorithm. The EPIs are of the Holstein and extended Holstein types, and a Hubbard interaction is also included. The number of phonons associated with the bipolaron, inverse mass, and radius are calculated and used to construct a phase diagram for bipolaron pairing (identifying the regions of pairing into intersite bipolarons and onsite bipolarons). From the inverse mass it is determined that for the extended interaction, there is a region of light pairing associated with intersite bipolarons formed in both BCC and FCC lattices. Intersite bipolarons in the extended model at large phonon frequency and large Coulomb repulsion become superlight due to first order hopping effects. The transition temperature of Bose--Einstein condensates of these pairs is estimated. It is determined that intersite bipolarons are associated with regions of high transition temperatures.

We have recently presented the multichannel Dyson equation as an alternative to the standard single-channel Dyson equation. While the latter involves a single many-body Green's function, the former uses a multichannel Green's function in which two or more many-body Green's functions are coupled. Quasiparticles and satellites are thus naturally treated on equal footing in the multichannel Dyson equation. To assess the accuracy of our approach we apply it here to the ground- and excited-state properties of the extended Hubbard dimer, an exactly solvable model for $H_2$. In particular, we focus on the potential energy surface as well as the corresponding spectral functions and HOMO-LUMO gaps, which are well-known challenges for many-body approximations such as second Born and $GW$. We show that the multichannel Dyson equation gives overall very good results for all properties considered and outperforms both $GW$ and second Born. In particular, the multichannel Dyson equation yields the correct ground-state energy and HOMO-LUMO gap in the dissociation limit contrary to $GW$.

1T-TiSe$_2$ exhibits a charge-density-wave (CDW) transition below 200 K, which is believed to be driven by a hybrid exciton-phonon mechanism, making it a versatile platform for investigating the interplay between electronic and lattice degrees of freedom. Although the corresponding band structure modifications below the CDW transition temperature are well established, only a few reports discuss the occurrence of the CDW fluctuating phase and its spectral features well above the transition temperature. Here, we report the direct observation of spectral features associated with CDW fluctuations at 295 K, using time-resolved extreme ultraviolet momentum microscopy. We investigated the transient melting and recovery of CDW fluctuations upon nonresonant ultrafast photoexcitation. Surprisingly, our results reveal that the coherent amplitude mode modulating the ultrafast CDW recovery persists at these elevated temperatures. The time-, energy- and momentum-resolved photoemission data supported by density functional perturbation theory further confirm that CDW fluctuations at these elevated temperatures are dominated by electron-phonon interaction. The analysis of these very localized microscopic fluctuations consequently provides new insights into the complex interplay between the electronic and lattice degrees of freedom at elevated temperatures and, therefore, on the nature of this quantum phase transition.

We investigate the nonequilibrium dynamics of the transverse field Ising chain (TFIC) using time-resolved resonant inelastic X-ray scattering (tr-RIXS) in a pump-probe setup within a condensed matter setting. The tr-RIXS spectra reveal distinctive features, including two-kink continua at high energies and oscillatory spectral weight at low energies, both strongly influenced by the pump dynamics. By systematically analyzing the nonequilibrium tr-RIXS cross section under various pump configurations, we demonstrate how the spectra encode detailed information about the system's transient states, capturing instantaneous Hamiltonian parameters and the traversal of equilibrium quantum critical points during dynamic evolution. Notably, we uncover a one-to-one correspondence between oscillatory features in the low-energy spectra and dynamical quantum phase transitions, offering a novel and experimentally accessible approach for their detection. These findings highlight the versatility of tr-RIXS as a powerful tool for studying quantum systems far from equilibrium and their critical phenomena.

We present a quantum Monte Carlo-based approach to detect and compute the most dominant correlations for many-body systems without prior knowledge. It is based on the measurement and analysis of the correlation density matrix between two (small) subsystems embedded in the full (large) sample. In order to benchmark this procedure, we investigate zero-temperature quantum phase transitions in one- and two-dimensional quantum Ising model as well as the two-dimensional bilayer Heisenberg antiferromagnet. The method paves the way for a systematic identification of unknown or exotic order parameters in unexplored phases on large systems accessible to quantum Monte Carlo methods.

The Hubbard model at $U\to\infty$ has recently been shown to have resonant valence bond (RVB) ground states on the corner-sharing sawtooth and pyrochlore lattices in the dilute doping limit of a single vacancy. In an effort to further generalize those results, I study how the ground state is modified when not all corners are shared between two tetrahedra as in the quasi-1D lattices of a pyrochlore stripe, and how to approach the problem in the case of finite doping. Using a non-Abelian version of the flux inequality, the tetrahedron chain is shown to have degenerate RVB-like ground states. The Bethe ansatz (BA) is adapted to solve the sawtooth chain with spinless or spin-polarized fermions and multiple holons, which is the first example of applying BA to a quasi-1D lattice.

A key to understand how electrons behave in crystalline solids is the band structure that connects the energy of electron waves to their wavenumber (k). Even in the phase of matter with only short-range order (liquid or amorphous solid), the coherent part of electron waves still possesses a band structure. Theoretical models for the band structure of liquid metals were formulated more than 5 decades ago, but thus far, bandstructure renormalization and pseudogap induced by resonance scattering have remained unobserved. Here, we report the observation of this unusual band structure at the interface of a crystalline insulator (black phosphorus) and disordered dopants (alkali metals). We find that a conventional parabolic band structure of free electrons bends back towards zero k with the pseudogap of 30-240 meV from the Fermi level. This is k renormalization caused by resonance scattering that leads to the formation of quasi-bound states in the scattering potential of alkali-metal ions. The depth of this potential tuned by different kinds of alkali metal (Na, K, Rb, and Cs) allows to classify the pseudogap of p-wave and d-wave resonance. Our results may provide a clue to the puzzling spectrum of various crystalline insulators doped by disordered dopants, such as the waterfall dispersion in cuprates.

The GdFeO$_3$-type perovskite antiferromagnet NaOsO$_3$, calculated here to be altermagnetic for all three typical collinear antiferromagnetic orders, was suggested early on to be a Slater-type insulator, due in large part to its continuous metal-insulator transition and its small energy gap. Below the N\'eel temperature, the gap opens along with ``weak magnetism'', accompanied by a discontinuity in the magnetic susceptibility. Most other properties remain continuous. Without explicit correlation in the band structure calculation, and neglecting spin-orbit coupling (SOC), already a gap opens. Inclusion of a modest on-site Coulomb repulsion ($U\sim$ 1 eV) is sufficient to eliminate a SOC-induced small band overlap, reproducing the experimentally observed gap of several tens meV as reported earlier. This evidence supports the viewpoint that NaOsO$_3$ lies in an unusual crossover region between Slater and Mott insulator. The unreported altermagnetism in NaOsO$_3$ is demonstrated and its consequences are considered. The origin of the very weak magnetism has been investigated using a combination of {\it ab initio} calculations and symmetry analysis of the magnetic space group, confirming the origin lying in the Dzyaloshinskii-Moriya spin-orbit coupling buttressed by altermagnetic order. After determining the easy axis, our calculation leads to an Os spin canting angle of about 3$^{\circ}$, accounting for the observed weak magnetism and the resulting discontinuity in the susceptibility. The altermagnetism spin-split bands (up to $\sim$100 meV) result in a chiral-split magnon spectrum in both acoustic and optical modes in the THz range, and lead to significant anomalous Hall conductivity upon hole doping.

Abinit is a widely used scientific software package implementing density functional theory and many related functionalities for excited states and response properties. This paper presents the novel features and capabilities, both technical and scientific, which have been implemented over the past 5 years. This evolution occurred in the context of evolving hardware platforms, high-throughput calculation campaigns, and the growing use of machine learning to predict properties based on databases of first principles results. We present new methodologies for ground states with constrained charge, spin or temperature; for density functional perturbation theory extensions to flexoelectricity and polarons; and for excited states in many-body frameworks including GW, dynamical mean field theory, and coupled cluster. Technical advances have extended abinit high-performance execution to graphical processing units and intensive parallelism. Second principles methods build effective models on top of first principles results to scale up in length and time scales. Finally, workflows have been developed in different community frameworks to automate \abinit calculations and enable users to simulate hundreds or thousands of materials in controlled and reproducible conditions.

The Landau-Lifshitz-Gilbert (LLG) equations are widely used to study spin dynamics in Mott insulators. However, because energy dissipation is typically introduced phenomenologically, its validity for describing nonequilibrium processes and long-time dynamics in real materials remains questionable. In this paper, we derive a generalized stochastic LLG equation from first principles for spin-orbital coupled Mott insulators, explicitly incorporating the coupling between electronic degrees of freedom and lattice vibrations. Our approach is based on the path-integral formalism formulated along the Keldysh contour, which naturally accounts for dissipation and thermal fluctuations through interactions with a phonon bath and emergent stochastic noise. We benchmark our theoretical framework by numerically integrating the equations of motion for a two-orbital spin chain coupled to Einstein phonons. The resulting energy relaxation mimics realistic cooling dynamics, exhibits nontrivial transient behavior during thermalization, and accurately reproduces thermodynamic properties upon equilibration. This demonstrates the robustness and reliability of our formalism for capturing both equilibrium and nonequilibrium aspects of dissipative spin dynamics in strongly correlated systems. We further show how electron-phonon coupling induces hybridization between electronic and phononic modes in the excitation spectrum. Finally, we demonstrate how the conventional LLG equations can be recovered as a limiting case of our microscopic theory, and discuss their applicability to real materials.