Quantum Marginal Problems David Gross (Colgone) Joint with: - - PowerPoint PPT Presentation

quantum marginal problems
SMART_READER_LITE
LIVE PREVIEW

Quantum Marginal Problems David Gross (Colgone) Joint with: - - PowerPoint PPT Presentation

Quantum Marginal Problems David Gross (Colgone) Joint with: Christandl, Doran, Lopes, Schilling, Walter Outline Overview: Marginal problems Overview: Entanglement Main Theme: Entanglement Polytopes Shortly: Beyond the Pauli


slide-1
SLIDE 1

Quantum Marginal Problems

David Gross (Colgone) Joint with: Christandl, Doran, Lopes, Schilling, Walter

slide-2
SLIDE 2

Outline

◮ Overview: Marginal problems ◮ Overview: Entanglement ◮ Main Theme: Entanglement Polytopes ◮ Shortly: Beyond the Pauli principle.

slide-3
SLIDE 3

Overview: Marginal Problems

slide-4
SLIDE 4

Marginals

◮ A marginal is obtained by

integrating out parts of high-dim

  • bject
slide-5
SLIDE 5

Marginals

◮ A marginal is obtained by

integrating out parts of high-dim

  • bject

◮ Not every set of marginals is

compatible

slide-6
SLIDE 6

Marginals

◮ A marginal is obtained by

integrating out parts of high-dim

  • bject

◮ Not every set of marginals is

compatible

◮ Deciding compatibility is the

marginal problem

slide-7
SLIDE 7

Marginals in classical probability

◮ Marginals are distributions of subsets of variables.

slide-8
SLIDE 8

Marginals in classical probability

◮ Marginals are distributions of subsets of variables.

One classical marginal prob well-known in quantum:

slide-9
SLIDE 9

Marginals in classical probability

◮ Marginals are distributions of subsets of variables.

One classical marginal prob well-known in quantum: Bell tests.

slide-10
SLIDE 10

Bell tests as marginal problems

◮ There are four random variables:

polarization along two axes, as seen by Alice/Bob

slide-11
SLIDE 11

Bell tests as marginal problems

◮ There are four random variables:

polarization along two axes, as seen by Alice/Bob

◮ Only certain pairs accessible ◮ Q: Are these marginals compatible with classical distribution?

slide-12
SLIDE 12

Bell tests as marginal problems

◮ There are four random variables:

polarization along two axes, as seen by Alice/Bob

◮ Only certain pairs accessible ◮ Q: Are these marginals compatible with classical distribution? ◮ Compatible marginals form convex

polytope

◮ Facets are Bell inequalities.

slide-13
SLIDE 13

Bell tests as marginal problems

◮ There are four random variables:

polarization along two axes, as seen by Alice/Bob

◮ Only certain pairs accessible ◮ Q: Are these marginals compatible with classical distribution? ◮ Compatible marginals form convex

polytope

◮ Facets are Bell inequalities. ◮ Testing locality NP-hard ⇒ so is classical marginal problem

slide-14
SLIDE 14

Marginals in quantum theory

◮ For subset Si specify state ρi. ◮ Q: Are these compatible:

ρi = tr\Si ρ for some global ρ?

slide-15
SLIDE 15

Marginals in quantum theory

◮ For subset Si specify state ρi. ◮ Q: Are these compatible:

ρi = tr\Si ρ for some global ρ? Would solve all finite-dim. few-body ground-state probs!

slide-16
SLIDE 16

Marginals in quantum theory

◮ For subset Si specify state ρi. ◮ Q: Are these compatible:

ρi = tr\Si ρ for some global ρ? Would solve all finite-dim. few-body ground-state probs! E.g.: For two-body Hamiltonian H =

n

  • i,j=1

hi,j, compute min

ρ trHρ = min ρ

  • i,j

trhi,j ρ = min

{ρi,j} comp.

  • i,j

trhi,j ρi,j.

slide-17
SLIDE 17

Marginals in quantum theory: Ground States

min

ρ trHρ =

min

{ρi,j} comp.

  • i,j

trhi,j ρi,j. Remarks:

◮ Left-hand side optimizes over O(dn) variables. ◮ R.h.s. over O(n2 d4). Exponential improvement!

slide-18
SLIDE 18

Marginals in quantum theory: Ground States

min

ρ trHρ =

min

{ρi,j} comp.

  • i,j

trhi,j ρi,j. Remarks:

◮ Left-hand side optimizes over O(dn) variables. ◮ R.h.s. over O(n2 d4). Exponential improvement! ◮ Optimization over convex set of compatible ρi,j.

slide-19
SLIDE 19

Marginals in quantum theory: Ground States

min

ρ trHρ =

min

{ρi,j} comp.

  • i,j

trhi,j ρi,j. Remarks:

◮ Left-hand side optimizes over O(dn) variables. ◮ R.h.s. over O(n2 d4). Exponential improvement! ◮ Optimization over convex set of compatible ρi,j.

General theory of convex optimization ⇒ Computational complexity of 2-RDM method (r.h.s) domi- nated by deciding compatibility of ρi,j’s.

slide-20
SLIDE 20

Marginals in quantum theory: Ground States

min

ρ trHρ =

min

{ρi,j} comp.

  • i,j

trhi,j ρi,j. Remarks:

◮ Left-hand side optimizes over O(dn) variables. ◮ R.h.s. over O(n2 d4). Exponential improvement! ◮ Optimization over convex set of compatible ρi,j.

General theory of convex optimization ⇒ Computational complexity of 2-RDM method (r.h.s) domi- nated by deciding compatibility of ρi,j’s. Two directions:

◮ Progress on q. marginal prob. ⇒ info about ground states ◮ Hardness of ground-states ⇒ hardness of q. marginals.

slide-21
SLIDE 21

Negative direction

Computational complexity of 2-RDM method (r.h.s) domi- nated by deciding compatibility of ρi,j’s.

◮ But finding two-body ground-states is NP-hard ◮ (. . . and even QMA-hard)

slide-22
SLIDE 22

Negative direction

Computational complexity of 2-RDM method (r.h.s) domi- nated by deciding compatibility of ρi,j’s.

◮ But finding two-body ground-states is NP-hard ◮ (. . . and even QMA-hard)

Thus: There is no efficient algorithm (quantum or classical) for the general two-body quantum marginal problem.

slide-23
SLIDE 23

Negative direction

Computational complexity of 2-RDM method (r.h.s) domi- nated by deciding compatibility of ρi,j’s.

◮ But finding two-body ground-states is NP-hard ◮ (. . . and even QMA-hard)

Thus: There is no efficient algorithm (quantum or classical) for the general two-body quantum marginal problem.

◮ Remains hard for Fermions. ◮ Argument works for classical marginal prob. (hardness of Ising) ◮ Leaves room for outer approximations → D. Mazziotti’s talk.

slide-24
SLIDE 24

Negative direction

Computational complexity of 2-RDM method (r.h.s) domi- nated by deciding compatibility of ρi,j’s.

◮ But finding two-body ground-states is NP-hard ◮ (. . . and even QMA-hard)

Thus: There is no efficient algorithm (quantum or classical) for the general two-body quantum marginal problem.

◮ Remains hard for Fermions. ◮ Argument works for classical marginal prob. (hardness of Ising) ◮ Leaves room for outer approximations → D. Mazziotti’s talk.

Natural Question: Is there subproblem with enough structure to be tractable?

slide-25
SLIDE 25

1-RDM marginal problem

1-RDM subproblem: marginals do not overlap, global state pure

slide-26
SLIDE 26

1-RDM marginal problem

1-RDM subproblem: marginals do not overlap, global state pure Classical version:

◮ Globally pure

⇔ no global randomness ⇒ no local randomness.

◮ . . . trivial.

slide-27
SLIDE 27

1-RDM marginal problem

1-RDM subproblem: marginals do not overlap, global state pure Classical version:

◮ Globally pure

⇔ no global randomness ⇒ no local randomness.

◮ . . . trivial.

Quantum version:

◮ Globally pure

⇒ no local randomness (in presence of entanglement).

◮ . . . seems non-trivial, but tractable!

slide-28
SLIDE 28

1-RDM marginal problem

Questions to be asked:

◮ Structure of set of 1-RDMs? ◮ What info about global ψ accessible from 1-RDM? ◮ Computational complexity of 1-RDM marginal prob.? ◮ Practical uses?

slide-29
SLIDE 29

Structure of 1-RDMs.

slide-30
SLIDE 30

Reduction to eigenvalues

◮ Local basis change does not affect compatibility ◮ ⇒ can assume ρi are diagonal ⇒ described by eigenvalues

  • λ(1), . . . ,

λ(n) ∈ ❘dn.

slide-31
SLIDE 31

Reduction to eigenvalues

◮ Local basis change does not affect compatibility ◮ ⇒ can assume ρi are diagonal ⇒ described by eigenvalues

  • λ(1), . . . ,

λ(n) ∈ ❘dn. Question becomes: Which set of ordered local eigenvalues λ(i) can occur?

slide-32
SLIDE 32

Reduction to eigenvalues

◮ Local basis change does not affect compatibility ◮ ⇒ can assume ρi are diagonal ⇒ described by eigenvalues

  • λ(1), . . . ,

λ(n) ∈ ❘dn. Question becomes: Which set of ordered local eigenvalues λ(i) can occur? Deep fact:

◮ Compatible spectra form convex polytope

slide-33
SLIDE 33

Reduction to eigenvalues

◮ Local basis change does not affect compatibility ◮ ⇒ can assume ρi are diagonal ⇒ described by eigenvalues

  • λ(1), . . . ,

λ(n) ∈ ❘dn. Question becomes: Which set of ordered local eigenvalues λ(i) can occur? Deep fact:

◮ Compatible spectra form convex polytope ◮ Highly non-trivial! (Global set is not convex).

slide-34
SLIDE 34

Reduction to eigenvalues

◮ Local basis change does not affect compatibility ◮ ⇒ can assume ρi are diagonal ⇒ described by eigenvalues

  • λ(1), . . . ,

λ(n) ∈ ❘dn. Question becomes: Which set of ordered local eigenvalues λ(i) can occur? Deep fact:

◮ Compatible spectra form convex polytope ◮ Highly non-trivial! (Global set is not convex). ◮ Several proofs, building on symplectic

geometry & asymptotic rep theory

[Klyachko, Kirwan, Christandl, Mitchison, Harrow, Daftuar, Hayden, . . . ]

slide-35
SLIDE 35

Reduction to eigenvalues

◮ Local basis change does not affect compatibility ◮ ⇒ can assume ρi are diagonal ⇒ described by eigenvalues

  • λ(1), . . . ,

λ(n) ∈ ❘dn. Question becomes: Which set of ordered local eigenvalues λ(i) can occur? Deep fact:

◮ Compatible spectra form convex polytope ◮ Highly non-trivial! (Global set is not convex). ◮ Several proofs, building on symplectic

geometry & asymptotic rep theory

[Klyachko, Kirwan, Christandl, Mitchison, Harrow, Daftuar, Hayden, . . . ] ◮ No conceptually simple proof known to me!

slide-36
SLIDE 36

Example: d = n = 2

Warm up: work out solution for two qubits.

slide-37
SLIDE 37

Example: d = n = 2

Warm up: work out solution for two qubits.

◮ Schmidt-decomposition:

|ψ =

  • λ(1)|e1 ⊗ |f1 +
  • λ(2)|e2 ⊗ |f2

◮ With

ρ1 = λ(1)|e1e1|+λ(2)|e2e2|, ρ2 = λ(1)|f1f1|+λ(2)|f2f2|.

slide-38
SLIDE 38

Example: d = n = 2

Warm up: work out solution for two qubits.

◮ Schmidt-decomposition:

|ψ =

  • λ(1)|e1 ⊗ |f1 +
  • λ(2)|e2 ⊗ |f2

◮ With

ρ1 = λ(1)|e1e1|+λ(2)|e2e2|, ρ2 = λ(1)|f1f1|+λ(2)|f2f2|.

◮ So eigenvalues must be equal:

λ1 = λ2.

slide-39
SLIDE 39

Example: d = n = 2

Warm up: work out solution for two qubits.

◮ Schmidt-decomposition:

|ψ =

  • λ(1)|e1 ⊗ |f1 +
  • λ(2)|e2 ⊗ |f2

◮ With

ρ1 = λ(1)|e1e1|+λ(2)|e2e2|, ρ2 = λ(1)|f1f1|+λ(2)|f2f2|.

◮ So eigenvalues must be equal:

λ1 = λ2. In terms of largest eigenvalue, get simple polytope:

slide-40
SLIDE 40

Further examples

Three qubits: Three fermions on 6 modes (“Dennis-Borland”): (c.f. M. Christandl’s talk)

slide-41
SLIDE 41

Summary: Structure of 1-RDM’s

◮ Compatible 1-RDMs described by convex polytopes of spectra.

slide-42
SLIDE 42

Summary: Structure of 1-RDM’s

◮ Compatible 1-RDMs described by convex polytopes of spectra. ◮ If this doesn’t surprise you, I’m terribly sad.

slide-43
SLIDE 43

Computational aspects.

slide-44
SLIDE 44

List inequalities?

[Klyachko, Altunbulak] ◮ Polytopes characterized by finitely

many linear ineqs D, x ≤ c.

slide-45
SLIDE 45

List inequalities?

[Klyachko, Altunbulak] ◮ Polytopes characterized by finitely

many linear ineqs D, x ≤ c.

◮ Ansatz so far: Compute all ineqs → Altunbulak’s talk ◮ Doesn’t seem to scale: too many ineqs

as n, d go up.

slide-46
SLIDE 46

List inequalities?

[Klyachko, Altunbulak] ◮ There might be better algorithm than

“checking all ineqs”. ❘

slide-47
SLIDE 47

List inequalities?

[Klyachko, Altunbulak] ◮ There might be better algorithm than

“checking all ineqs”.

◮ Ex.: ℓ1-unit ball in ❘n has 2n linear

ineqs, but membership equivalent to xℓ1 =

n

  • i=1

|xi| ≤ 1.

slide-48
SLIDE 48

Computational complexity

Thus, central open question: Q.: Is there a poly-time algorithm that decides the 1-RDM quantum marginal problem?

slide-49
SLIDE 49

Computational complexity

Thus, central open question: Q.: Is there a poly-time algorithm that decides the 1-RDM quantum marginal problem? Progress Nov. 2015 [Burgisser, Christandl, Mulmuley, Walter]: Problem in NP ∩ coNP

◮ Virtually guarantees that it can’t be proven hard ◮ Suggests it might be in P.

slide-50
SLIDE 50

Info about global state from 1-RDMs. Part 1: Selection rules.

slide-51
SLIDE 51

Selection rules

Selection rule, “Generalized Hartree-Fock”: If a state ψ maps to the boundary of the polytope, only few, special Slater determinants can appear in an expansion of ψ.

v

(a)

v

(b)

v

(c)

v

(d)

slide-52
SLIDE 52

Selection rules

Selection rule, “Generalized Hartree-Fock”: If a state ψ maps to the boundary of the polytope, only few, special Slater determinants can appear in an expansion of ψ.

◮ Stated by Klyachko (2009). He

didn’t feel proof was necessary.

◮ True for for general scenarios –

stated here for Fermions.

v

(a)

v

(b)

v

(c)

v

(d)

[Schilling, DG, Christandl, PRL ’13]

slide-53
SLIDE 53

Selection rules

◮ Consider n-Fermion system with modes {φ1, . . . , φd}.

slide-54
SLIDE 54

Selection rules

◮ Consider n-Fermion system with modes {φ1, . . . , φd}. ◮ Expand ψ ∈ ∧nCd in Slater dets:

ψ =

  • i1<···<in

ci1,...,in φi1 ∧ · · · ∧ φin. (1)

slide-55
SLIDE 55

Selection rules

◮ Consider n-Fermion system with modes {φ1, . . . , φd}. ◮ Expand ψ ∈ ∧nCd in Slater dets:

ψ =

  • i1<···<in

ci1,...,in φi1 ∧ · · · ∧ φin. (1)

◮ Assume (wlog) ρ(1)(ψ) is diagonal with eigenvalues λ

slide-56
SLIDE 56

Selection rules

◮ Consider n-Fermion system with modes {φ1, . . . , φd}. ◮ Expand ψ ∈ ∧nCd in Slater dets:

ψ =

  • i1<···<in

ci1,...,in φi1 ∧ · · · ∧ φin. (1)

◮ Assume (wlog) ρ(1)(ψ) is diagonal with eigenvalues λ ◮ Let D, x ≤ c be a face of the 1-RDM polytope.

slide-57
SLIDE 57

Selection rules

◮ Consider n-Fermion system with modes {φ1, . . . , φd}. ◮ Expand ψ ∈ ∧nCd in Slater dets:

ψ =

  • i1<···<in

ci1,...,in φi1 ∧ · · · ∧ φin. (1)

◮ Assume (wlog) ρ(1)(ψ) is diagonal with eigenvalues λ ◮ Let D, x ≤ c be a face of the 1-RDM polytope.

If the eigenvalues of ψ saturate the ineq. D, λ = c, then (1)

  • nly contains Slater dets whose eigenvalues do so as well.
slide-58
SLIDE 58

Selection rules: Proof

If the eigenvalues of ψ saturate the ineq. D, λ = c, then (1)

  • nly contains Slater dets whose eigenvalues do so as well.

Elementary proof [Alex Lopes, PhD thesis; Lopes, Schilling, DG, in eternal prep.]

slide-59
SLIDE 59

Selection rules: Proof

If the eigenvalues of ψ saturate the ineq. D, λ = c, then (1)

  • nly contains Slater dets whose eigenvalues do so as well.

Elementary proof [Alex Lopes, PhD thesis; Lopes, Schilling, DG, in eternal prep.] Trick:

◮ Introduce operator ˆ

D =

i Di a† i ai. ◮ Then

D, λ = tr ˆ Dρ(1) = tr ˆ D|ψψ|.

slide-60
SLIDE 60

Selection rules: Proof

If the eigenvalues of ψ saturate the ineq. D, λ = c, then (1)

  • nly contains Slater dets whose eigenvalues do so as well.

Elementary proof [Alex Lopes, PhD thesis; Lopes, Schilling, DG, in eternal prep.] Trick:

◮ Introduce operator ˆ

D =

i Di a† i ai. ◮ Then

D, λ = tr ˆ Dρ(1) = tr ˆ D|ψψ|. Selection rule equivalent to: If tr ˆ D|ψψ| = c, then ˆ D|ψ = c|ψ. (Non-trivial, as c need not be extremal eigenvalue of ˆ D).

slide-61
SLIDE 61

Selection rules: Proof

If the eigenvalues of ψ saturate the ineq. D, λ = c, then (1)

  • nly contains Slater dets whose eigenvalues do so as well.

Elementary proof [Alex Lopes, PhD thesis; Lopes, Schilling, DG, in eternal prep.] Trick:

◮ Introduce operator ˆ

D =

i Di a† i ai. ◮ Then

D, λ = tr ˆ Dρ(1) = tr ˆ D|ψψ|. Selection rule equivalent to: If tr ˆ D|ψψ| = c, then ˆ D|ψ = c|ψ. (Non-trivial, as c need not be extremal eigenvalue of ˆ D). Proof: Blackboard.

slide-62
SLIDE 62

Info about global state from 1-RDMs. Part 2: Entanglement.

slide-63
SLIDE 63

Entanglement

◮ Two pure states ψ, φ are in same entanglement class if they

can be converted into each other with finite probability of success using local operations and classical communication. ❈ ❈

slide-64
SLIDE 64

Entanglement

◮ Two pure states ψ, φ are in same entanglement class if they

can be converted into each other with finite probability of success using local operations and classical communication.

◮ Often referred to as SLOCC classes. But that sounds too

unpleasant. ❈ ❈

slide-65
SLIDE 65

Entanglement

◮ Two pure states ψ, φ are in same entanglement class if they

can be converted into each other with finite probability of success using local operations and classical communication.

◮ Often referred to as SLOCC classes. But that sounds too

unpleasant.

◮ Formally:

ψ ∼ φ ⇔ ψ = (g1 ⊗ · · · ⊗ gn)φ with gi local invertible matrices (filtering operations).

◮ Mathematically: We’re looking at SL(❈d)×n-orbits in

  • ❈dn.
slide-66
SLIDE 66

SLOCC, SLOCC! – Who’s There?

◮ For three qubits (d = 2, n = 3), equivalence classes known

since mid-1800s. Re-discovered in 2000 to great effect:

slide-67
SLIDE 67

Examples

Classes:

◮ Products ψ = φ1 ⊗ φ2 ⊗ φ3. ◮ Three classes of bi-separable states: ψ = φ1 ⊗ φ2,3. ◮ The W-class:

|W = |001 + |010 + |100.

◮ The GHZ-class:

|GHZ = |000 + |111.

slide-68
SLIDE 68

Further examples

4 qubits:

◮ Classification apparently first obtained in QI community

[Verstraete et al. (2002)].

◮ Nine families of four complex parameters each.

slide-69
SLIDE 69

Further examples

4 qubits:

◮ Classification apparently first obtained in QI community

[Verstraete et al. (2002)].

◮ Nine families of four complex parameters each.

Beyond:

◮ Number of parameters required to label orbits increases

exponentially.

◮ Only sporadic facts known.

slide-70
SLIDE 70

Desiderata

Can we come up with theory that

◮ is systematic

(any number of particles, local dimensions, symmetry constraints),

◮ is efficient

(only polynomial number of parameters have to be learned),

◮ experimentally feasible

(parameters easily accessible, robust to noise)?

slide-71
SLIDE 71

Desiderata

Can we come up with theory that

◮ is systematic

(any number of particles, local dimensions, symmetry constraints),

◮ is efficient

(only polynomial number of parameters have to be learned),

◮ experimentally feasible

(parameters easily accessible, robust to noise)? Claim: The single-site quantum marginal problem lives up to these standards.

slide-72
SLIDE 72

Entanglement Polytopes

slide-73
SLIDE 73

Central observation, entanglement polytopes

Set of allowed eigenvalues may depend on entanglement class

  • f global state.
slide-74
SLIDE 74

Central observation, entanglement polytopes

Set of allowed eigenvalues may depend on entanglement class

  • f global state.

Thus:

◮ To every class C, associated set ∆C of local eigenvalues of

states in (closure of) C.

slide-75
SLIDE 75

Central observation, entanglement polytopes

Set of allowed eigenvalues may depend on entanglement class

  • f global state.

Thus:

◮ To every class C, associated set ∆C of local eigenvalues of

states in (closure of) C.

◮ Turns out: ∆C is again polytope: the entanglement polytope

associated with C.

[Walter, Doran, Gross, Christandl, Science 2013]

slide-76
SLIDE 76

Central observation, entanglement polytopes

Set of allowed eigenvalues may depend on entanglement class

  • f global state.

Thus:

◮ To every class C, associated set ∆C of local eigenvalues of

states in (closure of) C.

◮ Turns out: ∆C is again polytope: the entanglement polytope

associated with C.

[Walter, Doran, Gross, Christandl, Science 2013] ◮ Clearly: the position of

λ(ψ) w.r.t. the entanglement polytopes contains all local information about global entanglement class.

slide-77
SLIDE 77

Examples re-visited: 3 qubit entanglement polytopes

For three qubits, polytopes resolve all 6 entanglement classes: [Hang et al. (2004), Sawicki et al. (2012), our paper]

slide-78
SLIDE 78

Examples re-visited: 3 qubit entanglement polytopes

For three qubits, polytopes resolve all 6 entanglement classes: [Hang et al. (2004), Sawicki et al. (2012), our paper] W-class corresponds to “upper pyramid”: λ(1)

max + λ(2) max + λ(3) max ≥ 2.

Any violation of that witnesses GHZ-type entanglement.

slide-79
SLIDE 79

Examples re-visited: 4 qubit entanglement polytopes

4 qubits:

◮ Entanglement classes:

9 families with up to four complex parameters each [Verstraete et al. (2002)].

slide-80
SLIDE 80

Examples re-visited: 4 qubit entanglement polytopes

4 qubits:

◮ Entanglement classes:

9 families with up to four complex parameters each [Verstraete et al. (2002)].

◮ Entanglement Polytopes:

13 polytopes, 7 of which are genuinely 4-party entangled.

slide-81
SLIDE 81

Examples re-visited: 4 qubit entanglement polytopes

4 qubits:

◮ Entanglement classes:

9 families with up to four complex parameters each [Verstraete et al. (2002)].

◮ Entanglement Polytopes:

13 polytopes, 7 of which are genuinely 4-party entangled.

◮ We feel: attractive balance between coarse-graining and

preserving structure. Example: 4-qubit W-class CW ∋ |0001 + |0010 + |0100 + |1000 again an “upper pyramid”: λ(1)

max + λ(2) max + λ(3) max + λ(4) max ≥ 3.

slide-82
SLIDE 82

Example: 4 qubit entanglement polytopes

slide-83
SLIDE 83

Example: Bosonic qubits

Consider n bosonic qubits: ψ ∈ Symn ❈2 .

slide-84
SLIDE 84

Example: Bosonic qubits

Consider n bosonic qubits: ψ ∈ Symn ❈2 .

◮ Symmetry ⇒ all local reductions are equal:

ρ(1)

i,j = ψ|a† i aj|ψ. ◮ ⇒ single number captures all: λmax ∈ [0.5, 1].

slide-85
SLIDE 85

Example: Bosonic qubits

Consider n bosonic qubits: ψ ∈ Symn ❈2 .

◮ Symmetry ⇒ all local reductions are equal:

ρ(1)

i,j = ψ|a† i aj|ψ. ◮ ⇒ single number captures all: λmax ∈ [0.5, 1].

Analyze polytopes:

◮ |0, . . . , 0 in all C’s ⇒ ∆C = [γC, 1].

slide-86
SLIDE 86

Example: Bosonic qubits

Consider n bosonic qubits: ψ ∈ Symn ❈2 .

◮ Symmetry ⇒ all local reductions are equal:

ρ(1)

i,j = ψ|a† i aj|ψ. ◮ ⇒ single number captures all: λmax ∈ [0.5, 1].

Analyze polytopes:

◮ |0, . . . , 0 in all C’s ⇒ ∆C = [γC, 1]. ◮ Turns out: Possible choices are

γC ∈ 1 2

N − k N : k = 0, 1, . . . , ⌊N/2⌋

  • . . .

◮ . . . with innermost point γ the image of W -type states.

slide-87
SLIDE 87

Example: No Solipsism

◮ A vector is genuinely n-partite entangled if it does not

factorize w.r.t. any bi-partition: ψ = ψ1 ⊗ ψ2.

slide-88
SLIDE 88

Example: No Solipsism

◮ A vector is genuinely n-partite entangled if it does not

factorize w.r.t. any bi-partition: ψ = ψ1 ⊗ ψ2. Observation: sometimes detectable from local spectra alone.

slide-89
SLIDE 89

Example: No Solipsism

◮ A vector is genuinely n-partite entangled if it does not

factorize w.r.t. any bi-partition: ψ = ψ1 ⊗ ψ2. Observation: sometimes detectable from local spectra alone. ⇔ spectra ( λ(1), . . . , λ(n)) compatible, but no bi-partition is. Example: (λ(1)

max, . . . , λ(n) max) =

1 2 + 1 n − 1, 1 − 1 n − 1, . . . , 1 − 1 n − 1

  • .
slide-90
SLIDE 90

Example: No Solipsism

◮ A vector is genuinely n-partite entangled if it does not

factorize w.r.t. any bi-partition: ψ = ψ1 ⊗ ψ2. Observation: sometimes detectable from local spectra alone. ⇔ spectra ( λ(1), . . . , λ(n)) compatible, but no bi-partition is. Example: (λ(1)

max, . . . , λ(n) max) =

1 2 + 1 n − 1, 1 − 1 n − 1, . . . , 1 − 1 n − 1

  • .

Interpretation:

◮ no solipsism: love needs a partner!

(And entangled qubits need their counter-parts).

slide-91
SLIDE 91

Example: Distillation

Entanglement measures from local information:

◮ (Linear) entropy of entanglement

E(ψ) = 1 − 1 N

  • i

trρ2

i

simple function of Euclidean distance of eigenvalue point to origin.

◮ “Closer to origin ⇒ more entanglement”.

slide-92
SLIDE 92

Example: Distillation

Entanglement measures from local information:

◮ (Linear) entropy of entanglement

E(ψ) = 1 − 1 N

  • i

trρ2

i

simple function of Euclidean distance of eigenvalue point to origin.

◮ “Closer to origin ⇒ more entanglement”. ◮ ⇒ can bound distillable entanglement from local information!

slide-93
SLIDE 93

Example: Distillation

Entanglement measures from local information:

◮ (Linear) entropy of entanglement

E(ψ) = 1 − 1 N

  • i

trρ2

i

simple function of Euclidean distance of eigenvalue point to origin.

◮ “Closer to origin ⇒ more entanglement”. ◮ ⇒ can bound distillable entanglement from local information! ◮ Can even give distillation procedure without need to know

state beyond local densities (generalizing [Verstraete et al. 2002]).

slide-94
SLIDE 94

Pure???

◮ Yeah, but no pure state exists in Nature.

slide-95
SLIDE 95

Pure???

◮ Yeah, but no pure state exists in Nature. ◮ Results are epsilonifiable: if distance d of spectrum to a

polytope ∆ exceeds 4N

  • 1 − p,

then ρ ∈ conv(∆).

◮ p = tr ρ2 is purity, which an be lower-bounded from local

information alone.

slide-96
SLIDE 96

Summary of Entanglement Polytopes

◮ Locally accessible info about global entanglement encoded in

entanglement polytopes – subpolytopes of the set of admissible local spectra.

◮ Provides a systematic and efficient way of obtaining

information about entanglement classes.

slide-97
SLIDE 97

Thank you for your attention!

David Gross (Uni Cologne) Oxford, April 2016