
2
sideration of quantum repeaters [33]. Thus, they answer
one corresponding open question in the review paper on
nonlocality in quantum networks [22]. We can also ap-
ply our methods to a part of the network instead of the
whole, which fits the sparse structure of the network in
the NISQ era and reduces the difficulty of computation.
GHZ state under decoherence.— As a warming-up ex-
ercise we discuss the Greenberger-Horne-Zeilinger (GHZ)
state of nqubits under decoherence,
ρ(α) = (1 −α)|GHZ+⟩⟨GHZ+|+α|GHZ−⟩⟨GHZ−|,(1)
where |GHZ±⟩= (|0···0⟩±|1···1⟩)/√2, and α∈[0,1/2]
describes the degree of decoherence. Despite its simplic-
ity, this example allows us to introduce our main ideas.
If all the nparties implement the same measurement
Z=|0⟩⟨0|−|1⟩⟨1|, then two possible combinations of
outcomes happen equally with probability 1/2, i.e., ei-
ther all of the outcomes are 0, or all of them are 1.
To simulate this statistical behaviour without genuine n-
partite entanglement, the state for simulation can only be
ρc= [|0···0⟩⟨0···0|+|1···1⟩⟨1···1|]/2, since no other
0−1 string appear as a combination of outcomes. Such
a simulation invalidates known methods with only sta-
tistical data [19, 26–28]. It costs at least one classical
bit of randomness, as the Shannon entropy or the Von
Neumann entropy of ρcis 1. However, the Von Neumann
entropy of the state ρ(α) is −[αlog α+(1−α) log(1−α)],
which is strictly less than 1 for α∈[0,1/2). This means
that we cannot simulate the statistical behaviour and the
Von Neumann entropy of ρ(α) simultaneously by a quan-
tum network with at most (n−1)-partite sources.
The Von Neumann entropy is one way to measure the
purity of the state, capturing partially the classical cor-
relations in the state. To continue, we examine firstly
different measures of purity and choose a suitable one
for our following methods. For a given state ρin the
d-dimensional space, the common measures of its pu-
rity [40–42] are R´enyi α-purity log2d−log2(Tr(ρα))/(1−
α), which converges to the Von Neumann entropy as
αtends to 1, and linear entropy purity Tr(ρ2)−1/d.
Through the whole text, we take τ= Tr(ρ2) to quantify
the purity, which determines R´enyi 2-purity and linear
entropy purity. The advantage of τover other quanti-
fiers, like the Von Neumann entropy, is that it fits the
covariance of experimental data well in our approach, as
both of them contain the information of ρ2. As for the
estimation of purity of a multipartite state with different
measures, it can be done efficiently with only local opera-
tions [38, 39], which are feasible in the network scenario.
Noisy quantum networks.— Noise is unavoidable for
the quantum network states in the NISQ era, either the
local noise or the global one. Quantum networks with dif-
ferent noise models can all be classified as CQNs. Firstly,
we develop the covariance matrix decomposition method
for CQN, where a key step is to separate the part related
to global classical correlation out in the whole covariance
matrix. For a given hypergraph G(V, E) and a state ρ
from CQN of G, the state ρcan be decomposed as
ρ=Xkpkρk, ρk=Oi∈VC(k)
iOe∈Eη(k)
e,(2)
where {pk}kwith Pkpk= 1 and pk>0 is the global clas-
sical correlation, C(k)
iis a local channel for the i-th party,
η(k)
eis an entangled state distributed from the source la-
beled by the hyperedge e.
For simplicity, we assume each party has only one mea-
surement, and denote Mithe measurement for the i-th
party. Then we introduce three kinds of covariance matri-
ces, Γ, Γ(k)and Γ(c), whose elements in the i-th row and
j-th column are Γij , Γ(k)
ij and Γ(c)
ij , respectively, where
Γij =⟨MiMj⟩−⟨Mi⟩⟨Mj⟩,⟨Mi⟩= Tr(ρMi),
Γ(k)
ij =⟨MiMj⟩k− ⟨Mi⟩k⟨Mj⟩k,⟨Mi⟩k= Tr(ρkMi),
Γ(c)
ij =Xkpk⟨Mi⟩k⟨Mj⟩k− ⟨Mi⟩⟨Mj⟩.(3)
The covariance matrix Γ is the one that can be observed
directly in experiments. The covariance matrices {Γ(k)}k
are hidden in the experimental data when we assume that
the randomness of the sampling {pk, ρk}kis inaccessible.
The covariance matrix Γ(c)can be viewed as a classical
covariance matrix, since it is only about the distribution
of classical data {⟨M1⟩k,...,⟨Mn⟩k}k. Throughout the
whole paper, we only consider the dichotomic measure-
ments with outcomes ±1. A pivotal observation is that
the classical covariance matrix Γ(c)can be separated out
from the observed one Γ perfectly, i.e.,
Γ = XkpkΓ(k)+ Γ(c),(4)
whose proof can be found in Sec. A in Supplemental Ma-
terial (SM) [43]. Since {Γ(k)}kare about network states
from IQN, the existing method in Ref. [19] can be em-
ployed to impose constraints on them. However, if there
is no limitation of Γ(c), the observed covariance matrix Γ
can still have arbitrary relation with the network topol-
ogy G. As it turns out, the purity of the state implies
a nontrivial condition on Γ(c), leading to a semi-definite
programming (SDP) to determine whether a state can
arise from CQN with a given topology.
Observation 1. For a given state ρfrom the CQN with
the network topology G(V, E), measurements {Mi}i∈V,
which result in the covariance matrix Γ, it holds that
Γ = Xe∈EΥe+T, ΠeΥeΠe= Υe⪰0,
T⪰0,maxi∈VTii ≤β, Tr(T)≤l1β, (5)
where l1is the maximal eigenvalue of Pi∈VMi⊗Mi,
β= 2√1−τ2,Tii is the i-th diagonal term of T,Πe=
Pi∈ePiwith Pito be the projection onto i-th row.