Publications

Order by: publication date number of citations

• Lattice T-duality from non-invertible symmetries in quantum spin chains

  Salvatore D. Pace, Arkya Chatterjee, Shu-Heng Shao

  SciPost Phys. 18, 121 (2025) · published 8 April 2025 | Toggle abstract · pdf

  Dualities of quantum field theories are challenging to realize in lattice models of qubits. In this work, we explore one of the simplest dualities, T-duality of the compact boson CFT, and its realization in quantum spin chains. In the special case of the XX model, we uncover an exact lattice T-duality, which is associated with a non-invertible symmetry that exchanges two lattice U(1) symmetries. The latter symmetries flow to the momentum and winding U(1) symmetries with a mixed anomaly in the CFT. However, the charge operators of the two U(1) symmetries do not commute on the lattice and instead generate the Onsager algebra. We discuss how some of the anomalies in the CFT are nonetheless still exactly realized on the lattice and how the lattice U(1) symmetries enforce gaplessness. We further explore lattice deformations preserving both U(1) symmetries and find a rich gapless phase diagram with special $\mathrm{Spin}(2k)_1$ WZW model points and whose phase transitions all have dynamical exponent $z>1$.

• Clock factorized quantum Monte Carlo method for long-range interacting systems

  Zhijie Fan, Chao Zhang, Youjin Deng

  SciPost Phys. Core 8, 036 (2025) · published 8 April 2025 | Toggle abstract · pdf

  Simulating long-range interacting systems is a challenging task due to its computational complexity that the computational effort for each local update is of order $\mathcal{O}(N)$, where $N$ is the size of the system. In this work, we introduce the clock factorized quantum Monte Carlo method, an efficient technique for simulating long-range interacting quantum systems. The method is developed by generalizing the clock Monte Carlo method for classical systems [Phys. Rev. E 99 010105 (2019)] to the path-integral representation of long-range interacting quantum systems, with some specific treatments for quantum cases and a few significant technical improvements in general. We first explain how the clock factorized quantum Monte Carlo method is implemented to reduce the computational overhead from $\mathcal{O}(N)$ to $\mathcal{O}(1)$. In particular, the core ingredients, including the concepts of bound probabilities and bound rejection events, the recursive sampling procedure, and the fast algorithms for sampling an extensive set of discrete and small probabilities, are elaborated. Next, we show how the clock factorized quantum Monte Carlo method can be flexibly implemented in various update strategies, like the Metropolis and worm-type algorithms. Finally, we demonstrate the high efficiency of the clock factorized quantum Monte Carlo algorithms using examples of three typical long-range interacting quantum systems, including the transverse field Ising model with long-range $z$-$z$ interaction, the extended Bose-Hubbard model with long-range density-density interactions, and the XXZ Heisenberg model with long-range spin interactions. We expect that the clock factorized quantum Monte Carlo method would find broad applications in statistical and condensed-matter physics.

• Photon counting statistics in Gaussian bosonic networks

  Kalle Sulo Ukko Kansanen, Pedro Portugal, Christian Flindt, Peter Samuelsson

  SciPost Phys. 18, 116 (2025) · published 1 April 2025 | Toggle abstract · pdf

  The statistics of transmitted photons in microwave cavities play a foundational role in microwave quantum optics and its technological applications. By utilizing quantum mechanical phase-space methods, we here develop a general theory of the photon counting statistics in Gaussian bosonic networks consisting of driven cavities with beamsplitter interactions and two-mode-squeezing. The dynamics of the network can be captured by a Lyapunov equation for the covariance matrix of the cavity fields, which generalizes to a Riccati equation, when counting fields are included. By solving the Riccati equation, we obtain the statistics of emitted and absorbed photons as well as the time-dependent correlations encoded in waiting time distributions and second-order coherence functions. To illustrate our theoretical framework, we first apply it to a simple linear network consisting of two coupled cavities, for which we evaluate the photon cross-correlations and discuss connections between the photon emission statistics and the entanglement between the cavities. We then consider a bosonic circulator consisting of three coupled cavities, for which we investigate how a synthetic flux may affect the direction of the photon flow, similarly to recent experiments. Our general framework paves the way for systematic investigations of the photon counting statistics in Gaussian bosonic networks.

• Neural quantum state study of fracton models

  Marc Machaczek, Lode Pollet, Ke Liu

  SciPost Phys. 18, 112 (2025) · published 28 March 2025 | Toggle abstract · pdf

  Fracton models host unconventional topological orders in three and higher dimensions and provide promising candidates for quantum memory platforms. Understanding their robustness against quantum fluctuations is an important task but also poses great challenges due to the lack of efficient numerical tools. In this work, we establish neural quantum states (NQS) as new tools to study phase transitions in these models. Exact and efficient parametrizations are derived for three prototypical fracton codes — the checkerboard and X-cube model, as well as Haah's code — both in terms of a restricted Boltzmann machine (RBM) and a correlation-enhanced RBM. We then adapt the correlation-enhanced RBM architecture to a perturbed checkerboard model and reveal its strong first-order phase transition between the fracton phase and a trivial field-polarizing phase. To this end, we simulate this highly entangled system on lattices of up to 512 qubits with high accuracy, representing a cutting-edge application of variational neural-network methods. In addition, we reproduce the phase transition of the X-cube model previously obtained with quantum Monte Carlo and high-order series expansion methods. Our work demonstrates the remarkable potential of NQS in studying complicated three-dimensional problems and highlights physics-oriented constructions of NQS architectures.

• Hierarchy of degenerate stationary states in a boundary-driven dipole-conserving spin chain

  Apoorv Srivastava, Shovan Dutta

  SciPost Phys. 18, 111 (2025) · published 27 March 2025 | Toggle abstract · pdf

  Kinetically constrained spin chains serve as a prototype for structured ergodicity breaking in isolated quantum systems. We show that such a system exhibits a hierarchy of degenerate steady states when driven by incoherent pump and loss at the boundary. By tuning the relative pump and loss and how local the constraints are, one can stabilize mixed steady states, noiseless subsystems, and various decoherence-free subspaces, all of which preserve large amounts of information. We also find that a dipole-conserving bulk suppresses current in steady state. These exact results based on the flow in Hilbert space hold regardless of the specific Hamiltonian or drive mechanism. Our findings show that a competition of kinetic constraints and local drives can induce different forms of ergodicity breaking in open systems, which should be accessible in quantum simulators.

• Many-body quantum dynamics with MCTDH-X

  Paolo Molignini, Sunayana Dutta, Elke Fasshauer

  SciPost Phys. Lect. Notes 94 (2025) · published 26 March 2025 | Toggle abstract · pdf

  The lecture notes on "Many-body Quantum Dynamics with MCTDH-X", adapted from the 2023 Heidelberg MCTDH Summer School, provide an in-depth exploration of the Multiconfigurational Time-Dependent Hartree approach for indistinguishable particles. They serve as a comprehensive guide for understanding and utilizing the MCTDH-X software for both bosonic and fermionic systems. The tutorial begins with an introduction to the MCTDH-X software, highlighting its capability to handle various quantum systems, including those with internal degrees of freedom and long-range interactions. The theoretical foundation is then laid out on how to solve the time-dependent and time-independent Schrödinger equations for many-body systems. The workflow section provides practical instructions on setting up and executing simulations using MCTDH-X. Detailed benchmarks against exact solutions are presented, showcasing the accuracy and reliability of the software in ground-state and dynamic simulations. The notes then delve into the dynamics of quantum systems, covering relaxation processes, time evolution, and the analysis of propagation for both bosonic and fermionic particles. The discussion includes the interpretation of various physical quantities such as energy, density distributions, and orbital occupations. Advanced features of MCTDH-X are also explored in the last section, including the calculation of correlation functions and the creation of visualizations through video tutorials. The notes conclude with a Linux/UNIX command cheat sheet, facilitating ease of use for users operating the software on different systems. Overall, these lecture notes provide a valuable resource for researchers and students in the field of quantum dynamics, offering both theoretical insights and practical guidance on the use of MCTDH-X for studying complex many-body systems.

• Quantum many-body simulations with PauliStrings.jl

  Nicolas Loizeau, J. Clayton Peacock, Dries Sels

  SciPost Phys. Codebases 54 (2025) · published 24 March 2025 | Toggle abstract · pdf

  We present the Julia package PauliStrings.jl for quantum many-body simulations, which performs fast operations on the Pauli group by encoding Pauli strings in binary. All of the Pauli string algebra is encoded into low-level logic operations on integers, and is made efficient by various truncation methods which allow for systematic extrapolation of the results. We illustrate the effectiveness of our package by (i) performing Heisenberg time evolution through direct numerical integration and (ii) by constructing a Liouvillian Krylov space. We benchmark the results against tensor network methods, and we find our package performs favorably. In addition, we show that this representation allows for easy encoding of any geometry. We present results for chaotic and integrable spin systems in 1D as well as some examples in 2D. Currently, the main limitations are the inefficiency of representing non-trivial pure states (or other low-rank operators), as well as the need to introduce dissipation to probe long-time dynamics.

  1 citation

• Codebase release 1.5 for PauliStrings.jl

  Nicolas Loizeau, J. Clayton Peacock, Dries Sels

  SciPost Phys. Codebases 54-r1.5 (2025) · published 24 March 2025 | Toggle abstract · src

  We present the Julia package PauliStrings.jl for quantum many-body simulations, which performs fast operations on the Pauli group by encoding Pauli strings in binary. All of the Pauli string algebra is encoded into low-level logic operations on integers, and is made efficient by various truncation methods which allow for systematic extrapolation of the results. We illustrate the effectiveness of our package by (i) performing Heisenberg time evolution through direct numerical integration and (ii) by constructing a Liouvillian Krylov space. We benchmark the results against tensor network methods, and we find our package performs favorably. In addition, we show that this representation allows for easy encoding of any geometry. We present results for chaotic and integrable spin systems in 1D as well as some examples in 2D. Currently, the main limitations are the inefficiency of representing non-trivial pure states (or other low-rank operators), as well as the need to introduce dissipation to probe long-time dynamics.

  1 citation

• Nonlocality of deep thermalization

  Harshank Shrotriya, Wen Wei Ho

  SciPost Phys. 18, 107 (2025) · published 21 March 2025 | Toggle abstract · pdf

  We study the role of global system topology in governing deep thermalization, the relaxation of a local subsystem towards a maximally-entropic, uniform distribution of post-measurement states, upon observing the complementary subsystem in a local basis. Concretely, we focus on a class of (1+1)d systems exhibiting 'maximally-chaotic' dynamics, and consider how the rate of the formation of such a universal wavefunction distribution depends on boundary conditions of the system. We find that deep thermalization is achieved exponentially quickly in the presence of either periodic or open boundary conditions; however, the rate at which this occurs is twice as fast for the former than for the latter. These results are attained analytically using the calculus of integration over unitary groups, and supported by extensive numerical simulations. Our findings highlight the nonlocal nature of deep thermalization, and clearly illustrates that the physics underlying this phenomenon goes beyond that of standard quantum thermalization, which only depends on the net build-up of entanglement between a subsystem and its complement.

• Conformal geometry from entanglement

  Isaac H. Kim, Xiang Li, Ting-Chun Lin, John McGreevy, Bowen Shi

  SciPost Phys. 18, 102 (2025) · published 19 March 2025 | Toggle abstract · pdf

  In a physical system with conformal symmetry, observables depend on cross-ratios, measures of distance invariant under global conformal transformations (conformal geometry for short). We identify a quantum information-theoretic mechanism by which the conformal geometry emerges at the gapless edge of a 2+1D quantum many-body system with a bulk energy gap. We introduce a novel pair of information-theoretic quantities $(\mathfrak{c}_{\textrm{tot}}, \eta)$ that can be defined locally on the edge from the wavefunction of the many-body system, without prior knowledge of any distance measure. We posit that, for a topological groundstate, the quantity $\mathfrak{c}_{\textrm{tot}}$ is stationary under arbitrary variations of the quantum state, and study the logical consequences. We show that stationarity, modulo an entanglement-based assumption about the bulk, implies (i) $\mathfrak{c}_{\textrm{tot}}$ is a non-negative constant that can be interpreted as the total central charge of the edge theory. (ii) $\eta$ is a cross-ratio, obeying the full set of mathematical consistency rules, which further indicates the existence of a distance measure of the edge with global conformal invariance. Thus, the conformal geometry emerges from a simple assumption on groundstate entanglement. We show that stationarity of $\mathfrak{c}_{\textrm{tot}}$ is equivalent to a vector fixed-point equation involving $\eta$, making our assumption locally checkable. We also derive similar results for 1+1D systems under a suitable set of assumptions.