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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

Equivariant Graph Neural Networks (eGNNs) trained on density-functional theory (DFT) data can potentially perform electronic structure prediction at unprecedented scales, enabling investigation of the electronic properties of materials with extended defects, interfaces, or exhibiting disordered phases. However, as interactions between atomic orbitals typically extend over 10+ angstroms, the graph representations required for this task tend to be densely connected, and the memory requirements to perform training and inference on these large structures can exceed the limits of modern GPUs. Here we present a distributed eGNN implementation which leverages direct GPU communication and introduce a partitioning strategy of the input graph to reduce the number of embedding exchanges between GPUs. Our implementation shows strong scaling up to 128 GPUs, and weak scaling up to 512 GPUs with 87% parallel efficiency for structures with 3,000 to 190,000 atoms on the Alps supercomputer.

Spintronic devices typically employ heterostructures with ferromagnets which break time-reversal symmetry and have non-vanishing magnetization. With the growing class of materials that support spin-polarized carriers, current, and excitations,it is possible to envision emerging spintronic applications that are not limited to magnetoresistance. Here we focus on chiral perovskites with no net magnetization where the space-inversion and mirror symmetries are broken to induce chiral structure. The known importance of these perovskites is further expanded by the demonstration of the chiral-induced spin selectivity (CISS). However, the generation of the spin-polarized carriers across the interface with these chiral perovskites remains to be fully understood. Our first-principles studies for two-dimensional PbBr$_4$-based chiral perovskites provide their electronic structure and an orbital-based symmetry analysis, which allows us to establish an effective Hamiltonian to elucidate the underlying origin of their chirality. We also use this analysis for the Edelstein effect, responsible for electrical generation of the nonequilibrium spin polarization in many materials, which in chiral perovskites could be a mechanism contributing to CISS. Furthermore, by examining optical properties of chiral perovskites and the opportunity to use them to realize tunable altermagnets, another class of zero-magnetization spintronic materials, we put forth a versatile materials platform for unconventional spintronics.

A quantum state of matter that is forbidden to interact with photons and is therefore undetectable by spectroscopic means is called a dark state. This basic concept can be applied to condensed matter where it suggests that a whole band of quantum states could be undetectable across a full Brillouin zone. Here we report the discovery of such condensed matter dark states in palladium diselenide as a model system that has two pairs of sublattices in the primitive cell. By using angle-resolved photoemission spectroscopy, we find valence bands that are practically unobservable over the whole Brillouin zone at any photon energy, polarisation, and scattering plane. Our model shows that two pairs of sublattices located at half-translation positions and related by multiple glide-mirror symmetries make their relative quantum phases polarised into only four kinds, three of which become dark due to double destructive interference. This mechanism is generic to other systems with two pairs of sublattices, and we show how the phenomena observed in cuprates, lead-halide perovskites, and density wave systems can be resolved by the mechanism of dark states. Our results suggest that the sublattice degree of freedom, which has been overlooked so far, should be considered in the study of correlated phenomena and optoelectronic characteristics.

Generalised charge pumps are topological obstructions to trivialising loops in the space of symmetric gapped Hamiltonians. We show that given mild conditions on such pumps, the associated loop has high-symmetry points which must be in distinct symmetry-protected topological (SPT) phases. To further elucidate the connection between pumps and SPTs, we focus on closed paths, `pivot loops', defined by two Hamiltonians, where the first is unitarily evolved by the second `pivot' Hamiltonian. While such pivot loops have been studied as entanglers for SPTs, here we explore their connection to pumps. We construct families of pivot loops which pump charge for various symmetry groups, often leading to SPT phases -- including dipole SPTs. Intriguingly, we find examples where non-trivial pumps do not lead to genuine SPTs but still entangle representation-SPTs (RSPTs). We use the anomaly associated to the non-trivial pump to explain the a priori `unnecessary' criticality between these RSPTs. We also find that particularly nice pivot families form circles in Hamiltonian space, which we show is equivalent to the Hamiltonians satisfying the Dolan-Grady relation -- known from the study of integrable models. This additional structure allows us to derive more powerful constraints on the phase diagram. Natural examples of such circular loops arise from pivoting with the Onsager-integrable chiral clock models, containing the aforementioned RSPT example. In fact, we show that these Onsager pivots underlie general group cohomology-based pumps in one spatial dimension. Finally, we recast the above in the language of equivariant families of Hamiltonians and relate the invariants of the pump to the candidate SPTs. We also highlight how certain SPTs arise in cases where the equivariant family is labelled by spaces that are not manifolds.

Understanding the confinement mechanism in gauge theories and the universality of effective string-like descriptions of gauge flux tubes remains a fundamental challenge in modern physics. We probe string modes of motion with dynamical matter in a digital quantum simulation of a (2+1) dimensional gauge theory using a superconducting quantum processor with up to 144 qubits, stretching the hardware capabilities with quantum-circuit depths comprising up to 192 two-qubit layers. We realize the $Z_2$-Higgs model ($Z_2$HM) through an optimized embedding into a heavy-hex superconducting qubit architecture, directly mapping matter and gauge fields to vertex and link superconducting qubits, respectively. Using the structure of local gauge symmetries, we implement a comprehensive suite of error suppression, mitigation, and correction strategies to enable real-time observation and manipulation of electric strings connecting dynamical charges. Our results resolve a dynamical hierarchy of longitudinal oscillations and transverse bending at the end points of the string, which are precursors to hadronization and rotational spectra of mesons. We further explore multi-string processes, observing the fragmentation and recombination of strings. The experimental design supports 300,000 measurement shots per circuit, totaling 600,000 shots per time step, enabling high-fidelity statistics. We employ extensive tensor network simulations using the basis update and Galerkin method to predict large-scale real-time dynamics and validate our error-aware protocols. This work establishes a milestone for probing non-perturbative gauge dynamics via superconducting quantum simulation and elucidates the real-time behavior of confining strings.

We propose a new method based on the compression of arbitrary multi-point correlators into quantics tensor trains and studying its so called bond dimension as an indicator of time scale entanglement in the system. This quantics tensor train diagnostics approach allows us to i) identify quantum phase transitions (ground state crossings) and thermal crossovers, ii) distinguish between these two by allowing flexibility through iii) using arbitrary available correlators and iv) the possibility to deal with low-accuracy input data through varying the quantics tensor train cutoffs. We verify the validity of the approach by studying the phase diagram of the Hubbard dimer, the four-site Hubbard ring with and without next-nearest neighbor hopping and the single-impurity Anderson model. We find time scale entanglement to be a system-inherent property becoming maximal at phase transitions and crossover, which translates to peaks in the quantics tensor train bond dimension of correlator functions. Additionally, the connection of the method to other physical quantities, such as various entanglement measures and susceptibilities, is shown.

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.

Fully-compensated ferrimagnets exhibit zero net magnetic moment yet display non-relativistic global spin splitting, making them highly advantageous for constructing high-performance spintronic devices. The general strategy is to break the inversion symmetry of conventional antiferromagnets or the rotational/mirror symmetry of altermagnets to achieve fully-compensated ferrimagnets. Here, we propose to induce fully-compensated ferrimagnetism by engineering the spin ordering rather than modifying the lattice structure. Bilayer stacking engineering offers a convenient platform to verify our proposal and readily enables switching between two distinct electronic states by tuning the $\mathrm{N\acute{e}el}$ vector of one layer. By the first-principles calculations, a bilayer system is constructed with monolayer $\mathrm{Cr_2C_2S_6}$ as the elementary building block to corroborate our proposal. This strategy can also be extended to inducing altermagnetism via spin ordering engineering. Our work offers an alternative route to realize non-relativistic spin splitting in zero-net-magnetization magnets, paving the way for the advancement and construction of low-power spintronic device.

Ruthenium dioxide (RuO$_2$) has been in the focus of contemporary condensed matter research as a prototypical candidate material for altermagnetism. In the face of daunting evidence for bulk magnetic order despite promising theoretical predictions, it naturally suggests the focus on thin films where Coulomb interactions are dimensionally quenched and may yield a more strongly correlated environment prone to magnetic ordering. Here, we combine scanning tunneling microscopy (STM), density functional theory (DFT), and density matrix renormalization group (DMRG) methods to investigate atomically ordered ultrathin RuO$_2$(110) grown on Ru(0001). Contrary to predictions of magnetic order, we observe a nonmagnetic charge density wave (CDW) instability that is driven by Fermi surface nesting within the flat-band surface state, and stabilized by the incommensurate moir\'e stacking with the substrate. We further identify a nonmagnetic and metastable 2 $\times$ 2 surface reconstruction that breaks in-plane mirror symmetry and can be reversibly toggled via STM tip manipulation. Our spin-polarized STM measurements find no sign of any magnetic instability at the RuO$_2$ surface. As much as our findings refute proposals for either bulk or surface magnetism in RuO$_2$, we establish ultrathin RuO$_2$(110) as an intriguing platform for exploring Moir\'e-assisted electronic surface order.

The sublattice current is recently proposed as a feasible scheme to realize the odd-parity altermagnetism(ALM) by breaking the nonmagnetic time reversal symmetry. We first adopt the spin group method to analyze why the broken nonmagnetic time reversal symmetry(TRS) is the sufficient condition for the odd-parity ALM, and find that it is the symmetry $[\bar{C}_2||\bar{E}]$ that allows the appearance of the odd-parity ALM where $\bar{C}_2$ represents a $180^\circ$ rotation around an axis perpendicular to the spins combined with the spin-space inversion, and $\bar{E}$ is the inversion in real space. As a representative example with the presence of sublattice currents, the optical conductivity of Haldane-Hubbard model is further studied using the Kubo formula. It is shown that there is no spin-polarized electrical current in both the longitudinal and traverse directions because of TRS in the odd-parity ALM, and they display a significant peak in the vicinity of the single-particle direct energy gap which reflects that both the longitudinal and traverse optical conductivities are dominated by the quasiparticle excitations around the Dirac points. We also study the hole doping dependence of the anomalous Hall conductivity $\sigma_{\rm Hall}$ as well as the staggered magnetization $M$, and find that they display opposite behaviors, i.e., $\sigma_{\rm Hall}$ in the ALM CI state is diminished monotonically from its initial value $C=2$ at half-filling as the hole doping is increased, while in the ALMI state exhibits a non-monotonous behavior.

We propose a framework for the design and analysis of optimization algorithms in variational quantum Monte Carlo, drawing on geometric insights into the corresponding function space. The framework translates infinite-dimensional optimization dynamics into tractable parameter-space algorithms through a Galerkin projection onto the tangent space of the variational ansatz. This perspective unifies existing methods such as stochastic reconfiguration and Rayleigh-Gauss-Newton, provides connections to classic function-space algorithms, and motivates the derivation of novel algorithms with geometrically principled hyperparameter choices. We validate our framework with numerical experiments demonstrating its practical relevance through the accurate estimation of ground-state energies for several prototypical models in condensed matter physics modeled with neural network wavefunctions.

Altermagnets, with spin splitting and vanishing magnetization, have been attributed to many fascinating phenomena and potential applications. In particular, integrating ferroelectricity with altermagnetism to enable magnetoelectric coupling and electric control of spin has drawn significant attention. However, its experimental realization and precise spin manipulation remain elusive. Here, by focusing on molecular ferroelectrics, the first discovered ferroelectrics renowned for their highly controllable molecular polarizations and structural flexibility, we reveal that these obstacles can be removed by an emergent multiferroic altermagnets with tunable spin polarization in a large class of fabricated organic materials. Using a symmetry-based design and a tight-binding model, we uncover the underlying mechanism of such molecular ferroelectric altermagnets and demonstrate how noncollinear molecular polarization can switch the spin polarization on and off and even reverse its sign. From the first-principles calculations, we verify the feasibility of these materials in a series of well-established hybrid organic-inorganic perovskites and metal-organic frameworks. Our findings bridge molecular ferroelectrics and altermagnetic spintronics, highlighting an unexplored potential of multifunctional organic multiferroics.

The polaron problem is a very old problem in condensed matter physics that dates back to the thirties, but still remain largely unsolved today, especially when electron-electron interaction is taken into consideration. The presence of both electron-electron and electron-phonon interactions in the problem invalidates most existing numerical methods, either computationally too expensive or simply intractable. The continuous time quantum Monte Carlo (CTQMC) methods could tackle this problem, but are only effective in the imaginary-time axis. In this work we present a method based on tensor network and the path integral formalism to solve polaron impurity problems. As both the electron and phonon baths can be integrated out via the Feynman-Vernon influence functional in the path integral formalism, our method is free of bath discretization error. It can also flexibly work on the imaginary, Keldysh, and the L-shaped Kadanoff contour. In addition, our method can naturally resolve several long-existing challenges: (i) non-diagonal hybridization function; (ii) measuring multi-time correlations beyond the single particle Green's functions. We demonstrate the effectiveness and accuracy of our method with extensive numerical examples against analytic solutions, exact diagonalization and CTQMC. We also perform full-fledged real-time calculations that have never been done before to our knowledge, which could be a benchmarking baseline for future method developments.

The magnetic properties of two-dimensional altermagnets can be obtained from a square lattice Heisenberg model with antiferromagetic nearest neighbor interaction and two types of next-nearest neighbor interactions arranged in a checkerboard pattern. Using nonlinear spin-wave theory we calculate for this model the corrections to the renormalized magnon spectrum and the staggered magnetization to first order in the inverse spin quantum number $1/S$. We also show that to order $1/S^2$ the ground state energy is not sensitive to the component of the interaction which is responsible for altermagnetism. At the $\Gamma$-point the $1/S$-correction to the magnon dispersion vanishes so that quantum fluctuations do not induce a gap in the magnon spectrum of altermagnets, as expected by Goldstone's theorem. We extract the leading $1/S$-corrections to the spin-wave velocity and the effective mass characterizing the curvature of the magnon dispersion in altermagnets.

We revisit the infinite variance problem in fermionic Monte Carlo simulations, which is widely encountered in areas ranging from condensed matter to nuclear and high-energy physics. The different algorithms, which we broadly refer to as determinantal quantum Monte Carlo (DQMC), are applied in many situations and differ in details, but they share a foundation in field theory, and often involve fermion determinants whose symmetry properties make the algorithm sign-problem-free. We show that the infinite variance problem arises as the observables computed in DQMC tend to form heavy-tailed distributions. To remedy this issue retrospectively, we introduce a tail-aware error estimation method to correct the otherwise unreliable estimates of confidence intervals. Furthermore, we demonstrate how to perform DQMC calculations that eliminate the infinite variance problem for a broad class of observables. Our approach is an exact bridge link method, which involves a simple and efficient modification to the standard DQMC algorithm. The method introduces no systematic bias and is straightforward to implement with minimal computational overhead. Our results establish a practical and robust solution to the infinite variance problem, with broad implications for improving the reliability of a variety of fundamental fermion simulations.

Terahertz two-dimensional coherent spectroscopy (THz-2DCS) is transforming our ability to probe, visualize, and control quantum materials far from equilibrium. This emerging technique brings multi-dimensional resolution to the ultrafast dynamics of nonequilibrium phases of matter, enabling new capabilities demanding precise coherent control and measurement of many-body dynamics and multi-order correlations. By mapping complex excitations across time and frequency dimensions, THz-2DCS delivers coherence tomography of driven quantum matter, thus revealing hidden excitation pathways, measuring higher order nonlinear response functions, disentangling various quantum pathways, capturing collective modes on ultrafast timescales and at terahertz frequencies. These experimental features frequently remain obscured in traditional single particle measurements, ultrafast spectroscopy techniques, and equilibrium-based probes. This Review traces the early development of THz-2DCS and showcases significant recent progress in leveraging this technique to probe and manipulate quantum material properties, including nonequilibrium superconductivity, nonlinear magnonics, dynamical topological phases, and the detection of novel excitations and exotic collective modes with potential technological impact. Looking forward, we identify critical opportunities in advancing THz-2DCS instrumentation and experimental strategies that are shaping future applications in THz optoelectronics, quantum information processing, and sensing.

Spontaneous magnon decays in an $S=\frac{1}{2}$ Heisenberg antiferromagnet on the equilateral triangular prism are investigated in terms of modified magnon Green's functions. In one dimension, the so-called infrared divergence prevents us from calculating any -- whether static or dynamic -- structure factor within the conventional spin-wave theory even at zero temperature. The well-known modified spin-wave theory initiated by Takahashi completely fails to treat anharmonicities to cause transverse-to-longitudinal coupling which are quite characteristic of noncollinear antiferromagnets. We propose imposing a double-constraint condition on spin waves to solve all these difficulties and get a full view of the nonlinear spin-wave dynamics in one-dimensional frustrated noncollinear antiferromagnets. We reveal a novel instability of the single-particle spectrum in the absence of any well-defined magnetically ordered ground state.

We study the ground state phase diagram of the extended Hubbard model in a half-filled 5/7 skewed ladder, which is topologically equivalent to a \{112\} grain boundary in graphene and related systems. Using the mean-field method, we identify various electronic and magnetic phases in the U-V plane, by calculating the site charge and spin densities. The electronic phases include partially charge-ordered metal or insulator, and fully charge-ordered insulator. The different magnetic phases of the model are non-magnet, spin density wave, spin split compensated ferrimagnet or partial antiferromagnet. Analysis of the electronic band structure reveals that the partially charge-ordered compensated ferrimagnetic phase exhibits spin polarisation, which can be quite interesting for spintronics applications. We also compute the polarisation as a function of $U$ using the Berry phase formalism and show that the system exhibits multiferroicity with coexisting compensated ferrimagnetic spin order alongside electronic polarisations.

Recent experiments observing femtosecond laser pulse (fsLP) exciting magnons in two-dimensional (2D) antiferromagnetic (AF) semiconductors -- such as CrSBr, NiPS$_3$, and MnPS$_3$, or their van der Waals heterostructures -- suggest exciton-mediation of such an effect. However, its microscopic details remain obscure as resonant coupling of magnons, living in the sub-meV energy range, to excitons, living in \mbox{$\sim 1$ eV} range, can hardly be operative. Here, we develop a quantum transport theory of this effect, in which time-dependent nonequilibrium Green's function (TDNEF) for electrons driven by fsLP are coupled self-consistently to the Landau-Lifshitz-Gilbert (LLG) equation describing classical dynamics of localized magnetic moments (LMMs) within 2D AF semiconductor. This theory explains how fsLP, of central frequency above the semiconductor gap, generates a photocurrent that subsequently exerts spin-transfer torque (STT) onto LMMs as a {\em nonequilibrium spintronic mechanism}. The collective motion of LMMs analyzed by windowed Fast Fourier transform (FFT) decodes frequencies of excited magnons, as well as their lifetime governed by {\em nonlocal} damping with the LLG equation due to, explicitly included via TDNEGF, electronic bath. The TDNEGF part of the loop is also used to include excitons via mean-field treatment, utilizing off-diagonal elements of the density matrix, of Coulomb interaction binding conduction-band electrons and valence-band holes. Finally, our theory predicts how excited magnons will {\em pump} time-dependent charge currents into the attached electrodes, or locally within AF semiconductor that will then emit electromagnetic radiation. The windowed FFT of these signals contains imprints of excited magnons, as well as their interaction with excitons, which could be exploited as a novel probe in future experiments.

Recent experiments observing femtosecond laser pulse (fsLP) exciting magnons within two-dimensional (2D) antiferromagnetic (AF) semiconductors -- such as CrSBr, NiPS$_3$, and MnPS$_3$, or their van der Waals heterostructures -- suggest exciton-mediation of such an effect. However, its microscopic details remain obscure as resonant coupling of magnons, living in the sub-meV energy range, to excitons, living in \mbox{$\sim 1$ eV} range, can hardly be operative. Here, we develop a quantum transport theory of this effect, in which time-dependent nonequilibrium Green's function (TDNEF) for electrons driven by classical vector potential of fsLP are coupled to the Landau-Lifshitz-Gilbert (LLG) equation describing classical dynamics of localized magnetic moments (LMMs) within 2D AF semiconductor. Our TDNEGF+LLG theory explains how fsLP, with central frequency above the semiconductor gap, generates a photocurrent that subsequently exerts spin-transfer torque (STT) onto LMMs as a genuine nonequilibrium spintronic effect. The collective motion of LMMs analyzed by windowed Fast Fourier transform (FFT) reveals frequencies of excited magnons, as well as their lifetime governed by nonlocal damping due to the bath of electrons. In addition, the TDNEGF part of our TDNEGF+LLG self-consistent loop computes a time-dependent density matrix whose off-diagonal elements are utilized to describe, at the mean-field level, inter-orbital Coulomb interaction binding electrons and holes into excitons. Our TDNEGF+LLG theory predicts how excited magnons {\em pump} charge current into the attached electrodes, or locally within AF semiconductor responsible for microwave emission. The windowed FFT of the former signal contains imprints of excited magnons, as well as their interaction with excitons, which could be exploited as a novel probe in future experiments.

We develop a tensor network framework based on the quantic tensor train (QTT) format to efficiently solve the Gross-Pitaevskii equation (GPE), which governs Bose-Einstein condensates under mean-field theory. By adapting time-dependent variational principle (TDVP) and gradient descent methods, we accurately handle the GPE's nonlinearities within the QTT structure. Our approach enables high-resolution simulations with drastically reduced computational cost. We benchmark ground states and dynamics of BECs--including vortex lattice formation and breathing modes--demonstrating superior performance over conventional grid-based methods and stable long-time evolution due to saturating bond dimensions. This establishes QTT as a powerful tool for nonlinear quantum simulations.

We develop a tensor network framework based on the quantic tensor train (QTT) format to efficiently solve the Gross-Pitaevskii equation (GPE), which governs Bose-Einstein condensates under mean-field theory. By adapting time-dependent variational principle (TDVP) and gradient descent methods, we accurately handle the GPE's nonlinearities within the QTT structure. Our approach enables high-resolution simulations with drastically reduced computational cost. We benchmark ground states and dynamics of BECs--including vortex lattice formation and breathing modes--demonstrating superior performance over conventional grid-based methods and stable long-time evolution due to saturating bond dimensions. This establishes QTT as a powerful tool for nonlinear quantum simulations.

Simple binary oxide glasses can exhibit a compression behavior distinct from that of their crystalline counterparts. In this study, we employed high-pressure X-ray absorption spectroscopy coupled to the diamond anvil cell to investigate in detail local structural changes around Ge in glassy $GeO_{2}$ up to 158 GPa. We conducted four independent runs, both with and without pressure-transmitting media. Up to 30 GPa, we observed no significant influence of the pressure medium on the pressure dependence of the $Ge-O$ bond length ($$). Between 10 and 30 GPa, the evolution of $$ shows substantial variability across our experiments and previous works. The measured values lie close to those reported for crystalline polymorphs, including the rutile- and $CaCl_{2}$-type phase of $GeO_{2}$. This finding suggests that the amorphous structure possesses considerable flexibility to transition among different atomic configurations. From 30 GPa to 158 GPa, our results for both $$ and the non-bonded cation-cation distance $$ demonstrate that edge-sharing octahedra remain the main structural motives in glassy $GeO_{2}$. Up to 100 GPa, compaction proceeds primarily via distortions of octahedral $O-Ge-O$ bond angles accompanied by octahedral bond shortening. Above 100 GPa, octahedral distortion becomes the prevailing mechanism. Compared to its crystalline analogues ($\alpha-PbO_{2}$ and pyrite-like phase), glassy $GeO_{2}$ exhibits a slightly less efficient compaction mechanism, likely due to kinetic constraints that inhibit reconstructive lattice rearrangements.

Incorporating zero-net-magnetization magnets that exhibit spin-splitting into spintronics delivers key advantages: faster switching dynamics, greater immunity to destabilizing fields, lower power consumption, and markedly improved overall efficiency. The collinear magnets with net-zero magnetization and spin-splitting mainly include altermagnet and fully-compensated ferrimagnet, which provide possibility to achieve hidden spin polarization (HSP) with net-zero spin polarization in total but non-zero local spin polarization. In addition to proposal of hidden altermagnetism, we hereby introduce this concept of hidden fully-compensated ferrimagnetism, where the total spin polarization is zero, but either of the two inversion-partner sectors possesses fully-compensated ferrimagnetism with non-zero local spin polarization in the real space. By the first-principle calculations, we predict that $PT$-bilayer $\mathrm{CrMoC_2S_6}$ is a possible hidden fully-compensated ferrimagnet, showing fully-compensated ferrimagnetic HSP, which can be separated and observed by an out-of-plane external electric field. Our works provide a class of hidden spin-polarized materials that facilitates the advancement of spintronics.

The fluctuation-dissipation theorem (FDT) is a fundamental result in statistical mechanics. It stipulates that, if perturbed out of equilibrium, a system responds at a rate proportional to a thermal-equilibrium property. Applications range from particle diffusion to electrical-circuit noise. To prove the FDT, one must prove that common thermal states obey a symmetry property, the Kubo-Martin-Schwinger (KMS) relation. Energy eigenstates of certain quantum many-body systems were recently proved to obey the KMS relation. The proof relies on the eigenstate thermalization hypothesis (ETH), which explains how such systems thermalize internally. This KMS relation contains a finite-size correction that scales as the inverse system size. Non-Abelian symmetries conflict with the ETH, so a non-Abelian ETH was proposed recently. Using it, we derive a KMS relation for SU(2)-symmetric quantum many-body systems' energy eigenstates. The finite-size correction scales as usual under certain circumstances but can be polynomially larger in others, we argue. We support the ordinary-scaling result numerically, simulating a Heisenberg chain of 16-24 qubits. The numerics, limited by computational capacity, indirectly support the larger correction. This work helps extend into nonequilibrium physics the effort, recently of interest across quantum physics, to identify how non-Abelian symmetries may alter conventional thermodynamics.

We develop a general framework based on the functional derivative to extract nonlinear dynamical response functions from the temporal evolution of physical quantities, without explicitly computing multipoint correlation functions. We validate our approach by calculating the second- and third-order optical responses in the Rice-Mele model and further apply it to a many-body interacting system using a tensor network method. This framework is broadly applicable to any method that can compute real-time dynamics, offering a powerful and versatile tool for investigating nonlinear responses in dynamical systems.

Quantum dots (QDs) exhibit unique properties arising from their reduced size and quantum confinement effects, including exceptionally bright and tunable photoluminescence. Among these, CuInS$_{2}$ QDs have gained significant attention owing to their remarkable broadband emission, making them highly desirable for various optoelectronic applications requiring efficient luminescent nanomaterials. However, maximizing radiative recombination in CuInS$_{2}$ QDs often necessitates minimizing intragap trap states. A common approach involves the introduction of Zn during the synthesis, which typically promotes the formation of a ZnS shell that passivates the QD surface. Despite its importance, the characterization and quantification of Zn incorporation using conventional techniques, such as optical spectroscopy or electron microscopy, remains challenging. In this study, we utilized X-ray absorption spectroscopy (XAS), in both X-ray absorption near-edge structure (XANES) and extended X-ray absorption fine structure (EXAFS) spectral ranges, to investigate Zn incorporation into CuInS$_{2}$ QDs with element-specific precision. This approach allowed us to detect the formation of a ZnS surface shell and to resolve the spatial distribution of Zn atoms within the QD lattice, distinguishing between Zn as a substituent, or as an interstitial defect. Additionally, we explored the optical and dynamical properties of CuInS$_{2}$ QDs using time-resolved optical spectroscopies, particularly in the presence of electron and hole acceptors. These results provide deeper insights into the role and effectiveness of the Zn-induced passivating layer, paving the way for optimizing QD performance in photoluminescence applications.

Materials with tailored quantum properties can be engineered from atomic scale assembly techniques, but existing methods often lack the agility and accuracy to precisely and intelligently control the manufacturing process. Here we demonstrate a fully autonomous approach for fabricating atomic-level defects using electron beams in scanning transmission electron microscopy (STEM) that combines advanced machine learning and automated beam control. As a proof of concept, we achieved controlled fabrication of MoS-nanowire (MoS-NW) edge structures by iterative and targeted exposure of $MoS_2$ monolayer to a focused electron beam to selectively eject sulfur atoms, utilizing high-angle annular dark-field (HAADF) imaging for feedback-controlled monitoring structural evolution of defects. A machine learning framework combining a random forest model and convolutional neural networks (CNN) was developed to decode the HAADF image and accurately identify atomic positions and species. This atomic-level information was then integrated into an autonomous decision-making platform, which applied predefined fabrication strategies to instruct beam control about atomic sites to be ejected. The selected sites were subsequently exposed to localized electron beam using an FPGA-controlled scan routine with precise control over beam positioning and duration. While the MoS-NW edge structures produced exhibit promising mechanical and electronic properties, the proposed autonomous fabrication framework is material-agnostic and can be extended to other 2D materials for the creation of diverse defect structures and heterostructures beyond $MoS_2$.

Significant challenges exist in efficient data analysis of most advanced experimental and observational techniques because the collected signals often include unwanted contributions--such as background and signal distortions--that can obscure the physically relevant information of interest. To address this, we have developed a self-supervised machine-learning approach for source separation using a dual implicit neural representation framework that jointly trains two neural networks: one for approximating distortions of the physical signal of interest and the other for learning the effective background contribution. Our method learns directly from the raw data by minimizing a reconstruction-based loss function without requiring labeled data or pre-defined dictionaries. We demonstrate the effectiveness of our framework by considering a challenging case study involving large-scale simulated as well as experimental momentum-energy-dependent inelastic neutron scattering data in a four-dimensional parameter space, characterized by heterogeneous background contributions and unknown distortions to the target signal. The method is found to successfully separate physically meaningful signals from a complex or structured background even when the signal characteristics vary across all four dimensions of the parameter space. An analytical approach that informs the choice of the regularization parameter is presented. Our method offers a versatile framework for addressing source separation problems across diverse domains, ranging from superimposed signals in astronomical measurements to structural features in biomedical image reconstructions.

New discoveries in chemistry and materials science, with increasingly expanding volume of requisite knowledge and experimental workload, provide unique opportunities for machine learning (ML) to take critical roles in accelerating research efficiency. Here, we demonstrate (1) the use of large language models (LLMs) for automated literature reviews, and (2) the training of an ML model to predict chemical knowledge (thermodynamic parameters). Our LLM-based literature review tool (LMExt) successfully extracted chemical information and beyond into a machine-readable structure, including stability constants for metal cation-ligand interactions, thermodynamic properties, and other broader data types (medical research papers, and financial reports), effectively overcoming the challenges inherent in each domain. Using the autonomous acquisition of thermodynamic data, an ML model was trained using the CatBoost algorithm for accurately predicting thermodynamic parameters (e.g., enthalpy of formation) of minerals. This work highlights the transformative potential of integrated ML approaches to reshape chemistry and materials science research.

Accurate grain orientation mapping is essential for understanding and optimizing the performance of polycrystalline materials, particularly in energy-related applications. Lithium nickel oxide (LiNiO$_{2}$) is a promising cathode material for next-generation lithium-ion batteries, and its electrochemical behaviour is closely linked to microstructural features such as grain size and crystallographic orientations. Traditional orientation mapping methods--such as manual indexing, template matching (TM), or Hough transform-based techniques--are often slow and noise-sensitive when handling complex or overlapping patterns, creating a bottleneck in large-scale microstructural analysis. This work presents a machine learning-based approach for predicting Euler angles directly from scanning transmission electron microscopy (STEM) diffraction patterns (DPs). This enables the automated generation of high-resolution crystal orientation maps, facilitating the analysis of internal microstructures at the nanoscale. Three deep learning architectures--convolutional neural networks (CNNs), Dense Convolutional Networks (DenseNets), and Shifted Windows (Swin) Transformers--are evaluated, using an experimentally acquired dataset labelled via a commercial TM algorithm. While the CNN model serves as a baseline, both DenseNets and Swin Transformers demonstrate superior performance, with the Swin Transformer achieving the highest evaluation scores and the most consistent microstructural predictions. The resulting crystal maps exhibit clear grain boundary delineation and coherent intra-grain orientation distributions, underscoring the potential of attention-based architectures for analyzing diffraction-based image data. These findings highlight the promise of combining advanced machine learning models with STEM data for robust, high-throughput microstructural characterization.

We calculate the dynamical current and kinetic-energy correlation functions for the first Mott lobe of the one-dimensional Bose-Hubbard model. We employ the strong-coupling expansion up to sixth order in $x=t/U$, and the dynamical density-matrix renormalization group method on rings with 64 sites. The correlation functions are finite above the single-particle gap with a square-root onset, as is also found from field theory close to the Mott transition. The correlation functions display a featureless superposition of the primary and tertiary Hubbard bands. We find very good agreement between all methods in the interaction/frequency regimes where they are applicable.

Studies relying on hydrodynamic theory and Kardar-Parisi-Zhang (KPZ) scaling have found that in the one-dimensional Hubbard model spin and charge transport are for all temperatures T > 0 anomalous superdiffusive at zero magnetic field, h = 0, and zero chemical potential, {\mu} = 0, respectively. However, this contradicts recent exact results that at very low temperature charge transport rather is normal diffusive. In this Letter we identify the mechanisms that control the different types of temperature dependence of the h = 0 spin and {\mu} = 0 charge transport and find that the latter is normal diffusive for all finite temperatures T > 0, in contrast to the hydrodynamic theory and KPZ scaling predictions.

Neural quantum states are a promising framework for simulating many-body quantum dynamics, as they can represent states with volume-law entanglement. As time evolves, the neural network parameters are typically optimized at discrete time steps to approximate the wave function at each point in time. Given the differentiability of the wave function stemming from the Schr\"odinger equation, here we impose a time-continuous and differentiable parameterization of the neural network by expressing its parameters as linear combinations of temporal basis functions with trainable, time-independent coefficients. We test this ansatz, referred to as the smooth neural quantum state ($s$-NQS) with a loss function defined over an extended time interval, under a sudden quench of a non-integrable many-body quantum spin chain. We demonstrate accurate time evolution using simply a restricted Boltzmann machine as the instantaneous neural network architecture. Furthermore, we demonstrate that the parameterization is efficient in the number of parameters and the smooth neural quantum state allows us to initialize and evaluate the wave function at times not included in the training set, both within and beyond the training interval.

Neural quantum states are a promising framework for simulating many-body quantum dynamics, as they can represent states with volume-law entanglement. As time evolves, the neural network parameters are typically optimized at discrete time steps to approximate the wave function at each point in time. Given the differentiability of the wave function stemming from the Schr\"odinger equation, here we impose a time-continuous and differentiable parameterization of the neural network by expressing its parameters as linear combinations of temporal basis functions with trainable, time-independent coefficients. We test this ansatz, referred to as the smooth neural quantum state ($s$-NQS) with a loss function defined over an extended time interval, under a sudden quench of a non-integrable many-body quantum spin chain. We demonstrate accurate time evolution using simply a restricted Boltzmann machine as the instantaneous neural network architecture. Furthermore, we demonstrate that the parameterization is efficient in the number of parameters and the smooth neural quantum state allows us to initialize and evaluate the wave function at times not included in the training set, both within and beyond the training interval.

The possibility that spatial geometry may emerge from the entanglement structure of a quantum many-body system is a subject of fundamental interest. Here, we propose and numerically test a candidate distance metric in 1D, d_E, defined purely from quantum mutual information (I) via the relation d_E = K_0 / sqrt(I). Using large-scale density-matrix renormalization group (DMRG) simulations, we compute this emergent distance for the ground state of the 1D spin-1/2 XXZ chain, a canonical model system. Our simulations show that in the quantum critical phase at anisotropy $\Delta = 1.0$, the mutual information exhibits a power-law decay consistent with the emergence of a valid metric space. In stark contrast, within the gapped, antiferromagnetic phase ($\Delta = 2.0$), where mutual information decays exponentially, the emergent distance grows exponentially, a behavior inconsistent with the triangle inequality. These results provide numerical evidence that this information-theoretic definition can yield a well-behaved geometry in critical systems, offering a quantitative tool for probing quantum phases and motivating further analytical investigation into the foundations of emergent space.

We present an instance of utility-grade quantum computation by calculating the ground-state energy of a 125-site flat Kagome lattice under the antiferromagnetic Heisenberg model (KAFH), using IBM's Falcon and Hummingbird quantum processors. For spin-1/2 KAFH, our best per-site ground-state energy estimate reaches -0.417J, and after applying open-boundary corrections, it closely approaches the established thermodynamic value of -0.438J. To achieve this, we propose a hybrid approach that splits the variational quantum eigensolver (VQE) into local (classical) and global (quantum) components for efficient hardware utilization. We further introduce a Hamiltonian engineering strategy that increases coupling on defect triangles to mimic loop-flip dynamics, allowing us to simplify the ansatz while retaining physical accuracy. Using a single-repetition, hardware-efficient ansatz, we entangle up to 103 qubits with high fidelity to determine the Hamiltonian's lowest eigenvalue. This work demonstrates the scalability of VQE for frustrated 2D systems and lays the foundation for future studies using deeper ansatz circuits and larger lattices.