Search NASA⌕ Search

SEARCH · Search NASA

Results for “Distribution System State Estimation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Pre-commercial Demonstration of Direct Non-iterative State Estimator for Operational Use of Synchrophasor Management Systems

The objective of this project was to demonstrate functionality and performance of a Direct Non-iterative State Estimator (DNSE) integrated with NYPA’s Energy Management System (EMS) and with enhanced Real Time Dynamics Monitoring System (RTDMS) synchrophasor platform from Electric Power Group (EPG). DNSE is designed to overcome a major obstacle to operational use of Synchro-Phasor Management Systems (SPMS) by providing to synchrophasor management systems (SPMS) applications a consistent and a complete synchrophasor data foundation in the same way that a traditional EMS State Estimator (SE) provides to EMS applications. Specifically, DNSE is designed to use synchrophasor measurements collected by a central PDC, Supervisory Control and Data Acquisition (SCADA) measurements, and Energy Management System (EMS) network model, to obtain the complete state of the utility’s operating model at rates that are close to the synchrophasor data rates. In this way, the system is comprehensive in that it does not only cover the part of the network that is “visible” via synchrophasors, but also the part that is only “visible” through the SCADA measurements. Visualization needs associated with the use of DNSE results are fulfilled through suitably enhanced Real Time Dynamics Monitoring System (RTDMS), with the enhancements implemented by EPG. This project had the following goals in mind: To advance the deployment of commercial grade DNSE software application that relies on synchrophasor and SCADA data ; Apply DNSE at other utilities, to address a generic and fundamental need for “clean” operational data for synchrophasor applications; Provide means for “live” estimated data access by control system operators; Enhance potential for situational awareness through full system operational model coverage; Sub-second execution rate of the Direct Non-iterative State Estimator, eventually at a near-phasor data rate execution speed, i.e. < 0.1 sec. Anticipated benefits from this projects are: Enhanced reliability and improvements in the economic efficiency of bulk power system planning and operations; Providing “clean” data to other synchrophasor applications; Enhancement of situational awareness by providing the full operational model updated at near synchrophasor rate; A production-grade software tool that incorporate synchrophasor and SCADA data; Provides a basis for development of next generation monitoring and control applications, based on both SCADA and PMU data. Quanta Technology (QT) team worked in collaboration with Electric Power Group (EPG) whose team has enhanced its commercial Real Time Dynamics Monitoring System (RTDMS) to accommodate the requirements posed by DNSE application. EPG also provided its ePDC and Model-less Data Conditioning (PDVC) software for integration with DNSE+. QT developed the system requirements for DNSE; developed system architecture and defined interfaces between internal DNSE components. The core DNSE algorithm with all surrounding interfaces was named DNSE+. Since the DNSE development was done in a simulated system environment, QT used its PMU simulator that was enhanced during this project for development and factory acceptance testing (FAT). SCADA data in this stage was simulated by commercial PSS/e software. The output of DNSE are estimates of System states in C37.118-2 format, sent to RTDMS for further processing and display. As the number of these states is large, it was necessary to expand the C37.111-2 standard to accommodate large data sets. This enhancement was implemented in RTDMS. The demonstration of pre-production DNSE technology was done at NYPA using streaming field data from NYPA PMUs and from its RTUs through their SCADA system. NYPA provided ICCP interface as well as Common Information Model (CIM). The relevance of the DNSE+ application is that it provides state estimation of the power system based on hybrid set of data, consisting of both available PMU data and SCADA measurements. As this is a direct, non-iterative method of calculation of the system states, if does not suffer from convergence issues which is potential problem for conventional state estimators. Also, it can take any available PMU measurements, so it does not need to have a high percentage of PMU coverage needed in the case of Linear State Estimator. As the DNSE calculates synchrophasors of the system states (both phase and absolute value) as sub-second rate, this application can provide a basis for development of next generation of applications based both on SCADA and PMU data.

24 POWER TRANSMISSION AND DISTRIBUTION↗
NASA NTRS1993

Recent advances in radar remote sensing of forest

On a global scale, forests represent most of the terrestrial standing biomass (80 to 90 percent). Thus, natural and anthropogenic change in forest covers can have major impacts not only on local ecosystems but also on global hydrologic, climatic, and biogeochemical cycles that involve exchange of energy, water, carbon, and other elements between the earth and atmosphere. Quantitative information on the state and dynamics of forest ecosystems and their interactions with the global cycles appear necessary to understand how the earth works as a natural system. The information required includes the lateral and vertical distribution of forest cover, the estimates of standing biomass (woody and foliar volume), the phenological and environmental variations and disturbances (clearcutting, fires, flood), and the longer term variations following deforestation (regeneration, successional stages). To this end, seasonal, annual, and decadal information is necessary in order to separate the long term effects in the global ecosystem from short term seasonal and interannual variations. Optical remote sensing has been used until now to study the forest cover at local, regional, and global scales. Radar remote sensing, which provides recent SAR data from space on a regular basis, represents an unique means of consistently monitoring different time scales, at all latitudes and in any atmospheric conditions. Also, SAR data have shown the potential to detect several forest parameters that cannot be inferred from optical data. The differences--and complementarity--lie in the penetration capabilities of SAR data and their sensitivity to dielectric and geometric properties of the canopy volume, whereas optical data are sensitive to the chemical composition of the external foliar layer of the vegetation canopy.

Letoan, Thuy↗
NASA NTRS1995

An approximate Kalman filter for ocean data assimilation: An example with an idealized Gulf Stream model

A practical method of data assimilation for use with large, nonlinear, ocean general circulation models is explored. A Kalman filter based on approximation of the state error covariance matrix is presented, employing a reduction of the effective model dimension, the error's asymptotic steady state limit, and a time-invariant linearization of the dynamic model for the error integration. The approximations lead to dramatic computational savings in applying estimation theory to large complex systems. We examine the utility of the approximate filter in assimilating different measurement types using a twin experiment of an idealized Gulf Stream. A nonlinear primitive equation model of an unstable east-west jet is studied with a state dimension exceeding 170,000 elements. Assimilation of various pseudomeasurements are examined, including velocity, density, and volume transport at localized arrays and realistic distributions of satellite altimetry and acoustic tomography observations. Results are compared in terms of their effects on the accuracies of the estimation. The approximate filter is shown to outperform an empirical nudging scheme used in a previous study. The examples demonstrate that useful approximate estimation errors can be computed in a practical manner for general circulation models.

Fukumori, Ichiro↗
DOE OSTI2015

Pre-Commercial Demonstration of Direct Non-Iterative State Estimator (DNSE+)

The objective of this DOE co-funded project is to demonstrate functionality and performance of a production-grade Direct Non-iterative State Estimator (DNSE) through its integration with NYPA’s Energy Management System (EMS) and an enhanced Real Time Dynamics Monitoring System (RTDMS) synchrophasor platform from Electric Power Group (EPG). A production-grade DNSE aims to overcome a major obstacle to the operational use of synchrophasor data in Synchro-Phasor Management Systems (SPMS) by providing SPMS applications with a consistent and a complete synchrophasor data foundation in the same way that a traditional EMS State Estimator (SE) provides to EMS applications. Specifically, the DNSE uses synchrophasor measurements and Supervisory Control and Data Acquisition (SCADA) measurements, together with the complete power system model, to simultaneously obtain the complete state of the entire utility’s operating model at rates that are close to the synchrophasor data rates. The full system coverage property is what distinguishes the DNSE from the linear synchrophasor SEs (which cover only parts of the operating model that are visible through synchrophasor measurements), and DNSE differs from traditional EMS SE in that the solution is obtained without iterations, without initial values (hence the “direct” in the name), and at a much higher execution rate. This presentation provides an overview of the project objectives, a brief description of the main tasks including the development work to productize and optimize DNSE from its current MATLAB code base, RTDMS enhancement for interfacing with the DNSE, factory acceptance testing, field installation, testing and demonstration, as well as the current status of the project execution.

24 POWER TRANSMISSION AND DISTRIBUTION↗
NASA NTRS1993

An atlas of monthly mean distributions of SSMI surface wind speed, ARGOS buoy drift, AVHRR/2 sea surface temperature, and ECMWF surface wind components during 1990

The following monthly mean global distributions for 1990 are proposed with a common color scale and geographical map: 10-m height wind speed estimated from the Special Sensor Microwave Imager (SSMI) on a United States (US) Air Force Defense Meteorological Satellite Program (DMSP) spacecraft; sea surface temperature estimated from the advanced very high resolution radiometer (AVHRR/2) on a U.S. National Oceanic and Atmospheric Administration (NOAA) spacecraft; Cartesian components of free drifting buoys which are tracked by the ARGOS navigation system on NOAA satellites; and Cartesian components on the 10-m height wind vector computed by the European Center for Medium-Range Weather Forecasting (ECMWF). Charts of monthly mean value, sampling distribution, and standard deviation values are displayed. Annual mean distributions are displayed.

Halpern, D.↗
NASA NTRS1993

An atlas of monthly mean distributions of SSMI surface wind speed, ARGOS buoy drift, AVHRR/2 sea surface temperature, and ECMWF surface wind components during 1991

The following monthly mean global distributions for 1991 are presented with a common color scale and geographical map: 10-m height wind speed estimated from the Special Sensor Microwave Imager (SSMI) on a United States Air Force Defense Meteorological Satellite Program (DMSP) spacecraft; sea surface temperature estimated from the advanced very high resolution radiometer (AVHRR/2) on a U.S. National Oceanic and Atmospheric Administration (NOAA) spacecraft; Cartesian components of free-drifting buoys which are tracked by the ARGOS navigation system on NOAA satellites; and Cartesian components of the 10-m height wind vector computed by the European Center for Medium-Range Weather Forecasting (ECMWF). Charts of monthly mean value, sampling distribution, and standard deviation value are displayed. Annual mean distributions are displayed.

Halpern, D.↗
DOE OSTI2016

A Bayesian Machine Learning Model for Estimating Building Occupancy from Open Source Data

Understanding building occupancy is critical to a wide array of applications including natural hazards loss analysis, green building technologies, and population distribution modeling. Due to the expense of directly monitoring buildings, scientists rely in addition on a wide and disparate array of ancillary and open source information including subject matter expertise, survey data, and remote sensing information. These data are fused using data harmonization methods which refer to a loose collection of formal and informal techniques for fusing data together to create viable content for building occupancy estimation. In this paper, we add to the current state of the art by introducing the Population Data Tables (PDT), a Bayesian based informatics system for systematically arranging data and harmonization techniques into a consistent, transparent, knowledge learning framework that retains in the final estimation uncertainty emerging from data, expert judgment, and model parameterization. PDT probabilistically estimates ambient occupancy in units of people/1000ft2 for over 50 building types at the national and sub-national level with the goal of providing global coverage. The challenge of global coverage led to the development of an interdisciplinary geospatial informatics system tool that provides the framework for capturing, storing, and managing open source data, handling subject matter expertise, carrying out Bayesian analytics as well as visualizing and exporting occupancy estimation results. We present the PDT project, situate the work within the larger community, and report on the progress of this multi-year project.Understanding building occupancy is critical to a wide array of applications including natural hazards loss analysis, green building technologies, and population distribution modeling. Due to the expense of directly monitoring buildings, scientists rely in addition on a wide and disparate array of ancillary and open source information including subject matter expertise, survey data, and remote sensing information. These data are fused using data harmonization methods which refer to a loose collection of formal and informal techniques for fusing data together to create viable content for building occupancy estimation. In this paper, we add to the current state of the art by introducing the Population Data Tables (PDT), a Bayesian model and informatics system for systematically arranging data and harmonization techniques into a consistent, transparent, knowledge learning framework that retains in the final estimation uncertainty emerging from data, expert judgment, and model parameterization. PDT probabilistically estimates ambient occupancy in units of people/1000 ft 2 for over 50 building types at the national and sub-national level with the goal of providing global coverage. The challenge of global coverage led to the development of an interdisciplinary geospatial informatics system tool that provides the framework for capturing, storing, and managing open source data, handling subject matter expertise, carrying out Bayesian analytics as well as visualizing and exporting occupancy estimation results. We present the PDT project, situate the work within the larger community, and report on the progress of this multi-year project.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
NASA NTRS1994

An Approximate Kalman Filter for Ocean Data Assimilation; An Example with an Idealized Gulf Stream Model

A practical method of data assimilation for use with large, nonlinear, ocean general circulation models is explored. A Kalman filter based on approximations of the state error covariance matrix is presented, employing a reduction of the effective model dimension, the error's asymptotic steady-state limit, and a time-invariant linearization of the dynamic model for the error integration. The approximations lead to dramatic computational savings in applying estimation theory to large complex systems. We examine the utility of the approximate filter in assimilating different measurement types using a twin experiment of an idealized Gulf Stream. A nonlinear primitive equation model of an unstable east-west jet is studied with a state dimension exceeding 170,000 elements. Assimilation of various pseudo measurements is examined, including velocity, density, and volume transport at localized arrays, and realistic distributions of satellite altimetry and acoustic tomography observations. Results are compared in terms of their effects on the accuracies of the estimation. The approximate filter is shown to outperform a previous study that used an empirical nudging scheme. The examples demonstrate that useful approximate estimation errors can be computed in a practical manner for general circulation models.

['Kalman↗
NASA NTRS2024

Real-Time Temperature Field State Observer and Digital Replica to Support Cyber-Physical Testing

Cyber-physical testing is a class of experimental methods in which a system is partitioned into physical and numerical components to study its overall behavior. Transfer systems are typically needed to capture relevant interactions by enforcing appropriate interface conditions between the physical and numerical portions. While cyber-physical testing is often performed using lumped-parameter systems, thermomechanical cyber-physical testing requires the use of thermal actuators as transfer systems to apply a thermal condition over a spatially distributed region or surface. Thus, the interactions between the numerical and physical components are of a distributed nature, i.e., a field rather than point values. In this paper, we develop and experimentally validate a state observer for temperature fields to enable thermomechanical cyber-physical testing. The state observer provides continuous estimates of temperature over a surface of the thermal actuator using only temperature measurements at discrete locations. The estimated temperature field is then leveraged to account for localized mechanical defects in the interface condition. The method consists of building a digital replica of one portion of the control plant. By imposing defects in the digital replica it is used to generate and output the response of a thermal load pattern that represents that of a physical system with mechanical defects.

Manuel Salmeron↗
DOE OSTI2024

NPFTURBULENCE: Best Estimate Aerosol Size Distribution by airborne measurements

The original data were collected during the field campaign of “Turbulent layers promoting New Particle Formation” experiment (NPFTURBULENCE; https://www.arm.gov/research/campaigns/aaf2024npfturbulence) over the Atmospheric Radiation Measurement (ARM) Southern Great Plains (SGP) Atmospheric Observatory (https://www.arm.gov/capabilities/observatories/sgp ) in north-central Oklahoma. The ARM Aerial Facility ArcticShark uncrewed aerial system (UAS) was based at Blackwell–Tonkawa Municipal Airport (IATA: BWL, ICAO: KBKN, FAA LID: BKN, 36.74475° N, 97.34918° W, 313.9 m MSL), for the field campaign from May 5 through May 29, 2024. The ArcticShark UAS performed 11 flights, including 10 research flights over the Central Facility of the ARM SGP to measure atmospheric state, turbulence, surface IT temperature and imagery, aerosol number concentration and size distribution. The current data set presents Best Estimate Aerosol Size Distribution: a merged aerosol size distribution composed of the data from 2 sensors: miniaturized Scanning Electrical Mobility Sizer (mSEMS) and Portable Optical Particle Spectrometer (POPS). The mSEMS data were interpolated to 1 second from “native” time resolution of about 15 second to match the other probe. The POPS data were converted from equivalent optical size into geometric size using value of aerosol refractive index of 1.477 from the HISCALE field campaign (same geographical area, altitudes, and time of year; http://www.arm.gov/campaigns/aaf2016hiscale ).

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1985

A summary of the UV, optical and IR properties of disks in CV's

The available observations of dwarf novae and novalike systems that are relevant to the study of accretion in cataclysmic variables are reviewed. This includes UV data from IUE and optical spectroscopy and infrared photometry. The flux distributions are compared to available steady state models to obtain an estimate of the mass accretion and the range of this parameter among the dwarf novae (U Gem, Z Cam and SU UMa types) and the novalikes (UX UMa, AM Her, DQ Her types). The normal and abnormal range of variation in the continuum and line fluxes is discussed.

Szkody, P.↗
NASA NTRS2009

Model Calibration in Watershed Hydrology

Hydrologic models use relatively simple mathematical equations to conceptualize and aggregate the complex, spatially distributed, and highly interrelated water, energy, and vegetation processes in a watershed. A consequence of process aggregation is that the model parameters often do not represent directly measurable entities and must, therefore, be estimated using measurements of the system inputs and outputs. During this process, known as model calibration, the parameters are adjusted so that the behavior of the model approximates, as closely and consistently as possible, the observed response of the hydrologic system over some historical period of time. This Chapter reviews the current state-of-the-art of model calibration in watershed hydrology with special emphasis on our own contributions in the last few decades. We discuss the historical background that has led to current perspectives, and review different approaches for manual and automatic single- and multi-objective parameter estimation. In particular, we highlight the recent developments in the calibration of distributed hydrologic models using parameter dimensionality reduction sampling, parameter regularization and parallel computing.

Yilmaz, Koray K.↗
DOE OSTI2014

Statistics of reversible transitions in two-state trajectories in force-ramp spectroscopy

A possible way to extract information about the reversible dissociation of a molecular adhesion bond from force fluctuations observed in force ramp experiments is discussed. For small loading rates the system undergoes a limited number of unbinding and rebinding transitions observable in the so-called force versus extension (FE) curves. The statistics of these transient fluctuations can be utilized to estimate the parameters for the rebinding rate. This is relevant in the experimentally important situation where the direct observation of the reversed FE-curves is hampered, e.g., due to the presence of soft linkers. I generalize the stochastic theory of the kinetics in two-state models to the case of time-dependent kinetic rates and compute the relevant distributions of characteristic forces. While for irreversible systems there is an intrinsic relation between the rupture force distribution and the population of the free-energy well of the bound state, the situation is slightly more complex if reversible systems are considered. For a two-state model, a “stationary” rupture force distribution that is proportional to the population can be defined and allows to consistently discuss quantities averaged over the transient fluctuations. While irreversible systems are best analyzed in the soft spring limit of small pulling device stiffness and large loading rates, here I argue to use the stiffness of the pulling device as a control parameter in addition to the loading rate.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2019

Biogeochemical Equation of State for the Sea-Air Interface

We have recently argued that marine interfacial surface tension must have a distinctive biogeography because it is mediated by fresh surfactant macromolecules released locally through the food web. Here we begin the process of quantification for associated climate flux implications. A low dimensionality (planar) equation of state is invoked at the global scale as our main analysis tool. For the reader’s convenience, fundamental surfactant physical chemistry principles are reviewed first, as they pertain to tangential forces that may alter oceanic eddy, ripple, and bubble fields. A model Prandtl (neutral) wind stress regime is defined for demonstration purposes. It is given the usual dependence on roughness, but then in turn on the tension reduction quantity known as surface pressure. This captures the main net influences of biology and detrital organics on global microlayer physics. Based on well-established surrogate species, tangent pressures are related to distributed ecodynamics as reflected by the current marine systems science knowledge base. Reductions to momentum and related heat-vapor exchange plus gas and salt transfer are estimated and placed on a coarse biogeographic grid. High primary production situations appear to strongly control all types of transfer, whether seasonally or regionally. Classic chemical oceanographic data on boundary state composition and behaviors are well reproduced, and there is a high degree of consistency with conventional micrometeorological wisdom. But although our initial best guesses are quite revealing, coordinated laboratory and field experiments will be required to confirm the broad hypotheses even partially. We note that if the concepts have large scale validity, they are super-Gaian. Biological control over key planetary climate-transfer modes may be accomplished through just a single rapidly renewed organic monolayer.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2011

Contributions of Precipitation and Soil Moisture Observations to the Skill of Soil Moisture Estimates in a Land Data Assimilation System

The contributions of precipitation and soil moisture observations to the skill of soil moisture estimates from a land data assimilation system are assessed. Relative to baseline estimates from the Modern Era Retrospective-analysis for Research and Applications (MERRA), the study investigates soil moisture skill derived from (i) model forcing corrections based on large-scale, gauge- and satellite-based precipitation observations and (ii) assimilation of surface soil moisture retrievals from the Advanced Microwave Scanning Radiometer for the Earth Observing System (AMSR-E). Soil moisture skill is measured against in situ observations in the continental United States at 44 single-profile sites within the Soil Climate Analysis Network (SCAN) for which skillful AMSR-E retrievals are available and at four CalVal watersheds with high-quality distributed sensor networks that measure soil moisture at the scale of land model and satellite estimates. The average skill (in terms of the anomaly time series correlation coefficient R) of AMSR-E retrievals is R=0.39 versus SCAN and R=0.53 versus CalVal measurements. The skill of MERRA surface and root-zone soil moisture is R=0.42 and R=0.46, respectively, versus SCAN measurements, and MERRA surface moisture skill is R=0.56 versus CalVal measurements. Adding information from either precipitation observations or soil moisture retrievals increases surface soil moisture skill levels by IDDeltaR=0.06-0.08, and root zone soil moisture skill levels by DeltaR=0.05-0.07. Adding information from both sources increases surface soil moisture skill levels by DeltaR=0.13, and root zone soil moisture skill by DeltaR=0.11, demonstrating that precipitation corrections and assimilation of satellite soil moisture retrievals contribute similar and largely independent amounts of information.

Reichle, Rolf H.↗
DOE OSTI2018

Resilient distributed state estimation with mobile agents: overcoming Byzantine adversaries, communication losses, and intermittent measurements

Applications in environmental monitoring, surveillance and patrolling typically require a network of mobile agents to collectively gain information regarding the state of a static or dynamical process evolving over a region. However, these networks of mobile agents also introduce various challenges, including intermittent observations of the dynamical process, loss of communication links due to mobility and packet drops, and the potential for malicious or faulty behavior by some of the agents. The main contribution of this paper is the development of resilient, fully-distributed, and provably correct state estimation algorithms that simultaneously account for each of the above considerations, and in turn, offer a general framework for reasoning about state estimation problems in dynamic, failure-prone and adversarial environments. Specifically, we develop a simple switched linear observer for dealing with the issue of time-varying measurement models, and resilient filtering techniques for dealing with worst-case adversarial behavior subject to time-varying communication patterns among the agents. Here, our approach considers both communication patterns that recur in a deterministic manner, and patterns that are induced by random packet drops. For each scenario, we identify conditions on the dynamical system, the patrols, the nominal communication network topology, and the failure models that guarantee applicability of our proposed techniques. Finally, we complement our theoretical results with detailed simulations that illustrate the efficacy of our algorithms in the presence of the technical challenges described above.

42 ENGINEERING↗
DOE OSTI2020

Optimization-Based Fast-Frequency Estimation and Control of Low-Inertia Microgrids

The lack of inertial response from non-synchronous, inverter-based generation in microgrids makes the power system vulnerable to a large rate of change of frequency (ROCOF) and frequency excursions. Energy storage systems (ESSs) can be utilized to provide fast-frequency support to prevent such large excursions in the system. However, fast-frequency support is a power-intensive application that has a significant impact on the ESS lifetime. In this paper, a framework that allows the ESS operator to provide fast-frequency support as a service is proposed. The framework maintains the desired quality-of-service (limiting the ROCOF and frequency) while taking into account the ESS lifetime and physical limits. The framework utilizes moving horizon estimation (MHE) to estimate the frequency deviation and ROCOF from noisy phase-locked loop (PLL) measurements. These estimates are employed by a model predictive control (MPC) algorithm that computes control actions by solving a finite-horizon, online optimization problem. Additionally, this approach avoids oscillatory behavior induced by delays that are common when using low-pass filters as with traditional derivative-based (virtual inertia) controllers. MATLAB/Simulink simulations on a test system from Cordova, Alaska, show the effectiveness of the MHE-MPC approach to reduce frequency deviations and ROCOF of a low-inertia microgrid.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2020

Inferring the Dynamics of the State Evolution During Quantum Annealing

To solve an optimization problem using a commercial quantum annealer, one has to represent the problem of interest as an Ising or a quadratic unconstrained binary optimization (QUBO) problem and submit its coefficients to the annealer, which then returns a user-specified number of low-energy solutions. It would be useful to know what happens in the quantum processor during the anneal process so that one could design better algorithms or suggest improvements to the hardware. However, existing quantum annealers are not able to directly extract such information from the processor. Hence, in this work we propose to use advanced features of D-Wave 2000Q to indirectly infer information about the dynamics of the state evolution during the anneal process. Specifically, D-Wave 2000Q allows the user to customize the anneal schedule, that is, the schedule with which the anneal fraction is changed from the start to the end of the anneal. Furthermore, using this feature, we design a set of modified anneal schedules whose outputs can be used to generate information about the states of the system at user-defined time points during a standard anneal. With this process, called "slicing", we obtain approximate distributions of lowest-energy anneal solutions as the anneal time evolves. We use our technique to obtain a variety of insights into the annealer, such as the state evolution during annealing, when individual bits in an evolving solution flip during the anneal process and when they stabilize, and we introduce a technique to estimate the freeze-out point of both the system as well as of individual qubits.

42 ENGINEERING↗