Write
153
edits
Full Quantum Process Tomography with Linear inversion: Difference between revisions
Jump to navigation
Jump to search
Full Quantum Process Tomography with Linear inversion (edit)
Revision as of 10:44, 11 September 2019
, 11 September 2019no edit summary
(Created page with "Quantum process tomography is a method used to characterize a physical process in an open quantum system. This method is similar to Quantum state tomography [link], but the go...") |
No edit summary |
||
| Line 21: | Line 21: | ||
* Measurement device. | * Measurement device. | ||
* Quantum computational resources to perform the operation. | * Quantum computational resources to perform the operation. | ||
==Notation== | ==Notation== | ||
| Line 37: | Line 36: | ||
==Properties== | ==Properties== | ||
* '''Figure of merit''': Density Matrix of the quantum process | |||
* Multiple copies of the quantum state are required in this method. | |||
* The process <math>\varepsilon</math> can be completely determined by <math>\chi</math> once the set of operators <math>\tilde{A_i}</math> has been fixed. | |||
* <math>\chi</math> will contain <math>N^4 - N^2</math> parameters | |||
* In the case of a “nonselective” quantum evolution, such as arises from uncontrolled interactions with an environment the <math>A_i</math> operators satisfy an additional completeness relation <math>\sum_i A_i^{\dagger} A_i = I</math>. This relation ensures that the trace factor <math>tr(\varepsilon(\rho))</math> is always equal to one, and thus the state change experienced by the system can be written <math>\rho \xrightarrow{} \varepsilon(\rho)</math> | |||
* Quantum state tomography has to be performed multiple times in this experiment. | |||
* This method is extremely resource intensive. | |||
* The main drawback of this technique is that sometimes the computed solution of the density matrix will not be a valid density matrix. This issue occurs when not enough measurements are made. Due to this drawback, the maximum likelihood estimation technique [link here] is preferred over this linear inversion technique. | |||
* Another drawback of this method is that an infinite number of measurement outcomes would be required to give the exact solution. | |||
* Entanglement fidelity can be evaluated with the knowledge of <math>A_i</math> operators. This quantity can be used to measure how closely the dynamics of the quantum system under consideration approximates that of some ideal quantum system. This quantity can be determined robustly, because of its linear dependence on the experimental errors. Suppose the target quantum operation is a unitary quantum operation <math>U(\rho) = U\rho U^\dagger</math> and the actual quantum operation implemented experimentally is <math>\varepsilon</math>. Then the entanglement fidelity is <math>F_e(\rho, U, \varepsilon) = \sum_i {|tr(U^{\dagger}A_i \rho )|}^2 = \sum_{mn} \chi_{mn} tr(U^{\dagger}\tilde{A_m}\rho) tr(\rho \tilde{A_n^{\dagger}} U)</math> | |||
* The minimum value of <math>F_e</math> over all possible states <math>\rho</math> is a single parameter which describes how well the experimental system implements the desired quantum logic gate | |||
* Minimum fidelity of the gate operation is <math>F = min_{|\psi \rangle} \langle\psi| U^{\dagger} \varepsilon (|\psi\rangle\langle\psi|) U|\psi\rangle</math> | |||
* Quantum channel capacity which is a measure of the amount of quantum information that can be sent reliably using a quantum communications channel which is described by a quantum operation <math>\varepsilon</math> can also be defined. | |||
* This procedure can also be used to determine the form of the [[Lindblad operator]] used in Markovian master equations of a certain form. | |||
* Quantum operations can also be used to describe measurements. For each measurement outcome, <math>i</math>, there is a associated quantum operation <math>\varepsilon_i</math>. The corresponding state change is given by <math>\rho \xrightarrow{} \varepsilon_i(\rho)/tr(\varepsilon_i(\rho))</math> where the probability of the measurement outcome occurring is <math>p_i = tr(\varepsilon_i(\rho))</math>. To determine the process the steps are the same, except the measurement has to be taken a large enough number of times that the probability pi can be reliably estimated. Next <math>\rho^{'}_{j}</math> is determined using tomography, thus obtaining <math>\varepsilon_i(\rho_j) = tr(\varepsilon_i(\rho_j))\rho^{'}_{j}</math> for each input <math>\rho_j</math>. The other steps remain to same to estimate <math>\varepsilon_i</math>. | |||
==Protocol Description== | ==Protocol Description== | ||
'''Input''': <math>p_j, j = 1, ..., N^2</math> | |||
'''Output''': Density matrix of the quantum process operator, <math>\varepsilon</math> | |||
* Select a fixed set of operators <math>\tilde{A_i}</math> such that <math>A_i = \sum_m a_{im} \tilde{A_m}</math> | |||
* The operator sum representation is <math>\varepsilon(\rho) = \sum_{i}A_i\rho A^{\dagger}_i</math> | |||
* Hence, <math>\varepsilon(\rho) = \sum_{mn}\tilde{A_m}\rho \tilde{A^{\dagger}_n}\chi_{mn}</math> | |||
* For <math>\rho_j = 1, 2, ..., N^2</math>: | |||
** The process <math>\varepsilon</math> is performed on <math>\rho_j</math> | |||
** The output state <math>\varepsilon(\rho_j)</math> is measured using quantum state tomography | |||
** <math>\varepsilon(\rho_j)</math> is expressed as a linear combination of basis states, <math>\varepsilon(\rho_j) = \sum_k \lambda_{jk} \rho_k</math> | |||
** Since <math>\tilde{A_m}\rho_j\tilde{A_n^{\dagger}} = \sum_k \beta^{mn}_{jk} \rho_k</math>, then <math>\sum_k \sum_{mn} \chi_{mn} \beta^{mn}_{jk} \rho_k = \sum_k \lambda_{jk} \rho_k</math> | |||
** Since <math>\rho_k</math> is independent, <math>\lambda_{jk} = \sum_{mn} \beta^{mn}_{jk} \chi_{mn}</math> | |||
* From there <math>\chi_{mn}</math> is defined as, <math>\chi_{mn} = \sum_{jk} \kappa^{mn}_{jk}</math> | |||
* Let the <math>U^{\dagger}</math> diagonalize <math>\chi</math>, <math>\chi_{mn} =\sum_{xy} U_{mn}d_x\delta_{xy}U^{*}_{ny}</math> | |||
* From there, <math>A_i = \sqrt{d_i} \sum_j U_{ij} \tilde{A_j}</math> | |||
* Thus <math>\varepsilon(\rho)</math> can be determined using <math>\varepsilon(\rho) = \sum_{i}A_i\rho A^{\dagger}_i</math> | |||
==Further Information== | ==Further Information== | ||
* [https://arxiv.org/abs/quant-ph/0402166 Maximum likelihood estimation] (also known as MLE or MaxLik) is a popular technique that overcomes the problem that the naive matrix inversion procedure in QPT, when performed on real (i.e., inherently noisy) experimental data, typically leads to an unphysical process matrix. | |||
* The standard approach of process tomography is grossly inaccurate in the case where the states and measurement operators used to interrogate the system are generated by gates that have some systematic error, a situation all but unavoidable in any practical setting. These errors in tomography cannot be fully corrected through oversampling or by performing a larger set of experiments. Hence [[gate set tomography]] was introduced. | |||
==Related Papers== | ==Related Papers== | ||
* Chuang et al arXiv:quant-ph/9610001v1: Prescription for experimental determination of | |||
the dynamics of a quantum black box | |||
* J. F. Poyatos, J. I. Cirac, and P. Zoller, PhysRevLett.78.390: Complete Characterization of a Quantum Process: The Two-Bit Quantum Gate | |||
* M.W. Mitchell, et al., Phys. Rev. Lett. 91, 120402 (2003) | |||
* J. L. O’Brien et al arXiv:quant-ph/0402166v2: Quantum process tomography of a controlled-not gate | |||
* S. T. Merkel, J. M. Gambetta, J. A. Smolin, S. Poletto, A. D. Corcoles, B. R. | |||
* Johnson, C. A. Ryan, and M. Steffen, Phys. Rev. A 87, 062119 (2013). | |||
<div style='text-align: right;'>''*contributed by Rhea Parekh''</div> | <div style='text-align: right;'>''*contributed by Rhea Parekh''</div> | ||