Protec'ng quantum gates from control noise Constantin Brif - PowerPoint PPT Presentation

Protec'ng ¡quantum ¡gates ¡ from ¡control ¡noise ¡ Constantin Brif Sandia National Laboratories Collaborators: Matthew Grace and Kevin Young (Sandia) David Hocker, Katharine Moore,Tak-San Ho, and Herschel Rabitz (Princeton University) Sandia National Laboratories is a multi-program laboratory managed and operated by Sandia Corporation, a wholly owned subsidiary of Lockheed Martin Corporation, for the U.S. Department of Energy's National Nuclear Security Administration under contract DE-AC04-94AL85000.

Basic ¡defini'ons: ¡Unitary ¡quantum ¡gates ¡ A unitary quantum gate is the basic functioning element of a quantum circuit. Some basic notation: number of qubits in the quantum gate system dimension of the system’s Hilbert space the target unitary transformation the actual evolution operator at the final time T The same unitary transformation is applied to any input state :

Controlled ¡quantum ¡gates ¡ An external classical control c ( t ) is necessary to operate the quantum gate. The Hamiltonian and evolution operator are functionals of the control: The gate fidelity is a measure of how well the target transformation is performed: It is convenient to use a normalized fidelity: or ¡ The gate fidelity is also a functional of the control:

Quantum ¡control ¡landscape ¡and ¡op'mality ¡ The functional dependence F = F [ c ( t )] is called the control landscape . Critical points of the control landscape satisfy: F A sufficient condition for optimality of a critical point is negative semidefiniteness of the Hessian: For a recent review, see C. Brif, R. Chakrabarti, and H. Rabitz, New J. Phys. 12 , 075008 (2010)

Op'mally ¡controlled ¡quantum ¡gates ¡ The analysis of regular critical points on the control landscape reveals that: l There is one maximum manifold: F = 1 l There is one minimum manifold: F = 0 l All other critical manifolds are saddles (can be avoided by a smart optimization algorithm) An optimal control solution c 0 ( t ) is perfect in ideal conditions: For quantum gate control with linear field coupling: The Hessian at the maximum: The Hessian at any optimal control solution has only non-positive eigenvalues. The “flatness” of the control landscape in the vicinity of an optimal solution depends on the number and magnitude of negative Hessian eigenvalues.

Op'mal ¡quantum ¡gates ¡with ¡noisy ¡controls ¡ All actual controls are noisy! Control noise Gate errors F < 1 Consider a quantum gate that operates in the vicinity of an optimal control: l additive noise : l multiplicative noise : Expanding for small noise: We consider a random noise process, so the error z ( t ) is a stochastic variable with the autocorrelation function : Expected value of the quantum gate fidelity:

Robustness ¡to ¡white ¡noise: ¡general ¡results ¡ White noise has zero correlation time: Expected value of the quantum gate fidelity: The diagonal elements of the optimal Hessian are time independent: The general expression for the expected value of the quantum gate fidelity in the presence of small white noise :

Robustness ¡to ¡addi've ¡white ¡noise ¡(AWN) ¡ Expected gate error is proportional to the total control time T : 0 10 2 = 10 − 1 σ a − 1 10 2 = 10 − 2 σ a − 2 10 Expected gate error 2 = 10 − 3 σ a − 3 10 2 = 10 − 4 σ a − 4 10 2 = 10 − 5 σ a − 5 10 2 = 10 − 6 σ a − 6 10 − 7 10 0.5 0.75 1 1.25 1.5 2 2.5 3 4 5 6 7 8 9 10 Control time For fidelities above 0.99, the discrepancy between the perturbative result and Monte Carlo averaging (over a sample of 20,000 noisy controls) is within the statistical error due to the finite sample size (~0.7%). The accuracy of the perturbative approximation is excellent for all fidelities above 0.9.

Robustness ¡to ¡mul'plica've ¡white ¡noise ¡(MWN) ¡ Expected gate error is proportional to the control fluence (energy): 0 10 2 = 10 − 1 σ m − 1 10 2 = 10 − 2 σ m − 2 10 2 = 10 − 3 Expected gate error σ m − 3 10 2 = 10 − 4 σ m − 4 10 2 = 10 − 5 σ m − 5 10 2 = 10 − 6 σ m − 6 10 − 7 10 0.8 1 1.2 1.5 2 2.5 3 4 5 6 7 8 9 10 12 15 18 Fluence of the optimal control field The perturbative result is once again in excellent agreement with Monte Carlo averaging (over a sample of 20,000 noise realizations) for all gate fidelities above 0.9.

Time-­‑op'mal ¡control: ¡improved ¡robustness ¡to ¡AWN ¡ Optimizing robustness to additive white noise in controls is equivalent to minimizing the total control time. For a given system Hamiltonian H and target gate W , there exists a critical value T * of control time, below which the target is no longer reachable. We developed a numerical procedure to identify T * and explore the Pareto front for two competing control objectives: gate fidelity maximization and control time minimization. 10 −2 Time-fidelity Pareto fronts for the CNOT T=6 gate in a two-qubit system, explored 10 −3 T=4.2 starting from different initial values of T . T=4.4 T=4.3 D (normalized) T=4.6 10 −4 For control times above T * , it is possible T=4.5 T=5 to decrease T (and thus improve T=4.4 10 −5 robustness) while keeping the nominal T=4.8 gate fidelity at 1. However, below T * , the T=4.2 nominal fidelity rapidly deteriorates as T 10 −6 decreases, i.e., one has to sacrifice fidelity to improve robustness. 10 −7 10 −8 K. W. Moore et al., arXiv:1112.0333 (2011) 3 3.2 3.4 3.6 3.8 4 4.2 T T*

Time-­‑op'mal ¡control: ¡improved ¡robustness ¡to ¡AWN ¡ −2 10 1.5 CNOT QFT QFT’ SWAP 1 random D (normalized) identity −4 10 log 10 ( T* ) 0.5 0 −6 10 W = SWAP, φ = 0 −0.5 W = QFT, φ = 0 W = QFT, φ = π /2 W = SWAP, φ = π /2 −1 −8 10 −1 −0.5 0 0.5 1 1.5 2 2.5 4 6 8 10 12 log 10 (coupling strength J (1,2) ) T Time-fidelity Pareto fronts for different The critical value T * of the control time target gates (SWAP and QFT) and also depends on the Hamiltonian of the different values of the gate’s global controlled system, in particular, on the phase (0 and π /2 ). strength of inter-qubit couplings. The critical value T * of the control time In the weak coupling regime, the dependence of the critical time on the depends both on the target gate and on coupling strength is its global phase. K. W. Moore et al., arXiv:1112.0333 (2011)

Fluence-­‑op'mal ¡control: ¡improved ¡robustness ¡to ¡MWN ¡ Optimizing robustness to multiplicative white noise in controls is equivalent to minimizing the fluence of the control field. Exploring the dependence of the optimal-field fluence on the initial-field fluence: Average fluence of 25 optimal control fields obtained by starting from random initial fields with the same fluence (the target gate is the Hadamard transform for a one-qubit system): 3 10 T = 1 To enhance the robustness to T = 3 MWN, one should use initial T = 4 2 fields with very small fluence. 10 T = 6 Average optimal field fluence T = 15 T = 20 1 10 0 10 − 1 10 − 6 − 5 − 4 − 3 − 2 − 1 0 1 2 3 10 10 10 10 10 10 10 10 10 10 Initial field fluence

Fluence-­‑op'mal ¡control: ¡improved ¡robustness ¡to ¡MWN ¡ The minimum value of the optimal-field fluence decreases with T , but not monotonically. Exploring the dependence of the minimum optimal-field fluence on control time T : Target: Hadamard transform Target: Pauli Y -gate 16 10 9 14 Fluence of the optimal control field Fluence of the optimal control field 8 12 7 10 6 8 5 4 6 3 4 2 2 1 0 0 2 4 6 8 10 12 14 2 4 6 8 10 12 14 Control time Control time The envelope function decays as 1/ T The envelope function decays as 1/ T The amplitude of oscillations is constant The amplitude of oscillations decays as exp(- γ T)

Robustness ¡to ¡colored ¡noise ¡ For small colored noise, the expected gate error (i.e., the expected fidelity decrease) is determined by an overlap integral involving the Hessian of the optimal control field and the noise autocorrelation function. We define this overlap as the robustness metric K . In particular, for additive colored noise , the robustness metric is For a wide-sense-stationary noise process: The Wiener-Khinchin theorem relates autocorrelation function and power spectral density : Thus, the robustness metric can be expressed as an overlap in the frequency domain: The goal is to minimize the expected gate error by searching for optimal controls with the Hessian “orthogonal” to the control noise (i.e., using the null space of the Hessian).