This paper was converted on www.awesomepapers.org from LaTeX by an anonymous user.
Want to know more? Visit the Converter page.

Generating Randomness from a Computable, Non-random Sequence of Qubits

Tejas Bhojraj Department of Mathematics
University of Wisconsin-Madison WI,USA bhojraj@wisc.edu
       
Abstract

Nies and Scholz [10] introduced the notion of a state to describe an infinite sequence of qubits and defined quantum-Martin-Löf randomness for states, analogously to the well known concept of Martin-Löf randomness for elements of Cantor space (the space of infinite sequences of bits). We formalize how ‘measurement’ of a state in a basis induces a probability measure on Cantor space. A state is ‘measurement random’ (mR) if the measure induced by it, under any computable basis, assigns probability one to the set of Martin-Löf randoms. Equivalently, a state is mR if and only if measuring it in any computable basis yields a Martin-Löf random with probability one. While quantum-Martin-Löf random states are mR, the converse fails: there is a mR state, ρ\rho which is not quantum-Martin-Löf random. In fact, something stronger is true. While ρ\rho is computable and can be easily constructed, measuring it in any computable basis yields an arithmetically random sequence with probability one. I.e., classical arithmetic randomness can be generated from a computable, non-quantum random sequence of qubits.

1 Introduction

Martin-Löf randomness for infinite sequences of bits, a central concept in algorithmic randomness[9], has recently been generalized to infinite sequences of qubits [10]. As is standard in the literature in algorithmic randomness and mathematical logic, let 2ω2^{\omega} denote the collection of infinite sequences of bits, let 2n2^{n} denote the set of bit strings of length nn, 2<ω:=n2n2^{<\omega}:=\bigcup_{n}2^{n} and let 2ω:=2<ω2ω2^{\leq\omega}:=2^{<\omega}\cup 2^{\omega}. X2ωX\in 2^{\omega} is said to be Martin-Löf random (MLR) if it is not in any effectively null set with respect to the uniform measure on 2ω2^{\omega} [9]. Nies and Scholz introduced the notion of a state to describe an infinite sequence of qubits [10] and then defined quantum-Martin-Löf randomness for states. We quickly sketch those parts of their work most relevant to the arguments in this paper. All definitions in the introduction section are from [10] or [9]. See [9] and [7] for a detailed introduction to Martin-Löf randomness and computability theory.

Definition 1.1.

A state, ρ=(ρn)n\rho=(\rho_{n})_{n\in\mathbb{N}} is an infinite sequence of density matrices such that ρn2n×2n\rho_{n}\in\mathbb{C}^{2^{n}\times 2^{n}} and n\forall n, PT2(ρn)=ρn1PT_{\mathbb{C}^{2}}(\rho_{n})=\rho_{n-1}.

Here, PT2PT_{\mathbb{C}^{2}} denotes the partial trace which ‘traces out’ the last qubit from 2n\mathbb{C}^{2^{n}}. ρ\rho is a infinite sequence of qubits whose first nn qubits are given by ρn\rho_{n}. For this notion to be meaningful, ρ\rho is required to be coherent in the following sense; for all nn, ρn\rho_{n}, when ‘restricted’ via the partial trace to it’s first n1n-1 qubits, has the same measurement statistics as the state on n1n-1 qubits given by ρn1\rho_{n-1}. We now sketch a few notions from computability theory pertinent to us. Let A2ωA\in 2^{\omega}. We define an AA-computable function to be a total function that can be realized by a Turing machine with AA as an oracle. By ‘computable’, we will refer to \emptyset-computable. The concept of an AA-computable sequence of natural numbers will come up frequently in our discussion.

Definition 1.2.

A sequence (an)n(a_{n})_{n\in\mathbb{N}} is said to be AA-computable if there is a AA-computable function φ\varphi such that φ(n)=an\varphi(n)=a_{n}

Definition 1.3.

A special projection is a hermitian projection matrix with complex algebraic entries.

Since the algebraic complex numbers have a computable presentation (see [10]), we may identify a special projection with a natural number and hence talk about computable sequences of special projections. Let II denote the two by two identity matrix.

Definition 1.4.

A quantum Σ10\Sigma_{1}^{0} set (or q-Σ10\Sigma_{1}^{0} set for short) G is a computable sequence of special projections G=(pi)iG=(p_{i})_{i\in\mathbb{N}} such that pip_{i} is 2i2^{i} by 2i2^{i} and range(piI)(p_{i}\otimes I)\subseteq range (pi+1),i(p_{i+1}),\forall i\in\mathbb{N}.

While a 2n2^{n} by 2n2^{n} special projection may be thought of as a computable projective measurement on a system of nn qubits, a q-Σ10\Sigma_{1}^{0} class corresponds to a computable sequence of projective measurements on longer and longer systems of qubits and mirrors the concept of a Σ10\Sigma_{1}^{0} class in computability theory. One of the many equivalent ways of defining a Σ10\Sigma_{1}^{0} class is as follows. If C2nC\subset 2^{n}, let C2ω\llbracket C\rrbracket\subseteq 2^{\omega} be the set of all X2ωX\in 2^{\omega} such that the initial segment of XX of length nn is in CC.

Definition 1.5.

A Σ10\Sigma_{1}^{0} class S2ωS\subseteq 2^{\omega} is any set of the form,

S=iAiS=\bigcup_{i\in\mathbb{N}}\llbracket A_{i}\rrbracket

where

  1. 1.

    Ai2i,iA_{i}\subseteq 2^{i},\forall i\in\mathbb{N}

  2. 2.

    The indices of AiA_{i} form a computable sequence. (Being a finite set, each AiA_{i} has a natural number coding it.

  3. 3.

    AiAi+1,i\llbracket A_{i}\rrbracket\subseteq\llbracket A_{i+1}\rrbracket,\forall i\in\mathbb{N}

A Σ10\Sigma_{1}^{0} class, S is coded (non-uniquely) by the index of the total computable function generating the sequence (Ai)i(A_{i})_{i\in\mathbb{N}} occurring in (2) in the definition of SS. Hence, the notion of a computable sequence of Σ10\Sigma_{1}^{0} classes makes sense (see [9], section 3.2). One sees that the special projections qiq_{i} in the definition of the q-Σ10\Sigma_{1}^{0}, GG play the role of the AiA_{i}s which generate a the Σ10\Sigma_{1}^{0} class, SS. The following notion is a quantum analog of the uniform measure of SS which equals limn(2n|An|)\lim_{n}(2^{-n}|A_{n}|), where |.||.| refers to the cardinality. (The uniform measure on 2ω2^{\omega} is the measure induced by letting the measure of τ\llbracket\tau\rrbracket to be 2|τ|2^{-|\tau|} for each τ2<ω\tau\in 2^{<\omega}. Here, |τ|:=n|\tau|:=n if τ2n\tau\in 2^{n}.)

Definition 1.6.

If G=(pn)nG=(p_{n})_{n\in\mathbb{N}} is a q-Σ10\Sigma_{1}^{0} class, define τ(G):=limn(2n|qn|)\tau(G):=\lim_{n}(2^{-n}|q_{n}|) where, |qn||q_{n}| is the rank of qnq_{n}.

Informally, a q-Σ10\Sigma_{1}^{0} class, G=(pn)nG=(p_{n})_{n\in\mathbb{N}} may be thought of as an observable whose expected value, when ‘measured’ on a state ρ=(ρn)n\rho=(\rho_{n})_{n\in\mathbb{N}} is ρ(G):=limn\rho(G):=\lim_{n} tr (ρnpn)(\rho_{n}p_{n}). The positive-semidefiniteness of density matrices, the coherence of the components of ρ\rho and condition (3) in definition 1.5 ensure that tr(ρnpn)tr(\rho_{n}p_{n}) is non decreasing in nn. As tr(ρnpn)1tr(\rho_{n}p_{n})\leq 1 for all nn, ρ(G):=limn\rho(G):=\lim_{n} tr (ρnpn)(\rho_{n}p_{n}) exists.

Definition 1.7.

A classical Martin-Löf test (MLT) is a computable sequence, (Sm)m(S_{m})_{m\in\mathbb{N}} of Σ01\Sigma_{0}^{1} classes such that the uniform measure of SmS_{m} is less than or equal to 2m2^{-m} for all m.

Definition 1.8.

A quantum Martin-Löf test (q-MLT) is a computable sequence, (Sm)m(S_{m})_{m\in\mathbb{N}} of q-Σ01\Sigma_{0}^{1} classes such that τ(Sm)\tau(S_{m}) is less than or equal to 2m2^{-m} for all m.

Having established the necessary prerequisites, we can define a quantum Martin-Löf random (q-MLR) state. Roughly speaking, a state is q-MLR if it cannot be ‘detected by projective measurements of arbitrarily small rank’.

Definition 1.9.

ρ\rho is q-MLR if for any q-MLT (Sm)m(S_{m})_{m\in\mathbb{N}}, infρm(Sm)=0{}_{m\in\mathbb{N}}\rho(S_{m})=0.

Definition 1.10.

ρ\rho is said to fail the q-MLT (Sm)m(S_{m})_{m\in\mathbb{N}}, at order δ\delta, if infρm(Sm)>δ{}_{m\in\mathbb{N}}\rho(S_{m})>\delta.

2 Measuring a state induces a measure on 2ω2^{\omega}

To fix notation, let X(n)X(n) denote the nnth bit of an X2ωX\in 2^{\leq\omega} , let p(E)p(E) stand for the probability of the event EE.

Definition 2.1.

An A-computable measurement system B=((b0n,b1n))nB=((b^{n}_{0},b^{n}_{1}))_{n\in\mathbb{N}} (or just ‘measurement system’ for short) is a sequence of orthonormal bases for 2\mathbb{C}^{2} such that each binb^{n}_{i} is complex algebraic and the sequence ((b0n,b1n))n((b^{n}_{0},b^{n}_{1}))_{n\in\mathbb{N}} is A-computable.

Let ρ=(ρn)n\rho=(\rho_{n})_{n\in\mathbb{N}} be a state and B=((b0n,b1n))nB=((b^{n}_{0},b^{n}_{1}))_{n\in\mathbb{N}} be a measurement system. We now work towards formalizing a notion of qubitwise measurement of ρ\rho in the bases in BB. A (probability) premeasure [7],pp (also called a measure representation [9]), is a function from the set of all finite bit strings to [0,1][0,1] satisfying n\forall n, τ2n,p(τ)=p(τ0)+p(τ1)\forall\tau\in 2^{n},p(\tau)=p(\tau 0)+p(\tau 1). pp induces a measure on 2ω2^{\omega} which is seen to be unique by Carathéodory’s extension theorem (See 6.12.1 in [7]). Flipping a 0,10,1 sided fair coin repeatedly induces a probability measure (which happens to be the uniform measure) on 2ω2^{\omega} as follows. Let the random variable Z(n)Z(n) denote the outcome of the the nnth coin flip. The sequence (Z(n))n(Z(n))_{n\in\mathbb{N}} induces a premeasure, pp, on 2<ω2^{<\omega} which extends to the uniform measure on 2ω2^{\omega}. Here, p(σ)=2np(\sigma)=2^{-n} is the probability that Z(i)=σ(i)Z(i)=\sigma(i) for all i|σ|i\leq|\sigma|. Similarly the act of measuring ρ\rho qubit by qubit in BB induces a premeasure on 2<ω2^{<\omega} which extends to a probability measure (denoted μρB\mu_{\rho}^{B}) on 2ω2^{\omega} as follows. Let the random variable X(n)X(n) be the 0,10,1 valued outcome of the measurement of the nnth qubit of ρ\rho. Let pp be the premeasure induced by the sequence (X(n))n(X(n))_{n\in\mathbb{N}} on 2<ω2^{<\omega}. pp extends to μρB\mu_{\rho}^{B} on 2ω2^{\omega}. For any A2ωA\subseteq 2^{\omega}, μρB(A)\mu_{\rho}^{B}(A) is the probability that XAX\in A where XX is the element of 2ω2^{\omega} obtained in the limit by the qubit by qubit measurement of ρ\rho in BB. The most conspicuous difference between the two situations is that while the (Z(n))n(Z(n))_{n\in\mathbb{N}} are independent, (X(n))n(X(n))_{n\in\mathbb{N}} need not be independent as the elements of ρ\rho can be entangled. We now formalize the above. The following calculations follow from standard results mentioned, for example, in [6].

We now define (X(n))n(X(n))_{n\in\mathbb{N}} and pp, the induced premeasure. Measure ρ1\rho_{1} by the measurement operators {|b01><b01|,|b11><b11|}\{|b^{1}_{0}\big{>}\big{<}b^{1}_{0}|,|b^{1}_{1}\big{>}\big{<}b^{1}_{1}|\} and define X(1):=iX(1):=i if bi1b^{1}_{i} was obtained by the above measurement. Let ρ2^\hat{\rho_{2}} be the density matrix corresponding to the post-measurement state of ρ2\rho_{2} given that ρ2\rho_{2} yields |bX(1)1><bX(1)1|I|b^{1}_{X(1)}\big{>}\big{<}b^{1}_{X(1)}|\otimes I if measured in the system

(|bi1><bi1|I)i{0,1}.(|b^{1}_{i}\big{>}\big{<}b^{1}_{i}|\otimes I)_{i\in\{0,1\}}.

I.e,

ρ2^=(|bX(1)1><bX(1)1|I)ρ2(|bX(1)1><bX(1)1|I)tr((|bX(1)1><bX(1)1|I)ρ2).\hat{\rho_{2}}=\dfrac{(|b^{1}_{X(1)}\big{>}\big{<}b^{1}_{X(1)}|\otimes I)\rho_{2}(|b^{1}_{X(1)}\big{>}\big{<}b^{1}_{X(1)}|\otimes I)}{tr((|b^{1}_{X(1)}\big{>}\big{<}b^{1}_{X(1)}|\otimes I)\rho_{2}\big{)}}.

To define X(2)X(2), measure ρ2^\hat{\rho_{2}} by the measurement operators

(I|bi2><bi2|)i{0,1},(I\otimes|b^{2}_{i}\big{>}\big{<}b^{2}_{i}|)_{i\in\{0,1\}},

and set X(2):=iX(2):=i if I|bi2><bi2|I\otimes|b^{2}_{i}\big{>}\big{<}b^{2}_{i}| is obtained. We use ρ2^\hat{\rho_{2}} instead of ρ2\rho_{2} to define X(2)X(2) to account for the previous measurement of the first qubit. X(n)X(n) is defined similarly. By the above,

p(ij):=p(X(1)=i,X(2)=j)=p(X(1)=i)p(X(2)=j|X(1)=i)=p(ij):=p(X(1)=i,X(2)=j)=p(X(1)=i)p(X(2)=j|X(1)=i)=
p(X(1)=i)tr[I|bj2><bj2|((|bi1><bi1|I)ρ2(|bi1><bi1|I)tr((|bi1><bi1|I)ρ2))].p(X(1)=i)tr\big{[}I\otimes|b^{2}_{j}\big{>}\big{<}b^{2}_{j}|(\dfrac{(|b^{1}_{i}\big{>}\big{<}b^{1}_{i}|\otimes I)\rho_{2}(|b^{1}_{i}\big{>}\big{<}b^{1}_{i}|\otimes I)}{tr((|b^{1}_{i}\big{>}\big{<}b^{1}_{i}|\otimes I)\rho_{2})})\big{]}.

Since PT2(ρ2)=ρ1PT_{\mathbb{C}^{2}}(\rho_{2})=\rho_{1}, p(X(1)=i)=tr((|bi1><bi1|I)ρ2)p(X(1)=i)=tr((|b^{1}_{i}\big{>}\big{<}b^{1}_{i}|\otimes I)\rho_{2}\big{)}. So,

p(ij)=tr[ρ2(|bi1bj2><bi1bj2|)].p(ij)=tr\big{[}\rho_{2}(|b^{1}_{i}b^{2}_{j}\big{>}\big{<}b^{1}_{i}b^{2}_{j}|)\big{]}.

Given τ2n\tau\in 2^{n}, similar calculations show that

p(τ):=p(X(1)=τ(1),,X(n)=τ(n))=tr[ρn(|i=1nbτ(i)i><i=1nbτ(i)i|].\displaystyle p(\tau):=p(X(1)=\tau(1),\dots,X(n)=\tau(n))=tr\big{[}\rho_{n}(|\bigotimes_{i=1}^{n}b^{i}_{\tau(i)}\big{>}\big{<}\bigotimes_{i=1}^{n}b^{i}_{\tau(i)}\big{|}\big{]}. (2.1)

This defines pp. The following lemma shows that p(.)p(.) is a premeasure. Define μρB\mu^{B}_{\rho} to be the unique probability measure induced by it.

Lemma 2.2.

n\forall n, τ2n,p(τ)=p(τ0)+p(τ1)\forall\tau\in 2^{n},p(\tau)=p(\tau 0)+p(\tau 1)

Proof.

Noting that for j{0,1}j\in\{0,1\},

ρn+1(|i=1nbτ(i)ibjn+1><i=1nbτ(i)ibjn+1|)=ρn+1(|i=1nbτ(i)i><i=1nbτ(i)i||bjn+1><bjn+1|),\rho_{n+1}(|\bigotimes_{i=1}^{n}b^{i}_{\tau(i)}\otimes b^{n+1}_{j}\big{>}\big{<}\bigotimes_{i=1}^{n}b^{i}_{\tau(i)}\otimes b^{n+1}_{j}|)=\rho_{n+1}(|\bigotimes_{i=1}^{n}b^{i}_{\tau(i)}\big{>}\big{<}\bigotimes_{i=1}^{n}b^{i}_{\tau(i)}|\otimes|b^{n+1}_{j}\big{>}\big{<}b^{n+1}_{j}|),

and letting A:=|i=1nbτ(i)i><i=1nbτ(i)i|A:=|\bigotimes_{i=1}^{n}b^{i}_{\tau(i)}\big{>}\big{<}\bigotimes_{i=1}^{n}b^{i}_{\tau(i)}|, the right hand side is

=tr[(A|b0n+1><b0n+1|)ρn+1+(A|b1n+1><b1n+1|)ρn+1]=tr\big{[}(A\otimes|b^{n+1}_{0}\big{>}\big{<}b^{n+1}_{0}|)\rho_{n+1}+(A\otimes|b^{n+1}_{1}\big{>}\big{<}b^{n+1}_{1}|)\rho_{n+1}]
=tr[(A(|b0n+1><b0n+1|+|b1n+1><b1n+1|))ρn+1]=tr[(AI)ρn+1]=tr[Aρn]=p(τ)=tr\big{[}(A\otimes(|b^{n+1}_{0}\big{>}\big{<}b^{n+1}_{0}|+|b^{n+1}_{1}\big{>}\big{<}b^{n+1}_{1}|))\rho_{n+1}]=tr\big{[}(A\otimes I)\rho_{n+1}]=tr[A\rho_{n}]=p(\tau)

Remark 2.3.

If BB is SS-computable and ρ\rho is TT-computable, then the sequence {μρB(σ)}σ\{\mu^{B}_{\rho}(\sigma)\}_{\sigma\in\mathbb{N}} is STS\oplus T-computable.

Here, STS\oplus T is obtained by putting SS on the even bits and TT on the odd bits [9].

3 Measurement Randomness

Let MLR2ωMLR\subset 2^{\omega} be the set of MLR bitstrings. If ρ\rho is a state and BB a measurement system, μρB(MLR)\mu^{B}_{\rho}(MLR) is the probability of getting a MLR bitstring by a qubit-wise measurement of ρ\rho as described in the previous section.

Definition 3.1.

ρ\rho is measurement random (mR) if for any computable measurement system, B, μρB(MLR)=1\mu^{B}_{\rho}(MLR)=1

Theorem 3.2.

All q-MLR states are also mR states.

Proof.

Let ρ=(ρn)n\rho=(\rho_{n})_{n\in\mathbb{N}} be q-MLR. Suppose towards a contradiction that there is a δ(0,1)\delta\in(0,1) and a computable B=((b0n,b1n))nB=((b^{n}_{0},b^{n}_{1}))_{n\in\mathbb{N}} such that μρB(2ω/MLR)>δ\mu^{B}_{\rho}(2^{\omega}/MLR)>\delta. Let (Sm)m(S^{m})_{m} be the universal MLT [9] and let for all mm,

Sm=miAim,\displaystyle S^{m}=\bigcup_{m\leq i}\llbracket A^{m}_{i}\rrbracket, (3.1)

where the AimA^{m}_{i}s satisfy the conditions of Definition 1.5. By the definition of a MLT, for all mm and all imi\geq m, we can write Aim={τ1m,i,,τkm,im,i}2iA^{m}_{i}=\{\tau^{m,i}_{1},\dots,\tau^{m,i}_{k^{m,i}}\}\subset 2^{i} for some 0km,i2im0\leq k^{m,i}\leq 2^{i-m}. Now define a q-MLT as follows. For all mm and imi\geq m, let τa=τam,i\tau_{a}=\tau_{a}^{m,i} for convenience and define the special projection:

pim=akm,i(|q=1ibτa(q)q><q=1ibτa(q)q|).\displaystyle p^{m}_{i}=\sum_{a\leq k^{m,i}}(|\bigotimes_{q=1}^{i}b^{q}_{\tau_{a}(q)}\big{>}\big{<}\bigotimes_{q=1}^{i}b^{q}_{\tau_{a}(q)}\big{|}\big{)}. (3.2)

Letting Pm:=(pim)miP^{m}:=(p^{m}_{i})_{m\leq i}, we see that (Pm)m(P^{m})_{m\in\mathbb{N}} is a q-MLT (For each mm, the sequence (pim)mi(p^{m}_{i})_{m\leq i} is computable since BB and (Aim)mi(A^{m}_{i})_{m\leq i} are computable. Condition 3 in Definition 1.5 implies that for all ii, range(pim)(p^{m}_{i})\subseteqrange(pi+1m)(p^{m}_{i+1}). So, PmP^{m} is a q-Σ01\Sigma_{0}^{1} class for all mm. km,i2imk^{m,i}\leq 2^{i-m} for all m,im,i implies that τ(Pm)2m\tau(P^{m})\leq 2^{-m} for all mm. Since (Sm)m(S^{m})_{m\in\mathbb{N}} is a MLT, (Pm)m(P^{m})_{m\in\mathbb{N}} is a computable sequence.) For all m, (2ω/MLR)Sm(2^{\omega}/MLR)\subseteq S^{m} holds by the definition of a universal MLT. Hence, since 3.1 is an increasing union and as μρB(2ω/MLR)>δ\mu^{B}_{\rho}(2^{\omega}/MLR)>\delta, for all mm there exists an i(m)>mi(m)>m such that

μρB(Ai(m)m)>δ.\displaystyle\mu^{B}_{\rho}(\llbracket A^{m}_{i(m)}\rrbracket)>\delta. (3.3)

Fix such an mm and corresponding i=i(m)i=i(m) and let Aim={τ1,,τkm,i}A^{m}_{i}=\{\tau_{1},\dots,\tau_{k^{m,i}}\} for some km,i2imk^{m,i}\leq 2^{i-m} as in 3.2. By 2.1 and 3.3, we have that

δ<akp(τa)=akm,itr[ρi(|q=1ibτa(q)q><q=1ibτa(q)q|)]=tr[ρiakm,i(|q=1ibτa(q)q><q=1ibτa(q)q|)]\displaystyle\delta<\sum_{a\leq k}p(\tau_{a})=\sum_{a\leq k^{m,i}}tr\big{[}\rho_{i}(|\bigotimes_{q=1}^{i}b^{q}_{\tau_{a}(q)}\big{>}\big{<}\bigotimes_{q=1}^{i}b^{q}_{\tau_{a}(q)}\big{|}\big{)}\big{]}=tr\big{[}\rho_{i}\sum_{a\leq k^{m,i}}(|\bigotimes_{q=1}^{i}b^{q}_{\tau_{a}(q)}\big{>}\big{<}\bigotimes_{q=1}^{i}b^{q}_{\tau_{a}(q)}\big{|}\big{)}\big{]} (3.4)

So, by 3.2 and 3.4, we see that for all mm there is an ii such that,

δ<tr[ρipim]ρ(Pm).\delta<tr[\rho_{i}p^{m}_{i}]\leq\rho(P^{m}).

So, inf(ρ(Pm))m>δ{}_{m}(\rho(P^{m}))>\delta, contradicting that ρ\rho is q-MLR. ∎

Definition 3.3.

ρ=(ρn)n\rho=(\rho_{n})_{n\in\mathbb{N}} is computable if the sequence (ρn)n(\rho_{n})_{n\in\mathbb{N}} is computable.

Theorem 3.4.

There is a computable state which is not q-MLR but is mR.

Proof.

All matrices in this proof are in the standard basis. Let ρ=n=5dn\rho=\bigotimes_{n=5}^{\infty}d_{n} and for N>5N>5, SN:=n=5NdnS_{N}:=\bigotimes_{n=5}^{N}d_{n}. where dnd_{n} is a 2n2^{n} by 2n2^{n} matrix with 2n2^{-n} along the diagonal and rn:=2n/nr_{n}:=\lfloor 2^{n}/n\rfloor many 2n2^{-n}s on the extreme ends of the anti-diagonal. Formally, define dnd_{n} to be the symmetric matrix such that: For irni\leq r_{n}, dn(i,j)=2nd_{n}(i,j)=2^{-n} if j=ij=i or j=2ni+1j=2^{n}-i+1 and dn(i,j)=0d_{n}(i,j)=0 otherwise. For rn<i<2nrnr_{n}<i<2^{n}-r_{n}, dn(i,j)=2nd_{n}(i,j)=2^{-n} if j=ij=i and dn(i,j)=0d_{n}(i,j)=0 otherwise. For example, r3=2r_{3}=2 and so,

d3=[2300000023023000023000230000000023000000002300000000230002300002302300000023]d_{3}=\begin{bmatrix}2^{-3}&0&0&0&0&0&0&2^{-3}\\ 0&2^{-3}&0&0&0&0&2^{-3}&0\\ 0&0&2^{-3}&0&0&0&0&0\\ 0&0&0&2^{-3}&0&0&0&0\\ 0&0&0&0&2^{-3}&0&0&0\\ 0&0&0&0&0&2^{-3}&0&0\\ 0&2^{-3}&0&0&0&0&2^{-3}&0\\ 2^{-3}&0&0&0&0&0&0&2^{-3}\\ \end{bmatrix}

Clearly, dnd_{n} is a density matrix. The theorem will be proved via the following lemmas.

Lemma 3.5.

ρ\rho is not q-MLR.

Proof.

It is easy to see that zero has multiplicity rnr_{n} as an eigenvalue of dnd_{n}. Hence, letting qn=2nrnq_{n}=2^{n}-r_{n}, the eigenpairs of dnd_{n} can be listed as {αin,vin}i=12n\{\alpha^{n}_{i},v^{n}_{i}\}_{i=1}^{2^{n}} where αin=0\alpha^{n}_{i}=0 if qn+1i2nq_{n}+1\leq i\leq 2^{n} and (vin)i=12n(v^{n}_{i})_{i=1}^{2^{n}} is a orthonormal basis of 2n\mathbb{C}^{2^{n}}.

Fix a N>5N>5. By properties of the Kronecker product, SNS_{N} has a orthonormal basis of eigenvectors:

{n=5Nvl(n)n:(l(n))n=5N is a sequence such that for all n,l(n)2n},\{\bigotimes_{n=5}^{N}v^{n}_{l(n)}:(l(n))_{n=5}^{N}\text{ is a sequence such that for all }n,l(n)\leq 2^{n}\},

and n=5Nvl(n)n\bigotimes_{n=5}^{N}v^{n}_{l(n)} has eigenvalue n=5Nαl(n)n\prod_{n=5}^{N}\alpha^{n}_{l(n)}. Letting MNM_{N} be those elements of the above eigenbasis having non-zero eigenvalues, we have that

MN={n=5Nvl(n)n:(l(n))n=5N is a sequence such that for all n,l(n)qn}.M_{N}=\{\bigotimes_{n=5}^{N}v^{n}_{l(n)}:(l(n))_{n=5}^{N}\text{ is a sequence such that for all }n,l(n)\leq q_{n}\}.

By the definition of qnq_{n},

|MN|=n=5N2n2n/nn=5N2n(2n/n)+1=n=5N2n(1n1+2n)=n=5N2nn=5N(1n1+2n).|M_{N}|=\prod_{n=5}^{N}2^{n}-\lfloor 2^{n}/n\rfloor\leq\prod_{n=5}^{N}2^{n}-(2^{n}/n)+1=\prod_{n=5}^{N}2^{n}(1-n^{-1}+2^{-n})=\prod_{n=5}^{N}2^{n}\prod_{n=5}^{N}(1-n^{-1}+2^{-n}).

Noting that n=5(1n1+2n)=0\prod_{n=5}^{\infty}(1-n^{-1}+2^{-n})=0, define a q-MLT (Tm)m(T_{m})_{m\in\mathbb{N}} as follows. Given mm, we describe the construction of TmT_{m}. Find N=N(m)N=N(m) such that n=5N(1n1+2n)<2m\prod_{n=5}^{N}(1-n^{-1}+2^{-n})<2^{-m}. Let γ(N):=n=5Nn\gamma(N):=\sum_{n=5}^{N}n and let

pγ(N)=vMN|v><v|.p_{\gamma(N)}=\sum_{v\in M_{N}}|v\big{>}\big{<}v|.

pγ(N)p_{\gamma(N)} is a special projection on 2γ(N)\mathbb{C}^{2^{{\gamma(N)}}} having rank equal to |MN||M_{N}|. Let pk=p_{k}=\emptyset for k<γ(N)k<\gamma(N) and

pk:=pγ(N)i=1kγ(N)Ip_{k}:=p_{\gamma(N)}\otimes\bigotimes_{i=1}^{k-\gamma(N)}I

for k>γ(N)k>\gamma(N). Using that ρ\rho is computable, it is easy to see that (pk)k(p_{k})_{k\in\mathbb{N}} is a q-Σ10\Sigma^{0}_{1} class. Let Tm:=(pk)kT_{m}:=(p_{k})_{k\in\mathbb{N}}. (Tm)m(T_{m})_{m\in\mathbb{N}} is a q-MLT since the choice of N(m)N(m) implies that τ(Tm)<2m\tau(T_{m})<2^{-m} and as N(m)N(m) can be computed from mm. (Tm)m(T_{m})_{m\in\mathbb{N}} demonstrates that ρ\rho is not q-MLR as follows. Fix mm arbitrarily and let N(m)N(m) be as above. Recalling that MNM_{N} is the set consisting of all eigenvectors of SNS_{N} with non-zero eigenvalue, we have that,

ρ(Tm)tr(ργ(N)pγ(N))=tr(SNpγ(N))=tr(SN)=1.\rho(T_{m})\geq tr(\rho_{\gamma(N)}p_{\gamma(N)})=tr(S_{N}p_{\gamma(N)})=tr(S_{N})=1.

Since mm was arbitrary, infm(ρ(Tm))=1inf_{m\in\mathbb{N}}(\rho(T_{m}))=1. ∎

The following technical lemma, although seems unmotivated at this juncture, is crucial at a later point in the proof.

Lemma 3.6.

Let {[ai,bi]T}i=1n\{[a_{i},b_{i}]^{T}\}_{i=1}^{n} be a set of unit column vectors in 2\mathbb{C}^{2}. Let V=i=1n[ai,bi]TV=\bigotimes_{i=1}^{n}[a_{i},b_{i}]^{T} be their Kronecker product. If V=[v1,v2,,v2n]TV=[v_{1},v_{2},\dots,v_{2^{n}}]^{T}, then for all k2n1k\leq 2^{n-1}, we have that

|vk||v2nk+1|=i=1n|ai||bi|.|v_{k}||v_{2^{n}-k+1}|=\prod_{i=1}^{n}|a_{i}||b_{i}|.
Proof.

For natural numbers uu and qq, let [u]q[u]_{q} denote the remainder obtained by dividing uu by qq. We use the following convention for the Kronecker product [11]:

[a1b1][a2b2]=[a1a2b1a2a1b2b1b2].\begin{bmatrix}a_{1}\\ b_{1}\\ \end{bmatrix}\otimes\begin{bmatrix}a_{2}\\ b_{2}\\ \end{bmatrix}=\begin{bmatrix}a_{1}a_{2}\\ b_{1}a_{2}\\ a_{1}b_{2}\\ b_{1}b_{2}\end{bmatrix}.

So, v1=i=1naiv_{1}=\prod_{i=1}^{n}a_{i} and v2n=i=1nbiv_{2^{n}}=\prod_{i=1}^{n}b_{i}. For any k2n1k\leq 2^{n-1}, vkv_{k} has the form vk=i=1ncikv_{k}=\prod_{i=1}^{n}c^{k}_{i}, for some cik{ai,bi}c^{k}_{i}\in\{a_{i},b_{i}\} and v2nk+1v_{2^{n}-k+1} has the form v2nk+1=i=1neikv_{2^{n}-k+1}=\prod_{i=1}^{n}e^{k}_{i}, for some eik{ai,bi}e^{k}_{i}\in\{a_{i},b_{i}\}. Note that c1k=a1c^{k}_{1}=a_{1} if and only if kk is odd if and only if e1k=b1e^{k}_{1}=b_{1}. Similarly, we have the following. c2k=a2c^{k}_{2}=a_{2} if and only if [k]22{1,2}[k]_{2^{2}}\in\{1,2\} if and only if e2k=b2e^{k}_{2}=b_{2}. c3k=a3c^{k}_{3}=a_{3} if and only if [k]23{1,,22}[k]_{2^{3}}\in\{1,\dots,2^{2}\} if and only if e3k=b3e^{k}_{3}=b_{3}. In general, for ini\leq n, for all k2n1k\leq 2^{n-1},

cik=ai[k]2i{1,,2i1}eik=bi.c^{k}_{i}=a_{i}\iff[k]_{2^{i}}\in\{1,\dots,2^{i-1}\}\iff e^{k}_{i}=b_{i}.

This proves the lemma. Intuitively, this happens for the following reason. Imagine moving from v1v_{1} to v2n1v_{2^{n-1}} (by incrementing kk) and keeping track of the values of cikc^{k}_{i} as you move along the vkv_{k}s. Also, imagine moving from v2nv_{2^{n}} to v2n1v_{2^{n-1}} and keeping track of the values of eike^{k}_{i} as you move along the v2nk+1v_{2^{n}-k+1}s. Both motions are in opposite directions since as kk is incremented, the first motion is from lower to higher indices and the second is from higher to lower indices. Consider the behavior of c1k,e1kc^{k}_{1},e^{k}_{1} as kk is incremented. At the ‘start’ point, c11=a1c^{1}_{1}=a_{1}, e11=b1e^{1}_{1}=b_{1}. Now, as you move (i.e as you increment kk), c1kc^{k}_{1} alternates between a1a_{1} and b1b_{1} equalling it’s starting value, a1a_{1} at odd kks and e1ke^{k}_{1} alternates between b1b_{1} and a1a_{1} equalling it’s starting value b1b_{1} for odd kks. Now, take any ini\leq n. cikc^{k}_{i} alternates between aia_{i} and bib_{i} in blocks of length 2i12^{i-1}. cik=aic^{k}_{i}=a_{i} when kk is in the first block, {1,2,,2i1}\{1,2,\dots,2^{i-1}\} (i.e, when [k]2i{1,2,,2i1}[k]_{2^{i}}\in\{1,2,\dots,2^{i-1}\}) and cik=bic^{k}_{i}=b_{i} when kk is in the second block, {2i1+1,,2i}\{2^{i-1}+1,\dots,2^{i}\}(i.e, when [k]2i{2i1+1,,0}[k]_{2^{i}}\in\{2^{i-1}+1,\dots,0\}) and so on. Similarly, eike^{k}_{i} alternates between bib_{i} and aia_{i} in blocks of length 2i12^{i-1}. ∎

Lemma 3.7.

Let nn\in\mathbb{N} and let {[ai,bi]T}i=1n\{[a_{i},b_{i}]^{T}\}_{i=1}^{n} be such that for all i, [ai,bi]T[a_{i},b_{i}]^{T} is unit column vector in 2\mathbb{C}^{2} and let W=i=1n[ai,bi]TW=\bigotimes_{i=1}^{n}[a_{i},b_{i}]^{T}. Then, |<W|dn|W>|[2n(12n1),2n(1+2n1)]|\big{<}W|d_{n}|W\big{>}|\in[2^{-n}(1-2n^{-1}),2^{-n}(1+2n^{-1})]

Proof.

Fix nn and VV as in the statement and write dnd_{n} as a block matrix with each block of size 2n12^{n-1} by 2n12^{n-1}.

dn=[ABBTA].d_{n}=\begin{bmatrix}A&B\\ B^{T}&A\\ \end{bmatrix}.

Letting V=i=1n1[ai,bi]TV=\bigotimes_{i=1}^{n-1}[a_{i},b_{i}]^{T}, in block form, W=[anVT,bnVT]TW=[a_{n}V^{T},b_{n}V^{T}]^{T}. Let V=[v1,v2,,v2n1]TV=[v_{1},v_{2},\dots,v_{2^{n-1}}]^{T}. It is easily checked that

<W|dn|W>=2n+anbnVBV+anbnVBTV.\big{<}W|d_{n}|W\big{>}=2^{-n}+a_{n}^{*}b_{n}V^{\dagger}BV+a_{n}b_{n}^{*}V^{\dagger}B^{T}V.

By the form of B we get,

VBV=2n[v1,v2,,v2n1][v2n1,v2n11,,v2n1rn+1,0,,0]T.V^{\dagger}BV=2^{-n}[v^{*}_{1},v^{*}_{2},\dots,v^{*}_{2^{n-1}}][v_{2^{n-1}},v_{2^{n-1}-1},\dots,v_{2^{n-1}-r_{n}+1},0,\dots,0]^{T}.
=2nk=1rnvkv2n1k+1.=2^{-n}\sum_{k=1}^{r_{n}}v^{*}_{k}v_{2^{n-1}-k+1}.

By the previous lemma,

|VBV|2nk=1rn|vk||v2n1k+1|=2nrni=1n1|ai||bi|=2nrni=1n1|ai|1|ai|2.|V^{\dagger}BV|\leq 2^{-n}\sum_{k=1}^{r_{n}}|v_{k}||v_{2^{n-1}-k+1}|=2^{-n}r_{n}\prod_{i=1}^{n-1}|a_{i}||b_{i}|=2^{-n}r_{n}\prod_{i=1}^{n-1}|a_{i}|\sqrt{1-|a_{i}|^{2}}.

Since x1x2x\sqrt{1-x^{2}} has a maximum value of 1/21/2 and recalling definition of rnr_{n},

|VBV|2n12n12nn=21nn.|V^{\dagger}BV|\leq 2^{-n}\dfrac{1}{2^{n-1}}\dfrac{2^{n}}{n}=\dfrac{2^{1-n}}{n}.

Similarly, |VBTV|21nn.|V^{\dagger}B^{T}V|\leq\dfrac{2^{1-n}}{n}. Noting that |anbn|,|anbn|1/2|a_{n}^{*}b_{n}|,|a_{n}b_{n}^{*}|\leq 1/2,

|<W|dn|W>|2n+|anbnVBV|+|anbnVBTV|2n+21nn,|\big{<}W|d_{n}|W\big{>}|\leq 2^{-n}+|a_{n}^{*}b_{n}V^{\dagger}BV|+|a_{n}b_{n}^{*}V^{\dagger}B^{T}V|\leq 2^{-n}+\dfrac{2^{1-n}}{n},

and

|<W|dn|W>|2n|anbnVBV||anbnVBTV|2n21nn.|\big{<}W|d_{n}|W\big{>}|\geq 2^{-n}-|a_{n}^{*}b_{n}V^{\dagger}BV|-|a_{n}b_{n}^{*}V^{\dagger}B^{T}V|\geq 2^{-n}-\dfrac{2^{1-n}}{n}.

Lemma 3.8.

ρ\rho is mR.

If pp is any measure on 2ω2^{\omega}, we can define Martin-Löf randomness with respect to pp exactly as we defined it for the uniform measure. Denote by MLR(p)MLR(p), the set of bitstrings Martin-Löf random with respect to pp [5].

Proof.

We use ideas similar to Theorem 196(a) in [5]. For convenience, for all i>5i>5, define

βi:=q=5i1q.\beta_{i}:=\sum_{q=5}^{i-1}q.

Let BB be any computable measurement system. We show that MLR(μρB)MLRMLR(\mu^{B}_{\rho})\subseteq MLR. Since μρB[MLR(μρB)]=1\mu^{B}_{\rho}[MLR(\mu^{B}_{\rho})]=1, this implies that μρB(MLR)=1\mu^{B}_{\rho}(MLR)=1. Denote μρB\mu^{B}_{\rho} by μ\mu for convenience. Let λ\lambda denote the uniform measure. We will abuse notation by writing μ(τ)\mu(\tau) instead of the more cumbersome μ(τ)\mu(\llbracket\tau\rrbracket) for τ2<ω\tau\in 2^{<\omega}. Let XMLR(μ)X\in MLR(\mu). Write XX as a concatenation of finite bitstrings : X=σ5σ6σnX=\sigma_{5}\sigma_{6}\dots\sigma_{n}\dots where σn2n\sigma_{n}\in 2^{n} for all nn\in\mathbb{N}. Let Sn:=σ5σ6σnS_{n}:=\sigma_{5}\sigma_{6}\dots\sigma_{n} be the concatenation upto nn. Let μi\mu_{i} be such that for all τ2i\tau\in 2^{i},

μi(τ):=tr[di(|q=1ibτ(q)q+βi><q=1ibτ(q)q+βi|)].\mu_{i}(\tau):=tr\big{[}d_{i}(|\bigotimes_{q=1}^{i}b^{q+\beta_{i}}_{\tau(q)}\big{>}\big{<}\bigotimes_{q=1}^{i}b^{q+\beta_{i}}_{\tau(q)}|)\big{]}.

By 2.1 and by the form of ρ\rho we see that,

μ(Sn)=i=5nμi(σi).\mu(S_{n})=\prod_{i=5}^{n}\mu_{i}(\sigma_{i}).

Note that μ\mu is computable [5] since ρ\rho and BB are. Since XMLR(μ)X\in MLR(\mu), by the Levin-Schnorr theorem (Theorem 90, section 5.6 in [5]) there is a C1C_{1} such that

n,log(μ(Sn))C1KM(Sn).\displaystyle\forall n,-\log(\mu(S_{n}))-C_{1}\leq KM(S_{n}).

By Theorem 89, section 5.6 in [5] fix a C2C_{2} such that

n,KM(Sn)log(λ(Sn))+C2.\displaystyle\forall n,KM(S_{n})\leq-\log(\lambda(S_{n}))+C_{2}.

By these inequalities and taking exponents, we see that there is a constant α>0\alpha>0 such that

n,μ(Sn)αλ(Sn).\forall n,\mu(S_{n})\geq\alpha\lambda(S_{n}).

Letting ri:=μi(σi)r_{i}:=\mu_{i}(\sigma_{i}) and δi:=λ(σi)ri\delta_{i}:=\lambda(\sigma_{i})-r_{i} in the above,

n,i=5nriαi=5nri+δi.\displaystyle\forall n,\prod_{i=5}^{n}r_{i}\geq\alpha\prod_{i=5}^{n}r_{i}+\delta_{i}. (3.5)

Let μ\mu^{\prime} be a probability measure on 2ω2^{\omega} such that for all σ2<ω,μ(σ):=2μ(σ)λ(σ)\sigma\in 2^{<\omega},\mu^{\prime}(\sigma):=2\mu(\sigma)-\lambda(\sigma). In particular, this implies that

n,μ(Sn)=i=5nriδi.\forall n,\mu^{\prime}(S_{n})=\prod_{i=5}^{n}r_{i}-\delta_{i}.

Note that μ\mu^{\prime} is computable since μ\mu and λ\lambda are. Applying the same argument which resulted in 3.5, we get that there is an ε>0\varepsilon>0 such that,

n,i=5nriεi=5nriδi.\displaystyle\forall n,\prod_{i=5}^{n}r_{i}\geq\varepsilon\prod_{i=5}^{n}r_{i}-\delta_{i}. (3.6)

By Lemma 3.7, for all i,ri[2i(12i1),2i(1+2i1)].i,r_{i}\in[2^{-i}(1-2i^{-1}),2^{-i}(1+2i^{-1})]. So, |δi|=|ri2i|[0,2i+1i1]|\delta_{i}|=|r_{i}-2^{-i}|\in[0,2^{-i+1}i^{-1}]. Hence,
ri+δi2i2i+1i12i+1i1=2i[14i1]>0r_{i}+\delta_{i}\geq 2^{-i}-2^{-i+1}i^{-1}-2^{-i+1}i^{-1}=2^{-i}[1-4i^{-1}]>0, since i5i\geq 5. Similarly, riδi0r_{i}-\delta_{i}\geq 0. By this, multiplying 3.5 and 3.6 gives,

n,i=5nri2αεi=5nri2δi2=αεi=5nri2i=5n(1δi2ri2).\displaystyle\forall n,\prod_{i=5}^{n}r^{2}_{i}\geq\alpha\varepsilon\prod_{i=5}^{n}r^{2}_{i}-\delta^{2}_{i}=\alpha\varepsilon\prod_{i=5}^{n}r^{2}_{i}\prod_{i=5}^{n}\big{(}1-\dfrac{\delta^{2}_{i}}{r^{2}_{i}}\big{)}. (3.7)

By the above,

|δi|ri2i+1i12i(12i1)=2(i2)1.\dfrac{|\delta_{i}|}{r_{i}}\leq\dfrac{2^{-i+1}i^{-1}}{2^{-i}(1-2i^{-1})}=2(i-2)^{-1}.

Letting F>0F>0 be the constant,

n,i=5n(1δi2ri2)i=5(1δi2ri2)i=5(14(i2)2)=F,\displaystyle\forall n,\prod_{i=5}^{n}\big{(}1-\dfrac{\delta^{2}_{i}}{r^{2}_{i}}\big{)}\geq\prod_{i=5}^{\infty}\big{(}1-\dfrac{\delta^{2}_{i}}{r^{2}_{i}}\big{)}\geq\prod_{i=5}^{\infty}\big{(}1-4(i-2)^{-2}\big{)}=F,

3.7 gives,

n,(αε)1i=5nri2i=5nri2δi2i=5nri2F.\displaystyle\forall n,(\alpha\varepsilon)^{-1}\prod_{i=5}^{n}r^{2}_{i}\geq\prod_{i=5}^{n}r^{2}_{i}-\delta^{2}_{i}\geq\prod_{i=5}^{n}r^{2}_{i}F. (3.8)

From 3.5, 3.6 and 3.8, it is easy to see that there is a G>0G>0 such that for all nn

i=5nri+δiGi=5nri.\prod_{i=5}^{n}r_{i}+\delta_{i}\geq G\prod_{i=5}^{n}r_{i}.

Recalling the definitions of rir_{i} and δi\delta_{i},

n,λ(Sn)Gμ(Sn).\forall n,\lambda(S_{n})\geq G\mu(S_{n}).

Letting D=C1log(G)D=C_{1}-\log(G) and recalling the definition of C1C_{1},

n,log(λ(Sn))log(μ(Sn))log(G)KM(Sn)+D.\forall n,-\log(\lambda(S_{n}))\leq-\log(\mu(S_{n}))-\log(G)\leq KM(S_{n})+D.

By Theorem 85 in [5], KM(.)K(.)+O(1)KM(.)\leq K(.)+O(1) and so there is a E>0E>0 such that

n,log(λ(Sn))K(Sn)+E.\forall n,-\log(\lambda(S_{n}))\leq K(S_{n})+E.

Noting that log(λ(Sn))=|Sn|=βn+n-\log(\lambda(S_{n}))=|S_{n}|=\beta_{n}+n, 3.2.14 from [9] implies that XX is MLR. ∎

The theorem is proved. ∎

4 Generalizations and future directions

We sketch some ways in which the previous section’s results generalize. Given S2ωS\in 2^{\omega}, we may relativize the notion of Martin-Löf randomness to define the set MLRS2ωMLR^{S}\subset 2^{\omega} of infinite bitstrings which are Martin-Löf random with respect to SS. The halting problem, denoted by \emptyset^{\prime}\subset\mathbb{N} is an incomputable set important in computability theory. Letting (n)\emptyset^{(n)} be the n1n-1th iterate of the halting problem, an element of Cantor space is said to be arithmetically random if it is in MLR(n)MLR^{\emptyset^{(n)}} for every nn (see 6.8.4 in [7]). Given S2ωS\in 2^{\omega}, relativizing the proof of Lemma 3.8 shows that MLRS(μρB)MLRSMLR^{S}(\mu^{B}_{\rho})\subseteq MLR^{S} as follows. Take an XMLRS(μρB)X\in MLR^{S}(\mu^{B}_{\rho}). Relativizing Theorems 85 and 90 from [5] and 3.2.14 from [9] to SS and noting that KMS(.)KM(.)KM^{S}(.)\leq KM(.) and following the proof of Lemma 3.8 shows that XMLRSX\in MLR^{S}. This shows that μρB(MLRS)=1\mu^{B}_{\rho}(MLR^{S})=1 holds for any S2ωS\in 2^{\omega} and any computable measurement system BB. In particular, this has an interesting application; if BB is any computable measurement system, μρB(MLR(n))=1\mu^{B}_{\rho}(MLR^{\emptyset^{(n)}})=1 for all nn. So,

μρB[n(MLR(n))]=1.\mu^{B}_{\rho}\big{[}\bigcap_{n\in\mathbb{N}}(MLR^{\emptyset^{(n)}})\big{]}=1.

So, measuring ρ\rho in any computable measurement system yields an arithmetically random infinite sequence of bits, with probability one. The above note naturally suggests a definition:

Definition 4.1.

ρ\rho is said to be strong measurement random (strong mR), if μρB(MLRS)=1\mu^{B}_{\rho}(MLR^{S})=1 holds for any S2ωS\in 2^{\omega} and any computable measurement system BB.

By Remark 2.3 and by the above discussion on relativizations, we can also consider measurement of a state in non-computable measurement systems by using an appropriate oracle. We do not explore this here.

One may ask if we can build other computable examples of ρ\rhos which are not q-MLR and are mR. We note that a straightforward modification of the proof of Theorem 3.4 yields a family of such ρ\rhos. We do not provide all the details here for lack of space. Let h:h:\mathbb{N}\longrightarrow\mathbb{N} and g:(0,1)g:\mathbb{N}\longrightarrow(0,1) be computable, satisfying the following for some constants δ(0,1)\delta\in(0,1) and F>0F>0:

n=5(1h(n)2n)=0,n=5(1h(n)[2ng(n)])=δ,\prod_{n=5}^{\infty}(1-h(n)2^{-n})=0,\prod_{n=5}^{\infty}(1-h(n)[2^{-n}-g(n)])=\delta,
n,g(n)2n and n=5[14g2(n)h2(n)(12g(n)h(n))2]=F.\forall n,g(n)\leq 2^{-n}\text{ and }\prod_{n=5}^{\infty}\big{[}1-\dfrac{4g^{2}(n)h^{2}(n)}{(1-2g(n)h(n))^{2}}\big{]}=F.

Let ρ\rho be defined as in the proof of the main Theorem but with rnr_{n} replaced by h(n)h(n) and with the h(n)h(n) many entries on the extreme ends of the anti-diagonal of dnd_{n} being equal to g(n)g(n) instead of 2n2^{-n}. Then, this ρ\rho is computable and mR (in fact, it is strong mR) and fails a q-MLT at order δ\delta.

We are working towards characterizing the set of states for which mR and q-MLR are equivalent. So far, we have shown that these notions are equivalent for states of the form n=1d\otimes_{n=1}^{\infty}d for some computable density matrix dd.

One may imagine using a sequence of POVMs in 2.1 instead of a sequence orthonormal bases. We use orthonormal bases to conform with the work [10] on which ours is based, which uses projective measurements and not POVMs. It would also be interesting to replicate the approach of [10] using POVMs instead of projective measurements.

5 Concluding comments

Measuring a finite dimensional quantum system or a composite system of finitely many qubits is a pivotal concept in quantum information theory [6]. It hence seems natural to consider defining a notion of measuring a state. Since measurement of a state yields an infinite sequence of bits, it is interesting to explore the relation between the randomness of the measured state and the randomness of the resulting sequence. This paper is motivated by these questions. The main result is that q-MLR is not equivalent to mR even for the computable states.

5.1 Remark

Intuitively, the non-equivalence of mR and q-MLR should not be surprising given that entanglement in composite systems cannot be detected by independent measurements of the subsystems. Let us elaborate on this remark. ρ\rho in 3.4 is built up from dnd_{n}s where each dnd_{n} has rnr_{n} many entangled eigenvectors with non-zero eigenvalue and rnr_{n} many entangled eigenvectors with zero eigenvalue. This inhomogeneity in the distribution of eigenvalues is solely due to these entangled eigenvectors (all the 2n2rn2^{n}-2r_{n} many non entangled eigenvectors of dnd_{n} have the same non-zero eigenvalues). A crucial part in showing that ρ\rho is non q-MLR was to use the inhomogeneous eigenvalue distribution to bound the size MNM_{N} (see 3 in the proof of 3.5). Heuristically speaking, the the non-quantum randomness of ρ\rho was a reflection of the non-uniform eigenvalue distribution of dnd_{n} which in turn was due to by the presence of entangled eigenvectors of dnd_{n}. It is hence reasonable to expect that the quantum non-randomness of ρ\rho, which stems from entanglement, cannot be captured by measurements in the sense of 2.1 using pure tensors (i.e. measuring each 2-dimensional subsystem independently).

5.2 Notes

We have shown how to extract classical arithmetic randomness from a computable, non-quantum random sequence of qubits. It seems plausible that our results may prove to be relevant to the construction of quantum random number generators [8]. Abbott, Calude and Svozil have also studied classical bit sequences resulting from measuring a quantum system[2, 4]. However, their notion of measurement is significantly different from ours. In contrast to our work which considers measurement of an infinite sequence of qubits, they studied the randomness of a sequence of bits generated by repeatedly measuring a finite dimensional quantum system. They go on to apply this to quantum random number generators and their certification [2, 3, 4].

5.3 Acknowledgements

I thank James Hanson for conjecturing that q-MLR is equivalent to mR when dd (see section 4) is four by four in the above. Joe Miller and Peter Cholak (independently) asked if there is a notion of ‘measuring a state’ and if one gets a MLR bitstring from measuring a q-MLR state with respect to such a notion. These questions were one of the factors which led me to explore this area. André Nies, whom I thank for introducing me to quantum algorithmic randomness, independently suggested that one might get a measure on Cantor space by ‘measuring’ a state. I thank Joe Miller, my thesis advisor, for his encouragement and guidance.

References

  • [1]
  • [2] Alastair A. Abbott (2015): Value Indefiniteness, Randomness and Unpredictability in Quantum Foundations. (De la Valeur Indéfinie aux Notions d’Aléatoire et d’Imprévisibilité Quantiques). Ph.D. thesis, École Normale Supérieure, Paris, France.
  • [3] Alastair A. Abbott, Cristian S. Calude & Karl Svozil (2014): A quantum random number generator certified by value indefiniteness. Mathematical Structures in Computer Science 24(3), 10.1103/PhysRevA.81.012109.
  • [4] Alastair A. Abbott, Cristian S. Calude & Karl Svozil (2015): On the Unpredictability of Individual Quantum Measurement Outcomes. In: Fields of Logic and Computation II, Lecture Notes in Computer Science 9300, Springer, pp. 69–86, 10.1038/438743a.
  • [5] A Shen et al. (2017): Kolmogorov Complexity and Algorithmic Randomness. Mathematical Surveys and Monographs : 220, AMS. ISBN: 978-1-4704-3182-2.
  • [6] Isaac Chuang & Michael Nielsen (2000): Quantum Computation and Quantum Information, second edition. Cambridge Univ. Press. ISBN-13: 978-1107002173.
  • [7] Rodney G. Downey & Denis R. Hirschfeldt (2010): Algorithmic Randomness and Complexity, 10.1007/978-0-387-68441-3.
  • [8] Miguel Herrero-Collantes & Juan Carlos Garcia-Escartin (2017): Quantum random number generators. Rev. Mod. Phys. 89, p. 015004, 10.1103/RevModPhys.89.015004.
  • [9] André Nies (2009): Computability and Randomness. ISBN: 0199230765, 9780199230761.
  • [10] André Nies & Volkher B. Scholz (2019-09): Martin-Löf random quantum states. Journal of Mathematical Physics 60(9), p. 092201, 10.1063/1.5094660.
  • [11] P. A. Regalia & S. K. Mitra (1989): Kronecker Products, Unitary Matrices, and Signal Processing Applications. SIAM Rev. 31(4), pp. 586–613, 10.1137/1031127.