This article is aimed at extending the framework of optimal control techniques to the situation where the control field values are restricted to a finite set. We propose generalizations of the standard GRAPE algorithm suited to this constraint. We test the validity and the efficiency of this approach for the inversion of an inhomogeneous ensemble of spin systems with different offset frequencies. It is shown that a remarkable efficiency can be achieved even for a very limited number of discrete values. Some applications in nuclear magnetic resonance are discussed.
This article is aimed at extending the framework of optimal control techniques to the situation where the control field values are restricted to a finite set. We propose generalizations of the standard GRAPE algorithm suited to this constraint. We test the validity and the efficiency of this approach for the inversion of an inhomogeneous ensemble of spin systems with different offset frequencies. It is shown that a remarkable efficiency can be achieved even for a very limited number of discrete values. Some applications in nuclear magnetic resonance are discussed.
Zero-field nuclear magnetic resonance (NMR) has been recently developed as a complementary tool for nondestructive structural investigations of matter. However, precise and efficient control in such systems is still at the beginning. Based on the controllability of spin systems, we theoretically study the issue of quantum control in zero-field NMR by using a numerical optimization method based on the gradient ascent pulse engineering (GRAPE) algorithm. A set of quantum gates, including single-qubit and two-qubit gates, are achieved as examples. In particular, we show that homonuclear spins can be individually manipulated in zero-field NMR by the GRAPE method, which extends the ability to manipulate spin systems at zero magnetic field. Quantum control realized by the numerical GRAPE method features the properties of high fidelity, robustness, and hardware-friendly characteristics. High-fidelity quantum control in zero-field NMR systems might provide promising applications in material science, chemical analysis, quantum information processing, and fundamental physics.
We analyze the simultaneous time-optimal control of two-spin systems. The two noncoupled spins, which differ in the value of their chemical offsets, are controlled by the same magnetic fields. Using an appropriate rotating frame, we restrict the study to the case of opposite shifts. We then show that the optimal solution of the inversion problem in a rotating frame is composed of a pulse sequence of maximum intensity and is similar to the optimal solution for inverting only one spin by using a nonresonant control field in the laboratory frame. An example is implemented experimentally using nuclear magnetic resonance techniques.
Quantum optimal control, a powerful toolbox for engineering an optimal control field modulation that most precisely implements a desired quantum operation in the best way possible, has evolved into one of the cornerstones for enabling quantum technologies. The gradient ascent pulse engineering (GRAPE) algorithm is a widely used method in quantum optimal control, which has achieved great success in different physical platforms. However, its computational complexity increases exponentially with the number of qubits, making it challenging to be implemented for large-scale quantum systems. To mitigate this issue, we present the iterative GRAPE algorithm (iGRAPE), which reduces the optimization problem into a series of lower-dimensional subproblems by incorporating disentanglement operations. Our numerical simulations on physical platforms such as nuclear magnetic resonance and superconducting quantum systems demonstrate that iGRAPE significantly enhances state preparation speed. Specifically, compared to GRAPE, iGRAPE achieves up to a five-fold acceleration in preparing Greenberger--Horne--Zeilinger states using a 12-qubit implementation, and up to a 13-fold acceleration for arbitrary state preparation with eight qubits. To further validate our findings, we conduct experimental validation of iGRAPE on a four-qubit nuclear magnetic resonance system. Overall, iGRAPE offers an efficient solution for implementing optimal control in large-scale quantum systems, holding great potential for advancing quantum technologies during the noisy intermediate-scale quantum era.
Quantum control aims to manipulate quantum systems toward specific quantum states or desired operations. Designing highly accurate and effective control steps is vitally important to various quantum applications, including energy minimization and circuit compilation. In this paper we focus on discrete binary quantum control problems and apply different optimization algorithms and techniques to improve computational efficiency and solution quality. Specifically, we develop a generic model and extend it in several ways. We introduce a squared L2 -penalty function to handle additional side constraints, to model requirements such as allowing at most one control to be active. We introduce a total variation (TV) regularizer to reduce the number of switches in the control. We modify the popular gradient ascent pulse engineering (GRAPE) algorithm, develop a new alternating direction method of multipliers (ADMM) algorithm to solve the continuous relaxation of the penalized model, and then apply rounding techniques to obtain binary control solutions. We propose a modified trust-region method to further improve the solutions. Our algorithms can obtain high-quality control results, as demonstrated by numerical studies on diverse quantum control examples.
Quadratic convergence throughout the active space is achieved for the gradient ascent pulse engineering (GRAPE) family of quantum optimal control algorithms. We demonstrate in this communication that the Hessian of the GRAPE fidelity functional is unusually cheap, having the same asymptotic complexity scaling as the functional itself. This leads to the possibility of using very efficient numerical optimization techniques. In particular, the Newton-Raphson method with a rational function optimization (RFO) regularized Hessian is shown in this work to require fewer system trajectory evaluations than any other algorithm in the GRAPE family. This communication describes algebraic and numerical implementation aspects (matrix exponential recycling, Hessian regularization, etc.) for the RFO Newton-Raphson version of GRAPE and reports benchmarks for common spin state control problems in magnetic resonance spectroscopy.
We analyze the saturation of a nuclear magnetic resonance (NMR) signal using optimal magnetic fields. We consider both the problems of minimizing the duration of the control and its energy for a fixed duration. We solve the optimal control problems by using geometric methods and a purely numerical approach, the grape algorithm, the two methods being based on the application of the Pontryagin maximum principle. A very good agreement is obtained between the two results. The optimal solutions for the energy-minimization problem are finally implemented experimentally with available NMR techniques.