DMRG Preprints

QwaveMPS is an open-source Python library for simulating one-dimensional quantum many-body waveguide systems using matrix product states (MPS). It provides a user-friendly interface for constructing, evolving, and analyzing quantum states and operators, facilitating studies in quantum physics and quantum information with waveguide QED systems. This approach enables efficient, scalable simulations by focusing computational resources on the most relevant parts of the quantum system. Thus, one can study a wide range of complex dynamical interactions, including time-delayed feedback effects in the non-Markovian regime and deeply non-linear systems, at a highly reduced computational cost compared to full Hilbert space approaches, making it both practical and convenient to model a variety of open waveguide-QED systems (in Markovian and non-Markovian regimes), treating quantized atoms and quantized photons on an equal footing.

We classify phases of one-dimensional matrix-product states (MPS) under symmetric circuits augmented with symmetric measurements and feedforward. Building on the framework introduced in Gunn et al., Phys. Rev. B 111, 115110 (2025), we extend the analysis from abelian and class-2 nilpotent groups to all finite nilpotent groups. For any such symmetry group $G$, we construct explicit protocols composed of $G$-symmetric circuits and measurements with feedforward that transform symmetry-protected topological (SPT) states into the trivial phase and vice versa using a finite number of measurement rounds determined by the nilpotency class of $G$. Although these transformations are approximate, we prove that their success probability converges to unity in the thermodynamic limit, establishing asymptotically deterministic equivalence. Consequently, all SPT phases protected by finite nilpotent groups collapse to a single phase once symmetric measurements and feedforward are allowed. We further show that the same holds for non-normal MPS with long-range correlations, including GHZ-type states. The central technical ingredient is a hierarchical structure of irreducible representations of nilpotent groups, which enables a recursive reduction of non-abelian components to abelian ones. Our results demonstrate that symmetric measurements lead to a complete collapse of both symmetry-protected and non-normal MPS phases for all finite nilpotent symmetry groups.

While nonstabilizerness (''magic'') is a key resource for universal quantum computation, its behavior in many-body quantum systems, especially near criticality, remains poorly understood. We develop a spectral transfer-matrix framework for the stabilizer Rényi entropy (SRE) in infinite matrix product states, showing that its spectrum contains universal subleading information. In particular, we identify an SRE correlation length -- distinct from the standard correlation length -- which diverges at continuous phase transitions and governs the spatial response of the SRE to local perturbations. We derive exact SRE expressions for the bond dimension $χ=2$ MPS ''skeleton'' of the cluster-Ising model, and we numerically probe its universal scaling along the $ℤ_2$ critical lines in the phase diagram. These results demonstrate that nonstabilizerness captures signatures of criticality and local perturbations, providing a new lens on the interplay between computational resources and emergent phenomena in quantum many-body systems.

We study tensor network varieties associated with the triangular graph, with a focus on the case where one of the physical dimensions is 2. This allows us to interpret the tensors as pencils of matrices. We provide a complete characterization of these varieties in terms of the Kronecker invariants of pencils. We determine their dimension, identifying the cases for which the dimension is smaller than the expected parameter count. We provide necessary conditions for membership in these varieties, in terms of the geometry of classical determinantal varieties, coincident root loci and plane cubic curves. We address some extensions to arbitrary graphs.

Previous studies of the $S=1/2$ triangular-lattice $J_1$--$J_2$ Heisenberg antiferromagnet have inferred the existence of a non-magnetic ground-state phase for an intermediate range of $J_2$, but disagree concerning whether it is a gapped $ℤ_2$ quantum spin liquid (QSL), a gapless (Dirac) QSL, or a weakly symmetry-broken phase. Using an improved dynamical density-matrix renormalization group method, we investigate the relevant intermediate $J_2$ regime for cylinders with circumferences from 6 to 9. Depending on the initial state and boundary conditions, we find two {\it distinct} variational states. The higher energy state is consistent with a Dirac QSL. In the lower-energy state, both the static and dynamical properties are qualitatively similar to the magnetically ordered state at $J_2=0$, suggestive of either a weakly magnetically ordered non-QSL or a gapped QSL proximate to a continuous transition to such an ordered state.

Quantum circuits have been shown to be a fertile ground for realizing long-range entangled phases of matter. While various quantum double models with non-chiral topological order have been theoretically investigated and experimentally implemented, the realization and characterization of chiral topological phases have remained less explored. Here we show that chiral topological phases in spin systems, i.e., chiral spin liquids, can be prepared in quantum circuits using the variational quantum eigensolver (VQE) framework. On top of the VQE ground state, signatures of the chiral topological order are revealed using the recently proposed tangent space excitation ansatz for quantum circuits. We show that, both topological ground state degeneracy and the chiral edge mode can be faithfully captured by this approach. We demonstrate our approach using the Kitaev honeycomb model, finding excellent agreement of low-energy excitation spectrum on quantum circuits with exact solution in all topological sectors. Further applying this approach to a non-exactly solvable chiral spin liquid model on square lattice, the results suggest this approach works well even when the topological sectors are not exactly known.

We present TensorCircuit-NG, a next-generation quantum software platform designed to bridge the gap between quantum physics, artificial intelligence, and high-performance computing. Moving beyond the scope of traditional circuit simulators, TensorCircuit-NG establishes a unified, tensor-native programming paradigm where quantum circuits, tensor networks, and neural networks fuse into a single, end-to-end differentiable computational graph. Built upon industry-standard machine learning backends (JAX, TensorFlow, PyTorch), the framework introduces comprehensive capabilities for approximate circuit simulation, analog dynamics, fermion Gaussian states, qudit systems, and scalable noise modeling. To tackle the exponential complexity of deep quantum circuits, TensorCircuit-NG implements advanced distributed computing strategies, including automated data parallelism and model-parallel tensor network slicing. We validate these capabilities on GPU clusters, demonstrating a near-linear speedup in distributed variational quantum algorithms. TensorCircuit-NG enables flagship applications, including end-to-end QML for CIFAR-100 computer vision, efficient pipelines from quantum states to neural networks via classical shadows, and differentiable optimization of tensor network states for many-body physics.

Accurately simulating the non-Markovian dynamics of open quantum systems remains a significant challenge. While the recently proposed time-evolving matrix product operator (TEMPO) algorithm based on path integrals successfully circumvents the exponential scaling associated with memory length, its reliance on layer-by-layer tensor contractions and compressions leads to steep scaling with respect to the system Hilbert space dimension. In this work, we introduce the effective Hamiltonian-based TEMPO (EH-TEMPO) algorithm, which reformulates the calculation of the Feynman-Vernon influence functional as an imaginary time evolution governed by an effective Hamiltonian. We demonstrate that this effective Hamiltonian admits a highly compact matrix product operator representation, enabling substantial compression with negligible loss of accuracy. Combining a one-shot global evolution with a backward retrieval approach, EH-TEMPO significantly reduces algorithmic complexity and is naturally suited for GPU acceleration. We benchmark the method against the process tensor TEMPO algorithm using the 7-site Fenna-Matthews-Olson complex model. The results demonstrate that EH-TEMPO achieves numerically exact accuracy with superior efficiency, delivering speedups of up to 17.5x on GPU hardware compared to standard CPU implementations.

Fractionalization of symmetry - exemplified by spin-charge separation in the 1D Hubbard model and fractional charges in the fractional quantum Hall effect - is a typical strongly correlated phenomena in quantum many-body systems. Despite the success in measuring velocity differences, however, it is still quite challenging in probing emergent excitations' anomalous dimensions experimentally. We propose a off-diagonal dissipative response protocol, leveraging dissipative response theory (DRT), to directly detect spin-charge separation. By selectively dissipating spin-$\downarrow$ particles and measuring the spin-$\uparrow$ response, we uncover a universal temporal signature: the off-diagonal response exhibits a crossover from cubic-in-time ($t^3$) growth at short times to linear-in-time ($t$) decay at long times. Crucially, the coefficients &acaron;rkappa^s (short-time) and &acaron;rkappa^l (long-time) encode the distinct anomalous dimensions and velocities of spinons and holons, providing unambiguous evidence of fractionalization. This signal vanishes trivially without spin-charge separation. Our predictions, verified numerically via tDMRG, with microscopic parameters linking with Luttinger parameters by Bethe ansatz, establish off-diagonal dissipative response as a probe of quantum fractionalization in synthetic quantum matter.