View source for PostBQP

In [[computational complexity theory]], '''PostBQP''' is a [[complexity class]] consisting of all of the [[computational problem]]s solvable in [[polynomial time]] on a [[quantum Turing machine]] with [[postselection]] and bounded error (in the sense that the algorithm is correct at least 2/3 of the time on all inputs).

Postselection is not considered to be a feature that a realistic computer (even a quantum one) would possess, but nevertheless postselecting machines are interesting from a theoretical perspective.

Removing either one of the two main features (quantumness, postselection) from '''PostBQP''' gives the following two complexity classes, both of which are subsets of '''PostBQP''':
* '''[[BQP]]''' is the same as '''PostBQP''' except without postselection
* '''BPP<sub>path</sub>''' is the same as '''PostBQP''' except that instead of quantum, the algorithm is a classical randomized algorithm (with postselection)<ref>{{cite journal |year=1997|title=Threshold computation and cryptographic security|journal= [[SIAM Journal on Computing]] | volume=26 | pages=59–78 | doi = 10.1137/S0097539792240467 | author=Y. Han and Hemaspaandra, L. and Thierauf, T.|citeseerx=10.1.1.23.510}}</ref>

The addition of [[postselection]] seems to make quantum Turing machines much more powerful: [[Scott Aaronson]] proved<ref name="PostBQP = PP">{{cite journal|last=Aaronson|first=Scott|year=2005|title=Quantum computing, postselection, and probabilistic polynomial-time|journal=[[Proceedings of the Royal Society A]] |volume=461|issue=2063|pages=3473–3482|doi=10.1098/rspa.2005.1546|arxiv=quant-ph/0412187|bibcode=2005RSPSA.461.3473A|s2cid=1770389}}.  Preprint available at [https://arxiv.org/abs/quant-ph/0412187]</ref><ref>{{cite web|url=http://weblog.fortnow.com/2004/01/complexity-class-of-week-pp-by-guest.html|title=Complexity Class of the Week: PP|last=Aaronson|first=Scott|date=2004-01-11|work=Computational Complexity Weblog|access-date=2008-05-02}}</ref> '''PostBQP''' is equal to '''[[PP (complexity)|PP]]''', a class which is believed to be relatively powerful, whereas '''[[BQP]]''' is not known even to contain the seemingly smaller class '''[[NP (complexity class)|NP]]'''. Using similar techniques, Aaronson also proved that small changes to the laws of quantum computing would have significant effects. As specific examples, under either of the two following changes, the "new" version of '''BQP''' would equal '''[[PP (complexity)|PP]]''':
* if we broadened the definition of 'quantum gate' to include not just unitary operations but linear operations, or
* if the probability of measuring a basis state <math>|x\rangle</math> was proportional to <math>|\alpha_x|^p</math> instead of <math>|\alpha_x|^2</math> for any even integer ''p > 2''.
{{TOClimit|limit=3}}

==Basic properties==
In order to describe some of the properties of '''PostBQP''' we fix a formal way of describing quantum postselection. Define a quantum algorithm to be a family of [[quantum circuit]]s (specifically, a [[uniform circuit family]]). We designate one qubit as the ''postselection qubit P'' and another as the ''output qubit Q''. Then '''PostBQP''' is defined by postselecting upon the event that the postselection qubit is <math>|1\rangle</math>. Explicitly, a language ''L'' is in PostBQP if there is a quantum algorithm ''A'' so that after running ''A'' on input ''x'' and measuring the two qubits ''P'' and ''Q'',
* ''P'' = 1 with nonzero probability
* if the input ''x'' is in ''L'' then {{math|1 = Pr[''Q'' = 1{{!}}''P'' = 1] ≥ 2/3}}
* if the input ''x'' is not in ''L'' then {{math|1 = Pr[''Q'' = 0{{!}}''P'' = 1] ≥ 2/3}}.

One can show that allowing a single postselection step at the end of the algorithm (as described above) or allowing intermediate postselection steps during the algorithm are equivalent.<ref name="PostBQP = PP"></ref><ref>{{cite journal|
    title = Quantum Complexity Theory|
    journal = SIAM Journal on Computing|
    year = 1997|
    volume = 26|
    issue = 5|
    pages = 11&ndash;20|author1=Ethan Bernstein  |author2=Umesh Vazirani
 |name-list-style=amp |doi=10.1137/s0097539796300921
|
    citeseerx = 10.1.1.144.7852}}</ref>

Here are three basic properties of '''PostBQP''' (which also hold for '''BQP''' via similar proofs):

# '''PostBQP''' is ''closed under complement''. Given a language ''L'' in '''PostBQP''' and a corresponding deciding circuit family, create a new circuit family by flipping the output qubit after measurement, then the new circuit family proves the complement of ''L'' is in '''PostBQP'''.
# You can do ''probability amplification'' in '''PostBQP'''. The definition of '''PostBQP''' is not changed if we replace the 2/3 value in its definition by any other constant strictly between 1/2 and 1. As an example, given a '''PostBQP''' algorithm ''A'' with success probability 2/3, we can construct another algorithm which runs three independent copies of ''A'', outputs a postselection bit equal to the [[logical conjunction|conjunction]] of the three "inner" ones, and outputs an output bit equal to the [[majority]] of the three "inner" ones; the new algorithm will be correct with conditional probability <math>(2/3)^3 + 3(1/3)(2/3)^2 = 20/27</math>, greater than the original 2/3.
# '''PostBQP''' is ''closed under [[intersection (set theory)|intersection]]''. Suppose we have '''PostBQP''' circuit families for two languages {{tmath|L_1}} and {{tmath|L_2}}, with respective postselection qubits and output qubits {{tmath|P_1, P_2, Q_1, Q_2}}. We may assume by probability amplification that both circuit families have success probability at least 5/6. Then we create a composite algorithm where the circuits for {{tmath|L_1}} and {{tmath|L_2}} are run independently, and we set ''P'' to the conjunction of {{tmath|P_1}} and {{tmath|P_2}}, and ''Q'' to the conjunction of {{tmath|Q_1}} and {{tmath|Q_2}}. It is not hard to see by a [[union bound]] that this composite algorithm correctly decides membership in <math>L_1 \cap L_2</math> with (conditional) probability at least 2/3.

More generally, combinations of these ideas show that '''PostBQP''' is closed under union and BQP truth-table reductions.

=={{sans-serif|PostBQP}} = {{sans-serif|PP}}==

[[Scott Aaronson]] showed<ref>{{cite journal | doi = 10.1098/rspa.2005.1546 | last1 = Aaronson | first1 = Scott|author9 = Aaronson, Scott | year = 2005 | title = Quantum computing, postselection, and probabilistic polynomial-time| journal = Proceedings of the Royal Society A | volume = 461 | issue = 2063 | pages = 3473&ndash;3482 | arxiv = quant-ph/0412187| bibcode = 2005RSPSA.461.3473A| s2cid = 1770389 }}</ref> that the complexity classes {{tmath|\mathsf{PostBQP} }} (postselected bounded error quantum polynomial time) and '''[[PP (complexity)|PP]]''' (probabilistic polynomial time) are equal. The result was significant because this [[PP (complexity)#PostBQP|quantum computation reformulation of {{tmath|\mathsf{PP} }}]] gave new insight and simpler proofs of properties of {{tmath|\mathsf{PP} }}.

The usual definition of a {{tmath|\mathsf{PostBQP} }} circuit family is one with two outbit qubits ''P'' (postselection) and ''Q'' (output) with a single measurement of ''P'' and ''Q'' at the end such that the probability of measuring {{math|1=''P'' = 1}} has nonzero probability, the conditional probability <span class="texhtml">Pr[''Q'' = 1|''P'' = 1] ≥ 2/3</span> if the input&nbsp;''x'' is in the language, and <span class="texhtml">Pr[''Q'' = 0|''P'' = 1] ≥ 2/3</span> if the input ''x'' is not in the language. For technical reasons we tweak the definition of {{tmath|\mathsf{PostBQP} }} as follows: we require that {{math|1=Pr[''P'' = 1] ≥ 2<sup>&minus;''n''<sup>''c''</sup></sup>}} for some constant ''c'' depending on the circuit family. Note this choice does not affect the [[PostBQP#Basic properties|basic properties of {{tmath|\mathsf{PostBQP} }}]], and also it can be shown that any computation consisting of typical gates (e.g. Hadamard, Toffoli) has this property whenever {{math|1=Pr[''P'' = 1] > 0}}.

== Proving {{sans-serif|PostBQP}} &sube; {{sans-serif|PP}} ==

Suppose we are given a {{tmath|\mathsf{PostBQP} }} family of circuits to decide a language ''L''. We assume without loss of generality (e.g. see the [[Quantum computer#EssentialProperties|inessential properties of quantum computers]]) that all gates have transition matrices that are represented with real numbers, at the expense of adding one more qubit.

Let {{math|&Psi;}} denote the final quantum state of the circuit before the postselecting measurement is made. The overall goal of the proof is to construct a {{tmath|\mathsf{PP} }} algorithm to decide ''L''. More specifically it suffices to have ''L'' correctly compare the squared amplitude of {{math|&Psi;}} in the states with {{math|1=''Q'' = 1, ''P'' = 1}} to the squared amplitude of {{math|&Psi;}} in the states with {{math|1=''Q'' = 0, ''P'' = 1}} to determine which is bigger. The key insight is that the comparison of these amplitudes can be transformed into comparing the acceptance probability of a {{tmath|\mathsf{PP} }} machine with 1/2.

=== Matrix view of PostBQP algorithms ===
Let ''n'' denote the input size, {{math|1=''B'' = ''B''(''n'')}} denote the total number of qubits in the circuit (inputs, ancillary, output and postselection qubits), and {{math|1=''G'' = ''G''(''n'')}} denote the total number of gates.
Represent the ''i''th gate by its transition matrix ''A''<sup>''i''</sup> (a real unitary <math>2^B \times 2^B</math> matrix) and let the initial state be <math>|x\rangle</math> (padded with zeroes). Then <math>\Psi = A^G A^{G-1}\dotsb A^2 A^1 |x\rangle</math>. Define ''S<sub>1</sub>'' (resp. ''S''<sub>0</sub>) to be the set of basis states corresponding to {{math|1=''P'' = 1, ''Q'' = 1}} (resp. {{math|1=''P'' = 1, ''Q'' = 0}}) and define the probabilities

: <math>\pi_1 := \text{Pr}[P=1,Q=1]=\sum_{\omega \in S_1} \Psi^2_\omega</math>

: <math>\pi_0 := \text{Pr}[P=1,Q=0]=\sum_{\omega \in S_0} \Psi^2_\omega.</math>

The definition of {{tmath|\mathsf{PostBQP} }} ensures that either <math>\pi_{1} \ge 2\pi_0</math> or <math>\pi_0 \ge 2\pi_1</math> according to whether ''x'' is in ''L'' or not.

Our {{tmath|\mathsf{PP} }} machine will compare <math>\pi_{1}</math> and <math>\pi_{0}</math>. In order to do this, we expand the definition of matrix multiplication:
:: <math>\Psi_\omega = \sum_{\alpha_1, \ldots, \alpha_{G}}A^{G}_{\omega,\alpha_{G}}A^{G-1}_{\alpha_G,\alpha_{G-1}}\dotsb A^2_{\alpha_3,\alpha_2} A^1_{\alpha_2,\alpha_1} x_{\alpha_1}</math>

where the sum is taken over all lists of ''G'' basis vectors <math>\alpha_i</math>. Now <math>\pi_1</math> and <math>\pi_0</math> can be expressed as a sum of pairwise products of these terms. Intuitively, we want to design a machine whose acceptance probability is something like <math>\tfrac{1}{2}(1+\pi_1-\pi_0)</math>, since then <math>x \in L</math> would imply that the acceptance probability is <math>\tfrac{1}{2}(1+\pi_{1}-\pi_{0})>\tfrac 1 2</math>, while <math>x \not\in L</math> would imply that the acceptance probability is <math>\tfrac{1}{2}(1+\pi_1-\pi_0)<\tfrac 1 2</math>.

====Technicality: we may assume entries of the transition matrices {{mvar|A<sup>i</sup>}} are rationals with denominator {{math|2{{sup|''f''(''n'')}}}} for some polynomial ''f''(''n'').====
The definition of {{tmath|\mathsf{PostBQP} }} tells us that <math>\pi_{1} \ge \tfrac{2}{3}(\pi_0+\pi_1)</math> if ''x'' is in ''L'', and that otherwise <math>\pi_{0} \ge \tfrac{2}{3}(\pi_0+\pi_1)</math>. Let us replace all entries of ''A'' by the nearest fraction with denominator <math>2^{f(n)}</math> for a large polynomial {{tmath|f(n)}} that we presently describe. What will be used later is that the ''new'' {{pi}} values satisfy <math>\pi_1 > \tfrac{1}{2}(\pi_0+\pi_1)</math> if ''x'' is in ''L'', and <math>\pi_0 > \tfrac{1}{2}(\pi_0+\pi_1)</math> if ''x'' is not in ''L''. Using the earlier technical assumption and by analyzing how the 1-norm of the computational state changes, this is seen to be satisfied if <math>(1+2^{-f(n)}2^{B})^{G}-1 < \tfrac{1}{6}2^{-n^c},</math> thus clearly there is a large enough ''f'' that is polynomial in ''n''.

===Constructing the PP machine===
Now we provide the detailed implementation of our {{tmath|\mathsf{PP} }} machine. Let {{mvar|&alpha;}} denote the sequence <math>\{\alpha_i\}_{i=1}^G</math> and define the shorthand notation
:: <math>\Pi(A, \omega, \alpha, x) := A^{G}_{\omega,\alpha_{G}}A^{G-1}_{\alpha_{G},\alpha_{G-1}}\dotsb A^2_{\alpha_3,\alpha_2} A^1_{\alpha_2,\alpha_1} x_{\alpha_1}</math>,
then 
:: <math>\pi_1 - \pi_0 = \sum_{\omega \in S_1} \sum_{\alpha,\alpha'} \Pi(A, \omega, \alpha, x)\Pi(A, \omega, \alpha', x) - \sum_{\omega \in S_0} \sum_{\alpha,\alpha'} \Pi(A, \omega, \alpha, x)\Pi(A, \omega, \alpha', x).</math>
We define our {{tmath|\mathsf{PP} }} machine to
* pick a basis state {{mvar|&omega;}} uniformly at random
* if <math>\omega \not\in S_0 \cup S_1</math> then STOP and accept with probability 1/2, reject with probability 1/2
* pick two sequences <math>\alpha,\alpha'</math> of ''G'' basis states uniformly at random
* compute <math>X = \Pi(A, \omega, \alpha, x)\Pi(A, \omega, \alpha', x)</math> (which is a fraction with denominator <math>2^{2f(n)G(n)}</math> such that <math>-1 \le X \le 1</math>)
* if <math>\omega \in S_1</math> then accept with probability <math>\tfrac{1+X}{2}</math> and reject with probability <math>\tfrac{1-X}{2}</math> (which takes at most {{tmath|2f(n)G(n)+1}} coin flips)
* otherwise (then <math>\omega \in S_0</math>) accept with probability <math>\tfrac{1-X}{2}</math> and reject with probability <math>\tfrac{1+X}{2}</math> (which again takes at most {{tmath|2f(n)G(n)+1}} coin flips)
Then it is straightforward to compute that this machine accepts with probability
<math>\frac{1}{2}+\frac{\pi_{1}-\pi_{0}}{2^{1+B(n)+2B(n)G(n)}},</math>
so this is a {{tmath|\mathsf{PP} }} machine for the language ''L'', as needed.

== Proving {{sans-serif|PP}} &sube;  {{sans-serif|PostBQP}} ==
Suppose we have a {{tmath|\mathsf{PP} }} machine with time complexity {{tmath|1=T:=T(n)}} on input ''x'' of length <math>n := |x|</math>. Thus the machine flips a coin at most ''T'' times during the computation. We can thus view the machine as a deterministic function ''f'' (implemented, e.g. by a classical circuit) which takes two inputs (''x, r'') where ''r'', a binary string of length ''T'', represents the results of the random coin flips that are performed by the computation, and the output of ''f'' is 1 (accept) or 0 (reject). The definition of {{tmath|\mathsf{PP} }} tells us that 
:: <math>x \in L \Leftrightarrow \#\{r \in \{0,1\}^T\mid f(x, r)=1\} > 2^{T-1}</math>
Thus, we want a {{tmath|\mathsf{PostBQP} }} algorithm that can determine whether the above statement is true.

Define ''s'' to be the number of random strings which lead to acceptance,
:: <math>s := \#\{r \in \{0,1\}^T\mid f(x, r)=1\}</math>
and so <math>2^T-s</math> is the number of rejected strings.
It is straightforward to argue that without loss of generality, <math>s \not\in \{0, 2^T/2, 2^T\}</math>; for details, see a similar [[without loss of generality]] assumption in the [[PP (complexity)#complement|proof that {{tmath|\mathsf{PP} }} is closed under complementation]].

=== Aaronson's algorithm ===
Aaronson's algorithm for solving this problem is as follows. For simplicity, we will write all quantum states as unnormalized. First, we prepare the state <math>|x \rangle \otimes \sum_{r \in \{0,1\}^T} |r \rangle |f(x,r) \rangle</math>. Second, we apply [[Hadamard gate]]s to the second register (each of the first ''T'' qubits), measure the second register and postselect on it being the all-zero string. It is easy to verify that this leaves the last register (the last qubit) in the residual state
:: <math> |\psi \rangle := (2^T-s)|0 \rangle + s|1 \rangle. </math>
Where ''H'' denotes the Hadamard gate, we compute the state
:: <math> H |\psi\rangle = (2^T |0\rangle + (2^T - 2s)|1 \rangle)/\sqrt{2} </math>.
Where <math>\alpha, \beta</math> are positive real numbers to be chosen later with <math>\alpha^2+\beta^2=1</math>, we compute the state <math>\alpha |0\rangle|\psi\rangle+\beta |1\rangle|H\psi\rangle</math> and measure the second qubit, postselecting on its value being equal to 1, which leaves the first qubit in a residual state depending on <math>\beta/\alpha</math> which we denote
:: <math> | \phi_{\beta/\alpha} \rangle := \alpha s |0\rangle + \frac{\beta}{\sqrt{2}}(2^T-2s)|1\rangle</math>.

Visualizing the possible states of a qubit as a circle, we see that if <math>s > 2^{T-1}</math>, (i.e. if <math>x \in L</math>) then <math>\phi_{\beta/\alpha}</math> lies in the open quadrant <math>Q_{acc} := (-|1\rangle, |0\rangle)</math> while if <math>s < 2^{T-1}</math>, (i.e. if <math>x \not\in L</math>) then <math>\phi_{\beta/\alpha}</math> lies in the open quadrant <math>Q_{rej} := (|0\rangle,|1\rangle)</math>. In fact for any fixed ''x'' (and its corresponding ''s''), as we vary the ratio <math>\beta/\alpha</math> in <math>(0, \infty)</math>, note that the image of <math>|\phi_{\beta/\alpha}\rangle </math> is precisely the corresponding open quadrant. In the rest of the proof, we make precise the idea that we can distinguish between these two quadrants.

===Analysis===
Let <math>|+\rangle = (|1\rangle+|0\rangle)/\sqrt{2}</math>, which is the center of <math>Q_{rej}</math>, and let <math>|-\rangle</math> be orthogonal to <math>|+\rangle</math>. Any qubit in <math>Q_{acc}</math>, when measured in the basis <math>\{|+\rangle, |-\rangle\}</math>, gives the value <math>|+\rangle</math> less than 1/2 of the time. On the other hand, if <math>x \not\in L</math> and we had picked <math>\beta/\alpha = r^* := \sqrt{2}s / (2^T-2s)</math> then measuring <math>| \phi_{\beta/\alpha} \rangle</math> in the basis <math>\{|+\rangle, |-\rangle\}</math> would give the value <math>|+\rangle</math> all of the time. Since we don't know ''s'' we also don't know the precise value of ''r*'', but we can try several (polynomially many) different values for <math>\beta/\alpha</math> in hopes of getting one that is "near" ''r*''.

Specifically, note <math>2^{-T} < r* < 2^T</math> and let us successively set <math>\beta/\alpha</math> to every value of the form <math>2^i</math> for <math>-T \leq i \leq T</math>. Then elementary calculations show that for one of these values of ''i'', the probability that the measurement of <math>| \phi_{2^i} \rangle</math> in the basis <math>\{|+\rangle, |-\rangle\}</math> yields <math>|+\rangle</math> is at least <math>(3+2\sqrt{2})/6 \approx 0.971.</math>

Overall, the {{tmath|\mathsf{PostBQP} }} algorithm is as follows. Let ''k'' be any constant strictly between 1/2 and <math>(3+2\sqrt{2})/6</math>.
We do the following experiment for each <math>-T \leq i \leq T</math>: construct and measure <math>| \phi_{2^i} \rangle</math> in the basis <math>\{|+\rangle, |-\rangle\}</math> a total of <math>C \log T</math> times where ''C'' is a constant. If the proportion of <math>|+\rangle</math> measurements is greater than ''k'', then reject. If we don't reject for any ''i'', accept. [[Chernoff bound]]s then show that for a sufficiently large universal constant ''C'', we correctly classify ''x'' with probability at least 2/3.

Note that this algorithm satisfies the technical assumption that the overall postselection probability is not too small: each individual measurement of <math>| \phi_{2^i} \rangle</math> has postselection probability <math>1/2^{O(T)}</math> and so the overall probability is <math>1/2^{O(T^2 \log T)}</math>.

== References ==
<references/>

{{quantum computing}}

{{DEFAULTSORT:Postbqp}}
[[Category:Articles containing proofs]]
[[Category:Quantum complexity theory]]
[[Category:Probabilistic complexity classes]]