arXiv Curator

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 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.

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.

La$_2$NiO$_4$ is a prototypical member of the Ruddlesden-Popper nickelate series that offers a valuable reference point for elucidating the key ingredients behind the intriguing properties of these systems. However, the structural and electronic properties of La$_2$NiO$_4$ under pressure and doping remain surprisingly underexplored. Here, we investigate these properties using density-functional-theory calculations. We find that its tetragonal $I4/mmm$ structure can be stabilized, not only under comparatively low pressures of $\sim$ 8 GPa, but also at ambient pressure via the partial substitution of La with Ba. Moreover, we show that the combined effects of Ba substitution and pressure leads to qualitative changes in the electronic structure towards the formal $d^{7.5}$ configuration of the superconducting bilayer nickelates. Further, while La$_2$NiO$_4$ can undergo a insulator-metal transition with pressure retaining G-type antiferromagnetic order, La$_{1.5}$Ba$_{0.5}$NiO$_4$ exhibits metallic behavior with an enhanced competition between different magnetic states. Our results thus offer new insights into the interplay of structure, doping, and magnetism across the Ruddlesden-Popper nickelate series.

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.

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.

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.

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.

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.

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.

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 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.

Time-dependent density functional theory (TDDFT) is a theory that describes the time evolution of quantum mechanical many-electron systems under the influence of external time-dependent electric and magnetic fields. INQ is a specially designed software to efficiently solve the real-time TDDFT equations on graphics processing units (GPUs), which aim to overcome the computational limitation of time and size scales of non-equilibrium quantum dynamics. In this work we will present an implementation of non-collinear TDDFT for the INQ code to simulate spin dynamics in real time and discuss the implementation of non-collinear magnetic effects into the code. We will discuss the implementation of exchange-correlation magnetic fields, spin-orbit coupling, and the interaction between the electronic system and external magnetic fields. We will then consider several prototypical examples of spin dynamics in magnetic clusters and solids after light excitation. Potential applications range from the study of real-time dynamics of magnons to ultrafast spin dynamics under linear and circularly polarized laser excitation, as well as spectroscopic signatures such as magnetic circular dichroism and pump-probe Kerr rotation.

Altermagnetism has recently emerged as a distinct and fundamental class of magnetic order. Exploring its interplay with quantum phenomena such as unconventional superconductivity, density-wave instabilities, and many-body effects represents a compelling frontier. In this work, we theoretically confirm the presence of high-temperature metallic altermagnetism in KV$_2$Se$_2$O. We demonstrate that the anomalous metal-insulator-metal transition arises from a Lifshitz transition associated with Fermi surface reconstruction. The previously reported spin-density wave gap is found to lie below the Fermi level in our study and is now recognized to be attributed to the V-shaped density of states, originating from orbital-selective and sublattice-resolved half-metal-like behavior on a specific V atom. Furthermore, we identify the instability from the nesting of spin-momentum-locked two-dimensional Fermi surfaces, which induces the SDW state. These findings position KV$_2$Se$_2$O as a promising platform for investigating the interplay among altermagnetism, unconventional superconductivity, and density-wave order.

While a predictive theory for unconventional superconductivity in Fe-based superconductors remains elusive, an extensively debated aspect is the interaction between phonons and strongly correlated electrons, and its potential role in the pairing mechanism. Here, through the combination of first principles dynamical mean field theory calculations and epitaxial growth of the single-layer FeX (X=Se, S, Te) on SrTiO$_3$ (STO)(001) substrate, which facilitates the controlled distortion of the FeX$_4$ tetrahedron, we demonstrate an unique superconducting dome where the superconducting gap peaks at a `magic' angle of the FeX$_4$ tetrahedron and the electron-phonon coupling (EPC) for the A$_{1g}$ mode is maximized for the FeSe film. Our findings uncover a significant role of electronic correlations in strengthening Cooper pairing in unconventional superconductors by enhancing EPC.

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 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.

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.

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 investigate the magnetotransport properties of a diluted half-filled one-band Hubbard model with second-nearest-neighbor hopping on a simple cubic lattice, aiming to explore the possibility of metallicity in diluted antiferromagnetic systems. Our semiclassical Monte Carlo (s-MC) calculations reveal an antiferromagnetic metallic regime in diluted correlated materials. This unexpected metallic regime naturally leads to a central question: how does the introduction of dilution into an antiferromagnetic material, especially with long-range magnetic interactions, induce metallicity -- a feature not commonly associated with antiferromagnets? To address this question, we demonstrate that when the on-site repulsion strength ($U$) is set to zero on a percentage of the sites (site dilution), the insulating state weakens due to percolative conduction among the diluted sites at low temperatures. Remarkably, this occurs without any significant alteration to the underlying long-range antiferromagnetic (AF) ordering in the system, thereby providing a pathway to realize antiferromagnetic metals. In addition, we show how the sublattice-dependent hopping can be exploited to engineer spin-polarized half-metallic antiferromagnets. Overall, our numerical results collectively provide a basis for understanding the combined effect of site dilution and competing interactions, which will assist in the design of new antiferromagnetic metals for future spintronic applications.

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.

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 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.

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.

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.

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.

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.

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.

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 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.

Stimulated by recent interest in altermagnets, a novel class of antiferromagnets with macroscopic time-reversal symmetry breaking, we investigate the coexistence of altermagnetism and superconductivity. By developing a Ginzburg--Landau theory based on microscopic models, we show that a phase-modulated Fulde--Ferrell superconducting state is stabilized via altermagnetic spin splitting, in contrast to the typical amplitude-modulated states that occur under a uniform Zeeman field. We apply our framework to different models: a two-sublattice model with altermagnetic order, a continuum model with an anisotropic Zeeman field mimicking altermagnetic spin splitting, and a conventional square-lattice model with two kinds of anisotropic Zeeman fields. Our results indicate that the multi-sublattice structure is crucial for realizing this exotic superconductivity, and highlight spin-split altermagnets as a promising platform for exploring phase-modulated superconductivity without external magnetic fields.

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$.

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 a dissipative version of the Anderson Impurity model, where an interacting impurity is coupled to a fermionic reservoir and exposed to Markovian dissipation in the form of two-body losses. Using a self-consistent hybridization expansion based on the Non-Crossing Approximation (NCA) we compute the dynamics of the impurity, its steady-state and spectral function. We show that the interplay between strong Coulomb repulsion and correlated dissipation gives rise to robust signatures of Kondo physics both at weak and strong losses. These include a strongly suppressed spin relaxation rate, displaying a characteristic Kondo-Zeno crossover and a spectral function where doublon band is quickly destroyed by dissipation while the coherent Kondo peak remains visible for weak losses, then disappears at intermediate values and finally re-emerge as the system enters in the Kondo-Zeno regime. As compared to the case of single particle losses we show that two-body dissipation protects Kondo physics. The picture obtained with NCA is confirmed by numerical simulations of exact dynamics on finite-size chains. We interpret these results using a dissipative Schrieffer-Wolff transformation, which leads to an effective Kondo model with residual impurity-bath losses which are suppressed by strong correlations or strong losses.

High magnetic fields have proven instrumental in exploring the physical properties of condensed matter, leading to groundbreaking discoveries such as the quantum Hall effect in 2D heterostructures and quantum oscillations in cuprate superconductors. The ability to conduct precise measurements at progressively higher magnetic fields continues to push the frontiers of knowledge and enable new discoveries. In this work, we present the development of a microwave technique for performing two-point transport measurements in semi-destructive pulsed magnetic fields (up to 125 T) and at low temperatures (down to 1.5 K) with unprecedented sensitivity. This new setup was tested on a variety of samples. We present results on the metal-insulator transition in InAs and we report notably the first observation of Shubnikov-de-Haas oscillations in WTe$_{2}$ at magnetic fields beyond 100 T.

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.

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.

This study provides a first-principles insight into half-Heusler CrMnSb to understand its deviation from the conventional Slater-Pauling semiconducting behavior. CrMnSb, having a valence electron count of 18, has been proposed to exhibit compensated ferrimagnetic character instead of the expected nonmagnetic semiconducting ground state. As half-Heusler systems with a valence electron count of 18 are not known to exhibit magnetic ordering, we have investigated the electronic and magnetic properties of CrMnSb using a combination of density functional theory and Green's function-based multiple-scattering theory. We show that, despite satisfying the 18 valence electron Slater-Pauling rule, CrMnSb does not exhibit ground-state nonmagnetic semiconducting behavior. Instead, it reveals a half-metallic, fully compensated ferrimagnetic ground state. This anomaly originates from the presence of localized sublattice moments, resulting from antiparallel alignment between Cr and Mn sublattices, which enforces half-metallic ferrimagnetism despite its ideal 18 valence electron count.

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.

Geometric frustration is a key ingredient in the emergence of exotic states of matter, such as the quantum spin liquid in Mott insulators. While there has been intense interest in experimentally tuning frustration in candidate materials, achieving precise and continuous control has remained a major hurdle -- particularly in accessing the properties of the ideally frustrated lattice. Here, we show that large, finely controlled anisotropic strains can effectively tune the degree of geometric frustration in the Mott insulating $\kappa$-(ET)$_2$Cu$_2$(CN)$_3$ -- a slightly anisotropic triangular-lattice quantum magnet. Using thermodynamic measurements of the elastocaloric effect, we experimentally map out a temperature-strain phase diagram that captures both the ground state of the isotropic lattice and the less frustrated parent state. Our results provide a new benchmark for calculations of the triangular-lattice Hubbard model as a function of frustration and highlight the power of lattice engineering as a route to realizing perfectly frustrated quantum materials.

Altermagnetism is a collinear compensated magnetically-ordered phase with a d, g or i-wave anisotropy and alternating spin polarization of the electronic structure in the position and momentum space. Its recent discovery was in part motivated by the research of compensated magnets towards highly scalable spintronic technologies. Simultaneously, altermagnetism shares the anisotropic higher-partial-wave nature of ordering with unconventional superfluid phases which have been at the forefront of research for the past several decades. These examples illustrate the interest in altermagnetism from a broad range of science and technology perspectives. After summarizing the diverse research context, we turn the focus of this review to the symmetry, microscopy and spectroscopy signatures of altermagnetism. We start from the description of spontaneously broken and retained symmetries which delineate the compensated altermagnetic ordering as a distinct magnetic phase. Next we focus on microscopic signatures and ordering mechanism of the altermagnetic phase. We highlight crystal-structure realizations of a characteristic ferroic order of anisotropic higher-partial-wave components of atomic-scale spin densities in altermagnets, ranging from weakly-interacting metals to strongly correlated insulators. The symmetry and microscopy signatures of altermagnetism are directly reflected in spin-dependent electronic spectra and responses. We review salient band-structure features originating from the altermagnetic ordering, and from its interplay with spin-orbit coupling and topological phenomena. Throughout the review we compare altermagnetism to traditional ferromagnetism and Neel antiferromagntism, and to the currently intensely explored magnetic phases with non-collinear symmetry-protected compensated spin orders. We accompany the theoretical discussions by references to relevant experiments.

We investigate the magnetotransport properties of a diluted half-filled one-band Hubbard model with second-nearest-neighbor hopping on a simple cubic lattice, aiming to explore the possibility of metallicity in diluted antiferromagnetic systems. Our semiclassical Monte Carlo calculations reveal an antiferromagnetic metallic regime in diluted correlated materials. This unexpected metallic regime naturally leads to a central question: how does the introduction of dilution into an antiferromagnetic material, especially with long-range magnetic interactions, induce metallicity -- a feature not commonly associated with antiferromagnets? To address this question, we demonstrate that when the on-site repulsive Hubbard interaction strength is set to zero on a percentage of the sites (site dilution), the insulating state weakens due to percolative conduction among the diluted sites at low temperatures. Remarkably, this occurs without any significant alteration to the underlying long-range antiferromagnetic ordering in the system, thereby providing a pathway to realize antiferromagnetic metals. In addition, we show how the sublattice-dependent hopping can be exploited to engineer spin-polarized half-metallic antiferromagnets. Overall, our numerical results collectively provide a basis for understanding the combined effect of site dilution and competing interactions, which will assist in the design of new antiferromagnetic metals for future spintronic applications.

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 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.

The interplay between symmetry and topology in magnetic materials makes it possible to engineer exotic phases and technologically useful properties. A key requirement for these pursuits is achieving control over local crystallographic and magnetic structure, usually through sample morphology (such as synthesis of bulk crystals versus thin-films) and application of magnetic or electric fields. Here we show that V$_{1/3}$NbS$_2$ can be crystallized in two ordered superlattices, distinguished by the periodicity of out-of-plane magnetic intercalants. Whereas one of these structures is metallic and displays the hallmarks of altermagnetism, the other superlattice, which has not been isolated before in this family of intercalation compounds, is a semimetallic noncollinear antiferromagnet that may enable access to topologically nontrivial properties. This observation of an unconventional superlattice structure establishes a powerful route for tailoring the tremendous array of magnetic and electronic behaviors hosted in related materials.

Excited eigenstates are crucial to understand the dynamics of quantum many-body systems. Tensor network states are one of the workhorses to compute ground states of many-body systems, yet the accurate computation of excited eigenstates is still challenging. Here, we develop a combination of the inexact Lanczos method, which aims at efficiently computing excited states, to tree tensor network states (TTNSs). We demonstrate our approach by computing excited vibrational states for three challenging problems: (1) 84 states in different energy intervals of acetonitrile (12-dimensional), (2) Fermi resonance states of the fluxional Zundel ion (15-dimensional), and (3) selected excited states of the fluxional and very correlated Eigen ion (33-dimensional). The proposed TTNS inexact Lanczos method is directly applicable to other quantum many-body systems.

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.

We find an efficient approach to approximately convert matrix product states (MPSs) into restricted Boltzmann machine wave functions consisting of a multinomial hidden unit through a canonical polyadic (CP) decomposition of the MPSs. This method allows us to generate well-behaved initial neural network quantum states for quantum many-body ground-state calculations in polynomial time of the number of variational parameters and systematically shorten the distance between the initial states and the ground states with increasing the rank of the CP decomposition. We demonstrate the efficiency of our method by taking the transverse-field Ising model as an example and discuss possible applications of our method to more general quantum many-body systems in which the ground-state wave functions possess complex nodal structures.

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.

Altermagnetism, featuring alternating spin structures in reciprocal space, has sparked growing interest. Here, we predict novel real-space alternative piezomagnetic and piezoelectric responses in an emerging altermagnetic family of Lieb lattices, specifically transition-metal chalcogenides M2WS4 (M = Mn, Fe, Co). The unique S4T crystal-spin symmetry leads to distinct magnetic and electric responses depending on the direction of applied stress. When subjected to axial stress, they exhibit a giant piezomagnetic response, which is about one to two orders of magnitude larger than that of most piezomagnetic materials, while the residual C2 symmetry suppresses the piezoelectric effect. In contrast, diagonal stress induces an imbalance of oppositely aligned electric dipole moments and a significant piezoelectric response, while in-plane mirror symmetry inhibits the piezomagnetic effect. This alternative piezoresponse offers an unprecedented opportunity to precisely control electric and magnetic properties independently, opening new avenues for altermagnetic materials in high-fidelity multifunctional memory and sensor applications.

Graph foundation models using graph neural networks promise sustainable, efficient atomistic modeling. To tackle challenges of processing multi-source, multi-fidelity data during pre-training, recent studies employ multi-task learning, in which shared message passing layers initially process input atomistic structures regardless of source, then route them to multiple decoding heads that predict data-specific outputs. This approach stabilizes pre-training and enhances a model's transferability to unexplored chemical regions. Preliminary results on approximately four million structures are encouraging, yet questions remain about generalizability to larger, more diverse datasets and scalability on supercomputers. We propose a multi-task parallelism method that distributes each head across computing resources with GPU acceleration. Implemented in the open-source HydraGNN architecture, our method was trained on over 24 million structures from five datasets and tested on the Perlmutter, Aurora, and Frontier supercomputers, demonstrating efficient scaling on all three highly heterogeneous super-computing architectures.

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.

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.

Using our developed thermal tensor-network approach, we investigate the spin Seebeck effect (SSE) of the triangular-lattice quantum antiferromagnet hosting spin supersolid phase. We focus on the low-temperature scaling of normalized spin current $\tilde{I}_S$ through the interface, and benchmark our approach on 1D Heisenberg chain, with the negative spinon spin current reproduced. In addition, we find an algebraic scaling $\tilde{I}_S \sim T^{\alpha}$ with varying exponent $\alpha$ in the Tomonaga-Luttinger liquid phase. On the triangular lattice, spin frustration dramatically enhances the low-temperature SSE, with characteristic spin-current signatures distinguishing different magnetic phases. Remarkably, we discover a persistent spin current $\tilde{I}_S$ in the spin supersolid phase, which saturates to a non-zero value in the low-temperature limit and can be ascribed to the Goldstone-mode-mediated spin supercurrents. Moreover, a universal scaling $\tilde{I}_S \sim T^{d/z}$ is found at the U(1)-symmetric polarization quantum critical points. These distinct quantum spin transport traits - particularly the sign reversal and characteristic temperature dependence in SSE - provide sensitive experimental probes for investigating spin-supersolid compounds such as Na$_2$BaCo(PO$_4$)$_2$. Moreover, our results establish spin supersolids as a tunable quantum platform for spin caloritronics in the ultralow-temperature regime.

Obtaining accurate representations of the eigenstates of an array of coupled superconducting qubits is a crucial step in the design of circuit quantum electrodynamics (QED)-based quantum processors. However, exact diagonalization of the device Hamiltonian is challenging for system sizes beyond tens of qubits. Here, we employ a variant of the density matrix renormalization group (DMRG) algorithm, DMRG-X, to efficiently obtain localized eigenstates of a 2D transmon array without the need to first compute lower-energy states. We also introduce MTDMRG-X, a new algorithm that combines DMRG-X with multi-target DMRG to efficiently compute excited states even in regimes with strong eigenstate hybridization. We showcase the use of these methods for the analysis of long-range couplings in a multi-transmon Hamiltonian including qubits and couplers, and we discuss eigenstate localization. These developments facilitate the design and parameter optimization of large-scale superconducting quantum processors.

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.

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.

The unique physical properties of altermagnets, when transplanted to photonic systems, are anticipated to offer a new degree of freedom for engineering electromagnetic waves. Here, we show that a photonic analogue of altermagnetism can be mimicked in photonic crystals, where engineered photonic crystals can host spin space group symmetries. Our approach allows for the creation of spin-split bands and the corresponding transport properties provide an effective platform for circularly polarized light isolation without the need of geometrodynamic spin-orbit interaction. Beyond the concurrent solid-state materials, we anticipate our work to offer photonic crystals as a versatile platform to test the spin-split band properties and inspire optical designs for photospintronic applications.

We introduce a holographic framework for the entanglement Hamiltonian in symmetry-protected topological (SPT) phases with area-law entanglement, whose reduced density matrix $\rho \propto e^{-H_e}$ can be treated as a lower-dimensional mixed state. By replicating $\rho$, we reconstruct the fixed-point SPT wavefunction, establishing an exact correspondence between the bulk strange correlator of the (d+1)-dimensional SPT state and the twisted R\'enyi-N operator of the d-dimensional reduced density matrix. Notably, the reduced density matrix exhibits long-range or quasi-long-range order along the replica direction, revealing a universal entanglement feature in SPT phases. As a colloary, we generalized the framework of twisted R\'enyi-N correlator to thermal states and open quantum systems, providing an alternative formulation of the Lieb-Schultz-Mattis theorem, applicable to both closed and open systems. Finally, we extend our protocol to mixed-state SPT phases and introduce new quantum information metrics -- twisted R\'enyi-N correlators of the surgery operator -- to characterize the topology of mixed states.

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.

Rare earth (RE)-free permanent magnets, as alternative substitutes for RE-containing magnets for sustainable energy technologies and modern electronics, have attracted considerable interest. We performed a comprehensive search for new hard magnetic materials in the ternary Fe-Co-Zr space by leveraging a scalable, machine learning-assisted materials discovery framework running on GPU-enabled exascale computing resources. This framework integrates crystal graph convolutional neural network (CGCNN) machine learning (ML) method with first-principles calculations to efficiently navigate the vast composition-structure space. The efficiency and accuracy of the ML approach enable us to reveal 9 new thermodynamically stable ternary Fe-Co-Zr compounds and 81 promising low-energy metastable phases with their formation energies within 0.1 eV/atom above the convex hull. The predicted compounds span a wide range of crystal symmetries and magnetic behaviors, providing a rich platform for tuning functional properties. Based on the analysis of site-specific magnetic properties, we show that the Fe6Co17Zr6 compound obtained from our ML discovery can be further optimized by chemical doping. Chemical substitutions lead to a ternary Fe5Co18Zr6 phase with a strong anisotropy of K1 = 1.1 MJ/m3, and a stable quaternary magnetic Fe5Co16Zr6Mn4 compound.

We introduce a modified 2D toric code Hamiltonian that exhibits explicit anyon confinement along a single spatial direction. By bounding the motion of these confined anyons, we obtain dipolar excitations with restricted mobility. We analyze the resulting logical operators, whose existence depends on the system size, as well as the structure of gapped boundaries and a tensor network representation of the ground state. Furthermore, when confinement is enforced in both directions, fractal-like excitations emerge, resulting in unpaired logical operators. We extend our construction to 3D models, such as the surface code and the X-cube model, leading to novel dipole-loop and dipole-planon excitations that arise from bounding confined excitations. These modifications are implemented through group cohomological twistings--projective representations of finite groups--with most examples based on Z2xZ2.

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.

The recent proposal of altermagnetism has drawn widespread attention to antiferromagnet (AFM) exhibiting spin splitting, extending beyond the realm of sign-alternating spin splitting in momentum space protected solely by rotational symmetry. Herrin, we propose a shear-strain strategy that enables significant modulation of d-wave altermagnets into an non-alter compensated magnets. A comprehensive analysis combining the magnetic moment compensation characteristics of opposite spin sublattices with the distribution of spin-resolved conduction channels in momentum space under the [001] crystal orientation reveals that shear strain breaks the rotational symmetry of alternatmagnets. To explore the application potential of non-alter compensated magnets, we designe RuO2/TiO2/RuO2 magnetic tunnel junctions (MTJ) with three crystallographic orientations ((001), (110), (100)) and investigated their transport properties under shear strain. This non-alter electronic structure not only enhances the tunneling magnetoresistance (TMR) in spin-split paths of intrinsic RuO2 (226% to 431%) but also enables substantial TMR in spin-degenerate paths (from 0-88%)). Our work provides guidelines for broadening magnetic materials and device platforms.

Variational tensor network optimization has become a powerful tool for studying classical statistical models in two dimensions. However, its application to three-dimensional systems remains limited, primarily due to the high computational cost associated with evaluating the free energy density and its gradient. This process requires contracting a triple-layer tensor network composed of a projected entangled pair operator and projected entangled pair states. In this paper, we employ a split corner-transfer renormalization group scheme tailored for the contraction of such triple-layer networks, which reduces the computational complexity from $\mathcal{O}(D^{12})$ in conventional methods to $\mathcal{O}(D^{9})$. Through numerical benchmarks on the three-dimensional classical Ising model, we demonstrate that the proposed scheme achieves numerical results comparable to the most recent Monte Carlo simulations, while also providing a substantial speedup over previous variational tensor network approaches. This makes this method well-suited for efficient gradient-based optimization in three-dimensional tensor network simulations.

This study investigates the size-dependent mechanical behavior of the HfNbTaTiZr refractory high-entropy alloy (RHEA) under uniaxial tension, with a focus on the effects of random solid-solution (RSS) and chemical short-range order (CSRO). A machine learning framework is developed to accelerate the parameterization of interatomic force fields (FFs), enabling molecular dynamics simulations of three nanocrystalline models: (i) a meta-atom (MA) mode representing the RHEA as a hypothetical sing-element system with averaged properties, (ii) a quinary RSS model with randomly distributed constituent atoms, and (iii) a Monte Carlo (MC) model with internal CSRO. The results reveal that RSS enhances strength, while CSRO reduces flow stress level but improves strain hardening and failure resistance. A transition from Hall-Petch (HP) strengthening to inverse Hall-Petch (IHP) softening is observed, with CSRO suppressing this transition. The underlying plastic mechanisms (i.e., dislocation slip, deformation twinning, phase transformation and grain boundary movements) are analyzed from both nanostructural and energetic perspectives. Theoretical models are established to describe the size-dependent yield strength and predict the critical grain size. Additionally, the contributions of different plastic mechanisms to the overall stress response are separately quantified. These findings provide new insights into the design and performance optimization of RHEAs through nanostructural engineering.

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.

Recent advances in Rydberg tweezer arrays bring novel opportunities for programmable quantum simulations beyond previous capabilities. In this work, we investigate a bosonic t-J-V model currently realized with Rydberg atoms. Through large-scale quantum Monte Carlo simulations, we uncover an emergent double supersolid (DSS) phase with the coexistence of two superfluids and crystalline order. Tunable long-range tunneling and repulsive hole-hole interactions enable a rich phase diagram featuring a double superfluid phase, a DSS phase, and an antiferromagnetic insulator. Intriguingly, within the DSS regime we observe an unconventional thermal enhancement of crystalline order. Our results establish the bosonic t-J-V model as a promising and experimentally accessible platform for exploring exotic quantum phases in Rydberg atom arrays.

This work focuses on optimizing the gates of a quantum circuit with a given topology to approximate the unitary time evolution governed by a Hamiltonian. Recognizing that unitary matrices form a mathematical manifold, we employ Riemannian optimization methods -- specifically the Riemannian trust-region algorithm -- which involves second derivative calculations with respect to the gates. Our key technical contribution is a matrix-free algorithmic framework that avoids the explicit construction and storage of large unitary matrices acting on the whole Hilbert space. Instead, we evaluate all quantities as sums over state vectors, assuming that these vectors can be stored in memory. We develop HPC-optimized kernels for applying gates to state vectors and for the gradient and Hessian computation. Further improvements are achieved by exploiting sparsity structures due to Hamiltonian conservation laws, such as parity conservation, and lattice translation invariance. We benchmark our implementation on the Fermi-Hubbard model with up to 16 sites, demonstrating a nearly linear parallelization speed-up with up to 112 CPU threads. Finally, we compare our implementation with an alternative matrix product operator-based approach.

In our previous article [arXiv:2307.12552], we introduced local topological order (LTO) axioms for quantum spin systems which allowed us to define a physical boundary manifested by a net of boundary algebras in one dimension lower. This gives a formal setting for topological holography, where the braided tensor category of DHR bimodules of the physical boundary algebra captures the bulk topological order. In this article, we extend the LTO axioms to quantum spin systems equipped with a topological boundary, again producing a physical boundary algebra for the bulk-boundary system, whose category of (topological) boundary DHR bimodules recovers the topological boundary order. We perform this analysis in explicit detail for Levin-Wen and Walker-Wang bulk-boundary systems. Along the way, we introduce a 2D braided categorical net of algebras built from a unitary braided fusion category (UBFC). Such nets arise as boundary algebras of Walker-Wang models. We consider the canonical state on this braided categorical net corresponding to the standard topological boundary for the Walker-Wang model. Interestingly, in this state, the cone von Neumann algebras are type I with finite dimensional centers, in contrast with the type II and III cone von Neumann algebras from the Levin-Wen models studied in [arXiv:2307.12552]. Their superselection sectors recover the underlying unitary category of our UBFC, and we conjecture the superselection category also captures the fusion and braiding.

The critical behavior of the classical Ising model on a three-dimensional fractal lattice with Hausdorff dimension $d_H = \ln32 / \ln4 = 2.5$ is investigated using the higher-order tensor renormalization group (HOTRG) method. We determine the critical temperature $T_c \approx 2.65231$ and the critical exponents for magnetization $\beta \approx 0.059$ and field response $\delta \approx 35$. Unlike a previously studied 2D fractal with $d_H \approx 1.792$, the specific heat for this 3D fractal exhibits a divergent singularity at $T_c$. The results are compared with those for regular lattices and other fractal structures to elucidate the role of dimensionality in critical phenomena.

We study Ruelle-Pollicott resonances of translationally invariant magnetization-conserving qubit circuits via the spectrum of the quasi-momentum-resolved truncated propagator of extensive observables. Diffusive transport of the conserved magnetization is reflected in the Gaussian quasi-momentum $k$ dependence of the leading eigenvalue (Ruelle-Pollicott resonance) of the truncated propagator for small $k$. This, in particular, allows us to extract the diffusion constant. For large $k$, the leading Ruelle-Pollicott resonance is not related to transport and governs the exponential decay of correlation functions. Additionally, we conjecture the existence of a continuum of eigenvalues below the leading diffusive resonance, which governs non-exponential decay, for instance, power-law hydrodynamic tails. We expect our conclusions to hold for generic systems with exactly one U(1) conserved quantity.

We report a scalable and cost-effective method for fabricating high-performance Bragg mirrors using a bottom-up approach that combines evaporation-induced self-assembly (EISA) and dip-coating. The photonic crystals are composed of alternating mesoporous SiO$_2$ and dense TiO$_2$ layers, providing a high refractive index contrast ($\sim$0.8). This enables strong reflectance (up to 96%) with as few as five bilayers and precise control of the photonic stop band across the visible spectrum by simply adjusting the deposition parameters. Integration of a thin film of the two-dimensional perovskite (PEA)$_2$PbI$_4$ leads to strong light--matter coupling at room temperature. Angle-resolved reflectance and photoluminescence measurements reveal the formation of upper and lower polariton branches, with a Rabi splitting of 90 meV. The observed polaritonic dispersion is well described by a two-level system and Green's function formalism. This work demonstrates an efficient strategy for constructing tunable optical cavities using simple solution-based methods. The combination of high optical quality, spectral tunability, and strong coupling performance positions this platform as a promising candidate for low-threshold polariton lasers, nonlinear optics, and integrated optoelectronic devices.

Using the developed thermal tensor-network approach, we investigate the spin Seebeck effect (SSE) of the triangular-lattice quantum antiferromagnet hosting spin supersolid phase. We focus on the low-temperature scaling behaviors of the normalized spin current across the interface. Using the 1D Heisenberg chain as a benchmark system, we observe a negative spinon spin current exhibiting algebraic scaling $T^{\alpha}$, with exponent $\alpha$, in the Tomonaga-Luttinger liquid phase. On the triangular lattice, spin frustration dramatically enhances the low-temperature SSE, with distinct spin-current signatures - particularly the sign reversal and characteristic temperature dependence - distinguishing different spin states. Remarkably, we discover a persistent, negative spin current in the spin supersolid phase. It saturates to a non-zero value in the low-temperature limit, and can be ascribed to the Goldstone-mode-mediated spin supercurrents. Moreover, a universal scaling $T^{d/z}$ is found at the U(1)-symmetric polarization quantum critical points. These distinct quantum spin transport traits provide sensitive probes for spin-supersolid compounds such as Na$_2$BaCo(PO$_4$)$_2$. Furthermore, our results establish spin supersolids as a tunable quantum platform for spin caloritronics in the ultralow-temperature regime.

We present an extended formulation of the non-Markovian stochastic Schr\"odinger equation with complex frequency modes (extended cNMSSE), enabling the treatment of open quantum system dynamics under general spectral densities. This extension employs complete non-exponential basis expansions for the bath correlation functions, thereby generalizing the applicability of the cNMSSE framework beyond environments with Debye-type spectral structures. By preserving the wavefunction-based description and favorable linear-scaling properties, the extended cNMSSE offers an efficient and flexible approach to simulate non-Markovian quantum dynamics. The method is implemented within a pseudo-Fock space using conventional ladder operators and solved numerically via matrix product state (MPS) techniques. Benchmark simulations on four representative cases, including discrete spectra, Ohmic spectra with exponential and algebraic cutoffs, and critically damped Brownian spectral densities, demonstrate excellent agreement with results of hierarchy of forward-backward stochastic Schr\"odinger equations (HFB-SSE) and extended hierarchical equation of motion (HEOM). The extended cNMSSE thus provides a robust and scalable framework for accurate simulations of non-Markovian open quantum systems with general environments.

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 consider the Blume--Capel model in the scaling limit to realize the thermal perturbation of the tricritical Ising fixed point. We develop a numerical scaling limit extrapolation for one-point functions and R\'enyi entropies in the ground state. In a mass quench scenario, we found long-lived oscillations despite the absence of explicit spin-flip symmetry breaking or confining potential. We construct form factors of branch-point twist fields in the paramagnetic phase. In the scaling limit of small quenches, we verify form factor predictions for the energy density and leading magnetic field using the dynamics of one-point functions, and branch-point twist fields using the dynamics of R\'enyi entropies.

Motivated by experiments on Rydberg atom arrays, we explore the properties of uniform quantum superpositions of kagome dimer configurations and construct an efficient algorithm for experimentally producing them. We begin by considering the thin cylinder limit, where these states have simple descriptions. We then develop a matrix product representation of the states on arbitrary cylinders, which leads to a natural protocol to efficiently grow them. We explain how our approach can be adapted to other quantum computing hardware.

Free theories are landmarks in the landscape of quantum field theories: their exact solvability serves as a pillar for perturbative constructions of interacting theories. Fuzzy sphere regularization, which combines quantum Hall physics with state-operator correspondence, has recently been proposed as a promising framework for simulating three-dimensional conformal field theories (CFTs), but so far it has not provided access to free theories. We overcome this limitation by designing a bilayer quantum Hall system that hosts an Ising tricritical point -- a nontrivial fixed point where first-order and second-order transitions meet -- which flows to the conformally coupled scalar theory in the infrared. The critical energy spectrum and operator structure match those at the Gaussian fixed point, providing nonperturbative evidence for the emergence of a free scalar CFT. Our results expand the landscape of CFTs realizable on the fuzzy sphere and demonstrate that even free bosonic theories -- previously inaccessible -- can emerge from interacting electrons in this framework.

False vacuum decay in scalar quantum field theory (QFT) plays a central role in cosmology and particle physics but has remained intractable for classical simulation due to its non-perturbative, highly entangled dynamics. We present a novel method for Hamiltonian simulation of false-vacuum decay in a self-interacting scalar QFT, based on a spatial discretisation into a lattice of bosonic modes. Each mode has a continuous or quasi-continuous local Hilbert space whose dimension is large enough to capture the self-interactions in the Hamiltonian. We refer to these local modes generically as ``qumodes'', echoing the ideal case in which the lattice would be realised on a continuous-variable quantum computer. Here we show that the same framework can be implemented on a tensor network, in which case the qumodes are quasi-continuous, corresponding to the local sites of the network having a very high dimensional Hilbert space. Using a time-evolving block decimation algorithm with Trotterised time evolution, we capture the real-time dynamics of the scalar field. In (1+1)-dimensions, we initialise a metastable vacuum and observe the coherent entanglement growth that seeds true-vacuum bubble nucleation. Our approach establishes the qumode network as a scalable framework for non-equilibrium scalar QFT phenomena and paves the way for higher-dimensional studies and continuous-variable quantum computing implementations.

We propose a quantum error mitigation method termed self-mitigation, which is comparable with zero-noise extrapolation, to achieve quantum utility on near-term, noisy quantum computers. We investigate the effectiveness of several quantum error mitigation strategies, including self-mitigation, by simulating quantum quench dynamics for Heisenberg spin chains with system sizes up to 104 qubits using IBM quantum processors. In particular, we discuss the limitations of zero-noise extrapolation and the advantages offered by self-mitigation at a large scale. The self-mitigation method shows stable accuracy with the large systems of 104 qubits with more than 3,000 CNOT gates. Also, we combine the discussed quantum error mitigation methods with practical entanglement entropy measuring methods, and it shows a good agreement with the theoretical estimation. Our study illustrates the usefulness of near-term noisy quantum hardware in examining the quantum quench dynamics of many-body systems at large scales, and lays the groundwork for surpassing classical simulations with quantum methods prior to the development of fault-tolerant quantum computers.

Photonic Integrated Circuits (PICs) provide superior speed, bandwidth, and energy efficiency, making them ideal for communication, sensing, and quantum computing applications. Despite their potential, PIC design workflows and integration lag behind those in electronics, calling for groundbreaking advancements. This review outlines the state of PIC design, comparing traditional simulation methods with machine learning approaches that enhance scalability and efficiency. It also explores the promise of quantum algorithms and quantum-inspired methods to address design challenges.

Quantum state tomography (QST) is a fundamental technique for estimating the state of a quantum system from measured data and plays a crucial role in evaluating the performance of quantum devices. However, standard estimation methods become impractical due to the exponential growth of parameters in the state representation. In this work, we address this challenge by parameterizing the state using a low-rank block tensor train decomposition and demonstrate that our approach is both memory- and computationally efficient. This framework applies to a broad class of quantum states that can be well approximated by low-rank decompositions, including pure states, nearly pure states, and ground states of Hamiltonians.

Existing quantum state tomography methods are limited in scalability due to their high computation and memory demands, making them impractical for recovery of large quantum states. In this work, we address these limitations by reformulating the maximum likelihood estimation (MLE) problem using the Burer-Monteiro factorization, resulting in a non-convex but low-rank parameterization of the density matrix. We derive a fully unconstrained formulation by analytically eliminating the trace-one and positive semidefinite constraints, thereby avoiding the need for projection steps during optimization. Furthermore, we determine the Lagrange multiplier associated with the unit-trace constraint a priori, reducing computational overhead. The resulting formulation is amenable to scalable first-order optimization, and we demonstrate its tractability using limited-memory BFGS (L-BFGS). Importantly, we also propose a low-memory version of the above algorithm to fully recover certain large quantum states with Pauli-based POVM measurements. Our low-memory algorithm avoids explicitly forming any density matrix, and does not require the density matrix to have a matrix product state (MPS) or other tensor structure. For a fixed number of measurements and fixed rank, our algorithm requires just $\mathcal{O}(d \log d)$ complexity per iteration to recover a $d \times d$ density matrix. Additionally, we derive a useful error bound that can be used to give a rigorous termination criterion. We numerically demonstrate that our method is competitive with state-of-the-art algorithms for moderately sized problems, and then demonstrate that our method can solve a 20-qubit problem on a laptop in under 5 hours.

QMetro++ is a Python package containing a set of tools dedicated to identifying optimal estimation protocols that maximize quantum Fisher information (QFI). Optimization can be performed for an arbitrary arrangement of input states, parameter encoding channels, noise correlations, control operations and measurements. The use of tensor networks and iterative see-saw algorithm allows for an efficient optimization even in the regime of large number of channel uses ($N\approx100$). Additionally, the package comes with an implementation of the recently developed methods for computing fundamental upper bounds on QFI, which serve as benchmarks of optimality of the outcomes of numerical optimization. All functionalities are wrapped up in a user-friendly interface which enables defining strategies at various levels of detail.