DMRG Preprints

Recently there has been considerable excitement surrounding the promising realization of high-spin Kitaev material, such as the quasi-2D compound CrI$_3$ and CrGeTe$_3$. However, the stability of quantum spin liquids (QSL) against single ion anisotropy (SIA) in these materials and the global quantum phase diagram of the extended spin-3/2 Kitaev model with finite SIA remain unclear. In this study, we perform large-scale density matrix renormalization group (DMRG) to explore the quantum phase diagram of the generalized spin-3/2 Kitaev-Heisenberg (K-H) model accompanied with SIA $A_c$. In the $A_c=0$ limit, the spin-3/2 K-H model exhibits a quantum phase diagram similar to that of a spin-1/2 system, including two QSLs around antiferromagnetic and ferromagnetic Kitaev models. For models with finite $A_c$, we map out the quantum phase diagram around two Kitaev points and observe distinct types of in-plane vortex orders developed from these two QSL phases. Interestingly, series of nearly degenerate vortex configurations are discovered in each vortex phases. Using linear spin-wave theory, we demonstrate that these vortex configurations can be understood as a consequence of the quantum correction on a continuous family of degenerate classical states.

Weak ergodicity breaking, particularly through quantum many-body scars (QMBS), has become a significant focus in many-body physics. Krylov state complexity quantifies the spread of quantum states within the Krylov basis and serves as a powerful diagnostic for analyzing nonergodic dynamics. In this work, we study spin-one XXZ magnets and reveal nonergodic behavior tied to QMBS. For the XY model, the nematic Néel state exhibits periodic revivals in Krylov complexity. In the generic XXZ model, we identify spin helix states as weakly ergodicity-breaking states, characterized by low entanglement and nonthermal dynamics. Across different scenarios, the Lanczos coefficients for scarred states display an elliptical pattern, reflecting a hidden SU(2) algebra that enables analytical results for Krylov complexity and fidelity. These findings, which exemplify the rare capability to characterize QMBS analytically, are feasible with current experimental techniques and offer deep insights into the nonergodic dynamics of interacting quantum systems.

We apply the fermionic tensor network (TN) state method to understand the strongly correlated nature in a doped Mott insulator. We conduct a comparative study of the $\sigma t$-$J$ model, in which the no-double-occupancy constraint remains unchanged but the quantum phase string effect associated with doped holes is precisely switched off. Thus, the ground state of the $\sigma t$-$J$ model can serve as a well-controlled reference state of the standard $t$-$J$ model. In the absence of phase string, the spin long-range antiferromagnetic (AFM) order is found to be essentially decoupled from the doped holes, and the latter contribute to a Fermi-liquid-like compressibility and a coherent single-particle propagation with a markedly reduced pairing tendency. In contrast, our TN calculations of the $t$-$J$ model indicate that the AFM order decreases much faster with doping and the single-particle propagation of doped holes gets substantially suppressed, concurrently with a much stronger charge compressibility at small doping and a significantly amplified Cooper pairing tendencies. These findings demonstrate that quantum many-body interference from phase strings plays a pivotal role in the $t$-$J$ model, mediating long-range entanglement between spin and charge degrees of freedom.

We construct a holographic tensor network for the double-scaled SYK model (DSSYK). The moment of the transfer matrix of DSSYK can be mapped to the matrix product state (MPS) of a spin chain. By adding the height direction as a holographic direction, we recast the MPS for DSSYK into the holographic tensor network whose building block is a 4-index tensor with the bond dimension three.

We propose a general strategy to develop quantum many-body approximations of primitives in linear algebra algorithms. As a practical example, we introduce a coupled-cluster inspired framework to produce approximate Hamiltonian moments, and demonstrate its application in various linear algebra algorithms for ground state estimation. Through numerical examples, we illustrate the difference between the ground-state energies arising from quantum many-body linear algebra and those from the analogous many-body perturbation theory. Our results support the general idea of designing quantum many-body approximations outside of perturbation theory, providing a route to new algorithms and approximations.

We introduce a novel numerical scheme for solving the Fokker-Planck equation of discretized Dean-Kawasaki models with a functional tensor network ansatz. The Dean-Kawasaki model describes density fluctuations of interacting particle systems, and it is a highly singular stochastic partial differential equation. By performing a finite-volume discretization of the Dean-Kawasaki model, we derive a stochastic differential equation (SDE). To fully characterize the discretized Dean-Kawasaki model, we solve the associated Fokker-Planck equation of the SDE dynamics. In particular, we use a particle-based approach whereby the solution to the Fokker-Planck equation is obtained by performing a series of density estimation tasks from the simulated trajectories, and we use a functional hierarchical tensor model to represent the density. To address the challenge that the sample trajectories are supported on a simplex, we apply a coordinate transformation from the simplex to a Euclidean space by logarithmic parameterization, after which we apply a sketching-based density estimation procedure on the transformed variables. Our approach is general and can be applied to general density estimation tasks over a simplex. We apply the proposed method successfully to the 1D and 2D Dean-Kawasaki models. Moreover, we show that the proposed approach is highly accurate in the presence of external potential and particle interaction.

The goal of multi-objective optimization is to understand optimal trade-offs between competing objective functions by finding the Pareto front, i.e., the set of all Pareto optimal solutions, where no objective can be improved without degrading another one. Multi-objective optimization can be challenging classically, even if the corresponding single-objective optimization problems are efficiently solvable. Thus, multi-objective optimization represents a compelling problem class to analyze with quantum computers. In this work, we use low-depth Quantum Approximate Optimization Algorithm to approximate the optimal Pareto front of certain multi-objective weighted maximum cut problems. We demonstrate its performance on an IBM Quantum computer, as well as with Matrix Product State numerical simulation, and show its potential to outperform classical approaches.

We show that the braiding of anyons in a quantum spin liquid leaves a distinct dynamical signature in the nonlinear pump-probe response. Using a combination of exact diagonalization and matrix product state techniques, we study the nonlinear pump-probe response of the toric code in a magnetic field, a model which hosts mobile electric $e$ and magnetic $m$ anyonic excitations. While the linear response signal oscillates and decays with time like $\sim t^{-1.3}$, the amplitude of the nonlinear signal for &hcedil;i^{(3)}_{XZZ} features a linear-in-time enhancement at early times. The comparison between &hcedil;i^{(3)}_{XZZ}, which involves the non-trivial braiding of $e$ and $m$ anyons, and &hcedil;i^{(3)}_{XXX} that involves the trivial braiding of the same types of anyons, serves to distinguish the braiding statistics of anyons. We support our analysis by constructing a hard-core anyon model with statistical gauge fields to develop further insights into the time dependence of the pump-probe response. Pump-probe spectroscopy provides a distinctive new probe of quantum spin liquid states, beyond the inconclusive broad features observed in single spin flip inelastic neutron scattering.

Moment polytopes of tensors, the study of which is deeply rooted in invariant theory, representation theory and symplectic geometry, have found relevance in numerous places, from quantum information (entanglement polytopes) and algebraic complexity theory (GCT program and the complexity of matrix multiplication) to optimization (scaling algorithms). Towards an open problem in algebraic complexity theory, we prove separations between the moment polytopes of matrix multiplication tensors and unit tensors. As a consequence, we find that matrix multiplication moment polytopes are not maximal, i.e. are strictly contained in the corresponding Kronecker polytope. As another consequence, we obtain a no-go result for a natural operational characterization of moment polytope inclusion in terms of asymptotic restriction. We generalize the separation and non-maximality to moment polytopes of iterated matrix multiplication tensors. Our result implies that tensor networks where multipartite entanglement structures beyond two-party entanglement are allowed can go beyond projected entangled-pair states (PEPS) in terms of expressivity. Our proof characterizes membership of uniform points in moment polytopes of tensors, and establishes a connection to polynomial multiplication tensors via the minrank of matrix subspaces. As a result of independent interest, we extend these techniques to obtain a new proof of the optimal border subrank bound for matrix multiplication.

Charge creation via quantum tunneling, i.e. dielectric breakdown, is one of the most fundamental and significant phenomena arising from strong light(field)-matter coupling. In this work, we conduct a systematic numerical analysis of quantum tunneling in one-dimensional Mott insulators described by the extended ($U$-$V$) Hubbard model. We discuss the applicability of the analytical formula for doublon-holon (DH) pair production, previously derived for the one-dimensional Hubbard model, which highlights the relationship between the tunneling threshold, the charge gap, and the correlation length. We test the formulas ability to predict both DH pair production and energy increase rate. Using tensor-network-based approaches, we demonstrate that the formula provides accurate predictions in the absence of excitonic states facilitated by the nearest-neighbor interaction $V$. However, when excitonic states emerge, the formula more accurately describes the rate of energy increase than the DH pair creation rate and in both cases gets improved by incorporating the exciton energy as the effective gap.