-
Preparing matrix product states via fusion: constraints and extensions
Authors:
David T. Stephen,
Oliver Hart
Abstract:
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 ca…
▽ More
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.
△ Less
Submitted 9 May, 2024; v1 submitted 25 April, 2024;
originally announced April 2024.
-
Non-onsite symmetries and quantum teleportation in split-index matrix product states
Authors:
David T. Stephen
Abstract:
We describe a class of spin chains with new physical and computational properties. On the physical side, the spin chains give examples of symmetry-protected topological phases that are defined by non-onsite symmetries, i.e. symmetries that are not a tensor product of single-site operators. These phases can be detected by string-order parameters, but notably do not exhibit entanglement spectrum deg…
▽ More
We describe a class of spin chains with new physical and computational properties. On the physical side, the spin chains give examples of symmetry-protected topological phases that are defined by non-onsite symmetries, i.e. symmetries that are not a tensor product of single-site operators. These phases can be detected by string-order parameters, but notably do not exhibit entanglement spectrum degeneracy. On the computational side, the spin chains represent a new class of states that can be used to deterministically teleport information across long distances, with the novel property that the necessary classical side processing is a non-linear function of the measurement outcomes. We also give examples of states that can serve as universal resources for measurement-based quantum computation, providing the first examples of such resources without entanglement spectrum degeneracy. The key tool in our analysis is a new kind of tensor network representation which we call split-index matrix product states (SIMPS). We develop the basic formalism of SIMPS, compare them to matrix product states, show how they are better equipped to describe certain kinds of non-onsite symmetries and quantum teleportation, and discuss how they are also well-suited to describing constrained spin chains.
△ Less
Submitted 24 April, 2024;
originally announced April 2024.
-
Fusion of one-dimensional gapped phases and their domain walls
Authors:
David T. Stephen,
Xie Chen
Abstract:
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 gappe…
▽ More
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.
△ Less
Submitted 27 March, 2024;
originally announced March 2024.
-
Subsystem Symmetry Fractionalization and Foliated Field Theory
Authors:
Po-Shen Hsin,
David T. Stephen,
Arpit Dua,
Dominic J. Williamson
Abstract:
Topological quantum matter exhibits a range of exotic phenomena when enriched by subdimensional symmetries. This includes new features beyond those that appear in the conventional setting of global symmetry enrichment. A recently discovered example is a type of subsystem symmetry fractionalization that occurs through a different mechanism to global symmetry fractionalization. In this work we exten…
▽ More
Topological quantum matter exhibits a range of exotic phenomena when enriched by subdimensional symmetries. This includes new features beyond those that appear in the conventional setting of global symmetry enrichment. A recently discovered example is a type of subsystem symmetry fractionalization that occurs through a different mechanism to global symmetry fractionalization. In this work we extend the study of subsystem symmetry fractionalization through new examples derived from the general principle of embedding subsystem symmetry into higher-form symmetry. This leads to new types of symmetry fractionalization that are described by foliation dependent higher-form symmetries. This leads to field theories and lattice models that support previously unseen anomalous subsystem symmetry fractionalization. Our work expands the range of exotic topological physics that is enabled by subsystem symmetry in field theory and on the lattice.
△ Less
Submitted 14 March, 2024;
originally announced March 2024.
-
Playing nonlocal games across a topological phase transition on a quantum computer
Authors:
Oliver Hart,
David T. Stephen,
Dominic J. Williamson,
Michael Foss-Feig,
Rahul Nandkishore
Abstract:
Many-body quantum games provide a natural perspective on phases of matter in quantum hardware, crisply relating the quantum correlations inherent in phases of matter to the securing of quantum advantage at a device-oriented task. In this paper we introduce a family of multiplayer quantum games for which topologically ordered phases of matter are a resource yielding quantum advantage. Unlike previo…
▽ More
Many-body quantum games provide a natural perspective on phases of matter in quantum hardware, crisply relating the quantum correlations inherent in phases of matter to the securing of quantum advantage at a device-oriented task. In this paper we introduce a family of multiplayer quantum games for which topologically ordered phases of matter are a resource yielding quantum advantage. Unlike previous examples, quantum advantage persists away from the exactly solvable point and is robust to arbitrary local perturbations, irrespective of system size. We demonstrate this robustness experimentally on Quantinuum's H1-1 quantum computer by playing the game with a continuous family of randomly deformed toric code states that can be created with constant-depth circuits leveraging mid-circuit measurements and unitary feedback. We are thus able to tune through a topological phase transition - witnessed by the loss of robust quantum advantage - on currently available quantum hardware. This behavior is contrasted with an analogous family of deformed GHZ states, for which arbitrarily weak local perturbations destroy quantum advantage in the thermodynamic limit. Finally, we discuss a topological interpretation of the game, which leads to a natural generalization involving an arbitrary number of players.
△ Less
Submitted 7 March, 2024;
originally announced March 2024.
-
Frobenius light cone and the shift unitary
Authors:
Chao Yin,
Andrew Lucas,
David T. Stephen
Abstract:
We bound the time necessary to implement the shift unitary on a one-dimensional ring, both using local Hamiltonians and those with power-law interactions. This time is constrained by the Frobenius light cone; hence we prove that shift unitaries are parametrically harder to implement than preparing long-range Bell pairs for sufficiently small power-law exponent. We note an intriguing similarity bet…
▽ More
We bound the time necessary to implement the shift unitary on a one-dimensional ring, both using local Hamiltonians and those with power-law interactions. This time is constrained by the Frobenius light cone; hence we prove that shift unitaries are parametrically harder to implement than preparing long-range Bell pairs for sufficiently small power-law exponent. We note an intriguing similarity between the proof of our results, and the hardness of preparing symmetry-protected topological states with symmetry-preserving Hamiltonians.
△ Less
Submitted 12 February, 2024;
originally announced February 2024.
-
Sequential Adiabatic Generation of Chiral Topological States
Authors:
Xie Chen,
Michael Hermele,
David T. Stephen
Abstract:
In previous work, it was shown that non-trivial gapped states can be generated from a product state using a sequential quantum circuit. Explicit circuit constructions were given for a variety of gapped states at exactly solvable fixed points. In this paper, we show that a similar generation procedure can be established for chiral topological states as well, despite the fact that they lack an exact…
▽ More
In previous work, it was shown that non-trivial gapped states can be generated from a product state using a sequential quantum circuit. Explicit circuit constructions were given for a variety of gapped states at exactly solvable fixed points. In this paper, we show that a similar generation procedure can be established for chiral topological states as well, despite the fact that they lack an exactly solvable form. Instead of sequentially applying local unitary gates, we sequentially evolve the Hamiltonian by changing local terms in one subregion and then the next. The Hamiltonian remains gapped throughout the process, giving rise to an adiabatic evolution mapping the ground state from a product state to a chiral topological state. We demonstrate such a sequential adiabatic generation process for free fermion chiral states like the Chern Insulator and the $p+ip$ superconductor. Moreover, we show that coupling a quantum state to a discrete gauge group can be achieved through a sequential quantum circuit, thereby generating interacting chiral topological states from the free fermion ones.
△ Less
Submitted 5 February, 2024;
originally announced February 2024.
-
Quantum teleportation implies symmetry-protected topological order
Authors:
Yifan Hong,
David T. Stephen,
Aaron J. Friedman
Abstract:
We constrain a broad class of teleportation protocols using insights from locality. In the "standard" teleportation protocols we consider, all outcome-dependent unitaries are Pauli operators conditioned on linear functions of the measurement outcomes. We find that all such protocols involve preparing a "resource state" exhibiting symmetry-protected topological (SPT) order with Abelian protecting s…
▽ More
We constrain a broad class of teleportation protocols using insights from locality. In the "standard" teleportation protocols we consider, all outcome-dependent unitaries are Pauli operators conditioned on linear functions of the measurement outcomes. We find that all such protocols involve preparing a "resource state" exhibiting symmetry-protected topological (SPT) order with Abelian protecting symmetry $\mathcal{G}_{k}= (\mathbb{Z}_2 \times \mathbb{Z}_2)^k$. The $k$ logical states are teleported between the edges of the chain by measuring the corresponding $2k$ string order parameters in the bulk and applying outcome-dependent Paulis. Hence, this single class of nontrivial SPT states is both necessary and sufficient for the standard teleportation of $k$ qubits. We illustrate this result with several examples, including the cluster state, variants thereof, and a nonstabilizer hypergraph state.
△ Less
Submitted 5 March, 2024; v1 submitted 18 October, 2023;
originally announced October 2023.
-
Sequential Quantum Circuits as Maps between Gapped Phases
Authors:
Xie Chen,
Arpit Dua,
Michael Hermele,
David T. Stephen,
Nathanan Tantivasadakarn,
Robijn Vanhove,
Jing-Yu Zhao
Abstract:
Finite-depth quantum circuits preserve the long-range entanglement structure in quantum states and map between states within a gapped phase. To map between states of different gapped phases, we can use Sequential Quantum Circuits which apply unitary transformations to local patches, strips, or other sub-regions of a system in a sequential way. The sequential structure of the circuit on the one han…
▽ More
Finite-depth quantum circuits preserve the long-range entanglement structure in quantum states and map between states within a gapped phase. To map between states of different gapped phases, we can use Sequential Quantum Circuits which apply unitary transformations to local patches, strips, or other sub-regions of a system in a sequential way. The sequential structure of the circuit on the one hand preserves entanglement area law and hence the gapped-ness of the quantum states. On the other hand, the circuit has generically a linear depth, hence it is capable of changing the long-range correlation and entanglement of quantum states and the phase they belong to. In this paper, we discuss systematically the definition, basic properties, and prototypical examples of sequential quantum circuits that map product states to GHZ states, symmetry-protected topological states, intrinsic topological states, and fracton states. We discuss the physical interpretation of the power of the circuits through connection to condensation, Kramers-Wannier duality, and the notion of foliation for fracton phases.
△ Less
Submitted 3 July, 2023;
originally announced July 2023.
-
Charting the space of ground states with tensor networks
Authors:
Marvin Qi,
David T. Stephen,
Xueda Wen,
Daniel Spiegel,
Markus J. Pflaum,
Agnès Beaudry,
Michael Hermele
Abstract:
We employ matrix product states (MPS) and tensor networks to study topological properties of the space of ground states of gapped many-body systems. We focus on families of states in one spatial dimension, where each state can be represented as an injective MPS of finite bond dimension. Such states are short-range entangled ground states of gapped local Hamiltonians. To such parametrized families…
▽ More
We employ matrix product states (MPS) and tensor networks to study topological properties of the space of ground states of gapped many-body systems. We focus on families of states in one spatial dimension, where each state can be represented as an injective MPS of finite bond dimension. Such states are short-range entangled ground states of gapped local Hamiltonians. To such parametrized families over $X$ we associate a gerbe, which generalizes the line bundle of ground states in zero-dimensional families (\emph{i.e.} in few-body quantum mechanics). The nontriviality of the gerbe is measured by a class in $H^3(X, \mathbb{Z})$, which is believed to classify one-dimensional parametrized systems. We show that when the gerbe is nontrivial, there is an obstruction to representing the family of ground states with an MPS tensor that is continuous everywhere on $X$. We illustrate our construction with two examples of nontrivial parametrized systems over $X=S^3$ and $X = \mathbb{R} P^2 \times S^1$. Finally, we sketch using tensor network methods how the construction extends to higher dimensional parametrized systems, with an example of a two-dimensional parametrized system that gives rise to a nontrivial 2-gerbe over $X = S^4$.
△ Less
Submitted 12 May, 2023;
originally announced May 2023.
-
Renormalization of Ising cage-net model and generalized foliation
Authors:
Zongyuan Wang,
Xiuqi Ma,
David T. Stephen,
Michael Hermele,
Xie Chen
Abstract:
A large class of type-I fracton models, including the X-cube model, have been found to be fixed points of the foliated renormalization group (RG). The system size of such foliated models can be changed by adding or removing decoupled layers of $2$D topological states and continuous deformation of the Hamiltonian. In this paper, we study a closely related model -- the Ising cage-net model -- and fi…
▽ More
A large class of type-I fracton models, including the X-cube model, have been found to be fixed points of the foliated renormalization group (RG). The system size of such foliated models can be changed by adding or removing decoupled layers of $2$D topological states and continuous deformation of the Hamiltonian. In this paper, we study a closely related model -- the Ising cage-net model -- and find that this model is not foliated in the same sense. In fact, we point out certain unnatural restrictions in the foliated RG, and find that removing these restrictions leads to a generalized foliated RG under which the Ising cage-net model is a fixed point, and which includes the original foliated RG as a special case. The Ising cage-net model thus gives a prototypical example of the generalized foliated RG, and its system size can be changed either by condensing / uncondensing bosonic planon excitations near a 2D plane or through a linear depth quantum circuit in the same plane. We show that these two apparently different RG procedures are closely related, as they lead to the same gapped boundary when implemented in part of a plane. Finally, we briefly discuss the implications for foliated fracton phases, whose universal properties will need to be reexamined in light of the generalized foliated RG.
△ Less
Submitted 26 July, 2023; v1 submitted 30 December, 2022;
originally announced January 2023.
-
Non-local finite-depth circuits for constructing SPT states and quantum cellular automata
Authors:
David T. Stephen,
Arpit Dua,
Ali Lavasani,
Rahul Nandkishore
Abstract:
Whether a given target state can be prepared by starting with a simple product state and acting with a finite-depth quantum circuit is a key question in condensed matter physics and quantum information science. It underpins classifications of topological phases, as well as the understanding of topological quantum codes, and has obvious relevance for device implementations. Traditionally, this ques…
▽ More
Whether a given target state can be prepared by starting with a simple product state and acting with a finite-depth quantum circuit is a key question in condensed matter physics and quantum information science. It underpins classifications of topological phases, as well as the understanding of topological quantum codes, and has obvious relevance for device implementations. Traditionally, this question assumes that the quantum circuit is made up of unitary gates that are geometrically local. Inspired by the advent of noisy intermediate-scale quantum devices, we reconsider this question with $k$-local gates, i.e. gates that act on no more than $k$ degrees of freedom, but are not restricted to be geometrically local. First, we construct explicit finite-depth circuits of symmetric $k$-local gates which create symmetry-protected topological (SPT) states from an initial a product state. Our construction applies both to SPT states protected by global symmetries and subsystem symmetries, but not to those with higher-form symmetries, which we conjecture remain nontrivial. Next, we show how to implement arbitrary translationally invariant quantum cellular automata (QCA) in any dimension using finite-depth circuits of $k$-local gates. These results imply that the topological classifications of SPT phases and QCA both collapse to a single trivial phase in the presence of $k$-local interactions. We furthermore argue that SPT phases are fragile to generic $k$-local symmetric perturbations. We conclude by discussing the implications for other phases, such as fracton phases, and surveying future directions. Our analysis opens a new experimentally motivated conceptual direction examining the stability of phases and the feasibility of state preparation without the assumption of geometric locality.
△ Less
Submitted 11 January, 2024; v1 submitted 13 December, 2022;
originally announced December 2022.
-
Universal measurement-based quantum computation in a one-dimensional architecture enabled by dual-unitary circuits
Authors:
David T. Stephen,
Wen Wei Ho,
Tzu-Chieh Wei,
Robert Raussendorf,
Ruben Verresen
Abstract:
A powerful tool emerging from the study of many-body quantum dynamics is that of dual-unitary circuits, which are unitary even when read `sideways', i.e., along the spatial direction. Here, we show that this provides the ideal framework to understand and expand on the notion of measurement-based quantum computation (MBQC). In particular, applying a dual-unitary circuit to a many-body state followe…
▽ More
A powerful tool emerging from the study of many-body quantum dynamics is that of dual-unitary circuits, which are unitary even when read `sideways', i.e., along the spatial direction. Here, we show that this provides the ideal framework to understand and expand on the notion of measurement-based quantum computation (MBQC). In particular, applying a dual-unitary circuit to a many-body state followed by appropriate measurements effectively implements quantum computation in the spatial direction. We show how the dual-unitary dynamics generated by the dynamics of the paradigmatic one-dimensional kicked Ising chain with certain parameter choices generate resource states for universal deterministic MBQC. Specifically, after $k$ time-steps, equivalent to a depth-$k$ quantum circuit, we obtain a resource state for universal MBQC on $\sim 3k/4$ encoded qubits. Our protocol allows generic quantum circuits to be `rotated' in space-time and gives new ways to exchange between resources like qubit number and coherence time in quantum computers. Beyond the practical advantages, we also interpret the dual-unitary evolution as generating an infinite sequence of new symmetry-protected topological phases with spatially modulated symmetries, which gives a vast generalization of the well-studied one-dimensional cluster state and shows that our protocol is robust to symmetry-respecting deformations.
△ Less
Submitted 16 May, 2024; v1 submitted 13 September, 2022;
originally announced September 2022.
-
Ergodicity breaking provably robust to arbitrary perturbations
Authors:
David T. Stephen,
Oliver Hart,
Rahul M. Nandkishore
Abstract:
We present a new route to ergodicity breaking via Hilbert space fragmentation that displays an unprecedented level of robustness. Our construction relies on a single emergent (prethermal) conservation law. In the limit when the conservation law is exact, we prove the emergence of Hilbert space fragmentation with an exponential number of frozen configurations. We further prove that every frozen con…
▽ More
We present a new route to ergodicity breaking via Hilbert space fragmentation that displays an unprecedented level of robustness. Our construction relies on a single emergent (prethermal) conservation law. In the limit when the conservation law is exact, we prove the emergence of Hilbert space fragmentation with an exponential number of frozen configurations. We further prove that every frozen configuration is absolutely stable to arbitrary perturbations, to all finite orders in perturbation theory. In particular, our proof is not limited to symmetric perturbations, or to perturbations with compact support, but also applies to perturbations with long-range tails, and even to arbitrary geometrically nonlocal $k$-body perturbations, as long as $k/L \rightarrow 0$ in the thermodynamic limit, where $L$ is linear system size. Additionally, we identify one-form $U(1)$ charges characterizing some non-frozen sectors, and discuss the dynamics starting from typical initial conditions, which we argue is best interpreted in terms of the magnetohydrodynamics of the emergent one-form symmetry.
△ Less
Submitted 5 October, 2022; v1 submitted 8 September, 2022;
originally announced September 2022.
-
Fractionalization of subsystem symmetries in two dimensions
Authors:
David T. Stephen,
Arpit Dua,
José Garre-Rubio,
Dominic J. Williamson,
Michael Hermele
Abstract:
The fractionalization of global symmetry charges is a striking hallmark of topological quantum order. Here, we discuss the fractionalization of subsystem symmetries in two-dimensional topological phases. In line with previous no-go arguments, we show that subsystem symmetry fractionalization is not possible in many cases due to the additional rigid geometric structure of the symmetries. However, w…
▽ More
The fractionalization of global symmetry charges is a striking hallmark of topological quantum order. Here, we discuss the fractionalization of subsystem symmetries in two-dimensional topological phases. In line with previous no-go arguments, we show that subsystem symmetry fractionalization is not possible in many cases due to the additional rigid geometric structure of the symmetries. However, we identify a new mechanism that allows fractionalization, involving global relations between macroscopically many symmetry generators. We find that anyons can fractionalize such relations, meaning that the total charge carried under all generators involved in the global relation is non-trivial, despite the fact that these generators multiply to the identity. We first discuss the general algebraic framework needed to characterize this new type of fractionalization, and then explore this framework using a number of exactly solvable models with $\mathbb{Z}_2$ topological order, including models having line and fractal symmetries. These models all showcase another necessary property of subsystem symmetry fractionalization: fractionalized anyons must have restricted mobility when the symmetry is enforced, such that they are confined to a single line or point in the case of line and fractal symmetries, respectively. Looking forward, we expect that our identification of the importance of global relations in fractionalization will hold significance for the classification of phases with subsystem symmetries in all dimensions.
△ Less
Submitted 24 March, 2022;
originally announced March 2022.
-
String order parameters for symmetry fractionalization in an enriched toric code
Authors:
José Garre-Rubio,
Mohsin Iqbal,
David T. Stephen
Abstract:
We study a simple model of symmetry-enriched topological order obtained by decorating a toric code model with lower-dimensional symmetry-protected topological states. We show that the symmetry fractionalization in this model can be characterized by string order parameters, and that these signatures are robust under the effects of external fields and interactions, up to the phase transition point.…
▽ More
We study a simple model of symmetry-enriched topological order obtained by decorating a toric code model with lower-dimensional symmetry-protected topological states. We show that the symmetry fractionalization in this model can be characterized by string order parameters, and that these signatures are robust under the effects of external fields and interactions, up to the phase transition point. This extends the recent proposal of [New Journal of Physics 21, 113016 (2019)] beyond the setting of fixed-point tensor network states, and solidifies string order parameters as a useful tool to characterize and detect symmetry fractionalization. In addition to this, we observe how the condensation of an anyon that fractionalizes a symmetry forces that symmetry to spontaneously break, and we give a proof of this in the framework of projected entangled pair states. This phenomenon leads to a notable change in the phase diagram of the toric code in parallel magnetic fields.
△ Less
Submitted 4 March, 2021; v1 submitted 5 November, 2020;
originally announced November 2020.
-
Subsystem symmetry enriched topological order in three dimensions
Authors:
David T. Stephen,
José Garre-Rubio,
Arpit Dua,
Dominic J. Williamson
Abstract:
We introduce a model of three-dimensional (3D) topological order enriched by planar subsystem symmetries. The model is constructed starting from the 3D toric code, whose ground state can be viewed as an equal-weight superposition of two-dimensional (2D) membrane coverings. We then decorate those membranes with 2D cluster states possessing symmetry-protected topological order under line-like subsys…
▽ More
We introduce a model of three-dimensional (3D) topological order enriched by planar subsystem symmetries. The model is constructed starting from the 3D toric code, whose ground state can be viewed as an equal-weight superposition of two-dimensional (2D) membrane coverings. We then decorate those membranes with 2D cluster states possessing symmetry-protected topological order under line-like subsystem symmetries. This endows the decorated model with planar subsystem symmetries under which the loop-like excitations of the toric code fractionalize, resulting in an extensive degeneracy per unit length of the excitation. We also show that the value of the topological entanglement entropy is larger than that of the toric code for certain bipartitions due to the subsystem symmetry enrichment. Our model can be obtained by gauging the global symmetry of a short-range entangled model which has symmetry-protected topological order coming from an interplay of global and subsystem symmetries. We study the non-trivial action of the symmetries on boundary of this model, uncovering a mixed boundary anomaly between global and subsystem symmetries. To further study this interplay, we consider gauging several different subgroups of the total symmetry. The resulting network of models, which includes models with fracton topological order, showcases more of the possible types of subsystem symmetry enrichment that can occur in 3D.
△ Less
Submitted 8 September, 2020; v1 submitted 8 April, 2020;
originally announced April 2020.
-
Inaccessible entanglement in symmetry protected topological phases
Authors:
Caroline de Groot,
David T. Stephen,
Andras Molnar,
Norbert Schuch
Abstract:
We study the entanglement structure of symmetry-protected topological (SPT) phases from an operational point of view by considering entanglement distillation in the presence of symmetries. We demonstrate that non-trivial SPT phases in one-dimension necessarily contain some entanglement which is inaccessible if the symmetry is enforced. More precisely, we consider the setting of local operations an…
▽ More
We study the entanglement structure of symmetry-protected topological (SPT) phases from an operational point of view by considering entanglement distillation in the presence of symmetries. We demonstrate that non-trivial SPT phases in one-dimension necessarily contain some entanglement which is inaccessible if the symmetry is enforced. More precisely, we consider the setting of local operations and classical communication (LOCC) where the local operations commute with a global onsite symmetry group $G$, which we call $G$-LOCC, and we define the inaccessible entanglement $E_{inacc}$ as the entanglement that cannot be used for distillation under $G$-LOCC. We derive a tight bound on $E_{inacc}$ which demonstrates a direct relation between inaccessible entanglement and the SPT phase, namely $\log(D_ω^2) \leq E_{inacc} \leq \log(|G|)$, where $D_ω$ is the topologically protected edge mode degeneracy of the SPT phase $ω$ with symmetry $G$. For particular phases such as the Haldane phase, $D_ω= \sqrt{|G|}$ so the bound becomes an equality. We numerically investigate the distribution of states throughout the bound, and show that typically the region near the upper bound is highly populated, and also determine the nature of those states lying on the upper and lower bounds. We then discuss the relation of $E_{inacc}$ to string order parameters, and also the extent to which it can be used to distinguish different SPT phases of matter.
△ Less
Submitted 15 March, 2020;
originally announced March 2020.
-
Detecting subsystem symmetry protected topological order via entanglement entropy
Authors:
David T. Stephen,
Henrik Dreyer,
Mohsin Iqbal,
Norbert Schuch
Abstract:
Subsystem symmetry protected topological (SSPT) order is a type of quantum order that is protected by symmetries acting on lower-dimensional subsystems of the entire system. In this paper, we show how SSPT order can be characterized and detected by a constant correction to the entanglement area law, similar to the topological entanglement entropy. Focusing on the paradigmatic two-dimensional clust…
▽ More
Subsystem symmetry protected topological (SSPT) order is a type of quantum order that is protected by symmetries acting on lower-dimensional subsystems of the entire system. In this paper, we show how SSPT order can be characterized and detected by a constant correction to the entanglement area law, similar to the topological entanglement entropy. Focusing on the paradigmatic two-dimensional cluster phase as an example, we use tensor network methods to give an analytic argument that almost all states in the phase exhibit the same correction to the area law, such that this correction may be used to reliably detect the SSPT order of the cluster phase. Based on this idea, we formulate a numerical method that uses tensor networks to extract this correction from ground-state wave functions. We use this method to study the fate of the SSPT order of the cluster state under various external fields and interactions, and find that the correction persists unless a phase transition is crossed, or the subsystem symmetry is explicitly broken. Surprisingly, these results uncover that the SSPT order of the cluster state persists beyond the cluster phase, thanks to a new type of subsystem time-reversal symmetry. Finally, we discuss the correction to the area law found in three-dimensional cluster states on different lattices, indicating rich behavior for general subsystem symmetries
△ Less
Submitted 5 September, 2019; v1 submitted 20 April, 2019;
originally announced April 2019.
-
Subsystem symmetries, quantum cellular automata, and computational phases of quantum matter
Authors:
David T. Stephen,
Hendrik Poulsen Nautrup,
Juani Bermejo-Vega,
Jens Eisert,
Robert Raussendorf
Abstract:
Quantum phases of matter are resources for notions of quantum computation. In this work, we establish a new link between concepts of quantum information theory and condensed matter physics by presenting a unified understanding of symmetry-protected topological (SPT) order protected by subsystem symmetries and its relation to measurement-based quantum computation (MBQC). The key unifying ingredient…
▽ More
Quantum phases of matter are resources for notions of quantum computation. In this work, we establish a new link between concepts of quantum information theory and condensed matter physics by presenting a unified understanding of symmetry-protected topological (SPT) order protected by subsystem symmetries and its relation to measurement-based quantum computation (MBQC). The key unifying ingredient is the concept of quantum cellular automata (QCA) which we use to define subsystem symmetries acting on rigid lower-dimensional lines or fractals on a 2D lattice. Notably, both types of symmetries are treated equivalently in our framework. We show that states within a non-trivial SPT phase protected by these symmetries are indicated by the presence of the same QCA in a tensor network representation of the state, thereby characterizing the structure of entanglement that is uniformly present throughout these phases. By also formulating schemes of MBQC based on these QCA, we are able to prove that most of the phases we construct are computationally universal phases of matter, in which every state is a resource for universal MBQC. Interestingly, our approach allows us to construct computational phases which have practical advantages over previous examples, including a computational speedup. The significance of the approach stems from constructing novel computationally universal phases of matter and showcasing the power of tensor networks and quantum information theory in classifying subsystem SPT order.
△ Less
Submitted 9 May, 2019; v1 submitted 22 June, 2018;
originally announced June 2018.
-
A computationally universal phase of quantum matter
Authors:
Robert Raussendorf,
Cihan Okay,
Dong-Sheng Wang,
David T. Stephen,
Hendrik Poulsen Nautrup
Abstract:
We provide the first example of a symmetry protected quantum phase that has universal computational power. Throughout this phase, which lives in spatial dimension two, the ground state is a universal resource for measurement based quantum computation.
We provide the first example of a symmetry protected quantum phase that has universal computational power. Throughout this phase, which lives in spatial dimension two, the ground state is a universal resource for measurement based quantum computation.
△ Less
Submitted 4 March, 2019; v1 submitted 28 February, 2018;
originally announced March 2018.
-
Computational Power of Symmetry-Protected Topological Phases
Authors:
David T. Stephen,
Dong-Sheng Wang,
Abhishodh Prakash,
Tzu-Chieh Wei,
Robert Raussendorf
Abstract:
We consider ground states of quantum spin chains with symmetry-protected topological (SPT) order as resources for measurement-based quantum computation (MBQC). We show that, for a wide range of SPT phases, the computational power of ground states is uniform throughout each phase. This computational power, defined as the Lie group of executable gates in MBQC, is determined by the same algebraic inf…
▽ More
We consider ground states of quantum spin chains with symmetry-protected topological (SPT) order as resources for measurement-based quantum computation (MBQC). We show that, for a wide range of SPT phases, the computational power of ground states is uniform throughout each phase. This computational power, defined as the Lie group of executable gates in MBQC, is determined by the same algebraic information that labels the SPT phase itself. We prove that these Lie groups always contain a full set of single-qubit gates, thereby affirming the long-standing conjecture that general SPT phases can serve as computationally useful phases of matter.
△ Less
Submitted 6 July, 2017; v1 submitted 23 November, 2016;
originally announced November 2016.
-
On the Circuit Diameter Conjecture
Authors:
Steffen Borgwardt,
Tamon Stephen,
Timothy Yusun
Abstract:
From the point of view of optimization, a critical issue is relating the combinatorial diameter of a polyhedron to its number of facets $f$ and dimension $d$. In the seminal paper of Klee and Walkup [KW67], the Hirsch conjecture of an upper bound of $f-d$ was shown to be equivalent to several seemingly simpler statements, and was disproved for unbounded polyhedra through the construction of a part…
▽ More
From the point of view of optimization, a critical issue is relating the combinatorial diameter of a polyhedron to its number of facets $f$ and dimension $d$. In the seminal paper of Klee and Walkup [KW67], the Hirsch conjecture of an upper bound of $f-d$ was shown to be equivalent to several seemingly simpler statements, and was disproved for unbounded polyhedra through the construction of a particular 4-dimensional polyhedron $U_4$ with 8 facets. The Hirsch bound for bounded polyhedra was only recently disproved by Santos.
We consider analogous properties for a variant of the combinatorial diameter called the circuit diameter. In this variant, the walks are built from the circuit directions of the polyhedron, which are the minimal non-trivial solutions to the system defining the polyhedron.
We are able to prove that circuit variants of the so-called non-revisiting conjecture and $d$-step conjecture both imply the circuit analogue of the Hirsch conjecture. For the equivalences in [KW67], the wedge construction was a fundamental proof technique. We exhibit why it is not available in the circuit setting, and what are the implications of losing it as a tool.
Further, we show the circuit analogue of the non-revisiting conjecture implies a linear bound on the circuit diameter of all unbounded polyhedra - in contrast to what is known for the combinatorial diameter. Finally, we give two proofs of a circuit version of the $4$-step conjecture. These results offer some hope that the circuit version of the Hirsch conjecture may hold in general.
A challenge in the circuit setting is that different realizations of polyhedra of the same combinatorial structure may have different diameters. We adapt the notion of simplicity to work with circuits in the form of C-simple and wedge-simple polyhedra. We show that it suffices to consider such polyhedra.
△ Less
Submitted 17 April, 2018; v1 submitted 23 November, 2016;
originally announced November 2016.
-
Symmetry-protected topological phases with uniform computational power in one dimension
Authors:
Robert Raussendorf,
Dongsheng Wang,
Abhishodh Prakash,
Tzu-Chieh Wei,
David Stephen
Abstract:
We investigate the usefulness of ground states of quantum spin chains with symmetry-protected topological order (SPTO) for measurement-based quantum computation. We show that, in spatial dimension one, if an SPTO phase supports quantum wire, then, subject to an additional symmetry condition that is satisfied in all cases so far investigated, it can also be used for quantum computation.
We investigate the usefulness of ground states of quantum spin chains with symmetry-protected topological order (SPTO) for measurement-based quantum computation. We show that, in spatial dimension one, if an SPTO phase supports quantum wire, then, subject to an additional symmetry condition that is satisfied in all cases so far investigated, it can also be used for quantum computation.
△ Less
Submitted 23 September, 2016;
originally announced September 2016.
-
Qudit quantum computation on matrix product states with global symmetry
Authors:
Dong-Sheng Wang,
David T. Stephen,
Robert Raussendorf
Abstract:
Resource states that contain nontrivial symmetry-protected topological order are identified for universal single-qudit measurement-based quantum computation. Our resource states fall into two classes: one as the qudit generalizations of the 1D qubit cluster state, and the other as the higher-symmetry generalizations of the spin-1 Affleck-Kennedy-Lieb-Tasaki (AKLT) state, namely, with unitary, orth…
▽ More
Resource states that contain nontrivial symmetry-protected topological order are identified for universal single-qudit measurement-based quantum computation. Our resource states fall into two classes: one as the qudit generalizations of the 1D qubit cluster state, and the other as the higher-symmetry generalizations of the spin-1 Affleck-Kennedy-Lieb-Tasaki (AKLT) state, namely, with unitary, orthogonal, or symplectic symmetry. The symmetry in cluster states protects information propagation (identity gate), while the higher symmetry in AKLT-type states enables nontrivial gate computation. This work demonstrates a close connection between measurement-based quantum computation and symmetry-protected topological order.
△ Less
Submitted 15 March, 2017; v1 submitted 22 September, 2016;
originally announced September 2016.
-
Algorithms for Colourful Simplicial Depth and Medians in the Plane
Authors:
Olga Zasenko,
Tamon Stephen
Abstract:
The colourful simplicial depth of a point x in the plane relative to a configuration of n points in k colour classes is exactly the number of closed simplices (triangles) with vertices from 3 different colour classes that contain x in their convex hull. We consider the problems of efficiently computing the colourful simplicial depth of a point x, and of finding a point, called a median, that maxim…
▽ More
The colourful simplicial depth of a point x in the plane relative to a configuration of n points in k colour classes is exactly the number of closed simplices (triangles) with vertices from 3 different colour classes that contain x in their convex hull. We consider the problems of efficiently computing the colourful simplicial depth of a point x, and of finding a point, called a median, that maximizes colourful simplicial depth.
For computing the colourful simplicial depth of x, our algorithm runs in time O(n log(n) + k n) in general, and O(kn) if the points are sorted around x. For finding the colourful median, we get a time of O(n^4). For comparison, the running times of the best known algorithm for the monochrome version of these problems are O(n log(n)) in general, improving to O(n) if the points are sorted around x for monochrome depth, and O(n^4) for finding a monochrome median.
△ Less
Submitted 25 August, 2016;
originally announced August 2016.
-
A Linear Time Algorithm for the $3$-neighbour Traveling Salesman Problem on Halin graphs and extensions
Authors:
Brad Woods,
Abraham Punnen,
Tamon Stephen
Abstract:
The Quadratic Travelling Salesman Problem (QTSP) is to find a least cost Hamilton cycle in an edge-weighted graph, where costs are defined on all pairs of edges contained in the Hamilton cycle. This is a more general version than the commonly studied QTSP which only considers pairs of adjacent edges. We define a restricted version of QTSP, the $k$-neighbour TSP (TSP($k$)), and give a linear time a…
▽ More
The Quadratic Travelling Salesman Problem (QTSP) is to find a least cost Hamilton cycle in an edge-weighted graph, where costs are defined on all pairs of edges contained in the Hamilton cycle. This is a more general version than the commonly studied QTSP which only considers pairs of adjacent edges. We define a restricted version of QTSP, the $k$-neighbour TSP (TSP($k$)), and give a linear time algorithm to solve TSP($k$) on a Halin graph for $k\leq 3$. This algorithm can be extended to solve TSP($k$) on any fully reducible class of graphs for any fixed $k$ in polynomial time. This result generalizes corresponding results for the standard TSP. TSP($k$) can be used to model various machine scheduling problems as well as an optimal routing problem for unmanned aerial vehicles (UAVs).
△ Less
Submitted 4 September, 2017; v1 submitted 8 April, 2015;
originally announced April 2015.
-
Circuit diameter and Klee-Walkup constructions
Authors:
Tamon Stephen,
Timothy Yusun
Abstract:
Consider a variant of the graph diameter of a polyhedron where each step in a walk between two vertices travels maximally in a circuit direction instead of along incident edges. Here circuit directions are non-trivial solutions to minimally-dependent subsystems of the presentation of the polyhedron. These can be understood as the set of all possible edge directions, including edges that may arise…
▽ More
Consider a variant of the graph diameter of a polyhedron where each step in a walk between two vertices travels maximally in a circuit direction instead of along incident edges. Here circuit directions are non-trivial solutions to minimally-dependent subsystems of the presentation of the polyhedron. These can be understood as the set of all possible edge directions, including edges that may arise from translation of the facets.
It is appealing to consider a circuit analogue of the Hirsch conjecture for graph diameter, as suggested by Borgwardt et al. [BFH15]. They ask whether the known counterexamples to the Hirsch conjecture give rise to counterexamples for this relaxed notion of circuit diameter. We show that the most basic counterexample to the unbounded Hirsch conjecture, the Klee-Walkup polyhedron, does have a circuit diameter that satisfies the Hirsch bound, regardless of representation. We also examine the circuit diameter of the bounded Klee-Walkup polytope.
△ Less
Submitted 17 March, 2015;
originally announced March 2015.
-
On the Numerical Treatment and Dependence of Thermohaline Mixing in Red Giants
Authors:
John Lattanzio,
Lionel Siess,
Ross Church,
George Angelou,
Richard Stancliffe,
Carolyn Doherty,
Thomas Stephen,
Simon Campbell
Abstract:
In recent years much interest has been shown in the process of thermohaline mixing in red giants. In low and intermediate mass stars this mechanism first activates at the position of the bump in the luminosity function, and has been identified as a likely candidate for driving the slow mixing inferred to occur in these stars. One particularly important consequence of this process, which is driven…
▽ More
In recent years much interest has been shown in the process of thermohaline mixing in red giants. In low and intermediate mass stars this mechanism first activates at the position of the bump in the luminosity function, and has been identified as a likely candidate for driving the slow mixing inferred to occur in these stars. One particularly important consequence of this process, which is driven by a molecular weight inversion, is the destruction of lithium. We show that the degree of lithium destruction, or in some cases production, is extremely sensitive to the numerical details of the stellar models. Within the standard 1D diffusion approximation to thermohaline mixing, we find that dfferent evolution codes, with their default numerical schemes, can produce lithium abundances that differ from one another by many orders of magnitude. This disagreement is worse for faster mixing. We perform experiments with four independent stellar evolution codes, and derive conditions for the spatial and temporal resolution required for a converged numerical solution. The results are extremely sensitive to the timesteps used. We find that predicted lithium abundances published in the literature until now should be treated with caution.
△ Less
Submitted 2 November, 2014; v1 submitted 23 October, 2014;
originally announced October 2014.
-
Computational Lower Bounds for Colourful Simplicial Depth
Authors:
Antoine Deza,
Tamon Stephen,
Feng Xie
Abstract:
The colourful simplicial depth problem in dimension d is to find a configuration of (d+1) sets of (d+1) points such that the origin is contained in the convex hull of each set (colour) but contained in a minimal number of colourful simplices generated by taking one point from each set. A construction attaining d^2+1 simplices is known, and is conjectured to be minimal. This has been confirmed up t…
▽ More
The colourful simplicial depth problem in dimension d is to find a configuration of (d+1) sets of (d+1) points such that the origin is contained in the convex hull of each set (colour) but contained in a minimal number of colourful simplices generated by taking one point from each set. A construction attaining d^2+1 simplices is known, and is conjectured to be minimal. This has been confirmed up to d=3, however the best known lower bound for d at least 4 is ((d+1)^2)/2.
A promising method to improve this lower bound is to look at combinatorial octahedral systems generated by such configurations. The difficulty to employing this approach is handling the many symmetric configurations that arise. We propose a table of invariants which exclude many of partial configurations, and use this to improve the lower bound in dimension 4.
△ Less
Submitted 29 October, 2012;
originally announced October 2012.
-
Counting inequivalent monotone Boolean functions
Authors:
Tamon Stephen,
Timothy Yusun
Abstract:
Monotone Boolean functions (MBFs) are Boolean functions $f: {0,1}^n \rightarrow {0,1}$ satisfying the monotonicity condition $x \leq y \Rightarrow f(x) \leq f(y)$ for any $x,y \in {0,1}^n$. The number of MBFs in n variables is known as the $n$th Dedekind number. It is a longstanding computational challenge to determine these numbers exactly - these values are only known for $n$ at most 8. Two mono…
▽ More
Monotone Boolean functions (MBFs) are Boolean functions $f: {0,1}^n \rightarrow {0,1}$ satisfying the monotonicity condition $x \leq y \Rightarrow f(x) \leq f(y)$ for any $x,y \in {0,1}^n$. The number of MBFs in n variables is known as the $n$th Dedekind number. It is a longstanding computational challenge to determine these numbers exactly - these values are only known for $n$ at most 8. Two monotone Boolean functions are inequivalent if one can be obtained from the other by renaming the variables. The number of inequivalent MBFs in $n$ variables was known only for up to $n = 6$. In this paper we propose a strategy to count inequivalent MBF's by breaking the calculation into parts based on the profiles of these functions. As a result we are able to compute the number of inequivalent MBFs in 7 variables. The number obtained is 490013148.
△ Less
Submitted 20 September, 2012;
originally announced September 2012.
-
Efficient Algorithms for Finding Tucker Patterns
Authors:
Cedric Chauve,
Tamon Stephen,
Maria Tamayo
Abstract:
The Consecutive Ones Property is an important notion for binary matrices, both from a theoretical and applied point of view. Tucker gave in 1972 a characterization of matrices that do not satisfy the Consecutive Ones Property in terms of forbidden submatrices, the Tucker patterns. We describe here a linear time algorithm to find a Tucker pattern in a non-C1P binary matrix, which allows to extract…
▽ More
The Consecutive Ones Property is an important notion for binary matrices, both from a theoretical and applied point of view. Tucker gave in 1972 a characterization of matrices that do not satisfy the Consecutive Ones Property in terms of forbidden submatrices, the Tucker patterns. We describe here a linear time algorithm to find a Tucker pattern in a non-C1P binary matrix, which allows to extract in linear time a certificate for the non-C1P. We also describe an output-sensitive algorithm to enumerate all Tucker patterns of a non-C1P binary matrix.
This paper had been withdrawn due to some missing cases in Algorithms 2 and 3.
△ Less
Submitted 30 June, 2012; v1 submitted 8 June, 2012;
originally announced June 2012.
-
A tight bound on the length of odd cycles in the incompatibility graph of a non-C1P matrix
Authors:
Mehrnoush Malekesmaeili,
Cedric Chauve,
Tamon Stephen
Abstract:
A binary matrix has the consecutive ones property (C1P) if it is possible to order the columns so that all 1s are consecutive in every row. In [McConnell, SODA 2004 768-777] the notion of incompatibility graph of a binary matrix was introduced and it was shown that odd cycles of this graph provide a certificate that a matrix does not have the consecutive ones property. A bound of (k+2) was claimed…
▽ More
A binary matrix has the consecutive ones property (C1P) if it is possible to order the columns so that all 1s are consecutive in every row. In [McConnell, SODA 2004 768-777] the notion of incompatibility graph of a binary matrix was introduced and it was shown that odd cycles of this graph provide a certificate that a matrix does not have the consecutive ones property. A bound of (k+2) was claimed for the smallest odd cycle of a non-C1P matrix with k columns. In this note we show that this result can be obtained simply and directly via Tucker patterns, and that the correct bound is (k+2) when k is even, but (k+3) when k is odd.
△ Less
Submitted 2 September, 2011;
originally announced September 2011.
-
Embedding a pair of graphs in a surface, and the width of 4-dimensional prismatoids
Authors:
Francisco Santos,
Tamon Stephen,
Hugh Thomas
Abstract:
A prismatoid is a polytope with all its vertices contained in two parallel facets, called its bases. Its width is the number of steps needed to go from one base to the other in the dual graph. The first author recently showed that the existence of counter-examples to the Hirsch conjecture is equivalent to that of $d$-prismatoids of width larger than $d$, and constructed such prismatoids in dimensi…
▽ More
A prismatoid is a polytope with all its vertices contained in two parallel facets, called its bases. Its width is the number of steps needed to go from one base to the other in the dual graph. The first author recently showed that the existence of counter-examples to the Hirsch conjecture is equivalent to that of $d$-prismatoids of width larger than $d$, and constructed such prismatoids in dimension five. Here we show that the same is impossible in dimension four. This is proved by looking at the pair of graph embeddings on a 2-sphere that arise from the normal fans of the two bases.
△ Less
Submitted 13 April, 2011;
originally announced April 2011.
-
An Euler characteristic proof that 4-prismatoids have width at most 4
Authors:
Tamon Stephen,
Hugh Thomas
Abstract:
We show that there is no 4-dimensional analogue of the 5-prismatoids used in Santos' recent counterexample to the Hirsch conjecture.
We show that there is no 4-dimensional analogue of the 5-prismatoids used in Santos' recent counterexample to the Hirsch conjecture.
△ Less
Submitted 15 April, 2011; v1 submitted 16 January, 2011;
originally announced January 2011.
-
Expected Crossing Numbers
Authors:
Bojan Mohar,
Tamon Stephen
Abstract:
The expected value for the weighted crossing number of a randomly weighted graph is studied. A variation of the Crossing Lemma for expectations is proved. We focus on the case where the edge-weights are independent random variables that are uniformly distributed on [0,1].
The expected value for the weighted crossing number of a randomly weighted graph is studied. A variation of the Crossing Lemma for expectations is proved. We focus on the case where the edge-weights are independent random variables that are uniformly distributed on [0,1].
△ Less
Submitted 12 December, 2010;
originally announced December 2010.
-
More Colourful Simplices
Authors:
Antoine Deza,
Tamon Stephen,
Feng Xie
Abstract:
We show that any point in the convex hull of each of (d+1) sets of (d+1) points in general position in \R^d is contained in at least (d+1)^2/2 simplices with one vertex from each set. This improves the known lower bounds for all d >= 4.
We show that any point in the convex hull of each of (d+1) sets of (d+1) points in general position in \R^d is contained in at least (d+1)^2/2 simplices with one vertex from each set. This improves the known lower bounds for all d >= 4.
△ Less
Submitted 31 May, 2010; v1 submitted 25 January, 2010;
originally announced January 2010.
-
Minimal Conflicting Sets for the Consecutive Ones Property in ancestral genome reconstruction
Authors:
Cedric Chauve,
Utz-Uwe Haus,
Tamon Stephen,
Vivija P. You
Abstract:
A binary matrix has the Consecutive Ones Property (C1P) if its columns can be ordered in such a way that all 1's on each row are consecutive. A Minimal Conflicting Set is a set of rows that does not have the C1P, but every proper subset has the C1P. Such submatrices have been considered in comparative genomics applications, but very little is known about their combinatorial structure and efficie…
▽ More
A binary matrix has the Consecutive Ones Property (C1P) if its columns can be ordered in such a way that all 1's on each row are consecutive. A Minimal Conflicting Set is a set of rows that does not have the C1P, but every proper subset has the C1P. Such submatrices have been considered in comparative genomics applications, but very little is known about their combinatorial structure and efficient algorithms to compute them. We first describe an algorithm that detects rows that belong to Minimal Conflicting Sets. This algorithm has a polynomial time complexity when the number of 1's in each row of the considered matrix is bounded by a constant. Next, we show that the problem of computing all Minimal Conflicting Sets can be reduced to the joint generation of all minimal true clauses and maximal false clauses for some monotone boolean function. We use these methods on simulated data related to ancestral genome reconstruction to show that computing Minimal Conflicting Set is useful in discriminating between true positive and false positive ancestral syntenies. We also study a dataset of yeast genomes and address the reliability of an ancestral genome proposal of the Saccahromycetaceae yeasts.
△ Less
Submitted 21 December, 2009;
originally announced December 2009.
-
Computing knock out strategies in metabolic networks
Authors:
Utz-Uwe Haus,
Steffen Klamt,
Tamon Stephen
Abstract:
Given a metabolic network in terms of its metabolites and reactions, our goal is to efficiently compute the minimal knock out sets of reactions required to block a given behaviour. We describe an algorithm which improves the computation of these knock out sets when the elementary modes (minimal functional subsystems) of the network are given. We also describe an algorithm which computes both the…
▽ More
Given a metabolic network in terms of its metabolites and reactions, our goal is to efficiently compute the minimal knock out sets of reactions required to block a given behaviour. We describe an algorithm which improves the computation of these knock out sets when the elementary modes (minimal functional subsystems) of the network are given. We also describe an algorithm which computes both the knock out sets and the elementary modes containing the blocked reactions directly from the description of the network and whose worst-case computational complexity is better than the algorithms currently in use for these problems. Computational results are included.
△ Less
Submitted 29 December, 2007;
originally announced January 2008.
-
A quadratic lower bound for colourful simplicial depth
Authors:
Tamon Stephen,
Hugh Thomas
Abstract:
We show that any point in the convex hull of each of (d+1) sets of (d+1) points in R^d is contained in at least floor((d+2)^2/4) simplices with one vertex from each set.
We show that any point in the convex hull of each of (d+1) sets of (d+1) points in R^d is contained in at least floor((d+2)^2/4) simplices with one vertex from each set.
△ Less
Submitted 23 July, 2007; v1 submitted 16 December, 2005;
originally announced December 2005.
-
The Colourful Feasibility Problem
Authors:
Antoine Deza,
Sui Huang,
Tamon Stephen,
Tamás Terlaky
Abstract:
We study a colourful generalization of the linear programming feasibility problem, comparing the algorithms introduced by Barany and Onn with new methods. We perform benchmarking on generic and ill-conditioned problems, as well as as recently introduced highly structured problems. We show that some algorithms can lead to cycling or slow convergence, but we provide extensive numerical experiments…
▽ More
We study a colourful generalization of the linear programming feasibility problem, comparing the algorithms introduced by Barany and Onn with new methods. We perform benchmarking on generic and ill-conditioned problems, as well as as recently introduced highly structured problems. We show that some algorithms can lead to cycling or slow convergence, but we provide extensive numerical experiments which show that others perform much better than predicted by complexity arguments. We conclude that the most efficient method for all but the most ill-conditioned problems is a proposed multi-update algorithm.
△ Less
Submitted 30 November, 2005;
originally announced November 2005.
-
Colourful Simplicial Depth
Authors:
Antoine Deza,
Sui Huang,
Tamon Stephen,
Tamás Terlaky
Abstract:
Inspired by Barany's colourful Caratheodory theorem, we introduce a colourful generalization of Liu's simplicial depth. We prove a parity property and conjecture that the minimum colourful simplicial depth of any core point in any d-dimensional configuration is d^2+1 and that the maximum is d^(d+1)+1. We exhibit configurations attaining each of these depths and apply our results to the problem o…
▽ More
Inspired by Barany's colourful Caratheodory theorem, we introduce a colourful generalization of Liu's simplicial depth. We prove a parity property and conjecture that the minimum colourful simplicial depth of any core point in any d-dimensional configuration is d^2+1 and that the maximum is d^(d+1)+1. We exhibit configurations attaining each of these depths and apply our results to the problem of bounding monochrome (non-colourful) simplicial depth.
△ Less
Submitted 8 January, 2006; v1 submitted 1 June, 2005;
originally announced June 2005.
-
A majorization bound for the eigenvalues of some graph Laplacians
Authors:
Tamon Stephen
Abstract:
It is conjectured that the Laplacian spectrum of a graph is majorized by its conjugate degree sequence. In this paper, we prove that this majorization holds for a class of graphs including trees. We also show that a generalization of this conjecture to graphs with Dirichlet boundary conditions is equivalent to the original conjecture.
It is conjectured that the Laplacian spectrum of a graph is majorized by its conjugate degree sequence. In this paper, we prove that this majorization holds for a class of graphs including trees. We also show that a generalization of this conjecture to graphs with Dirichlet boundary conditions is equivalent to the original conjecture.
△ Less
Submitted 7 November, 2004;
originally announced November 2004.
-
The distribution of values in the quadratic assignment problem
Authors:
Alexander Barvinok,
Tamon Stephen
Abstract:
We obtain a number of results regarding the distribution of values of a quadratic function f on the set of nxn permutation matrices (identified with the symmetric group S_n) around its optimum (minimum or maximum). In particular, we estimate the fraction of permutations sigma such that f(sigma) lies within a given neighborhood of the optimal value of f. We identify some ``extreme'' functions f (…
▽ More
We obtain a number of results regarding the distribution of values of a quadratic function f on the set of nxn permutation matrices (identified with the symmetric group S_n) around its optimum (minimum or maximum). In particular, we estimate the fraction of permutations sigma such that f(sigma) lies within a given neighborhood of the optimal value of f. We identify some ``extreme'' functions f (there are 4 of those for n even and 5 for n odd) such that the distribution of every quadratic function around its optimum is a certain ``mixture'' of the distributions of the extremes and describe a natural class of functions (which includes, for example, the objective function in the Traveling Salesman Problem) with a relative abundance of near-optimal permutations. In particular, we identify a large class of functions f with the property that permutations in the vicinity of the optimal permutation (in the Hamming metric of S_n) tend to produce near optimal values of f (such is, for example, the objective function in the symmetric Traveling Salesman Problem) and show that for general f, just the opposite behavior may take place: an average permutation in the vicinity of the optimal permutation may be much worse than an average permutation in the whole group S_n.
△ Less
Submitted 24 April, 2001;
originally announced April 2001.