Download Is a random state entangled ?

Survey
yes no Was this document useful for you?
   Thank you for your participation!

* Your assessment is very important for improving the workof artificial intelligence, which forms the content of this project

Document related concepts

Hydrogen atom wikipedia , lookup

Copenhagen interpretation wikipedia , lookup

Quantum dot wikipedia , lookup

Quantum field theory wikipedia , lookup

Quantum decoherence wikipedia , lookup

Path integral formulation wikipedia , lookup

Topological quantum field theory wikipedia , lookup

Coherent states wikipedia , lookup

Delayed choice quantum eraser wikipedia , lookup

Scalar field theory wikipedia , lookup

Max Born wikipedia , lookup

Measurement in quantum mechanics wikipedia , lookup

Quantum fiction wikipedia , lookup

Bell test experiments wikipedia , lookup

Quantum electrodynamics wikipedia , lookup

Many-worlds interpretation wikipedia , lookup

Quantum computing wikipedia , lookup

Orchestrated objective reduction wikipedia , lookup

Bell's theorem wikipedia , lookup

Density matrix wikipedia , lookup

T-symmetry wikipedia , lookup

Canonical quantization wikipedia , lookup

Interpretations of quantum mechanics wikipedia , lookup

EPR paradox wikipedia , lookup

Probability amplitude wikipedia , lookup

Symmetry in quantum mechanics wikipedia , lookup

Quantum machine learning wikipedia , lookup

History of quantum field theory wikipedia , lookup

Quantum group wikipedia , lookup

Quantum teleportation wikipedia , lookup

Quantum state wikipedia , lookup

Quantum key distribution wikipedia , lookup

Hidden variable theory wikipedia , lookup

Quantum entanglement wikipedia , lookup

Transcript
Is a random state entangled ?
Guillaume AUBRUN
(joint work with Stanislaw Szarek and Deping Ye)
Université Lyon 1, France
Guillaume Aubrun (Lyon)
Random quantum states
Roscoff, June 2012
1 / 18
The probabilistic method in quantum information theory
The probabilistic method: use random techniques to show existence of
objects for which explicit constructions are not known (Erdős–Rényi ;
Shannon). Major revolution in combinatorics since 1950s !
A priori, Quantum Information Theory is a good candidate for applying the
probabilistic method, because the objects have very very large dimensions.
While the “curse of dimensionality” makes numerical method impractical,
high dimensions usually boost random techniques (“blessom of
dimensionality”).
In QIT, the probabilistic method has been applied successfully by Hastings
(counterexample to the additivity conjecture). An important pioneering
work by Hayden–Leung–Winter: “Aspects of generic entanglement”.
This motivates the study the properties of random states or random
channels, which might become a basic tool in future years.
Guillaume Aubrun (Lyon)
Random quantum states
Roscoff, June 2012
2 / 18
The probabilistic method in quantum information theory
The probabilistic method: use random techniques to show existence of
objects for which explicit constructions are not known (Erdős–Rényi ;
Shannon). Major revolution in combinatorics since 1950s !
A priori, Quantum Information Theory is a good candidate for applying the
probabilistic method, because the objects have very very large dimensions.
While the “curse of dimensionality” makes numerical method impractical,
high dimensions usually boost random techniques (“blessom of
dimensionality”).
In QIT, the probabilistic method has been applied successfully by Hastings
(counterexample to the additivity conjecture). An important pioneering
work by Hayden–Leung–Winter: “Aspects of generic entanglement”.
This motivates the study the properties of random states or random
channels, which might become a basic tool in future years.
Guillaume Aubrun (Lyon)
Random quantum states
Roscoff, June 2012
2 / 18
The probabilistic method in quantum information theory
The probabilistic method: use random techniques to show existence of
objects for which explicit constructions are not known (Erdős–Rényi ;
Shannon). Major revolution in combinatorics since 1950s !
A priori, Quantum Information Theory is a good candidate for applying the
probabilistic method, because the objects have very very large dimensions.
While the “curse of dimensionality” makes numerical method impractical,
high dimensions usually boost random techniques (“blessom of
dimensionality”).
In QIT, the probabilistic method has been applied successfully by Hastings
(counterexample to the additivity conjecture). An important pioneering
work by Hayden–Leung–Winter: “Aspects of generic entanglement”.
This motivates the study the properties of random states or random
channels, which might become a basic tool in future years.
Guillaume Aubrun (Lyon)
Random quantum states
Roscoff, June 2012
2 / 18
The probabilistic method in quantum information theory
The probabilistic method: use random techniques to show existence of
objects for which explicit constructions are not known (Erdős–Rényi ;
Shannon). Major revolution in combinatorics since 1950s !
A priori, Quantum Information Theory is a good candidate for applying the
probabilistic method, because the objects have very very large dimensions.
While the “curse of dimensionality” makes numerical method impractical,
high dimensions usually boost random techniques (“blessom of
dimensionality”).
In QIT, the probabilistic method has been applied successfully by Hastings
(counterexample to the additivity conjecture). An important pioneering
work by Hayden–Leung–Winter: “Aspects of generic entanglement”.
This motivates the study the properties of random states or random
channels, which might become a basic tool in future years.
Guillaume Aubrun (Lyon)
Random quantum states
Roscoff, June 2012
2 / 18
States in quantum information theory
Definition
A (mixed) state on Cn is a positive, self-adjoint, trace 1 operator on Cn .
A pure state is a state with rank 1, i.e. a one-dimensional projector. A
pure state is denoted |ψihψ|, where ψ is a unit vector in the range.
Pure states are the extreme points of the set D = D(Cn ) of all states.
D(Cn ) is a convex body in the (affine) space of trace 1 self-adjoint
operators on Cn . The Euclidean structure is given by the Hilbert–Schmidt
inner product hA, Bi = Tr AB, and dim D(Cn ) = n2 − 1.
The role of the origin is played by Id/n, the maximally mixed state. We
write for example
kId/nkD = 0,
where k · kD stands for the gauge of convex body D (which is not centrally
symmetric).
Guillaume Aubrun (Lyon)
Random quantum states
Roscoff, June 2012
3 / 18
States in quantum information theory
Definition
A (mixed) state on Cn is a positive, self-adjoint, trace 1 operator on Cn .
A pure state is a state with rank 1, i.e. a one-dimensional projector. A
pure state is denoted |ψihψ|, where ψ is a unit vector in the range.
Pure states are the extreme points of the set D = D(Cn ) of all states.
D(Cn ) is a convex body in the (affine) space of trace 1 self-adjoint
operators on Cn . The Euclidean structure is given by the Hilbert–Schmidt
inner product hA, Bi = Tr AB, and dim D(Cn ) = n2 − 1.
The role of the origin is played by Id/n, the maximally mixed state. We
write for example
kId/nkD = 0,
where k · kD stands for the gauge of convex body D (which is not centrally
symmetric).
Guillaume Aubrun (Lyon)
Random quantum states
Roscoff, June 2012
3 / 18
States in quantum information theory
Definition
A (mixed) state on Cn is a positive, self-adjoint, trace 1 operator on Cn .
A pure state is a state with rank 1, i.e. a one-dimensional projector. A
pure state is denoted |ψihψ|, where ψ is a unit vector in the range.
Pure states are the extreme points of the set D = D(Cn ) of all states.
D(Cn ) is a convex body in the (affine) space of trace 1 self-adjoint
operators on Cn . The Euclidean structure is given by the Hilbert–Schmidt
inner product hA, Bi = Tr AB, and dim D(Cn ) = n2 − 1.
The role of the origin is played by Id/n, the maximally mixed state. We
write for example
kId/nkD = 0,
where k · kD stands for the gauge of convex body D (which is not centrally
symmetric).
Guillaume Aubrun (Lyon)
Random quantum states
Roscoff, June 2012
3 / 18
Random induced states
Mixed states on D(Cn ) can be obtained as partial traces of pure states on
a larger system.
Definition (random “induced” states)
Let n, s be integers, |ψihψ| be a random Haar-distributed pure state on
Cn ⊗ Cs , and ρ = TrCs |ψihψ| (partial trace over the environment). Then
ρ is a random state on Cn . We denote by µn,s be its distribution.
Miracle: µn,n is the normalized Lebesgue measure on D(Cn )
(Życzkowski–Sommers).
More generally, for s > n, the measure µn,s has a log-concave density,
proportionnal to (det ρ)s−n .
Guillaume Aubrun (Lyon)
Random quantum states
Roscoff, June 2012
4 / 18
Random induced states
Mixed states on D(Cn ) can be obtained as partial traces of pure states on
a larger system.
Definition (random “induced” states)
Let n, s be integers, |ψihψ| be a random Haar-distributed pure state on
Cn ⊗ Cs , and ρ = TrCs |ψihψ| (partial trace over the environment). Then
ρ is a random state on Cn . We denote by µn,s be its distribution.
Miracle: µn,n is the normalized Lebesgue measure on D(Cn )
(Życzkowski–Sommers).
More generally, for s > n, the measure µn,s has a log-concave density,
proportionnal to (det ρ)s−n .
Guillaume Aubrun (Lyon)
Random quantum states
Roscoff, June 2012
4 / 18
Random induced states
Mixed states on D(Cn ) can be obtained as partial traces of pure states on
a larger system.
Definition (random “induced” states)
Let n, s be integers, |ψihψ| be a random Haar-distributed pure state on
Cn ⊗ Cs , and ρ = TrCs |ψihψ| (partial trace over the environment). Then
ρ is a random state on Cn . We denote by µn,s be its distribution.
Miracle: µn,n is the normalized Lebesgue measure on D(Cn )
(Życzkowski–Sommers).
More generally, for s > n, the measure µn,s has a log-concave density,
proportionnal to (det ρ)s−n .
Guillaume Aubrun (Lyon)
Random quantum states
Roscoff, June 2012
4 / 18
Wishart random matrices
Uniform measure on the sphere = normalized Gaussian vector
Fact (Induced states = normalized Wishart matrices)
If W is a n × s random matrix with independent NC (0, 1) entries, then
WW ∗
Tr(WW ∗ )
has distribution µn,s .
What happens when n if fixed and s → ∞ ?
Law of large numbers : the measure µn,s converge to the state Id/n.
Central limit approximation : the measure µn,s can be approximated by
Id
1
+ √ G
n
n s
where G is a GUE0 matrix (GUE conditionned to have trace 0).
Guillaume Aubrun (Lyon)
Random quantum states
Roscoff, June 2012
5 / 18
Wishart random matrices
Uniform measure on the sphere = normalized Gaussian vector
Fact (Induced states = normalized Wishart matrices)
If W is a n × s random matrix with independent NC (0, 1) entries, then
WW ∗
Tr(WW ∗ )
has distribution µn,s .
What happens when n if fixed and s → ∞ ?
Law of large numbers : the measure µn,s converge to the state Id/n.
Central limit approximation : the measure µn,s can be approximated by
Id
1
+ √ G
n
n s
where G is a GUE0 matrix (GUE conditionned to have trace 0).
Guillaume Aubrun (Lyon)
Random quantum states
Roscoff, June 2012
5 / 18
Wishart random matrices
Uniform measure on the sphere = normalized Gaussian vector
Fact (Induced states = normalized Wishart matrices)
If W is a n × s random matrix with independent NC (0, 1) entries, then
WW ∗
Tr(WW ∗ )
has distribution µn,s .
What happens when n if fixed and s → ∞ ?
Law of large numbers : the measure µn,s converge to the state Id/n.
Central limit approximation : the measure µn,s can be approximated by
Id
1
+ √ G
n
n s
where G is a GUE0 matrix (GUE conditionned to have trace 0).
Guillaume Aubrun (Lyon)
Random quantum states
Roscoff, June 2012
5 / 18
Wishart random matrices
Uniform measure on the sphere = normalized Gaussian vector
Fact (Induced states = normalized Wishart matrices)
If W is a n × s random matrix with independent NC (0, 1) entries, then
WW ∗
Tr(WW ∗ )
has distribution µn,s .
What happens when n if fixed and s → ∞ ?
Law of large numbers : the measure µn,s converge to the state Id/n.
Central limit approximation : the measure µn,s can be approximated by
Id
1
+ √ G
n
n s
where G is a GUE0 matrix (GUE conditionned to have trace 0).
Guillaume Aubrun (Lyon)
Random quantum states
Roscoff, June 2012
5 / 18
Central limit approximation made rigorous
Proposition
Let K ⊂ D(Cn ) be a convex body. Let ρ be a state with distribution µn,s
and G a n × n GUE0 matrix. Then, when n and s/n tend to infinity,
1
E kρkK ∼ √ E kId/n + G kK .
n s
Two ingredients in the proof:
1
2
Both sides converge to the semicircle distribution (Marčenko–Pastur
distributions degenerate to semicircles when s/n become large), and
extreme eigenvalues are under control (Bai–Yin).
Majorization theory: A ≺ B if A ∈ conv{UBU ∗ : U unitary}.
A ≺ B ⇐⇒ ∀t ∈ R, Tr |A − tId| 6 Tr |B − tId|.
Can be checked from an approximate knowledge of spectra.
Guillaume Aubrun (Lyon)
Random quantum states
Roscoff, June 2012
6 / 18
Central limit approximation made rigorous
Proposition
Let K ⊂ D(Cn ) be a convex body. Let ρ be a state with distribution µn,s
and G a n × n GUE0 matrix. Then, when n and s/n tend to infinity,
1
E kρkK ∼ √ E kId/n + G kK .
n s
Two ingredients in the proof:
1
2
Both sides converge to the semicircle distribution (Marčenko–Pastur
distributions degenerate to semicircles when s/n become large), and
extreme eigenvalues are under control (Bai–Yin).
Majorization theory: A ≺ B if A ∈ conv{UBU ∗ : U unitary}.
A ≺ B ⇐⇒ ∀t ∈ R, Tr |A − tId| 6 Tr |B − tId|.
Can be checked from an approximate knowledge of spectra.
Guillaume Aubrun (Lyon)
Random quantum states
Roscoff, June 2012
6 / 18
Central limit approximation made rigorous
Proposition
Let K ⊂ D(Cn ) be a convex body. Let ρ be a state with distribution µn,s
and G a n × n GUE0 matrix. Then, when n and s/n tend to infinity,
1
E kρkK ∼ √ E kId/n + G kK .
n s
Two ingredients in the proof:
1
2
Both sides converge to the semicircle distribution (Marčenko–Pastur
distributions degenerate to semicircles when s/n become large), and
extreme eigenvalues are under control (Bai–Yin).
Majorization theory: A ≺ B if A ∈ conv{UBU ∗ : U unitary}.
A ≺ B ⇐⇒ ∀t ∈ R, Tr |A − tId| 6 Tr |B − tId|.
Can be checked from an approximate knowledge of spectra.
Guillaume Aubrun (Lyon)
Random quantum states
Roscoff, June 2012
6 / 18
Central limit approximation made rigorous
Proposition
Let K ⊂ D(Cn ) be a convex body. Let ρ be a state with distribution µn,s
and G a n × n GUE0 matrix. Then, when n and s/n tend to infinity,
1
E kρkK ∼ √ E kId/n + G kK .
n s
Two ingredients in the proof:
1
2
Both sides converge to the semicircle distribution (Marčenko–Pastur
distributions degenerate to semicircles when s/n become large), and
extreme eigenvalues are under control (Bai–Yin).
Majorization theory: A ≺ B if A ∈ conv{UBU ∗ : U unitary}.
A ≺ B ⇐⇒ ∀t ∈ R, Tr |A − tId| 6 Tr |B − tId|.
Can be checked from an approximate knowledge of spectra.
Guillaume Aubrun (Lyon)
Random quantum states
Roscoff, June 2012
6 / 18
Separable states
We now replace Cn by Cd ⊗ Cd (n = d 2 ) and introduce the entanglement
vs separability dichotomy.
Definition (entanglement vs. separability)
A product state on Cd ⊗ Cd is a state of the form ρ1 ⊗ ρ2 , where ρ1 , ρ2
are states on Cd .
A state on Cd ⊗ Cd is separable if it can be written as a convex
combination of product states — otherwise it is entangled.
Let S be the set of separable states on Cd ⊗ Cd .
S = conv{|ψ1 ⊗ ψ2 ihψ1 ⊗ ψ2 | : ψ1 , ψ2 ∈ Cd , |ψ1 | = |ψ2 | = 1}.
Fact
The dimension of S is d 4 − 1, and equals the dimension of D.
Even more: D and S have the same inradius (Barnum–Gurvits).
Guillaume Aubrun (Lyon)
Random quantum states
Roscoff, June 2012
7 / 18
Separable states
We now replace Cn by Cd ⊗ Cd (n = d 2 ) and introduce the entanglement
vs separability dichotomy.
Definition (entanglement vs. separability)
A product state on Cd ⊗ Cd is a state of the form ρ1 ⊗ ρ2 , where ρ1 , ρ2
are states on Cd .
A state on Cd ⊗ Cd is separable if it can be written as a convex
combination of product states — otherwise it is entangled.
Let S be the set of separable states on Cd ⊗ Cd .
S = conv{|ψ1 ⊗ ψ2 ihψ1 ⊗ ψ2 | : ψ1 , ψ2 ∈ Cd , |ψ1 | = |ψ2 | = 1}.
Fact
The dimension of S is d 4 − 1, and equals the dimension of D.
Even more: D and S have the same inradius (Barnum–Gurvits).
Guillaume Aubrun (Lyon)
Random quantum states
Roscoff, June 2012
7 / 18
Separable states
We now replace Cn by Cd ⊗ Cd (n = d 2 ) and introduce the entanglement
vs separability dichotomy.
Definition (entanglement vs. separability)
A product state on Cd ⊗ Cd is a state of the form ρ1 ⊗ ρ2 , where ρ1 , ρ2
are states on Cd .
A state on Cd ⊗ Cd is separable if it can be written as a convex
combination of product states — otherwise it is entangled.
Let S be the set of separable states on Cd ⊗ Cd .
S = conv{|ψ1 ⊗ ψ2 ihψ1 ⊗ ψ2 | : ψ1 , ψ2 ∈ Cd , |ψ1 | = |ψ2 | = 1}.
Fact
The dimension of S is d 4 − 1, and equals the dimension of D.
Even more: D and S have the same inradius (Barnum–Gurvits).
Guillaume Aubrun (Lyon)
Random quantum states
Roscoff, June 2012
7 / 18
The problem: (when) is a random state entangled ?
The dichotomy entangled vs separable is fundamental in QIT (useful vs
useless). This motivates the following question.
Question
Let ρ be a random state on Cd ⊗ Cd with distribution µd 2 ,s . For which
values of d, s is ρ typically separable ?
1
If s = 1, ρ is entangled with probability 1.
2
If s ≫ d, ρ is separable with large probability (law of large numbers).
So we expect a threshold between typical entanglement and typical
separability.
Guillaume Aubrun (Lyon)
Random quantum states
Roscoff, June 2012
8 / 18
The problem: (when) is a random state entangled ?
The dichotomy entangled vs separable is fundamental in QIT (useful vs
useless). This motivates the following question.
Question
Let ρ be a random state on Cd ⊗ Cd with distribution µd 2 ,s . For which
values of d, s is ρ typically separable ?
1
If s = 1, ρ is entangled with probability 1.
2
If s ≫ d, ρ is separable with large probability (law of large numbers).
So we expect a threshold between typical entanglement and typical
separability.
Guillaume Aubrun (Lyon)
Random quantum states
Roscoff, June 2012
8 / 18
The problem: (when) is a random state entangled ?
The dichotomy entangled vs separable is fundamental in QIT (useful vs
useless). This motivates the following question.
Question
Let ρ be a random state on Cd ⊗ Cd with distribution µd 2 ,s . For which
values of d, s is ρ typically separable ?
1
If s = 1, ρ is entangled with probability 1.
2
If s ≫ d, ρ is separable with large probability (law of large numbers).
So we expect a threshold between typical entanglement and typical
separability.
Guillaume Aubrun (Lyon)
Random quantum states
Roscoff, June 2012
8 / 18
The problem: (when) is a random state entangled ?
The dichotomy entangled vs separable is fundamental in QIT (useful vs
useless). This motivates the following question.
Question
Let ρ be a random state on Cd ⊗ Cd with distribution µd 2 ,s . For which
values of d, s is ρ typically separable ?
1
If s = 1, ρ is entangled with probability 1.
2
If s ≫ d, ρ is separable with large probability (law of large numbers).
So we expect a threshold between typical entanglement and typical
separability.
Guillaume Aubrun (Lyon)
Random quantum states
Roscoff, June 2012
8 / 18
Our main theorem
We prove that the threshold is sharp and holds somewhere between d 3
and d 3 log2 d.
Theorem ( A+Szarek+Ye )
There is a function s0 (d), satisfying d 3 . s0 (d) . d 3 log2 d, such that if ρ
is a random state on Cd ⊗ Cd with distribution µd 2 ,s , and ε > 0,
1
If s 6 (1 − ε)s0 (d), then ρ is entangled w.h.p.
2
If s > (1 + ε)s0 (d), then ρ is separable w.h.p.
w.h.p. = with high probability, i.e. P(·) > 1 − C (ε) exp(−c(ε)s).
Known before: ρ entangled for s 6 d 2 and ρ separable for s > Cd 4
(Hayden–Leung–Winter).
Guillaume Aubrun (Lyon)
Random quantum states
Roscoff, June 2012
9 / 18
Our main theorem
We prove that the threshold is sharp and holds somewhere between d 3
and d 3 log2 d.
Theorem ( A+Szarek+Ye )
There is a function s0 (d), satisfying d 3 . s0 (d) . d 3 log2 d, such that if ρ
is a random state on Cd ⊗ Cd with distribution µd 2 ,s , and ε > 0,
1
If s 6 (1 − ε)s0 (d), then ρ is entangled w.h.p.
2
If s > (1 + ε)s0 (d), then ρ is separable w.h.p.
w.h.p. = with high probability, i.e. P(·) > 1 − C (ε) exp(−c(ε)s).
Known before: ρ entangled for s 6 d 2 and ρ separable for s > Cd 4
(Hayden–Leung–Winter).
Guillaume Aubrun (Lyon)
Random quantum states
Roscoff, June 2012
9 / 18
Convex geometry
Let k · kS be the gauge of the convex body S. Then
ρ is separable ⇐⇒ kρkS 6 1
ρ is entangled ⇐⇒ kρkS > 1
If ρ is a random state with distribution µd 2 ,s , we have to decide whether
the random variable kρkS is typically larger or smaller than 1.
Guillaume Aubrun (Lyon)
Random quantum states
Roscoff, June 2012
10 / 18
Convex geometry
Let k · kS be the gauge of the convex body S. Then
ρ is separable ⇐⇒ kρkS 6 1
ρ is entangled ⇐⇒ kρkS > 1
If ρ is a random state with distribution µd 2 ,s , we have to decide whether
the random variable kρkS is typically larger or smaller than 1.
Guillaume Aubrun (Lyon)
Random quantum states
Roscoff, June 2012
10 / 18
The strategy of proof
We proceed via two independent steps.
Decide for which values of d, s, one has
1
EkρkS < 1 − ε or EkρkS > 1 + ε.
Show that the function ρ 7→ kρkS concentrate arounds its mean.
2
Step 2 is relatively routine, except for the following simple idea.
Idea
If f : S N−1 → R is a 1-Lipschitz function, and Ω ⊂ S N−1 a subset of large
measure such that f|Ω is ε-Lipschitz, one may consider
f˜ = any ε-Lipschitz extension of f|Ω to S N−1 .
and apply standard concentration inequalities to f˜.
Guillaume Aubrun (Lyon)
Random quantum states
Roscoff, June 2012
11 / 18
The strategy of proof
We proceed via two independent steps.
Decide for which values of d, s, one has
1
EkρkS < 1 − ε or EkρkS > 1 + ε.
Show that the function ρ 7→ kρkS concentrate arounds its mean.
2
Step 2 is relatively routine, except for the following simple idea.
Idea
If f : S N−1 → R is a 1-Lipschitz function, and Ω ⊂ S N−1 a subset of large
measure such that f|Ω is ε-Lipschitz, one may consider
f˜ = any ε-Lipschitz extension of f|Ω to S N−1 .
and apply standard concentration inequalities to f˜.
Guillaume Aubrun (Lyon)
Random quantum states
Roscoff, June 2012
11 / 18
Mean width
Definition
If K ⊂ RN is a convex body containing the origin, its mean width is
Z
Z
w (K ) =
maxhx, θidθ =
kθkK ◦ dθ.
S N−1 x∈K
S N−1
We have to compute EkρkS . By our proposition about Gaussian
approximation, if s d 2 , and G a d 2 × d 2 GUE0 random matrix.
1
√ EkId/d 2 + G kS .
s
The right-hand side is the
√ Gaussian mean width of the dual (polar) convex
body S ◦ , which equals dim ≈ d 2 times the usual mean width. Therefore
1
EkρkS ∼ √ w (S ◦ ).
s
EkρkS ∼
d2
This shows that a sharp threshold occurs at the value s0 = w (S ◦ )2 .
Guillaume Aubrun (Lyon)
Random quantum states
Roscoff, June 2012
12 / 18
Mean width
Definition
If K ⊂ RN is a convex body containing the origin, its mean width is
Z
Z
w (K ) =
maxhx, θidθ =
kθkK ◦ dθ.
S N−1 x∈K
S N−1
We have to compute EkρkS . By our proposition about Gaussian
approximation, if s d 2 , and G a d 2 × d 2 GUE0 random matrix.
1
√ EkId/d 2 + G kS .
s
The right-hand side is the
√ Gaussian mean width of the dual (polar) convex
body S ◦ , which equals dim ≈ d 2 times the usual mean width. Therefore
1
EkρkS ∼ √ w (S ◦ ).
s
EkρkS ∼
d2
This shows that a sharp threshold occurs at the value s0 = w (S ◦ )2 .
Guillaume Aubrun (Lyon)
Random quantum states
Roscoff, June 2012
12 / 18
Mean width
Definition
If K ⊂ RN is a convex body containing the origin, its mean width is
Z
Z
w (K ) =
maxhx, θidθ =
kθkK ◦ dθ.
S N−1 x∈K
S N−1
We have to compute EkρkS . By our proposition about Gaussian
approximation, if s d 2 , and G a d 2 × d 2 GUE0 random matrix.
1
√ EkId/d 2 + G kS .
s
The right-hand side is the
√ Gaussian mean width of the dual (polar) convex
body S ◦ , which equals dim ≈ d 2 times the usual mean width. Therefore
1
EkρkS ∼ √ w (S ◦ ).
s
EkρkS ∼
d2
This shows that a sharp threshold occurs at the value s0 = w (S ◦ )2 .
Guillaume Aubrun (Lyon)
Random quantum states
Roscoff, June 2012
12 / 18
Computing the mean width of S ◦
To complete the proof it remains to compute the mean width of S ◦ .
Proposition
For absolute constants c, C , we have
cd 3/2 6 w (S ◦ ) 6 Cd 3/2 log d
Challenge: remove the log(d) factor.
It turns out that w (S) is simple to estimate, instead of w (S ◦ ).
Lemma
For absolute constants c, C , we have
cd −3/2 6 w (S) 6 Cd −3/2 .
Guillaume Aubrun (Lyon)
Random quantum states
Roscoff, June 2012
13 / 18
Computing the mean width of S ◦
To complete the proof it remains to compute the mean width of S ◦ .
Proposition
For absolute constants c, C , we have
cd 3/2 6 w (S ◦ ) 6 Cd 3/2 log d
Challenge: remove the log(d) factor.
It turns out that w (S) is simple to estimate, instead of w (S ◦ ).
Lemma
For absolute constants c, C , we have
cd −3/2 6 w (S) 6 Cd −3/2 .
Guillaume Aubrun (Lyon)
Random quantum states
Roscoff, June 2012
13 / 18
Estimating w (S)
Fact (union bound)
If P ⊂ B2N is a polytope with v vertices, then
r
log v
w (P) 6 C
.
N
Here N = d 4 − 1, and the convex body S can be approximated by a
polytope with 100d vertices. Hence
w (S) 6 Cd −3/2 .
For the lower bound, we use the Urysohn inequality w (S) > vrad(S). The
volume of S is large because
S the symmetrized body (with one extra
dimension) Σ = conv(−S S) has a large inradius. By the
Rogers–Shephard inequality, the volumes of S and Σ are comparable.
Guillaume Aubrun (Lyon)
Random quantum states
Roscoff, June 2012
14 / 18
Estimating w (S)
Fact (union bound)
If P ⊂ B2N is a polytope with v vertices, then
r
log v
w (P) 6 C
.
N
Here N = d 4 − 1, and the convex body S can be approximated by a
polytope with 100d vertices. Hence
w (S) 6 Cd −3/2 .
For the lower bound, we use the Urysohn inequality w (S) > vrad(S). The
volume of S is large because
S the symmetrized body (with one extra
dimension) Σ = conv(−S S) has a large inradius. By the
Rogers–Shephard inequality, the volumes of S and Σ are comparable.
Guillaume Aubrun (Lyon)
Random quantum states
Roscoff, June 2012
14 / 18
Duality argument: the MM ∗ -estimate
One trivially has w (K )w (K ◦ ) > 1. Conversely, one has
Theorem (The MM ∗ -estimate; Pisier, Figiel–Tomczak-Jaegermann)
If K ⊂ RN is a symmetric convex body in the `-position, then
1 6 w (K )w (K ◦ ) 6 C log N.
Every symmetric convex body has an linear image which is in the
`-position. A guarantee for K being in the `-position is that the group of
isometries of K acts irreducibly.
Guillaume Aubrun (Lyon)
Random quantum states
Roscoff, June 2012
15 / 18
Applying the MM ∗ -estimate
We cannot apply the MM ∗ -estimate to S directly because
1
S is not symmetric,
2
Its group of isometries does not act irreducibly.
We introduce another symmetrization, Ssym = S ∩ −S. It is checked that
S and Ssym have comparable volumes, mean widths and dual mean widths.
For the second point, the isometry group contains “local unitaries”
ρ 7→ (U ⊗ V )ρ(U ∗ ⊗ V ∗ )
and has an irreducible subspace of small codimension: the space
span{σ1 ⊗ σ2 : Tr σ1 = Tr σ2 = 0}.
We can apply the MM ∗ -estimate into this subspace, and this is enough.
Guillaume Aubrun (Lyon)
Random quantum states
Roscoff, June 2012
16 / 18
Applying the MM ∗ -estimate
We cannot apply the MM ∗ -estimate to S directly because
1
S is not symmetric,
2
Its group of isometries does not act irreducibly.
We introduce another symmetrization, Ssym = S ∩ −S. It is checked that
S and Ssym have comparable volumes, mean widths and dual mean widths.
For the second point, the isometry group contains “local unitaries”
ρ 7→ (U ⊗ V )ρ(U ∗ ⊗ V ∗ )
and has an irreducible subspace of small codimension: the space
span{σ1 ⊗ σ2 : Tr σ1 = Tr σ2 = 0}.
We can apply the MM ∗ -estimate into this subspace, and this is enough.
Guillaume Aubrun (Lyon)
Random quantum states
Roscoff, June 2012
16 / 18
Applying the MM ∗ -estimate
We cannot apply the MM ∗ -estimate to S directly because
1
S is not symmetric,
2
Its group of isometries does not act irreducibly.
We introduce another symmetrization, Ssym = S ∩ −S. It is checked that
S and Ssym have comparable volumes, mean widths and dual mean widths.
For the second point, the isometry group contains “local unitaries”
ρ 7→ (U ⊗ V )ρ(U ∗ ⊗ V ∗ )
and has an irreducible subspace of small codimension: the space
span{σ1 ⊗ σ2 : Tr σ1 = Tr σ2 = 0}.
We can apply the MM ∗ -estimate into this subspace, and this is enough.
Guillaume Aubrun (Lyon)
Random quantum states
Roscoff, June 2012
16 / 18
The case of N qubits
A harder problem: (C2 )⊗N instead of (Cd )⊗2 .
Definition
A state ρ on (C2 )⊗N is a product state if it is of the form ρ1 ⊗ · · · ⊗ ρN ,
where (ρi ) are states on C2 . A state is separable if it is a convex
combination of product states, otherwise it is entangled.
Question
Consider a random state on (C2 )⊗N with distribution µ2N ,s . It is typically
separable or entangled ? What is the threshold ?
The volume and mean width of S can be estimated (not so directly, one
has to put our convex body in the Löwner position). However, the same
approach fails at the very last step, because the set of separable states on
N qubits has “few” symmetries, and one cannot apply efficiently the
MM ∗ -estimate.
Guillaume Aubrun (Lyon)
Random quantum states
Roscoff, June 2012
17 / 18
The case of N qubits
A harder problem: (C2 )⊗N instead of (Cd )⊗2 .
Definition
A state ρ on (C2 )⊗N is a product state if it is of the form ρ1 ⊗ · · · ⊗ ρN ,
where (ρi ) are states on C2 . A state is separable if it is a convex
combination of product states, otherwise it is entangled.
Question
Consider a random state on (C2 )⊗N with distribution µ2N ,s . It is typically
separable or entangled ? What is the threshold ?
The volume and mean width of S can be estimated (not so directly, one
has to put our convex body in the Löwner position). However, the same
approach fails at the very last step, because the set of separable states on
N qubits has “few” symmetries, and one cannot apply efficiently the
MM ∗ -estimate.
Guillaume Aubrun (Lyon)
Random quantum states
Roscoff, June 2012
17 / 18
The case of N qubits
A harder problem: (C2 )⊗N instead of (Cd )⊗2 .
Definition
A state ρ on (C2 )⊗N is a product state if it is of the form ρ1 ⊗ · · · ⊗ ρN ,
where (ρi ) are states on C2 . A state is separable if it is a convex
combination of product states, otherwise it is entangled.
Question
Consider a random state on (C2 )⊗N with distribution µ2N ,s . It is typically
separable or entangled ? What is the threshold ?
The volume and mean width of S can be estimated (not so directly, one
has to put our convex body in the Löwner position). However, the same
approach fails at the very last step, because the set of separable states on
N qubits has “few” symmetries, and one cannot apply efficiently the
MM ∗ -estimate.
Guillaume Aubrun (Lyon)
Random quantum states
Roscoff, June 2012
17 / 18
The end
THANK YOU !
Guillaume Aubrun (Lyon)
Random quantum states
Roscoff, June 2012
18 / 18