Search NASA⌕ Search

SEARCH · Search NASA

Results for “determining the coordination number”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Effects of Al:Si and (Al+Na):Si Ratios on the Properties of the International Simple Glass, Part II: Structure

High-alumina containing high-level waste (HLW) will be vitrified at the Waste Treatment Plant at the Hanford Site. The resulting glasses, high in alumina, will have distinct composition-structure-property (C-S-P) relationships compared to previously studied HLW glasses. These C-S-P relationships determine the processability and product durability of glasses and therefore must be understood. The main purpose of this study is to understand the detailed structural changes caused by Al:Si and (Al+Na):Si substitutions in a simplified nuclear waste model glass (ISG, International Simple Glass) by combining experimental structural characterizations and molecular dynamics (MD) simulations. The structures of these two series of glasses were characterized by neutron total scattering and 27Al, 23Na, 29Si, and 11B solid-state nuclear magnetic resonance (NMR) spectroscopy. Additionally, MD simulations were used to generate atomistic structural models of the borosilicate glasses and simulation results were validated by the experimental structural data. Short-range (e.g., bond distance, coordination number, etc.) and medium-range (e.g., oxygen speciation, network connectivity, polyhedral linkages) structural features of the borosilicate glasses were systematically investigated as a function of the degree of substitution. The results show that bond distance and coordination number of the cation-oxygen pairs are relatively insensitive to Al:Si and (Al+Na):Si substitutions with the exception of the B-O pair. Additionally, the Al:Si substitution results in an increase of tri-bridging oxygen species, while (Al+Na):Si substitution creates non-bridging oxygen species. Charge compensator preferences were found for Si-[NBO] (Na+), [3]B-[NBO] (Na+), [4]B (mostly Ca2+), [4]Al (nearly equally split Na+ and Ca2+), and [6]Zr (mostly Ca2+). The network former-BO-network former linkages preferences were also tabulated; Si-O-Al and Al-O-Al were preferred at the expense of lower Si-O-[3]B and [3]B-O-[3]B linkages. These results provide insights on the structural origins of property changes such as glass transition temperature caused by the substitutions, providing a basis for future improvements of theoretical and computer simulation models.

Lu, Xiaonan↗
DOE OSTI2022

Z-Contrast Enhancement in Au–Pt Nanocatalysts by Correlative X-ray Absorption Spectroscopy and Electron Microscopy: Implications for Composition Determination

The properties of bimetallic nanoparticles (BNPs) vary widely as a function of their composition and size distributions. X-ray absorption fine structure analysis is commonly used to characterize their structure, but its application to elements that are close to each other in the periodic table is hampered by poor Z-contrast. For this study, we trained an artificial neural network to recognize the partial coordination numbers in AuPt NPs synthesized via peptide templating using their X-ray absorption near-edge structure spectra. This approach, combined with scanning transmission electron microscopy analysis, revealed unique details of this prototype catalytic system that has different forms of heterogeneities at both the intra- and inter-particle levels. Our method based on the enhancement of Z-contrast of metal species will have implications for compositional studies of BNPs.

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

Active learning of ternary alloy structures and energies

Abstract Machine learning models with uncertainty quantification have recently emerged as attractive tools to accelerate the navigation of catalyst design spaces in a data-efficient manner. Here, we combine active learning with a dropout graph convolutional network (dGCN) as a surrogate model to explore the complex materials space of high-entropy alloys (HEAs). We train the dGCN on the formation energies of disordered binary alloy structures in the Pd-Pt-Sn ternary alloy system and improve predictions on ternary structures by performing reduced optimization of the formation free energy, the target property that determines HEA stability, over ensembles of ternary structures constructed based on two coordinate systems: (a) a physics-informed ternary composition space, and (b) data-driven coordinates discovered by the Diffusion Maps manifold learning scheme. Both reduced optimization techniques improve predictions of the formation free energy in the ternary alloy space with a significantly reduced number of DFT calculations compared to a high-fidelity model. The physics-based scheme converges to the target property in a manner akin to a depth-first strategy, whereas the data-driven scheme appears more akin to a breadth-first approach. Both sampling schemes, coupled with our acquisition function, successfully exploit a database of DFT-calculated binary alloy structures and energies, augmented with a relatively small number of ternary alloy calculations, to identify stable ternary HEA compositions and structures. This generalized framework can be extended to incorporate more complex bulk and surface structural motifs, and the results demonstrate that significant dimensionality reduction is possible in thermodynamic sampling problems when suitable active learning schemes are employed.

Chemistry↗
DOE OSTI2022

Radar Retrieval Evaluation and Investigation of Dendritic Growth Layer Polarimetric Signatures in a Winter Storm

Abstract This study evaluates ice particle size distribution and aspect ratio φ Multi-Radar Multi-Sensor (MRMS) dual-polarization radar retrievals through a direct comparison with two legs of observational aircraft data obtained during a winter storm case from the Investigation of Microphysics and Precipitation for Atlantic Coast-Threatening Snowstorms (IMPACTS) campaign. In situ cloud probes, satellite, and MRMS observations illustrate that the often-observed K dp and Z DR enhancement regions in the dendritic growth layer can either indicate a local number concentration increase of dry ice particles or the presence of ice particles mixed with a significant number of supercooled liquid droplets. Relative to in situ measurements, MRMS retrievals on average underestimated mean volume diameters by 50% and overestimated number concentrations by over 100%. IWC retrievals using Z DR and K dp within the dendritic growth layer were minimally biased relative to in situ calculations where retrievals yielded −2% median relative error for the entire aircraft leg. Incorporating φ retrievals decreased both the magnitude and spread of polarimetric retrievals below the dendritic growth layer. While φ radar retrievals suggest that observed dendritic growth layer particles were nonspherical (0.1 ≤ φ ≤ 0.2), in situ projected aspect ratios, idealized numerical simulations, and habit classifications from cloud probe images suggest that the population mean φ was generally much higher. Coordinated aircraft radar reflectivity with in situ observations suggests that the MRMS systematically underestimated reflectivity and could not resolve local peaks in mean volume diameter sizes. These results highlight the need to consider particle assumptions and radar limitations when performing retrievals. significance statement Developing snow is often detectable using weather radars. Meteorologists combine these radar measurements with mathematical equations to study how snow forms in order to determine how much snow will fall. This study evaluates current methods for estimating the total number and mass, sizes, and shapes of snowflakes from radar using images of individual snowflakes taken during two aircraft legs. Radar estimates of snowflake properties were most consistent with aircraft data inside regions with prominent radar signatures. However, radar estimates of snowflake shapes were not consistent with observed shapes estimated from the snowflake images. Although additional research is needed, these results bolster understanding of snow-growth physics and uncertainties between radar measurements and snow production that can improve future snowfall forecasting.

Meteorology & Atmospheric Sciences↗
DOE OSTI2026

The Stability of Nd Hydroxyl Complexes at Near Neutral to Alkaline pH from 25 to 75 °C: Implications for Monazite Solubility in Hydrothermal Aqueous Fluids

Speciation of rare earth elements (REE) in aqueous fluids controls their mobilization during fluid-rock interaction. Thermodynamic modeling provides important insights into the factors controlling REE mobility and fractionation in aqueous fluids, in which aqueous REE hydroxyl complexes are particularly important across a broad pH range. However, the stability of these REE hydroxyl complexes remains poorly determined due to limited experimental data above ambient temperature. In this work, UV–vis spectrophotometry experiments were conducted to determine the cumulative hydrolysis constants for neodymium (Nd) from 25 to 75 °C in near-neutral to alkaline solutions. The color indicator m -cresol purple was used to determine in situ pH. Alkaline NaOH-bearing aqueous solutions were doped with varying initial Nd concentrations (0 to ∼0.155 mmol/kg), which resulted in the release of protons (H + ) and pH decrease during hydrolysis of Nd 3+ according to Nd 3+ + n H 2 O = Nd(OH) n 3– n + H + . The average OH – ligand number coordinated to Nd 3+ was found to increase from 1.0 to 1.6 at 25 °C and from 1.2 to 2.8 at 75 °C over a pH range from 6.3 to 9.0. The measured speciation shows an increased predominance of Nd(OH) 3 0 over the Nd(OH) 2 + and Nd(OH) 2+ species with increased temperature and pH. The increasing cumulative formation constants (log*β° n , n = 1 to 3) retrieved from 25 to 75 °C differ by 0.1 to 1.9 logarithmic units in comparison to existing thermodynamic databases. These updated thermodynamic data have important implications for geochemical modeling of the speciation of Nd hydroxyl complexes and the solubility of monazite as a function of pH.

58 GEOSCIENCES↗
DOE OSTI2023

Spatio-Temporal Surrogates for Interaction of a Jet with High Explosives: Part II - Clustering Extremely High-Dimensional Grid-Based Data

Building an accurate surrogate model for the spatio-temporal outputs of a computer simulation is a challenging task. A simple approach to improve the accuracy of the surrogate is to cluster the outputs based on similarity and build a separate surrogate model for each cluster. This clustering is relatively straightforward when the output at each time step is of moderate size. However, when the spatial domain is represented by a large number of grid points, numbering in the millions, the clustering of the data becomes more challenging. In this report, we consider output data from simulations of a jet interacting with high explosives. These data are available on spatial domains of different sizes, at grid points that vary in their spatial coordinates, and in a format that distributes the output across multiple files at each time step of the simulation. We first describe how we bring these data into a consistent format prior to clustering. Borrowing the idea of random projections from data mining, we reduce the dimension of our data by a factor of thousand, making it possible to use the iterative k-means method for clustering. We show how we can use the randomness of both the random projections, and the choice of initial centroids in k-means clustering, to determine the number of clusters in our data set. Our approach makes clustering of extremely high dimensional data tractable, generating meaningful cluster assignments for our problem, despite the approximation introduced in the random projections.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Quantifying a partial polyamorphic transition in a cerium-based metallic glass during cooling

Cerium-based metallic glasses are prototype polyamorphous systems with pressure-induced polyamorphic transitions extensively reported. Cooling typically has a similar effect on materials as compression with regard to reducing volume. However, previous studies show dramatically different behavior of Ce-based metallic glasses between cooling and compression, whose origin remains unclear. Here, using in situ low-temperature synchrotron high-energy x-ray diffraction, the structural evolution of a Ce68Al10Cu20Co2 metallic glass is accurately determined and analyzed by a structure factor and a reduced pair distribution function (PDF) during cooling from 298 to 83 K. An unusually large linear thermal expansion coefficient is revealed, which is associated with both continuous but inconsistent structural changes between the two subpeaks of the first atomic shell in terms of average bond lengths and coordination numbers. These phenomena are suggested to be attributed to a gradual 4f electron delocalization of only a minimal amount (∼2.6% at 83 K) of Ce atoms by quantitative analysis of the PDF data. However, a previously expected global polymorphic transition from a low-density amorphous state to a high-density amorphous state with an abrupt volume collapse is not observed. Moreover, electrical resistivity also shows a continuous increase during cooling without any sharp change. It is clarified that cryogenic temperatures could facilitate but are not powerful enough alone to trigger a global polymorphic transition in the Ce68Al10Cu20Co2 metallic glass, suggesting a wide distribution of its local atomic environment.

Chen, Zhi↗
DOE OSTI2021

Distributed directory of named data elements in coordination namespace

An approach is described that provides a distributed directory structure within a storage of an information handling system (a local node). A request is received with the request corresponding to a shared virtual address. The shared virtual address that is shared amongst a number of nodes that includes the local node and some remote nodes. A Global Address Space Directory (GASD) is retrieved that corresponds to a global virtual address space. The GASD is stored in a Coordination Namespace that is stored in a memory that is distributed amongst the nodes. A mapping that is included in the GASD is used to determine the node where the shared virtual address currently resides. The shared virtual address is then accessed from the node where it currently resides.

Johns, Charles Ray↗
DOE OSTI2021

Molecular dynamics simulation of the effect of cooling rate on the structure and properties of lithium disilicate glass

The effects of cooling rate on the structure and properties of lithium disilicate (LS2) glass are investigated using molecular dynamics (MD) computer simulations. The evolution of structural features such as pair distribution function, bond angle distribution, and Li coordination distribution are determined, and correlated with dynamic and static properties to elucidate the effects of cooling rate. The density, elastic moduli, and diffusion coefficient are found to be highly sensitive to cooling rate, whereas the Si pair distribution function and bond angle distribution are weakly affected by the cooling rate. Additionally, the changes of Si-O-Si bond angle and Li coordination number suggest the formation of Li cluster at lower cooling rates. Furthermore, by comparing results from other simulations and reported experiments, we confirm that the increase of cooling rate leads to an increase of conductivity and a decrease of density. Finally, at very high cooling rates, we find that all the atoms in LS2-glass do not relax simultaneously, but do so in two distinct configurations.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Structure–Activity Relationships for Ethanol Dehydrogenation to Acetaldehyde by Silica-Supported Zinc Oxide Catalysts

Silica-supported ZnO efficiently catalyzes the nonoxidative dehydrogenation of ethanol to acetaldehyde, which is relevant for production of 1,3-butadiene from bioethanol. Characterization with in situ spectroscopies under dehydrated conditions (high sensitivity-low energy ion scattering (HS-LEIS), diffuse reflectance (DR) UV–vis, X-ray absorption spectroscopy (XAS), diffuse reflectance Fourier transform infrared spectroscopy (DRIFTS), inelastic neutron scattering (INS), and UV Raman), and ammonia adsorption probed by temperature-programmed desorption followed by DRIFTS and mass spectrometry (DRIFTS-MS NH 3 -TPD), and DFT calculations revealed that the supported ZnO x phase was present as isolated surface ZnO x sites on SiO 2 , with the vast majority coordinated by two siloxane bonds and one silicon atom with two nonbridging oxygens ((≡SiO) 2 Zn 2+ O 2 Si=), anchored at 4-, 5-, and 6-membered siloxane rings. A minor fraction of surface ZnO x sites possessed Lewis acidity, and even fewer sites possessed a Bro̷nsted acidic Zn(OH) + Si moiety. Ethanol temperature-programmed surface reaction-mass spectrometry (TPSR-MS) with various oxidative or ethanol reaction pretreatments indicated that only sites with Lewis and Bro̷nsted acidic character (Zn(OH) + Si) were active for ethanol dehydrogenation, while the majority surface (≡SiO) 2 Zn 2+ O 2 Si= sites were inactive. Greater heterogeneity among all surface ZnO x sites, as assessed by in situ DR UV–vis spectroscopy, was associated with a greater number of ZnO x sites that were active for ethanol dehydrogenation as well as lower enthalpic barriers for acetaldehyde production among the most active surface ZnO x sites. Turnover frequencies and the apparent activation energy for ethanol dehydrogenation were determined from steady-state kinetics. Together, these findings suggested that anchoring inactive surface (≡SiO) 2 Zn 2+ O 2 Si= sites on the silica support caused a greater number of active surface ZnO x sites to adopt a more strained configuration, promoting ethanol dehydrogenation catalysis. Pretreatments and catalysts that promoted desorption of ethanol during TPSR, taken as a marker of surface dehydroxylation, were associated with an increased number of the most active surface (Zn(OH) + Si) sites. Such findings suggested that inactive surface ZnO x sites were activated for ethanol dehydrogenation by dehydroxylation of the support and/or decreased coordination to hemilabile siloxane ligands.

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

Ionic Contraction across the Lanthanide Series Decreases the Temperature-Induced Disorder of the Water Coordination Sphere

In liquid, temperature affects the structures of lanthanide complexes in multiple ways that depend upon complex interactions between ligands, anions, and solvent molecules. The relative simplicity of lanthanide aqua ions (Ln 3+ ) make them well suited to determine how temperature induces structural changes in lanthanide complexes. In this work, we performed a combination of ab initio molecular dynamics (AIMD) simulations and extended X-ray absorption fine structure (EXAFS) measurements, both at 25 and 90 °C, to determine how temperature affects the first- and second-coordination spheres of three Ln 3+ (Ce 3+ , Sm 3+ , and Lu 3+ ) aqua ions. AIMD simulations show first lanthanide coordination spheres that are similar at 25 and 90 °C, more so for the Lu 3+ ion that remains as eight-coordinate than for the Ce 3+ and Sm 3+ ions that change their preferred coordination number from nine (at 25 °C) to eight (at 90 °C). The measured EXAFS spectra are very similar at 25 and 90 °C, for the Ce 3+ , Sm 3+ , and Lu 3+ ions, suggesting that the dynamical disorder of the Ln 3+ ions in liquid water is sufficient such that temperature-induced changes do not clearly manifest changes in the structure of the three ions. Both AIMD simulations and EXAFS measurements show very similar structures of the first coordination sphere of the Lu 3+ ion at 25 and 90 °C.

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

Determination of the structure of lithium niobosilicate glasses by molecular dynamics simulation with a new Nb-O potential

A partial charge empirical potential for the Nb-O pair has been developed based on existing parameters of other pairs in a Teter potential to enable molecular dynamics simulations of both lithium niobate crystal and lithium niobosilicate glasses. The developed potential is capable of describing structural features of lithium niobosilicate glasses in a wide composition range, including niobium coordination number and bond length, density, bond angle distribution, polyhedral distribution and fraction of bridging oxygen. Furthermore, the results obtained using this new potential show good agreement with experimental data of density and structure from Raman spectroscopy of lithium niobosilicate glasses and X-ray absorption spectroscopy near Nb K-edge.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Cooperative effect of local active stresses on the macroscopic contractility of elastic fiber networks

The collective action of actively contractile units embedded in elastic biopolymer networks plays a crucial role in regulating the network's macroscopic mechanical response. Here, in this study, we investigate how the macroscopic boundary stress in model elastic fiber networks depends on the number and nature of embedded contractile units, each exerting an isotropic force dipole, as well as on the bending stiffness of fibers. We find that the macroscopic stress increases nonlinearly with the number of dipoles due to mutual stiffening of initially soft, bending-dominated networks. Using effective medium theory, we relate this enhanced contractility to an increase in the effective average network coordination number due to constraints imposed by the force dipoles. By comparing three distinct force dipole models that differ in their local structures, we demonstrate that the specific manner in which an active unit constrains the network strongly influences the onset and nature of the stiffening transition. Our results highlight that not only the quantity but also the local geometry of force-generating units critically determines the macroscopic mechanical behavior. This framework provides a physical basis for understanding how biological systems—such as molecular motors in the cytoskeleton, or adherent cells in the extracellular matrix—can modulate network-scale nonlinear elastic properties through local tuning of active force-generating units.

Biological and medical sciences↗
DOE OSTI2023

Structurally Constrained Evolutionary Algorithm for the Discovery and Design of Metastable Phases

Metastable materials are abundant in nature and technology, showcasing remarkable properties that inspire innovative materials design. However, traditional crystal structure prediction methods, which rely solely on energetic factors to determine a structure’s fitness, are not suitable for predicting the vast number of potentially synthesizable phases that represent a local minimum corresponding to a state in thermodynamic equilibrium. Here, we present a new approach for the prediction of metastable phases with specific structural features, and interface this method with the XTALOPT evolutionary algorithm. Our method relies on structural features that include the local crystalline order (e.g., the coordination number or chemical environment), and symmetry (e.g., Bravais lattice and space group) to filter the breeding pool of an evolutionary crystal structure search. The effectiveness of this approach is benchmarked on three known metastable systems: XeN 8 , with a two-dimensional polymeric nitrogen sublattice, brookite TiO 2 , and a high pressure BaH 4 phase that was recently characterized. Additionally, a newly predicted metastable melaminate salt, P1¯WC 3 N 6 , was found to possess an energy that is lower than two phases proposed in a recent computational study. Here, the method presented here could help in identifying the structures of compounds that have already been synthesized, and developing new synthesis targets with desired properties.

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

Kβ X-ray Emission Spectroscopy as a Probe of Cu(I) Sites: Application to the Cu(I) Site in Preprocessed Galactose Oxidase

Cu(I) active sites in metalloproteins are involved in O 2 activation, but their O 2 reactivity is difficult to study due to the Cu(I) d 10 closed shell which precludes the use of conventional spectroscopic methods. Kβ X-ray emission spectroscopy (XES) is a promising technique for investigating Cu(I) sites as it detects photons emitted by electronic transitions from occupied orbitals. Here, we demonstrate the utility of Kβ XES in probing Cu(I) sites in model complexes and a metalloprotein. Using Cu(I)Cl, emission features from double-ionization (DI) states are identified using varying incident X-ray photon energies, and a reasonable method to correct the data to remove DI contributions is presented. Kβ XES spectra of Cu(I) model complexes, having biologically relevant N/S ligands and different coordination numbers, are compared and analyzed, with the aid of density functional theory (DFT) calculations, to evaluate the sensitivity of the spectral features to the ligand environment. While the low-energy Kβ 2,5 emission feature reflects the ionization energy of ligand np valence orbitals, the high-energy Kβ 2,5 emission feature corresponds to transitions from molecular orbitals (MOs) having mainly Cu 3d character with the intensities determined by ligand-mediated d–p mixing. A Kβ XES spectrum of the Cu(I) site in preprocessed galactose oxidase (GO pre ) supports the 1Tyr/2His structural model that was determined by our previous X-ray absorption spectroscopy and DFT study. Finally, the high-energy Kβ 2,5 emission feature in the Cu(I)-GO pre data has information about the MO containing mostly Cu 3d x 2 –y 2 character that is the frontier molecular orbital (FMO) for O 2 activation, which shows the potential of Kβ XES in probing the Cu(I) FMO associated with small-molecule activation in metalloproteins.

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

GPS Spoofing Mitigation and Timing Risk Analysis in Networked Phasor Measurement Units via Stochastic Reachability

To address phasor measurement unit (PMU) vulnerability to spoofing, we propose the use of a set-valued state estimation technique known as stochastic reachability (SR)-based distributed Kalman filter (DKF) that computes secure global positioning system (GPS) timing across a network of receivers. Utilizing SR, we estimate not only GPS time but also its stochastic reachable set, which is parameterized by probabilistic zonotope (p-Zonotope). While requiring known measurement error bounds in only non-spoofed conditions, we designed a two-tiered approach. We first performed measurement-level spoofing mitigation via deviation of a measurement innovation from its expected p-Zonotope. We then performed state-level timing risk analysis via a determination of the intersection probability of the estimated p-Zonotope with an unsafe set that violates IEEE C37.118.1a-2014 standards. Finally, we validated our SR-DKF algorithm by subjecting it to a simulated receiver network to coordinate signal-level spoofing. We demonstrate improved timing accuracy and successful spoofing mitigation via the use of our SR-DKF algorithm. We also validated the robustness of the estimated timing risk as the number of receivers were varied.

47 OTHER INSTRUMENTATION↗
DOE OSTI2020

A Quantum Analog Coprocessor for Correlated Electron Systems Simulation

Analog quantum simulation is an approach for studying physical systems that might otherwise be computationally intractable to simulate on classical high-performance computing (HPC) systems. The key idea behind analog quantum simulation is the realization of a physical system with a low-energy effective Hamiltonian that is the same as the low-energy effective Hamiltonian of some target system to be studied. Purpose-built nanoelectronic devices are a natural candidate for implementing the analog quantum simulation of strongly correlated materials that are otherwise challenging to study using classical HPC systems. However, realizing devices that are sufficiently large to study the properties of a non-trivial material system (e.g., those described by a Fermi-Hubbard model) will eventually require the fabrication, control, and measurement of at least 0(10) quantum dots, or other engineered quantum impurities. As a step toward large-scale analog or digital quantum simulation platforms based on nanoelectronic devices, we propose a new approach to analog quantum simulation that makes use of the large Hilbert space dimension of the electronic baths that are used to adjust the occupancy of one or a few engineered quantum impurities. This approach to analog quantum simulation allows us to study a wide array of quantum impurity models. We can further augment the computational power of such an approach by combining it with a classical computer to facilitate dynamical mean-field theory (DMFT) calculations. DMFT replaces the solution of a lattice impurity problem with the solution of a family of localized impurity problems with bath couplings that are adjusted to satisfy a self-consistency condition between the two models. In DMFT, the computationally challenging task is the high-accuracy solution of an instance of a quantum impurity model that is determined self-consistently in coordination with a mean-field calculation. We propose using one or a few engineered quantum impurities with adjustable couplings to baths to realize an analog quantum coprocessor that effects the solution of such a model through measurements of a physical quantum impurity, operating in coordination with a classical computer to achieve a self-consistent solution to a DMFT calculation. We focus on implementation details relevant to a number of technologies for which Sandia has design, fabrication, and measurement expertise. The primary technical advances outlined in this report concern the development of a supporting modeling capability. As with all analog quantum simulation platforms, the successful design and operation of individual devices depends critically on one's ability to predict the effective low-energy Hamiltonian governing its dynamics Our project has made this possible and lays the foundation for future experimental implementations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Effect of outer-cylinder rotation on the radially heated Taylor–Couette flow

A Taylor–Couette setup with radial heating is considered where a Boussinesq fluid is sheared in the annular region between two concentric, independently rotating cylinders maintained at different temperatures. Linear stability analysis is performed to determine the Taylor number for the onset of instability. Two radius ratios corresponding to wide and thin gaps with several rotation rate ratios are considered. The rotation of the outer cylinder is found to have a general stabilizing effect on the stability threshold as compared to pure inner-cylinder rotation, with a few exceptions. The radial heating sets up an axial flow which breaks the reflection symmetry of isothermal Taylor–Couette flow in the axial coordinate. This symmetry breaking separates linear stability thresholds, and we find the fastest growing modes with both positive and negative azimuthal numbers for different parameters. Another important finding of the current study is the discovery of unstable modes in the Rayleigh-stable regime. Furthermore, closed disconnected neutral curves (CDNCs) are observed for both wide and thin gaps which can separate from or merge into open neutral–stability curves. Alternatively, CDNCs can also morph into open neutral stability curves as the rotation rate ratio is changed. CDNCs are observed to be sensitive to changes in control parameters, and their appearance/disappearance is shown to induce discontinuous jumps in the critical Taylor number. Finally, for both wide and thin gaps, the fastest-growing modes found in the pure corotation case are shown to have their origins in the instability islands at smaller values of rotation rate ratios.

97 MATHEMATICS AND COMPUTING↗