The recently discovered nickelate superconductor ${\mathrm{La}}_3{\mathrm{Ni}}_2{\mathrm{O}}_7$ has a high transition temperature near 80 K under pressure, providing an additional avenue for exploring unconventional superconductivity. Here, with state-of-the-art tensor-network methods, we study a bilayer ${t\text{−}J\text{−}J}_{\perp }$ model for ${\mathrm{La}}_3{\mathrm{Ni}}_2{\mathrm{O}}_7$ and find a robust $s$-wave superconductive (SC) order mediated by interlayer magnetic couplings. Large-scale density matrix renormalization group calculations find algebraic pairing correlations with Luttinger parameter $K_{\mathrm{SC}}\lesssim 1$. Infinite projected entangled-pair state method obtains a nonzero SC order directly in the thermodynamic limit, and estimates a strong pairing strength ${\overline{\mathrm{\Delta }}}_z\sim \mathcal{O}(0.1)$. Tangent-space tensor renormalization group simulations elucidate the temperature evolution of SC pairing and further determine a high SC temperature $T_c^{*}/J\sim \mathcal{O}(0.1)$. Because of the intriguing orbital selective behaviors and strong Hund’s rule coupling in the compound, ${t\text{−}J\text{−}J}_{\perp }$ model has strong interlayer spin exchange (while negligible interlayer hopping), which greatly enhances the SC pairing in the bilayer system. Such a magnetically mediated pairing has also been observed recently in the optical lattice of ultracold atoms. Our accurate and comprehensive tensor-network calculations reveal a robust SC order in the bilayer ${t\text{−}J\text{−}J}_{\perp }$ model and shed light on the pairing mechanism of the high-$T_c$ nickelate superconductor.

We speed up thermal simulations of quantum many-body systems in both one- (1D) and two-dimensional (2D) models in an exponential way by iteratively projecting the thermal density matrix $\widehat{\rho }=e^{-\beta \widehat{H}}$ onto itself. We refer to this scheme of doubling $\beta$ in each step of the imaginary time evolution as the exponential tensor renormalization group (XTRG). This approach is in stark contrast to conventional Trotter-Suzuki-type methods which evolve $\widehat{\rho }$ on a linear quasicontinuous grid in inverse temperature $\beta \equiv 1/T$. As an aside, the large steps in XTRG allow one to swiftly jump across finite-temperature phase transitions, i.e., without the need to resolve each singularly expensive phase-transition point right away, e.g., when interested in low-energy behavior. A fine temperature resolution can be obtained, nevertheless, by using interleaved temperature grids. In general, XTRG can reach low temperatures exponentially fast and, thus, not only saves computational time but also merits better accuracy due to significantly fewer truncation steps. For similar reasons, we also find that the series expansion thermal tensor network approach benefits in both efficiency and precision, from the logarithmic temperature scale setup. We work in an (effective) 1D setting exploiting matrix product operators (MPOs), which allows us to fully and uniquely implement non-Abelian and Abelian symmetries to greatly enhance numerical performance. We use our XTRG machinery to explore the thermal properties of Heisenberg models on 1D chains and 2D square and triangular lattices down to low temperatures approaching ground-state properties. The entanglement properties, as well as the renormalization-group flow of entanglement spectra in MPOs, are discussed, where logarithmic entropies (approximately $\ln \beta$) are shown in both spin chains and square-lattice models with gapless towers of states. We also reveal that XTRG can be employed to accurately simulate the Heisenberg $XXZ$ model on the square lattice which undergoes a thermal phase transition. We determine its critical temperature based on thermal physical observables, as well as entanglement measures. Overall, we demonstrate that XTRG provides an elegant, versatile, and highly competitive approach to explore thermal properties, including finite-temperature thermal phase transitions as well as the different ordering tendencies at various temperature scales for frustrated systems.

Frustrated magnets hold the promise of material realizations of exotic phases of quantum matter, but direct comparisons of unbiased model calculations with experimental measurements remain very challenging. Here we design and implement a protocol of employing many-body computation methodologies for accurate model calculations—of both equilibrium and dynamical properties—for a frustrated rare-earth magnet TmMgGaO₄ (TMGO), which explains the corresponding experimental findings. Our results confirm TMGO is an ideal realization of triangular-lattice Ising model with an intrinsic transverse field. The magnetic order of TMGO is predicted to melt through two successive Kosterlitz–Thouless (KT) phase transitions, with a floating KT phase in between. The dynamical spectra calculated suggest remnant images of a vanishing magnetic stripe order that represent vortex–antivortex pairs, resembling rotons in a superfluid helium film. TMGO therefore constitutes a rare quantum magnet for realizing KT physics, and we further propose experimental detection of its intriguing properties.

A linearized tensor renormalization group algorithm is developed to calculate the thermodynamic properties of low-dimensional quantum lattice models. This new approach employs the infinite time-evolving block decimation technique, and allows for treating directly the transfer-matrix tensor network that makes it more scalable. To illustrate the performance, the thermodynamic quantities of the quantum XY spin chain as well as the Heisenberg antiferromagnet on a honeycomb lattice are calculated by the linearized tensor renormalization group method, showing the pronounced precision and high efficiency.

Magic-angle twisted bilayer graphene has recently become a thriving material platform realizing correlated electron phenomena taking place within its topological flat bands. Several numerical and analytical methods have been applied to understand the correlated phases therein, revealing some similarity with the quantum Hall physics. In this work, we provide a Mott-Hubbard perspective for the TBG system. Employing the large-scale density matrix renormalization group on the lattice model containing the projected Coulomb interactions only, we identify a first-order quantum phase transition between the insulating stripe phase and the quantum anomalous Hall state with the Chern number of ±1. Our results not only shed light on the mechanism of the quantum anomalous Hall state discovered at three-quarters filling, but also provide an example of the topological Mott insulator, i.e., the quantum anomalous Hall state in the strong coupling limit.