Definitions for Steering, Entanglement, Nonlocality, and the EPR Paradox
DefinitionQAlgorithms_Papers_Independent_Intakethird100_P62The objects the formalized statements of Steering, Entanglement, Nonlocality, and the EPR Paradox (arXiv:quant-ph/0612147) are stated in. Each was checked against the paper by a blind read-back and a moderator.
ProjectiveResolution: PDF p.2: a complete orthogonal spectral resolution, with zero projectors permitted for padding unused outcome labels. This includes degenerate observables.IsResponse: PDF p.2, (3),(5),(6): for each setting, a measurable stochastic kernel on the hidden variable, with no strict positivity requirement. No measurability in the external setting parameter is imposed by the source.IsHiddenEnsemble: PDF p.2: one normalized ensemble of arbitrary Bob density matrices. Integrability excludes the Bochner integral's default value on bad inputs.assemblage: PDF p.2, before (5): unnormalized conditional Bob operator. In particular, a zero-probability outcome is kept as a zero operator, never renormalized.jointProbability: PDF p.2, (2): joint Born outcome probability.HasProjectiveLHS: PDF p.2, (5): one ensemble and one response family precede every setting and outcome. The sample space is arbitrary and may support continuous mixtures.HasHybrid: PDF p.2, (6): hybrid factorization against every Bob projective observable, with the same ensemble and Alice kernel for all Bob settings and outcomes.HasBellLocal: PDF p.2, (3): one local hidden-variable model for all local projective settings. Each response sees only its own setting and the shared variable.Steerable: PDF pp.2--3: state-level Alice-to-Bob steering, allowing all projective observables. A model for all settings also models every restricted strategy.BellNonlocal: PDF p.2: state-level Bell nonlocality for single-copy projective measurements.wernerEta: PDF p.3, (8): convert the paper's eta to the antisymmetric weight p=(d-1+(d+1)eta)/(2d) of the frozen Werner family. For d>=2 this is exactly Eq.(8), including the physical negative-eta interval. The frozen family uses a different parameter, so it cannot be applied directly at eta.oneProjective: Sanity: the one-dimensional projective measurement has a single certain outcome.
import Mathlib
import Definitions.Def_QAlgorithms_Defs_BipartiteSeparability
import Definitions.Def_QAlgorithms_Papers_Independent_Information100_P60
import Definitions.Def_QAlgorithms_Papers_Independent_Information100_P04
/-!
Wiseman--Jones--Doherty, arXiv:quant-ph/0612147v3, PDF pp.1--4.
Finite-dimensional, single-copy, local projective measurements only. Steering is
from Alice to Bob. All hidden-variable integrals below use genuine probability
measures; hidden states may be mixed. The general-Hilbert-space hierarchy and hybrid equivalence, and Gaussian
Theorem 2, are explicitly deferred in the unit's HARD.md. The finite model
below is used to state only the Werner and isotropic steering thresholds.
A d-dimensional projective measurement has at most d nonzero spectral
projectors. Fin d labels accommodate every such measurement by zero padding;
relabeling the distinct eigenvalues does not change either hidden-model test.
-/
noncomputable section
open MeasureTheory
open scoped BigOperators Matrix.Norms.Elementwise
open QAlgorithms.Papers.Independent.Information100.P60 (Mat)
namespace QAlgorithms.Papers.Independent.Intakethird100.P62
attribute [local instance] Classical.propDecidable Classical.decEq
/-- PDF p.2: continuous extended norm for finite-dimensional operator integrals.
The elementwise matrix norm and topology are those of the underlying function space. -/
local instance matrixContinuousENorm (d : ℕ) : ContinuousENorm (Mat (Fin d)) :=
inferInstanceAs (ContinuousENorm (Fin d → Fin d → ℂ))
/-- Sanity: the selected matrix extended norm is continuous in every finite dimension. -/
example (d : ℕ) : Continuous (enorm : Mat (Fin d) → ENNReal) :=
ContinuousENorm.continuous_enorm
/-- PDF p.2: Borel structure for finite-dimensional density-operator ensembles. -/
local instance matrixMeasurable (d : ℕ) : MeasurableSpace (Mat (Fin d)) := borel _
/-- PDF p.2: a complete orthogonal spectral resolution, with zero projectors
permitted for padding unused outcome labels. This includes degenerate observables. -/
structure ProjectiveResolution (d : ℕ) where
projector : Fin d → Mat (Fin d)
hermitian : ∀ a, (projector a).IsHermitian
idempotent : ∀ a, projector a * projector a = projector a
orthogonal : ∀ a b, a ≠ b → projector a * projector b = 0
complete : ∑ a, projector a = 1
/-- PDF p.2, (3),(5),(6): for each setting, a measurable stochastic kernel
on the hidden variable, with no strict positivity requirement. No measurability
in the external setting parameter is imposed by the source. -/
def IsResponse {Omega S : Type} [MeasurableSpace Omega] {d : ℕ}
(r : S → Fin d → Omega → ℝ) : Prop :=
(∀ s a, Measurable (r s a)) ∧
(∀ s a w, 0 ≤ r s a w) ∧
(∀ s w, ∑ a, r s a w = 1)
/-- PDF p.2: one normalized ensemble of arbitrary Bob density matrices.
Integrability excludes the Bochner integral's default value on bad inputs. -/
def IsHiddenEnsemble {Omega : Type} [MeasurableSpace Omega] {b : ℕ}
(mu : Measure Omega) (rho : Omega → Mat (Fin b)) : Prop :=
IsProbabilityMeasure mu ∧ Measurable rho ∧ Integrable rho mu ∧
∀ w, IsDensityMatrix (rho w)
/-- PDF p.2, before (5): unnormalized conditional Bob operator. In particular,
a zero-probability outcome is kept as a zero operator, never renormalized. -/
def assemblage {a b : ℕ} (W : Mat (Fin a × Fin b)) (A : ProjectiveResolution a) (x : Fin a) : Mat (Fin b) :=
traceLeft (W * Matrix.kronecker (A.projector x) (1 : Mat (Fin b)))
/-- PDF p.2, (2): joint Born outcome probability. -/
def jointProbability {a b : ℕ} (W : Mat (Fin a × Fin b))
(A : ProjectiveResolution a) (B : ProjectiveResolution b) (x : Fin a) (y : Fin b) : ℝ :=
effectProb (Matrix.kronecker (A.projector x) (B.projector y)) W
/-- PDF p.2, (5): one ensemble and one response family precede every setting
and outcome. The sample space is arbitrary and may support continuous mixtures. -/
def HasProjectiveLHS {a b : ℕ} (W : Mat (Fin a × Fin b)) (M : Set (ProjectiveResolution a)) : Prop :=
∃ (Omega : Type) (m : MeasurableSpace Omega),
letI : MeasurableSpace Omega := m
∃ (mu : Measure Omega) (rho : Omega → Mat (Fin b))
(r : M → Fin a → Omega → ℝ),
IsHiddenEnsemble mu rho ∧ IsResponse r ∧
(∀ A x, Integrable (fun w => (r A x w : ℂ) • rho w) mu) ∧
∀ A x, assemblage W A.val x = ∫ w, (r A x w : ℂ) • rho w ∂mu
/-- PDF p.2, (6): hybrid factorization against every Bob projective observable,
with the same ensemble and Alice kernel for all Bob settings and outcomes. -/
def HasHybrid {a b : ℕ} (W : Mat (Fin a × Fin b)) (M : Set (ProjectiveResolution a)) : Prop :=
∃ (Omega : Type) (m : MeasurableSpace Omega),
letI : MeasurableSpace Omega := m
∃ (mu : Measure Omega) (rho : Omega → Mat (Fin b))
(r : M → Fin a → Omega → ℝ),
IsHiddenEnsemble mu rho ∧ IsResponse r ∧
(∀ A x (B : ProjectiveResolution b) y,
Integrable (fun w => r A x w * effectProb (B.projector y) (rho w)) mu) ∧
∀ A x (B : ProjectiveResolution b) y, jointProbability W A.val B x y =
∫ w, r A x w * effectProb (B.projector y) (rho w) ∂mu
/-- PDF p.2, (3): one local hidden-variable model for all local projective
settings. Each response sees only its own setting and the shared variable. -/
def HasBellLocal {a b : ℕ} (W : Mat (Fin a × Fin b)) : Prop :=
∃ (Omega : Type) (m : MeasurableSpace Omega),
letI : MeasurableSpace Omega := m
∃ (mu : Measure Omega)
(r : ProjectiveResolution a → Fin a → Omega → ℝ)
(s : ProjectiveResolution b → Fin b → Omega → ℝ),
IsProbabilityMeasure mu ∧ IsResponse r ∧ IsResponse s ∧
(∀ A B x y, Integrable (fun w => r A x w * s B y w) mu) ∧
∀ A B x y, jointProbability W A B x y = ∫ w, r A x w * s B y w ∂mu
/-- PDF pp.2--3: state-level Alice-to-Bob steering, allowing all projective
observables. A model for all settings also models every restricted strategy. -/
def Steerable {a b : ℕ} (W : Mat (Fin a × Fin b)) : Prop :=
¬ HasProjectiveLHS W Set.univ
/-- PDF p.2: state-level Bell nonlocality for single-copy projective measurements. -/
def BellNonlocal {a b : ℕ} (W : Mat (Fin a × Fin b)) : Prop := ¬ HasBellLocal W
/-- PDF p.3, (8): convert the paper's eta to the antisymmetric weight
p=(d-1+(d+1)eta)/(2d) of the frozen Werner family. For d>=2 this is
exactly Eq.(8), including the physical negative-eta interval. The frozen
family uses a different parameter, so it cannot be applied directly at eta. -/
def wernerEta (d : ℕ) (eta : ℝ) : Mat (Fin d × Fin d) :=
QAlgorithms.Papers.Independent.Information100.P04.werner d
(((d : ℝ) - 1 + ((d : ℝ) + 1) * eta) / (2 * (d : ℝ)))
/-- Sanity: the one-dimensional projective measurement has a single certain outcome. -/
def oneProjective : ProjectiveResolution 1 where
projector := fun _ => 1
hermitian := by intro x; simp
idempotent := by intro x; simp
orthogonal := by
intro x y h
exact (h (Subsingleton.elim x y)).elim
complete := by simp
/-- Sanity: stochastic kernels can have constant deterministic responses. -/
example {Omega S : Type} [MeasurableSpace Omega] :
IsResponse (fun (_ : S) (_ : Fin 1) (_ : Omega) => (1 : ℝ)) := by
refine ⟨?_, ?_, ?_⟩
· intro s x; exact measurable_const
· intro s x w; norm_num
· intro s w; simp
/-- Sanity: impossible outcomes never invoke a normalized conditional state. -/
example {a b : ℕ} (A : ProjectiveResolution a) (x : Fin a) :
assemblage (0 : Mat (Fin a × Fin b)) A x = 0 := by
ext i j
simp [assemblage, traceLeft]
/-- Sanity: eta=0 in (8) is the maximally mixed product density. -/
example : wernerEta 2 0 = (1 / 4 : ℂ) • (1 : Mat (Fin 2 × Fin 2)) := by
norm_num [wernerEta, QAlgorithms.Papers.Independent.Information100.P04.werner,
QAlgorithms.Papers.Independent.Information100.P04.wernerP1] <;> module
/-- Sanity: the two-qubit eta=1 endpoint has antisymmetric weight one. -/
example : wernerEta 2 1 =
QAlgorithms.Papers.Independent.Information100.P04.werner 2 1 := by
norm_num [wernerEta]
/-- Sanity: at d=2 the finite harmonic number is 3/2. -/
example : (harmonic 2 : ℝ) = 3 / 2 := by norm_num [harmonic, Finset.sum_range_succ]
end QAlgorithms.Papers.Independent.Intakethird100.P62
Read-back
What the Lean code literally says, in plain math · gpt-6.1-sol
QAlgorithms.Papers.Independent.Intakethird100.P62.BellNonlocal
For any natural numbers and any complex matrix on , being Bell nonlocal means that there do not exist a measurable space , a probability measure , and real responses for all projective resolutions , , such that each response is measurable in , all responses are pointwise nonnegative, their outcome sums are pointwise , and, for every , the product is integrable and
A projective resolution in dimension consists of Hermitian idempotent matrices with pairwise zero products for distinct indices and sum ; individual projectors may be zero. No density condition on is required. The definition includes zero dimensions, where such responses cannot be normalized on a probability space, so the negated existence condition holds.
QAlgorithms.Papers.Independent.Intakethird100.P62.HasBellLocal
For any natural numbers and any complex matrix on , this condition means that there exist a type with a measurable-space structure, a probability measure , and real response functions and defined for every projective resolution and , respectively. For every fixed resolution and outcome, the response as a function of must be measurable; all responses must be pointwise nonnegative, and and for every resolution and every . For all , the product response must be integrable and satisfy
Each projective resolution in dimension consists of Hermitian idempotent matrices with distinct members multiplying to zero and sum ; zero projectors are allowed. No hidden matrix ensemble, measurability in the resolution, or density assumption on is required. If either or , the condition is impossible: the empty projective resolution exists, but its response normalization cannot hold on a probability space.
QAlgorithms.Papers.Independent.Intakethird100.P62.HasHybrid
For any natural numbers , complex matrix on , and set of projective resolutions on , this condition means that there exist a measurable space , a probability measure , a measurable Bochner-integrable function , and real functions , such that every is positive semidefinite with trace , every is measurable and pointwise nonnegative, and for every . For every , every projective resolution on , and every , the scalar function must be integrable, and
A projective resolution in dimension consists of Hermitian idempotent matrices with pairwise zero products for distinct indices and sum ; zero projectors are allowed. No density assumption on is imposed. If is empty, only the hidden-ensemble requirements remain; is impossible, and with is impossible because response normalization would require an empty sum to equal .
QAlgorithms.Papers.Independent.Intakethird100.P62.HasProjectiveLHS
For any natural numbers , complex matrix on , and set of projective resolutions on , this condition means that there exist a type , a measurable-space structure on it, a probability measure , a measurable Bochner-integrable matrix function , and real functions for , such that every is positive semidefinite with trace , every is measurable, , , and, for every , the matrix function is integrable and
Here each supplies Hermitian idempotent matrices , with pairwise zero products for distinct indices and sum , and . All response and density conditions are pointwise; no density condition on , or measurability in , is required. An empty removes the response and assemblage requirements but retains the hidden-ensemble requirements. In particular, is impossible; if and , response normalization is impossible.
QAlgorithms.Papers.Independent.Intakethird100.P62.IsHiddenEnsemble
For any measurable space , natural number , measure on , and function assigning a complex matrix to each , the hidden-ensemble condition means that is a probability measure, that is measurable and Bochner integrable with respect to , and that every is positive semidefinite with trace . The density-matrix condition is required at every point, including points of measure zero. The condition cannot hold when is empty, because no probability measure exists there, or when , because a zero-dimensional matrix has trace .
QAlgorithms.Papers.Independent.Intakethird100.P62.IsResponse
For any types and , any specified measurable-space structure on , any natural number , and any function , the response condition means that is measurable for every , that for every , and that for every . No measurable-space structure or measurability in is required. These requirements hold pointwise, rather than merely almost everywhere. If or is empty, all requirements are vacuous; if and both and are nonempty, the normalization requirement is impossible.
QAlgorithms.Papers.Independent.Intakethird100.P62.ProjectiveResolution
For every natural number , a projective resolution consists of complex matrices , indexed by , satisfying , , whenever , and . No rank or nonzero condition is imposed on any individual projector. The case is included: there are no projectors, and completeness holds because the empty sum and the identity matrix on the zero-dimensional space are the same empty matrix.
QAlgorithms.Papers.Independent.Intakethird100.P62.Steerable
For any natural numbers and any complex matrix on , being steerable means that there do not exist a measurable space , a probability measure , a measurable Bochner-integrable matrix function whose values are all positive semidefinite with trace , and responses for every projective resolution on , such that each is measurable, all responses are pointwise nonnegative, , every is integrable, and
for every . Here a projective resolution is a family of Hermitian idempotent matrices with pairwise zero products for distinct indices and sum , and . No density assumption on is part of this definition. Zero dimensions are included: if , response normalization is impossible, and if , the hidden density matrices are impossible; hence the negated existence condition holds in either case.
QAlgorithms.Papers.Independent.Intakethird100.P62.assemblage
For any natural numbers , any complex matrix indexed by , any family of Hermitian idempotent matrices with for and , and any outcome , the assemblage is the matrix , where . Explicitly, its -entry is . No positivity, Hermiticity, or trace normalization of is assumed. When , there is no possible outcome ; is allowed and gives an empty output matrix.
QAlgorithms.Papers.Independent.Intakethird100.P62.jointProbability
For any natural numbers , any complex matrix on the index set , projective resolutions and , and outcomes , the joint-probability expression is . Here each resolution is a family of Hermitian idempotent matrices, with distinct members multiplying to zero and with sum equal to the identity; zero projectors are permitted. No density-matrix assumption is imposed on , so this definition alone imposes no probability bounds on the resulting real number. If either dimension is zero, the corresponding outcome cannot be supplied.
QAlgorithms.Papers.Independent.Intakethird100.P62.lhs_correlation_equivalence
For every pair of natural numbers with and , every positive semidefinite matrix of trace on , and every set of projective measurements on , the following two existence statements are equivalent. A projective measurement on here consists of exactly Hermitian idempotents , with for and ; zero projectors are allowed. The first existence statement asks for a measurable space , a probability measure , a measurable and -integrable family of density matrices on , and measurable real functions for , with and for every , such that every matrix-valued function is integrable and
where . The second existence statement asks for a measurable space, probability measure, density-matrix family, and response functions with those same measurability, integrability of the density-matrix family, nonnegativity, and normalization requirements, such that, for every , every , every projective measurement on , and every , the scalar function is integrable and
The witnesses in the two existence statements need not be identical. The set may be empty, in which case all requirements indexed by are vacuous.
QAlgorithms.Papers.Independent.Intakethird100.P62.oneProjective
This defines the projective resolution in dimension whose sole projector is the identity matrix . It is Hermitian and idempotent, and its sum is the identity. The requirement that projectors with distinct indices multiply to zero is vacuous because the outcome index set has only one member.
QAlgorithms.Papers.Independent.Intakethird100.P62.strict_hierarchy
For positive subsystem dimensions , call a matrix a density matrix when it is positive semidefinite and has trace . A projective measurement on means a family of Hermitian idempotent matrices satisfying for and ; zero projectors are allowed. Say that has a projective hidden-state model if there exist a measurable space , a probability measure , a measurable, -integrable family of density matrices on , and response functions , indexed by every such projective measurement , which are measurable, nonnegative, and satisfy for every , such that each is integrable and
Here . “Steerable” means that no such model exists. Say that is Bell local if there exist a measurable probability space and measurable, nonnegative response functions , normalized pointwise over their respective outcomes, for all projective measurements on the two subsystems, such that every product is integrable and
“Bell nonlocal” means the absence of this model. Say that is separable positive if for some natural number , nonnegative real weights , and density matrices ; this definition imposes no separate normalization on the weights and allows an empty sum. The declaration asserts the conjunction of three statements: for every and every density matrix on , Bell nonlocality implies steerability, and steerability implies that is not separable positive; there exists a real with for which , where , is a density matrix that is not separable positive and is not steerable; and there exists a real with for which this same formula gives a density matrix that is steerable and Bell local. The universal statement includes or , corresponding to a one-dimensional subsystem.
QAlgorithms.Papers.Independent.Intakethird100.P62.wernerEta
For every natural number and every real number , set , with the arithmetic in this expression taking place in . Let be the complex matrix indexed by pairs from , with when and , and otherwise, and set , . The defined matrix is
The real coefficients are regarded as complex scalars, and in the displayed denominator is real subtraction. There is no restriction on and no density-matrix requirement. All divisions use the total convention ; thus gives the empty matrix, and gives the matrix with entry .
Confirmed by the mission captain (proposal self-audit).
Confirmed by the moderator at approval.