Search NASA⌕ Search

SEARCH · Search NASA

Results for “horizon matching”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Standardized Data to Improve Understanding and Modeling of Soil Nitrogen at Continental Scale

Abstract Nitrogen (N) is a key limiting nutrient in terrestrial ecosystems, but there remain critical gaps in our ability to predict and model controls on soil N cycling. This may be in part due to lack of standardized sampling across broad spatial–temporal scales. Here, we introduce a continentally distributed, publicly available data set collected by the National Ecological Observatory Network (NEON) that can help fill these gaps. First, we detail the sampling design and methods used to collect and analyze soil inorganic N pool and net flux rate data from 47 terrestrial sites. We address methodological challenges in generating a standardized data set, even for a network using uniform protocols. Then, we evaluate sources of variation within the sampling design and compare measured net N mineralization to simulated fluxes from the Community Earth System Model 2 (CESM2). We observed wide spatiotemporal variation in inorganic N pool sizes and net transformation rates. Site explained the most variation in NEON’s stratified sampling design, followed by plots within sites. Organic horizons had larger pools and net N transformation rates than mineral horizons on a sample weight basis. The majority of sites showed some degree of seasonality in N dynamics, but overall these temporal patterns were not matched by CESM2, leading to poor correspondence between observed and modeled data. Looking forward, these data can reveal new insights into controls on soil N cycling, especially in the context of other environmental data sets provided by NEON, and should be leveraged to improve predictive modeling of the soil N cycle.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Quantum fate of the Choptuik naked singularity

Classical critical collapse provides a dynamical route from smooth initial data to a naked singularity, representing a sharper violation of predictability than ordinary black hole singularities. We argue that this distinction is erased by quantum backreaction. Building on the semiclassical interior analysis, where quantum self-energy of the collapsing matter generates a universal growing mode and a finite mass gap, we study the exterior naked singularity region that determines global visibility in the Einstein-scalar system. We analyze controlled exterior models in both $2+1$ and $3+1$ dimensions. In the former, smooth matching and physical boundary conditions analytically select a vacuum polarization state, whose backreaction cloaks the classically naked region by a quantum trapped branch. In the latter, numerical horizon tracing shows that near a quantum-shifted threshold the exterior develops finite-mass marginally trapped surfaces rather than a zero-mass naked endpoint. These results suggest a global quantum picture in which the Choptuik naked singularity shares the fate of an ordinary black hole singularity: quantum effects push the putative Cauchy horizon behind a quantum-generated horizon, thereby reducing the loss of predictability to the standard black hole evaporation problem.

FOS: Physical sciences↗
NASA NTRS2018

Remote Sensing of the Ocean Surface Refractive Index via Short-Wave Infrared Polarimetry

The linear polarization of sunlight reflected by ocean surfaces in the shortwave infrared (SWIR), at geometries where specular reflection dominates the signal, is a direct function of the refractive index of the surface microlayer(SML). This simple physical concept is at the base of a novel technique presented in this study. We invert observations obtained by the airborne Research Scanning Polarimeter (RSP) in the sun glint region, where each pixel's radiance is dominated by the signal originating from the wave slopes oriented precisely to cause specular reflection. The SWIR wavelength ensures minimization of aerosol interference when radiance travels through the atmosphere; strong absorption by the water body then limits the penetration depth to the first micrometer or so, effectively probing the SML. The resulting Degree of Linear Polarization (DoLP) is then governed by the refractive index via the Fresnel law for the specfic pixel geometry, independently of the wind speed. The selected dataset concerns several field deployments from both low- and high-altitude aircraft, including total reflectance measurements with the sole purpose of accounting for the residual aerosol effect.Stable retrievals from transects above pure seawater yielded values of refractive index that match the values published in the literature within an accuracy of 5 x 10 [superscript 4]. Flying over the oil spill caused by the explosion of the Deepwater Horizon platform, detected variations were found compatible with the presence of an oil slick. The robustness of the results, guaranteedby the high RSP polarimetric accuracy (less than or equal to 0:2%), opens the possibility for remote-sensing detection of other entities that similarly affect the refractive index including whitecaps, microplastics, biological gels, seaweed and grass mats.

Surface Microlayer↗
DOE OSTI2025

Enhancing Autonomous Control of Microreactors Using Multi-Agent Reinforcement Learning

In order for microreactors to be economically competitive, operation costs will need to be minimized through some degree of autonomous control. Previous work has demonstrated the effectiveness of reinforcement learning (RL) for load-following control in a drum-controlled microreactor. This study extends that work by exploring the potential of RL to independently control each of the reactor’s drums. We compare a single-agent RL approach with a multi-agent RL (MARL) framework, testing them for generalization across different load-following power profiles and control timescales, and for robustness in cases of randomly disabled control drums. Since the point kinetics simulation environment used in this study cannot resolve spatial effects, we assume that in the absence of spatially localized disturbances, optimal drum movements should be symmetrical. We demonstrate that single-agent RL is able to achieve accurate performance only when symmetric actions are ignored; otherwise, it fails to train a useful controller. Meanwhile, the MARL framework performs symmetric actions by design and trains a robust, accurate agent, as evidenced by mean absolute errors in power matching of 0.41% for the training power profile, 0.68% for a profile with half the drums disabled, and 0.21% for a profile on a realistic load-following time horizon.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

M-STAR: Magnetism second target advanced reflectometer at the Spallation Neutron Source

M-STAR is a next generation polarized neutron reflectometer with advanced capabilities. A new focusing guide concept is optimized for samples with dimensions down to a millimeter range. A proposed hybrid pulse-skipping chopper will enable experiments at constant geometry at one incident angle in a broad range of wavevector transfer Q up to 0.3 A –1 for specular, off-specular, and GISANS measurements. M-STAR will empower nanoscience and spintronics studies routinely on small samples (~2 × 2 mm 2 ) and of atomic-scale thickness using versatile experimental conditions of magnetic and/or electric fields, light, and temperature applied in situ to novel complex device-like nanosystems with multiple buried interfaces. M-STAR will enable improved grazing incidence diffraction measurements, as a surface-sensitive depth-resolved probe of, e.g., the out-of-plane component of atomic magnetic moments in ferromagnetic, antiferromagnetic, and more complex structures as well as in-plane atomic-scale structures inaccessible with contemporary diffractometry and reflectometry. New horizons will be opened by the development of an option to probe near-surface dynamics with inelastic grazing incidence scattering in the time-of-flight mode. These novel options in combination with ideally matched parameters of the second target station will place M-STAR in the world’s leading position for high resolution polarized reflectometry.

47 OTHER INSTRUMENTATION↗
NASA NTRSDate not supplied

An Updated Secondary Lunar Meteoroid Ejecta Model for Engineering Design

Introduction: The surface of the Moon is constantly being bombarded by a flux of meteoroids of various sizes. Impacts due to these meteoroids produce secondary ejecta material at much lower speeds but with a total mass larger than the original impactor. Details about the secondary ejecta are important for planning missions on the lunar surface. In this work, an updated ejecta model is presented called the Meteoroid Model of Secondary Ejecta (MeMoSeE), to re-place the Apollo-era ejecta model, NASA SP-8013 [1], in the SLS-SPEC 159 Design Specification for Natural Environments (DSNE) [2]. The model pro-duces secondary ejecta flux environments for a user-specified location on the lunar surface, and sorts the incoming secondary flux by angular direction and speed. Methods: MeMoSeE is separated into three parts: the inputs, the conversion step, and the integration of fluxes. Inputs to the model include the primary meteoroid fluxes and the primary near-Earth object (NEO) fluxes. Meteoroid fluxes, both asteroidal and cometary, are calculated using the Meteoroid Engineering Model (MEM3) [3] for different locations on the Moon. For each surface location, an ephemeris is generated using the JPL HORIZONS System [4] that feeds into MEM3. The NEO fluxes are approximated by the high-density population of MEM3 (i.e., only the directionality), where the speed distribution of the NEO fluxes is renormalized to match observations [5]. The regolith properties are used as defined in the DSNE [2]. The conversion step utilizes scaling laws given by Housen & Holsapple 2011 [6] to convert the primary impactor flux to the total mass of secondary ejecta. The ejecta distribution, at the point-of-impact (POI), is separated into a zenith angle and azimuthal angle distribution. The zenith angle distribution follows a beta distribution where the peak depends on the impact altitude angle and the impact azimuth [7]. We employ an ejecta azimuth distribution that is based on Rival & Mandeville 1999 [8] which focuses ejecta in the downstream direction for more oblique impacts. Finally, during the integration step, we sum secondary ejecta number fluxes at a particular region-of-interest (ROI) that originated from many POI locations over the entire surface of the Moon. We keep track of both altitude and azimuth angle bins as well as a range of speed bins, following the igloo gridding as done in MEM3 [9]. The ejecta particle size distribution and density is assumed to be the same as the lunar regolith [10, 2]. Results: The primary fluxes are computed for one Metonic cycle (19 years) for various locations over the lunar surface with a fixed orientation. Both the angular and speed distributions of the primary fluxes are dependent on the latitude and longitude. The over-all primary fluxes show a roughly 13% increase from the eastern limb to the western limb. In general, the speed distribution of the secondary fluxes span from a user-defined minimum speed to the escape speed of the Moon (2.38 km/s), roughly following a power-law relation [6]. Different parts of the speed distribution come from different primary impact locations on the Moon. The secondary ejecta is dominated by the slowest speeds, where these particles originate nearby the ROI. For speeds around 71% of the escape speed, the secondary ejecta originates from locations near the antipodal point. On the other hand, for secondary ejecta speeds that exceed roughly 90% the escape speed, the ejecta particles come from all over the lunar surface to the ROI. Comparing the secondary ejecta fluxes from MeMoSeE with NASA SP-8013 [1], there is a reduction by about 2-3 orders of magnitude for secondary ejecta particles greater than 1 µg. These estimates agree with recent findings from Bjorkman & Christiansen 2019 [11]. The secondary ejecta fluxes are also compared with the primary fluxes, where the ejecta fluxes are roughly an order of magnitude greater than the primary fluxes. References: [1] Cour-Palais, B. G., (1969) NASA SP-8013. [2] NASA SLS-SPEC-159 Rev. H (2020). [3] Moorhead, A. V., et al. (2019) JS&R, 1-17. [4] Giorgini, J. D., (2015) IAUGA, 29, 2256293. [5] Moorhead, A. V., (2020) Memo OSMA/MEO/Lunar-001. [6] Housen, K. R., and Holsapple, K. A., (2011) Icarus, 211(1), 856-875. [7] Gault, D. E., and Wedekind, J. A., (1978) L&PSCP, 9, 3843-3875. [8] Rival, M., and Mandeville, J., (1999) Space Debris, 1(1), 45-57. [9] Moorhead, A. V., (2019) MEM3 User Guide. [10] Carrier, W. D., (1973) The Moon, 6(3-4), 250-263. [11] Bjorkman, M. D. & Christiansen, E. L., (2019) ODC

impacts↗
DOE OSTI2024

Absorptive effects and classical black hole scattering

Abstract We describe an approach to incorporating the physical effects of the absorption of energy by the event horizon of black holes in the scattering amplitudes based post-Minkowskian, point-particle effective description. Absorptive dynamics are incorporated in a model-independent way by coupling the usual point-particle description to an invisible sector of gapless internal degrees-of-freedom. The leading order dynamics of this sector are encoded in the low-energy expansion of a spectral density function obtained by matching an absorption cross section in the ultraviolet description. This information is then recycled using the scattering amplitudes based Kosower-Maybee-O’Connell in-in formalism to calculate the leading absorptive contribution to the impulse and change in rest mass of a Schwarzschild black hole scattering with a second compact body sourcing a massless scalar, electromagnetic or gravitational field. The results obtained are in complete agreement with previous worldline Schwinger-Keldysh calculations and provide an alternative on-shell scattering amplitudes approach to incorporating horizon absorption effects in the gravitational two-body problem.

Physics↗
NASA NTRS1977

Ultraviolet spectroscopy of the zodiacal light at 20-deg elongation

The zodiacal light at 20-deg elongation and 10-deg inclination was observed by rocket ultraviolet spectrometers at 10-15-A resolution in the spectral range 1200-3200 A during an experiment designed to observe comet Kohoutek (1973 XII). The data were obtained above 180 km when scattered horizon light in the startracker caused a loss of tracking on the comet. Airglow emission due to NO and O(+), identified spectroscopically and by its variation with altitude, is significant between 1900 and 2500 A. Longward of 2600 A, the spectrum matches that of the sun, and the derived value of the color ratio, relative to the visible, is 0.90 + or - 0.20. At 1600 A, an upper limit on the zodiacal-light emission of 0.07 R per A or 7 hundred-millionths erg/s per sq cm/sterad per A is obtained.

Feldman, P. D.↗
DOE OSTI2025

Kerr effective black hole geometries in supergravity

We derive the explicit embedding of the effective Kerr spacetimes, which are pertinent to the vanishing of static Love numbers, soft hair descriptions of Kerr black holes, and low-frequency scalar-Kerr scattering amplitudes, as solutions within 𝑁 = 2 supergravity. These spacetimes exhibit a hidden 𝑆⁢𝐿⁡(2,𝑅) × 𝑈⁡(1) or 𝑆⁢𝑂⁡(4,2) symmetry resembling the so called subtracted geometries with 𝑆⁢𝐿⁡(2,𝑅) × 𝑆⁢𝐿⁡(2,𝑅) symmetry, which accurately represent the near-horizon geometry of Kerr black holes and, as we will argue most accurately represents the internal structure of the Kerr black hole. To quantify the differences among the effective Kerr spacetimes, we compare their physical quantities, internal structures, and geodesic equations. Although their thermodynamic properties, including entropy, match those of Kerr, our study uncovers significant differences in the interiors of these effective Kerr solutions. A careful examination of the internal structure of the spacetimes highlights the distinctions between various effective Kerr geometries and their quasinormal spectra.

quantum aspects of black holes↗
DOE OSTI2025

Generalized entropy of gravitational fluctuations

The corrections to holographic entanglement entropy from bulk quantum fields in a classical gravitational background are now well understood. They lead, in particular, to unitary Page curves for evaporating black holes. However, the correct treatment of quantum fluctuations of the metric, including graviton excitations, is a longstanding problem. We provide a gauge-invariant prescription for the generalized entropy of gravitons in anti-de Sitter space in terms of areas and bulk entanglement entropy, generalizing the quantum extremal surface prescription to accommodate fluctuations in the semiclassical spacetime geometry. This task requires a careful treatment of the area operator on the graviton Hilbert space and the definition of a “quantum extremal gauge” in which the extremal surface is unperturbed. It also requires us to determine the correct vacuum modular Hamiltonian for the graviton field, which we fix by requiring that it doesn’t contain a boundary term in extremal gauge. We check our prescription with an explicit computation of the vacuum-subtracted generalized entropy of states containing a graviton in an AdS-Rindler background. Our results exactly match vacuum-subtracted von Neumann entropies for stress-tensor excited states in holographic conformal field theory with d > 2 dimensions. We also use covariant phase space techniques to give a partial proof of our prescription when the entanglement wedge for the background spacetime has a bifurcate Killing horizon. Along the way, we identify a class of perturbative graviton states that have parametrically larger generalized entropy, in the small G N expansion, than any low-energy excitations of an ordinary quantum field.

1/N expansion↗
NASA NTRS2022

Electromagnetic signatures from supermassive binary black holes approaching merger

We present fully relativistic predictions for the electromagnetic emission produced by accretion disks surrounding spinning and nonspinning supermassive binary black holes on the verge of merging. We use the code Bothros to post-process data from 3D General Relativistic Magnetohydrodynamic (GRMHD) simulations via ray-tracing calculations. These simulations model the dynamics of a circumbinary disk and the mini-disks that form around two equal-mass black holes orbiting each other at an initial separation of 20 gravitational radii, and evolve the system for more than 10 orbits in the inspiral regime. We model the emission as the sum of thermal blackbody radiation emitted by an optically thick accretion disk and a power-law spectrum extending to hard X-rays emitted by a hot optically thin corona. We generate time-dependent spectra, images, and light curves at various frequencies to investigate intrinsic periodic signals in the emission, as well as the effects of the black hole spin. We find that prograde black hole spin makes mini-disks brighter since the smaller ISCO angular momentum demands more dissipation before matter plunges to the horizon. However, compared to mini-disks in larger separation binaries with spinning black holes, our mini-disks are less luminous: unlike those systems, their mass accretion rate is lower than in the circumbinary disk, and they radiate with lower efficiency because their inflow times are shorter. Compared to a single black hole system matched in mass and accretion rate, these binaries have spectra noticeably weaker and softer in the UV. Finally, we discuss the implications of our findings for the potential observability of these systems.

Eduardo M. Gutiérrez↗
DOE OSTI2020

AdS 6 /CFT 5 with O7-planes

Type IIB AdS6 solutions with orientifold 7-planes are constructed. The geometry is a warped product of AdS 6 and S 2 over a Riemann surface Σ and the O7-planes correspond to a particular type of puncture on Σ. The solutions are identified as near-horizon limits of (p,q) 5-brane webs with O7-planes. The dual 5d SCFTs have relevant deformations to linear quiver gauge theories which have SO(•) or USp(•) nodes or SU(•) nodes with hypermultiplets in symmetric or antisymmetric representations, in addition to SU(•) nodes with fundamental hypermultiplets. The S 5 free energies are obtained holographically and matched to field theory computations using supersymmetric localization to support the proposed dualities.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

From microscopic black hole entropy toward Hawking radiationc

The AdS/CFT correspondence has recently provided a novel approach for counting the microstates of black holes impressively matching the macroscopic Bekenstein–Hawking entropy formula of rotating electrically charged asymptotically AdS black holes in four to seven dimensions. This approach is designed for supersymmetric extremal black holes, but can also be extended to nonsupersymmetric, near-extremal black holes. Besides the dual higher-dimensional boundary CFT, an effective 2D CFT emerges in a certain near-horizon limit accounting for both the extremal and the near-extremal black hole entropies. Furthermore, this effective 2D description is universal across dimensions and comes naturally equipped with an approach to quantitatively tackle aspects of Hawking radiation.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1998

Comparison of MPEG digital video with super VHS tape for diagnostic echocardiographic readings

BACKGROUND: Digital recording of echocardiographic studies is on the clinical horizon. However, full digital capture of complete echocardiographic studies in traditional video format is impractical, given current storage capacity and network bandwidth. To overcome these constraints, we evaluated the diagnostic image quality of digital video by using MPEG (Motion Picture Experts Group) compression. METHODS AND RESULTS: Fifty-eight complete, consecutive studies were recorded simultaneously with the use of MPEG-1 and sVHS videotape. Each matched MPEG and sVHS study pair was reviewed by two from a total of six readers, and findings were recorded with the use of a detailed, computerized reporting tool. Intrareader and interreader discrepancies were characterized as major or minor and analyzed in total and for specific subgroups of findings (left and right ventricular parameters, valvular insufficiency, and left ventricular regional wall motion). Intrareader discrepancies were reviewed by a consensus panel for agreement with either MPEG or sVHS findings. There was an exact concordance between MPEG and sVHS readings in 83% of findings. The majority of discrepancies were minor, with major discrepancies in only 2.7% of findings. There was no difference in the rate of consensus panel agreement with MPEG or sVHS for instances of intrareader discrepancy, either in total or for any subgroup of findings. Interreader discrepancy rates were nearly identical for both MPEG and sVHS. CONCLUSIONS: MPEG-1 digital video is equivalent to sVHS videotape for diagnostic echocardiography. MPEG increases the range of practical options for digital echocardiography and offers, for the first time, the advantages of digital recording in a familiar video format.

Non-NASA Center↗
NASA NTRS2020

Integrated Demand Management: Concepts and Procedures

This report provides a comprehensive description of the Integrated Demand Management concept. Motivation: NASA’s Integrated Demand Management (IDM) research explores the idea that, under certain conditions, time-based flow management (TBFM) arrival operations can benefit from the coordinated use of a strategic traffic management initiative (TMI) to “precondition” the inbound demand. The research was motivated by the observation that TBFM was usually turned off during convective weather, even in facilities where it was routinely used. Our hypothesis was that strategic adjustments to the inbound traffic so that it provided a better match to the off-nominal changes in capacity observed during these conditions could enable TBFM scheduling to continue to provide effective support for arrival traffic management. Concept: IDM proposes that a TMI (e.g., a Collaborative Trajectory Options Program, or CTOP) be used to adjust the rate and/or geographic distribution across flows of the traffic inbound to a high-demand, TBFM-managed airport before that traffic reaches the TBFM planning horizon. After this strategic preconditioning, TBFM can then tactically fine-tune the demand to deliver a well-managed, orderly feed to the destination airport. Coordinated use of these two flow management capabilities is intended to improve system performance in terms of: • Equity of ground delay assignment, avoiding excessive ground delay for TBFM-scheduled departures, without penalizing longer flights; • Throughput, by distributing traffic to maximize use of available capacity; • Predictability for operators, providing advance notice about the impact on individual flights; • Increased flexibility, supporting operator mitigation strategies such as slot swapping or trajectory options; • Efficiency of flight operations, using ground delay more effectively and reducing airborne delay. The operational description in this document highlights how the IDM concept builds upon already existing tools and procedures, and also indicates where tool enhancements could facilitate conduct of IDM operations. However, enhanced tools are not a requirement for concept introduction. In fact, initial deployment that focused on training procedures and rationale for coordinated use of TFMS and TBFM, without changes to existing tools, might be a simpler way to introduce and to familiarize traffic managers with the idea of preconditioning. The concept and procedures described in this document can hopefully provide useful guidance for introduction of IDM into field operations..

IDM↗
NASA NTRS1987

Environmental modeling and recognition for an autonomous land vehicle

An architecture for object modeling and recognition for an autonomous land vehicle is presented. Examples of objects of interest include terrain features, fields, roads, horizon features, trees, etc. The architecture is organized around a set of data bases for generic object models and perceptual structures, temporary memory for the instantiation of object and relational hypotheses, and a long term memory for storing stable hypotheses that are affixed to the terrain representation. Multiple inference processes operate over these databases. Researchers describe these particular components: the perceptual structure database, the grouping processes that operate over this, schemas, and the long term terrain database. A processing example that matches predictions from the long term terrain model to imagery, extracts significant perceptual structures for consideration as potential landmarks, and extracts a relational structure to update the long term terrain database is given.

Lawton, D. T.↗
NASA NTRS2014

The Nustar Spectrum of Mrk 335: Extreme Relativistic Effects Within Two Gravitational Radii of the Event Horizon?

We present 3-50 keV NuSTAR observations of the active galactic nuclei Mrk 335 in a very low flux state. The spectrum is dominated by very strong features at the energies of the iron line at 5-7 keV and Compton hump from 10-30 keV. The source is variable during the observation, with the variability concentrated at low energies, which suggesting either a relativistic reflection or a variable absorption scenario. In this work, we focus on the reflection interpretation, making use of new relativistic reflection models that self consistently calculate the reflection fraction, relativistic blurring and angle-dependent reflection spectrum for different coronal heights to model the spectra. We find that the spectra can be well fitted with relativistic reflection, and that the lowest flux state spectrum is described by reflection alone, suggesting the effects of extreme light-bending occurring within approx. 2 gravitational radii (RG) of the event horizon. The reflection fraction decreases sharply with increasing flux, consistent with a point source moving up to above 10 RG as the source brightens. We constrain the spin parameter to greater than 0.9 at the 3(sigma) confidence level. By adding a spin-dependent upper limit on the reflection fraction to our models, we demonstrate that this can be a powerful way of constraining the spin parameter, particularly in reflection dominated states. We also calculate a detailed emissivity profile for the iron line, and find that it closely matches theoretical predictions for a compact source within a few RG of the black hole.

galaxies: active aEuro" galaxies: individual: Mrk ↗
DOE OSTI2022

Model predictive control for robust quantum state preparation

A critical engineering challenge in quantum technology is the accurate control of quantum dynamics. Model-based methods for optimal control have been shown to be highly effective when theory and experiment closely match. Consequently, realizing high-fidelity quantum processes with model-based control requires careful device characterization. In quantum processors based on cold atoms, the Hamiltonian can be well-characterized. For superconducting qubits operating at milli-Kelvin temperatures, the Hamiltonian is not as well-characterized. Unaccounted for physics (i.e., mode discrepancy), coherent disturbances, and increased noise compromise traditional model-based control. This work introduces model predictive control (MPC) for quantum control applications. MPC is a closed-loop optimization framework that (i) inherits a natural degree of disturbance rejection by incorporating measurement feedback, (ii) utilizes finite-horizon model-based optimizations to control complex multi-input, multi-output dynamical systems under state and input constraints, and (iii) is flexible enough to develop synergistically alongside other modern control strategies. We show how MPC can be used to generate practical optimized control sequences in representative examples of quantum state preparation. Specifically, we demonstrate for a qubit, a weakly-anharmonic qubit, and a system undergoing crosstalk, that MPC can realize successful model-based control even when the model is inadequate. These examples showcase why MPC is an important addition to the quantum engineering control suite.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗