bound-entangled 0.1.0__tar.gz

bound_entangled-0.1.0/LICENSE

@@ -0,0 +1,21 @@
+MIT License
+
+Copyright (c) 2026 Anton Van Cleemput
+
+Permission is hereby granted, free of charge, to any person obtaining a copy
+of this software and associated documentation files (the "Software"), to deal
+in the Software without restriction, including without limitation the rights
+to use, copy, modify, merge, publish, distribute, sublicense, and/or sell
+copies of the Software, and to permit persons to whom the Software is
+furnished to do so, subject to the following conditions:
+
+The above copyright notice and this permission notice shall be included in all
+copies or substantial portions of the Software.
+
+THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
+IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
+FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
+AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
+LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
+OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE
+SOFTWARE.

bound_entangled-0.1.0/PKG-INFO

@@ -0,0 +1,112 @@
+Metadata-Version: 2.4
+Name: bound-entangled
+Version: 0.1.0
+Summary: Families of bound entangled states for Quantum Information
+Author-email: AVlyx <anton@vancleem.com>
+License: MIT
+Project-URL: Homepage, https://github.com/AVlyx/bound-entangled
+Requires-Python: >=3.9
+Description-Content-Type: text/markdown
+License-File: LICENSE
+Requires-Dist: toqito>=1.3.0
+Requires-Dist: numpy>=2.4.1
+Provides-Extra: test
+Requires-Dist: pytest; extra == "test"
+Provides-Extra: docs
+Requires-Dist: mkdocs; extra == "docs"
+Requires-Dist: mkdocs-material; extra == "docs"
+Requires-Dist: mkdocstrings[python]; extra == "docs"
+Dynamic: license-file
+
+# bound-entangled
+
+Reference implementations of several families of **bound entangled states** from
+the quantum information literature, as plain NumPy density matrices.
+
+A bound entangled state is entangled yet has a positive partial transpose (PPT),
+so no pure entanglement can be distilled from it. This library collects
+constructions that are otherwise scattered across papers into a single, tested
+package, each factory returning the density matrix `rho` of the state.
+
+## Installation
+
+```bash
+pip install -e .
+# with the test dependencies:
+pip install -e ".[test]"
+```
+
+Requires Python ≥ 3.9. Depends on [`toqito`](https://toqito.readthedocs.io) and
+NumPy.
+
+## Usage
+
+Packages are organized by the bipartite Hilbert space the states live in
+(`c3_otimes_c3` = C³ ⊗ C³, `cm_otimes_cn` = C^m ⊗ C^n, etc.).
+
+```python
+from c3_otimes_c3 import chessboard_extremal_PPT, tiles_upb, pyramid_upb
+from cd_otimes_cd import yu_oh
+from c5_otimes_c5 import sn3_grid_state
+
+rho = chessboard_extremal_PPT()   # 9x9 PPT-entangled chessboard state
+rho = yu_oh(d=3, x=0.5, y=0.5)    # Yu-Oh nonlocal bound entangled state
+```
+
+## States
+
+### `c3_otimes_c3` — C³ ⊗ C³
+
+| Factory | State | Reference |
+|---|---|---|
+| `chessboard`, `chessboard_extremal_PPT` | Bruß–Peres chessboard states | [quant-ph/9911056](https://arxiv.org/abs/quant-ph/9911056) |
+| `cross_hatch` | 3×3 "cross-hatch" grid state (CCNR-detected) | [1705.09261](https://arxiv.org/abs/1705.09261) |
+| `steering_state` | Steerable bound entangled state (counterexample to the stronger Peres conjecture) | [1405.0262](https://arxiv.org/abs/1405.0262) |
+| `tiles_upb`, `pyramid_upb` | States from the Tiles / Pyramid unextendible product bases | [quant-ph/9808030](https://arxiv.org/abs/quant-ph/9808030) |
+| `parametrized_upb` | Six-parameter family of UPBs generalizing Tiles / Pyramid | [quant-ph/9908070](https://arxiv.org/abs/quant-ph/9908070) |
+| `horodecki` | 3x3 Horodecki bound entangled state | [quant-ph/9703004](https://arxiv.org/abs/quant-ph/9703004) |
+
+
+### `c4_otimes_c4` — C⁴ ⊗ C⁴
+
+| Factory | State | Reference |
+|---|---|---|
+| `pianni` | 4×4 Benatti–Floreanini–Piani state | [quant-ph/0411095](https://arxiv.org/abs/quant-ph/0411095) |
+
+### `c5_otimes_c5` — C⁵ ⊗ C⁵
+
+| Factory | State | Reference |
+|---|---|---|
+| `sn3_grid_state` | Smallest known Schmidt-number-3 PPT bound entangled state | [2402.12966](https://arxiv.org/abs/2402.12966) |
+
+### `cd_otimes_cd` — C^d ⊗ C^d
+
+| Factory | State | Reference |
+|---|---|---|
+| `yu_oh`, `is_valid_yu_oh_input` | Yu–Oh family of nonlocal bound entangled states | [1509.08991](https://arxiv.org/abs/1509.08991) |
+| `gen_tiles1` | GenTiles1 UPB generalizing Tiles to d⊗d, even d≥4 | [quant-ph/9908070](https://arxiv.org/abs/quant-ph/9908070) |
+
+### `cm_otimes_cn` — C^m ⊗ C^n (parametric constructions)
+
+| Factory | State | Reference |
+|---|---|---|
+| `grid_state` | Quantum grid states from graph edges | [1705.09261](https://arxiv.org/abs/1705.09261) |
+| `generalized_grid_state` | Grid states generalized to hyperedges | [2402.12966](https://arxiv.org/abs/2402.12966) |
+| `random_NPT`, `random_PPT`, `random_PPT_close_to_the_PPT_edge` | Random density matrices by PPT class (rejection sampling) | — |
+| `gen_tiles2` | GenTiles2 UPB generalizing Tiles to m⊗n, n>3, m≥3, n≥m | [quant-ph/9908070](https://arxiv.org/abs/quant-ph/9908070) |
+
+
+## Documentation
+
+Full API docs (built with [MkDocs](https://www.mkdocs.org/) +
+[mkdocstrings](https://mkdocstrings.github.io/)) are published to GitHub Pages
+on every push to `main`. To preview locally:
+
+```bash
+pip install -e ".[docs]"
+mkdocs serve
+```
+
+## License
+
+MIT

bound_entangled-0.1.0/README.md

@@ -0,0 +1,92 @@
+# bound-entangled
+
+Reference implementations of several families of **bound entangled states** from
+the quantum information literature, as plain NumPy density matrices.
+
+A bound entangled state is entangled yet has a positive partial transpose (PPT),
+so no pure entanglement can be distilled from it. This library collects
+constructions that are otherwise scattered across papers into a single, tested
+package, each factory returning the density matrix `rho` of the state.
+
+## Installation
+
+```bash
+pip install -e .
+# with the test dependencies:
+pip install -e ".[test]"
+```
+
+Requires Python ≥ 3.9. Depends on [`toqito`](https://toqito.readthedocs.io) and
+NumPy.
+
+## Usage
+
+Packages are organized by the bipartite Hilbert space the states live in
+(`c3_otimes_c3` = C³ ⊗ C³, `cm_otimes_cn` = C^m ⊗ C^n, etc.).
+
+```python
+from c3_otimes_c3 import chessboard_extremal_PPT, tiles_upb, pyramid_upb
+from cd_otimes_cd import yu_oh
+from c5_otimes_c5 import sn3_grid_state
+
+rho = chessboard_extremal_PPT()   # 9x9 PPT-entangled chessboard state
+rho = yu_oh(d=3, x=0.5, y=0.5)    # Yu-Oh nonlocal bound entangled state
+```
+
+## States
+
+### `c3_otimes_c3` — C³ ⊗ C³
+
+| Factory | State | Reference |
+|---|---|---|
+| `chessboard`, `chessboard_extremal_PPT` | Bruß–Peres chessboard states | [quant-ph/9911056](https://arxiv.org/abs/quant-ph/9911056) |
+| `cross_hatch` | 3×3 "cross-hatch" grid state (CCNR-detected) | [1705.09261](https://arxiv.org/abs/1705.09261) |
+| `steering_state` | Steerable bound entangled state (counterexample to the stronger Peres conjecture) | [1405.0262](https://arxiv.org/abs/1405.0262) |
+| `tiles_upb`, `pyramid_upb` | States from the Tiles / Pyramid unextendible product bases | [quant-ph/9808030](https://arxiv.org/abs/quant-ph/9808030) |
+| `parametrized_upb` | Six-parameter family of UPBs generalizing Tiles / Pyramid | [quant-ph/9908070](https://arxiv.org/abs/quant-ph/9908070) |
+| `horodecki` | 3x3 Horodecki bound entangled state | [quant-ph/9703004](https://arxiv.org/abs/quant-ph/9703004) |
+
+
+### `c4_otimes_c4` — C⁴ ⊗ C⁴
+
+| Factory | State | Reference |
+|---|---|---|
+| `pianni` | 4×4 Benatti–Floreanini–Piani state | [quant-ph/0411095](https://arxiv.org/abs/quant-ph/0411095) |
+
+### `c5_otimes_c5` — C⁵ ⊗ C⁵
+
+| Factory | State | Reference |
+|---|---|---|
+| `sn3_grid_state` | Smallest known Schmidt-number-3 PPT bound entangled state | [2402.12966](https://arxiv.org/abs/2402.12966) |
+
+### `cd_otimes_cd` — C^d ⊗ C^d
+
+| Factory | State | Reference |
+|---|---|---|
+| `yu_oh`, `is_valid_yu_oh_input` | Yu–Oh family of nonlocal bound entangled states | [1509.08991](https://arxiv.org/abs/1509.08991) |
+| `gen_tiles1` | GenTiles1 UPB generalizing Tiles to d⊗d, even d≥4 | [quant-ph/9908070](https://arxiv.org/abs/quant-ph/9908070) |
+
+### `cm_otimes_cn` — C^m ⊗ C^n (parametric constructions)
+
+| Factory | State | Reference |
+|---|---|---|
+| `grid_state` | Quantum grid states from graph edges | [1705.09261](https://arxiv.org/abs/1705.09261) |
+| `generalized_grid_state` | Grid states generalized to hyperedges | [2402.12966](https://arxiv.org/abs/2402.12966) |
+| `random_NPT`, `random_PPT`, `random_PPT_close_to_the_PPT_edge` | Random density matrices by PPT class (rejection sampling) | — |
+| `gen_tiles2` | GenTiles2 UPB generalizing Tiles to m⊗n, n>3, m≥3, n≥m | [quant-ph/9908070](https://arxiv.org/abs/quant-ph/9908070) |
+
+
+## Documentation
+
+Full API docs (built with [MkDocs](https://www.mkdocs.org/) +
+[mkdocstrings](https://mkdocstrings.github.io/)) are published to GitHub Pages
+on every push to `main`. To preview locally:
+
+```bash
+pip install -e ".[docs]"
+mkdocs serve
+```
+
+## License
+
+MIT

bound_entangled-0.1.0/pyproject.toml

@@ -0,0 +1,37 @@
+[build-system]
+requires = ["setuptools>=68.0"]
+build-backend = "setuptools.build_meta"
+
+[project]
+name = "bound-entangled"
+version = "0.1.0"
+description = "Families of bound entangled states for Quantum Information"
+readme = "README.md"
+requires-python = ">=3.9"
+license = {text = "MIT"}
+authors = [
+    {name = "AVlyx", email = "anton@vancleem.com"}
+]
+dependencies = [
+    "toqito>=1.3.0",
+    "numpy>=2.4.1"
+]
+
+[project.urls]
+Homepage = "https://github.com/AVlyx/bound-entangled"
+
+[project.optional-dependencies]
+test = ["pytest"]
+docs = ["mkdocs", "mkdocs-material", "mkdocstrings[python]"]
+
+[tool.setuptools.packages.find]
+where = ["src"]
+
+[tool.pytest.ini_options]
+testpaths = ["tests"]
+# src: import the library the same flat way the source does (`from utils import ...`)
+#      even without the editable install. tests: makes the shared `qi` helper
+#      module importable from every test file.
+pythonpath = ["src", "tests"]
+# importlib mode lets the two identically named `test_upb.py` modules coexist.
+addopts = "--import-mode=importlib"

bound_entangled-0.1.0/setup.cfg

@@ -0,0 +1,4 @@
+[egg_info]
+tag_build = 
+tag_date = 0
+

bound_entangled-0.1.0/src/bound_entangled.egg-info/PKG-INFO

@@ -0,0 +1,112 @@
+Metadata-Version: 2.4
+Name: bound-entangled
+Version: 0.1.0
+Summary: Families of bound entangled states for Quantum Information
+Author-email: AVlyx <anton@vancleem.com>
+License: MIT
+Project-URL: Homepage, https://github.com/AVlyx/bound-entangled
+Requires-Python: >=3.9
+Description-Content-Type: text/markdown
+License-File: LICENSE
+Requires-Dist: toqito>=1.3.0
+Requires-Dist: numpy>=2.4.1
+Provides-Extra: test
+Requires-Dist: pytest; extra == "test"
+Provides-Extra: docs
+Requires-Dist: mkdocs; extra == "docs"
+Requires-Dist: mkdocs-material; extra == "docs"
+Requires-Dist: mkdocstrings[python]; extra == "docs"
+Dynamic: license-file
+
+# bound-entangled
+
+Reference implementations of several families of **bound entangled states** from
+the quantum information literature, as plain NumPy density matrices.
+
+A bound entangled state is entangled yet has a positive partial transpose (PPT),
+so no pure entanglement can be distilled from it. This library collects
+constructions that are otherwise scattered across papers into a single, tested
+package, each factory returning the density matrix `rho` of the state.
+
+## Installation
+
+```bash
+pip install -e .
+# with the test dependencies:
+pip install -e ".[test]"
+```
+
+Requires Python ≥ 3.9. Depends on [`toqito`](https://toqito.readthedocs.io) and
+NumPy.
+
+## Usage
+
+Packages are organized by the bipartite Hilbert space the states live in
+(`c3_otimes_c3` = C³ ⊗ C³, `cm_otimes_cn` = C^m ⊗ C^n, etc.).
+
+```python
+from c3_otimes_c3 import chessboard_extremal_PPT, tiles_upb, pyramid_upb
+from cd_otimes_cd import yu_oh
+from c5_otimes_c5 import sn3_grid_state
+
+rho = chessboard_extremal_PPT()   # 9x9 PPT-entangled chessboard state
+rho = yu_oh(d=3, x=0.5, y=0.5)    # Yu-Oh nonlocal bound entangled state
+```
+
+## States
+
+### `c3_otimes_c3` — C³ ⊗ C³
+
+| Factory | State | Reference |
+|---|---|---|
+| `chessboard`, `chessboard_extremal_PPT` | Bruß–Peres chessboard states | [quant-ph/9911056](https://arxiv.org/abs/quant-ph/9911056) |
+| `cross_hatch` | 3×3 "cross-hatch" grid state (CCNR-detected) | [1705.09261](https://arxiv.org/abs/1705.09261) |
+| `steering_state` | Steerable bound entangled state (counterexample to the stronger Peres conjecture) | [1405.0262](https://arxiv.org/abs/1405.0262) |
+| `tiles_upb`, `pyramid_upb` | States from the Tiles / Pyramid unextendible product bases | [quant-ph/9808030](https://arxiv.org/abs/quant-ph/9808030) |
+| `parametrized_upb` | Six-parameter family of UPBs generalizing Tiles / Pyramid | [quant-ph/9908070](https://arxiv.org/abs/quant-ph/9908070) |
+| `horodecki` | 3x3 Horodecki bound entangled state | [quant-ph/9703004](https://arxiv.org/abs/quant-ph/9703004) |
+
+
+### `c4_otimes_c4` — C⁴ ⊗ C⁴
+
+| Factory | State | Reference |
+|---|---|---|
+| `pianni` | 4×4 Benatti–Floreanini–Piani state | [quant-ph/0411095](https://arxiv.org/abs/quant-ph/0411095) |
+
+### `c5_otimes_c5` — C⁵ ⊗ C⁵
+
+| Factory | State | Reference |
+|---|---|---|
+| `sn3_grid_state` | Smallest known Schmidt-number-3 PPT bound entangled state | [2402.12966](https://arxiv.org/abs/2402.12966) |
+
+### `cd_otimes_cd` — C^d ⊗ C^d
+
+| Factory | State | Reference |
+|---|---|---|
+| `yu_oh`, `is_valid_yu_oh_input` | Yu–Oh family of nonlocal bound entangled states | [1509.08991](https://arxiv.org/abs/1509.08991) |
+| `gen_tiles1` | GenTiles1 UPB generalizing Tiles to d⊗d, even d≥4 | [quant-ph/9908070](https://arxiv.org/abs/quant-ph/9908070) |
+
+### `cm_otimes_cn` — C^m ⊗ C^n (parametric constructions)
+
+| Factory | State | Reference |
+|---|---|---|
+| `grid_state` | Quantum grid states from graph edges | [1705.09261](https://arxiv.org/abs/1705.09261) |
+| `generalized_grid_state` | Grid states generalized to hyperedges | [2402.12966](https://arxiv.org/abs/2402.12966) |
+| `random_NPT`, `random_PPT`, `random_PPT_close_to_the_PPT_edge` | Random density matrices by PPT class (rejection sampling) | — |
+| `gen_tiles2` | GenTiles2 UPB generalizing Tiles to m⊗n, n>3, m≥3, n≥m | [quant-ph/9908070](https://arxiv.org/abs/quant-ph/9908070) |
+
+
+## Documentation
+
+Full API docs (built with [MkDocs](https://www.mkdocs.org/) +
+[mkdocstrings](https://mkdocstrings.github.io/)) are published to GitHub Pages
+on every push to `main`. To preview locally:
+
+```bash
+pip install -e ".[docs]"
+mkdocs serve
+```
+
+## License
+
+MIT

bound_entangled-0.1.0/src/bound_entangled.egg-info/SOURCES.txt

@@ -0,0 +1,32 @@
+LICENSE
+README.md
+pyproject.toml
+src/bound_entangled.egg-info/PKG-INFO
+src/bound_entangled.egg-info/SOURCES.txt
+src/bound_entangled.egg-info/dependency_links.txt
+src/bound_entangled.egg-info/requires.txt
+src/bound_entangled.egg-info/top_level.txt
+src/c3_otimes_c3/__init__.py
+src/c3_otimes_c3/chessboard.py
+src/c3_otimes_c3/cross_hatch.py
+src/c3_otimes_c3/horodecki.py
+src/c3_otimes_c3/steering.py
+src/c3_otimes_c3/upb/__init__.py
+src/c3_otimes_c3/upb/parametrized_UPB.py
+src/c3_otimes_c3/upb/pyramid_UPB.py
+src/c3_otimes_c3/upb/tiles_UPB.py
+src/c4_otimes_c4/__init__.py
+src/c4_otimes_c4/pianni.py
+src/c5_otimes_c5/__init__.py
+src/c5_otimes_c5/sn3_grid_state.py
+src/cd_otimes_cd/__init__.py
+src/cd_otimes_cd/gen_tiles.py
+src/cd_otimes_cd/yu_oh.py
+src/cm_otimes_cn/__init__.py
+src/cm_otimes_cn/gen_tiles.py
+src/cm_otimes_cn/generalized_grid_state.py
+src/cm_otimes_cn/grid_state.py
+src/cm_otimes_cn/random_not_bound_entangled.py
+src/utils/__init__.py
+src/utils/ketbra.py
+src/utils/upb.py

bound_entangled-0.1.0/src/bound_entangled.egg-info/dependency_links.txt

@@ -0,0 +1 @@
+

bound_entangled-0.1.0/src/bound_entangled.egg-info/requires.txt

@@ -0,0 +1,10 @@
+toqito>=1.3.0
+numpy>=2.4.1
+
+[docs]
+mkdocs
+mkdocs-material
+mkdocstrings[python]
+
+[test]
+pytest

bound_entangled-0.1.0/src/bound_entangled.egg-info/top_level.txt

@@ -0,0 +1,6 @@
+c3_otimes_c3
+c4_otimes_c4
+c5_otimes_c5
+cd_otimes_cd
+cm_otimes_cn
+utils

bound_entangled-0.1.0/src/c3_otimes_c3/__init__.py

@@ -0,0 +1,7 @@
+from c3_otimes_c3.chessboard import chessboard, chessboard_extremal_PPT
+from c3_otimes_c3.cross_hatch import cross_hatch
+from c3_otimes_c3.steering import steering_state
+from c3_otimes_c3.upb.tiles_UPB import tiles_upb
+from c3_otimes_c3.upb.pyramid_UPB import pyramid_upb
+from c3_otimes_c3.upb.parametrized_UPB import parametrized_upb
+from c3_otimes_c3.horodecki import horodecki

bound_entangled-0.1.0/src/c3_otimes_c3/chessboard.py

@@ -0,0 +1,57 @@
+"""https://arxiv.org/abs/quant-ph/9911056
+
+D. Bruß, A. Peres, "Construction of quantum states with bound entanglement",
+Phys. Rev. A 61, 030301(R) (2000). Named "chessboard states" because their
+8x8 matrix form looks like a chessboard (see also
+https://arxiv.org/abs/2010.08372, Appendix C2).
+"""
+
+from toqito.states import chessboard as chessboard_toqito
+from typing import Optional
+import numpy as np
+
+
+def chessboard(
+    a: complex, b: complex, c: complex, d: complex, m: complex, n: complex, s: Optional[complex] = None, t: Optional[complex] = None
+) -> np.ndarray:
+    """A chessboard state.
+
+    Args:
+        a: free real parameter.
+        b: free real parameter.
+        c: free real parameter.
+        d: free real parameter.
+        m: free real parameter.
+        n: free real parameter.
+        s: defaults to c/n if None.
+        t: defaults to a*d/m if None.
+
+    Returns:
+        np.ndarray: a chessboard quantum state.
+    """
+    return chessboard_toqito([a, b, c, d, m, n], s, t)  # type: ignore
+
+
+def chessboard_extremal_PPT() -> np.ndarray:
+    """Extremal PPT chessboard state.
+
+    The state from the chessboard family with m = n = b = -3/5, a = 3/5,
+    c = -d = 6/5, shown to be an extremal point of the PPT-entangled states
+    in https://arxiv.org/abs/2010.08372 (Fig. 2 / Fig. 9, page 25).
+
+    Returns:
+        np.ndarray: the extremal PPT chessboard state.
+    """
+    return chessboard(
+        # m = n = b = -3/5
+        m=-3 / 5,
+        n=-3 / 5,
+        b=-3 / 5,
+        # a = 3/5
+        a=3 / 5,
+        # c = -d = 6/5
+        c=6 / 5,
+        d=-6 / 5,
+        s=(3 / 5) * (6 / 5) / (-3 / 5),
+        t=(3 / 5) * (-6 / 5) / (-3 / 5),
+    )

bound_entangled-0.1.0/src/c3_otimes_c3/cross_hatch.py

@@ -0,0 +1,25 @@
+"""https://arxiv.org/abs/1705.09261
+
+J. Lockhart, O. Gühne, S. Severini, "Entanglement properties of quantum grid
+states", Phys. Rev. A 97, 062340 (2018).
+"""
+
+from cm_otimes_cn.grid_state import grid_state
+
+
+def cross_hatch():
+    """The 3x3 "cross-hatch" grid state.
+
+    A bound entangled example of a quantum grid state, detected by the CCNR
+    (realignment) criterion.
+
+    Returns:
+        np.ndarray: the cross-hatch grid state.
+    """
+    return grid_state(
+        (3, 3),
+        ((0, 0), (1, 2)),
+        ((1, 0), (2, 2)),
+        ((0, 1), (2, 0)),
+        ((0, 2), (2, 1)),
+    )

bound_entangled-0.1.0/src/c3_otimes_c3/horodecki.py

@@ -0,0 +1,23 @@
+"""https://arxiv.org/abs/quant-ph/9703004
+
+P. Horodecki, "Separability criterion and inseparable mixed states with
+positive partial transposition", Phys. Lett. A 232, 333 (1997), Section 4.1.
+One of the first known families of bound entangled states: PPT for every
+a_param in [0, 1], separable only at the endpoints a_param = 0 or 1.
+"""
+
+from toqito.states import horodecki as horodecki_toqito
+import numpy as np
+
+
+def horodecki(a_param: float) -> np.ndarray:
+    """The 3x3 Horodecki state.
+
+    Args:
+        a_param: free real parameter in [0, 1]. The state is PPT for every
+            value in this range, and separable only at a_param = 0 or 1.
+
+    Returns:
+        np.ndarray: the 3x3 Horodecki density matrix.
+    """
+    return horodecki_toqito(a_param, dim=[3, 3])  # type: ignore

bound_entangled-0.1.0/src/c3_otimes_c3/steering.py

@@ -0,0 +1,42 @@
+"""https://arxiv.org/abs/1405.0262
+
+T. Moroder, O. Gittsovich, M. Huber, O. Gühne, "Steering Bound Entangled
+States: A Counterexample to the Stronger Peres Conjecture", Phys. Rev. Lett.
+113, 050404 (2014).
+"""
+
+import numpy as np
+from math import sqrt
+from utils import ketbra
+
+
+def steering_state(m1: float, m2: float) -> np.ndarray:
+    """A 3x3 bound entangled state that is steerable, parametrized by m1, m2.
+
+    Counterexample to the (stronger) Peres conjecture that bound entangled
+    states cannot be used to demonstrate EPR steering.
+
+    Args:
+        m1: first free real parameter.
+        m2: second free real parameter (with m1**2 + m2**2 <= 1).
+
+    Returns:
+        np.ndarray: the steerable bound entangled state.
+    """
+    m3 = sqrt((1 - m1**2 - m2**2) / 2)
+
+    psi_1 = np.zeros(9, dtype=np.float64)
+    psi_2 = np.zeros(9, dtype=np.float64)
+    psi_3 = np.zeros(9, dtype=np.float64)
+    psi_t_3 = np.zeros(9, dtype=np.float64)
+
+    psi_1[1 * 3 + 2], psi_1[2 * 3 + 1] = 1 / sqrt(2), 1 / sqrt(2)
+    psi_2[0], psi_2[1 * 3 + 1], psi_2[2 * 3 + 2] = 1 / sqrt(3), 1 / sqrt(3), -1 / sqrt(3)
+    psi_3[1], psi_3[1 * 3], psi_3[1 * 3 + 1], psi_3[2 * 3 + 2] = m1, m2, m3, m3
+    psi_t_3[2], psi_t_3[2 * 3], psi_t_3[2 * 3 + 1], psi_t_3[1 * 3 + 2] = m1, -m2, m3, -m3
+
+    lambda1 = 1 - (2 + 3 * m1 * m2) / (4 - 2 * m1**2 + m1 * m2 - 2 * m2**2)
+    lambda3 = 1 / (4 - 2 * m1**2 + m1 * m2 - 2 * m2**2)
+    lambda2 = 1 - lambda1 - 2 * lambda3
+
+    return lambda1 * ketbra(psi_1, psi_1) + lambda2 * ketbra(psi_2, psi_2) + lambda3 * (ketbra(psi_3, psi_3) + ketbra(psi_t_3, psi_t_3))