Search

During the last two decades many new insights into the dynamics of strongly coupled quantum many-body systems have been obtained using gauge/gravity duality, with black holes often playing a universal role. In this white paper we summarize the results obtained and offer some outlook for future developments, including the ongoing mutually beneficial feedback loop with the study of more general, not necessarily holographic, quantum many-body systems., Comment: 18+24 pages

Black holes provide a window into the microscopic structure of spacetime in quantum gravity. Recently the quantum information contained in Hawking radiation has been calculated, verifying a key aspect of the consistency of black hole evaporation with quantum mechanical unitarity. This calculation relied crucially on recent progress in understanding the emergence of bulk spacetime from a boundary holographic description. Spacetime wormholes have played an important role in understanding the underpinnings of this result, and the precision study of such wormholes, in this and other contexts, has been enabled by the development of low-dimensional models of holography. In this white paper we review these developments and describe some of the deep open questions in this subject. These include the nature of the black hole interior, potential applications to quantum cosmology, the gravitational explanation of the fine structure of black holes, and the development of further connections to quantum information and laboratory quantum simulation., Comment: 16 + 17 pages. v2: references added

Modulated symmetries are internal symmetries that act in a non-uniform, spatially modulated way and are generalizations of, for example, dipole symmetries. In this paper, we systematically study the gauging of finite Abelian modulated symmetries in ${1+1}$ dimensions. Working with local Hamiltonians of spin chains, we explore the dual symmetries after gauging and their potential new spatial modulations. We establish sufficient conditions for the existence of an isomorphism between the modulated symmetries and their dual, naturally implemented by lattice reflections. For instance, in systems of prime qudits, translation invariance guarantees this isomorphism. For non-prime qudits, we show using techniques from ring theory that this isomorphism can also exist, although it is not guaranteed by lattice translation symmetry alone. From this isomorphism, we identify new Kramers-Wannier dualities and construct related non-invertible reflection symmetry operators using sequential quantum circuits. Notably, this non-invertible reflection symmetry exists even when the system lacks ordinary reflection symmetry. Throughout the paper, we illustrate these results using various simple toy models., Comment: 71 pages. v2: minor changes + added references

Space-time dependence of imaginary-time propagators, vital for \textit{ab initio} and many-body calculations based on quantum field theories, has been revealed to be compressible using Quantum Tensor Trains (QTTs) [Phys. Rev. X {\bf 13}, 021015 (2023)]. However, the impact of system parameters, like temperature, on data size remains underexplored. This paper provides a comprehensive numerical analysis of the compactness of local imaginary-time propagators in QTT for one-time/-frequency objects and two-time/-frequency objects, considering truncation in terms of the Frobenius and maximum norms. To study worst-case scenarios, we employ random pole models, where the number of poles grows logarithmically with the inverse temperature and coefficients are random. The Green's functions generated by these models are expected to be more difficult to compress than those from physical systems. The numerical analysis reveals that these propagators are highly compressible in QTT, outperforming the state-of-the-art approaches such as intermediate representation and discrete Lehmann representation. For one-time/-frequency objects and two-time/-frequency objects, the bond dimensions saturate at low temperatures, especially for truncation in terms of the Frobenius norm. We provide counting-number arguments for the saturation of bond dimensions for the one-time/-frequency objects, while the origin of this saturation for two-time/-frequency objects remains to be clarified. This paper's findings highlight the critical need for further research on the selection of truncation methods, tolerance levels, and the choice between imaginary-time and imaginary-frequency representations in practical applications., Comment: 23 pages, 11 figures

The entanglement negativity $\mathcal{E}(A:B)$ is a useful measure of quantum entanglement in bipartite mixed states. In random tensor networks (RTNs), which are related to fixed-area states, it was found in [arXiv:2101.11029] that the dominant saddles computing the even R\'enyi negativity $\mathcal{E}^{(2k)}$ generically break the $\mathbb{Z}_{2k}$ replica symmetry. This calls into question previous calculations of holographic negativity using 2D CFT techniques that assumed $\mathbb{Z}_{2k}$ replica symmetry and proposed that the negativity was related to the entanglement wedge cross section. In this paper, we resolve this issue by showing that in general holographic states, the saddles computing $\mathcal{E}^{(2k)}$ indeed break the $\mathbb{Z}_{2k}$ replica symmetry. Our argument involves an identity relating $\mathcal{E}^{(2k)}$ to the $k$-th R\'enyi entropy on subregion $AB^*$ in the doubled state $|{\rho_{AB}}\rangle_{AA^*BB^*}$, from which we see that the $\mathbb{Z}_{2k}$ replica symmetry is broken down to $\mathbb{Z}_{k}$. For $k<1$, which includes the case of $\mathcal{E}(A:B)$ at $k=1/2$, we use a modified cosmic brane proposal to derive a new holographic prescription for $\mathcal{E}^{(2k)}$ and show that it is given by a new saddle with multiple cosmic branes anchored to subregions $A$ and $B$ in the original state. Using our prescription, we reproduce known results for the PSSY model and show that our saddle dominates over previously proposed CFT calculations near $k=1$. Moreover, we argue that the $\mathbb{Z}_{2k}$ symmetric configurations previously proposed are not gravitational saddles, unlike our proposal. Finally, we contrast holographic calculations with those arising from RTNs with non-maximally entangled links, demonstrating that the qualitative form of backreaction in such RTNs is different from that in gravity., Comment: 29 pages, 10 figures

An important result in the theory of quantum control is the "universality" of $2$-local unitary gates, i.e. the fact that any global unitary evolution of a system of $L$ qudits can be implemented by composition of $2$-local unitary gates. Surprisingly, recent results have shown that universality can break down in the presence of symmetries: in general, not all globally symmetric unitaries can be constructed using $k$-local symmetric unitary gates. This also restricts the dynamics that can be implemented by symmetric local Hamiltonians. In this paper, we show that obstructions to universality in such settings can in general be understood in terms of superoperator symmetries associated with unitary evolution by restricted sets of gates. These superoperator symmetries lead to block decompositions of the operator Hilbert space, which dictate the connectivity of operator space, and hence the structure of the dynamical Lie algebra. We demonstrate this explicitly in several examples by systematically deriving the superoperator symmetries from the gate structure using the framework of commutant algebras, which has been used to systematically derive symmetries in other quantum many-body systems. We clearly delineate two different types of non-universality, which stem from different structures of the superoperator symmetries, and discuss its signatures in physical observables. In all, our work establishes a comprehensive framework to explore the universality of unitary circuits and derive physical consequences of its absence., Comment: 20+11 pages, 5 figures

Modern quantum physics is very modular: we first understand basic building blocks (``XXZ Hamiltonian'' ``Jaynes-Cummings'' etc.) and then combine them to explore novel effects. A typical example is placing known systems inside an optical cavity. The Schrieffer-Wolff perturbation method is particularly suited for dealing with these problems, since it casts the perturbation expansion in terms of operator corrections to a Hamiltonian, which is more intuitive than energy level corrections, as in traditional time-independent perturbation theory. However, the method lacks a systematic approach.% and has largely remained a niche topic. In these notes we discuss how \emph{eigenoperator decompositions}, a concept largely used in open quantum systems, can be employed to construct an intuitive and systematic formulation of Schrieffer-Wolff perturbation theory. To illustrate this we revisit various papers in the literature, old and new, and show how they can instead be solved using eigenoperators. Particular emphasis is given to perturbations that couple two systems with very different transition frequencies (highly off-resonance), leading to the so-called dispersive interactions.

Out-of-equilibrium fermionic quantum impurity models (QIM), describing a small interacting system coupled to a continuous fermionic bath, play an important role in condensed matter physics. Solving such models is a computationally demanding task, and a variety of computational approaches are based on finding approximate representations of the bath by a finite number of modes. In this paper, we formulate the problem of finding efficient bath representations as that of approximating a kernel of the bath's Feynman-Vernon influence functional by a sum of complex exponentials, with each term defining a fermionic pseudomode. Under mild assumptions on the analytic properties of the bath spectral density, we provide an analytic construction of pseudomodes, and prove that their number scales polylogarithmically with the maximum evolution time $T$ and the approximation error $\varepsilon$. We then demonstrate that the number of pseudomodes can be significantly reduced by an interpolative matrix decomposition (ID). Furthermore, we present a complementary approach, based on constructing rational approximations of the bath's spectral density using the ``AAA'' algorithm, followed by compression with ID. The combination of two approaches yields a pseudomode count scaling as $N_\text{ID} \sim \log(T)\log(1/\varepsilon)$, and the agreement between the two approches suggests that the result is close to optimal. Finally, to relate our findings to QIM, we derive an explicit Liouvillian that describes the time evolution of the combined impurity-pseudomodes system. These results establish bounds on the computational resources required for solving out-of-equilibrium QIMs, providing an efficient starting point for tensor-network methods for QIMs.

Two quantum theories which look different but are secretly describing the same low-energy physics are said to be dual to each other. When realized in the Topological Holography formalism, duality corresponds to changing the gapped boundary condition on the top boundary of a topological field theory, which determines the symmetry of the system, while not affecting the bottom boundary where all the dynamics take place. In this paper, we show that duality in the Topological Holography formalism can be realized with a Sequential Quantum Circuit applied to the top boundary. As a consequence, the Hamiltonians before and after the duality mapping have exactly the same spectrum in the corresponding symmetry sectors, and the entanglement in the corresponding low-energy eigenstates differs by at most an area law term., Comment: 22 pages, 10 figures, 5 tables

Non-Hermitian quantum dynamics lie in an intermediate regime between unitary Hamiltonian dynamics and trace-preserving non-unitary open quantum system dynamics. Given differences in the noise tolerance of unitary and non-unitary dynamics, it is interesting to consider implementing non-Hermitian dynamics on a noisy quantum computer. In this paper, we do so for a non-Hermitian Ising Floquet model whose many-body dynamics gives rise to persistent temporal oscillations, a form of time crystallinity. In the simplest two qubit case that we consider, there is an infinitely long-lived periodic steady state at certain fine-tuned points. These oscillations remain reasonably long-lived over a range of parameters in the ideal non-Hermitean dynamics and for the levels of noise and imperfection expected of modern day quantum devices. Using a generalized Floquet analysis, we show that infinitely long-lived oscillations are generically lost for arbitrarily weak values of common types of noise and compute corresponding damping rate. We perform simulations using IBM's Qiskit platform to confirm our findings; however, experiments on a real device (ibmq-lima) do not show remnants of these oscillations., Comment: v2

Achieving practical quantum speedup with limited resources is a crucial challenge in both academic and industrial communities. To address this, a partially fault-tolerant quantum computing architecture called ``space-time efficient analog rotation quantum computing architecture (STAR architecture)'' has been recently proposed. This architecture focuses on minimizing resource requirements while maximizing the precision of non-Clifford gates, essential for universal quantum computation. However, non-deterministic processes such as the repeat-until-success (RUS) protocol and state injection can introduce significant computational overhead. Therefore, optimizing the logical circuit to minimize this overhead by using efficient fault-tolerant operations is essential. This paper presents an efficient method for simulating the time evolution of the 2D Hubbard model Hamiltonian, a promising application of the STAR architecture. We present two techniques, parallel injection protocol and adaptive injection region updating, to reduce unnecessary time overhead specific to our architecture. By integrating these with the existing fSWAP technique, we develop an efficient Trotter-based time evolution operation for the 2D Hubbard model. Our analysis reveals an acceleration of over 10 times compared to naive serial compilation. This optimized compilation enables us to estimate the computational resources required for quantum phase estimation of the 2D Hubbard model. For devices with a physical error rate of $p_{\rm phys} = 10^{-4}$, we estimate that approximately $6.5 \times 10^4$ physical qubits are required to achieve faster ground state energy estimation of the $8\times8$ Hubbard model compared to classical computation., Comment: 22 pages, 33 figures

The simulation of quantum many-body systems, relevant for quantum chemistry and condensed matter physics, is one of the most promising applications of near-term quantum computers before fault-tolerance. However, since the vast majority of quantum computing technologies are built around qubits and discrete gate-based operations, the translation of the physical problem into this framework is a crucial step. This translation will often be device specific, and a suboptimal implementation will be punished by the exponential compounding of errors on real devices. The importance of an efficient mapping is already revealed for models of spinful fermions in two or three dimensions, which naturally arise when the relevant physics relates to electrons. Using the most direct and well-known mapping, the Jordan-Wigner transformation, leads to a non-local representation of local degrees of freedom, and necessities efficient decompositions of non-local unitary gates into a sequence of hardware accessible local gates. In this paper, we provide a step-by-step recipe for simulating the paradigmatic two-dimensional Fermi-Hubbard model on a quantum computer using only local operations. To provide the ingredients for such a recipe, we briefly review the plethora of different approaches that have emerged recently but focus on the Derby-Klassen compact fermion mapping in order to make our discussion concrete. We provide a detailed recipe for an end-to-end simulation including embedding on a physical device, preparing initial states such as ground states, simulation of unitary time evolution, and measurement of observables and spectral functions. We explicitly compute the resource requirements for simulating a global quantum quench and conclude by discussing the challenges and future directions for simulating strongly-correlated fermionic matter on quantum computers., Comment: 15 pages, 7 figures (+6 pages appendices)

Chirality is a property of a gapped phase of matter in two spatial dimensions that can be manifested through non-zero thermal or electrical Hall conductance. In this paper, we prove two no-go theorems that forbid such chirality for a quantum state in a finite dimensional local Hilbert space with strict area law entanglement entropies. As a crucial ingredient in the proofs, we introduce a new quantum information-theoretic primitive called instantaneous modular flow, which has many other potential applications., Comment: 5+9 pages, 4 figures

Recently, subdimensional particles including fractons have attracted much attention from various areas. Notable features of such matter phases are mobility constraints and subextensive ground state degeneracies (GSDs). In this paper, we propose a BF-like theory motivated by the Godbillon-Vey invariant, which is a mathematical invariant of the foliated manifold. Our theory hosts subsystem higher form symmetries which manifestly ensure the mobility constraint and subextensive GSD through the spontaneous symmetry breaking. We also discuss some lattice spin models which realize the same low energy behaviours as the BF-like theory. Furthermore, we explore dynamical matter theories which are coupled to the BF-like theory., Comment: 50 pages, 15 figures

Projected Entangled Pair States (PEPS) are a class of quantum many-body states that generalize Matrix Product States for one-dimensional systems to higher dimensions. In recent years, PEPS have advanced understanding of strongly correlated systems, especially in two dimensions, e.g., quantum spin liquids. Typically described by tensor networks on regular lattices (e.g., square, cubic), PEPS have also been adapted for irregular graphs, however, the computational cost becomes prohibitive for dense graphs with large vertex degrees. In this paper, we present a PEPS algorithm to simulate low-energy states and dynamics defined on arbitrary, fluctuating, and densely connected graphs. We introduce a cut-off, $\kappa \in \mathbb{N}$, to constrain the vertex degree of the PEPS to a set but tunable value, which is enforced in the optimization by applying a simple edge-deletion rule, allowing the geometry of the PEPS to change and adapt dynamically to the system's correlation structure. We benchmark our flexible PEPS algorithm with simulations of classical spin glasses and quantum annealing on densely connected graphs with hundreds of spins, and also study the impact of tuning $\kappa$ when simulating a uniform quantum spin model on a regular (square) lattice. Our work opens the way to apply tensor network algorithms to arbitrary, even fluctuating, background geometries., Comment: 4 pages (main text), 7 pages (appendix), 15 figures

In this paper, we study a stochastically driven non-equilibrium quantum system where the driving protocols consist of hopping and waiting processes. The waiting times between two hopping processes satisfy a heavy-tailed distribution. By calculating the squared width of the wavepackets, our findings demonstrate the emergence of various anomalous transport phenomenons when the system remains unchanged within the heavy-tailed regime, including superdiffusive, subdiffusive, and standard diffusive motion. Only subdiffusion occurs when the system has evolved during the waiting process. All these transport behaviors are accompanied by a breakdown of ergodicity, highlighting the complex dynamics induced by the stochastic driving mechanism., Comment: 7 pages, 8 figures

Simulation of quantum many-body systems is a promising application of quantum computers. However, implementing the time-evolution operator as a quantum circuit efficiently on near-term devices with limited resources is challenging. Standard approaches like Trotterization often require deep circuits, making them impractical. This paper proposes a hybrid quantum-classical algorithm called Local Subspace Variational Quantum Compilation (LSVQC) for compiling the time-evolution operator. The LSVQC uses variational optimization to reproduce the action of the target time-evolution operator within a physically reasonable subspace. Optimization is performed on small local subsystems based on the Lieb-Robinson bound, allowing for cost function evaluation using small-scale quantum devices or classical computers. Numerical simulations on a spin-lattice model and an $\mathit{\text{ab initio}}$ effective model of strongly correlated material Sr$_2$CuO$_3$ demonstrate the algorithm's effectiveness. It is shown that the LSVQC achieves a 95% reduction in circuit depth compared to Trotterization while maintaining accuracy. The subspace restriction also reduces resource requirements and improves accuracy. Furthermore, we estimate the gate count needed to execute the quantum simulations using the LSVQC on near-term quantum computing architectures in the noisy intermediate-scale or early fault-tolerant quantum computing era. Our estimation suggests that the acceptable physical gate error rate for the LSVQC can be significantly larger than for Trotterization., Comment: 29 pages, 14 figures

In quantum information science, Clifford operators and stabilizer codes play a central role for systems of qubits (or qudits). In this paper, we study the analogous objects for systems of Majorana fermions. A crucial role is played by fermion parity symmetry, which is an unbreakable symmetry present in any system in which the fundamental degrees of freedom are fermionic. We prove that the subgroup of parity-preserving fermionic Cliffords can be represented by the orthogonal group over the binary field $\mathbb{F}_2$, and we show how it can be generated by braiding operators and used to construct any (even-parity) Majorana stabilizer code. We also analyze the frame potential for this so-called p-Clifford group, proving that it is equivalent to the frame potential of the ordinary Clifford group acting on a fixed-parity sector of the Hilbert space., Comment: 16 pages, small errors fixed, comments welcome

The bulk-boundary correspondence of topological phases suggests strong connections between the topological features in a d+1-dimensional bulk and the potentially gapless theory on the (d-1)+1-dimensional boundary. In 2+1D topological phases, a direct correspondence can exist between anyonic excitations in the bulk and the topological point defects/primary fields in the boundary 1+1D conformal field theory. In this paper, we study how line excitations in 3+1D topological phases become line defects in the boundary 2+1D theory using the Topological Holography/Symmetry Topological Field Theory framework. We emphasize the importance of "descendent" line excitations and demonstrate in particular the effect of the Majorana chain defect: it leads to a distinct loop condensed gapped boundary state of the 3+1D fermionic Z2 topological order, and leaves signatures in the 2+1D Majorana-cone critical theory that describes the transition between the two types of loop condensed boundaries. Effects of non-invertible line excitations, such as Cheshire strings, are also discussed in bosonic 3+1D topological phases and the corresponding 2+1D critical points., Comment: 19 pages, 13 figures. Comments are welcome

In this paper we show how to measure in the setting of digital quantum simulations the reflection and transmission amplitudes of the one-dimensional scattering of a particle with a short-ranged potential. The main feature of the protocol is the coupling between the particle and an ancillary spin-1/2 degree of freedom. This allows us to reconstruct tomographically the scattering amplitudes, which are in general complex numbers, from the readout of one qubit. Applications of our results are discussed., Comment: 9+4 pages, 8+6 figures

The non-equilibrium process where the system does not evolve to the featureless state is one of the new central objects in the non-equilibrium phenomena. In this paper, starting from the short-range entangled state in the two-dimensional conformal field theories ($2$d CFTs), the boundary state with a regularization, we evolve the system with the inhomogeneous Hamiltonians called M\"obius/SSD ones. Regardless of the details of CFTs considered in this paper, during the M\"obius evolution, the entanglement entropy exhibits the periodic motion called quantum revival. During SSD time evolution, except for some subsystems, in the large time regime, entanglement entropy and mutual information are approximated by those for the vacuum state. We argue the time regime for the subsystem to cool down to vacuum one is $t_1 \gg \mathcal{O}(L\sqrt{l_A})$, where $t_1$, $L$, and $l_A$ are time, system, and subsystem sizes. This finding suggests the inhomogeneous quench induced by the SSD Hamiltonian may be used as the preparation for the approximately-vacuum state. We propose the gravity dual of the systems considered in this paper, furthermore, and generalize it. In addition to them, we discuss the relation between the inhomogenous quenches and continuous multi-scale entanglement renormalization ansatz (cMERA)., Comment: 32+4 pages, 11 figures

Some topological lattice models in two spatial dimensions exhibit intricate lattice size dependence in their ground state degeneracy (GSD). This and other features such as the position-dependent anyonic excitations are manifestations of UV/IR mixing. In the first part of the paper we exactly calculate the topological entanglement entropy (TEE) of one such model, the rank-2 toric code, after identifying the minimum entropy states for a non-contractible boundary around the torus. The resulting TEE, as with the GSD, shows intricate dependence on the lattice size. In the latter part of the paper we focus on the fact that the rank-2 toric code is an example of a translation symmetry-enriched topological phase, and show that viewing distinct lattice size as a consequence of different translation symmetry defects can explain both our TEE results and the GSD of the rank-2 toric code. Our work establishes the translation symmetry defect framework as a robust description of the UV/IR mixing in topological lattice models., Comment: 17 pages, 6 figures

Transformations between the plateau states of the quantum Hall effect (QHE) are an archetypical example of quantum phase transitions (QPTs) between phases with non-trivial topological order. These transitions appear to be well-described by the single-particle network theories. The long-standing problem with this approach is that it does not account for Coulomb interactions. In this paper, we show that experimental data in the quantum critical regime for both integer and fractional QHEs can be quantitatively explained by the recently developed phenomenological model of QPTs in interacting systems. This model assumes that all effects of interactions are contained in the life-time of fluctuations as set by the Planckian time $\tau_P=\hbar/k_BT$. The dephasing length is taken as the distance traveled by a non-interacting particle along the bulk edge state over this time. We show that the model also provides quantitative description of QPTs between the ground states of anomalous QHE and axion and Chern insulators. These analyzed systems are connected in that the QPTs occur via quantum percolation. Combining the presented results with the results of two companion papers, we conclude that the Planckian time is the encompassing characteristic of QPTs in interacting systems, independent of dimensionality and microscopic physics., Comment: 6 pages, 3 figures

We carry out a field-theoretical renormalization group procedure based on the Callan-Symanzik equation to calculate the detailed flow for the (multi) two-channel Kondo model and its compactified versions. In doing so, we go beyond the universal terms in the beta function we obtained using poor man's scaling (see arXiv:2308.03590 (companion paper II)) and culminate our analysis of how the compactified versions of the model fare against the original one. Among other results, we explore the large-channel-number limit and extend our considerations to the finite temperature crossover region. Moreover, we gain insights into the contradistinction between the consistent vs. conventional bosonization-debosonization formalisms, thereby advancing our understanding on multiple fronts. In particular, we make use of renormalization-flow arguments to further justify the consistent refermionization of the parallel Kondo interaction we presented earlier (see arXiv:2308.03569 (companion paper I)), Comment: 28 pages, 14 figures, 2 tables

For moir\'e bilayer TMD superlattices, full-configuration-interaction (FCI) calculations are presented that take into account both the intra-moir\'e-quantum-dot (MQD) charge-carrier Coulombic interactions, as well as the crystal-field effect from the surrounding moir\'e pockets (inter-moir\'e-QD interactions). Such FCI calculations enable an effective computational embedding strategy and allow for a complete interpretation of the counterintuitive experimental observations reported recently in the context of moir\'e TMD superlattices at integer fillings $\nu=2$ and 4. Two novel states of matter are reported: (i) a genuinely quantum-mechanical supercrystal of {\it sliding\/} Wigner molecules (WMs) for unstrained moir\'e TMD materials (when the crystal field is commensurate with the trilobal symmetry of the confining potential in each embedded MQD) and (ii) a supercrystal of {\it pinned\/} Wigner molecules when strain is involved and the crystal field is incommensurate with the trilobal symmetry of the confining potential in each embedded MQD. The case of $\nu=3$ is an exception, in that both unstrained and strained cases produce a supercrystal of pinned WMs, which is due to the congruence of intrinsic (that of the WM) and external (that of the confining potential of the MQD) $C_3$ point-group symmetries. Furthermore, it is shown that the unrestricted Hartree-Fock approach fails to describe the supercrystal of sliding WMs in the unstrained case, providing a qualitative agreement only in the case of a supercrystal of pinned WMs, Comment: 7 pages, 4 color figures. For similar papers, see https://sites.gatech.edu/cyannouleas/. arXiv admin note: text overlap with arXiv:2403.12262

The spectral properties of electrons confined in a wire-like quasi-one-dimensional (1D) elongated quantum dot (EQD) coupler between silicon qubits, are investigated with a newly developed valley-augmented unrestricted Hartree-Fock (va-UHF) method, generalized to include the valley degree of freedom treated as an isospin, allowing calculations for a large number of electrons. The lower energy symmetry-broken solutions of the self-consistent generalized Pople-Nesbet equations exhibit, for a confinement that has been modeled after an experimentally fabricated one in silicon, formation of Wigner-molecular polymeric (longitudinal) chains, initiating through charge accumulation at the edges of the finite-length quasi-1D wire. An increasing number of parallel zig-zag chains form as the number of electrons loaded into the confinement is increased, with the formation of newly added chains determined by the strength of the transverse harmonic confinement. The broken-symmetry va-UHF solutions, subsequently augmented by the quantum-mechanically required parity-restoration, go beyond the va-UHF single-determinant solution, predicting formation of entangled Wigner-molecular chains whose charge distributions obliterate the zig-zag organization of the broken-symmetry solutions. The symmetry-restored va-UHF methodology enables systematic investigations of multi-electron complex nano-scale confined structures that could be targeted for future imaging microscopy experiments in silicon and other materials (e.g., 1D domain walls in TMD materials), and quantum information utilization., Comment: 10 pages, 10 color figures. Accepted for publication in Physical Review Applied. For related papers, see https://sites.gatech.edu/cyannouleas/

To simulate thermalizing systems at long times, the most straightforward approach is to calculate the thermal properties at the corresponding energy. In a quantum many-body system of size $N$, for local observables and many initial states, this approach has an error of $O(1/N)$, which is reminiscent of the finite-size error of the equivalence of ensembles. In this paper, we propose a simple and efficient numerical method so that the simulation error is of higher order in $1/N$. This finite-size error scaling is proved by assuming the eigenstate thermalization hypothesis.

The qubits in a D-wave quantum annealer (D-wave QA) are designed on a Pegasus graph that is different from structure of a combinatorial optimization problem. This situation requires embedding with the chains connected by ferromagnetic (FM) coupling $J_c$ between the qubits. Weak and strong $J_c$ values induce chain breaking and enforcement of chain energy, which reduce the accuracy of quantum annealing (QA) measurements, respectively. In addition, we confirmed that even though the D-Wave Ocean package provides a default coupling $J_c^{\text{default}}$, it is not an optimal coupling $J_c^{\text{optimal}}$ that maximizes the possible correct rate of QA measurements. In this paper, we present an algorithm how $J_c^{\text{optimal}}$ with the maximum probability $p$ for observing the possible lowest energy is determined. Finally, we confirm that the extracted $J_c^{\text{optimal}}$ show much better $p$ than $J_c^{\text{default}}$ in QA measurements of various parameters of frustrated and fully connected combinatorial optimization problems. The open code is available in \textit{https://github.com/HunpyoLee/OptimizeChainStrength}., Comment: 5 pages, 6 figures

In this paper, we theoretically study a class of 3D non-liquid states that show exotic boundary phenomena in the thermodynamical limit. More concretely, we focus on a class of 3D fracton topological orders formed via stacking 2D twisted \(\mathbb{Z}_N\) topologically ordered layers along \(z\)-direction. Nearby layers are coupled while maintaining translation symmetry along \(z\) direction. The effective field theory is given by infinite-component Chern-Simons theory, with an integer-valued symmetric block-tridiagonal Toeplitz \(K\)-matrix whose size is thermodynamically large. With open boundary conditions (OBC) along \(z\), certain choice of \(K\)-matrices exhibits exotic boundary ``Toeplitz braiding'', where the mutual braiding phase angle between two anyons at opposite boundaries oscillates and remains non-zero in the thermodynamic limit. In contrast, in trivial case, the mutual braiding phase angle decays exponentially to zero in the thermodynamical limit. As a necessary condition, this phenomenon requires the existence of boundary zero modes in the \(K\)-matrix spectrum under OBC. We categorize nontrivial \(K\)-matrices into two distinct types. Each type-I possesses two boundary zero modes, whereas each type-II possesses only one boundary zero mode. Interestingly, the integer-valued Hamiltonian matrix of the familiar 1D ``Su-Schrieffer-Heeger model'' can be used as a non-trivial $K$-matrix. Importantly, since large-gauge-invariance ensures integer quantized \(K\)-matrix entries, global symmetries are not needed to protect these zero modes. We also present numerical simulation as well as finite size scaling, further confirming the above analytical results. Motivated by the present field-theoretical work, it will be interesting to construct 3D lattice models for demonstrating Toeplitz braiding, which is left to future investigation., Comment: 23p

Rydberg atom array experiments have demonstrated the ability to act as powerful quantum simulators, preparing strongly-correlated phases of matter which are challenging to study for conventional computer simulations. A key direction has been the implementation of interactions on frustrated geometries, in an effort to prepare exotic many-body states such as spin liquids and glasses. In this paper, we apply two-dimensional recurrent neural network (RNN) wave functions to study the ground states of Rydberg atom arrays on the kagome lattice. We implement an annealing scheme to find the RNN variational parameters in regions of the phase diagram where exotic phases may occur, corresponding to rough optimization landscapes. For Rydberg atom array Hamiltonians studied previously on the kagome lattice, our RNN ground states show no evidence of exotic spin liquid or emergent glassy behavior. In the latter case, we argue that the presence of a non-zero Edwards-Anderson order parameter is an artifact of the long autocorrelations times experienced with quantum Monte Carlo simulations. This result emphasizes the utility of autoregressive models, such as RNNs, to explore Rydberg atom array physics on frustrated lattices and beyond., Comment: 13 pages, 5 figures, 3 tables. Link to GitHub repository: https://github.com/mhibatallah/RNNWavefunctions

It is a challenging problem to correctly characterize the symmetry-protected topological (SPT) phases in open quantum systems. As the measurement-based quantum computation (MBQC) utilizes non-trivial edge states of the SPT phases as the logical qubit, its computational power is closely tied to the non-trivial topological nature of the phases. In this paper, we propose to use the gate fidelity which is a measure of the computational power of the MBQC to identify the SPT phases in mixed-state settings. Specifically, we investigate the robustness of the Haldane phase by considering the MBQC on the Affleck-Kennedy-Lieb-Tasaki state subject to different types of noises. To illustrate how our criterion works, we analytically and numerically calculated the gate fidelity to find that its behavior depends crucially on whether the noises satisfy a certain symmetry condition with respect to the on-site $\mathbb{Z}_2 \times \mathbb{Z}_2$ symmetry. In particular, the fidelity for the identity gate, which is given by the sum of the non-local string order parameters, plays an important role. Furthermore, we demonstrate that a stronger symmetry conditions are required to be able to perform other (e.g., the $Z$-rotation gate) gates with high fidelity. By examining which unitary gates can be implemented with the MBQC on the decohered states, we can gain a useful insight into the richer structure of noisy SPT states that cannot be captured solely by the string order parameters., Comment: 17 pages, 8 figures

We address dissipative dynamics of the one-dimensional nearest-neighbour $XX$ spin-$1/2$ chain governed by the Gorini-Kossakowski-Sudarshan-Lindblad (GKSL) equation. In the absence of dissipation the model is integrable. We identify a broad class of dissipative terms that generically destroy integrability but leave the operator space of the model fragmented into an extensive number of dynamically disjoint subspaces of varying dimensions. In sufficiently small subspaces the GKSL equation in the Heisenberg representation can be easily solved, sometimes in a closed analytical form. We provide an example of such an exact solution for a specific choice of dissipative terms. It is found that observables experience the Wannier-Stark localization in the corresponding operator subspace. As a result, the expectation values of the observables are linear combinations of essentially a few discrete decay modes, the long time dynamics being governed by the slowest mode. We examine the complex Liouvillian eigenvalue corresponding to this latter mode as a function of the dissipation strength. We find an exceptional point at a critical dissipation strength that separates oscillating and non-oscillating decay. We also describe a different type of dissipation that leads to a single decay mode in the whole operator subspace. Finally, we point out that our exact solutions of the GKSL equation entail exact solutions of the Schr\"odinger equation describing the quench dynamics in closed spin ladders dual to the dissipative spin chains., Comment: Originally this work appeared as a part of arXiv:2304.03155v1. Later this preprint was split into two separate papers. The first one retained the original arXiv identifier. The second one is presented here

Recently, it was proposed that the chiral central charge of a gapped, two-dimensional quantum many-body system is proportional to a bulk ground state entanglement measure known as the modular commutator. While there is significant evidence to support this relation, we show in this paper that it is not universal. We give examples of lattice systems that have vanishing chiral central charge which nevertheless give nonzero "spurious" values for the modular commutator for arbitrarily large system sizes, in both one and two dimensions. Our examples are based on cluster states and utilize the fact that they can generate nonlocal modular Hamiltonians., Comment: 8 pages, 5 figures

In this paper, we study the twisted gauging on the (1+1)d lattice and construct various non-local mappings on the lattice operators. To be specific, we define the twisted Gauss law operator and implement the twisted gauging of the finite group on the lattice motivated by the orbifolding procedure in the conformal field theory, which involves the data of non-trivial element in the second cohomology group of the gauge group. We show the twisted gauging is equivalent to the two-step procedure of first applying the SPT entangler and then untwisted gauging. We use the twisted gauging to construct the triality (order 3) and $p$-ality (order $p$) mapping on the $\mathbb{Z}_p\times \mathbb{Z}_p$ symmetric Hamiltonians, where $p$ is a prime. Such novel non-local mappings generalize Kramers-Wannier duality and they preserve the locality of symmetric operators but map charged operators to non-local ones. We further construct quantum process to realize these non-local mappings and analyze the induced mappings on the phase diagrams. For theories that are invariant under these non-local mappings, they admit the corresponding non-invertible symmetries. The non-invertible symmetry will constrain the theory at the multicritical point between the gapped phases. We further give the condition when the non-invertible symmetry can have symmetric gapped phase with a unique ground state., Comment: 49 pages, 12 figures, 3 tables. v2 updates references and minor changes. Comments are welcome

It is known that particles with exotic properties can emerge in systems made of simple constituents such as qubits, due to long-range quantum entanglement. In this paper, we provide quantitative characterizations of entanglement necessary for emergent anyons and fermions by using the geometric entanglement measure (GEM) which quantifies the maximal overlap between a given state and any short-range entangled states. For systems with emergent anyons, based on the braiding statistics, we show that the GEM scales linearly in the system size regardless of microscopic details. The phenomenon of emergent anyons can also be understood within the framework of quantum error correction (QEC). Specifically, we show that the GEM of any 2D stabilizer codes must be at least quadratic in the code distance. Our proof is based on a generic prescription for constructing string operators, establishing a rigorous and direct connection between emergent anyons and QEC. For systems with emergent fermions, despite that the ground state subspaces could be exponentially huge and their coding properties could be rather poor, we show that the GEM also scales linearly in the system size. Our results also establish an intriguing link between quantum anomaly and entanglement: a quantum state respecting anomalous $1$-form symmetries, be it pure or mixed, must be long-range entangled and have large GEM, offering a non-trivial class of intrinsically mixed state phases.

Recovering quantum information from a noisy quantum system is one of the central challenges in quantum information science and fundamental physics. The key to this goal is explicitly constructing a decoder. In this paper, we provide two explicit decoding quantum circuits that are both capable of recovering quantum information when a decoupling condition is satisfied, i.e., when quantum information is in principle recoverable. The decoders are constructed by using the fixed-point amplitude amplification algorithm based on the quantum singular value transformation, which significantly extends an approach by Yoshida and Kitaev in a specific noise model to general situations. We also show that the proposed decoding circuits reduce the computational cost compared to a previously known explicit decoder. Our constructions not only show an intriguing intersection between decoders and quantum algorithms but also reveal the power of an algorithmic approach to recovering quantum information., Comment: 23 pages, 12 figures, 2 tables

In one spatial dimension, families of short-range entangled many-body quantum states, parameterized over some parameter space, can be topologically distinguished and classified by topological invariants built from the higher Berry phase -- a many-body generalization of the Berry phase. Previous works identified the underlying mathematical structure (the gerbe structure) and introduced a multi-wavefunction overlap, a generalization of the inner product in quantum mechanics, which allows for the extraction of the higher Berry phase and topological invariants. In this paper, building on these works, we introduce a connection, the higher Berry connection, for a family of parameterized Matrix Product States (MPS) over a parameter space. We demonstrate the use of our formula for simple non-trivial models., Comment: 34 pages, 6 figures

We consider families of invertible many-body quantum states in $d$ spatial dimensions that are parameterized over some parameter space $X$. The space of such families is expected to have topologically distinct sectors classified by the cohomology group $\mathrm{H}^{d+2}(X;\mathbb{Z})$. These topological sectors are distinguished by a topological invariant built from a generalization of the Berry phase, called the higher Berry phase. In the previous work, we introduced a generalized inner product for three one-dimensional many-body quantum states, (``triple inner product''). The higher Berry phase for one-dimensional invertible states can be introduced through the triple inner product and furthermore the topological invariant, which takes its value in $\mathrm{H}^{3}(X;\mathbb{Z})$, can be extracted. In this paper, we introduce an inner product of four two-dimensional invertible quantum many-body states. We use it to measure the topological nontriviality of parameterized families of 2d invertible states. In particular, we define a topological invariant of such families that takes values in $\mathrm{H}^{4}(X;\mathbb{Z})$. Our formalism uses projected entangled pair states (PEPS). We also construct a specific example of non-trivial parameterized families of 2d invertible states parameterized over $\mathbb{R}P^4$ and demonstrate the use of our formula. Applications for symmetry-protected topological phases are also discussed., Comment: 31pages, 10 figures

Symmetry in mixed quantum states can manifest in two distinct forms: \textit{strong symmetry}, where each individual pure state in the quantum ensemble is symmetric with the same charge, and \textit{weak symmetry}, which applies only to the entire ensemble. This paper explores a novel type of spontaneous symmetry breaking (SSB) where a strong symmetry is broken to a weak one. While the SSB of a weak symmetry is measured by the long-ranged two-point correlation function $\mathrm{Tr}(O_xO^{\dagger}_y\rho)$, the strong-to-weak SSB (SW-SSB) is measured by the fidelity $F(\rho, O_xO^{\dagger}_y\rho O_yO^{\dagger}_x)$, dubbed the \textit{fidelity correlator}. We prove that SW-SSB is a universal property of mixed-state quantum phases, in the sense that the phenomenon of SW-SSB is robust against symmetric low-depth local quantum channels. { We also show that the symmetry breaking is "spontaneous " in the sense that the effect of a local symmetry-breaking measurement cannot be recovered locally.} We argue that a thermal state at a nonzero temperature in the canonical ensemble (with fixed symmetry charge) should have spontaneously broken strong symmetry. Additionally, we study non-thermal scenarios where decoherence induces SW-SSB, leading to phase transitions described by classical statistical models with bond randomness. In particular, the SW-SSB transition of a decohered Ising model can be viewed as the "ungauged" version of the celebrated toric code decodability transition. We confirm that, in the decohered Ising model, the SW-SSB transition defined by the fidelity correlator is the only physical transition in terms of channel recoverability. We also comment on other (inequivalent) definitions of SW-SSB, through correlation functions with higher R\'enyi indices., Comment: 17+6 pages, 4 figures

In this paper, we investigate the relationship between entanglement and non-stabilizerness (also known as magic) in matrix product states (MPSs). We study the relation between magic and the bond dimension used to approximate the ground state of a many-body system in two different contexts: full state of magic and mutual magic (the non-stabilizer analogue of mutual information, thus free of boundary effects) of spin-1 anisotropic Heisenberg chains. Our results indicate that obtaining converged results for non-stabilizerness is typically considerably easier than entanglement. For full state magic at critical points and at sufficiently large volumes, we observe convergence with $1/\chi^2$, with $\chi$ being the MPS bond dimension. At small volumes, magic saturation is so quick that, within error bars, we cannot appreciate any finite-$\chi$ correction. Mutual magic also shows a fast convergence with bond dimension, whose specific functional form is however hindered by sampling errors. As a by-product of our study, we show how Pauli-Markov chains (originally formulated to evaluate magic) resets the state of the art in terms of computing mutual information for MPS. We illustrate this last fact by verifying the logarithmic increase of mutual information between connected partitions at critical points. By comparing mutual information and mutual magic, we observe that, for connected partitions, the latter is typically scaling much slower - if at all - with the partition size, while for disconnected partitions, both are constant in size., Comment: 12 pages, 13 figures

Confinement/deconfinement, captivating attributes of high-energy elementary particles, have recently garnered wide attention in quantum simulations based on cold atoms. Yet, the partial confinement, an intermediate state between the confinement and deconfinement, remains underexplored. The partial confinement encapsulates the phenomenon that the confining behavior of charged particles is contingent upon their relative positions. In this paper, we demonstrate that the spin-1 quantum link model provides an excellent platform for exploring partial confinement. We conduct a comprehensive investigation of the physics emerging from partial confinement in both the context of equilibrium and non-equilibrium dynamics. Potential experimental setups using cold atoms are also discussed. Our work offers a simple and feasible routine for the study of confinement-related physics in the state-of-the-art artificial quantum systems subject to gauge symmetries.

In the era of noisy, intermediate-scale quantum (NISQ) devices, the efficient preparation of many-body resource states is a task of paramount importance. In this paper we focus on the deterministic preparation of matrix-product states (MPS) in constant depth by utilizing measurements and classical communication to fuse smaller states into larger ones. We place strong constraints on the MPS that can be prepared using this method, which we refer to as MPS fusion. Namely, we establish that it is necessary for the MPS to have a flat entanglement spectrum. Using the recently introduced split-index MPS (SIMPS) representation, we then introduce a family of states that belong to interesting phases of matter protected by non-onsite symmetries and serve as resources for long-range quantum teleportation, but which lie beyond the scope of ordinary MPS fusion. It is shown constructively that these states can be prepared in constant depth using a broader class of measurement-assisted protocols, which we dub SIMPS fusion. Even in cases when MPS fusion is possible, using SIMPS fusion can give rise to significantly reduced resource overhead. Our results therefore simultaneously establish the boundaries of conventional MPS fusion and push the envelope of which states can be prepared using measurement-assisted protocols., Comment: V2: Updated references

Efficient characterization of higher dimensional many-body physical states presents significant challenges. In this paper, we propose a new class of Project Entangled Pair State (PEPS) that incorporates two isometric conditions. This new class facilitates the efficient calculation of general local observables and certain two-point correlation functions, which have been previously shown to be intractable for general PEPS, or PEPS with only a single isometric constraint. Despite incorporating two isometric conditions, our class preserves the rich physical structure while enhancing the analytical capabilities. It features a large set of tunable parameters, with only a subleading correction compared to that of general PEPS. Furthermore, we analytically demonstrate that this class can encode universal quantum computations and can represent a transition from topological to trivial order., Comment: 5+7 pages, 1 figure

A (2+1)D topological ordered phase with U(1) symmetry may or may not have a symmetric gapped edge state, even if both thermal and electric Hall conductivity are vanishing. It is recently discovered that there are "higher" versions of Hall conductivity valid for fermionic fractional quantum Hall (FQH) states, which obstructs symmetry-preserving gapped edge state beyond thermal and electric Hall conductivity. In this paper, we show that one can extract higher Hall conductivity from a single wave function of an FQH state, by evaluating the expectation value of the "partial rotation" unitary which is a combination of partial spatial rotation and a U(1) phase rotation. This result is verified numerically with the fermionic Laughlin state with $\nu=1/3$, $1/5$, as well as the non-Abelian Moore-Read state. Together with topological entanglement entropy, we prove that the expectation values of the partial rotation completely determines if a bosonic/fermionic Abelian topological order with U(1) symmetry has a symmetry-preserving gappable edge state or not. We also show that thermal and electric Hall conductivity of Abelian topological order can be extracted by partial rotations. Even in non-Abelian FQH states, partial rotation provides the Lieb-Schultz-Mattis type theorem constraining the low-energy spectrum of the bulk-boundary system. The generalization of higher Hall conductivity to the case with Lie group symmetry is also presented., Comment: 17 pages, 4 figures, minor edits

An intriguing feature of type II$_1$ von Neumann algebra is that the entropy of the mixed states is negative. Although the type classification of von Neumann algebra and its consequence in holography have been extensively explored recently, there has not been an explicit calculation of entropy in some physically interesting models with type II$_1$ algebra. In this paper, we study the entanglement entropy $S_n$ of the fixed length state $\{|n\rangle\}$ in Double-Scaled Sachdev-Ye-Kitaev model, which has been recently shown to exhibit type II$_1$ von Neumann algebra. These states furnish an orthogonal basis for 0-particle chord Hilbert space. We systematically study $S_n$ and its R\'enyi generalizations $S_n^{(m)}$ in various limit of DSSYK model, ranging $q\in[0,1]$. We obtain exotic analytical expressions for the scaling behavior of $S_n^{(m)}$ at large $n$ for random matrix theory limit ($q=0$) and SYK$_2$ limit ($q=1$), for the former we observe highly non-flat entanglement spectrum. We then dive into triple scaling limits where the fixed chord number states become the geodesic wormholes with definite length connecting left/right AdS$_2$ boundary in Jackiw-Teitelboim gravity. In semi-classical regime, we match the boundary calculation of entanglement entropy with the dilaton value at the center of geodesic, as a nontrivial check of the Ryu-Takayanagi formula.

In a recent paper (arXiv:2206.05152v4), using the exact diagonalization technique, I calculated the energy and other physical properties (electron density, pair correlation function) of a system of $N\le 7$ two-dimensional electrons at the Landau level filling factor $\nu=1/3$, and showed that the variational many-body wave function proposed for this filling factor by Laughlin is far from the true ground state. In this paper I continue to study exact properties of a small ($N\le 7$) system of two-dimensional electrons lying on the lowest Landau level. I analyze the energies and electron densities of the systems with $N\le 7$ electrons continuously as a function of the magnetic field in the range $1/4\lesssim\nu<1$. The physical mechanisms of the appearance of energy gaps in many-particle electron spectra are elucidated. The results obtained clarify the true nature of the ground and excited states of the considered systems., Comment: 22 pages, 15 figures, 10 tables. Part II of the theory of the fractional quantum Hall effect; for the first part see arXiv:2206.05152v4

Fractonic superfluids are exotic phases of matter in which bosons are subject to mobility constraints, resulting in features beyond those of conventional superfluids. These exotic phases arise from the spontaneous breaking of higher-rank symmetry (HRS) in many-body systems with higher-moment conservation, such as dipoles, quadrupoles, and angular moments. The aim of this paper is to introduce exciting developments on the theory of spontaneous symmetry breaking in such systems, which we refer to as ``many-fracton systems''. More specifically, we introduce exciting progress on general aspects of HRS, minimal model construction, realization of symmetry-breaking ground states, order parameter, off-diagonal long-range order (ODLRO), Noether currents with continuity equations, Gross-Pitaevskii equations, quantum fluctuations, Goldstone modes, specific heat, generalized Mermin-Wagner theorem, critical current, Landau criterion, symmetry defects, and Kosterlitz-Thouless (KT)-like physics, hydrodynamics, and dipolar Hubbard model realization. This paper is concluded with several future directions., Comment: Title changed, references updated. A short review on recent progress on higher rank symmetry breaking, fractonic superfluids, dipole (and other higher moments) conservation, and related topics

Non-invertible symmetries of quantum field theories and many-body systems generalize the concept of symmetries by allowing non-invertible operations in addition to more ordinary invertible ones described by groups. The aim of this paper is to point out that these non-invertible symmetries act on local operators by quantum operations, i.e. completely positive maps between density matrices, which form a natural class of operations containing both unitary evolutions and measurements and play an important role in quantum information theory. This observation will be illustrated by the Kramers--Wannier duality of the one-dimensional quantum Ising chain, which is a prototypical example of non-invertible symmetry operations., Comment: 5 pages, 2 figures

The development of powerful numerical techniques has drastically improved our understanding of quantum matter out of equilibrium. Inspired by recent progress in the area of noisy intermediate-scale quantum devices, this paper highlights hybrid Schr\"odinger-Feynman techniques as an innovative approach to efficiently simulate certain aspects of many-body quantum dynamics on classical computers. To this end, we explore the nonequilibrium dynamics of two large subsystems, which interact sporadically in time, but otherwise evolve independently from each other. We consider subsystems with tunable disorder strength, relevant in the context of many-body localization, where one subsystem can act as a bath for the other. Importantly, studying the full interacting system, we observe that signatures of thermalization are enhanced compared to the reference case of having two independent subsystems. Notably, with the here proposed Schr\"odinger-Feynman method, we are able to simulate the pure-state survival probability in systems significantly larger than accessible by standard sparse-matrix techniques., Comment: 7 pages, 3 figures

Finite depth quantum circuits provide an equivalence relation between gapped phases. Moreover, there can be nontrivial domain walls either within the same gapped phase or between different gapped phases, whose equivalence relations are given by finite depth quantum circuits in one lower dimension. In this paper, we use such unitary equivalence relations to study the fusion of one-dimensional gapped phases. In particular, we use finite depth circuits to fuse two gapped phases, local unitaries to fuse two domain walls, and a combination of both to fuse gapped phases with domain walls. This provides a concrete illustration of some simple aspects of the `higher-category' structure of gapped defects in a higher-dimensional trivial gapped bulk state.