arXiv:1703.10118v2 [quant-ph] 18 Jun 2018 otx fapiueapicto,w yial ee otecomputa the to refer typically algor we search n Grover amplification, the amplitude of of version context generalized a describe further eso fGoe’ loih r itdi pedxA. resu Appendix Necessary in obtained. listed results are am the algorithm Grover’s of of of scenarios summary version a model with some o paper terms within in the exemplified coherence are of relations entropy relative trade-off conte the probability on quantifier success algorithmic coherence estimate the the two-sided on and material coherence in quantum required utilized between the relations review be betwee we will relations II, coherence general Sect. formulate geometric will the F we Third, addressed. Second, mainly are considered. aspects be following The amplification. amplitude n nageetmngm ntedsrt nlgeo nlgGrov analog [25]. of in inc analogue addressed discrete for first were Grover the resource original computations quantu in a the monogamy performing in be entanglement depletion in coherence may and coherence examined coherence [23] quantum is of whether of level authors ask quantum role the they the at particular, of coherence Studies th quantifying of of 21]. quan states problem about which The think in should bases, basis. we of impossib speed-up, number computational be reg quantum limited to the quantum only a seems with of entanglement connected states without are entangled speed-up use quantum to that algo assumes Grover [15] the Deutsch of of versions generalized of Details operatio compu quantum rotation amplitudes. addition, th the the In in to of one. phases unitary th Due changing arbitrary of with by blocks transform [11]. modified building algorithm be formulate may Grover to algorithm the attempted have database by researchers whereas searchin inspired [11], for item, algorithms optimal any quantum clo evaluate is of more to algorithm are Grover’s oracle algorithms algorithms that the Shor’s quantum recognized and numerous soon Grover’s to was that led arguments have gave [4] [8] results Shor’s ebrated nti eto,w ealdfiiin fqatmchrnequantifi coherence quantum of definitions recall we section, this In uis hntebssi omdby formed is basis the Then qubits. ti elkonta nageeti e eorei unu infor quantum in resource key a is entanglement that known well is It discoverie fundamental of one is [1–3] algorithm search Grover The h i fti oki osuytaeo eain ewe quantum between relations trade-off study to is work this of aim The nterl fdaigwt unu oeec napiueam amplitude in coherence quantum with dealing of role the On eateto hoeia hsc,IktkSaeUniversi State Irkutsk Physics, Theoretical of Department ewrs rvrsagrtm unu erh eaieen relative search, quantum algorithm, Grover’s Keywords: illus are presented processes. relations amplification general amplitude wit The relations coheren expository. complementarity of most coherence, entropy relative of amplification entropy the amplitude ative and in coherence b probability geometric relations success the trade-off the addition examine and We In ence prepared. one. be ope unitary may rotation al arbitrary amplitudes the search with in Grover phases transform changing the Hadamard with of connected ideas are key tions generalizes technique This mltd mlfiaini n fpiaytosi uliga building in tools primary of one is amplification Amplitude N 2 = I PRELIMINARIES II. .INTRODUCTION I. n lxyE Rastegin E. Alexey rhnra es ahket Each kets. orthonormal e se per sntrpeetdepiil.Tdy hr xssaclass a exists there Today, explicitly. represented not is rp,goerccoherence geometric tropy, tm tdigterl fqatmchrnei the in coherence quantum of role the Studying ithm. n ehiuso mltd mlfiain General amplification. amplitude of techniques and s y aai v 0 rus 603 Russia 664003, Irkutsk 20, Bv. Gagarin ty, r osdrdi et I.I atclr edrv a derive we particular, In III. Sect. in considered are h ucs rbblt.I et V h presented the IV, Sect. In probability. success the f se.Terslso h aes[6 7 manifested 17] [16, papers the of results The ister. u orltoswt epc otecomputational the to respect with correlations tum igmysatwt nabtayiiildistribution initial arbitrary an with start may ting h ucs rbblt emt ethe be to seem probability success the h unu oeec n h ucs probability. success the and coherence quantum n t fsligrcrineutosfrageneralized a for equations recursion solving of lts rtdwti eea oe cnro of scenarios model several within trated loih sgnrlya osbe h original The possible. as generally as algorithm e loih.Teppr[4 sdvtdt coherence to devoted is [24] paper The algorithm. esn h oe fqatmagrtm.The algorithms. quantum of power the reasing ih eedsrbdi [12–14]. in described were rithm eaecniee.I em fterel- the of terms In considered. are ce ain n elcn h intermediate the replacing and rations e nteohrhn,qatmcomputations quantum hand, other the On le. oih.Ptnilyueu modifica- useful Potentially gorithm. yqeist rce[,1] ee einvoke we Here, 10]. [9, oracle to queries by g rcse.A esrso coherence, of measures As processes. eyrltdta n ih xeta rt It first. at expect might one than related sely ltd mlfiain nSc.V econclude we V, Sect. In amplification. plitude inlbss e h unu eitrcontain register quantum the Let basis. tional rirr nta itiuino the of distribution initial arbitrary , opttoswr eotdi 2–4.In [22–24]. in reported were computations m urnl h ujc fatv eerhs[18– researches active of subject the currently r htwl eue hog h ae.We paper. the through used be will that ers rsac.Sc nlge fdgtlquantum digital of analogues Such search. er eitraerpeetd nodrt analyze to order In represented. are register e htmtvt unu opttos Cel- computations. quantum motivate that s o leri rbes[–] h uhr of authors The [5–7]. problems algebraic for tenmaue fqatmcoher- quantum of measures etween rt oegnrlseaiso optn will computing of scenarios general more irst, sadrpaigteitreit Hadamard intermediate the replacing and ns ainpoesn.Teqatmparallelism quantum The processing. mation grtm o unu computers. quantum for lgorithms t h ae sognzda olw.In follows. as organized is paper The xt. oeec n h ucs rbblt in probability success the and coherence ra plcblt fsac problems search of applicability broad e | x i sidxdb binary by indexed is plification n srn with -string 2 xj 0, 1 . A rigorous framework for the quantification of coherence was proposed in [18]. We consider the set of all∈{ diagonal} density matrices written as I
N−1 δ = δ x x . (1) x | ih | x=0 X One further asks how far the given state is from states of the form (1). The authors of [18] listed general conditions for quantifiers of coherence. From a general perspective, coherence measures are discussed in [19–21]. In the present paper, we will mainly use the relative entropy of coherence and the geometric coherence. Using the quantum relative entropy as a measure of distinguishability, we arrive at the relative entropy of coherence. The quantum relative entropy of ρ with respect to σ is defined as [26, 27]
tr(ρ ln ρ ρ ln σ) , if ran(ρ) ran(σ) , D (ρ σ) := − ⊆ (2) 1 || + , otherwise . ( ∞ By ran(ρ), we mean here the range of ρ. The role of the quantum relative entropy in quantum information science is considered in [26–28]. The relative-entropy-based measure of coherence is defined as [18]
C1(ρ) := min D1(ρ δ) . (3) δ∈I || The minimization finally gives [18]
C (ρ)= S (ρ ) S (ρ) , (4) 1 1 diag − 1 where S (ρ) = tr(ρ ln ρ) is the von Neumann entropy, and ρ is obtained from ρ by vanishing all off-diagonal 1 − diag elements calculated in the basis x . The first term in the right-hand side of (4) is equal to the Shannon entropy of the discrete distribution with probabilities| i p(x)= x ρ x , viz. h | | i N−1 S (ρ )= H (p) := p(x) ln p(x) . (5) 1 diag 1 − x=0 X The maximal value of (5) is equal to ln N, so that the quantity (4) cannot exceed ln N. Properties of the relative entropy of coherence are discussed in [18, 20, 21]. Uncertainty relations in terms of the relative entropy of coherence were derived in [29–31]. Coherence monotones of the Tsallis type were examined in [32]. Relative R´enyi entropies of coherence were considered in [33, 34]. In contrast to (4), these quantifiers do not reduce to a simple analytical expression. Several distance-based quantifiers of coherence were considered [18, 20]. The ℓ1-norm of coherence is often used as intuitively natural quantifier. For the given state ρ, we define
Cℓ1 (ρ) := min ρ δ = x ρ y . (6) δ∈I − ℓ1 h | | i x=6 y X
Complementarity and uncertainty relations in terms of the quantity (6) were studied in [35, 36]. A factorization relation for the ℓ1-norm of coherence was proved in [37]. Another intuitively attractable way is to use the squared ℓ2-norm. However, the corresponding coherence quantifier does not satisfy the monotonicity requirement [18]. The trace norm also induces an interesting candidate to quantify the amount of coherence [38, 39]. We will use the geometric coherence which is introduced in terms of the quantum fidelity [40, 41]. Following Jozsa [41], the fidelity of density matrices ρ and σ is expressed as
ρ σ ρ √σ 2 F ( , )= √ 1 , (7)
† where A 1 = tr √A A is the trace norm. Another known definition is written as the square root of (7) [26, 27]. The fidelityk k ranges between 0 and 1 taking the value 1 for two identical states. Using the unity minus fidelity as a distance measure, one defines the geometric coherence by [20]
Cg(ρ) := 1 max F (ρ, δ) . (8) − δ∈I 3
Properties of this coherence quantifier are summarized in subsection III.C.3 of [20]. For pure states, the formula (8) reduces to 2 Cg ψ =1 max x ψ . (9) | i − x h | i We will also use some bounds on the geometric coherence. To each state ρ, we assign the index of coincidence
N−1 I(ρ) := p(x)2 , (10) x=0 X where p(x)= x ρ x . The authors of [42] have proved that h | | i N 1 N − 1 1 tr(ρ2) I(ρ) C (ρ) . (11) N − − N 1 − ≤ g ( r − ) Using the unity minus the square root of (7), we also have a distance measure. The corresponding coherence quantifier was studied in [18, 38]. Duality relations between the coherence and path information were examined in [43–45]. The authors of [46] proposed the concept of relative quantum coherence. Let us recall some results concerning generalizations of the Grover algorithm. In our calculations, we will follow the scheme of analysis developed in [12, 13]. Suppose that the search space contains N =2n items indexed by binary n-string x = (x1 xn) with xj 0, 1 so that x 0, 1,...,N 1 . The problem is to find one of marked items which form the set··· . In amplitude∈ { amplification,} ∈ we{ try to increase− } maximally amplitudes of states x just for x . Without lossM of generality, the number of marked items is assumed to obey 1 N/2. | i ∈MEach concrete item is tested by means of oracle-computable Boolean function x ≤f |M|(x) such≤ that f(x) = 1 for x and f(x) = 0 for x C. The original Grover algorithm starts with initializing7→ n-qubit register to 0 and applying∈ M the Hadamard transform∈ M to get a uniform amplitude distribution | i
− 1 N 1 H 0 = x . (12) | i √ | i N x=0 X Following [12, 13], we will use phase rotations of two kinds. Rotations of the first kind non-trivially act on unknown marked states. For an arbitrary phase γ, we write
N−1 J (γ)= exp iγf(x) x x . (13) f | ih | x=0 X Thus, the operation (13) rotates all the marked states by the phase γ. For γ = π, we have x x for f(x) = 1, as in the original formulation. Its realization with invoking the oracle black box is well explained| i 7→ − in | i the literature. Rotations of the second class rotate some prescribed state s of the basis. For any phase β, one reads | i J (β)= 11 1 exp(iβ) s s , (14) s − − | ih | where 11 is the identity operator. So, the ket s is rotated by β, whereas other basis kets remain unchanged. One of the well known ingredients of Grover’s algorithm| i can be written as
J (π)=2 0 0 11 . (15) − 0 | ih |− Clothing (15) by the Hadamard transforms, we obtain the inversion about the average. So, the original Grover iteration is written as G = H J (π) H J (π) . (16) 0 − 0 f For many reasons, generalizations of Grover’s original algorithm have been developed. In the generalized formula- tion, the iteration (16) is replaced by G = U J (β) U† J (γ) , (17) − s f with arbitrary β and γ. Here, the Hadamard transform on n qubits is replaced by an arbitrary unitary operator. To the given s and U, we assign η = U s , whence | i | i | i U J (β) U† = 1 exp(iβ) η η 11 . (18) − s − | ih |− 4
In this way, the terms U and s are both taken into account by η solely. After t iterations, the state of the quantum register is represented as | i | i
g(t) = k (t) x + ℓ (t) y . (19) | i x | i y | i ∈M C xX yX∈M Iterations of the form (17) will be applied to an arbitrary initial distribution of amplitudes g(0) . It must be stressed that initializing the state of the quantum register can be quite challenging [47]. Grover’s search| algorithmi for a mixed initial state of the register was analyzed in [14]. Any single iteration changes amplitudes according to the equations
k (t +1)= x G g(t) , ℓ (t +1)= y G g(t) , x h | | i y h | | i C where x and y . These equation are accompanied by initial amplitudes kx(0) and ℓy(0) of marked and unmarked∈ items, M respectively.∈ M The authors of [13] have examined the above recursion equations. Their results required for us are summarized in Appendix A.
III. GENERAL RELATIONS BETWEEN QUANTUM COHERENCE AND THE SUCCESS PROBABILITY
In this section, we focus on general relations concerning coherence changes in amplitude amplification. In practice, computing devices are inevitably exposed to noise. In this situation, the state of a quantum register after t steps is described by density matrix ρ(t). Then the probability to measure one of the marked state is written as
P (t)= x ρ(t) x . (20) suc h | | i ∈M xX Amplitude amplification processes aim to magnify amplitudes of desired states unknown a priori. In order to increase the success probability Psuc(t), a coherence of the register should be used somehow. For original Grover’s formulation, this issue was addressed in [23]. The writers of [23] also concerned other quantum correlations such as pairwise or multipartite entanglement and discord. We shall analyze a coherence behavior in more general setting. We begin with a simple situation, where the geometric coherence immediately links to the success probability. For pure states, the geometric coherence is expressed by (9). We first suppose that only one amplitude should be maximized in a particular case. This situation is sufficiently typical, but not general. For a time, we also assume that states of a quantum register are pure during the performance of algorithm. After a proper number of iterations, amplitudes in some superposition ψ will be sufficiently small, except for the unique one. We then have | i C ψ + P ψ =1 . (21) g | i suc | i Under the described circumstances, an increase of the success p robability implies a decrease of the geometric coherence, and vice versa. In the form of inequality, we can exceed the above relation to mixed states of a quantum register and arbitrary number of marked states.
Proposition 1 Let ρ be a density matrix normalized as tr(ρ)=1, and let Psuc(ρ) be defined according to (20). The geometric coherence of ρ satisfies
N 1 N P − 1 1 tr(ρ2) P 2 (1 P )2 C (ρ) 1 suc . (22) N − − N 1 − suc − − suc ≤ g ≤ − M ( r − ) Proof. To prove the right-hand side of (22), we consider the diagonal state 1 δ = x x . 0 M | ih | ∈M xX Using the property P4(b) of [41], we can write the inequalities
Psuc max F (δ, ρ) F (δ0, ρ) tr(δ0ρ)= . (23) δ∈I ≥ ≥ M Combining (8) with (23) completes the proof of the upper bound. 5
C Let us proceed to the left-hand side of (22). For all x and y , one has p(x) Psuc and p(y) 1 Psuc, so that ∈M ∈M ≤ ≤ −
p(x)2 P 2 , p(y)2 (1 P )2 . (24) ≤ suc ≤ − suc ∈M C xX yX∈M 2 2 Hence, we obtain I(ρ) Psuc + (1 Psuc) . Combining the latter with (11) immediately gives the left-hand side of (22). ≤ − The statement of Proposition 1 is a two-sided estimate of the geometric coherence in terms of the success probability Psuc(ρ). During the proof, we have seen another relation that deserves to be given explicitly. It follows from (23) that
Psuc(ρ) M max F (δ, ρ) . (25) ≤ δ∈I Operationally, this relation shows that the success probability is bounded by the number of marked objects multiplied by the maximum overlap of the state with an incoherent state. It should be noted that the left-hand side of (22) provides a non-trivial lower bound only when
P 2 + (1 P )2 tr(ρ2) . (26) suc − suc ≤ This condition always holds for pure states, since tr(ρ2) = 1 for ρ = ψ ψ . Returning to the case of amplitude amplification, we see the following. If the success probability is determined| ih | by the unique state of the calculation basis, i.e., M = 1, then
C (ρ)+ P (ρ) 1 . (27) g suc ≤ In reality, states of the quantum register are inevitably exposed to noise. With some amount even small, they will become mixed. So, the trade-off relation between quantum coherence and the success probability may only be restricted here. Note that the geometric coherence is not a coherence measure in the sense that is proposed by the authors of [20]. They added the list of axioms for coherence quantifiers by two items called the uniqueness for pure states and the additivity. In this regard, we are also interested in relations of the success probability with other coherence quantifiers. The coherence quantifier based on the relative entropy is one of the most natural measure for these purposes. At the same time, this measure is not connected with Psuc so immediately as the geometric coherence. Here, the following statement holds.
Proposition 2 Let ρ be a density matrix normalized as tr(ρ)=1. For the given Psuc(ρ) and S1(ρ), the relative entropy of coherence satisfies the inequalities
h (P ) S (ρ) C (ρ) 1 suc − 1 ≤ 1 M N M P ln + (1 P ) ln − S (ρ) , (28) ≤ suc P − suc 1 P − 1 suc − suc where h (P ) is the binary Shannon entropy. Using the argument of max p(x) = max x ρ x , we further define 1 suc h | | i P , arg max p(x) , Ω := suc ∈M (29) 1 P , otherwise . ( − suc Then the relative entropy of coherence satisfies
ln Ω S (ρ) C (ρ) . (30) − − 1 ≤ 1 Proof. To prove the right-hand side of (28), we will use the results of [48]. The authors of [48] addressed the situation, when independent events are somehow collected into nonempty and pairwise disjoint subsets. Upper bounds on the Shannon entropy of the original probability distribution is then expressed in terms of new probabilities C corresponding to these subsets. In our case, we have the subsets and with probabilities Psuc and 1 Psuc, respectively. Applying theorem 3 of [48] to the probabilities p(x), weM writeM −
C H (p) P ln |M| + (1 P ) ln |M | . 1 ≤ suc P − suc 1 P suc − suc 6
Combining this with the definition of C1(ρ) completes the proof of the right-hand side of (28). Let us proceed to the left-hand side of (28). We first observe that 1 1 p(x)=1 , p(y)=1 . P 1 P suc ∈M suc C xX − yX∈M Applying Jensen’s inequality to convex function ξ ln ξ then gives 7→ − 1 1 H (p) P ln p(x)2 (1 P ) ln p(y)2 . (31) 1 ≥ − suc P − − suc 1 P suc ∈M ! suc C ! xX − yX∈M Combining (24) with (31) and decreasing of the function ξ ln ξ finally gives H1(p) h1(Psuc). The proof of (30) begins with applying Jensen’s inequality7→ to − convex function ξ ≥ln ξ, so that 7→ − N−1 p(x) ln p(x) ln p(x)2 . (32) − ≥− ∈M x=0 ! xX X If arg max p(x) , then we have p(x) max p(x) P for all x 0, 1,...,N 1 , whence ∈M ≤ ≤ suc ∈{ − } N−1 p(x)2 P . (33) ≤ suc x=0 X If arg max p(x) , the right-hand side of (33) is replaced with 1 Psuc. Together with (32) and decreasing of ξ ln ξ, these6∈ facts M provide (30). − 7→ − If we know S1(ρ), then Proposition 2 gives a two-sided estimate of the relative entropy of coherence in terms of the success probability. Of course, productive computations should keep states of the register closely to the pure ones with small values of the von Neumann entropy. When Psuc approaches 1, a band of allowed values of C1(ρ) becomes more and more narrow. This picture characterizes good algorithms, which ensure high chances for the success. In general, we see some complementarity between coherence and the success probability. For the fixed Psuc(ρ) and S1(ρ), a difference between the lower and upper bounds depends on the ratio M/N. The less this ratio is, the less width of coherence variations is given by (28). When other terms are fixed, the right-hand side of (28) decreases with Psuc in the interval (M/N; 1). The value M/N gives the success probability in the trivial algorithm of random choice of items, when no amplification actually takes place. Over this value, any increase of the success probability will lead to decreasing range of allowed changes of quantum coherence with respect to the computational basis. To reach sufficiently high values of Psuc, an algorithm of amplitude amplification should somehow provide a coherence depletion. Note also that the right-hand side of (28) is saturated in some examples of amplitude amplification. Thus, this upper bound cannot be improved without considering additional parameters. We also obtained two inequalities for estimating the relative entropy of coherence from below. The formula (30) can be used, when we know in which of the sets and C the argument of max p(x) lies. This inequality may sometimes give a stronger lower bound on the relativeM entropyM of coherence than (28). In principle, the condition arg max p(x) is sufficiently natural for a good technique of amplitude amplification. We aim to build algorithms ∈M that provide values Psuc > 1/2 and even as closely to 1 as possible. Since M N M, the above condition is very plausible. Initially, this condition may be recorded a priori in the initial amplitude≤ − distribution. When we do not have data to apply the definition (29), the left-hand side of (28) should solely be used as an estimate from below. Among distance-based quantifiers, the ℓ1-norm of coherence is one of most intuitive [18]. On the other hand, it does not fulfill additional axioms proposed in [20]. It seems that trade-off relations between Cℓ1 (ρ) and Psuc(ρ) are enough complicated to formulate. Nevertheless, certain conclusions can be obtained in some particular cases. In this regard, we recall the corresponding discussion given in [23]. N−1 Suppose that amplitudes in the superposition ψ = x=0 cx x has the same absolute values separately for labels in and C, namely | i | i M M P C c = α x , c = α′ y . (34) | x| ∀ ∈M | y| ∀ ∈M Such amplitude distribution is typical during a performance of the standard Grover algorithm. In this case, we easily obtain α2 = P /M and (α′)2 = (1 P )/(N M). Combining the latter with the expression suc − suc − N−1 1+ C ψ = c c , ℓ1 | i | x y| x,y=0 X 7 we finally obtain
2 C ψ = MP + (N M)(1 P ) 1 . (35) ℓ1 | i suc − − suc − At first glance, a complementarity between p the coherencep quantifier and the success probability is not obvious. Following [23], we now take natural assumption M N. Then we rewrite (35) in the form ≪ C ψ + NP = N 1+ O M/N . (36) ℓ1 | i suc p If Psuc is close to 1 with 1 Psuc = O(M/N ), then the right-hand side of (36) can be replaced with N 1+ O(M/N) . Rescaling by the denominator− N, this relation becomes very similar to (21) and (27). This rescaling could be expected here, since the ℓ1-norm of coherence can increase up to N 1 [35], whereas the geometric coherence cannot exceed 1. In the context of original Grover’s algorithm, the relation− (36) was presented in [23]. We only note that the result (36) reflects no more than a “boxcar” distribution of amplitudes. In a similar manner, we could consider other situations of amplitude amplification with a priori knowledge. In general, however, complementarity relations between the ℓ1-norm of coherence and the success probability are not easy to formulate. Even the simplest non-trivial choice (34) has lead to (35). Despite of a simple structure of the density matrix, the formula (35) is complicated enough. Other coherence quantifiers should be preferred, when we focus on quantum coherence as a potential resource in amplitude amplification. The statements of Propositions 1 and 2 are formulated for an arbitrary quantum state. In this sense, they may be not optimal for analyzing the change of coherence just in amplitude amplification processes. Nevertheless, we can use them whenever the actual density matrix or some of its characteristics are known. We will present examples that demonstrate how quantum coherence may be anti-correlated with the success probability in amplitude amplification.
IV. SOME MODEL EXAMPLES OF COHERENCE CHANGES IN AMPLITUDE AMPLIFICATIONS
In this section, we will illustrate relations between coherence and the success probability within explicit model examples. We wish to study possible corollaries of the use of generalized blocks in amplitude amplification from the viewpoint of their bearing on quantum coherence. As we have recalled above, effects of generalized blocks can be described by means of the single ket η . In the following, we will make some model assumptions about η without a discussion of their origin. We rather| tryi to understand, whether quantum coherence is a key resource in| amplitudei amplification. These studies allow us to emphasize distinctions between cases, when marked and unmarked states are dealt with consistently or inconsistently. For convenience, we begin with the original Grover formulation [1].
A. Original Grover’s formulation
Let us set the initial amplitude distribution (12) and the rotation angles are β = γ = π. By M = , we mean the number of marked states, then N M = C . For all x and y C, the initial amplitudes|M| appear as − |M | ∈M ∈M 1 kx(0) = ℓy(0) = . (37) √N For the original algorithm, we have η = H 0 and η = η =1/√N, whence | i | i x y M N M W = , W = − . (38) k N ℓ N ′ ′ ′ ′ Further, one gives kx(t)= √N kx(t) and ℓy(t)= √Nℓy(t), including kx(0) = ℓx(0) = 1. In line with (A3) and (A4), the initial weighted averages appear as k ′(0) = ℓ ′(0) = 1, so that
′ ′ e e ∆kx = ∆ℓy =0 . (39) ′ ′ ′ ′ Then the formulas (A12) and (A13) merely say that kx(t) = kx(t) and ℓx(t) = ℓx(t) for all t. In the original formulation, we actually deal only with two different values of amplitudes. Substituting β = γ = π together with (38) into (A6), we get e e N 2M cos ω = W W = − . (40) ℓ − k N 8
Recall that we assume 1 M N/2. So, the parameter ω is defined by (A6) and ω (0; π/2]. It will be convenient to remember the formulas≤ ≤ ∈ M M sin2 ω/2= , cos2 ω/2=1 . (41) N − N
For β = γ = π, we further have ω± =2π ω due to (A5), whence the eigenvalues (A5) read as ± ±iω λ± = e = cos ω i sin ω . ± Further calculations lead to the following expressions, i exp(+iω/2) i exp( iω/2) ξ = − , ξ = − − . (42) 1 2 sin ω/2 2 2 sin ω/2 Using (A9) together with b = 1+cos ω and (42), we obtain exp(+iω/2) exp( iω/2) ξ = , ξ = − − . (43) 3 2cos ω/2 4 2cos ω/2
Since t is integer, we may take ω± = ω instead of ω± =2π ω. Due to (A7), one has ± ± sin ω(t +1/2) cos ω(t +1/2) k ′(t)= , ℓ ′(t)= . (44) sin ω/2 cos ω/2 ′ ′ ′ ′ As was already mentioned, kx(te)= k (t) and ℓy(t)= ℓ (t). Thus,e we finally write 1 sin ω(t +1/2) 1 cos ω(t +1/2) kx(t)= e e , ℓy(t)= , √N sin ω/2 √N cos ω/2 for all x and y C. The probabilities of interest appear as ∈M ∈M 2 M 2 2 Psuc(t)= kx(t) = sin ω(t +1/2) = sin ω(t +1/2) , (45) | | N sin2 ω/2 x∈M X N M 1 P (t)= ℓ (t) 2 = − cos2 ω(t +1/2) = cos2 ω(t +1/2) . (46) − suc | y | N cos2 ω/2 y∈MC X We wish to relate these probabilities with the coherence quantifiers. Since the state of the register is pure, its von Neumann entropy is zero. Further, the diagonal part of the density matrix reads as M entries ρ = diag k (t) 2,..., k (t) 2, ℓ (t) 2,..., ℓ (t) 2 . (47) diag | x | | x | | y | | y | z }| { N − M entries Hence, the geometric coherence obeys the complementarity relation| {z } P (t) C g(t) =1 suc , (48) g | i − M 2 2 whenever kx(t) ℓy(t) . Otherwise, the geometric coherence is equal to the 1 minus 1 Psuc(t) /(N M). Further, the| relative| ≥ entropy| | of coherence is calculated as − − C g(t) = M k (t) 2 ln k (t) 2 (N M) ℓ (t) 2 ln ℓ (t) 2 . (49) 1 | i − | x | | x | − − | y | | y | Due to (41), we can write M k (t) 2 = sin2 ω(t +1/2) = P (t) , (50) | x | suc (N M) ℓ (t) 2 = cos2 ω(t +1/2) =1 P (t) . (51) − | y | − suc Combining (49) with the last two formulas gives P (t) 1 P (t) C g(t) = P (t) ln suc 1 P (t) ln − suc . (52) 1 | i − suc M − − suc N M − The latter coincides with the right-hand side of (28). In contrast to (48), the formula (52) holds for all t. Therefore, the upper bound given by (28) is saturated in the original formulation of Grover’s search. For the given Psuc, the relative entropy of coherence reaches the maximal value approved by the right-hand side of (28). Grover’s search algorithm works so that any coherence decreasing is used in the most efficient way. 9
B. Marked and unmarked states are consistently amplified or decayed
Let us consider the case, when the terms U and s are such that marked states are amplified or decayed consistently. In this particular situation, some prior knowledge| i is available to users. We assume that amplitudes of the state η read as | i
M N M η = η , η = − η , (53) x MN y (N M)N r s − where x and y C. The angles β and γ are the same as in the original formulation. When amplitudes are balanced∈ with M respect∈ to M both and C, we describe them by a single non-integer parameter M > 0. An efficiency M M η is increased with Mη >M and decreased with Mη Combining (53) with (A3) and (A4) directly gives k ′ (t)= k ′(t) for x and ℓ ′ (t)= ℓ ′(t) for y C. Rewriting x ∈M x ∈M (A6) with Mη instead of M, we have the parameter ωη such that e e M M sin2 ω /2= η , cos2 ω /2=1 η . (55) η N η − N The latter should be used simultaneously with (41). The initial amplitude distribution (37) is the same as in the original formulation. Correspondingly to (53), we have MN (N M)N k ′ (t)= k (t) , ℓ ′ (t)= − ℓ (t) . (56) x M x y N M y s η s − η Together with kx(0) = ℓy(0) = 1/√N, one then obtains ′ ′ M sin ω/2 kx(0) = k (0) = = , (57) sMη sin ωη/2 e N M cos ω/2 ℓ ′ (0) = ℓ ′(0) = − = , (58) x N M cos ω /2 s − η η e where the parameter ω is again defined by (40). For β = γ = π, we further have ω± =2π ω . As calculations show, the formulas (42) and (43) are replaced with ± η i exp(+iω/2) i exp( iω/2) ξ1 = − , ξ2 = − − , (59) 2 sin ωη/2 2 sin ωη/2 exp(+iω/2) exp( iω/2) ξ3 = , ξ4 = − − . (60) 2cos ωη/2 2cos ωη/2 Similarly to (44), we obtain the averaged amplitudes sin ωηt + ω/2 cos ωηt + ω/2 k ′(t)= , ℓ ′(t)= . (61) sin ωη/2 cos ωη/2 e ′ ′ e ′ ′ C These formulas also represent the amplitudes kx(t) = k (t) for x and ℓx(t) = ℓ (t) for y . Substituting t = 0, they reduce to (57) and (58). Combining (56) with (61) immediat∈ M ely gives ∈ M e e 2 2 sin ωηt + ω/2 cos ωηt + ω/2 k (t) 2 = , ℓ (t) 2 = . | x | M | y | N M − Calculating the probabilities then results in P (t) = sin2 ω t + ω/2 , 1 P (t) = cos2 ω t + ω/2 . (62) suc η − suc η 10 2 2 When kx(t) ℓy(t) , the geometric coherence again obeys (48). We also note that relations (50) and (51) are still valid. The| relative| ≥ | entropy| of coherence is again connected with the probabilities by (52). The success probability is first maximized, when ωηt + ω/2 becomes as close to π/2 as possible. As t is integer, we take one of the two numbers π ω π ω − , − . (63) 2 ω 2 ω η η With growth of Mη >M in (53), the key parameter ωη >ω also increases. It is natural that both M and Mη are very small in comparison with N. Then we approximately write ωη 2 Mη/N. Hence, the optimal measurement time ≈ −1/2 can be estimated. When Mη grows and other parameters are fixed, thep integers (63) decrease proportionally to Mη . By ωη, one characterizes a rate of amplitude amplification. Increasing Mη due to conducive prior knowledge, one reduces the number of iterations required for maximizing the success probability. For the fixed M and N, we herewith accelerate an algorithm work. When Mη