The recent discovery of superconductivity with a transition temperature $T_c$ over 40 K in La$_3$Ni$_2$O$_7$ and (La,Pr)$_{3}$Ni$_2$O$_7$ thin films at ambient pressure marks an important step in the field of nickelate superconductors. Here, we perform a renormalized mean-field theory study of the superconductivity in $\mathrm{La_3Ni_2O_7}$ thin films, using a bilayer two-orbital $t-J$ model. Our result reveals an $s_\pm$-wave pairing symmetry driven by the strong interlayer superexchange coupling of $d_{z^2}$ orbital, resembling the pressurized bulk case. Also, we roughly reproduce the experimentally reported nodeless shape of the superconducting gap at the $\beta$ pocket and the superconducting $T_c$. To gain insight into the orbital feature of superconductivity, we explore the projection of different pairing bonds on Fermi surface. We find that the nodeless gap at $\beta$ pocket is related to the interlayer pairing within both $d_{z^2}$ and $d_{x^2-y^2}$ orbitals, in which the latter is triggered by the former through hybridization, and both hold the same sign.

Employing combination of first-principles calculations, low-energy model construction, and variational Monte Carlo solution of the ab-initio derived Hubbard model, we study the effect of hydrogenation in the electronic structure and superconducting properties of infinite-layer nickelate, NdNiO$_2$. We find that the introduction of hydrogen at the apical oxygen vacancy position strongly influences the Wannier function corresponding to the effective interstitial orbital at the Ni site bound to the hydrogen. This results in the near disappearance of the electron pocket at the $k_z$ = $\pi$ Fermi surface, keeping that of $k_z$ = 0 unchanged, compared to the dehydrogenated case. The two-band model thus remains valid even in the presence of H. The calculated superconducting order parameters both in absence and presence of H, show a two-hump superconductivity arising the two overlapping domes, one arising from $d_{x^{2}-y^{2}}$ and another arising from interstitial orbital degree of freedom. Hydrogenation strengthens the latter, marginally affecting the former.

High-temperature superconductivity in pressurized and strained bilayer nickelates (La,Pr)$_3$Ni$_2$O$_7$ has emerged as a new frontier. One of the key unresolved issues concerns the fermiology that underlies superconductivity. On both theoretical and experimental sides, no general consensus has been reached, and conflicting results exist regarding whether the relevant Fermi surface involves a hole $\gamma$ pocket, or not. Here, we address this issue by unveiling a Janus-faced role of the $\gamma$ pocket in spin-fluctuation-mediated superconductivity. We show that this pocket simultaneously induces dominant pair-breaking and pair-forming fluctuations for the leading $s_\pm$-wave pairing. Consequently, an optimal superconducting transition temperature $T_\mathrm{c}$ is achieved when the $\gamma$ pocket surfaces at the Fermi level, placing the system near a Lifshitz transition. This suggests that superconductivity can emerge provided the maximum energy level of the $\gamma$ pocket lies sufficiently close to the Fermi level, either from below or above. Our finding not only reconciles two opposing viewpoints on the fermiology, but also naturally explains recent experiments on (La,Pr)$_3$Ni$_2$O$_7$ thin films, including the superconductivity under compressive strain, two conflicting measurements on the Fermi surface, and the dome shape of $T_\mathrm{c}$ as a function of hole doping.

We theoretically study quasi-two-dimensional nickel compounds, where the nickel ions assume Ni$^{2+}$ $d^8$ valence state and feature a low-spin $S=0$ ground state quasidegenerate with $S=1$ ionic excitations. Such a level structure is supported by square-planar coordination of nickel ions or a suitable substitution of apical oxygens. We construct the corresponding singlet-triplet exchange model and explore its phase diagram and excitation spectrum. By hole doping, we further introduce mobile Ni$^{3+}$ $d^7$ ionic configurations with effective spins $S=1/2$, and analyze their interactions with the $d^8$ singlet-triplet background. The interplay with the triplet excitations in the $d^8$ sector is found to have a deep impact on the propagation of the doped hole-like charge carriers and is identified as a powerful source of Cooper pairing among them.

The mysterious quantum-critical Planckian bad metal phase with perfect T-linear resistivity persisting beyond the quasi-particle limit and universal T-linear scattering rate has been observed in various high-Tc cuprate superconductors. Here, we develop a realistic theoretical approach to this phase in an analytically solvable large-N multi-channel Kondo lattice model, derived from a heavy-fermion formulated conventionaL t-J model, known for qualitatively describing cuprates. This phase is originated from critical charge Kondo fluctuations where disordered local bosonic charge fluctuations couple to spinon and heavy conduction-electron Fermi surfaces near a charge-Kondo-breakdown local quantum critical point associated with pseudogap-to-Fermi liquid transition. Our results show excellent agreement with experiments and offer broad implications for other unconventional superconductors.

The ghost Gutzwiller ansatz (gGut) embedding technique was shown to achieve comparable accuracy to the gold standard dynamical mean-field theory method in simulating real material properties, yet at a much lower computational cost. Despite that, gGut is limited by the algorithmic bottleneck of computing the density matrix of the underlying effective embedding model, a quantity which must be converged within a self-consistent embedding loop. We develop a hybrid quantum-classical gGut technique which computes the ground state properties of embedding Hamiltonians with the help of a quantum computer, using the sample-based quantum-selected configuration interaction (QSCI) algorithm. We study the applicability of SCI-based methods to the evaluation of the density of states for single-band Anderson impurity models within gGut and find that such ground states of interest become sufficiently sparse in the CI basis as the number of ghost orbitals is increased. Further, we investigate the performance of QSCI using local unitary cluster Jastrow (LUCJ) variational quantum states in combination with a circuit cutting technique, prepared on IQM's quantum hardware for system sizes of up to 11 ghost orbitals, equivalent to 24 qubits. We report converged gGut calculations which correctly capture the metal-to-insulator phase transition in the Fermi-Hubbard model on the Bethe lattice by using quantum samples to build an SCI basis with as little as $1\%$ of the total CI basis states.

A supercurrent is well recognized as being of prime importance within mean-field theory, but remains largely unexplored in strongly correlated electron systems (SCES) and the quantum critical region. To clarify the impact of the supercurrent on magnetism and superconductivity near an antiferromagnetic quantum critical point, we study the two-dimensional Hubbard model based on a fluctuation exchange approximation for a current-carrying superconducting state. We show a supercurrent-induced antiferromagnetism and emergence of spin-triplet Cooper pairs. The former results from Bogoliubov Fermi surfaces, suppression in the superconducting gap, and strong correlation effects beyond the mean-field theory. Our results suggest that the supercurrent can bring out rich phenomena of superconductivity in SCES.

Inspired by the recent discovery of high-temperature superconductivity in bilayer nickelates, we investigate the role of magnetism emerging from a hypothetical insulating $d^8$ parent state. We demonstrate that due to the interplay of superexchange and Hund's coupling, the system can be in a high-spin, low-spin or spin-gapped state. The low-spin state has singlets across the bilayer in the $d_{z^2}$ orbital, with charge carriers in the $d_{x^2-y^2}$ orbital. Thus, at low energy scales, it behaves as an effective one band system when hole doped. By contrast, the high-spin state is a more robust, spin-1 antiferromagnet. Using Hartree-Fock methods, we find that for fixed interaction strength and doping, high-spin magnetism remains more robust than the low-spin counterpart. Whether this implies that the high spin state provides a stronger pairing glue, or more strongly competes with superconductivity remains an open question. Our analysis therefore underscores the importance of identifying the spin state for understanding superconductivity in nickelates.

The superconducting gap is a characteristic feature of high-T$_c$ superconductors and provides crucial information on the pairing mechanism underlying high-temperature superconductivity. Here, we employ high-resolution resonant inelastic X-ray scattering (RIXS) at the Cu $L_3$-edge to investigate the superconducting gap in the overdoped cuprate $\mathrm{Bi}_2\mathrm{Sr}_2\mathrm{Ca}_2\mathrm{Cu}_3\mathrm{O}_{10+\delta}$ ($T_\mathrm{c}$ = 107 K). By analyzing antisymmetrized, temperature-dependent RIXS spectra over a range of in-plane momentum transfers, we observe a clear suppression of low-energy spectral weight below T$_c$, indicative of superconducting gap formation. This suppression is most pronounced at small momentum transfers ($|\boldsymbol{q}_\parallel| \leq 0.18$ r.l.u.) and corresponds to a gap size of approximately 2$\Delta_0 \sim$ 130 meV. Comparison with theoretical calculations of the momentum-dependent charge susceptibility supports a d-wave symmetry of the superconducting gap, while an isotropic s-wave gap fails to reproduce key experimental features. These findings establish RIXS as a powerful, bulk-sensitive probe of superconducting gap symmetry and highlight its utility for studying materials beyond the reach of surface-sensitive techniques such as ARPES and STM.

We propose a theory for bilayer nickelate materials, where a large tetragonal field - intrinsic or induced by epitaxial strain - lifts the orbital degeneracy and localizes the $3z^2-r^2$ orbital states. These states host local spins $S=1/2$ bound into singlets by strong interlayer coupling, and their dynamics is described by weakly dispersive singlet-triplet excitations ("triplons"). The charge carriers occupy the wide bands of $x^2-y^2$ symmetry, and their Cooper pairing is mediated by the high-energy triplon excitations. As the $x^2-y^2$ band filling increases, i.e. moving further away from the Ni$^{3+}$ valence state, the indirect Ruderman-Kittel-Kasuya-Yosida interactions between local spins induce spin-density-wave order via triplon condensation. Implications of the model for compressively strained La$_3$Ni$_2$O$_7$ films and electron doped oxychloride Sr$_3$Ni$_2$O$_5$Cl$_2$ are discussed.

We compute high-resolution angle-resolved photoemission spectroscopy of the Hubbard model using the unbiased determinant quantum Monte Carlo algorithm, revealing an asymmetry between electron and hole doping. Electron doping exhibits more coherent quasiparticles and stronger antiferromagnetic correlations compared to hole doping. At low doping, a nodal-antinodal dichotomy on the Fermi surface is observed, similar to cuprate experiments. The dichotomy reflects the momentum dependence of the Mott gap, as manifested in both the spectral function and the self-energy. For hole doping, we observe a transition towards the pseudogap, without signature of pocket formation. The simulated nuclear magnetic resonance pseudogap temperatures do not necessarily agree with the temperature determined by spectroscopy. These findings collectively suggest the pseudogap is a smooth crossover driven by strong correlations.

We investigate the ground-state properties of the two-dimensional Hubbard model, based on the off-diagonal wave function variational Monte Carlo method. We use an optimized wave function that is improved from an initial one-body wave function by multiplying by multiple correlation operators that are given by the form of exp($-S$)-type where $S$ is a suitable operator. We examine the inhomogeneous ground state at near 1/8 doping in the strongly correlated region where the on-site Coulomb interaction is larger than the bandwidth. The $d$-wave superconductivity with a spatially oscillating gap function can coexist with the charge ordering in the ground state without magnetic ordering, and superconducting condensation energy increases by this coexistence. We also show the $d$-wave pair correlation function as a function of lattice sites. The correlation function indicates that the long-range superconducting order indeed exists, and also that the strong electron correlation suppresses the pair correlation function.

Antiferromagnetism and superconductivity are often viewed as competing orders in correlated electron systems. Here, we demonstrate that kinetic frustration in hole motion facilitates their coexistence within the square-lattice repulsive Hubbard model. Combining exact analytical solutions on tailored geometries with large-scale numerical simulations, we reveal a robust pairing mechanism: holes on opposite sublattices behave as if they carry opposite effective charges due to spin singlet formation from kinetic frustration. This emergent property suppresses phase separation and fosters a coherent $d$-wave superconducting channel embedded within a long-range antiferromagnetic background. Our findings establish a minimal yet broadly applicable framework for stabilizing strong-coupling superconductivity in doped Mott insulators.

We explore the finite-density phase diagram of the single-flavour Gross-Neveu-Wilson (GNW) model using matrix product state (MPS) simulations. At zero temperature and along the symmetry line of the phase diagram, we find a sequence of inhomogeneous ground states that arise through a real-space version of the mechanism of Hilbert-space fragmentation. For weak interactions, doping the symmetry-protected topological (SPT) phase of the GNW model leads to localized charges or holes at periodic arrangements of immobile topological defects separating the fragmented subchains: a topological crystal. Increasing the interactions, we observe a transition into a parity-broken phase with a pseudoscalar condensate displaying a modulated periodic pattern. This soliton lattice is a sequence of topological charges corresponding to anti-kinks, which also bind the doped fermions at their respective centers. Out of this symmetry line, we show that quasi-spiral profiles appear with a characteristic wavevector set by the density $k = 2{\pi}{\rho}$, providing non-perturbative evidence for chiral spirals beyond the large-N limit. These results demonstrate that various exotic inhomogeneous phases can arise in lattice field theories, and motivate the use of quantum simulators to confirm such QCD-inspired phenomena in future experiments.

The multipoint numerical renormalization group (mpNRG) is a powerful impurity solver that provides accurate spectral data useful for computing local, dynamic correlation functions in imaginary or real frequencies non-perturbatively across a wide range of interactions and temperatures. It gives access to a local, non-perturbative four-point vertex in imaginary and real frequencies, which can be used as input for subsequent computations such as diagrammatic extensions of dynamical mean--field theory. However, computing and manipulating the real-frequency four-point vertex on large, dense grids quickly becomes numerically challenging when the density and/or the extent of the frequency grid is increased. In this paper, we compute four-point vertices in a strongly compressed quantics tensor train format using quantics tensor cross interpolation, starting from discrete partial spectral functions obtained from mpNRG. This enables evaluations of the vertex on frequency grids with resolutions far beyond the reach of previous implementations. We benchmark this approach on the four-point vertex of the single-impurity Anderson model across a wide range of physical parameters, both in its full form and its asymptotic decomposition. For imaginary frequencies, the full vertex can be represented to an accuracy on the order of $2\cdot 10^{-3}$ with maximum bond dimensions not exceeding 120. The more complex full real-frequency vertex requires maximum bond dimensions not exceeding 170 for an accuracy of $\lesssim 2\%$. Our work marks another step toward tensor-train-based diagrammatic calculations for correlated electronic lattice models starting from a local, non-perturbative mpNRG vertex.

The Sn/Si(111)-({\sqrt}3{\times}{\sqrt}3)R30{\deg} surface, a 2D Mott insulator, has long been predicted and then found experimetally to metallize and even turn superconducting upon boron doping. In order to clarify the structural, spectroscopic and theoretical details of that evolution, here we present ARPES data supplementing morphology and scanning tunneling measurements. These combined experimental results are compared with predictions from a variety of electronic structure approaches, mostly density functional DFT+U, but not neglecting Mott-Hubbard models, both ordered and disordered. These theoretical pictures address different spectroscopic aspects, including the 2D Fermi surface, the Hubbard bands, etc. While no single picture account for all observations at once,the emergent hypothesis compatible with all data is that metallization arises from sub-subsurface boron doping, additional to the main standard subsurface boron geometry, that would leave the surface insulating. These results advance the indispensable frame for the further understanding of this fascinating system.

Non-Hermitian (NH) quantum systems have emerged as a powerful framework for describing open quantum systems, non-equilibrium dynamics, and engineered quantum optical materials. However, solving the ground-state properties of NH systems is challenging due to the exponential scaling of the Hilbert space, and exotic phenomena such as the emergence of exceptional points. Another challenge arises from the limitations of traditional methods like exact diagonalization (ED). For the past decade, neural networks (NN) have shown promise in approximating many-body wavefunctions, yet their application to NH systems remains largely unexplored. In this paper, we explore different NN architectures to investigate the ground-state properties of a parity-time-symmetric, one-dimensional NH, transverse field Ising model with a complex spectrum by employing a recurrent neural network (RNN), a restricted Boltzmann machine (RBM), and a multilayer perceptron (MLP). We construct the NN-based many-body wavefunctions and validate our approach by recovering the ground-state properties of the model for small system sizes, finding excellent agreement with ED. Furthermore, for larger system sizes, we demonstrate that the RNN outperforms both the RBM and MLP. However, we show that the accuracy of the RBM and MLP can be significantly improved through transfer learning, allowing them to perform comparably to the RNN for larger system sizes. These results highlight the potential of neural network-based approaches--particularly for accurately capturing the low-energy physics of NH quantum systems.

Nonlinear dynamics govern a wide array of natural phenomena and are essential for understanding nonequilibrium behaviors in condensed matter systems. In magnetically ordered materials, magnons - the quanta of spin waves - exhibit intrinsic nonlinearities that are of great interest in fundamental research and practical applications. Despite progress in the nonlinear control of magnon modes in antiferromagnetic materials, the transition from perturbative to non-perturbative regimes of magnon coherences has remained elusive. Here, we explore the nonlinear dynamics of a magnon mode in an antiferromagnet using two-dimensional terahertz spectroscopy with waveguide-enhanced terahertz fields. By driving the magnon mode far from equilibrium, we demonstrate the emergence of high-order magnon coherences and delineate a distinct transition into non-perturbative magnon nonlinearities. This behavior originates from the intrinsic anharmonicity of the magnetic potential and marks a regime dominated by magnon self-interactions at large spin deflection angles. These findings provide fundamental mechanistic insights that might be exploited for ultrafast switching and other advanced magnonic applications.

We propose a Cooper pair analogue of Dzyaloshinskii-Moriya's anisotropic coupling for $d$-vectors in an inhomogeneous unconventional superconductor, which can arise when itinerant electrons couple to a local exchange field of frustrated spin moments. Using perturbative Green's function methods, we examine an $s$-$d$ model on a geometrically frustrated lattice with noncollinear spin texture that breaks time-reversal and parity symmetries, leading to spin-triplet pairing correlations. The effective Josephson tunneling in the presence of such an exchange field leads to an anisotropic Josephson coupling between $d$-vectors, which prefers a spatially inhomogeneous pairing order. Furthermore, we discuss anomalous vortices arising from the spatially inhomogeneous $d$-vector in the absence of a magnetic field for nonunitary triplet pairing order. We also propose a Josephson diode effect in which the diode efficiency is proportional to the spin chirality of the local exchange field. Our results can help understand the experiments in superconducting proximitized Mn$_3$Ge and in 4Hb-TaS$_2$.

We analyze the quantum geometric contribution to the intrinsic second-order nonlinear Hall effect (NLHE) for a general multiband Hamiltonian. The nonlinear conductivity, obtained in Green's function formalism, is decomposed into its quantum geometric constituents using a projector-based approach. In addition to the previously identified Berry curvature and interband quantum metric dipoles, we obtain a third term of quantum geometric origin, given by the momentum derivative of the $intraband$ quantum metric. This contribution, which we term the intraband quantum metric dipole, provides substantial corrections to the NLHE in topological magnets and becomes the dominant geometric term in topological antiferromagnets with gapped Dirac cones. Considering generalized 2D and 3D Weyl/Dirac Hamiltonians, describing a large class of topological band crossings with sizable quantum geometry, we derive analytical expressions of the NLHE, thereby revealing the individual contributions of the three quantum geometric terms. Combined with an exhaustive symmetry classification of all magnetic space groups, this analysis leads to the identification of several candidate materials expected to exhibit large intrinsic NLHE, including the antiferromagnets $\text{Yb}_3\text{Pt}_4$, $\text{CuMnAs}$, and $\text{CoNb}_3\text{S}_6$, as well as the nodal-plane material $\text{MnNb}_3\text{S}_6$. Finally, our projector-based approach yields a compact expression for the NLHE in terms of momentum derivatives of the Bloch Hamiltonian matrix alone, enabling efficient numerical evaluation of each contribution in the aforementioned materials.

We demonstrate that the local nonfreeness, an unbiased measure of correlation between electrons at a single lattice site, can be computed as the mutual information between local natural spin orbitals. Using this concept, we prove that local electron correlations in the Hubbard model are fully classical: in the natural basis the local reduced density matrix is separable and no quantum correlations beyond entanglement are present. Finally, we compare different theoretical descriptions of magnetic and nonmagnetic states, showing that local classical correlations are drastically influenced by nonlocal processes. Our results confirm the relation between local classical correlations within an open system and the nonlocal correlations with its quantum environment.

SU(N) Hubbard models exhibit a rich variety of phases, which may be realized through quantum simulation with ultracold atomic gases in optical lattices. In this work we study the Mott insulating phases of the SU(3) Hubbard model at 1/3-filling using infinite projected entangled-pair states, optimized with both imaginary time evolution and variational optimization. In the limit of strong interactions we reproduce the antiferromagnetic 3-sublattice ordered state previously identified in the SU(3) Heisenberg model. At intermediate interaction strength we find antiferromagnetic states exhibiting zigzag patterns of different lengths, in agreement with previous Hartree-Fock and constrained-path auxiliary-field quantum Monte Carlo calculations. We study the color order parameter and energy anisotropy, which are discontinuous across the phase transitions. Finally, we analyze the different energy contributions in two competing phases, identifying low-energy bonds at the corners of the zigzag that help stabilize the zigzag states.

We study an extension of the 2D Fermi--Hubbard model, which was recently introduced in [Das et al., Phys. Rev. Lett. 132, 263402 (2024)] and shown to describe altermagnetism that can be studied in cold atom systems. Using an updated Hartree--Fock method that can detect instabilities towards phase separation, we show that the model is in a mixed phase in large parts of the parameter regime at half-filling. We argue that the occurrence of a mixed phase is an indication of exotic physics which, in this model, occurs in parameter regimes accessible in cold atom experiments.

The development of machine-learning models for atomic-scale simulations has benefited tremendously from the large databases of materials and molecular properties computed in the past two decades using electronic-structure calculations. More recently, these databases have made it possible to train universal models that aim at making accurate predictions for arbitrary atomic geometries and compositions. The construction of many of these databases was however in itself aimed at materials discovery, and therefore targeted primarily to sample stable, or at least plausible, structures and to make the most accurate predictions for each compound - e.g. adjusting the calculation details to the material at hand. Here we introduce a dataset designed specifically to train machine learning models that can provide reasonable predictions for arbitrary structures, and that therefore follows a different philosophy. Starting from relatively small sets of stable structures, the dataset is built to contain massive atomic diversity (MAD) by aggressively distorting these configurations, with near-complete disregard for the stability of the resulting configurations. The electronic structure details, on the other hand, are chosen to maximize consistency rather than to obtain the most accurate prediction for a given structure, or to minimize computational effort. The MAD dataset we present here, despite containing fewer than 100k structures, has already been shown to enable training universal interatomic potentials that are competitive with models trained on traditional datasets with two to three orders of magnitude more structures. We describe in detail the philosophy and details of the construction of the MAD dataset. We also introduce a low-dimensional structural latent space that allows us to compare it with other popular datasets and that can be used as a general-purpose materials cartography tool.

A quantitative first-principles description of complex substitutional materials like alloys is challenging due to the vast number of configurations and the high computational cost of solving the quantum-mechanical problem. Therefore, materials properties must be modeled. The Cluster Expansion (CE) method is widely used for this purpose, but it struggles with properties that exhibit non-linear dependencies on composition, often failing even in a qualitative description. By looking at CE through the lens of machine learning, we resolve this severe problem and introduce a non-linear CE approach, yielding extremely accurate and computationally efficient results as demonstrated by distinct examples.

We study thermodynamic consequences of the hidden Bose-Einstein singularities, which have been predicted to cause a pseudogap phase, based on quantum field theory of ordered phases. Starting from the Luttinger-Ward functional for the grand thermodynamic potential, we derive expressions of the Helmholtz free energy, internal energy, entropy, and heat capacity in a form suitable for numerical studies. They are applied to the weakly attractive Hubbard model in three dimensions to calculate the thermodynamic potentials numerically and continuously for the normal, pseudogap, and superconducting phases on the same footing. It is shown that the entry into the pseudogap phase is detectable as singularities of thermodynamic potentials, especially a discontinuity in the heat capacity.

Large-scale foundation models, including neural network interatomic potentials (NIPs) in computational materials science, have demonstrated significant potential. However, despite their success in accelerating atomistic simulations, NIPs face challenges in directly predicting electronic properties and often require coupling to higher-scale models or extensive simulations for macroscopic properties. Machine learning (ML) offers alternatives for structure-to-property mapping but faces trade-offs: feature-based methods often lack generalizability, while deep neural networks require significant data and computational power. To address these trade-offs, we introduce HackNIP, a two-stage pipeline that leverages pretrained NIPs. This method first extracts fixed-length feature vectors (embeddings) from NIP foundation models and then uses these embeddings to train shallow ML models for downstream structure-to-property predictions. This study investigates whether such a hybridization approach, by ``hacking" the NIP, can outperform end-to-end deep neural networks, determines the dataset size at which this transfer learning approach surpasses direct fine-tuning of the NIP, and identifies which NIP embedding depths yield the most informative features. HackNIP is benchmarked on Matbench, evaluated for data efficiency, and tested on diverse tasks including \textit{ab initio}, experimental, and molecular properties. We also analyze how embedding depth impacts performance. This work demonstrates a hybridization strategy to overcome ML trade-offs in materials science, aiming to democratize high-performance predictive modeling.

We calculate the spin correlations that result from the motion of a single dopant in the hard-core Fermi-Hubbard model, as the geometry evolves from a square to a triangular lattice. In particular, we consider the square lattice with an additional hopping along one diagonal, whose strength is continuously varied. We use a high-temperature expansion which expresses the partition function as a sum over closed paths taken by the dopant. We sample thousands of diagrams in the space of closed paths using the quantum Monte Carlo approach of Raghavan and Elser [1,2], which is free of finite-size effects and allows us to simulate temperatures as low as $T \sim 0.3|t|$, even in cases where there is a sign problem. For the case of a hole dopant, we find a crossover from kinetic ferromagnetism to kinetic antiferromagnetism as the geometry is tuned from square to triangular, which can be observed in current quantum gas microscopes.

In this work we investigate the ground state properties of a candidate quantum spin liquid using a superconducting Noisy Intermediate-Scale Quantum (NISQ) device. Specifically, we study the antiferromagnetic Heisenberg model on a Kagome lattice, a geometrically frustrated structure that gives rise to a highly degenerate energy spectrum. To successfully simulate this system, we employ a qubit reduction strategy leveraging the Contextual Subspace methodology, significantly reducing the problem size prior to execution on the quantum device. We improve the quality of these subspaces by using the wavefunctions obtained from low bond dimension Density Matrix Renormalization Group (DMRG) calculations to bias the subspace stabilizers through a symplectic approximate symmetry generator extraction algorithm. Reducing the Hamiltonian size allows us to implement tiled circuit ensembles and deploy the Variational Quantum Eigensolver (VQE) to estimate the ground state energy. We adopt a hybrid quantum error mitigation strategy combining Readout Error Mitigation (REM), Symmetry Verification (SV) and Zero Noise Extrapolation (ZNE). This approach yields high-accuracy energy estimates, achieving error rates on the order of 0.01% and thus demonstrating the potential of near-term quantum devices for probing frustrated quantum materials.

In recent years, the variational Monte Carlo (VMC) calculations of projected entangled pair states (PEPS) has emerged as a competitive method for computing the ground states of many-body quantum systems. This method is particularly important for fermion systems where sign problems are abundant. We derive and explain the algorithms for the VMC calculations of fermion PEPS in the swap gates formulation. It is the purpose of this paper to be as concise and precise as possible. Diagramatic tensor notation is used whenever possible to serve this purpose. As a separate result, we prove the detailed balance of sequential sampling of tensor networks.

Raman spectroscopy is a powerful experimental technique for characterizing molecules and materials that is used in many laboratories. First-principles theoretical calculations of Raman spectra are important because they elucidate the microscopic effects underlying Raman activity in these systems. These calculations are often performed using the canonical harmonic approximation which cannot capture certain thermal changes in the Raman response. Anharmonic vibrational effects were recently found to play crucial roles in several materials, which motivates theoretical treatments of the Raman effect beyond harmonic phonons. While Raman spectroscopy from molecular dynamics (MD-Raman) is a well-established approach that includes anharmonic vibrations and further relevant thermal effects, MD-Raman computations were long considered to be computationally too expensive for practical materials computations. In this perspective article, we highlight that recent advances in the context of machine learning have now dramatically accelerated the involved computational tasks without sacrificing accuracy or predictive power. These recent developments highlight the increasing importance of MD-Raman and related methods as versatile tools for theoretical prediction and characterization of molecules and materials.

Mn-rich disordered rocksalt (DRX) cathode materials exhibit a phase transformation from a disordered to a partially disordered spinel-like structure ($\delta$-phase) during electrochemical cycling. In this computational study, we used charge-informed molecular dynamics with a fine-tuned CHGNet foundation potential to investigate the phase transformation in Li$_{x}$Mn$_{0.8}$Ti$_{0.1}$O$_{1.9}$F$_{0.1}$. Our results indicate that transition metal migration occurs and reorders to form the spinel-like ordering in an FCC anion framework. The transformed structure contains a higher concentration of non-transition metal (0-TM) face-sharing channels, which are known to improve Li transport kinetics. Analysis of the Mn valence distribution suggests that the appearance of tetrahedral Mn$^{2+}$ is a consequence of spinel-like ordering, rather than the trigger for cation migration as previously believed. Calculated equilibrium intercalation voltage profiles demonstrate that the $\delta$-phase, unlike the ordered spinel, exhibits solid-solution signatures during the 0-TM to Li$_{\text{tet}}$ conversion reaction. A higher Li capacity is obtained than in the DRX phase. This study provides atomic insights into solid-state phase transformation and its relation to experimental electrochemistry, highlighting the potential of charge-informed machine learning interatomic potentials for understanding complex oxide materials.

The optical spectra of two-dimensional (2D) semiconductors are dominated by tightly bound excitons and trions. In the low doping limit, trions are often described as three-body quasiparticles consisting of two electrons and one hole or vice versa. However, trions are more rigorously understood as quasiparticles arising from the interaction between an exciton and excitation of the Fermi sea - referred to as exciton Fermi polaron. Here we employ pump-probe spectroscopy to directly observe the formation of exciton Fermi polarons in a model system composed of a WSe2 monolayer adjacent to twisted bilayer graphene (tBLG). Following the pump-injection of Rydberg excitons in WSe2, a time-delayed probe pulse tracks the development of Rydberg exciton Fermi polarons as interactions with localized carriers in the tBLG moire superlattice evolve. Both the exciton Fermi polaron relaxation rate and binding energy are found to increase with electron or hole density. Our findings provide insight into the optical response of fundamental excitations in 2D Van der Waals systems and reveal how many-body interactions give rise to emergent quasiparticles.

We present the first proof of principle that normalizing flows can accurately learn the Boltzmann distribution of the fermionic Hubbard model - a key framework for describing the electronic structure of graphene and related materials. State-of-the-art methods like Hybrid Monte Carlo often suffer from ergodicity issues near the time-continuum limit, leading to biased estimates. Leveraging symmetry-aware architectures as well as independent and identically distributed sampling, our approach resolves these issues and achieves significant speed-ups over traditional methods.

Advancing complementary metal-oxide-semiconductor (CMOS) technology into the sub-1-nm angstr\"om-scale technology nodes is expected to involve alternative semiconductor channel materials, as silicon transistors encounter severe performance degradation at physical gate lengths below 10 nm. Two-dimensional (2D) semiconductors have emerged as strong candidates for overcoming short-channel effects due to their atomically thin bodies, which inherently suppress electrostatic leakage and improve gate control in aggressively scaled field-effect transistors (FETs). Among the growing library of 2D materials, the MoSi$_2$N$_4$ family -- a synthetic septuple-layered materials -- has attracted increasing attention for its remarkable ambient stability, suitable bandgaps, and favorable carrier transport characteristics, making it a promising platform for next-generation transistors. While experimental realization of sub-10-nm 2D FETs remains technologically demanding, computational device simulation using first-principles density functional theory combined with nonequilibrium Green's function transport simulations provide a powerful and cost-effective route for exploring the performance limits and optimal design of ultrascaled FET. This review consolidates the current progress in the computational design of MoSi$_2$N$_4$ family FETs. We review the physical properties of MoSi$_2$N$_4$ that makes them compelling candidates for transistor applications, as well as the simulated device performance and optimization strategy of MoSi$_2$N$_4$ family FETs. Finally, we identify key challenges and research gaps, and outline future directions that could accelerate the practical deployment of MoSi$_2$N$_4$ family FET in the angstr\"om-scale CMOS era.

We integrate concepts from topological band theory and strong-correlation physics with the principle of orbital symmetry conservation, which characterizes reactions as symmetry allowed or forbidden based on whether molecular orbitals cross along the reaction coordinate. Using a $4\pi$ electrocyclization as an example, we show how Green's functions generalize the concepts of symmetry allowed and symmetry forbidden reactions even in the presence of strong multi-reference correlations. We demonstrate how symmetry forbidden reactions are characterized by strong multi-reference correlations, resulting in crossings of Green's function zeros rather than poles as MO-theory would predict. Taking into account Green's functions zeros, a topological invariant is introduced that captures symmetry protected crossings of poles or zeros. We discuss the effects of symmetry breaking and outline generalizations of our approach to reactions without any conserved spatial symmetries along the reaction path. Our work lays the groundwork for systematic application of modern topological methods to chemical reactions and can be generalized to reactions involving different spin states or to excited states.

The Czochralski (Cz) method is a widely used process for growing high-quality single crystals, critical for applications in semiconductors, optics, and advanced materials. Achieving optimal growth conditions requires precise control of process and furnace design parameters. Still, data scarcity -- especially for new materials -- limits the application of machine learning (ML) in predictive modeling and optimization. This study proposes a transfer learning approach to overcome this limitation by adapting ML models trained on a higher data volume of one source material (Si) to a lower data volume of another target material (Ge and GaAs). The materials were deliberately selected to assess the robustness of the transfer learning approach in handling varying data similarity, with Cz-Ge being similar to Cz-Si, and GaAs grown via the liquid encapsulated Czochralski method (LEC), which differs from Cz-Si. We explore various transfer learning strategies, including Warm Start, Merged Training, and Hyperparameters Transfer, and evaluate multiple ML architectures across two different materials. Our results demonstrate that transfer learning significantly enhances predictive accuracy with minimal data, providing a practical framework for optimizing Cz growth parameters across diverse materials.

Weak first-order pseudocriticality with approximate scale invariance has been observed in a variety of settings, including the intriguing case of deconfined criticality in 2+1 dimensions. Recently, this has been interpreted as extremely slow flows ("walking behavior") for real-valued couplings in proximity to a bona fide critical point with complex-valued couplings, described by a complex conformal field theory (CFT). Here we study an SU($N$) generalization of the the Heisenberg antiferromagnet, which is a familiar model for deconfined criticality in 2+1 dimensions. We show that in 1+1 dimensions the model is located near a complex CFT, whose proximity can be tuned as a function of $N$. We employ state-of-the-art quantum Monte Carlo simulations for continuous $N$ along with an improved loop estimator for the R\'{e}nyi entanglement entropy based on a nonequilibrium work protocol. These techniques allow us to track the central charge of this model in detail as a function of $N$, where we observe excellent agreement with CFT predictions. Notably, this includes the region $N>2$, where the CFT moves into the complex plane and pseudocritical drifts enable us to recover the real part of the complex central charge with remarkable accuracy. Since the present model with $N=3$ is also equivalent to the spin-1 biquadratic model, our work sheds new light on the dimerized phase of the spin-1 chain, demonstrating that it is pseudocritical and proximate to a complex CFT.

Ultrafast pump-probe imaging reveals that the efficiency of optical excitation of coherent spins waves in epitaxial iron garnet films can be effectively controlled by an external electric field at room temperature. Although a femtosecond laser pulse alone does not excite any pronounced coherent spin oscillations, an electrical gating with the field of 0.5 MV/m dramatically changes the outcome in a laser-induced launching of spin waves. The effect, demonstrated under room-temperature conditions, is estimated to be orders of magnitude larger than in magnetic van der Waals semiconductors observed at 10 K. This electrical gating of laser-induced spin dynamics enriches opto-magnonics with a new tool and thus opens up a new avenue in fundamental and applied magnonics research.

Machine learning has significantly advanced the understanding and application of structural materials, with an increasing emphasis on integrating existing data and quantifying uncertainties in predictive modeling. This study presents a comprehensive methodology utilizing a mixed density network (MDN) model, trained on extensive experimental data from literature. This approach uniquely predicts the distribution of dislocation density, inferred as a latent variable, and the resulting stress distribution at the grain level. The incorporation of statistical parameters of those predicted distributions into a dislocation-mediated plasticity model allows for accurate stress-strain predictions with explicit uncertainty quantification. This strategy not only improves the accuracy and reliability of mechanical property predictions but also plays a vital role in optimizing alloy design, thereby facilitating the development of new materials in a rapidly evolving industry.

Electronic and optoelectronic applications of two-dimensional (2D) semiconductors demand precise control over material quality, including thickness, composition, doping, and defect density. Conventional benchmarking methods (e.g., charge transport, confocal mapping, electron or scanning probe microscopy) are slow, perturb sample quality, or involve trade-offs between speed, resolution, and scan area. To accelerate assessment of 2D semiconductors, we demonstrate a broadband, wide-field hyperspectral optical microscope for 2D materials (2D-HOM) that rapidly captures a spatial-spectral data cube within seconds. The data cube includes x-y spatial coordinate (a 300 * 300 $\mu\mathrm{m}^2$ field, with ~ 1 $\mu\mathrm{m}$ resolution) and a selectable wavelength range between 1100 to 200 nm at each pixel. Using synthesized films and heterostructures of transition metal dichalcogenides ($\mathrm{MoS}_{2}$, $\mathrm{WS}_{2}$, $\mathrm{V}_{x}\mathrm{W}_{1-x}\mathrm{S}_{2}$, and $\mathrm{WSe}_{2}$), we show that this cost-effective technique detects spectral fingerprints of material identity, doping, grain boundaries, and alloy composition, and enables advanced analysis, including unsupervised machine learning for spatial segmentation.

Polycrystalline materials have numerous applications due to their unique properties, which are often determined by the grain boundaries. Hence, quantitative characterization of grain as well as interface orientation is essential to optimize these materials, particularly energy materials. Using scanning transmission electron microscopy, matter can be analysed in an extremely fine grid of scan points via electron diffraction patterns at each scan point. By matching the diffraction patterns to a simulated database, the crystal orientation of the material as well as the orientation of the grain boundaries at each scan point can be determined. This pattern matching approach is highly time intensive. Artificial intelligence promises to be a very powerful tool for pattern recognition. In this work, we train convolutional neural networks (CNNs) on dynamically simulated diffraction patterns of LiNiO2, an important cathode-active material for Lithium-ion batteries, to predict the orientation of grains in terms of three Euler angles for the complete fundamental orientation region. Results demonstrate that these networks outperform the conventional pattern matching algorithm with increased accuracy and efficiency. The increased accuracy of the CNN models can be attributed to the fact that these models are trained by data incorporating dynamical effects. This work is the first attempt to apply deep learning for analysis of electron diffraction data and enlightens the great potential of ML to accelerate the analysis of electron microscopy data, toward high-throughput characterization technique.

Herein, we present a novel analysis framework for grain size profile of ferroelectric hafnia to tackle critical shortcomings inherent in the current microstructural analysis. We vastly enhanced visibility of grains with ion beam treatment and performed accurate grain segmentation using deep neural network (DNN). By leveraging our new method, we discovered unexpected discrepancies that contradict previous results, such as deposition temperature (Tdep) and post-metallization annealing (PMA) dependence of grain size statistics, prompting us to reassess earlier interpretations. Combining microstructural analysis with electrical tests, we found that grain size reduction had both positive and negative outcomes: it caused significant diminishing of die-to-die variation (~68 % decrease in standard deviation) in coercive field (Ec), while triggering an upsurge in leakage current. These uncovered results signify robustness of our method in characterization of ferroelectric hafnia for in-depth examination of both device variability and reliability.

DFT is vital for materials discovery, and at the base of extensive molecular and materials databases, chemical reaction predictions, and machine learning potentials. The widespread use of DFT in chemistry and materials science aims for "chemical accuracy," but this is limited by the unknown exchange and correlation (XC) functional. A meta-GGA, the restored regularized strongly constrained and appropriately normed, r r2SCAN XC functional, fulfils 17 exact constraints of the XC energy. r2SCAN still appears inadequate at predicting material properties of strongly correlated compounds. Inaccuracies of r2SCAN arise from functional and density-driven errors, linked to the self-interaction error. We introduce a new method, r2SCANY@r2SCANX, for simulating transition metal oxides accurately. r2SCANY@r2SCANX utilizes different fractions of exact exchange: X to set the electronic density, and Y to set the energy density functional approximation. r2SCANY@r2SCANX addresses functional-driven and density-driven inaccuracies. Using just one or two universal parameters, r2SCANY@r2SCANX enhances the r2SCAN predictions of the properties of 18 correlated oxides, outperforming the highly parameterized DFT+U method. The O2 overbinding in r2SCAN (~0.3 eV/O2) reduces to just ~0.03 eV/O$_2$ with any X in r2SCAN10@r2SCANX. Uncertainties for oxide oxidation energies and magnetic moments are reduced by r2SCAN10@r2SCAN50, minimizing r2SCAN density-driven errors. The computationally efficient r2SCAN10@r2SCAN is nearly as accurate as the hybrid r2SCAN10 for oxidation energies. Thus, accurate energy differences can be achieved by rate-limiting self-consistent iterations and geometry optimizations with the efficient r2SCAN. Subsequently, expensive hybrid functionals can be applied in a fast-to-execute single post-self-consistent calculation, as in r2SCAN10@r2SCAN, which is 10 to 300 times faster than r2SCAN10.

Free energy calculations are at the heart of physics-based analyses of biochemical processes. They allow us to quantify molecular recognition mechanisms, which determine a wide range of biological phenomena from how cells send and receive signals to how pharmaceutical compounds can be used to treat diseases. Quantitative and predictive free energy calculations require computational models that accurately capture both the varied and intricate electronic interactions between molecules as well as the entropic contributions from motions of these molecules and their aqueous environment. However, accurate quantum-mechanical energies and forces can only be obtained for small atomistic models, not for large biomacromolecules. Here, we demonstrate how to consistently link accurate quantum-mechanical data obtained for substructures to the overall potential energy of biomolecular complexes by machine learning in an integrated algorithm. We do so using a two-fold quantum embedding strategy where the innermost quantum cores are treated at a very high level of accuracy. We demonstrate the viability of this approach for the molecular recognition of a ruthenium-based anticancer drug by its protein target, applying traditional quantum chemical methods. As such methods scale unfavorable with system size, we analyze requirements for quantum computers to provide highly accurate energies that impact the resulting free energies. Once the requirements are met, our computational pipeline FreeQuantum is able to make efficient use of the quantum computed energies, thereby enabling quantum computing enhanced modeling of biochemical processes. This approach combines the exponential speedups of quantum computers for simulating interacting electrons with modern classical simulation techniques that incorporate machine learning to model large molecules.

We show that two Yukawa-SYK models with a weak tunneling contact can have an exotic hybrid superconducting thermofield-double-like state that is holographically dual to a traversable wormhole connecting two black holes with charged scalar hair. The hybrid superconducting thermo-field-double/wormhole state is distinguishable by anomalous scaling of revival oscillations in the fermionic Green's function, but also in a unique Andreev-revival in the anomalous Green's function. The existence of this TFD/wormhole state surprisingly shows that the some quantum critical effects can survive the phase transition to superconductivity. This Andreev-revival is in principle an accessible signature of the transition to the TFD/wormhole phase detectable in the ac-Josephson current.

The Mott-Hubbard metal-insulator transition is a paradigmatic phenomenon where Coulomb interactions between electrons drive a metal-insulator phase transition. It has been extensively studied within the Hubbard model, where a quantum critical transition occurs at a finite temperature second-order critical point. This work investigates the Mott metal-insulator transition in a modified periodic Anderson model that may be viewed as a three-orbital lattice model including an interacting, localized orbital coupled to a delocalized conduction orbital via a second conduction orbital. This model could also be viewed as a bilayer model involving a conventional periodic Anderson model layer coupled to a metallic layer. Within the dynamical mean field theory, this model possesses a strictly zero temperature quantum critical point separating a Fermi liquid and a Mott insulating phase. By employing a simplified version of the dynamical mean field theory, namely, the two-site or linearized dynamical mean field theory, we provide an analytical estimate of the critical parameter strengths at which the transition occurs at zero temperature. We also provide an analytical estimate of the single-site von Neumann entanglement entropy. This measure can be used as a robust identifier for the phase transition. We extend these calculations to their cluster version to incorporate short-range, non-local spatial correlations and discuss their effects on the transition observed in this model.

In strongly correlated electron systems, superconductivity and charge density waves often coexist in close proximity, suggesting a deeper relationship between these competing phases. Recent research indicates that these orders can intertwine, with the superconducting order parameter coupling to modulations in the electronic density. To elucidate this interplay, we study a two-dimensional XY model with a periodic modulation of the coupling strength in one spatial direction. Using a combination of tensor network methods and Monte Carlo simulations, we reveal a non-monotonic, dome-like dependence of $T_c$ on the modulation wavelength, with the peak $T_c$ shifting to longer wavelengths as the modulation strength grows. The origin of this phenomenon is traced back to an effective pinning of vortices in the valleys of the modulation, confirmed by a comparison to modulated $q$-state clock models. These findings shed new light on the phase behavior of intertwined superconducting and charge-ordered states, offering a deeper understanding of their complex interactions.

Motivated by recent experiments, we investigate the spin-1/2 XXZ model on the triangular lattice with strong Ising anisotropy, combining large-scale numerical simulations and analytical methods to uncover unconventional spin dynamics. First, we compute the dynamical spin structure factor using density matrix renormalization group (DMRG) simulations and find excellent agreement with inelastic neutron scattering data on the layered compound $\text{K}_2\text{Co}(\text{SeO}_3)_2$. The low-energy spectrum reveals a roton-like minimum at the $M$ point, absent in linear spin-wave theory, accompanied by peak intensity and a broad continuum above it. Near the $\Gamma$ point, we observe an approximately linear dispersion with vanishing spectral weight. Second, we compare two analytical frameworks that reproduce the observed features. The first is a hard-core boson approach, which includes: (i) an effective staggered boson model (ESBM) at zero magnetic field, (ii) perturbation theory applied to the one-third magnetization plateau, and (iii) a self-consistent mean-field Schwinger boson theory (SBT). The second framework is based on a variational supersolid quantum dimer model (QDM) ansatz, combined with a single-mode approximation. The SBT captures the broad continuum, the $M$-point minimum, and linear dispersion at $\Gamma$, whereas the QDM reproduces the roton minimum and linear dispersion at finite momentum near $\Gamma$. Remarkably, both the QDM wavefunction and the DMRG ground state exhibit nearly identical structure factors with pronounced transverse photon-like excitations. Together, our comprehensive theoretical and numerical analysis elucidates the microscopic origin of supersolid excitations in the XXZ triangular lattice model and their proximity to a spin liquid phase observed experimentally.

We study the dynamical structure factor of the frustrated spin-$3/2$ $J_1$-$J_2$ Heisenberg chains, with particular focus on the partially dimerized phase that emerges between two Kosterlitz-Thouless transitions. Using a valence bond solid ansatz corroborated by density matrix renormalization group simulations, we investigate the nature of magnon and spinon excitations through the single-mode approximation. We show that the magnon develops an incommensurate dispersion at $J_2 \approx 0.32J_1$, while the spinons, viewed as domain walls between degenerate valence bond solid states, become incommensurate at $J_2 \approx 0.4J_1$ beyond the Lifshitz point ($J_2 \approx 0.388J_1$). The dynamical structure factor exhibits rich spectral features shaped by the interplay between these excitations, with magnons appearing as resonances embedded in the spinon continuum. The spinon gap shows a nonmonotonic behavior, reaching a peak near the center of the partially dimerized phase and closing at the boundaries, suggesting the appearance of a floating phase as a result of the condensation of incommensurate spinons. Comparative analysis with the spin-$5/2$ case confirms the universality of these phenomena across half-integer higher-spin systems. Our results provide detailed insight into how fractionalization and incommensurate condensation govern the spectral properties of frustrated spin chains, offering a unified picture across different spin magnitudes.

In quantum many-body physics, one aims to understand emergent phenomena and effects of strong interactions, ideally by developing a simple theoretical picture. Recently, progress in quantum simulators has enabled the measurement of site resolved snapshots of Fermi-Hubbard systems at finite doping on square as well as triangular lattice geometries. These experimental advances pose the quest for theorists to analyze the ensuing data in order to gain insights into these prototypical, strongly correlated many-body systems. Here we employ machine learning techniques to optimize the mean-field parameters of a resonating valence bond (RVB) state through comparison with experimental data, thus determining a possible underlying simple model that is physically motivated and fully interpretable. We find that the resulting RVB states are capable of capturing two- as well as three-point correlations measured in experiments, even when they are not specifically used in the optimization. The analysis of the mean-field parameters and their doping dependence can be used to obtain physical insights and shed light on the nature of possible underlying quantum spin liquid states. Our results show that finite temperature data from Fermi-Hubbard quantum simulators can be well captured by RVB states. This work paves the way for a new, systematic analysis of data from numerical as well as quantum simulation of strongly correlated quantum many-body systems.