Non-Boolean Quantum Amplitude Ampification and Quantum Mean Estimation Systems and Methods

Abstract

Generalizations of quantum amplitude amplification and amplitude estimation algorithms work with non-boolean oracles (by way of definition, the action of a non-boolean oracle Uφ on an eigenstate |x is to apply a state-dependent phase-shift φ(x); unlike boolean oracles, the eigenvalues exp(iφ(x)) of a non-boolean oracle are not restricted to be ±1). The non-boolean amplitude amplification algorithm preferentially amplifies the amplitudes of the eigenstates based on the value of φ(x). Starting from a given initial superposition state |ψ0, the basis states with lower values of cos(φ) are amplified at the expense of the basis states with higher values of cos(φ). The non-boolean quantum mean estimation algorithm uses quantum phase estimation to estimate the expectation ψ0|Uφ|ψ0 (i.e., the expected value of exp(iφ(x)) for a random x sampled by making a measurement on |ψ0). The quantum mean estimation algorithm offers a quadratic speedup over its counterpart boolean algorithm known in the art.

Description

FIELD OF THE INVENTION

The present invention relates generally to quantum computing technology. More particularly, this invention pertains to devices, systems, and associated methods for achieving computational speed increases in quantum algorithms.

BACKGROUND OF THE INVENTION

Grover's algorithm is a quantum search algorithm for finding the unique input x_good that satisfies Equation (1), as follows:

ƒ_bool(x_good)=1,  (1)

for a given boolean function ƒ_bool: {0, 1, . . . , N−1} →{0,1}. Such an input x_good satisfying this boolean function is referred to as the “winning” input of ƒ_bool. Grover's algorithm has also been adapted to work with boolean functions with multiple winning inputs, where the goal is to find any one of the winning inputs.

An important generalization of Grover's algorithm is the amplitude amplification algorithm in which the function ƒ_bool is accessed through a boolean quantum oracle Û_ƒbool that acts on the orthonormal basis states |0, . . . , |N−1 as follows (Equation (2)):

$\begin{array}{cc}{\hat{U}}_{f_{\mathrm{bool}}}❘x\rangle =\{\begin{array}{cc}-❘x\rangle ,& \mathrm{if}{f}_{\mathrm{bool}}\left(x\right)=1,\\ +❘x\rangle ,& \mathrm{if}{f}_{\mathrm{bool}}\left(x\right)=0.\end{array}& \left(2\right)\end{array}$

In this way, the oracle marks the winning states by flipping their phase (that is, shifting the phase by π). Given a superposition state |ψ₀, the goal of the amplitude amplification algorithm is to amplify the amplitudes (in the superposition state) of the winning states. The algorithm accomplishes this iteratively by initializing a quantum system in the state |ψ₀ and performing the operation S_ψ0 Û_ƒbool on the system during each iteration, where

S _{ψ 0} ≡2|ψ₀ψ₀|−I.  (3)

As shown in Equation (3), I is the identity operator. Performing a measurement on the system after the iterative amplification process results in one of winning states with high probability. Grover's original algorithm is a special case of the amplitude amplification algorithm, where a) the uniform superposition state |s given by Equation (4):

$\begin{array}{cc}❘s\rangle =\frac{1}{\sqrt{N}}\sum _{x=0}^{N-1}❘x\rangle & \left(4\right)\end{array}$

is used as the initial state |ψ₀ of the system, and b) the number of winning inputs is exactly one.

Closely related to the amplitude amplification algorithm is the amplitude estimation algorithm, which combines features from the amplitude amplification algorithm and a quantum phase estimation (QPE) algorithm to estimate the probability that making a measurement on the initial state |ψ₀ will yield a winning input. If the uniform superposition state |s is used as |ψ₀, the amplitude estimation algorithm can help estimate the number of winning inputs of ƒ_bool (Note: This special case is also referred to as the quantum counting algorithm).

The amplitude amplification algorithm and the amplitude estimation algorithm have a wide range of applications and are important primitives that feature as subroutines employed by other quantum algorithms. The amplitude amplification algorithm can be used to find a winning input to ƒ_bool with (√{square root over (N)}) queries of the quantum oracle, regardless of whether the number of winning states is a priori known or unknown (Note: The (√{square root over (N)}) and (N) scalings for the quantum and classical algorithms, respectively, hold assuming that the number of winning states does not scale with N). This represents a quadratic speedup over classical algorithms, which typically require (N) evaluations of the function ƒ_bool. Similarly, the amplitude estimation algorithm also offers a quadratic speedup over the corresponding classical approaches. The quadratic speedup due to the amplitude amplification algorithm has been shown to be optimal for oracular quantum search algorithms.

A limitation of known amplitude amplification and estimation algorithms is that they work only with boolean oracles, which classify the basis states as good and bad. In situations where one is interested in using these algorithms in the context of a non-boolean function of the input x, a typical approach is to create a boolean oracle from the non-boolean function by using a threshold value of the function as a decision boundary. That is, the winning states are the ones for which the value of the function is, say, less than the chosen threshold value. In this way, the problem at hand may be adapted to work with the standard amplitude amplification and estimation algorithms.

Accordingly, a need exists for a solution to at least one of the aforementioned challenges in increasing the computation speed of widely applicable quantum algorithms. For instance, an established need exists for adaptation of certain primitive quantum algorithms to work directly with non-boolean functions.

This background information is provided to reveal information believed by the applicant to be of possible relevance to the present invention. No admission is necessarily intended, nor should be construed, that any of the preceding information constitutes prior art against the present invention.

SUMMARY OF THE INVENTION

With the above in mind, embodiments of the present invention are related to quantum amplitude amplification and amplitude estimation algorithms to work with non-boolean oracles. By way of definition, the action of a non-boolean oracle U_φ on an eigenstate |x is apply a state-dependent, real-valued phase-shift φ(x). Unlike boolean oracles, the eigenvalues exp(iφ(x)) of a non-boolean oracle are not restricted to be ±1.

In one embodiment of the present invention, a non-boolean amplitude amplification algorithm, starting from an initial superposition state |ψ₀, preferentially amplifies the amplitudes of the basis states based on the value of cos(φ). An objective of the algorithm is to preferentially amplify the states with lower values of cos(φ(x)). The algorithm is iterative in nature such that, after K iterations, the probability for a measurement of the system to yield x (namely p_K(x)) differs from the original probability p₀(x) by a factor that is linear in cos(φ(x)). The coefficient −λ_K of this linear dependence controls the degree (and direction) of the preferential amplification.

More specifically, embodiments of the present invention (in one or more of method, system, and/or device form) may include the following steps:

(1) Initialize a two-register system in the |Ψ₀ state:

(2) Perform K iterations: During odd iterations, act on an input basis state using a selective phase-flip unitary operator circuit S_Ψ0 and a controlled unitary operator circuit U_φ. During even iterations, act on an input basis state using the selective phase-flip unitary operator circuit S_Ψ0 and a controlled inverse unitary operator circuit U_φ^†; and

(3) After the K iterations, measure the ancilla in the 0/1 basis.

Up to a certain number of iterations, the iterative steps may be designed to amplify the amplitude of the basis states |0,x and |1,x with lower values of cos(φ(x)). The measurement of the ancilla at the end of the iterations may be performed simply to ensure that the two registers are not entangled in the final state.

In another embodiment of the present invention, a quantum mean estimation algorithm uses QPE as a subroutine in order to estimate the expectation of U_φ under |ψ₀ (i.e., ψ₀|ψ_φ|ψ₀). The algorithm offers a quadratic speedup over the classical approach of estimating the expectation, as a sample mean over randomly sampled inputs.

More specifically, embodiments of the present invention (in one or more of method, system, and/or device form) may include the following steps:

(1) Perform the QPE algorithm with a) the two-register unitary operator under consideration, and b) the superposition state |Ψ₀ in place of the eigenstate required by the QPE algorithm as input. Let the output of this step, appropriately scaled to be an estimate of the phase angle in the range [0,2π), be {circumflex over (ω)}.

(2) Return cos({circumflex over (ω)}) as the estimate for cos(θ) (i.e., the real part of E_ψ0[e^iφ]).

These and other objects, features, and advantages of the present invention will become more readily apparent from the attached drawings and the detailed description of the preferred embodiments, which follow.

BRIEF DESCRIPTION OF THE DRAWINGS

The preferred embodiments of the invention will hereinafter be described in conjunction with the appended drawings provided to illustrate and not to limit the invention, where like designations denote like elements, and in which:

FIG. 1 is a schematic diagram depicting an exemplary quantum circuit implementing a selective phase-flip unitary operator for non-boolean amplitude amplification according to an embodiment of the present invention;

FIG. 2 is a schematic diagram depicting an exemplary quantum circuit implementing a first conditional oracle call unitary operator for a non-boolean amplitude amplification algorithm according to an embodiment of the present invention;

FIG. 3 is a schematic diagram depicting an exemplary quantum circuit implementing a second conditional oracle call unitary operator for a non-boolean amplitude amplification algorithm according to an embodiment of the present invention;

FIG. 4 is schematic diagram depicting a quantum circuit for a first exemplary non-boolean amplification algorithm according to an embodiment of the present invention;

FIG. 5 is pseudocode for a non-boolean amplification algorithm according to FIG. 4;

FIG. 6 is a state diagram depicting stage evolution of a two-register system through iterations of the non-boolean amplitude amplification algorithm of FIG. 5;

FIG. 7 is schematic diagram depicting a quantum circuit for a second exemplary non-boolean amplification algorithm according to an embodiment of the present invention;

FIG. 8 is a non-boolean amplification algorithm according to FIG. 7;

FIG. 9 is a schematic diagram depicting a quantum circuit for a non-boolean mean estimation algorithm according to an embodiment of the present invention;

FIG. 10 depicts a graph of oracle function and initial state of a toy example for quantum simulation;

FIG. 11 depicts a graph of quantum simulation results from application of the non-boolean amplitude amplification algorithm of FIG. 4 to the toy example of FIG. 10; and

FIG. 12 depicts graphs of quantum simulation results from application of the non-boolean mean estimation algorithm of FIG. 9 to the toy example of FIG. 10.

Like reference numerals refer to like parts throughout the several views of the drawings.

DETAILED DESCRIPTION OF THE INVENTION

The present invention will now be described more fully hereinafter with reference to the accompanying drawings, in which preferred embodiments of the invention are shown. This invention may, however, be embodied in many different forms and should not be construed as limited to the embodiments set forth herein. Rather, these embodiments are provided so that this disclosure will be thorough and complete, and will fully convey the scope of the invention to those skilled in the art.

Although the following detailed description contains many specifics for the purposes of illustration, anyone of ordinary skill in the art will appreciate that many variations and alterations to the following details are within the scope of the invention. Accordingly, the following embodiments of the invention are set forth without any loss of generality to, and without imposing limitations upon, the claimed invention.

As used herein, the word “exemplary” or “illustrative” means “serving as an example, instance, or illustration.” Any implementation described herein as “exemplary” or “illustrative” is not necessarily to be construed as preferred or advantageous over other implementations. All of the implementations described below are exemplary implementations provided to enable persons skilled in the art to make or use the embodiments of the disclosure and are not intended to limit the scope of the disclosure, which is defined by the claims.

Furthermore, in this detailed description, a person skilled in the art should note that quantitative qualifying terms such as “generally,” “substantially,” “mostly,” and other terms are used, in general, to mean that the referred to object, characteristic, or quality constitutes a majority of the subject of the reference. The meaning of any of these terms is dependent upon the context within which it is used, and the meaning may be expressly modified.

Referring initially to FIGS. 1-9, a non-boolean quantum amplitude amplification algorithm, as also a quantum mean estimation algorithm based in the non-boolean quantum amplitude amplification algorithm, both according to embodiments of the present invention are now described in detail. Throughout this disclosure, the present invention may be referred to as a family of non-boolean quantum amplitude amplification algorithms, a family of non-boolean quantum algorithms, a non-boolean quantum algorithm, a non-boolean quantum method, a non-boolean quantum oracle, a non-boolean quantum system, a method, an oracle, and/or a system. Those skilled in the art will appreciate that this terminology is only illustrative and does not affect the scope of the invention. For instance, the present invention may just as easily relate to an instantiation of an object from a library of non-boolean quantum oracles.

Generally speaking, the present invention is a generalization of amplitude amplification and estimation algorithms to work with quantum oracles for non-boolean functions. Hereinafter, the qualifiers “boolean” and “non-boolean” will be used to distinguish embodiments of the present invention from known boolean versions of the amplitude amplification algorithm and their related applications.

Oracle for a Non-Boolean Function:

The behavior of the boolean quantum oracle Û_ƒbool of Equation (2) may be generalized to non-boolean functions by allowing the oracle to perform arbitrary phase-shifts on the different basis states. More concretely, let φ: {0, 1, . . . , N−1} → be a real-valued function, and let U_φ be a quantum oracle given by Equation (5), as follows:

$\begin{array}{cc}{U}_{\phi }\equiv \sum _{x=0}^{N-1}{e}^{i\phi \left(x\right)}❘x\rangle \langle x❘.& \left(5\right)\end{array}$

The actions of the oracle U_φ and its inverse U_φ^† on the basis states |0, . . . , |N−1 may be given by Equations (6) and (7), as follows:

U _φ |x)=e^+iφ(x)|x=[cos(φ(x))+i sin(φ(x))]|x,  (6)

U _φ ^† |x)=e^−iφ(x)|x=[cos(φ(x))−i sin(φ(x))]|x.  (7)

Goal of the Non-Boolean Amplitude Amplification Algorithm:

Given an oracle U_φ and an initial state |ψ₀, a goal of the non-boolean amplitude amplification algorithm of the present invention may be to preferentially amplify the amplitudes of the basis states |x with lower values of cos(φ(x)), at the expense of the amplitude of states with higher values of cos(φ(x)). Depending on the context in which the algorithm is to be used, a different function of interest ƒ (which is intended to guide the amplification) may be appropriately mapped onto the function cp. For example, and without limitation, if the range of ƒ is [0,1] and one intends to amplify the states with higher values of ƒ, then options for formulating the problem in terms of co include the following:

φ(x)=πƒ(x) or φ(x)=arccos(1−2ƒ(x)).  (8)

In both cases of Equations (8), cos(φ) is monotonically decreasing in ƒ.

The connection between the boolean and non-boolean amplitude amplification algorithms may be seen as follows: If either of the two options in Equations (8) is used to map a boolean function ƒ_bool onto φ_bool, then

$\begin{array}{cc}{\phi }_{\mathrm{bool}}\left(x\right)=\{\begin{array}{cc}\pi ,& \mathrm{if}{f}_{\mathrm{bool}}\left(x\right)=1,\\ 0,& \mathrm{if}{f}_{\mathrm{bool}}\left(x\right)=0.\end{array}& \left(9\right)\end{array}$

In the case of Equation (9), the oracle U_φ and its inverse U_φ^† both reduce to a boolean oracle as follows (Equation 10):

U _{φ bool} =U _{φ bool} ^† =Û _{ƒ bool} .  (10)

Congruently, the task of amplifying (the amplitude of) the states with lower values of cos(φ) may align with the task of amplifying the winning states |x with ƒ_bool(x)=1.

Goal of the Quantum Mean Estimation Algorithm:

Given a generic unitary operator U and a state |ψ₀, a goal of the quantum mean estimation algorithm may be to estimate the quantity of Equation (11):

ψ₀|U|ψ₀.  (11)

The task of estimating the quantity of Equation (11) may be phrased in terms of the oracle U_φ, as estimating the expectation of the eigenvalue e^iφ(x) for a state |x chosen randomly by making a measurement on the superposition state |ψ₀. The connection between the two tasks may be seen in Equation (12), as follows:

$\begin{array}{cc}\langle {\psi }_0❘U_{\phi }❘{\psi }_0\rangle =\sum _{x=0}^{N-1}\langle {\psi }_0❘e^{i\phi \left(x\right)}❘x\rangle \langle x❘{\psi }_0\rangle =\sum _{x=0}^{N-1}❘\langle x❘{\psi }_0\rangle {❘}^2{e}^{i\phi \left(x\right)},& \left(12\right)\end{array}$

wherein |x|ψ₀|² is the probability for a measurement on |ψ₀ to yield x.

The only difference between (a) estimating ψ₀|U_φ|ψ₀ for an oracle U_φ, and (b) estimating ψ₀|U|ψ₀ for a generic unitary operator U is that {|ψ₀|U_φ|ψ₀, . . . , |N−1)} is known beforehand to be an eigenbasis of U_φ. On the other hand, the eigenstates of a generic unitary operator U may be a priori unknown. However, as described in detail hereinbelow, the mean estimation algorithm does not use the knowledge of the eigenstates and, therefore, may be applicable for generic unitary operators U as well.

As described hereinabove, the mean estimation algorithm of the present invention is a generalization of known boolean amplitude estimation algorithm(s). Regarding the connection between the respective tasks of these algorithms, note that the eigenvalues of a boolean oracle may be either +1 or −1, and the expectation of the eigenvalue under |ψ₀ is directly related to the probability of a measurement yielding a winning state with eigenvalue −1. This probability is precisely the quantity estimated by the non-boolean amplitude estimation algorithm of the present invention.

Non-Boolean Amplitude Amplification Algorithm:

Setup and Notation: Various embodiments of a non-boolean amplitude amplification algorithm as described herein may include not only a quantum system, or qubit(s), as input to a quantum oracle, but also may employ one extra ancilla qubit. For example, and without limitation, let a quantum system used in certain embodiments of the present algorithm comprise two quantum registers. The first register may contain the lone ancilla qubit, and the second register (input qubit) may be acted upon by the quantum oracle.

The notations |a⊗|b and |a,b may both refer to a state where the two registers are unentangled, with the first register in state |a and the second register in state |b. The tensor product notation ⊗ may also be used to combine operators that act on the individual registers into operators that simultaneously act on both registers. Such two-register operators may be represented by boldface symbols (e.g., S_Ψ0, U_φ, I). Likewise, boldface symbols may be used to represent the states of the two-register system in the bra-ket notation (e.g., |Ψ₀). As used herein, any state written in the bra-ket notation (e.g., |ψ) will be unit normalized (i.e., normalized to 1). The dagger notation (†) may be used to denote the Hermitian conjugate of an operator, which is also the inverse for a unitary operator.

Unless otherwise specified herein, {|0, |1, . . . , |N−1} may be used as the basis for (that is, as the state space of) the second register. Any measurement of the second register may refer to measurement in this basis. Likewise, unless otherwise specified,

{|0,0,|0,1, . . . ,|0,N−1}∪{|1,0,|1,1, . . . ,|1,N−1}  (13)

may be used as the basis for the two-register system.

In Equation (14) below, let |ψ₀ be the initial state of the second register from which the amplification process is to begin, and let A₀ be the unitary operator that changes the state of the second register from |0 to |ψ₀ (Note: Assumed here, only for notational convenience, is that there exists a state |0| which is simultaneously an eigenstate of U_φ, as well as a special, easy-to-prepare state of the second register; a person of skill in the art will immediately recognize the algorithms described herein may be modified to work even without this assumption):

$\begin{array}{cc}|{\psi }_0\rangle \equiv A_0❘0\rangle \equiv \sum _{x=0}^{N-1}{a}_0\left(x\right)❘x\rangle ,& \left(14\right)\end{array}$ $\mathrm{such}\mathrm{that}$ $\sum _{x=0}^{N-1}|a_0\left(x\right){|}^2=1,$

where a₀(x) is the initial amplitude of the basis state |x.

The algorithm introduced hereinbelow may initialize the ancilla (i.e., first register) in the |+ state given by Equation (15), as follows:

$\begin{array}{cc}❘+\rangle =\frac{❘0\rangle +❘1\rangle }{\sqrt{2}}.& \left(15\right)\end{array}$

Anticipating this initialization, let the two-register state |Ψ₀ be defined as in Equation (16), as follows:

$\begin{array}{cc}\begin{array}{c}❘{\Psi }_0\rangle \equiv ❘+,{\psi }_0\rangle =\frac{❘0,{\psi }_0\rangle +❘1,{\psi }_0\rangle }{\sqrt{2}}\\ =\frac{1}{\sqrt{2}}\sum _{x=0}^{N-1}{a}_0\left(x\right)[❘0,x\rangle +❘1,x\rangle ]\end{array}.& \left(16\right)\end{array}$

Required Unitary Operations: The following unitary operations may be used in certain embodiments of a generalized amplitude amplification algorithm of the present invention:

Selective Phase-Flip Operator. Let the two-register unitary operator S_Ψ0 be defined as in Equation (17):

$\begin{array}{cc}\begin{array}{c}{S}_{{\Psi }_0}\equiv 2❘{\Psi }_0\rangle \langle {\Psi }_0❘-I\\ =2❘+,{\psi }_0\rangle \langle +,{\psi }_0❘-I\end{array},& \left(17\right)\end{array}$

where I is the two-register identity operator. S_Ψ0 may leave the state |Ψ₀ unchanged and may flip the phase of any state orthogonal to |Ψ₀. S_Ψ0 may be the two-register generalization of S_ψ0 used in counterpart boolean amplitude amplification algorithm(s). From Equations (14) and (16), it follows (in Equation (18)) that

|Ψ₀=[H⊗A₀]|0,0,  (18)

where H is the Hadamard transform. Thus, S_Ψ0 may be expressed as in Equation (19):

S _{Ψ 0} =[H⊗A₀][2|0,00,0|−I][H⊗A₀^†].  (19)

Referring initially to FIG. 1, an exemplary quantum circuit 100 implementing non-boolean amplitude amplification according to an embodiment of the present invention will now be described in detail. A quantum system of circuit 100 may comprise a first register (ancilla qubit) 110 and a second register (input qubit) 120. The second register 120 may be configured to be acted upon by a quantum oracle 130. As illustrated, Equation (19) is used to drive an implementation of circuit 130 for S_Ψ0, provided one has access to the quantum circuits that implement A₀ (unitary operator 142) and A₀^† (inverse unitary operator 144). Hadamard transforms 150 are denoted as H.

Conditional Oracle Calls: Let the two-register unitary operator U be defined as in Equation (20):

U _φ≡|0|⊗U_φ+|11|⊗U_φ^†.  (20)

This operator's action on the basis states of the two-register system is given by Equations (21) and (22):

U _φ|0,x)=e^+iφ(x)|0,x,  (21)

U _φ|1,x)=e^−iφ(x)|1,x.  (22)

If the ancilla is in state |0, U_φ acts U_φ on the second register. On the other hand, if the ancilla is in state |1, U_φ acts U_φ^† on the second register. The inverse of U_φ is given by Equation (23):

U _φ ^†=|00|⊗U_φ^†+|11|⊗U_φ,  (23)

and the action of U_φ^† on the basis states is given by Equations (24) and (25):

U _φ ^†|0,x=e^−iφ(x)|0,x,  (24)

U _φ ^†|1,x=e^+iφ(x)|1,x.  (25)

The amplitude amplification algorithm for non-boolean functions may involve calls to both U_φ and U_φ^†.

Referring now to FIGS. 2 and 3, exemplary quantum circuits 200 and 300 each implementing non-boolean amplitude amplification according to an embodiment of the present invention will now be described in detail. FIG. 2 depicts a circuit implementation 200 of U_φ using a) ancilla 210 serving as a control qubit, b) bit-flip (or Pauli-X) gates 240 denoted as X, and c) second register 220 acted upon by controlled U_φ and U_φ^† operations (i.e., oracles) 232, 234, respectively. FIG. 3 depicts a circuit implementation 330 of U_φ^† using a) ancilla 310 serving as a control qubit, b) bit-flip (or Pauli-X) gates 340 denoted as X, and c) second register 320 acted upon by controlled U_φ and U_φ^† operations 332, 334, respectively.

Algorithm Description:

Referring now to FIGS. 4 and 5, an exemplary quantum circuit 400 and pseudocode 500 implementing a non-boolean amplitude amplification quantum algorithm according to an embodiment of the present invention will now be described in detail. For example, and without limitation, the amplitude amplification algorithm for non-boolean functions is iterative and may comprise the following steps:

(1) Initialize a two-register system in the |Ψ₀ state;

(2) Perform K iterations: During odd iterations 402, apply operations 430, 432 denoted S_Ψ0 U_φ on the system. During even iterations 404, apply operations 430, 434 denoted S_Ψ0 U_φ^† on the system; and

(3) After the K iterations, measure 440 the ancilla (first register 410) in the 0/1 basis.

Up to a certain number of iterations, the iterative steps 402, 404 may be designed to amplify the amplitude of the basis states |0,x and |1,x with lower values of cos(φ(x)). The measurement 440 of the ancilla 410 at the end of the algorithm may be performed simply to ensure that the two registers 410, 420 are not entangled in the final state 450 of the system. The specification of K (i.e., the number of iterations to perform) is included in an analysis of the present algorithm described hereinbelow.

Connection to Known Boolean Amplitude Amplification Algorithm:

From Equations (20) and (23), a person of skill in the art will immediately recognize that for the boolean oracle case given by φ(x)=πƒ_bool(x), U_φ and U_φ^† both reduce to I⊗Û_ƒbool, where Û_ƒbool is the oracle used in known boolean amplitude amplification algorithm(s). Furthermore, if a first register is in the |+ state, from Equations (3) and (17), the action of S_Ψ0 is given by Equation (26):

S _{Ψ 0} [|+⊗|ψ]=|+⊗[S_ψ0|ψ].  (26)

Note that the first register is unaffected here. Thus, for the boolean oracle case, the algorithm 400, 500 of FIGS. 4 and 5 may reduce to simply acting S_ψ0 Û_ƒbool on the second register 420 during each iteration 402, 404; the ancilla qubit 410 remains untouched and unentangled from the second register 420. In this way, algorithm 400, 500 is a generalization of known boolean amplitude amplification algorithm(s) described hereinabove.

The two key differences of the generalized algorithm from the boolean algorithm, apart from the usage of a non-boolean oracle, may be as follows:

(1) The addition of the ancilla 410, which doubles the dimension of the state space of the system, and

(2) Alternating between using U_φ and U_φ^† (operations 432 and 434, respectively) during the odd iterations 402 and even iterations 404.

Analysis of the Non-Boolean Amplitude Amplification Algorithm:

Still referring to FIGS. 4 and 5, let |Ψ_k be the state of the two-register system after k=0, 1, . . . , K iterations of the amplitude amplification algorithm 400 (but before the measurement 440 of the ancilla 410). For k>0, |Ψ_k may be recursively written as in Equation (27):

$\begin{array}{cc}❘{\Psi }_k\rangle \equiv \{\begin{array}{cc}{S}_{{\Psi }_0}{U}_{\phi }❘{\Psi }_{k-1}\rangle ,& \mathrm{if}k\mathrm{is}\mathrm{odd},\\ S_{{\Psi }_0}{U}_{\phi }^{†}❘{\Psi }_{k-1}\rangle ,& \mathrm{if}k\mathrm{is}\mathrm{even}\end{array}.& \left(27\right)\end{array}$

Let ã_k (0,x) and ã_k(1,x) be the normalized amplitudes of the basis states |0,x and |1,x, respectively, in the superposition |Ψ_k.

$\begin{array}{cc}❘{\Psi }_k\rangle \equiv \sum _{x=0}^N[{\tilde{a}}_k\left(0,x\right)❘0,x\rangle +{\tilde{a}}_k\left(1,x\right)❘1,x\rangle ].& \left(28\right)\end{array}$

In the initial state |Ψ₀, the amplitudes ã₀(0,x) and ã₀(1,x) are both given, from Equation (16), by the following (Equation (29)):

$\begin{array}{cc}{\tilde{a}}_0\left(0,x\right)={\tilde{a}}_0\left(1,x\right)=\frac{a_0\left(x\right)}{\sqrt{2}}.& \left(29\right)\end{array}$

Let the parameter θ∈[0,π] be implicitly defined by Equation (30):

$\begin{array}{cc}\cos \left(\theta \right)\equiv \sum _{x=0}^{N-1}{❘a_0\left(x\right)❘}^2\cos \left(\phi \left(x\right)\right).& \left(30\right)\end{array}$

cos(θ) is the expected value of cos(φ(x)) over bitstrings x sampled by measuring the state |ψ₀.

Let the two-register states |αand |β be defined as follows (Equations (31) and (32)):

|α≡U_φ|Ψ₀,  (31)

|β≡U_φ^†|Ψ₀.  (32)

These register states may be used to track the evolution of the system through the iterative steps 402, 404 of algorithm 400, 500. Using Equations (16), (20), and (23), |α and |0 may be written as follows (Equations (33) and (34)):

$\begin{array}{cc}❘\alpha \rangle =\frac{1}{\sqrt{2}}\sum _{x=0}^{N-1}{a}_0\left(x\right)[e^{i\phi \left(x\right)}❘0,x\rangle +e^{-i\phi \left(x\right)}❘1,x\rangle ],& \left(33\right)\end{array}$ $\begin{array}{cc}❘\beta \rangle =\frac{1}{\sqrt{2}}\sum _{x=0}^{N-1}{a}_0\left(x\right)[e^{-i\phi \left(x\right)}❘0,x\rangle +e^{i\phi \left(x\right)}❘1,x\rangle ].& \left(34\right)\end{array}$

Note that θ, |α, and |β are all implicitly dependent on the function φ and the initial state |ψ₀. For notational convenience, these dependencies are not explicitly indicated.

The First Iteration:

After one iterative step, the system may be in state |Ψ₁ given by Equation (35):

|Ψ₁=S_Ψ0U_φ|Ψ₀.  (35)

Using Equations (17) and (31), this state may be written as follows (Equation (36)):

|Ψ₁=S_Ψ0|α=2Ψ₀|αΨ₀−|α.  (36)

From Equations (16), (30), and (33), it follows that

$\begin{array}{cc}\begin{array}{c}\langle {\Psi }_0❘\alpha \rangle =\frac{1}{2}\sum _{x=0}^{N-1}{❘a_0\left(x\right)❘}^2\left[e^{i\phi \left(x\right)}\langle 0,x❘0,x\rangle +e^{-i\phi \left(x\right)}\langle 1,x❘1,x\rangle \right]\\ =\frac{1}{2}\sum _{x=0}^{N-1}{❘a_0\left(x\right)❘}^2\left[e^{i\phi \left(x\right)}+e^{-i\phi \left(x\right)}\right]\\ =\frac{1}{2}\sum _{x=0}^{N-1}{❘a_0\left(x\right)❘}^{2}2\cos \left(\phi \left(x\right)\right)=\cos \left(\theta \right)\end{array}.& \left(37\right)\end{array}$

Note from Equation (37) that Ψ₀|α is real-valued. Key to the functioning of the algorithm 400, 500, the motivation behind adding an ancilla qubit 410 (effectively doubling the number of basis states) is precisely to make Ψ₀|α real-valued.

From Equations (36) and (37), |Ψ₁ may be written as follows (Equation (38)):

|Ψ₁=2 cos(θ)|Ψ₀−|α.  (38)

From Equations (16) and (33), the amplitude ã₁(0,x) of the basis state |0,x in the superposition |Ψ₁ may be written as follows (Equation (39)):

$\begin{array}{cc}{\tilde{a}}_1\left(0,x\right)=\langle 0,x❘{\Psi }_1\rangle =\frac{a_0\left(x\right)}{\sqrt{2}}\left[2\cos \left(\theta \right)-e^{i\phi \left(x\right)}\right].& \left(39\right)\end{array}$

Likewise, the amplitude ã₁(1,x) of the basis state |1,x may be written as follows (Equation (40)):

$\begin{array}{cc}{\tilde{a}}_1\left(1,x\right)=\langle 1,x❘{\Psi }_1\rangle =\frac{a_0\left(x\right)}{\sqrt{2}}\left[2\cos \left(\theta \right)-e^{-i\phi \left(x\right)}\right].& \left(40\right)\end{array}$

Equations (39) and (40) show that, after one iterative step, the amplitudes of |0,x and |1,x have acquired factors of [2 cos(θ)−e^iφ(x)] and [2 cos(θ)−e^−iφ(x)], respectively. Now, Equation (41) shows that, if cos(θ) is positive, the magnitude of the “amplitude amplification factor” after one iteration is monotonically decreasing in cos(φ):

$\begin{array}{cc}\frac{{❘{\tilde{a}}_1\left(0,x\right)❘}^2}{{❘{\tilde{a}}_0\left(0,x\right)❘}^2}=\frac{{❘{\tilde{a}}_1\left(1,x\right)❘}^2}{{❘{\tilde{a}}_0\left(1,x\right)❘}^2}=4{\cos }^2\left(\theta \right)+1-4\cos \left(\theta \right)\cos \left(\phi \left(x\right)\right).& \left(41\right)\end{array}$

Such preferential amplification of states with lower values of cos(φ) is precisely what the algorithm 400, 500 set out to do.

Note that this monotonicity property relies on Ψ₀|α being real valued in Equation (36). If Ψ₀|α is complex, with a phase δ∉{0,π}, then the amplification will be monotonic in cos(φ−δ), which does not meet the present goal of the algorithm 400, 500. The case where Ψ₀|α is not real-valued is explored further hereinbelow.

Identities to Track the Subsequent Iterations:

Equations (16), (17), (20), (23), (31), and (32) may be used to derive the identities of Equations (42), which capture the actions of the operators S_Ψ0, U_φ, and U_φ^† on the states IWO), |α, and |β:

$\begin{array}{cc}{S}_{{\Psi }_0}❘{\Psi }_0\rangle =❘{\Psi }_0\rangle ,& \left(42\right)\end{array}$ $S_{{\Psi }_0}❘\alpha \rangle =2\cos \left(\theta \right)❘{\Psi }_0\rangle -❘\alpha \rangle ,$ $S_{{\Psi }_0}❘\beta \rangle =2\cos \left(\theta \right)❘{\Psi }_0\rangle -❘\beta \rangle ,$ $U_{\phi }❘{\Psi }_0\rangle =❘\alpha \rangle ,$ $U_{\phi }❘\alpha \rangle =\frac{1}{\sqrt{2}}\sum _{x=0}^{N-1}{a}_0\left(x\right)[e^{2i\phi \left(x\right)}❘0,x\rangle +e^{-2i\phi \left(x\right)}❘1,x\rangle ],$ $U_{\phi }❘\beta \rangle =❘{\Psi }_0\rangle ,$ $U_{\phi }^{†}❘{\Psi }_0\rangle =❘\beta \rangle ,$ $U_{\phi }^{†}❘\alpha \rangle =❘{\Psi }_0\rangle ,$ $U_{\phi }^{†}❘\beta \rangle =\frac{1}{\sqrt{2}}\sum _{x=0}^{N-1}{a}_0\left(x\right)[e^{-2i\phi \left(x\right)}❘0,x\rangle +e^{2i\phi \left(x\right)}❘1,x\rangle ].$

A person of skill in the art will immediately recognize that the subspace spanned by the states |Ψ₀, |α, and |β is almost stable under the action of S_Ψ0, U_φ, and U_φ^†. Only the actions of U_φ on |α and |β may take the state of the system out of this subspace. Continuing to refer to FIGS. 4 and 5, the motivation behind alternating between using U_φ, and U_φ^† during the odd iterations 402 and the even iterations 404 is to keep the state of the system within the subspace spanned by |Ψ₀, |α, and |β.

From the identities above, the expressions of Equations (43) may be written capturing the relevant actions of the odd iteration operator 430, 432 denoted S_Ψ0 U_φ and the even iteration operator 430, 434 denoted S_Ψ0 U_φ^†:

S _{Ψ 0} U _φ|Ψ₀)=2 cos(θ)|Ψ₀−|α,

S _{Ψ 0} U _φ|β)=|Ψ₀,

S _{Ψ 0} U _φ ^†|Ψ₀)=2 cos(θ)|Ψ₀−|β,

S _{Ψ 0} U _φ ^†|α)=|Ψ₀.  (43)

Referring now to FIG. 6, and continuing to refer to FIGS. 4 and 5, state diagram 600 illustrates how, from the expressions of Equation (43), the odd iteration operator 430, 432 denoted S_Ψ0 U_φ may map any (first) state 610 in the space spanned by |Ψ₀ and |β to a (second) state 620 in the space spanned by |Ψ₀ and |α. Conversely, the even iteration operator 430, 434 denoted S_Ψ0 U_φ^† may map any state 620 in the space spanned by |Ψ₀ and |α to a state 610 in the space spanned by |Ψ₀ and |β. Because the algorithm 400, 500 begins with the system initialized in the state |Ψ₀, the state of the system may oscillate between the two subspaces 610, 620 during the odd and even iterations 402, 404.

State After k Iterations:

Using Equations (27) and (43), the state |Ψ_k of the two-register system after k≥0 iterations may be written, in matrix multiplication notation, as shown in Equation (44):

$\begin{array}{cc}❘{\Psi }_k\rangle =\{\begin{array}{cc}{[❘{\Psi }_0\rangle ❘\alpha \rangle ]\left[\begin{array}{cc}2\cos \left(\theta \right)& 1\\ -1& 0\end{array}\right]}^k\left[\begin{array}{c}1\\ 0\end{array}\right],& \mathrm{if}k\mathrm{is}\mathrm{odd},\\ {[❘{\Psi }_0\rangle ❘\beta \rangle ]\left[\begin{array}{cc}2\cos \left(\theta \right)& 1\\ -1& 0\end{array}\right]}^k\left[\begin{array}{c}1\\ 0\end{array}\right],& \mathrm{if}k\mathrm{is}\mathrm{even}\end{array}.& \left(44\right)\end{array}$

To simplify Equation (44), let the matrix M_θ be defined as follows (Equation (45)):

$\begin{array}{cc}{M}_{\theta }\equiv \left[\begin{array}{cc}2\cos \left(\theta \right)& 1\\ -1& 0\end{array}\right].& \left(45\right)\end{array}$

Substituting Equation (45) into Equation (44) yields Equation (46), as follows:

$\begin{array}{cc}❘{\Psi }_k\rangle =\{\begin{array}{cc}{\left[\begin{array}{c}❘{\Psi }_0\rangle \\ ❘\alpha \rangle \end{array}\right]}^T{M}_{\theta }^k\left[\begin{array}{c}1\\ 0\end{array}\right],& \mathrm{if}k\mathrm{is}\mathrm{odd},\\ {\left[\begin{array}{c}❘{\Psi }_0\rangle \\ ❘\beta \rangle \end{array}\right]}^T{M}_{\theta }^k\left[\begin{array}{c}1\\ 0\end{array}\right],& \mathrm{if}k\mathrm{is}\mathrm{even}\end{array}.& \left(46\right)\end{array}$

where the superscript T denotes transposition. M_θ may be diagonalized as

$\begin{array}{cc}{M}_{\theta }=S_{\theta }\left[\begin{array}{cc}{e}^{-i\theta }& 0\\ 0& e^{i\theta }\end{array}\right]S_{\theta }^{-1},& \left(47\right)\end{array}$

and where the matrix S_θ and its inverse S_θ¹ may be given by

$\begin{array}{cc}\begin{array}{cc}{S}_{\theta }=\frac{-i}{2\sin \left(\theta \right)}\left[\begin{array}{cc}{e}^{‐i\theta }& e^{i\theta }\\ -1& -1\end{array}\right],& S_{\theta }^{-1}=\left[\begin{array}{cc}-1& -e^{i\theta }\\ 1& e^{-i\theta }\end{array}\right].\end{array}& \left(48\right)\end{array}$

Now, M_θ^k may be written as follows (Equation (49)):

$\begin{array}{cc}{M}_{\theta }^k=S_{\theta }\left[\begin{array}{cc}{e}^{-i\theta }& 0\\ 0& e^{i\theta }\end{array}\right]S_{\theta }^{-1}=S_{\theta }\left[\begin{array}{cc}{e}^{-ik\theta }& 0\\ 0& e^{ik\theta }\end{array}\right]S_{\theta }^{-1}.& \left(49\right)\end{array}$

From Equation (48) and Equation (49), it follows that

$\begin{array}{cc}\begin{array}{c}{M}_{\theta }^k\left[\begin{array}{c}1\\ 0\end{array}\right]=S_{\theta }\left[\begin{array}{cc}{e}^{-ik\theta }& 0\\ 0& e^{ik\theta }\end{array}\right]\left[\begin{array}{c}-1\\ 1\end{array}\right]\\ =S_{\theta }\left[\begin{array}{c}-e^{-\mathrm{ik}\theta }\\ e^{ik\theta }\end{array}\right]=\frac{1}{\sin \left(\theta \right)}\left[\begin{array}{c}\sin (\left(k+1\right)\theta \\ -s\mathrm{in}\left(k\theta \right)\end{array}\right]\end{array}.& \left(50\right)\end{array}$

Plugging (50) back into (46) leads to the following (Equation (51)):

$\begin{array}{cc}❘{\Psi }_k\rangle =\{\begin{array}{cc}\frac{1}{\sin \left(\theta \right)}[\sin \left(\left(k+1\right)\theta \right)❘{\Psi }_0\rangle -\sin \left(k\theta \right)❘\alpha \rangle ],& \mathrm{if}k\mathrm{is}\mathrm{odd},\\ \frac{1}{\sin \left(\theta \right)}[\sin \left(\left(k+1\right)\theta \right)❘{\Psi }_0\rangle -\sin \left(k\theta \right)❘\beta \rangle ],& \mathrm{if}k\mathrm{is}\mathrm{even}.\end{array}& \left(51\right)\end{array}$

Basis State Amplitudes after k Iterations:

From Equations (16), (33), (34), and (51), the amplitudes ã_k(0,x) of the basis states |0,x after k≥0 iterations may be written as follows (Equation (52)):

$\begin{array}{cc}{\tilde{a}}_k\left(0,x\right)=\langle 0,x❘{\Psi }_k\rangle =\{\begin{array}{cc}\frac{a_0\left(x\right)}{\sqrt{2}\sin \left(\theta \right)}\left[\sin \left(\left(k+1\right)\theta \right)-\sin \left(k\theta \right)e^{i\phi \left(x\right)}\right],& \mathrm{if}k\mathrm{is}\mathrm{odd},\\ \frac{a_0\left(x\right)}{\sqrt{2}\sin \left(\theta \right)}\left[\sin \left(\left(k+1\right)\theta \right)-\sin \left(k\theta \right)e^{-i\phi \left(x\right)}\right],& \mathrm{if}k\mathrm{is}\mathrm{even}.\end{array}& \left(52\right)\end{array}$

Similarly, the amplitudes ã_k(1,x) of the basis states |1,x may be written as follows (Equation (53)):

$\begin{array}{cc}{\tilde{a}}_k\left(1,x\right)=\langle 1,x❘{\Psi }_k\rangle =\{\begin{array}{cc}\frac{a_0\left(x\right)}{\sqrt{2}\sin \left(\theta \right)}\left[\sin \left(\left(k+1\right)\theta \right)-\sin \left(k\theta \right)e^{i\phi \left(x\right)}\right],& \mathrm{if}k\mathrm{is}\mathrm{odd},\\ \frac{a_0\left(x\right)}{\sqrt{2}\sin \left(\theta \right)}\left[\sin \left(\left(k+1\right)\theta \right)-\sin \left(k\theta \right)e^{-i\phi \left(x\right)}\right],& \mathrm{if}k\mathrm{is}\mathrm{even}.\end{array}& \left(53\right)\end{array}$

These expressions may be summarized, for b∈{0,1}, as follows (Equation (54)):

$\begin{array}{cc}{\tilde{a}}_k\left(b,x\right)=\{\begin{array}{cc}\frac{a_0\left(x\right)}{\sqrt{2}\sin \left(\theta \right)}\left[\sin \left(\left(k+1\right)\theta \right)-\sin \left(k\theta \right)e^{i\phi \left(x\right)}\right],& \mathrm{if}k+b\mathrm{is}\mathrm{odd},\\ \frac{a_0\left(x\right)}{\sqrt{2}\sin \left(\theta \right)}\left[\sin \left(\left(k+1\right)\theta \right)-\sin \left(k\theta \right)e^{-i\phi \left(x\right)}\right],& \mathrm{if}k+b\mathrm{is}\mathrm{even}.\end{array}& \left(54\right)\end{array}$

Amplitudes After Ancilla Measurement. Note that the magnitudes of the amplitudes of the states |0,x and |1,x are equal; that is,

|ã₁(0,x)|=|ã_k(1,x)|,  (55)

for all k≥0 and x∈{0, 1, . . . , N−1}. Therefore, a measurement 440 of the ancilla qubit in the first register 410 after K iterations may yield a value of either 0 or 1 with equal probability. Let |ψ_K,b be the normalized state of the second register 420 after performing K iterations, followed by a measurement 440 of the ancilla qubit 410, which may yield a value b∈{0,1}. |ψ_K,b may be written as follows (Equation (56)):

$\begin{array}{cc}❘{\psi }_{K,b}\rangle =\sum _{x=0}^{N-1}{a}_{K,b}\left(x\right)❘x\rangle ,& \left(56\right)\end{array}$

where a_K,b(x) are the normalized amplitudes of the basis states of the second register 420 (after performing K iterations and the measurement 440 of the ancilla 410). a_K,b(x) may simply be given by Equation (57):

a _K,b(x)=√{square root over (2)}/ã_K(b,x).  (57)

Much of the definition hereinbelow holds a) regardless of whether a measurement 440 is performed on the ancilla qubit 410 after the K iterations, and b) regardless of the value yielded by the ancilla measurement 440 (if performed). The primary goal of measuring 440 the ancilla 410 is to make the two registers 410, 420 unentangled from each other.

Basis State Probabilities After K Iterations:

Continuing to refer to FIGS. 4 and 5, let p_K(x) be the probability for a measurement of the second register after K≥0 iterations to yield x. This probability may be written, in terms of the amplitudes in the “Basis State Amplitudes After K Iterations” section hereinabove, as follows (Equation (58)):

p _K(x)=[|ã_K(0,x)|²+|ã_K(1,x)|²]=|a_K,0(x)|²=|a_K,1(x)|².   (58)

This expression shows that the probability p_K(x) depends neither on whether the ancilla was measured, nor on the result of the ancilla measurement (if performed). From Equation (54), p_K(x) may be written as follows (Equation (59)):

$\begin{array}{cc}\begin{array}{c}{p}_K\left(x\right)=\frac{p_0\left(x\right)}{{\sin }^2\left(\theta \right)}{❘\sin \left(\left(K+1\right)\theta \right)-\sin \left(K\theta \right)e^{i\phi \left(x\right)}❘}^2\\ =\frac{p_0\left(x\right)}{{\sin }^2\left(\theta \right)}\left[{\sin }^2\left(K\theta \right)+{\sin }^2\left(\left(K+1\right)\theta \right)-2\sin \left(K\theta \right)\sin \left(\left(K+1\right)\theta \right)\cos \left(\phi \left(x\right)\right)\right].\end{array}& \left(59\right)\end{array}$

A person of skill in the art will immediately recognize that the probability amplification factor p_K(x)/p₀(x) is monotonic in cos(φ(x)) for all K≥0. The following trigonometric identities (Equation (60)) help elucidate the K dependence of this amplification factor:

sin²(C)+sin²(C+D)=sin²(D)+2 sin(C)sin(C+D)cos(D),

2 sin(C)sin(C+D)=cos(D)−cos(2C+D).  (60)

Setting C=Kθ and D=θ, these identities may be used to rewrite Equation (59) as follows (Equation (61)):

p _K(x)=p₀(x){1−λ_K[cos(φ(x))−cos(θ)]},  (61)

where the K-dependent factor λ_K is given by Equation (62):

$\begin{array}{cc}{\lambda }_K\equiv \frac{2\sin \left(K\theta \right)\sin \left(\left(K+1\right)\theta \right)}{{\sin }^2\left(\theta \right)}=\frac{\cos \left(\theta \right)-\cos \left(\left(2K+1\right)\theta \right)}{{\sin }^2\left(\theta \right)}.& \left(62\right)\end{array}$

For notational convenience, the fact that λ_K depends on θ is not explicitly indicated. The result in Equation (61) may be summarized as follows:

(1) Applying K iterations of the non-boolean amplitude amplification algorithm changes the probability of measuring x by a factor that is a linear function of cos(φ(x)). If the second register is initially in an equiprobable state (i.e., if p₀(x)=constant), then the probability p_K(x) after K iterations is itself a linear function of cos(φ(x)).

(2) If cos(φ(x))=cos(θ) for some x, the probability of a measurement of the second register yielding that x is unaffected by the algorithm.

(3) The slope of the linear dependence is −λ_K. If λ_K is positive, the states with cos(φ)<cos(θ) are amplified. Conversely, if λ_K is negative, states with cos(φ)>cos(θ) are amplified. The magnitude of λ_K controls the degree to which the preferential amplification has been performed.

From Equation (62), λ_K is an oscillatory function of K, centered around cos(θ)/sin²(θ) with an amplitude of 1/sin²(θ) and a period of π/θ. Recalling from Equation (30) that

$\begin{array}{cc}\sum _{x=0}^{N-1}{p}_0\left(x\right)\cos \left(\phi \left(x\right)\right)=\cos \left(\theta \right)=\sum _{x=0}^{N-1}{p}_0\left(x\right)\cos \left(\theta \right),& \left(63\right)\end{array}$

one can verify that for any K≥0, the probabilities p_K(x) from Equation (61) add up to 1.

From the definition of λ_K in Equation (62), for all K, λ_K is bounded from above by λ_optimal defined as follows (Equation (64)):

$\begin{array}{cc}{\lambda }_K\le {\lambda }_{\mathrm{optimal}}\equiv \frac{\cos \left(\theta \right)+1}{{\sin }^2\left(\theta \right)}=\frac{1}{1-\cos \left(\theta \right)}.& \left(64\right)\end{array}$

The λ_K=λ_optimal case represents the maximal preferential amplification of lower values of cos(φ) achievable by the algorithm. Let p_optimal(x) be the state probability function corresponding to λ_K=λ_optimal. From Equation (61),

$\begin{array}{cc}{p}_{\mathrm{optimal}}\left(x\right)=p_0\left(x\right)\left[1-\frac{\cos \left(\phi \left(x\right)\right)-\cos \left(\theta \right)}{1-\cos \left(\theta \right)}\right].& \left(65\right)\end{array}$

Note that p_optimal(x)=0 for inputs x with the highest possible value of cos(φ(x)), namely 1. In other words, p_optimal reaches the limit set by the non-negativity of probabilities, in the context of the non-boolean amplitude amplification algorithm 400, 500 described hereinabove.

Number of Iterations to Perform:

In the non-boolean amplitude amplification algorithm 400, 500 described hereinabove, the number of iterations K to perform is left unspecified. Armed with Equation (62), this aspect of the algorithm 400, 500 will now be described in detail. Higher values of λ_K are advantageous to preferentially amplify lower values of cos(φ). Equation (62) illustrates that λ_K may be monotonically increasing for K=0, 1, . . . as long as 0≤(2K+1)θ≤π+θ or, equivalently, for

$\begin{array}{cc}0\le K\le \lfloor \frac{\pi }{2\theta }\rfloor ,& \left(66\right)\end{array}$

where └v┘ denotes the floor of v. A good approach is to stop the algorithm 400, 500 just before the first iteration that, if performed, would cause value of λ_K to decrease (that is, from its value after the previous iteration). This approach leads to the choice K for the number of iterations to perform, given by

$\begin{array}{cc}\tilde{K}=\lfloor \frac{\pi }{2\theta }\rfloor .& \left(67\right)\end{array}$

The corresponding value of λ_K for K={tilde over (K)} is given by Equation (68):

$\begin{array}{cc}{\lambda }_{\tilde{K}}=\frac{1}{{\sin }^2\left(\theta \right)}\left[\cos \left(\theta \right)-\cos \left(\lfloor \frac{\pi }{2\theta }\rfloor 2\theta +\theta \right)\right].& \left(68\right)\end{array}$

The choice {tilde over (K)} in Equation (67) for the number of iterations offers an amplification iff π>2θ>0 or, equivalently, iff 0<cos(θ)<1. At one of the extremes, namely θ=π/2, one has λ_{{tilde over (K)}}=0. The other extreme, namely cos(θ)=1, corresponds to every state x with a non zero amplitude in the initial state |ψ₀ having cos(φ(x))=1. There is no scope for preferential amplification in this case.

Equation (68) illustrates that λ_K exactly equals λ_optimal defined in Equation (64) if π/(2θ) is a half-integer. In terms of θ, this condition may be written as follows (Equation (69)):

$\begin{array}{cc}\theta \in \left\{\frac{\pi }{3},\frac{\pi }{5},\frac{\pi }{7},\dots \left(\mathrm{harmonic}\mathrm{progression}\right)\right\}\Rightarrow {\lambda }_{\tilde{K}}={\lambda }_{\mathrm{optimal}}.& \left(69\right)\end{array}$

For generic values of θ, Equation (68) illustrates that AR satisfies the following (Equation 70)):

$\begin{array}{cc}\frac{2\cos \left(\theta \right)}{{\sin }^2\left(\theta \right)}\le {\lambda }_{\tilde{K}}\le {\lambda }_{\mathrm{optimal}}=\frac{\cos \left(\theta \right)+1}{{\sin }^2\left(\theta \right)}.& \left(70\right)\end{array}$

This equation may be rewritten as follows (Equation (71)):

λ_optimal[1−tan²(θ/2)]≤λ_{{tilde over (K)}}≤λ_optimal,  (71)

using the following identity:

$\begin{array}{cc}\frac{2\cos \left(\theta \right)}{\cos \left(\theta \right)+1}=\frac{2\left[{\cos }^2\left(\theta /2\right)-{\sin }^2\left(\theta /2\right)\right]}{2{\cos }^2\left(\theta /2\right)-1+1}=1-{\tan }^2\left(\theta /2\right).& \left(72\right)\end{array}$

Equation (71) illustrates that for small θ, λ_R is approximately equal to λ_optimal, (θ²) relative error.

Mean and Higher Moments of Cos(φ) after K Iterations:

Let μ_K⁽ⁿ⁾ be the n-th raw moment of cos(φ(x)) for a random value of x sampled by measuring the second register after K iterations.

$\begin{array}{cc}{\mu }_K^{\left(n\right)}\equiv \sum _{x=0}^{N-1}{p}_K\left(x\right){\cos }^n\left(\phi \left(x\right)\right).& \left(73\right)\end{array}$

Under the notation of Equation (73), μ₀⁽¹⁾ is simply cos(θ). From Equation (61),μ_K⁽ⁿ⁾ may be written in terms of the initial moments (K=0) as follows (Equation (74)):

$\begin{array}{cc}\begin{array}{c}{\mu }_K^{\left(n\right)}=\sum _{x=0}^{N-1}{p}_0\left(x\right){\cos }^n\left(\phi \left(x\right)\right)\left[1-{\lambda }_K\left[\cos \left(\phi \left(x\right)\right)-\cos \left(\theta \right)\right]\right]\\ ={\mu }_0^{\left(n\right)}-{\lambda }_K\left[{\mu }_0^{\left(n+1\right)}-{\mu }_0^{\left(n\right)}{\mu }_0^{\left(1\right)}\right]\end{array}.& \left(74\right)\end{array}$

In particular, let μ_K and σ_K² represent the expected value and variance, respectively, of cos(φ(x)) after K iterations.

μ_K≡μ_K⁽¹⁾,

σ_K²≡μ_K⁽²⁾−[μ_K⁽¹⁾]².  (75)

Now, the result in Equation (74) for n=1 may be written as follows (Equation (76)):

μ_K−μ₀=−λ_Kσ₀².  (76)

For λ_K>0, Equation (76) captures the reduction in the expected value of cos(φ(x)) resulting from K iterations of the algorithm 400, 500.

Cumulative Distribution Function of Cos(φ) after K Iterations:

Let F_K^cos(y) be the probability that cos(φ(x))≤y for an x sampled as per the probability distribution p_K. F_K^cos is the cumulative distribution function of cos(φ) for a measurement after K iterations, and may be written as follows (Equations (77) and (78)):

$\begin{array}{cc}{F}_K^{\cos }\left(y\right)=\sum _{x=0}^{N-1}\left[p_K\left(x\right)1_{[0,\infty )}\left[y-\cos \left(\phi \left(x\right)\right)\right]\right],& \left(77\right)\end{array}$ $\begin{array}{cc}1-F_K^{\cos }\left(y\right)=\sum _{x=0}^{N-1}\left[p_K\left(x\right)\left(1-1_{[0,\infty )}\left[y-\cos \left(\phi \left(x\right)\right)\right]\right)\right],& \left(78\right)\end{array}$

where 1_|0,∞) is the Heaviside step function, which equals 0 when its argument is negative and 1 when its argument is non-negative. Using the expression for p_K in Equation (61), these equations may be written as follows (Equations (79) and (80)):

$\begin{array}{cc}{F}_K^{\cos }\left(y\right)=F_0^{\cos }\left(y\right)\left(1+{\lambda }_K{\mu }_0\right)-{\lambda }_K\sum _{x=0}^{N-1}\left\{p_0\left(x\right)\cos \left(\phi \left(x\right)\right)1_{[0,\infty )}\left[y-\cos \left(\phi \left(x\right)\right)\right]\right\},& \left(79\right)\end{array}$ $\begin{array}{cc}1-F_K^{\cos }\left(y\right)=\left(1-F_0^{\cos }\left(y\right)\right)\left(1+{\lambda }_K{\mu }_0\right)-{\lambda }_K\sum _{x=0}^{N-1}\left\{p_0\left(x\right)\cos \left(\phi \left(x\right)\right)\left[1-1_{[0,\infty )}\left[y-\cos \left(\phi \left(x\right)\right)\right]\right]\right\}.& \left(80\right)\end{array}$

Every x that provides a non-zero contribution to the summation in Equation (79) satisfies cos(φ(x))≤y. This fact may be used to write Equation (81), as follows:

λ_K≥0⇒F_K^cos(y)≥F₀^cos(y)(1+λ_K(μ₀−y)).  (81)

Likewise, every x that provides a non-zero contribution to the summation in Equation (80) satisfies cos(φ(x))>y. This fact may be used to write Equation (82), as follows:

λ_K≥0⇒1−F_K^cos(y)≤(1−F₀^cos(y))(1+λ_K(μ₀−y)).  (82)

The inequalities in Equations (81) and (82) may be summarized as follows (Equation (83)):

λ_K≥0⇒F_k^cos(y)≥F₀^cos(y)+λ_K max{F₀^cos(y)(μ₀−y),(1−F₀^cos(y))(y−μ₀)},   (83)

where the max function represents the maximum of its two arguments. This equation provides a lower bound on the probability that a measurement after K iterations yields a state whose cos(φ) value is no higher than y. Derivation of stronger bounds (or even the exact expression) for F_K^cos(y) may be possible if additional information is known about the initial distribution of cos(φ). For y≤μ₀, the first argument of the max function in Equation (83) will be active, and for y≥μ₀, the second argument will be active.

For the λ_K≤0 case, it can similarly be shown that

λ_K≤0⇒F_K^cos(y)≤F₀^cos(y)+λ_K min{F₀^cos(y)(μ₀−y),(1−F₀^cos(y))(y−μ₀)},   (84)

where the min function represents the minimum of its two arguments.

Boolean Oracle Case:

As described in the “Number of Iterations to Perform” section hereinabove, the heuristic choice for the number of iterations for the non-boolean amplitude amplification algorithm of the present invention (namely {tilde over (K)}=└π/(2θ)┘) is different from, but analogous to, the result for known boolean amplitude amplification algorithm of the prior art (namely └π/(4θ_a)┘). Consider the parameter θ in the boolean oracle case, say θ_bool. Let P₀^good be the probability for a measurement on the initial state |ψ₀ to yield a winning state. From Equation (30),

$\begin{array}{cc}\cos \left({\theta }_{\mathrm{bool}}\right)=\left[-1\times P_0^{\mathrm{good}}\right]+\left[1\times \left(1-P_0^{\mathrm{good}}\right)\right]=1-2P_0^{\mathrm{good}},& \left(85\right)\end{array}$ $\begin{array}{cc}\Rightarrow {\sin }^2\left({\theta }_{\mathrm{bool}}/2\right)=P_0^{\mathrm{good}}.& \left(86\right)\end{array}$

Thus, in the boolean oracle case, sin²(θ/2) reduces to the initial probability of “success” (that is, measuring a winning state), which is captured by sin²(θ_a) in known boolean amplitude amplification algorithm(s). The parameter θ used herein reduces to the parameter 2θ_a used in known boolean amplitude amplification algorithm(s), and └π/(2θ)┘ reduces to └π/(4θ_a)┘.

In this way, the results of the “Analysis of the Non-Boolean Amplitude Amplification Algorithm” section hereinabove in general, and the “Number of Iterations to perform” section hereinabove in particular, may be seen as generalizations of the corresponding results in known boolean amplitude amplification algorithm(s) as described hereinabove.

Alternative Formulation of the Algorithm:

Referring again to FIGS. 4 and 5, in the formulation of the non-boolean amplitude amplification algorithm 400, 500 of the present invention, an ancilla qubit (first register) 410 may be included for the purpose of making the quantity Ψ₀|U_φ|Ψ₀=Ψ₀|α real valued. If, in a particular use case, it is guaranteed that ψ₀|U_φ|ψ₀ will be real-valued (or have a negligible imaginary part, for example, and without limitation, by replacing the function φ(x) with φ′(x)=r(x)φ(x), where r: {0, 1, . . . N−1}→{−1,+1} is a random function independent of φ(x), with mean 0 (for x sampled by measuring |ψ₀) even without introducing the ancilla 410, then the algorithm 400, 500 described above may be used without the ancilla 410: that is, may alternate between applying operator 430, 432 denoted S_ψ0 U_φ during the odd iterations and applying operator 430, 434 denoted S_ψ0 U_φ^† during the even iterations. In other words, the properties and structure of two-register system are not exploited in the algorithm description, beyond making Ψ₀|U_φ|Ψ₀ real-valued.

However, Equations (33) and (34) illustrate that the states |α=U_φ|Ψ₀ and |β=U_φ^†|Ψ₀ are related by

|β=[X⊗I]|α,|α=[X⊗I]|β,  (87)

where X is the bit-flip or the Pauli-X operator 340. This relation may be exploited to avoid having two separate cases—k being odd and even—in the final expression for |Ψ_k in Equation (51). The expression for both cases may be made the same by acting the Pauli-X gate 340 on the ancilla 410, once at the end, if the total number of iterations is even.

More interestingly, the relationship between |α and |β may be used to avoid having two different operations in the first place, for the odd and even iterations. This leads to the following alternative formulation of the non-boolean amplitude amplification algorithm: During each iteration, odd or even, act the same operator Q_iter defined by

Q _iter ≡S _{Ψ 0} U _φ[X⊗I].  (88)

This alternative formulation is depicted as a circuit 700 in FIG. 7 and as a pseudocode 800 in FIG. 8.

From Equations (43) and (87), the action of Q_iter on and |α may be derived as follows (Equations (89) and (90)):

$\begin{array}{cc}{Q}_{iter}❘{\Psi }_0\rangle =S_{{\Psi }_0}{U}_{\phi }❘{\Psi }_0\rangle =2\cos \left(\theta \right)❘{\Psi }_0\rangle -❘\alpha \rangle ,& \left(89\right)\end{array}$ $\begin{array}{cc}{Q}_{iter}❘\alpha \rangle =S_{{\Psi }_0}{U}_{\phi }❘\beta \rangle =❘{\Psi }_0\rangle .& \left(90\right)\end{array}$

Let |Ψ_k^alt be the state of the two-register system after k iterations under this alternative formulation 700, 800 (before any measurement of the ancilla 410), as follows (Equation (91)):

|Ψ_k^alt)≡Q_iter^k|Ψ₀.  (91)

Using similar manipulations as those leading to Equation (51) hereinabove, |Ψ_k^alt) may be expressed, for all k≥0, as follows (Equation (92)):

$\begin{array}{cc}❘{\Psi }_k^{\mathrm{alt}}\rangle =\frac{1}{\sin \left(\theta \right)}[\sin \left(\left(k+1\right)\theta \right)❘{\Psi }_0\rangle -\sin \left(k\theta \right)❘\alpha \rangle ].& \left(92\right)\end{array}$

Note that this expression for |Ψ_k^alt) is almost identical to the expression for |Ψ_k in Equation (51), but without two separate cases for the odd and even values of k. Much of the analysis of the original formulation of the non-boolean amplitude amplification algorithm 400, 500 described hereinabove holds for the alternative formulation 700, 800 described hereinabove as well, including the expressions for the state probabilities p_K(x), mean μ_K, raw moments λ_K⁽ⁿ⁾, and the cumulative distribution function F_K^cos.

In addition to simplifying the amplification algorithm (by using the same operation for every iteration), the Q_iter operator used in algorithms 700, 800 allows for a clearer presentation of the quantum mean estimation algorithm, as described hereinbelow.

Quantum Mean Estimation Algorithm:

A quantum mean estimation algorithm according to certain embodiments of the present invention may operate to estimate the expected value of e^iφ(x) for x sampled by measuring a given superposition state |ψ₀. Let E_ψ0[e^iφ] denote this expected value, as follows (Equation (93)):

$\begin{array}{cc}{E}_{{\psi }_0}\left[e^{i\phi }\right]\equiv \langle {\psi }_0❘U_{\phi }❘{\psi }_0\rangle =\sum _{x=0}^{N-1}{❘a_0\left(x\right)❘}^2{e}^{i\phi \left(x\right)}.& \left(93\right)\end{array}$

This equation may be written as follows (Equation (94)):

E _{ψ 0} [e^iφ]=Re[E_ψ0[e^iφ]]+i Im[E_ψ0[e^iφ]],  (94)

where the real and imaginary parts are given by Equations (95) and (96), respectively, as follows:

$\begin{array}{cc}\mathrm{Re}\left[E_{{\psi }_0}\left[e^{i\phi }\right]\right]=\sum _{x=0}^{N-1}{❘a_0\left(x\right)❘}^2\cos \left(\phi \left(x\right)\right),& \left(95\right)\end{array}$ $\begin{array}{cc}\mathrm{Im}\left[E_{{\psi }_0}\left[e^{i\phi }\right]\right]=\sum _{x=0}^{N-1}{❘a_0\left(x\right)❘}^2\sin \left(\phi \left(x\right)\right).& \left(96\right)\end{array}$

The mean estimation may therefore be performed in two parts: one for estimating the mean of cos(φ), and the other for estimating the mean of sin(φ). Note that the expectation of cos(φ) under the state |ψ₀ is precisely cos(θ) defined in Equation (30).

Estimating the Mean of cos(φ):

The connection shown in the “Boolean Oracle Case” section hereinabove between the parameter θ and the parameter θ_a serves as the intuition behind the quantum mean estimation algorithm of the present invention. In the known boolean amplitude estimation algorithm described hereinabove, the parameter θ_a is estimated using QPE (Note: The estimation of θ_Q is only (needed to be) performed up to a two-fold ambiguity of {θ_a,π−θ_a}). The estimate for θ_a is then turned into an estimate for the initial winning probability. The non-boolean quantum mean estimation algorithm of the present invention may operate to estimate the parameter θ defined in Equation (30) using QPE (Note: The estimation of θ is only (needed to be) performed up to a two fold ambiguity of {θ,2π−θ}). The estimate for θ may then be translated into an estimate for the initial expected value of cos(φ(x)), namely cos(θ).

Toward the end of actualizing the abovementioned intuition into a working algorithm, a key observation is that |Ψ may be written as follows (Equation (97)):

$\begin{array}{cc}❘{\Psi }_0\rangle =\frac{❘{\eta }_{+}\rangle -❘{\eta }_{-}\rangle }{\sqrt{2}},& \left(97\right)\end{array}$

where |η₊ and |η₋ are given by Equation (98):

$\begin{array}{cc}❘{\eta }_{\pm }\rangle =\frac{e^{\pm i\theta }❘{\Psi }_0\rangle -❘\alpha \rangle }{i\sqrt{2}\sin \left(\theta \right)}.& \left(98\right)\end{array}$

Note: The forms of the eigenstates |η₊ and |η₋ in Equation (98) may be determined from the form of the matrix S_θ in Equation (48). These forms may also be determined from Equation (92) by rewriting the sin functions in terms of complex exponential functions.

The expressions for |η₊ and |η₋ in Equation (98) may be used to verify Equation (97). Crucially, |η₊ and |η₋ are unit normalized eigenstates of the unitary operator Q_iter, with eigenvalues e^iθ and e^−iθ, respectively.

Q _iter|η_±=e^±iθ|η_±,  (99)

η_±|η_±=1.  (100)

The properties of |η₊ and |η₋ in Equation (99) and Equation (100) may be verified using Equation (89), Equation (90), and Equation (98). The observations in Equation (97) and Equation (99) lead to the following algorithm for estimating cos(θ):

(1) Perform the QPE algorithm with a) the two-register operator Q_iter serving the role of the unitary operator under consideration, and b) the superposition state |Ψ₀) in place of the eigenstate required by the QPE algorithm as input. Let the output of this step, appropriately scaled to be an estimate of the phase angle in the range [0,2π), be {circumflex over (ω)}.

(2) Return cos({circumflex over (ω)}) as the estimate for cos(θ) (i.e., the real part of E_ψ0[e^iφ]). Note: If the circuit implementation of Q_iter is wrong by an overall (state independent) phase ϕ_err, then the estimate for cos(θ) is cos({circumflex over (ω)}−ϕ_err). This is important, for example, if the operation [2|0,00,0|−I] is only implemented up to a factor of −1 (i.e., with ϕ_err=π). Note that the final state probabilities under the non-boolean amplitude amplification algorithm are unaffected by such an overall phase error.

Proof of correctness of the algorithm: |Ψ₀ is a superposition of the eigenstates |η₊ and |η₋ of the unitary operator Q_iter. This implies that {circumflex over (ω)} may either be an estimate for the phase angle of |η₊, namely θ, or an estimate for the phase angle of |η₋, namely 2π−θ. (Note: If θ=0, the phase angle being estimated is 0 for both |η₊ and |η₋. For θ≠0, {circumflex over (ω)} will be an estimate for either θ or 2π−θ with equal probability. Since θ lies in [0,π] and 2π−θ lies in [π,2π], the output {circumflex over (ω)} can be converted into an estimate for θ alone (although doing so is not necessary)). Since, cos(2π−θ)=cos(θ), it follows that cos({circumflex over (ω)}) is an estimate for cos(θ).

Estimating the Mean of e^iφ:

The algorithm for estimating the expected value of cos(φ) in the section hereinabove may be re-purposed to estimate the expected value of sin(φ) by using the fact that

sin(φ)=cos(φ−π/2).  (101)

In other words, the imaginary part of E_ψ0[e^iφ] is the real part of E_ψ0[e^i(φ−π/2)]. By using the oracle U_φ−π/2 (for the function φ−π/2), instead of U_φ, in the mean estimation algorithm of the preceding section, the imaginary part of E_ψ0[e^iφ] may also be estimated. This completes the estimation of E_ψ0[e^iφ].

For concreteness, U_φ−π2 may be explicitly written as

U _φ−π/2 =e ^−iπ/2|00|⊗U_φ+e^iπ/2|11|⊗U_φ^†.  (102)

Referring now to FIG. 9, exemplary circuit 900 illustrates implementation of U_φ−π/2 using the oracle 932 denoted U_φ, the bit-flip operator 940 denoted X, and the phase-shift operator R_ϕ (shown as elements 950, 952).

Note that the algorithm does not use the knowledge that {|0, . . . , |N−1} is an eigenbasis of U_φ. Therefore, this algorithm may be used to estimate ψ₀|U|ψ₀ for any unitary operator U.

Quantum Speedup:

The speedup offered by the quantum mean estimation algorithm 900 over classical methods will now be discussed in detail, in the context of estimating the mean of cos(φ) alone. This discussion may be extended in a straightforward way to estimation of E_ψ0[e^iφ].

Classical Approaches to Estimating the Mean: For an arbitrary function φ and a known sampling distribution p₀(x) for the inputs x, one classical approach to finding the mean of cos(φ(x)) is to sequentially query the value of φ(x) for all the inputs and use the query results to compute the mean. Let the permutation (x₀, x₁, . . . , x_N-1) of the inputs (0, 1, . . . , N−1) be the order in which the inputs are queried. The range of allowed values of cos(θ), based only on results for the first q inputs, is given by the following (Equation (103)):

$\begin{array}{cc}\left[\sum _{j=0}^{q-1}{p}_0\left(x\right)\cos \left(\phi \left(x_j\right)\right)-\sum _{j=q}^{N-1}{p}_0\left(x\right)\right]\le \cos \left(\theta \right)\le \left[\sum _{j=0}^{q-1}{p}_0\left(x\right)\cos \left(\phi \left(x_j\right)\right)+\sum _{j=q}^{N-1}{p}_0\left(x\right)\right].& \left(103\right)\end{array}$

These bounds are derived by setting the values of cos(φ) for all the unqueried inputs to their highest and lowest possible values, namely +1 and −1. The range of allowed values of cos(θ) shrinks as more and more inputs are queried. In particular, if p₀(x) is equal for all the inputs x, the width of the allowed range (based on q queries) is given by 2(N−q)/N. This strategy may take (N) queries before the width of the allowed range reduces to even, say, 1. Thus, this strategy will not be feasible for large values of N.

An alternative classical approach is to probabilistically estimate the expected value as follows:

(1) Independently sample q random inputs (x₁, . . . x_q) as per the distribution p₀.

(2) Return the sample mean of cos(φ) over the random inputs as an estimate for cos(θ).

Under this approach, the standard deviation of the estimate scales as ˜σ₀/√{square root over (q)}, where σ₀ is the standard deviation of cos(φ) under the distribution p₀.

Precision Vs Number of Queries for the Quantum Algorithm: Note that one call to the operator Q_iter corresponds to (1) calls to A₀ and U_φ (and their inverses). Let q be the number of times the (controlled) Q_iter operation is performed during the QPE subroutine. As q increases, the uncertainty on the estimate for the phase-angle θ (up to a two-fold ambiguity) falls at the rate of (1/q). Consequently, the uncertainty on cos(θ) also falls at the rate of (1/q). This represents a quadratic speedup over the classical, probabilistic approach, under which the error falls as (1/√{square root over (q)}).

Note that the variance of the estimate for cos(θ) is independent of a) the size of input space N, and b) the variance σ₀² of cos(φ(x)) under the distribution p₀(x). It only depends on the true value of cos(θ) and the number of queries q performed during the QPE subroutine.

Demonstrating the Algorithms Using a Toy Example

Referring again to FIGS. 4, 5, and 9, the non-boolean amplitude amplification algorithm 400, 500 and the mean estimation algorithm 900 will both be demonstrated using a toy example. Let the input to the oracle U_φ (i.e., the second register 420), contain 8 qubits. This leads to 2⁸=256 basis input states, namely |0, . . . , |255. Let the toy function φ(x) be

$\begin{array}{cc}\begin{array}{cc}\phi \left(x\right)=\frac{x}{255}\frac{\pi }{4},& \mathrm{for}x=0,1,\dots ,255.\end{array}& \left(104\right)\end{array}$

The largest phase-shift applied by the corresponding oracle 432 denoted U_φ on any basis state is π/4, for the state |255. Since, cos(φ(x)) is monotonically decreasing in x, the goal of the amplitude amplification algorithm 400, 500 is to amplify the probabilities of higher values of x.

Let the initial state, from which the amplification is performed, be the uniform superposition state |s.

$\begin{array}{cc}❘{\psi }_0\rangle =❘s\rangle =\frac{1}{\sqrt{256}}\sum _{x=0}^{255}❘x\rangle .& \left(105\right)\end{array}$

Such simple forms for the oracle function and the initial state allow for a good demonstration of the algorithms of the present invention.

For this toy example, from Equation (30), cos(θ) and θ are given by

$\begin{array}{cc}\begin{array}{c}\cos \left(\theta \right)=\frac{1}{256}\sum _{x=0}^{255}\cos \left(\frac{x}{255}\frac{\pi }{4}\right)\approx 0.9001,\\ \theta \approx 0.457.\end{array}& \left(106\right)\end{array}$

Referring now to FIG. 10, graph 1000 shows the value 1010 of λ_K, from Equation (62), for the first few of a range 1020 of values of K (more specifically, for K=0, 1, . . . , 14). The dots 1030 (only a subset labeled, for clarity) correspond to the different integer values 1020 of K. The black solid curve 1040 depicts the sinusoidal dependence of the value 1010 of λ_K on the range 1040 of K. The dotted lines 1050 indicate that the value 1010 of λ_K oscillates around a reference 1060 of cos(θ)/sin²(θ) with an amplitude 1070 of 1/sin²(θ) and a period 1080 of π/θ (in K values). The heuristic choice for the total number of iterations {tilde over (K)}=[π/(2θ)] is 3 for this example.

Amplitude Amplification:

Continuing with the toy example, and referring now to chart 1100 of FIG. 11, the quantum circuit for the non-boolean amplitude amplification algorithm 400, 500 was implemented in a quantum computing framework for three different values of the total number of iterations K, namely K=1, 2, and 3. The solid histograms 1130, 1140, 1150 show the observed measurement frequencies 1110 of the different values 1120 of x∈{0, 1 . . . , 255} after performing the non-boolean amplitude amplification algorithm 400, 500. The solid histograms 1130, 1140, 1150 correspond to the total number of iterations K being 1, 2, and 3, respectively. In each case, the observed frequencies are based on simulating (and measuring) the circuit for the algorithm 10⁶ times (i.e., 10⁶ shots). The dashed curves 1132, 1142, 1152, in good agreement (almost coincident) with their corresponding solid histograms 1130, 1140, 1150, respectively, show the predictions p_K(x) (for the measurement frequencies) computed using Equation (61). While p_K(x) may be technically defined only for the integer values of x, in FIG. 11 the dashed curves 1132, 1142, 1152 are interpolated for non-integer values of x using Equation (61).

As can be seen from FIG. 11, in each case, the algorithm 400, 500 preferentially amplifies lower values of cos(φ(x)) or, equivalently, higher values of x. This result is expected from the fact that λ_K>0 for all three values of K. Furthermore, as K increases from 0 to {tilde over (K)}=3, the preferential amplification grows stronger. Note that the probabilities of the x-s for which cos (Φ(x))≈cos(θ) are left approximately unchanged by the algorithm, as indicated by the dotted crosshair 1160.

Mean Estimation:

Continuing with the toy example, and referring now to graphs 1200, 1250 of FIG. 12, results of demonstrating the quantum circuit for the non-boolean mean estimation algorithm 900 will now be described in detail. Only the estimation of cos(θ) (i.e., the real part of E_ψ0[e^iφ]) is demonstrated herein. The imaginary part may also be estimated using the same technique, as described in “Estimating the Mean of e^iφ” hereinabove.

Let M be the number of qubits used in the QPE subroutine of the mean estimation algorithm 900, to contain the phase information. This corresponds to performing the (controlled) Q_iter operation 2^M−1 times during the QPE subroutine. Note that the estimated phase {circumflex over (ω)} may only take the following discrete values:

$\begin{array}{cc}\hat{\omega }\in \left\{\frac{2\pi j}{2^M}❘j\in \left\{0,\dots ,2^M-1\right\}\right\}.& \left(107\right)\end{array}$

In this way, the value of M controls the precision of the estimated phase and, by extension, the precision of the estimate for cos(θ); the higher the value of M, the higher the precision.

For the demonstration, two different quantum circuits were implemented, again using a quantum computing framework, for the mean estimation algorithm 900; one with M=4 and the other with M=8. Each circuit was simulated (and measured) 10⁶ times, to get a sample of {circumflex over (ω)} values all in the range [0,2π).

The observed frequencies 1210 (scaled by 1/bin-width) of the different values 1220 of {circumflex over (ω)} are shown as histograms on a linear scale in top graph 1200, and on a logarithmic scale in the bottom graph 1250. Here the bin-width of the histograms is given by 2π/2^M, which is the difference between neighboring allowed values of {circumflex over (ω)}.

The dashed histogram 1230 and solid histogram 1240 correspond to the circuits with 4 and 8 phase measurement qubits, respectively. The exact values of θ and 2π−θ for the toy example are indicated with vertical dotted lines 1232, 1242, respectively. In both cases (M=4 and M=8), the observed frequencies peak near the exact values of θ and 2π−θ, demonstrating that {circumflex over (ω)} is a good estimate for them, up to a two-fold ambiguity. (Note: The upward trends near the left (−1) and right (+1) edges of the plots in are artifacts caused by the Jacobian determinant for the map from {circumflex over (ω)} to cos({circumflex over (ω)})). Furthermore, as expected, using more qubits for estimating the phase leads to a more precise estimate.

Ancilla-Free Versions of the Algorithms:

Both algorithms introduced this in this paper so far, namely (1) the amplitude algorithm of FIGS. 4 and 5 and its alternative formulation in FIG. 8, and (2) the mean estimation algorithm of FIG. 9, use an ancilla qubit to make the quantity Ψ₀|U_φ|Ψ₀ real valued. This is important for achieving the respective goals of the algorithms. However, the same algorithms may be performed without the ancilla, albeit to achieve different goals, which may be relevant in some use cases.

Ancilla-Free Non-Boolean Amplitude Amplification:

The ancilla-free version of the amplitude amplification algorithm is almost identical to the algorithm introduced hereinabove at FIGS. 4 and 5. The only difference is that, in the ancilla-free version, the single-register operators S_ψ0, U_φ, and U_φ^† are used in place of the two-register operators S_Ψ0, U_φ, and U_φ^†, respectively. For concreteness, the ancilla-free algorithm may proceed as follows:

(1) Initialize a system in the state |ψ₀.

(2) Act the operation S_ψ0U_φ during the odd iterations and S_ψ0U_φ^† during the even iterations.

Analogous to the two-register states |α and |β in Equations (31) and (32), let the single-register states |α′ and |β′ be defined as

$\begin{array}{cc}\begin{array}{cc}2❘{\alpha }^{\prime }\rangle \equiv U_{\phi }❘{\psi }_0\rangle & =\sum _{x=0}^{N-1}{e}^{i\phi \left(x\right)}{a}_0\left(x\right)❘x\rangle ,\\ ❘{\beta }^{\prime }\rangle \equiv U_{\phi }^{†}❘{\psi }_0\rangle & =\sum _{x=0}^{N-1}{e}^{-i\phi \left(x\right)}{a}_0\left(x\right)❘x\rangle .\end{array}& \left(108\right)\end{array}$

Analogous to θ defined in Equation (30), let θ′∈[0,π/2] and δ∈[0,2π) be implicitly defined by

$\begin{array}{cc}\cos \left({\theta }^{\prime }\right)e^{i\delta }\equiv \langle {\psi }_0❘{\alpha }^{\prime }\rangle =\sum _{x=0}^{N-1}{❘a_0\left(x\right)❘}^2{e}^{i\phi \left(x\right)}.& \left(109\right)\end{array}$

Variables cos(θ′) and δ are the magnitude and phase, respectively, of the initial (i.e., x sampled from |ψ₀) expected value of e^iφ. An important difference between θ′ and θ is that cos(θ′) is restricted to be non-negative, unlike cos(θ), which may be positive, negative, or zero.

Note that cos(θ′) may be written as

$\begin{array}{cc}\cos \left({\theta }^{\prime }\right)=\sum _{x=0}^{N-1}{❘a_0\left(x\right)❘}^2{e}^{i\mathrm{\phi \prime }\left(x\right)},& \left(110\right)\end{array}$

where φ′(x) is given by

φ′(x)≡φ(x)−δ.  (111)

Acting the oracle U_φ for the function φ may be thought of as acting the oracle U_φ′ for the function φ′, followed by performing a global, state independent phase-shift of δ. Furthermore, from Equation (120), it can be seen that ψ₀|U_φ′|ψ₀ is real-valued. This observation may be used to re-purpose the analysis in the “Analysis of the Non-Boolean Amplitude Amplification Algorithm” section hereinabove for the ancilla-free version; the corresponding results are presented here without explicit proofs.

Let |ψ_k be the state of the system of the after k≥0 iterations of the ancilla-free algorithm. Analogous to Equation (51), |ψ′_k may be written as

$\begin{array}{cc}❘{\psi }_k^{\prime }\rangle =\{\begin{array}{cc}\frac{e^{i\delta }}{\sin \left(\theta \text{'}\right)}[\sin \left(\left(k+1\right){\theta }^{\prime }\right)❘{\psi }_0\rangle -\sin \left(k{\theta }^{\prime }\right)e^{-i\delta }❘{\alpha }^{\prime }\rangle ],& \mathrm{if}k\mathrm{is}\mathrm{odd},\\ \frac{1}{\sin \left({\theta }^{\prime }\right)}[\sin \left(\left(k+1\right){\theta }^{\prime }\right)❘{\psi }_0\rangle -\sin \left(k{\theta }^{\prime }\right)e^{i\delta }❘{\beta }^{\prime }\rangle ],& \mathrm{if}k\mathrm{is}\mathrm{even}.\end{array}& \left(112\right)\end{array}$

Let p′_K(x) be the probability of measuring the system in state x after K iterations. Analogous to Equation (61), p′_K(x) can be written as follows (Equation (113)):

p′ _K(x)=p₀(x){1−λ′_K[cos(φ(x)−δ)−cos(θ′)]},  (113)

where the λ′_K, the ancilla-free analogue of λ_K, is given by

$\begin{array}{cc}{\lambda }_K^{\prime }=\frac{2\sin \left(K{\theta }^{\prime }\right)\sin \left(\left(K+1\right){\theta }^{\prime }\right)}{{\sin }^2\left({\theta }^{\prime }\right)}=\frac{\cos \left({\theta }^{\prime }\right)-\cos \left(\left(2K+1\right){\theta }^{\prime }\right)}{{\sin }^2\left({\theta }^{\prime }\right)}.& \left(114\right)\end{array}$

In this case, the probability amplification factor p′_K/p₀ is linear in cos(φ−δ).

Ancilla-Free Mean Estimation:

The ancilla-free mean estimation algorithm described in this subsection may estimate the magnitude of ψ₀|U|ψ₀ for a given unitary operator U. Hereinbelow the algorithm is presented in terms of the oracle U_φ, and the goal of the algorithm is to estimate cos(θ′) from Equation (109) (i.e., the magnitude of E_ψ0[e^iφ]≡ψ₀|U_φ|ψ₀).

Let the unitary operator Q_evenodd be defined as follows (Equation 115)):

Q _evenodd ≡S _{ψ 0} U _φ S _{ψ 0} U _φ ^†.  (115)

Its action corresponds to performing the (ancilla-free) even-iteration operation once, followed by the odd-iteration operation. Analogous to Equations (97) and (98), the state |ψ₀ can be written as

$\begin{array}{cc}❘{\psi }_0\rangle =\frac{❘{\eta }_{+}^{\prime }\rangle -|{\eta }_{-}^{\prime }\rangle }{\sqrt{2}},& \left(116\right)\end{array}$

where |η′₊ and |η″₋ are given by

$\begin{array}{cc}❘{\eta }_{\pm }^{\prime }\rangle =\frac{e^{\pm i\theta \prime }❘{\psi }_0\rangle -e^{-i\delta }❘{\alpha }^{\prime }\rangle }{i\sqrt{2}\sin \left({\theta }^{\prime }\right)}.& \left(117\right)\end{array}$

|η′₊ and |η′₋ are unit-normalized eigenstates of Q_evenodd with eigenvalues e^2iθ, and e^−2iθ, respectively.

$\begin{array}{cc}{Q}_{e\mathrm{ven}\mathrm{odd}}❘{\eta }_{\pm }^{\prime }\rangle =e^{\pm 2i\theta \prime }❘{\eta }_{\pm }^{\prime }\rangle ,& \left(118\right)\end{array}$ $\begin{array}{cc}\langle {\eta }_{\pm }^{\prime }❘{\eta }_{\pm }^{\prime }\rangle =1.& \left(119\right)\end{array}$

The observations in Equations (118) and (119) lead to the following algorithm for estimating cos(θ′):

(1) Perform the QPE algorithm with Q_evenodd serving the role of the unitary operator under consideration, and the superposition state |ψ₀ in place of the eigenstate required by the QPE algorithm as input. Let the output of this step, appropriately scaled to be an estimate of the phase angle in the range [0,2π), be {circumflex over (ω)}.

(2) Return |cos({circumflex over (ω)}/2)| as the estimate for cos(θ′).

In this version of the algorithm, {circumflex over (ω)} will be an estimate for either 2θ or 2π−2θ. (Note: If θ′=0, the phase angle being estimated is 0 for both |η′₊ and |η′₋.) So, {circumflex over (ω)}/2 will be an estimate for either θ′ or π−θ′. Since, a) cos(π−θ′)=−cos(θ′), and b) cos(θ′) is a non-negative number, it follows that |cos({circumflex over (ω)}/2)| is an estimate for cos(θ′).

Potential advantageous applications of the algorithms of the present invention may include the following:

Approximate Optimization:

A straightforward application of the non-boolean amplitude amplification algorithm of the present invention is in the optimization of objective functions defined over a discrete input space. The objective function to be optimized needs to be mapped onto the function φ of the oracle U_φ, with the basis states of the oracle corresponding to the different discrete inputs of the objective function. After performing an appropriate number of iterations of the algorithm, measuring the state of the system will yield “good” states with amplified probabilities. Multiple repetitions of the algorithm (multiple shots) may be performed to probabilistically improve the quality of the optimization.

Note that the technique is not guaranteed to yield the true optimal input, and the performance of the technique will depend crucially on factors like a) the map from the objective function to the oracle function φ, b) the number of iterations K, c) the initial superposition |ψ₀, and, in particular, d) the initial distribution of φ under the superposition |ψ₀. This approach joins known quantum optimization techniques.

Simulating Probability Distributions:

The amplitude amplification algorithm may be useful for simulating certain probability distributions. By choosing the initial state |ψ₀, oracle U_φ, and the number of iterations K, one can control the final sampling probabilities p_K(x) of the basis states; the exact expression for p_K(x) in terms of these factors is given in Equation (61).

Estimating the Overlap Between Two States:

Let |ψ and |ϕ be two different states produced by acting the unitary operators A and B, respectively, on the state |0.

|ψ=A|0,|ϕ=B|0.  (120)

Estimating the overlap |ψ|ϕ| between the two states is an important task with several applications, including in Quantum Machine Learning (QML). Several algorithms, including the Swap test, can be used for estimating this overlap. For the Swap test, the uncertainty in the estimated value of |ψ|ϕ| falls as (1/√{square root over (q)}) in the number of queries q to the unitaries A and B (used to the create the states |ψand |ϕ).

On the other hand, the mean estimation algorithm of the present invention also be used to estimate ψ|ϕby noting that

ψ|ϕ=0|A^†B|0.  (121)

So, by setting

U≡A ^† B,  (122)

|ψ₀≡|0,  (123)

values ψ|ϕ may be estimated as ψ₀|U|ψ₀ using the mean estimation algorithm of FIG. 9. If one is only interested in the magnitude of ψ|ϕ, the ancilla-free version of the mean estimation algorithm described hereinabove will also suffice. Since, for the mean estimation algorithm, the uncertainty of the estimate falls as (1/q) in the number of queries q to the unitaries A and B (or their inverses), this approach offers a quadratic speedup over the Swap test. Furthermore, the (1/q) scaling of the error achieved by this approach matches the performance of the known optimal quantum algorithm for the overlap-estimation task.

Meta-Oracles to Evaluate the Superposition |ψ and Unitary U:

Recall from Equations (97) and (98) that

$\begin{array}{cc}{Q}_{iter}❘{\Psi }_0\rangle =\frac{e^{\dot{t}\theta }❘{\eta }_{+}\rangle -e^{-i\theta }❘{\eta }_{-}\rangle }{\sqrt{2}},& \left(124\right)\end{array}$

where cos(θ) is the real part of ψ₀|U_φ|ψ₀. Note that the parameter θ depends on the superposition |ψ₀ and the unitary U_φ. The action of Q_iter on |Ψ₀) is to apply a phase-shift of θ on the projection along |η₊ and a phase-shift of −θ on the projection along |η₋. This property may be used to create a meta-oracle which evaluates the superposition |Ψ₀ and/or the unitary U_φ (or a generic unitary U) based on the corresponding value of θ. More specifically, if the circuit A₀ for producing |Ψ₀ and/or the circuit for U are additionally parameterized using “control” quantum registers (provided as inputs to the circuits), then a meta-oracle may be created using Equation (24) to evaluate the states of the control registers. Such meta-oracles may be used with quantum optimization algorithms, including the non-boolean amplitude amplification algorithm of the present invention, to find “good” values (or states) of the control registers.

Variational quantum circuits (i.e., quantum circuits parameterized by (classical) free parameters) have several applications, including in QML. Likewise, quantum circuits parameterized by quantum registers may also have applications (e.g., in QML and quantum statistical inference).

Some of the illustrative aspects of the present invention may be advantageous in solving the problems herein described and other problems not discussed which are discoverable by a skilled artisan.

While the above description contains much specificity, these should not be construed as limitations on the scope of any embodiment, but as exemplifications of the presented embodiments thereof. Many other ramifications and variations are possible within the teachings of the various embodiments. While the invention has been described with reference to exemplary embodiments, it will be understood by those skilled in the art that various changes may be made and equivalents may be substituted for elements thereof without departing from the scope of the invention. In addition, many modifications may be made to adapt a particular situation or material to the teachings of the invention without departing from the essential scope thereof. Therefore, it is intended that the invention not be limited to the particular embodiment disclosed as the best or only mode contemplated for carrying out this invention, but that the invention will include all embodiments falling within the scope of the appended claims. Also, in the drawings and the description, there have been disclosed exemplary embodiments of the invention and, although specific terms may have been employed, they are unless otherwise stated used in a generic and descriptive sense only and not for purposes of limitation, the scope of the invention therefore not being so limited. Moreover, the use of the terms first, second, etc. do not denote any order or importance, but rather the terms first, second, etc. are used to distinguish one element from another. Furthermore, the use of the terms a, an, etc. do not denote a limitation of quantity, but rather denote the presence of at least one of the referenced item.

Thus, the scope of the invention should be determined by the appended claims and their legal equivalents, and not by the examples given.

Claims

1. A method of performing quantum calculation on an oracle Uφ for a non-boolean function φ, comprising: initializing an ancilla qubit in a |+ state and an input qubit in a |ψ0 state of a plurality of eigenstates |x to define a two-register state |Ψ0≡|+,ψ0;for each of a plurality K of iterations receiving, using the input qubit, a respective one of the plurality of eigenstates |x defining an input basis state,for odd iterations, acting on the input basis state using a selective phase-flip unitary operator circuit SΨ0 and a controlled unitary operator circuit Uφ, andfor even iterations, acting on the input basis state using the selective phase-flip unitary operator circuit SΨ0 and a controlled inverse unitary operator circuit Uφ†.

2. The method according to claim 1, further comprising measuring, after the plurality K of iterations, the ancilla qubit in a 0/1 basis.

3. The method according to claim 1, wherein the selective phase-flip unitary operator circuit SΨ0 further comprises: a Hadamard transform H,a unitary operator A0,a unitary operator I, andan inverse unitary operator A0†;wherein the selective phase-flip unitary operator SΨ0=[H⊗A0][2|0,00,0|−I][H⊗A0†].

4. The method according to claim 1, wherein the acting on the input basis state using the selective phase-flip unitary operator circuit SΨ0 and the controlled unitary operator circuit Uφ further comprises, for the ancilla qubit in a state |0, acting on the input qubit as Uφ|0,x=e+iφ(x)|0,x, Uφ|1,x=e−iφ(x)|1,x.

5. The method according to claim 1, wherein the acting on the input basis state using the selective phase-flip unitary operator circuit SΨ0 and the controlled inverse unitary operator circuit Uφ† further comprises, for the ancilla qubit in a state |1, acting on the input qubit as Uφ†|0,x=e+iφ(x)|0,x, Uφ†|1,x=e−iφ(x)|1,x.

6. The method according to claim 1, further comprising: receiving, using the input qubit, the |ψ0 state defining an input random state;acting on the input random state using a controlled estimation unitary operator circuit Uφ−π/2.

7. The method according to claim 6, wherein the controlled estimation unitary operator circuit Uφ−π/2 further comprises: the controlled unitary operator circuit Uφ,at least one bit-flip operator X, andat least one phase-shift operator Rϕ;wherein the controlled estimation unitary operator circuit Uφ−π/2=e−iπ/2|00|⊗Uφ+eiπ/2|11|⊗Uφ†.

8. A quantum computing device for performing quantum calculation on an oracle Uφ for a non-boolean function φ, comprising: a two-register quantum system comprising an ancilla qubit, andan input qubit; anda non-boolean quantum oracle comprising a selective phase-flip unitary operator circuit SΨ0,a controlled unitary operator circuit Uφ, anda controlled inverse unitary operator circuit Uφ†;wherein the quantum computing device is configured to initialize the ancilla qubit in a |+ state and the input qubit in a |ψ0 state of a plurality of eigenstates |x to define a two-register state |Ψ0≡|+,ψ0;for each of a plurality K of iterations receive, using the input qubit, a respective one of the plurality of eigenstates |x defining an input basis state,for odd iterations of the plurality K of iterations, act on the input basis state using the selective phase-flip unitary operator circuit SΨ0 and the controlled unitary operator circuit Uφ, andfor even iterations of the plurality K of iterations, act on the input basis state using the selective phase-flip unitary operator circuit SΨ0 and the controlled inverse unitary operator circuit Uφ†.

9. The quantum computing device according to claim 8, further configured to measure, after the plurality K of iterations, the ancilla qubit in a 0/1 basis.

10. The quantum computing device according to claim 8, wherein the selective phase-flip unitary operator circuit SΨ0 further comprises: a Hadamard transform H,a unitary operator A0,a unitary operator I, andan inverse unitary operator A0†;wherein the selective phase-flip unitary operator SΨ0=[H⊗A0][2|0,00,0|−I][H⊗A0†].

11. The quantum computing device according to claim 8, further configured, for the ancilla qubit in a state |0, to act on the input qubit as Uφ|0,x=e+iφ(x)|0,x, Uφ|1,x=e−iφ(x)|1,x.

12. The quantum computing device according to claim 8, further configured, for the ancilla qubit in a state |1, to act on the input qubit as Uφ†|0,x=e+iφ(x)|0,x, Uφ†|1,x=e−iφ(x)|1,x.

13. The quantum computing device according to claim 8, further comprising a controlled estimation unitary operator circuit Uφ−π/2, and configured to: receive, using the input qubit, the |ψ0 state defining an input random state; andact on the input random state using the controlled estimation unitary operator circuit Uφ−π/2.

14. The quantum computing device according to claim 13, wherein the controlled estimation unitary operator circuit Uφ−π/2 further comprises: the controlled unitary operator circuit Uφ,at least one bit-flip operator X, andat least one phase-shift operator Rϕ;wherein the controlled estimation unitary operator circuit Uφ−π/2=e−iπ/2|00|⊗Uφ+eiπ/2|11|⊗Uφ†.

15. A system of quantum circuits for implementing an oracle U q for a non-boolean function φ, the system configured to: initialize an ancilla qubit in a |+ state and an input qubit in a |ψ0 state of a plurality of eigenstates |x to define a two-register state |Ψ0≡|+,ψ0;for each of a plurality K of iterations receive, using the input qubit, a respective one of the plurality of eigenstates |x defining an input basis state,for odd iterations of the plurality K of iterations, act on the input basis state using a selective phase-flip unitary operator circuit SΨ0 and a controlled unitary operator circuit Uφ, andfor even iterations of the plurality K of iterations, act on the input basis state using the selective phase-flip unitary operator circuit SΨ0 and a controlled inverse unitary operator circuit Uφ†.

16. The system according to claim 15, further configured to measure, after the plurality K of iterations, the ancilla qubit in a 0/1 basis.

17. The system according to claim 15, wherein the selective phase-flip unitary operator circuit SΨ0 further comprises: a Hadamard transform H,a unitary operator A0,a unitary operator I, andan inverse unitary operator A0†;wherein the selective phase-flip unitary operator SΨ0=[H⊗A0][2|0,00,0|−I][H⊗A0†].

18. The system according to claim 15, further configured, for the ancilla qubit in a state |0, to act on the input qubit as Uφ|0,x=e+iφ(x)|0,x, Uφ|1,x=e−iφ(x)|1,x.

19. The system according to claim 15, further configured, for the ancilla qubit in a state |1, to act on the input qubit as Uφ†|0,x=e+iφ(x)|0,x, Uφ†|1,x=e−iφ(x)|1,x.

20. The system according to claim 15, further configured to: receive, using the input qubit, the |ψ0 state defining an input random state; andact on the input random state using a controlled estimation unitary operator circuit Uφ−π/2 comprising: the controlled unitary operator circuit Uφ,at least one bit-flip operator X, andat least one phase-shift operator Rϕ;wherein the controlled estimation unitary operator circuit Uφ−π/2=e−iπ/2|00|⊗Uφ+eiπ/2|11|⊗Uφ†.