Second order gradient ascent pulse engineering

Controlling NMR spin systems for quantum computation

Nuclear magnetic resonance is arguably both the best available quantum technology for implementing simple quantum computing experiments and the worst technology for building large scale quantum computers that has ever been seriously put forward. After a few years of rapid growth, leading to an implementation of Shor's quantum factoring algorithm in a seven-spin system, the field started to reach its natural limits and further progress became challenging. Rather than pursuing more complex algorithms on larger systems, interest has now largely moved into developing techniques for the precise and efficient manipulation of spin states with the aim of developing methods that can be applied in other more scalable technologies and within conventional NMR. However, the user friendliness of NMR implementations means that they remain popular for proof-of-principle demonstrations of simple quantum information protocols.

Optimal control pulses for the 1.2-GHz (28.2-T) NMR spectrometers

The ability to measure nuclear magnetic resonance (NMR) spectra with a large sample volume is crucial for concentration-limited biological samples to attain adequate signal-to-noise (S/N) ratio. The possibility to measure with a 5-mm cryoprobe is currently absent at the 1.2-GHz NMR instruments due to the exceedingly high radio frequency power demands, which is four times compared to 600-MHz instruments. Here, we overcome the high-power demands by designing optimal control (OC) pulses with up to 20 times lower power requirements than currently necessary at a 1.2-GHz spectrometer. We show that multidimensional biomolecular NMR experiments constructed using these OC pulses can bestow improvement in the S/N ratio of up to 26%. With the expected power limitations of a 5-mm cryoprobe, we observe an enhancement in the S/N ratio of more than 240% using our OC sequences. This motivates the development of a cryoprobe with a larger volume than the current 3-mm cryoprobes.

Simulation and design of shaped pulses beyond the piecewise-constant approximation

Response functions of resonant circuits create ringing artefacts if their input changes rapidly. When physical limits of electromagnetic spectroscopies are explored, this creates two types of problems. Firstly, simulation: the system must be propagated accurately through every response transient, this may be computationally expensive. Secondly, optimal control: circuit response must be taken into account; it may be advantageous to design pulses that are resilient to such distortions. At the root of both problems is the popular piecewise-constant approximation for control sequences in the rotating frame; in magnetic resonance it has persisted since the earliest days and has become entrenched in the commercially available hardware. In this paper, we report an implementation and benchmarks of recent Lie-group methods that can efficiently simulate and optimise smooth control sequences.

Adaptive optimal control of entangled qubits

Developing fast, robust, and accurate methods for optimal control of quantum systems comprising interacting particles is one of the most active areas of current science. Although a valuable repository of algorithms is available for numerical applications in quantum control, the high computational cost is somewhat overlooked. Here, we present a fast and accurate optimal control algorithm for systems of interacting qubits, QOALA (quantum optimal control by adaptive low-cost algorithm), which is predicted to offer [Formula: see text](M2) speedup for an M-qubit system, compared to the state-of-the-art exact methods, without compromising overall accuracy of the optimal solution. The method is general and compatible with diverse Hamiltonian structures. The proposed approach uses inexpensive low-accuracy approximations of propagators far from the optimum, adaptively switching to higher accuracy, higher-cost propagators when approaching the optimum. In addition, the utilization of analytical Lie algebraic derivatives that do not require computationally expensive matrix exponential brings even better performance.

Band-selective universal 90° and 180° rotation pulses covering the aliphatic carbon chemical shift range for triple resonance experiments on 1.2 GHz spectrometers

Biomolecular NMR spectroscopy requires large magnetic field strengths for high spectral resolution. Today's highest fields comprise proton Larmor frequencies of 1.2 GHz and even larger field strengths are to be expected in the future. In protein triple resonance experiments, various carbon bandwidths need to be excited by selective pulses including the large aliphatic chemical shift range. When the spectrometer field strength is increased, the length of these pulses has to be decreased by the same factor, resulting in higher rf-amplitudes being necessary in order to cover the required frequency region. Currently available band-selective pulses like Q3/Q5 excite a narrow bandwidth compared to the necessary rf-amplitude. Because the maximum rf-power allowed in probeheads is limited, none of the selective universal rotation pulses reported so far is able to cover the full [Formula: see text]C aliphatic region on 1.2 GHz spectrometers. In this work, we present band-selective 90° and 180° universal rotation pulses (SURBOP90 and SURBOP180) that have a higher ratio of selective bandwidth to maximum rf-amplitude than standard pulses. Simulations show that these pulses perform better than standard pulses, e. g. Q3/Q5, especially when rf-inhomogeneity is taken into account. The theoretical and experimental performance is demonstrated in offset profiles and by implementing the SURBOP pulses in an HNCACB experiment at 1.2 GHz.

Quantum-enhanced radiometry via approximate quantum error correction

Quantum sensing based on exotic quantum states is appealing for practical metrology applications and fundamental studies. However, these quantum states are vulnerable to noise and the resulting quantum enhancement is weakened in practice. Here, we experimentally demonstrate a quantum-enhanced sensing scheme with a bosonic probe, by exploring the large Hilbert space of the bosonic mode and developing both the approximate quantum error correction and the quantum jump tracking approaches. In a practical radiometry scenario, we attain a 5.3 dB enhancement of sensitivity, which reaches 9.1 × 10⁻⁴ Hz^−1/2 when measuring the excitation population of a receiver mode. Our results demonstrate the potential of quantum sensing with near-term quantum technologies, not only shedding new light on the quantum advantage of sensing, but also stimulating further efforts on bosonic quantum technologies.

Exotic quantum states can be advantageous for sensing, but are very fragile, so that some form of quantum error correction is needed. Here, the authors show how approximate QEC helps overcoming decoherence due to noise when measuring the excitation population of a receiver mode in a superconducting circuit.

SORDOR pulses: expansion of the Böhlen-Bodenhausen scheme for low-power broadband magnetic resonance

A novel type of efficient broadband pulse, called second-order phase dispersion by optimised rotation (SORDOR), has recently been introduced. In contrast to adiabatic excitation, SORDOR-90 pulses provide effective transverse 90∘ rotations throughout their bandwidth, with a quadratic offset dependence of the phase in the x , y plane. Together with phase-matched SORDOR-180 pulses, this enables the Böhlen-Bodenhausen broadband refocusing approach for linearly frequency-swept pulses to be extended to any type of 90∘ /180∘ pulse-delay sequence. Example pulse shapes are characterised in theory and experiment, and an example application is given with a 19 F -PROJECT experiment for measuring relaxation times with reduced distortions due to J -coupling evolution.

Dynamical Uncertainty Propagation with Noisy Quantum Parameters

Many quantum technologies rely on high-precision dynamics, which raises the question of how these are influenced by the experimental uncertainties that are always present in real-life settings. A standard approach in the literature to assess this is Monte Carlo sampling, which suffers from two major drawbacks. First, it is computationally expensive. Second, it does not reveal the effect that each individual uncertainty parameter has on the state of the system. In this Letter, we evade both these drawbacks by incorporating propagation of uncertainty directly into simulations of quantum dynamics, thereby obtaining a method that is orders of magnitude faster than Monte Carlo simulations and directly provides information on how each uncertainty parameter influences the system dynamics. Additionally, we compare our method to experimental results obtained using the IBM quantum computers.

Concurrent J-evolving refocusing pulses

Conventional refocusing pulses are optimised for a single spin without considering any type of coupling. However, despite the fact that most couplings will result in undesired distortions, refocusing in delay-pulse-delay-type sequences with desired heteronuclear coherence transfer might be enhanced considerably by including coupling evolution into pulse design. We provide a proof of principle study for a Hydrogen-Carbon refocusing pulse sandwich with inherent J-evolution following the previously reported ICEBERG-principle with improved performance in terms of refocusing performance and/or overall effective coherence transfer time. Pulses are optimised using optimal control theory with a newly derived quality factor and z-controls as an efficient tool to speed up calculations. Pulses are characterised in theory and experiment and compared to conventional concurrent refocusing pulses, clearly showing an improvement for the J-evolving pulse sandwich. As a side-product, also efficient J-compensated resfocusing pulse sandwiches - termed BUBU pulses following the nomenclature of previous J-compensated BUBI and BEBE^tr pulse sandwiches - have been optimised.

Optimal control gradient precision trade-offs: Application to fast generation of DeepControl libraries for MRI

We have recently demonstrated supervised deep learning methods for rapid generation of radiofrequency pulses in magnetic resonance imaging (https://doi.org/10.1002/mrm.27740, https://doi.org/10.1002/mrm.28667). Unlike the previous iterative optimization approaches, deep learning methods generate a pulse using a fixed number of floating-point operations - this is important in MRI, where patient-specific pulses preferably must be produced in real time. However, deep learning requires vast training libraries, which must be generated using the traditional methods, e.g., iterative quantum optimal control methods. Those methods are usually variations of gradient descent, and the calculation of the gradient of the performance metric with respect to the pulse waveform can be the most numerically intensive step. In this communication, we explore various ways in which the calculation of gradients in quantum optimal control theory may be accelerated. Four optimization avenues are explored: truncated commutator series expansions at zeroth and first order, a novel midpoint truncation scheme at first order, and the exact complex-step method. For the spin systems relevant to MRI, the first-order midpoint truncation is found to be sufficiently accurate, but also significantly faster than the machine precision gradient. This makes the generation of training databases for the machine learning methods considerably more realistic.

Maximizing efficiency of dipolar recoupling in solid-state NMR using optimal control sequences

Dipolar recoupling is a central concept in the nuclear magnetic resonance spectroscopy of powdered solids and is used to establish correlations between different nuclei by magnetization transfer. The efficiency of conventional cross-polarization methods is low because of the inherent radio frequency (rf) field inhomogeneity present in the magic angle spinning (MAS) experiments and the large chemical shift anisotropies at high magnetic fields. Very high transfer efficiencies can be obtained using optimal control–derived experiments. These sequences had to be optimized individually for a particular MAS frequency. We show that by adjusting the length and the rf field amplitude of the shaped pulse synchronously with sample rotation, optimal control sequences can be successfully applied over a range of MAS frequencies without the need of reoptimization. This feature greatly enhances their applicability on spectrometers operating at differing external fields where the MAS frequency needs to be adjusted to avoid detrimental resonance effects.

Quantum control of bosonic modes with superconducting circuits

Bosonic modes have wide applications in various quantum technologies, such as optical photons for quantum communication, magnons in spin ensembles for quantum information storage and mechanical modes for reversible microwave-to-optical quantum transduction. There is emerging interest in utilizing bosonic modes for quantum information processing, with circuit quantum electrodynamics (circuit QED) as one of the leading architectures. Quantum information can be encoded into subspaces of a bosonic superconducting cavity mode with long coherence time. However, standard Gaussian operations (e.g., beam splitting and two-mode squeezing) are insufficient for universal quantum computing. The major challenge is to introduce additional nonlinear control beyond Gaussian operations without adding significant bosonic loss or decoherence. Here we review recent advances in universal control of a single bosonic code with superconducting circuits, including unitary control, quantum feedback control, driven-dissipative control and holonomic dissipative control. Various approaches to entangling different bosonic modes are also discussed.

Conservation-Law-Based Global Bounds to Quantum Optimal Control

Active control of quantum systems enables diverse applications ranging from quantum computation to manipulation of molecular processes. Maximum speeds and related bounds have been identified from uncertainty principles and related inequalities, but such bounds utilize only coarse system information and loosen significantly in the presence of constraints and complex interaction dynamics. We show that an integral-equation-based formulation of conservation laws in quantum dynamics leads to a systematic framework for identifying fundamental limits to any quantum control scenario. We demonstrate the utility of our bounds in three scenarios-three-level driving, decoherence suppression, and maximum-fidelity gate implementations-and show that in each case our bounds are tight or nearly so. Global bounds complement local-optimization-based designs, illuminating performance levels that may be possible, as well as those that cannot be surpassed.

High-Efficiency Arbitrary Quantum Operation on a High-Dimensional Quantum System

The ability to manipulate quantum systems lies at the heart of the development of quantum technology. The ultimate goal of quantum control is to realize arbitrary quantum operations (AQUOs) for all possible open quantum system dynamics. However, the demanding extra physical resources impose great obstacles. Here, we experimentally demonstrate a universal approach of AQUO on a photonic qudit with the minimum physical resource of a two-level ancilla and a log_{2}d-scale circuit depth for a d-dimensional system. The AQUO is then applied in a quantum trajectory simulation for quantum subspace stabilization and quantum Zeno dynamics, as well as incoherent manipulation and generalized measurements of the qudit. Therefore, the demonstrated AQUO for complete quantum control would play an indispensable role in quantum information science.

Experimental Simulation of Open Quantum System Dynamics via Trotterization

Digital quantum simulators provide a diversified tool for solving the evolution of quantum systems with complicated Hamiltonians and hold great potential for a wide range of applications. Although much attention is paid to the unitary evolution of closed quantum systems, dissipation and noise are vital in understanding the dynamics of practical quantum systems. In this work, we experimentally demonstrate a digital simulation of an open quantum system in a controllable Markovian environment with the assistance of a single ancillary qubit. By Trotterizing the quantum Liouvillians, the continuous evolution of an open quantum system is effectively realized, and its application in error mitigation is demonstrated by adjusting the simulated noise intensities. High-order Trotter for open quantum dynamics is also experimentally investigated and shows higher accuracy. Our results represent a significant step toward hardware-efficient simulation of open quantum systems and error mitigation in quantum algorithms in noisy intermediate-scale quantum systems.

Optimization of phase dispersion enables broadband excitation without homonuclear coupling artifacts

In NMR spectroscopy, many specialized shaped pulses are available for broadband excitation, beyond the bandwidth of conventional high-powered hard pulses. These shaped pulses typically have long duration. However, long-duration pulses are unsuitable for spectra containing significant homonuclear couplings, such as polyfluorinated compounds in 19F NMR. J-coupling evolution during the excitation pulse leads to spectral artifacts and incorrect peak integrals. Here, we report an approach to optimal control pulse design which significantly reduces the pulse length required to excite large bandwidths of chemical shift frequencies. The target state phase is not chosen beforehand but is instead only constrained to be linearly dependent on offset frequency. The first-order phase of the target state is then treated as a free-variable, to be optimized at the same time as the RF waveform itself. The resulting spectra are easily phased using standard NMR processing software. We observe that the required pulse length is significantly shorter than for currently available in-phase excitation schemes. Spectral artifacts from homonuclear couplings are avoided. We also demonstrate that pure in-phase excitation can be obtained over the same bandwidth by appending two inversion pulses, at the expense of increased overall duration.

High-Fidelity Software-Defined Quantum Logic on a Superconducting Qudit

We present an efficient approach to achieving arbitrary, high-fidelity control of a multilevel quantum system using optimal control techniques. As an demonstration, we implement a continuous, software-defined microwave pulse to realize a 0↔2 SWAP gate that achieves an average gate fidelity of 99.4%. We describe our procedure for extracting the system Hamiltonian, calibrating the quantum and classical hardware chain, and evaluating the gate fidelity. Our work represents an alternative, fully generalizable route towards achieving universal quantum control by leveraging optimal control techniques.

Universal nonselective excitation and refocusing pulses with improved robustness to off-resonance for Magnetic Resonance Imaging at 7 Tesla with parallel transmission

PURPOSE

In MRI at ultra-high field, the k T -point and spiral nonselective (SPINS) pulse design techniques can be advantageously combined with the parallel transmission (pTX) and universal pulse techniques to create uniform excitation in a calibration-free manner. However, in these approaches, pulse duration is typically increased as compared to standard hard pulses, and excitation quality in regions exhibiting large resonance frequency offsets often suffer. This limitation is inherent to structure of k T -point or SPINS pulse, and likely can be mitigated using parameterization-free pulse design approaches.

METHODS

The Gradient Ascent Pulse Engineering (GRAPE) algorithm was used to design parameterization-free RF and magnetic field gradient (MFG) waveforms for creating 8 ∘ excitation, up to 105 ∘ scalable refocusing and inversion, nonselectively across the brain. Simulations were performed to provide flip angle normalized root-mean-squares error (FA-NRMSE) estimations for the 8 ∘ and the 180 ∘ k T -point, SPINS, and GRAPE pulses. GRAPE pulses were tested experimentally with anatomical head scans at 7T.

RESULTS

As compared to k T -points and SPINS, GRAPE provided substantial improvement of excitation, refocusing, and inversion quality at off-resonance while at least preserving the same global FA-NRMSE performance. As compared to k T -points, GRAPE allowed for a substantial reduction of the pulse duration for the 8 ∘ excitation and the 105 ∘ refocusing.

CONCLUSIONS

Parameterization-free universal nonselective pTX-pulses were successfully computed using GRAPE. Performance gains as compared to k T -points were validated numerically and experimentally for three imaging protocols. In its current implementation, the computational burden of GRAPE limits its use to applications where pulse computations are not subject to time constraints.

Large-scale magnetic resonance simulations: A tutorial

Computational modeling is becoming an essential tool in magnetic resonance to design and optimize experiments, test the performance of theoretical models, and interpret experimental data. Recent theoretical research and software development made possible simulations of large spin systems, for example, proteins with thousands of spins, in reasonable time. In the last few years, the Fokker-Planck formalism also re-emerged due to its ability to handle spatial dynamics. The purpose of this tutorial is to describe advantages and disadvantages of the most common formalisms, the latest developments and strategies to improve the computational efficiency, and to guide users in the setting up of a simulation using the Spinach software.

Demonstration of Controlled-Phase Gates between Two Error-Correctable Photonic Qubits

To realize fault-tolerant quantum computing, it is necessary to store quantum information in logical qubits with error correction functions, realized by distributing a logical state among multiple physical qubits or by encoding it in the Hilbert space of a high-dimensional system. Quantum gate operations between these error-correctable logical qubits, which are essential for implementation of any practical quantum computational task, have not been experimentally demonstrated yet. Here we demonstrate a geometric method for realizing controlled-phase gates between two logical qubits encoded in photonic fields stored in cavities. The gates are realized by dispersively coupling an ancillary superconducting qubit to these cavities and driving it to make a cyclic evolution depending on the joint photonic state of the cavities, which produces a conditional geometric phase. We first realize phase gates for photonic qubits with the logical basis states encoded in two quasiorthogonal coherent states, which have important implications for continuous-variable-based quantum computation. Then we use this geometric method to implement a controlled-phase gate between two binomially encoded logical qubits, which have an error-correctable function.

Efficient Quantum Error Correction of Dephasing Induced by a Common Fluctuator

Quantum error correction is expected to be essential in large-scale quantum technologies. However, the substantial overhead of qubits it requires is thought to greatly limit its utility in smaller, near-term devices. Here we introduce a new family of special-purpose quantum error-correcting codes that offer an exponential reduction in overhead compared to the usual repetition code. They are tailored for a common and important source of decoherence in current experiments, whereby a register of qubits is subject to phase noise through coupling to a common fluctuator, such as a resonator or a spin defect. The smallest instance encodes one logical qubit into two physical qubits, and corrects decoherence to leading-order using a constant number of one- and two-qubit operations. More generally, while the repetition code on n qubits corrects errors to order t^{O(n)}, with t the time between recoveries, our codes correct to order t^{O(2^{n})}. Moreover, they are robust to model imperfections in small- and intermediate-scale devices, where they already provide substantial gains in error suppression. As a result, these hardware-efficient codes open a potential avenue for useful quantum error correction in near-term, pre-fault tolerant devices.

Quantum Adiabatic Doping with Incommensurate Optical Lattices

Quantum simulations of Fermi-Hubbard models have been attracting considerable effort in the optical lattice research, with the ultracold antiferromagnetic atomic phase reached at half filling in recent years. An unresolved issue is to dope the system while maintaining the low thermal entropy. Here we propose to achieve the low temperature phase of the doped Fermi-Hubbard model using incommensurate optical lattices through adiabatic quantum evolution. In this theoretical proposal, we find that one major problem about the adiabatic doping is atomic localization in the incommensurate lattice, potentially causing an exponential slowing down of the adiabatic procedure. We study both one- and two-dimensional incommensurate optical lattices, and find that the localization prevents efficient adiabatic doping in the strong lattice regime for both cases. With density matrix renormalization group calculation, we further show that the slowing down problem in one dimension can be circumvented by considering interaction induced many-body delocalization, which is experimentally feasible using Feshbach resonance techniques. This protocol is expected to be efficient as well in two dimensions where the localization phenomenon is less stable.

On the use of frequency-swept pulses and pulses designed with optimal control theory for the acquisition of ultra-wideline NMR spectra

Frequency-swept (FS) pulses, such as wideband uniform-rate smooth-truncation (WURST) pulses, have found much success for the acquisition of ultra-wideline (UW) solid-state NMR spectra. In this preliminary study, new pulses and pulse sequences are explored in simulation and experimentally for several nuclei exhibiting UWNMR powder patterns under static conditions, including ¹¹⁹Sn (I = 1/2), ¹⁹⁵Pt (I = 1/2), ²H (I = 1), and ⁷¹Ga (I = 3/2). First, hyperbolic secant (HS) and tanh/tan (THT) pulses are tested and implemented as excitation and refocusing pulses in spin-echo and Carr-Purcell/Meiboom Gill (CPMG)-type sequences, and shown to have comparable performances to analogous WURST pulses. Second, optimal control theory (OCT) is utilized for the design of new Optimal Control Theory Optimized Broadband Excitation and Refocusing (OCTOBER) pulses, using carefully parameterized WURST, THT, and HS pulses as starting points. Some of the new OCTOBER pulses used in spin-echo sequences are capable of efficient broadband excitation and refocusing, in some cases resulting in spectra with increased signal enhancements over those obtained in experiments using conventional FS pulses. Finally, careful consideration of the spin dynamics of several systems, by monitoring of the time evolution of the density matrix via the Liouville-von Neumann equation and analysis of the time-resolved Fourier transforms of the pulses, lends insight into the underlying mechanisms of the FS and OCTOBER pulses. This is crucial for understanding their performance in terms of generating uniformly excited patterns of high signal intensity, and for identifying trends that may offer pathways to generalized parameterization and/or new pulse shapes.

Witnessing Quantum Resource Conversion within Deterministic Quantum Computation Using One Pure Superconducting Qubit

Deterministic quantum computation with one qubit (DQC1) is iconic in highlighting that exponential quantum speedup may be achieved with negligible entanglement. Its discovery catalyzed a heated study of general quantum resources, and various conjectures regarding their role in DQC1's performance advantage. Coherence and discord are prominent candidates, respectively, characterizing nonclassicality within localized and correlated systems. Here we realize DQC1 within a superconducting system, engineered such that the dynamics of coherence and discord can be tracked throughout its execution. We experimentally confirm that DQC1 acts as a resource converter, consuming coherence to generate discord during its operation. Our results highlight superconducting circuits as a promising platform for both realizing DQC1 and related algorithms, and experimentally characterizing resource dynamics within quantum protocols.

Flexible and efficient optimization of quantitative sequences using automatic differentiation of Bloch simulations

PURPOSE

To investigate a computationally efficient method for optimizing the Cramér-Rao Lower Bound (CRLB) of quantitative sequences without using approximations or an analytical expression of the signal.

METHODS

Automatic differentiation was applied to Bloch simulations and used to optimize several quantitative sequences without the need for approximations or an analytical expression. The results were validated with in vivo measurements and comparisons to prior art. Multi-echo spin echo and DESPO T 1 were used as benchmarks to verify the CRLB implementation. The CRLB of the Magnetic Resonance Fingerprinting (MRF) sequence, which has a complicated analytical formulation, was also optimized using automatic differentiation.

RESULTS

The sequence parameters obtained for multi-echo spin echo and DESPO T 1 matched results obtained using conventional methods. In vivo, MRF scans demonstrate that the CRLB optimization obtained with automatic differentiation can improve performance in presence of white noise. For MRF, the CRLB optimization converges in 1.1 CPU hours for N TR = 400 and has O ( N TR ) asymptotic runtime scaling for the calculation of the CRLB objective and gradient.

CONCLUSIONS

Automatic differentiation can be used to optimize the CRLB of quantitative sequences without using approximations or analytical expressions. For MRF, the runtime is computationally efficient and can be used to investigate confounding factors as well as MRF sequences with a greater number of repetitions.

Heisenberg-limited single-mode quantum metrology in a superconducting circuit

Two-mode interferometers lay the foundations for quantum metrology. Instead of exploring quantum entanglement in the two-mode interferometers, a single bosonic mode also promises a measurement precision beyond the shot-noise limit (SNL) by taking advantage of the infinite-dimensional Hilbert space of Fock states. Here, we demonstrate a single-mode phase estimation that approaches the Heisenberg limit (HL) unconditionally. Due to the strong dispersive nonlinearity and long coherence time of a microwave cavity, quantum states of the form [Formula: see text] can be generated, manipulated and detected with high fidelities, leading to an experimental phase estimation precision scaling as ∼N^-0.94. A 9.1 dB enhancement of the precision over the SNL at N = 12 is achieved, which is only 1.7 dB away from the HL. Our experimental architecture is hardware efficient and can be combined with quantum error correction techniques to fight against decoherence, and thus promises quantum-enhanced sensing in practical applications.

Defeating the Matrix

These are personal recollections and musings, written for the 50th Anniversary of the Journal of Magnetic Resonance. They are distilled from twenty years of writing simulation code, and filtered through hindsight and sarcastic intransigence. To me, the biggest recent achievements of the magnetic resonance community in the field of theory and simulation have been the successful war on the exponential scaling, the powerful and general simulation software, and the return to elegant notation. It appears that our future will be defined by computers. Three aspects are pertinent: simulation as the experiment is designed, optimal control as the experiment proceeds, and machine learning at the data processing stage.

Hybrid-state free precession in nuclear magnetic resonance

The dynamics of large spin-1/2 ensembles are commonly described by the Bloch equation, which is characterized by the magnetization's non-linear response to the driving magnetic field. Consequently, most magnetic field variations result in non-intuitive spin dynamics, which are sensitive to small calibration errors. Although simplistic field variations result in robust spin dynamics, they do not explore the richness of the system's phase space. Here, we identify adiabaticity conditions that span a large experiment design space with tractable dynamics. All dynamics are trapped in a one-dimensional subspace, namely in the magnetization's absolute value, which is in a transient state, while its direction adiabatically follows the steady state. In this hybrid state, the polar angle is the effective drive of the spin dynamics. As an example, we optimize this drive for robust and efficient quantification of spin relaxation times and utilize it for magnetic resonance imaging of the human brain.

Optimized quantification of spin relaxation times in the hybrid state

PURPOSE

The optimization and analysis of spin ensemble trajectories in the hybrid state-a state in which the direction of the magnetization adiabatically follows the steady state while the magnitude remains in a transient state.

METHODS

Numerical optimizations were performed to find spin ensemble trajectories that minimize the Cramér-Rao bound for T 1 -encoding, T 2 -encoding, and their weighted sum, respectively, followed by a comparison between the Cramér-Rao bounds obtained with our optimized spin-trajectories, Look-Locker sequences, and multi-spin-echo methods. Finally, we experimentally tested our optimized spin trajectories with in vivo scans of the human brain.

RESULTS

After a nonrecurring inversion segment on the southern half of the Bloch sphere, all optimized spin trajectories pursue repetitive loops on the northern hemisphere in which the beginning of the first and the end of the last loop deviate from the others. The numerical results obtained in this work align well with intuitive insights gleaned directly from the governing equation. Our results suggest that hybrid-state sequences outperform traditional methods. Moreover, hybrid-state sequences that balance T 1 - and T 2 -encoding still result in near optimal signal-to-noise efficiency for each relaxation time. Thus, the second parameter can be encoded at virtually no extra cost.

CONCLUSIONS

We provided new insights into the optimal encoding processes of spin relaxation times in order to guide the design of robust and efficient pulse sequences. We found that joint acquisitions of T 1 and T 2 in the hybrid state are substantially more efficient than sequential encoding techniques.

Quantum generative adversarial learning in a superconducting quantum circuit

Generative adversarial learning is one of the most exciting recent breakthroughs in machine learning. It has shown splendid performance in a variety of challenging tasks such as image and video generation. More recently, a quantum version of generative adversarial learning has been theoretically proposed and shown to have the potential of exhibiting an exponential advantage over its classical counterpart. Here, we report the first proof-of-principle experimental demonstration of quantum generative adversarial learning in a superconducting quantum circuit. We demonstrate that, after several rounds of adversarial learning, a quantum-state generator can be trained to replicate the statistics of the quantum data output from a quantum channel simulator, with a high fidelity (98.8% on average) so that the discriminator cannot distinguish between the true and the generated data. Our results pave the way for experimentally exploring the intriguing long-sought-after quantum advantages in machine learning tasks with noisy intermediate-scale quantum devices.

Experimental repetitive quantum channel simulation

Universal control of quantum systems is a major goal to be achieved for quantum information processing, which demands thorough understanding of fundamental quantum mechanics and promises applications of quantum technologies. So far, most studies concentrate on ideally isolated quantum systems governed by unitary evolutions, while practical quantum systems are open and described by quantum channels due to their inevitable coupling to environment. Here, we experimentally simulate arbitrary quantum channels for an open quantum system, i.e. a single photonic qubit in a superconducting quantum circuit. The arbitrary channel simulation is achieved with minimum resource of only one ancilla qubit and measurement-based adaptive control. By repetitively implementing the quantum channel simulation, we realize an arbitrary Liouvillian for a continuous evolution of an open quantum system for the first time. Our experiment provides not only a testbed for understanding quantum noise and decoherence, but also a powerful tool for full control of practical open quantum systems.

Feedback control optimisation of ESR experiments

Numerically optimised microwave pulses are used to increase excitation efficiency and modulation depth in electron spin resonance experiments performed on a spectrometer equipped with an arbitrary waveform generator. The optimisation procedure is sample-specific and reminiscent of the magnet shimming process used in the early days of nuclear magnetic resonance - an objective function (for example, echo integral in a spin echo experiment) is defined and optimised numerically as a function of the pulse waveform vector using noise-resilient gradient-free methods. We found that the resulting shaped microwave pulses achieve higher excitation bandwidth and better echo modulation depth than the pulse shapes used as the initial guess. Although the method is theoretically less sophisticated than simulation based quantum optimal control techniques, it has the advantage of being free of the linear response approximation; rapid electron spin relaxation also means that the optimisation takes only a few seconds. This makes the procedure fast, convenient, and easy to use. An important application of this method is at the final stage of the implementation of theoretically designed pulse shapes: compensation of pulse distortions introduced by the instrument. The performance is illustrated using spin echo and out-of-phase electron spin echo envelope modulation experiments. Interface code between Bruker SpinJet arbitrary waveform generator and Matlab is included in versions 2.2 and later of the Spinach library.

A simplified framework to optimize MRI contrast preparation

PURPOSE

This article proposes a rigorous optimal control framework for the design of preparation schemes that optimize MRI contrast based on relaxation time differences.

METHODS

Compared to previous optimal contrast preparation schemes, a drastic reduction of the optimization parameter number is performed. The preparation scheme is defined as a combination of several block pulses whose flip angles, phase terms and inter-pulse delays are optimized to control the magnetization evolution.

RESULTS

The proposed approach reduces the computation time of <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"><mml:msub><mml:mi>B</mml:mi> <mml:mn>0</mml:mn></mml:msub> </mml:math> -robust preparation schemes to around a minute (whereas several hours were required with previous schemes), with negligible performance loss. The chosen parameterization allows to formulate the total preparation duration as a constraint, which improves the overall compromise between contrast performance and preparation time. Simulation, in vitro and in vivo results validate this improvement, illustrate the straightforward applicability of the proposed approach, and point out its flexibility in terms of achievable contrasts. Major improvement is especially achieved for short-T₂ enhancement, as shown by the acquisition of a non-trivial contrast on a rat brain, where a short-T₂ white matter structure (corpus callosum) is enhanced compared to surrounding gray matter tissues (hippocampus and neocortex).

CONCLUSIONS

This approach proposes key advances for the design of optimal contrast preparation sequences, that emphasize their ability to generate non-standard contrasts, their potential benefit in a clinical context, and their straightforward applicability on any MR system.

Constant gradient elastography with optimal control RF pulses

This article presents a new motion encoding strategy to perform magnetic resonance elastography (MRE). Instead of using standard motion encoding gradients, a tailored RF pulse is designed to simultaneously perform selective excitation and motion encoding in presence of a constant gradient. The RF pulse is designed with a numerical optimal control algorithm, in order to obtain a magnetization phase distribution that depends on the displacement characteristics inside each voxel. As a consequence, no post-excitation encoding gradients are required. This offers numerous advantages, such as reducing eddy current artifacts, and relaxing the constraint on the gradients maximum switch rate. It also allows to perform MRE with ultra-short TE acquisition schemes, which limits T₂ decay and optimizes signal-to-noise ratio. The pulse design strategy is developed and analytically analyzed to clarify the encoding mechanism. Finally, simulations, phantom and ex vivo experiments show that phase-to-noise ratios are improved when compared to standard MRE encoding strategies.

Rotation relaxation splitting for optimizing parallel RF excitation pulses with T1- and T2-relaxations in MRI

Exact solutions of the Bloch equations with T₁- and T₂-relaxation terms for piecewise constant magnetic fields are numerically challenging. We therefore investigate an approximation for the achieved magnetization in which rotations and relaxations are split into separate operations. We develop an estimate for its accuracy and explicit first and second order derivatives with respect to the complex excitation radio frequency voltages. In practice, the deviation between an exact solution of the Bloch equations and this rotation relaxation splitting approximation seems negligible. Its computation times are similar to exact solutions without relaxation terms. We apply the developed theory to numerically optimize radio frequency excitation waveforms with T₁- and T₂-relaxations in several examples.

NMRCloudQ: a quantum cloud experience on a nuclear magnetic resonance quantum computer

Cloud-based quantum computing is anticipated to be the most useful and reachable form for public users to experience with the power of quantum. As initial attempts, IBM Q has launched influential cloud services on a superconducting quantum processor in 2016, but no other platforms has followed up yet. Here, we report our new cloud quantum computing service - NMRCloudQ (http://nmrcloudq.com/zh-hans/), where nuclear magnetic resonance, one of the pioneer platforms with mature techniques in experimental quantum computing, plays as the role of implementing computing tasks. Our service provides a comprehensive software environment preconfigured with a list of quantum information processing packages, and aims to be freely accessible to either amateurs that look forward to keeping pace with this quantum era or professionals that are interested in carrying out real quantum computing experiments in person. In our current version, four qubits are already usable with in average 99.10% single-qubit gate fidelity and 97.15% two-qubit fidelity via randomized benchmaking tests. Improved control precisions as well as a new seven-qubit processor are also in preparation and will be available later.

Cooperative broadband spin echoes through optimal control

The Hahn echo sequence is one of the most common building blocks in magnetic resonance, consisting of an excitation pulse and a refocusing pulse. Conventional approaches to improve the performance of echo experiments focused on the optimization of individual pulses, compensating their own imperfections. Here we present an approach to concurrently design both pulses such that they also compensate each others imperfections. The fact that for such cooperative pulses the individual pulses do not need to be perfect provides additional degrees of freedom, resulting in improved overall Hahn echo performance. Single-scan cooperative pulses are compared to conventional approaches by simulations as well as experiments.

Bang-bang optimal control of large spin systems: Enhancement of 13C-13C singlet-order at natural abundance

Using a bang-bang optimal control technique, we transfer polarization from a set of abundant high-γ nuclei directly to singlet order of a low-γ spin-pair. This approach is analogous to algorithmic cooling (AC) procedure used in quantum state purification. Specifically, we apply this method for enhancing the singlet order in a natural abundant ¹³C- ¹³C spin pair by exploiting nine equivalent protons of an 11-spin system. Compared to the standard method not involving polarization transfer, we find an enhancement of singlet order by about 3.4 times. In addition, since the singlet magnetization is contributed by the faster relaxing protons, the recycle delay is halved. Thus effectively we observe a reduction in the overall experimental time by a factor of 23. We also discuss a possible extension of AC, known as heat-bath algorithmic cooling (HBAC).

An optimal control approach to design entire relaxation dispersion experiments

A general approach is introduced to optimize experiments for the analysis of spin systems in the presence of chemical exchange. Rather than optimizing individual pulse sequence elements, such as refocusing pulses, entire relaxation dispersion sequences are optimized in the form of a single shaped pulse. This is achieved by defining a performance index that is only based on the remaining signal after the relaxation dispersion sequence for a range of exchange, relaxation, offset, and rf inhomogeneity parameters. The approach is demonstrated by optimizing energy-limited broadband relaxation dispersion sequences that closely approach the overall effect of ideal CPMG sequences. As illustrated both theoretically and experimentally, significant improvements are found compared to standard amplitude or energy-limited CPMG sequences.

Active control of the spatial MRI phase distribution with optimal control theory

This paper investigates the use of Optimal Control (OC) theory to design Radio-Frequency (RF) pulses that actively control the spatial distribution of the MRI magnetization phase. The RF pulses are generated through the application of the Pontryagin Maximum Principle and optimized so that the resulting transverse magnetization reproduces various non-trivial and spatial phase patterns. Two different phase patterns are defined and the resulting optimal pulses are tested both numerically with the ODIN MRI simulator and experimentally with an agar gel phantom on a 4.7T small-animal MR scanner. Phase images obtained in simulations and experiments are both consistent with the defined phase patterns. A practical application of phase control with OC-designed pulses is also presented, with the generation of RF pulses adapted for a Magnetic Resonance Elastography experiment. This study demonstrates the possibility to use OC-designed RF pulses to encode information in the magnetization phase and could have applications in MRI sequences using phase images.

Implementing a universal gate set on a logical qubit encoded in an oscillator

A logical qubit is a two-dimensional subspace of a higher dimensional system, chosen such that it is possible to detect and correct the occurrence of certain errors. Manipulation of the encoded information generally requires arbitrary and precise control over the entire system. Whether based on multiple physical qubits or larger dimensional modes such as oscillators, the individual elements in realistic devices will always have residual interactions, which must be accounted for when designing logical operations. Here we demonstrate a holistic control strategy which exploits accurate knowledge of the Hamiltonian to manipulate a coupled oscillator-transmon system. We use this approach to realize high-fidelity (98.5%, inferred), decoherence-limited operations on a logical qubit encoded in a superconducting cavity resonator using four-component cat states. Our results show the power of applying numerical techniques to control linear oscillators and pave the way for utilizing their large Hilbert space as a resource in quantum information processing.A logical qubit is a two-dimensional subspace of a higher dimensional system, whose manipulation requires precise control over the whole system. Here the authors demonstrate a control strategy which exploits precise knowledge of the Hamiltonian to manipulate a coupled oscillator-transmon system.

Perspectives of shaped pulses for EPR spectroscopy

This article describes current uses of shaped pulses, generated by an arbitrary waveform generator, in the field of EPR spectroscopy. We show applications of sech/tanh and WURST pulses to dipolar spectroscopy, including new pulse schemes and procedures, and discuss the more general concept of optimum-control-based pulses for applications in EPR spectroscopy. The article also describes a procedure to correct for experimental imperfections, mostly introduced by the microwave resonator, and discusses further potential applications and limitations of such pulses.

Optimal control design of preparation pulses for contrast optimization in MRI

This work investigates the use of MRI radio-frequency (RF) pulses designed within the framework of optimal control theory for image contrast optimization. The magnetization evolution is modeled with Bloch equations, which defines a dynamic system that can be controlled via the application of the Pontryagin Maximum Principle (PMP). This framework allows the computation of optimal RF pulses that bring the magnetization to a given state to obtain the desired contrast after acquisition. Creating contrast through the optimal manipulation of Bloch equations is a new way of handling contrast in MRI, which can explore the theoretical limits of the system. Simulation experiments carried out on-resonance quantify the contrast improvement when compared to standard T₁ or T₂ weighting strategies. The use of optimal pulses is also validated for the first time in both in vitro and in vivo experiments on a small-animal 4.7T MR system. Results demonstrate their robustness to static field inhomogeneities as well as the fact that they can be embedded in standard imaging sequences without affecting standard parameters such as slice selection or echo type. In vivo results on rat and mouse brains illustrate the ability of optimal contrast pulses to create non-trivial contrasts on well-studied structures (white matter versus gray matter).

Hybrid Quantum-Classical Approach to Quantum Optimal Control

A central challenge in quantum computing is to identify more computational problems for which utilization of quantum resources can offer significant speedup. Here, we propose a hybrid quantum-classical scheme to tackle the quantum optimal control problem. We show that the most computationally demanding part of gradient-based algorithms, namely, computing the fitness function and its gradient for a control input, can be accomplished by the process of evolution and measurement on a quantum simulator. By posing queries to and receiving answers from the quantum simulator, classical computing devices update the control parameters until an optimal control solution is found. To demonstrate the quantum-classical scheme in experiment, we use a seven-qubit nuclear magnetic resonance system, on which we have succeeded in optimizing state preparation without involving classical computation of the large Hilbert space evolution.

Local SAR, global SAR, and power-constrained large-flip-angle pulses with optimal control and virtual observation points

PURPOSE

To present a constrained optimal-control (OC) framework for designing large-flip-angle parallel-transmit (pTx) pulses satisfying hardware peak-power as well as regulatory local and global specific-absorption-rate (SAR) limits. The application is 2D and 3D spatial-selective 90° and 180° pulses.

THEORY AND METHODS

The OC gradient-ascent-pulse-engineering method with exact gradients and the limited-memory Broyden-Fletcher-Goldfarb-Shanno method is proposed. Local SAR is constrained by the virtual-observation-points method. Two numerical models facilitated the optimizations, a torso at 3 T and a head at 7 T, both in eight-channel pTx coils and acceleration-factors up to 4.

RESULTS

The proposed approach yielded excellent flip-angle distributions. Enforcing the local-SAR constraint, as opposed to peak power alone, reduced the local SAR 7 and 5-fold with the 2D torso excitation and inversion pulse, respectively. The root-mean-square errors of the magnetization profiles increased less than 5% with the acceleration factor of 4.

CONCLUSION

A local and global SAR, and peak-power constrained OC large-flip-angle pTx pulse design was presented, and numerically validated for 2D and 3D spatial-selective 90° and 180° pulses at 3 T and 7 T. Magn Reson Med 77:374-384, 2017. © 2015 Wiley Periodicals, Inc.

Efficient simulation of ultrafast magnetic resonance experiments

We present a convenient and powerful simulation formalism for ultrafast NMR spectroscopy. The formalism is based on the Fokker–Planck equation that supports systems with complicated combinations of classical spatial dynamics and quantum mechanical spin dynamics.