Showing 1–19 of 19 results for author: Goodwin, D

A universal neutral-atom quantum computer with individual optical addressing and non-destructive readout

Authors: A. G. Radnaev, W. C. Chung, D. C. Cole, D. Mason, T. G. Ballance, M. J. Bedalov, D. A. Belknap, M. R. Berman, M. Blakely, I. L. Bloomfield, P. D. Buttler, C. Campbell, A. Chopinaud, E. Copenhaver, M. K. Dawes, S. Y. Eubanks, A. J. Friss, D. M. Garcia, J. Gilbert, M. Gillette, P. Goiporia, P. Gokhale, J. Goldwin, D. Goodwin, T. M. Graham , et al. (33 additional authors not shown)

Abstract: Quantum computers must achieve large-scale, fault-tolerant operation to deliver on their promise of transformational processing power [1-4]. This will require thousands or millions of high-fidelity quantum gates and similar numbers of qubits [5]. Demonstrations using neutral-atom qubits trapped and manipulated by lasers have shown that this modality can provide high two-qubit gate (CZ) fidelities… ▽ More Quantum computers must achieve large-scale, fault-tolerant operation to deliver on their promise of transformational processing power [1-4]. This will require thousands or millions of high-fidelity quantum gates and similar numbers of qubits [5]. Demonstrations using neutral-atom qubits trapped and manipulated by lasers have shown that this modality can provide high two-qubit gate (CZ) fidelities and scalable operation [6-10]. However, the gates in these demonstrations are driven by lasers that do not resolve individual qubits, with universal computation enabled by physical mid-circuit shuttling of the qubits. This relatively slow operation will greatly extend runtimes for useful, large-scale computation. Here we demonstrate a universal neutral-atom quantum computer with gate rates limited by optical switching times, rather than shuttling, by individually addressing tightly focused laser beams at an array of single atoms. We achieve CZ fidelity of 99.35(4)% and local single qubit RZ gate fidelity of 99.902(8)%. Moreover, we demonstrate non-destructive readout of alkali-atom qubits with sub-percent loss, which boosts operational speed. This technique also enables us to measure 99.73(3)% CZ fidelity with atom-loss events excluded, which is a record among long lived neutral-atom qubits and highlights the path to higher fidelity and error correction. Our results represent a critical step towards large-scale, fault-tolerant neutral-atom quantum computers that can execute computations on practical timescales. △ Less

Submitted 19 August, 2024; v1 submitted 15 August, 2024; originally announced August 2024.

Seeing Through the Clouds: Cloud Gap Imputation with Prithvi Foundation Model

Authors: Denys Godwin, Hanxi Li, Michael Cecil, Hamed Alemohammad

Abstract: Filling cloudy pixels in multispectral satellite imagery is essential for accurate data analysis and downstream applications, especially for tasks which require time series data. To address this issue, we compare the performance of a foundational Vision Transformer (ViT) model with a baseline Conditional Generative Adversarial Network (CGAN) model for missing value imputation in time series of mul… ▽ More Filling cloudy pixels in multispectral satellite imagery is essential for accurate data analysis and downstream applications, especially for tasks which require time series data. To address this issue, we compare the performance of a foundational Vision Transformer (ViT) model with a baseline Conditional Generative Adversarial Network (CGAN) model for missing value imputation in time series of multispectral satellite imagery. We randomly mask time series of satellite images using real-world cloud masks and train each model to reconstruct the missing pixels. The ViT model is fine-tuned from a pretrained model, while the CGAN is trained from scratch. Using quantitative evaluation metrics such as structural similarity index and mean absolute error as well as qualitative visual analysis, we assess imputation accuracy and contextual preservation. △ Less

Submitted 30 April, 2024; originally announced April 2024.

Low-loss, compact, fibre-integrated cell for quantum memories

Authors: Cameron McGarry, Kerrianne Harrington, Daniel J. Goodwin, Charles Perek-Jennings, Tim A. Birks, Kristina Rusimova, Peter J. Mosley

Abstract: We present a low-loss, compact, hollow core optical fibre (HCF) cell integrated with single mode fibre (SMF). The cell is designed to be filled with atomic vapour and used as a component in photonic quantum technologies, with applications in quantum memory and optical switching. We achieve a total insertion loss of 0.6(2) dB at 780 nm wavelength via graded index fibre to ensure efficient mode matc… ▽ More We present a low-loss, compact, hollow core optical fibre (HCF) cell integrated with single mode fibre (SMF). The cell is designed to be filled with atomic vapour and used as a component in photonic quantum technologies, with applications in quantum memory and optical switching. We achieve a total insertion loss of 0.6(2) dB at 780 nm wavelength via graded index fibre to ensure efficient mode matching coupled with anti-reflection coatings to minimise loss at the SMF-HCF interfaces. We also present numerical modelling of these interfaces, which can be undertaken efficiently without the need for finite element simulation. We encapsulate the HCF core by coupling to the SMF inside a support capillary, enhancing durability and facilitating seamless integration into existing fibre platforms. △ Less

Submitted 19 January, 2024; originally announced January 2024.

Foundation Models for Generalist Geospatial Artificial Intelligence

Authors: Johannes Jakubik, Sujit Roy, C. E. Phillips, Paolo Fraccaro, Denys Godwin, Bianca Zadrozny, Daniela Szwarcman, Carlos Gomes, Gabby Nyirjesy, Blair Edwards, Daiki Kimura, Naomi Simumba, Linsong Chu, S. Karthik Mukkavilli, Devyani Lambhate, Kamal Das, Ranjini Bangalore, Dario Oliveira, Michal Muszynski, Kumar Ankur, Muthukumaran Ramasubramanian, Iksha Gurung, Sam Khallaghi, Hanxi, Li , et al. (8 additional authors not shown)

Abstract: Significant progress in the development of highly adaptable and reusable Artificial Intelligence (AI) models is expected to have a significant impact on Earth science and remote sensing. Foundation models are pre-trained on large unlabeled datasets through self-supervision, and then fine-tuned for various downstream tasks with small labeled datasets. This paper introduces a first-of-a-kind framewo… ▽ More Significant progress in the development of highly adaptable and reusable Artificial Intelligence (AI) models is expected to have a significant impact on Earth science and remote sensing. Foundation models are pre-trained on large unlabeled datasets through self-supervision, and then fine-tuned for various downstream tasks with small labeled datasets. This paper introduces a first-of-a-kind framework for the efficient pre-training and fine-tuning of foundational models on extensive geospatial data. We have utilized this framework to create Prithvi, a transformer-based geospatial foundational model pre-trained on more than 1TB of multispectral satellite imagery from the Harmonized Landsat-Sentinel 2 (HLS) dataset. Our study demonstrates the efficacy of our framework in successfully fine-tuning Prithvi to a range of Earth observation tasks that have not been tackled by previous work on foundation models involving multi-temporal cloud gap imputation, flood mapping, wildfire scar segmentation, and multi-temporal crop segmentation. Our experiments show that the pre-trained model accelerates the fine-tuning process compared to leveraging randomly initialized weights. In addition, pre-trained Prithvi compares well against the state-of-the-art, e.g., outperforming a conditional GAN model in multi-temporal cloud imputation by up to 5pp (or 5.7%) in the structural similarity index. Finally, due to the limited availability of labeled data in the field of Earth observation, we gradually reduce the quantity of available labeled data for refining the model to evaluate data efficiency and demonstrate that data can be decreased significantly without affecting the model's accuracy. The pre-trained 100 million parameter model and corresponding fine-tuning workflows have been released publicly as open source contributions to the global Earth sciences community through Hugging Face. △ Less

Submitted 8 November, 2023; v1 submitted 28 October, 2023; originally announced October 2023.

Accelerated Newton-Raphson GRAPE methods for optimal control

Authors: David L. Goodwin, Mads Sloth Vinding

Abstract: A Hessian based optimal control method is presented in Liouville space to mitigate previously undesirable polynomial scaling of computation time. This new method, an improvement to the state-of-the-art Newton-Raphson GRAPE method, is derived with respect to two exact time-propagator derivative techniques: auxiliary matrix and ESCALADE methods. We observed that compared to the best current implemen… ▽ More A Hessian based optimal control method is presented in Liouville space to mitigate previously undesirable polynomial scaling of computation time. This new method, an improvement to the state-of-the-art Newton-Raphson GRAPE method, is derived with respect to two exact time-propagator derivative techniques: auxiliary matrix and ESCALADE methods. We observed that compared to the best current implementation of Newton-Raphson GRAPE method, for an ensemble of 2-level systems, with realistic conditions, the new auxiliary matrix and ESCALADE Hessians can be 4-200 and 70-600 times faster, respectively. △ Less

Submitted 18 August, 2022; v1 submitted 20 July, 2022; originally announced July 2022.

Concurrent J-Evolving Refocusing Pulses

Authors: Sebastian Ehni, Martin R. M. Koos, Tony Reinsperger, Jens D. Haller, David L. Goodwin, Burkhard Luy

Abstract: 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 stud… ▽ More 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 detail 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 refocusing pulse sandwiches -- termed BUBU pulses following the nomenclature of previous J-compensated BUBI and BEBE(tr) pulse sandwiches -- have been optimised. △ Less

Submitted 12 September, 2021; originally announced September 2021.

Comments: 17 pages, 5 figures

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

Authors: Mads Sloth Vinding, David L. Goodwin, Ilya Kuprov, Torben Ellegaard Lund

Abstract: 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-spe… ▽ More 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 fidelity 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 fidelity 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 truncation is found to be sufficiently accurate, but also up to five times faster than the machine precision gradient. This makes the generation of training databases for the machine learning methods considerably more realistic. △ Less

Submitted 2 July, 2021; originally announced July 2021.

Comments: The following article has been submitted to the Journal of Chemical Physics. After it is published, it will be found at https://aip.scitation.org/journal/jcp

Serverless inferencing on Kubernetes

Authors: Clive Cox, Dan Sun, Ellis Tarn, Animesh Singh, Rakesh Kelkar, David Goodwin

Abstract: Organisations are increasingly putting machine learning models into production at scale. The increasing popularity of serverless scale-to-zero paradigms presents an opportunity for deploying machine learning models to help mitigate infrastructure costs when many models may not be in continuous use. We will discuss the KFServing project which builds on the KNative serverless paradigm to provide a s… ▽ More Organisations are increasingly putting machine learning models into production at scale. The increasing popularity of serverless scale-to-zero paradigms presents an opportunity for deploying machine learning models to help mitigate infrastructure costs when many models may not be in continuous use. We will discuss the KFServing project which builds on the KNative serverless paradigm to provide a serverless machine learning inference solution that allows a consistent and simple interface for data scientists to deploy their models. We will show how it solves the challenges of autoscaling GPU based inference and discuss some of the lessons learnt from using it in production. △ Less

Submitted 24 July, 2020; v1 submitted 14 July, 2020; originally announced July 2020.

Comments: 4 pages, 1 figure, presented at workshop on "Challenges in Deploying and Monitoring Machine Learning System" at ICML 2020

Intra-model Variability in COVID-19 Classification Using Chest X-ray Images

Authors: Brian D Goodwin, Corey Jaskolski, Can Zhong, Herick Asmani

Abstract: X-ray and computed tomography (CT) scanning technologies for COVID-19 screening have gained significant traction in AI research since the start of the coronavirus pandemic. Despite these continuous advancements for COVID-19 screening, many concerns remain about model reliability when used in a clinical setting. Much has been published, but with limited transparency in expected model performance. W… ▽ More X-ray and computed tomography (CT) scanning technologies for COVID-19 screening have gained significant traction in AI research since the start of the coronavirus pandemic. Despite these continuous advancements for COVID-19 screening, many concerns remain about model reliability when used in a clinical setting. Much has been published, but with limited transparency in expected model performance. We set out to address this limitation through a set of experiments to quantify baseline performance metrics and variability for COVID-19 detection in chest x-ray for 12 common deep learning architectures. Specifically, we adopted an experimental paradigm controlling for train-validation-test split and model architecture where the source of prediction variability originates from model weight initialization, random data augmentation transformations, and batch shuffling. Each model architecture was trained 5 separate times on identical train-validation-test splits of a publicly available x-ray image dataset provided by Cohen et al. (2020). Results indicate that even within model architectures, model behavior varies in a meaningful way between trained models. Best performing models achieve a false negative rate of 3 out of 20 for detecting COVID-19 in a hold-out set. While these results show promise in using AI for COVID-19 screening, they further support the urgent need for diverse medical imaging datasets for model training in a way that yields consistent prediction outcomes. It is our hope that these modeling results accelerate work in building a more robust dataset and a viable screening tool for COVID-19. △ Less

Submitted 30 April, 2020; originally announced May 2020.

Comments: 7 pages, 5 figures; Writing, analysis, and design carried out by authors Brian and Corey; experiments carried out by authors Can and Herick; results and code located at https://github.com/synthetaic/COVID19-IntraModel-Variability and https://covidresearch.ai/datasets/dataset?id=2

MSC Class: 68T01 ACM Class: I.2.0

Dogfooding: use IBM Cloud services to monitor IBM Cloud infrastructure

Authors: William Pourmajidi, Andriy Miranskyy, John Steinbacher, Tony Erwin, David Godwin

Abstract: The stability and performance of Cloud platforms are essential as they directly impact customers' satisfaction. Cloud service providers use Cloud monitoring tools to ensure that rendered services match the quality of service requirements indicated in established contracts such as service-level agreements. Given the enormous number of resources that need to be monitored, highly scalable and capable… ▽ More The stability and performance of Cloud platforms are essential as they directly impact customers' satisfaction. Cloud service providers use Cloud monitoring tools to ensure that rendered services match the quality of service requirements indicated in established contracts such as service-level agreements. Given the enormous number of resources that need to be monitored, highly scalable and capable monitoring tools are designed and implemented by Cloud service providers such as Amazon, Google, IBM, and Microsoft. Cloud monitoring tools monitor millions of virtual and physical resources and continuously generate logs for each one of them. Considering that logs magnify any technical issue, they can be used for disaster detection, prevention, and recovery. However, logs are useless if they are not assessed and analyzed promptly. Thus, we argue that the scale of Cloud-generated logs makes it impossible for DevOps teams to analyze them effectively. This implies that one needs to automate the process of monitoring and analysis (e.g., using machine learning and artificial intelligence). If the automation will witness an anomaly in the logs --- it will alert DevOps staff. The automatic anomaly detectors require a reliable and scalable platform for gathering, filtering, and transforming the logs, executing the detector models, and sending out the alerts to the DevOps staff. In this work, we report on implementing a prototype of such a platform based on the 7-layered architecture pattern, which leverages micro-service principles to distribute tasks among highly scalable, resources-efficient modules. The modules interact with each other via an instance of the Publish-Subscribe architectural pattern. The platform is deployed on the IBM Cloud service infrastructure and is used to detect anomalies in logs emitted by the IBM Cloud services, hence the dogfooding. △ Less

Submitted 13 July, 2019; originally announced July 2019.

Journal ref: Proceedings of the 29th Annual International Conference on Computer Science and Software Engineering (CASCON'19), 2019, pp. 344-353

Second order phase dispersion by optimised rotation pulses

Authors: David L. Goodwin, Martin R. M. Koos, Burkhard Luy

Abstract: We show that the duration of broadband universal control pulses can be halved by choosing control targets with a quadratic function of phase dispersion. This class of control pulses perform a broadband universal rotation around an axis, in the Bloch sphere representation of two-level systems, given by this phase dispersion function. We present an effective optimal control method to avoid the probl… ▽ More We show that the duration of broadband universal control pulses can be halved by choosing control targets with a quadratic function of phase dispersion. This class of control pulses perform a broadband universal rotation around an axis, in the Bloch sphere representation of two-level systems, given by this phase dispersion function. We present an effective optimal control method to avoid the problem of convergence to local extrema traps. △ Less

Submitted 2 August, 2020; v1 submitted 20 April, 2019; originally announced April 2019.

Comments: 8 pages, 4 figures

Journal ref: Phys. Rev. Research 2, 033157 (2020)

Optimal control of mirror pulses for atom interferometry

Authors: Jack Saywell, Ilya Kuprov, David Goodwin, Max Carey, Tim Freegarde

Abstract: Atom matterwave interferometry requires mirror and beamsplitter pulses that are robust to inhomogeneities in field intensity, magnetic environment, atom velocity and Zeeman sub-state. Pulse shapes determined using quantum control methods offer significantly improved interferometer performance by allowing broader atom distributions, larger interferometer areas and higher contrast. We have applied g… ▽ More Atom matterwave interferometry requires mirror and beamsplitter pulses that are robust to inhomogeneities in field intensity, magnetic environment, atom velocity and Zeeman sub-state. Pulse shapes determined using quantum control methods offer significantly improved interferometer performance by allowing broader atom distributions, larger interferometer areas and higher contrast. We have applied gradient ascent pulse engineering (GRAPE) to optimise the design of phase-modulated mirror pulses for a Mach-Zehnder light-pulse atom interferometer, with the aim of increasing fringe contrast when averaged over atoms with an experimentally relevant range of velocities, beam intensities, and Zeeman states. Pulses were found to be highly robust to variations in detuning and coupling strength, and offer a clear improvement in robustness over the best established composite pulses. The peak mirror fidelity in a cloud of $\sim 80\ μ$K ${}^{85}$Rb atoms is predicted to be improved by a factor of 2 compared with standard rectangular $π$ pulses. △ Less

Submitted 19 April, 2018; v1 submitted 12 April, 2018; originally announced April 2018.

Journal ref: Phys. Rev. A 98, 023625 (2018)

Advanced Optimal Control Methods for Spin Systems

Authors: David L. Goodwin

Abstract: Work within this thesis advances optimal control algorithms for application to magnetic resonance systems. Specifically, presenting a quadratically convergent version of the gradient ascent pulse engineering method. The work is formulated in a superoperator representation of the Liouville-von Neumann equation. A Newton-grape method is developed using efficient calculation of analytical second di… ▽ More Work within this thesis advances optimal control algorithms for application to magnetic resonance systems. Specifically, presenting a quadratically convergent version of the gradient ascent pulse engineering method. The work is formulated in a superoperator representation of the Liouville-von Neumann equation. A Newton-grape method is developed using efficient calculation of analytical second directional derivatives. The method is developed to scale with the same complexity as methods that use only first directional derivatives. Algorithms to ensure a well-conditioned and positive definite matrix of second directional derivatives are used so the sufficient conditions of gradient-based numerical optimisation are met. A number of applications of optimal control in magnetic resonance are investigated: solid-state nuclear magnetic resonance, magnetisation-to-singlet pulses, and electron spin resonance experiments. △ Less

Submitted 26 April, 2018; v1 submitted 28 March, 2018; originally announced March 2018.

Comments: 250 pages, 56 figures

Journal ref: David L. Goodwin. Advanced Optimal Control Methods for Spin Systems. PhD thesis, University of Southampton, 2017

Feedback control optimisation of ESR experiments

Authors: David L. Goodwin, William K. Myers, Christiane R. Timmel, Ilya Kuprov

Abstract: 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… ▽ More 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. △ Less

Submitted 18 October, 2018; v1 submitted 13 July, 2017; originally announced July 2017.

Journal ref: J. Magn. Reson. 297, 9-16 (2018)

Database Engines: Evolution of Greenness

Authors: Andriy V. Miranskyy, Zainab Al-zanbouri, David Godwin, Ayse Basar Bener

Abstract: Context: Information Technology consumes up to 10\% of the world's electricity generation, contributing to CO2 emissions and high energy costs. Data centers, particularly databases, use up to 23% of this energy. Therefore, building an energy-efficient (green) database engine could reduce energy consumption and CO2 emissions. Goal: To understand the factors driving databases' energy consumption a… ▽ More Context: Information Technology consumes up to 10\% of the world's electricity generation, contributing to CO2 emissions and high energy costs. Data centers, particularly databases, use up to 23% of this energy. Therefore, building an energy-efficient (green) database engine could reduce energy consumption and CO2 emissions. Goal: To understand the factors driving databases' energy consumption and execution time throughout their evolution. Method: We conducted an empirical case study of energy consumption by two MySQL database engines, InnoDB and MyISAM, across 40 releases. We examined the relationships of four software metrics to energy consumption and execution time to determine which metrics reflect the greenness and performance of a database. Results: Our analysis shows that database engines' energy consumption and execution time increase as databases evolve. Moreover, the Lines of Code metric is correlated moderately to strongly with energy consumption and execution time in 88% of cases. Conclusions: Our findings provide insights to both practitioners and researchers. Database administrators may use them to select a fast, green release of the MySQL database engine. MySQL database-engine developers may use the software metric to assess products' greenness and performance. Researchers may use our findings to further develop new hypotheses or build models to predict greenness and performance of databases. △ Less

Submitted 9 January, 2017; originally announced January 2017.

Modified Newton-Raphson GRAPE methods for optimal control of spin systems

Authors: D. L. Goodwin, Ilya Kuprov

Abstract: 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… ▽ More 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. △ Less

Submitted 28 April, 2016; v1 submitted 8 October, 2015; originally announced October 2015.

Comments: submitted for publication

Auxiliary matrix formalism for interaction representation transformations, optimal control and spin relaxation theories

Authors: D. L. Goodwin, Ilya Kuprov

Abstract: Auxiliary matrix exponential method is used to derive simple and numerically efficient general expressions for the following, historically rather cumbersome and hard to compute, theoretical methods: (1) average Hamiltonian theory following interaction representation transformations; (2) Bloch-Redfield-Wangsness theory of nuclear and electron relaxation; (3) gradient ascent pulse engineering versio… ▽ More Auxiliary matrix exponential method is used to derive simple and numerically efficient general expressions for the following, historically rather cumbersome and hard to compute, theoretical methods: (1) average Hamiltonian theory following interaction representation transformations; (2) Bloch-Redfield-Wangsness theory of nuclear and electron relaxation; (3) gradient ascent pulse engineering version of quantum optimal control theory. In the context of spin dynamics, the auxiliary matrix exponential method is more efficient than methods based on matrix factorizations and also exhibits more favourable complexity scaling with the dimension of the Hamiltonian matrix. △ Less

Submitted 24 July, 2015; v1 submitted 29 May, 2015; originally announced June 2015.

Impedance Spectra of Mixed Conductors: a 2D Study of Ceria

Authors: Francesco Ciucci, Yong Hao, David G. Goodwin

Abstract: In this paper we develop an analytical framework for the study of electrochemical impedance of mixed ionic and electronic conductors (MIEC). The framework is based on first-principles and it features the coupling of electrochemical reactions, surface transport and bulk transport processes. We utilize this work to analyze two dimensional systems relevant for fuel cell science via finite element m… ▽ More In this paper we develop an analytical framework for the study of electrochemical impedance of mixed ionic and electronic conductors (MIEC). The framework is based on first-principles and it features the coupling of electrochemical reactions, surface transport and bulk transport processes. We utilize this work to analyze two dimensional systems relevant for fuel cell science via finite element method (FEM). Alternate current Impedance Spectroscopy (IS) of a ceria symmetric cell is simulated near equilibrium condition (zero bias) for a wide array of working conditions including variations of temperature and $H_2$ partial pressure on a two dimensional fuel cell sample with patterned metal electrodes. The model shows agreement of IS curves with the experimental literature with the relative error on the impedance being consistently below 2%. Important two-dimensional effects such the effects of thickness decrease and the influence of variable electronic and ionic diffusivities on the impedance spectra are also explored. △ Less

Submitted 16 April, 2009; v1 submitted 20 March, 2009; originally announced March 2009.

Two-dimensional electrochemical model for mixed conductors: a study of ceria

Authors: Francesco Ciucci, William C. Chueh, Sossina M. Haile, David G. Goodwin

Abstract: A two-dimensional small bias model has been developed for a patterned metal current collector $|$ mixed oxygen ion and electronic conductor (MIEC) $|$ patterned metal current collector electrochemical cell in a symmetric gas environment. Specifically, we compute the electrochemical potential distributions of oxygen vacancies and electrons in the bulk and near the surface for… ▽ More A two-dimensional small bias model has been developed for a patterned metal current collector $|$ mixed oxygen ion and electronic conductor (MIEC) $|$ patterned metal current collector electrochemical cell in a symmetric gas environment. Specifically, we compute the electrochemical potential distributions of oxygen vacancies and electrons in the bulk and near the surface for $\text{Pt} | \text{Sm}_{0.15}\text{Ce}_{0.85}\text{O}_{1.925} | \text{Pt}$ symmetric cell in a $\text{H}_2-\text{H}_2\text{O}-\text{Ar}$ (reducing) atmosphere from 500 to $650^o C$. Using a two-dimensional finite-element model, we show that two types of electronic current exist within the cell: an in-plane drift-diffusion current that flows between the gas $|$ ceria chemical reaction site and the metal current collector, and a cross-plane current that flows between the two metal electrodes on the opposite side of the cell. By fitting the surface reaction constant $\tilde k_f^0$ to experimental electrode resistance values while fixing material properties such as bulk ionic and electronic equilibrium defect concentrations and mobilities, we are able to separate the electrode polarization into the surface reaction component and the in-plane electron drift-diffusion component. We show that for mixed conductors with a low electronic conductivity (a function of oxygen partial pressure) or a high surface reaction rate constant, the in-plane electron drift-diffusion resistance can become rate-limiting in the electrode reaction. △ Less

Submitted 18 March, 2009; originally announced March 2009.