Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “qualitative methods”

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 541 records · Page 30
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Simulation Center for Runaway Electron Avoidance and Mitigation (SCREAM SciDAC) (Technical Final Report)

Runaway electrons can severely damage the plasma facing components on ITER during a major disruption and pose a major risk for tokamak fusion. It has been recognized that an adequate disruption mitigation system (DMS) is essential for the safe operation of ITER. The United States is responsible for the design and implementation of the disruption mitigation system on ITER, and in July 2016 the Simulation Center for Runaway Electron Avoidance and Mitigation (SCREAM) was launched by DOE, in a joint Fusion Energy Sciences (FES) and Advanced Scientific Computing Research (ASCR) collaboration. SCREAM was a comprehensive theory and simulation SciDAC center that provided physics guidance in the avoidance and mitigation of runaway electrons, and in tandem with domestic and international experiments, helped establish the qualitative and quantitative bases for safe operational scenarios and viable mitigation techniques. The SCREAM center assembled a national team of experts in runaway electron physics, tokamak disruptions, magnetohydrodynamic (MHD) simulation, and advanced algorithms and computing. The team combined advanced simulation and analysis capability facilitated by direct participation of ASCR SciDAC institutes with theoretical models and code development by FES scientists to focus on the runaway risk for ITER and tokamaks in general. The research scope was focussed on integrated simulations of kinetic runaway electrons, including MHD and fluid models of impurity transport, within a research plan guided by theory. The specific research tasks were (1) establish the fundamental physics of runaway generation, saturation, and dynamical evolution in a tokamak; (2) examine the critical path toward runaway avoidance; and (3) investigate the viability and effectiveness of the leading candidate schemes for runaway mitigation. In all three areas, members of the team carried out scoping studies that established the readiness for rapid and critical advances, especially in the deployment and further development of large-to extreme-scale simulation tools. Our multi-pronged computational approach included (1) relativistic Fokker-Planck solvers with discretization in phase space, (2) self-consistent particle-in-cell techniques, (3) particle-based Monte-Carlo, and (4) MHD-particle hybrid simulations. Cross-check between these different methods provided an additional means for verification and further bolstered the fidelity of our physics prediction. Validation against experimental results brings confidence to the predictive capability for ITER and frequently leads to new ideas for understanding and mitigating the thermal quench driven runaway electron phenomenon.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Polysulfide entrapment and retardation in gel electrolyte Li–S batteries: experiments and modeling

Gelling the electrolyte in lithium–sulfur (Li–S) batteries provides advantages in terms of safety and mitigating polysulfide shuttling leading to reduced side reactions and dendrite growth on the anode side. In the current study, organic and ceramic crosslinkers have been used to improve the entrapment of polysulfides and prevent them from shuttling to the anode side. However, the performance of gel electrolytes is still far from being on par with liquid electrolytes in terms of raw capacity. The discharge curves exhibit an extended first plateau for gel electrolyte batteries compared to liquid electrolytes, showing that polysulfides are better trapped in the gel system. Yet, in the conversion of soluble polysulfides into insoluble polysulfides, the second plateau was delayed. Qualitative and quantitative analysis of the polysulfide diffusion in the liquid electrolyte and gel electrolyte systems were also conducted to predict the relative diffusivity ratio between the liquid electrolyte and gel electrolyte for simulation. To probe what happens at the end of the discharge, where the mass transport in the gel electrolyte appears to be slow, a modified discharge procedure was used where the cell is first paused after discharge and then discharged a second time before recharging normally. Coin cell tests for two gel electrolytes in addition to liquid electrolyte show up to 7% capacity recovery after the pause, which indicates that the diffusion of species during the pause is responsible for the failure recovery. Additionally, the recovery is higher for gel electrolytes compared to liquid, for higher C-rates, and for longer pause times. To better understand this behavior, a Li–S numerical model with mass transport limited reactions was used to examine the polysulfide diffusion coefficients expected in different electrolyte systems. The model is able to reproduce the trends seen in experiments and yields a higher recovery for smaller Li 2 S 4 diffusion coefficients, suggesting that insufficient Li 2 S 4 mass transport is responsible for failure at the end of discharge. By identifying the effect that gel electrolytes have on the state of the battery during discharge using both experimental and simulation methods, this work is a step towards leveraging the advantages of gel electrolytes and mitigating their weaknesses.

25 ENERGY STORAGE↗
DOE OSTI2021

Collective Nanoparticle Dynamics Associated with Bridging Network Formation in Model Polymer Nanocomposites

The addition of nanoparticles (NPs) to polymers is a powerful method to improve the mechanical and other properties of macromolecular materials. Such hybrid polymer–particle systems are also rich in fundamental soft matter physics. Among several factors contributing to mechanical reinforcement, a polymer-mediated NP network is considered to be the most important in polymer nanocomposites (PNCs). Here, we present an integrated experimental–theoretical study of the collective NP dynamics in model PNCs using X-ray photon correlation spectroscopy and microscopic statistical mechanics theory. Silica NPs dispersed in unentangled or entangled poly(2-vinylpyridine) matrices over a range of NP loadings are used. Static collective structure factors of the NP subsystems at temperatures above the bulk glass transition temperature reveal the formation of a network-like microstructure via polymer-mediated bridges at high NP loadings above the percolation threshold. The NP collective relaxation times are up to 3 orders of magnitude longer than the self-diffusion limit of isolated NPs and display a rich dependence with observation wavevector and NP loading. A mode-coupling theory dynamical analysis that incorporates the static polymer-mediated bridging structure and collective motions of NPs is performed. It captures well both the observed scattering wavevector and NP loading dependences of the collective NP dynamics in the unentangled polymer matrix, with modest quantitative deviations emerging for the entangled PNC samples. Here, we identify an unusual and weak temperature dependence of collective NP dynamics, in qualitative contrast with the mechanical response. Hence, the present study has revealed key aspects of the collective motions of NPs connected by polymer bridges in contact with a viscous adsorbing polymer medium and identifies some outstanding remaining challenges for the theoretical understanding of these complex soft materials.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Robustness of neutron star merger simulations to changes in neutrino transport and neutrino-matter interactions

Binary neutron star mergers play an important role in nuclear astrophysics: Their gravitational wave and electromagnetic signals carry information about the equation of state of cold matter above nuclear saturation density, and they may be one of the main sources of r-process elements in the Universe. Neutrino-matter interactions during and after merger impact the properties of these electromagnetic signals and the relative abundances of the produced r-process elements. Existing merger simulations are, however, limited in their ability to realistically model neutrino transport and neutrino-matter interactions. Here, we perform a comparison of the impact of the use of state-of-the art two-moment or Monte Carlo transport schemes on the outcome of merger simulations, for a single binary neutron star system with a short-lived neutron star remnant [(5–10) ms]. We also investigate the use of different reaction rates in the simulations. While the best transport schemes generally agree well on the qualitative impact of neutrinos on the system, differences in the behavior of the high-density regions can significantly impact the collapse time and the properties of the hot tidal arms in this metastable merger remnant. The chosen interaction rates and transport algorithm as well as recent improvements by Radice et al. to the two-moment algorithms can all contribute to changes at the (10–30)% level in the global properties of the merger remnant and outflows. Furthermore, the limitations of previous moment schemes fixed by Radice et al. also appear sufficient to explain the large difference that we observed in the production of heavy-lepton neutrinos in a previous comparison of Monte Carlo and moment schemes in the context of a low-mass binary neutron star system.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Comparing GGA, GGA+ U , and meta-GGA functionals for redox-dependent binding at open metal sites in metal–organic frameworks

Metal–organic frameworks (MOFs) with open metal sites have been widely investigated for the selective adsorption of small molecules via redox mechanisms where charge transfer can take place between the binding site and the adsorbate of interest. Quantum-chemical screening methods based on density functional theory have emerged as a promising route to accelerate the discovery of MOFs with enhanced binding affinities toward various adsorbates. However, the success of this approach is linked to the accuracy of the underlying density functional approximations (DFAs). In this work, we compare commonly used generalized gradient approximation (GGA), GGA+U, and meta-GGA exchange-correlation functionals in modeling redox-dependent binding at open metal sites in MOFs using O 2 and N 2 as representative small molecules. We find that the self-interaction error inherent to the widely used Perdew, Burke, and Ernzerhof (PBE) GGA predicts metal sites that are artificially redox-active, as evidenced by their strong binding affinities, short metal–adsorbate bond distances, and large degree of charge transfer. The incorporation of metal-specific, empirical Hubbard U corrections based on the transition metal oxide literature systematically reduces the redox activity of the open metal sites, often improving agreement with experiment. Additionally, the binding behavior shifts from strong chemisorption to weaker physisorption as a function of U. The M06-L meta-GGA typically predicts binding energies between those of PBE-D3(BJ) and PBE-D3(BJ)+U when using empirically derived U values from the transition metal oxide literature. Despite the strong sensitivity of the binding affinities toward a given DFA, the GGA, GGA+U, and meta-GGA approaches often yield the same qualitative trends and structure–property relationships.

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

Biofouling and Corrosion Study for a Novel Linear Guided Wave Energy Converter (Post Access Report)

The overall objective of this project was to examine the reliability and performance of antibiofouling coatings used for a wave energy converter (WEC) developed by E-Wave Technologies. The particular coatings were selected for their low toxicity and potential compatibility with aquaculture. The aim of this work was to 1) test coating solutions to prevent biofouling growth and saltwater corrosion on the static (paddle and attachment frame surface) components of the WEC that are submerged, 2) determine adhesion of the coatings to system components, and 3) assess the ease and effectiveness of biofouling cleaning to insure long term performance of the system. An analysis of commercial coatings was performed using methods to examine the prevention of biofouling and coating adhesion properties on two key materials of the WEC, which were 316L low carbon marine grade stainless steel (SS) and Ultra High Molecular Weight Polyethylene (PE). Three marine antifouling paints were selected based on their unique properties to test how different paint styles perform on different materials. The selected paints were ePaint Ecominder self-polishing paint with Zinc Omadine for slime control, Pettit ECO HRT Copper-Free ablative antifouling with Econea biocide, and Intersleek 1100SR foul release. Pacific Northwest National Laboratory (PNNL) prepared PE and SS substrates coated with the three paints and compared the performance against uncoated substrates when submerged in raw seawater for 3-, 6-, and 9-month (m) time periods. Results in adhesion testing indicated that Pettit and ePaint materials clearly bonded strongly to SS, but did not bond comparably well to PE. It was noted during adhesion testing that the Intersleek surfaces were especially difficult to test as the paint highly resists bonding to the epoxy adhesives used with the adherence testing platform. The wear rate of the coatings was not measured under this study; however, based on adhesion testing, coatings in the sliding regions of the Ewave device are expected to wear rapidly. Sandia National Laboratories (SNL) evaluated the adhesion of three different paints to PE and SS substrates which were exposed to a marine environment for time intervals of 0, 3, 6, and 9 months. From qualitative visual analysis of the 3 in2 coupons when pulled from the tank, the 3 in2 coupons generally only appeared to have biofouling consisting of filamentous algae or diatoms, which all have relatively low mass and can be easily wiped from the surface of coupons. Qualitative visual analysis indicated that ECO HRT and Unpainted were consistently worse than Intersleek and Ecominder at all time points. Results provide insight to aid with down-selection of commercial coatings under static conditions to support reliability of the WEC and potential maintenance schedules. This investigation was conducted using small coupon samples suspended in seawater and the development of testing rigs for dynamic component level testing is needed for future work. In addition, the potential toxicity of these commercial coatings on aquaculture has not been determined by this study. One recommendation is to conduct toxicity investigations at the Environmental Toxicity Laboratory at Oak Ridge National Laboratory.

16 TIDAL AND WAVE POWER↗
DOE OSTI2023

A Gridded Solar Irradiance Ensemble Prediction System Based on WRF-Solar EPS and the Analog Ensemble

The WRF-Solar Ensemble Prediction System (WRF-Solar EPS) and a calibration method, the analog ensemble (AnEn), are used to generate calibrated gridded ensemble forecasts of solar irradiance over the contiguous United States (CONUS). Global horizontal irradiance (GHI) and direct normal irradiance (DNI) retrievals, based on geostationary satellites from the National Solar Radiation Database (NSRDB) are used for both calibrating and verifying the day-ahead GHI and DNI predictions (GDIP). A 10-member ensemble of WRF-Solar EPS is run in a re-forecast mode to generate day-ahead GDIP for three years. The AnEn is used to calibrate GDIP at each grid point independently using the NSRDB as the “ground truth”. Performance evaluations of deterministic and probabilistic attributes are carried out over the whole CONUS. The results demonstrate that using the AnEn calibrated ensemble forecast from WRF-Solar EPS contributes to improving the overall quality of the GHI predictions with respect to an AnEn calibrated system based only on the deterministic run of WRF-Solar. In fact, the calibrated WRF-Solar EPS’s mean exhibits a lower bias and RMSE than the calibrated deterministic WRF-Solar. Moreover, using the ensemble mean and spread as predictors for the AnEn allows a more effective calibration than using variables only from the deterministic runs. Finally, it has been shown that the recently introduced algorithm of correction for rare events is of paramount importance to obtain the lowest values of GHI from the calibrated ensemble (WRF-Solar EPS AnEn), qualitatively consistent with those observed from the NSRDB.

14 SOLAR ENERGY↗
DOE OSTI2025

Neutron generation dynamics inside a MA-class dense plasma focus Z-pinch

Dense plasma focii (DPFs) are appealing as energy efficient sources of short pulses of ions, neutrons, and x rays. The output of these sources is expected to scale with input current (I4), but has been shown to drop at the MA level [S. K. H. Auluck, “On the failure of neutron yield scaling in the dense plasma focus,” Phys. Plasmas 30, 080701 (2023)]. New results on the MegaJOuLe Neutron Imaging Radiography DPF showed neutron yield production in agreement with the input current scaling beyond the previously observed drop. This work provides insight into the pinch formation on a DPF and reports on the two different mechanisms leading to neutron generation inside a DPF using a combination of kinetic simulations and experimental data. A combination of particle-in-cell (PIC) and 1D shock theory results are used to describe the pinch formation and disassembly and the corresponding thermonuclear and beam-target mechanisms. The temporal evolution of the pinch column predicted by the PIC simulations shows qualitative agreement with the experimental data from plasma photon emission as well as temporal neutron pulse shapes. In MJ-class DPFs, both thermonuclear and beam-target mechanisms can occur over the course of the implosion and contribute to the total neutron production. Hence the neutron source size of a DPF will change throughout the implosion. Experimental neutron radiographs show the increase in source size as the pinch breaks apart, in agreement with simulation's prediction.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Mineralogical, magnetic and geochemical data constrain the pathways and extent of weathering of mineralized sedimentary rocks

The oxidative weathering of sulfidic rock can profoundly impact watersheds through the resulting export of acidity and metals. Weathering leaves a record of mineral transformation, particularly involving minor redox-sensitive phases, that can inform the development of conceptual and quantitative models. In sulfidic sedimentary rocks, however, variations in depositional history, diagenesis and mineralization can change or overprint the distributions of these trace minerals, complicating the interpretation of weathering signatures. Here we show that a combination of bulk mineralogical and geochemical techniques, micrometer-resolution X-ray fluorescence microprobe analysis and rock magnetic measurements, applied to drill core samples and single weathered fractures, can provide data that enable the development of a geochemically consistent weathering model. This work focused on one watershed in the Upper Colorado River Basin sitting within the Mesaverde Formation, a sedimentary sandstone bedrock with disseminated sulfide minerals, including pyrite and sphalerite, that were introduced during diagenesis and subsequent magmatic-hydrothermal mineralization. Combined analytical methods revealed the pathways of iron (Fe), carbonate and silicate mineral weathering and showed how pH controls element retention or release from the actively weathering fractured sandstone. Drill core logging, whole rock X-ray diffraction, and geochemical measurements document the progression from unweathered rock at depth to weathered rock at the surface. X-ray microprobe analyses of a 1-cm size weathering profile along a fracture surface are consistent with the mobilization of Fe(II) and Fe(III) into acidic pore water from the dissolution of primary pyrite, Fe-sphalerite, chlorite, and minor siderite and pyrrhotite. These reactions are followed by the precipitation of secondary minerals such as of goethite and jarosite, a Fe-(oxyhydr)oxide and hydrous Fe(III) sulfate, respectively. Microscale analyses also helped explain the weathering reactions responsible for the mineralogical transformations observed in the top and most weathered section of the drill core. For example, dissolution of feldspar and chlorite neutralizes the acidity generated by Fe and sulfide mineral oxidation, oversaturating the solution in both Fe-oxides. The combination of X-ray spectromicroscopy and magnetic measurements show that the Fe(III) product is goethite, mainly present either as a coating on fracture surfaces in the actively weathering region of the core or more homogeneously contained within the unconsolidated regolith at the top of the core. Low-temperature magnetic data reveal the presence of ferromagnetic Fe-sulfide pyrrhotite that, although it occurs at trace concentrations, could provide a qualitative proxy for unweathered sulfide minerals because the loss of pyrrhotite is associated with the onset of oxidative weathering. Pyrrhotite loss and goethite formation are detectable through room-temperature magnetic coercivity changes, suggesting that rock magnetic measurements can determine weathering intensity in rock samples at many scales. In conclusion, this work contributes evidence that the weathering of sulfidic sedimentary rocks follows a geochemical pattern in which the abundance of sulfide minerals controls the generation of acidity and dissolved elements, and the pH-dependent mobility of these elements controls their export to the ground- and surface-water.

58 GEOSCIENCES↗
DOE OSTI2025

Laboratory time series moisture manipulative experiment from sediment across San Antonio, Texas: time series aerobic respiration and geochemistry

This dataset supports a broader study examining the effects of wetting and drying on hyporheic zone respiration. The dataset provides data generated from a laboratory moisture manipulation experiment. The contents include time series aerobic respiration and moisture; dissolved oxygen; sediment geochemistry data; and field metadata (including qualitative information on instream and river corridor characteristics). Samples were collected as part of the WHONDRS Allison Veach collaboration (AV1). The data package associated with the AV1 study is available at https://data.ess-dive.lbl.gov/datasets/doi:10.15485/2529428. AV1 sampling occurred across 7 perennial and 7 intermittent streams in San Antonio, Texas. Each stream/site was visited both in summer during base flow (July-September 2023) and winter during peak flow (January-February 2024). This study uses subsamples from a subset of AV1 samples. The original field samples were labeled as AV1_###. Subsequent subsamples for this study were labeled as EV_###. The labels from the field samples and the EV subsamples can be mapped directly based on the digits following the prefix and underscore (i.e., EV_001 is a subsample from AV1_001). See the critical details section below for more details on sample naming. For details on how to navigate data packages generated by this project, see https://data.ess-dive.lbl.gov/portals/PNNLRiverCorridorSFA/About. This dataset is comprised of one main data folder containing (1) file-level metadata; (2) data dictionary; (3) field metadata; (4) readme; (5) field protocol; and a (6) a subfolder with sediment sample data from the incubation experiment. The sample data subfolder contains (1) effect size; (2) iron (II); (3) gravimetric moisture; (4) respiration rates; (5) raw dissolved oxygen values and plots; (6) specific conductance; (7) pH; (8) temperature; (9) a summary containing mean, median, and standard deviation values of each data type for each treatment (wet and dry); and (10) methods codes. All files are .csv or.pdf.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Towards Reconstruction of Complex Flow Fields Using Unit Flows

Many complex turbulent flows in nature and engineering can be qualitatively regarded as being constituted of multiple simpler unit flows. The objective of this work is to characterize the coherent structures in such complex flows as a combination of constituent unitary flow structures for the purpose of reduced-order representation. While turbulence is clearly a non-linear phenomenon, we aim to establish the degree to which the optimally weighted superposition of unitary flow structures can represent the complex flow structures. The rationale for investigating such superposition stems from the fact that the large-scale coherent structures are generated by underlying flow instabilities that may be reasonably described using linear analysis. Clearly, the degree of validity of superposition will depend on the flow under consideration. In this work, we take the first step toward establishing a procedure for investigating superposition. Experimental data of single and triple tandem jets in crossflow are used to demonstrate the procedure. A composite triple tandem jet flow field is generated from optimal superposition of single jet data and compared against ‘true’ triple jet data. Direct comparisons between the true and composite fields are made for spatial, temporal, and kinetic energy content. The large-scale features (obtained from proper orthogonal decomposition or POD) of true and composite tandem jet wakes exhibit nearly 70% agreement in terms of modal eigenvector correlation. Corresponding eigenvalues reveal that the kinetic energy of the flow is also emulated with only a slight overprediction. Temporal frequency features are also examined in an effort to completely characterize POD modes. The proposed method serves as a foundation for more rigorous and robust dimensional reduction in complex flows based on unit flow modes.

Kristo, Paul J.↗
DOE OSTI2020

Evaluation of a visual interpretation method for tau-PET with 18 F-flortaucipir

Introduction: Positron emission tomography targeting tau (tau-PET) is a promising diagnostic tool for the identification of Alzheimer's disease (AD). Currently available data rely on quantitative measures, and a visual interpretation method, critical for clinical translation, is needed. Methods: In this work, we developed a visual interpretation method for 18 F-flortaucipir tau-PET and tested it on 274 individuals (cognitively normal controls, patients with mild cognitive impairment [MCI], AD dementia, and non-AD diagnoses). Two readers interpreted 18 F-flortaucipir PET using two complementary indices: a global visual score and a visual distribution pattern. Results: Global visual scores were reliable, correlated with global cortical 18 F-flortaucipir standardized uptake value ratio (SUVR) and were associated with clinical diagnosis and amyloid status. The AD-like 18 F-flortaucipir pattern had good sensitivity and specificity to identify amyloid-positive patients with AD dementia or MCI. Discussion: This 18 F-flortaucipir visual rating scheme is associated with SUVR quantification, clinical diagnosis, and amyloid status, and constitutes a promising approach to tau measurement in clinical settings.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

A Meta-Level Framework for Evaluating Resilience in Net-Zero Carbon Power Systems with Extreme Weather Events in the United States

Important changes are underway in the U.S. power industry in the way that electricity is sourced, transported, and utilized. Disruption from extreme weather events and cybersecurity events is bringing new scrutiny to power-system resilience. Recognizing the complex social and technical aspects that are involved, this article provides a meta-level framework for coherently evaluating and making decisions about power-system resilience. It does so by examining net-zero carbon strategies with quantitative, qualitative, and integrative dimensions across discrete location-specific systems and timescales. The generalizable framework is designed with a flexibility and logic that allows for refinement to accompany stakeholder review processes and highly localized decision-making. To highlight the framework’s applicability across multiple timescales, processes, and types of knowledge, power system outages are reviewed for extreme weather events, including 2021 and 2011 winter storms that impacted Texas, the 2017 Hurricane Maria that affected Puerto Rico, and a heatwave/wildfire event in California in August 2020. By design, the meta-level framework enables utility decision-makers, regulators, insurers, and communities to analyze and track levels of resilience safeguards for a given system. Future directions to advance an integrated science of resilience in net-zero power systems and the use of this framework are also discussed.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2024

Solid State Quantum Refrigeration Superconducting, Absorption and Measurement Based (Final Technical Report)

During this DOE grant, DE-SC0017890, in place for the past six years, all proposed research was carried out and published in peer-reviewed papers, as well as other projects that emerged during the research. In that effort the research team accomplished all proposed research, as well as many closely related research projects discovered and conceived of during the grant. These works included “Efficient Quantum Measurement Engines”, a work published in Physical Review Letters, giving a theory of quantum measurement-based engines, which uses quantum measurement as a resource. These engines are designed to efficiently convert energy from the stochastic quantum measurement process into useful work. Further publications include “Experimental Realization of a Quantum Dot Energy Harvester”, a joint theory and experimental work in collaboration with the group of Charles Smith in Cambridge, UK, as well as long time theoretical collaborators, Rafael Sánchez and Björn Sothmann. This work, featured as an Editor’s Suggestion in Physical Review Letters, realized an earlier theoretical proposal of ours, whereby two resonant tunneling quantum dots are connected to a central electronic cavity that is heated by a hot energy source. We also published “Superconducting Quantum Refrigerator: Breaking and Rejoining Cooper Pairs with Magnetic Field Cycles” a work done in collaboration with experimentalist Francesco Giazotto from ENS Pisa, Italy, which also resulted in a patent. This paper, published in Phys. Rev. Applied, advanced the concept of a cyclic fridge based on the normal/superconducting phase transition together with layered materials separated by tunnel junctions. We also completed the proposed research on a heat transistor, publishing “Thermal transistor and thermometer based on Coulomb-coupled conductors”, carried out as a collaboration between my group and theorists Splettstoesser (Lund U., Sweden), Sothmann (U. Duisburg-Essen, Germany), and Sánchez (U. Autónoma de Madrid, Spain). We carried out an analysis of a quantum coupled to a quantum point contact as a sensitive thermometer and heat transistor. We found the optimal statistical estimator for the temperature and compared it with experiments on the same type of devices. We also investigated autonomous quantum absorption refrigerators using quantum dots to cool by using a very hot thermal reservoir to drive heat between two other reservoirs. In the article “Quantifying the quantum heat contribution from a driven superconducting circuit”, we demonstrated that for a driven superconducting circuit, we showed heat flow provided by a hot source to the qubit can be switched on and off by varying external parameters, the frequency and the intensity of the driving. In the work “Stochastic thermodynamic cycles of a mesoscopic thermoelectric engine”, we reconsidered the autonomous thermoelectric heat engine in terms of underlying cycles. Rather than periodic behavior, the cycles were stochastic in nature. Nevertheless, by undertaking a graph theoretical analysis of the elementary transport processed, great quantitative and qualitative insight could be found. We also considered the quantum measurement process and showed that a quantum version of Maxwell’s demon could be related to the work extraction of a quantum system, closely related to arrow-of-time measures for quantum measurement, as described in our article “Thermodynamics of quantum measurement and Maxwell's demon's arrow of time”. This work was selected in Phys. Rev. A as an Editor’s Suggestion. A recent preprint titled “Cyclic Superconducting Quantum Refrigerators Using Guided Fluxon Propagation” accomplished an important piece of this grant: to propose a new kind of quantum refrigerator using the dynamics of fluxons in a type II superconductor. This invention envisioned a race-track type geometry where fluxons are confined. By applying a gradient of magnetic field together with electric current in a Corbino geometry, the circulating fluxons can actively cool a cold reservoir, realizing a new type of cyclic superconducting refrigerator. We also investigated the possibility of thermal control from different points of view. The application of quantum measurement to the system gives a new kind of control on the system of interest – we have pioneered this approach and shown that measurement can boost the thermal power of quantum engines as described in “Continuous measurement boosted adiabatic quantum thermal machines”. The ability to have heat flows on demand is an outstanding challenge, and we have provided new solutions to this problem in Thermal control across a chain of electronic nanocavities” for a chain of electron cavities using gating voltage control. The control methods using qubit/qubit coupling to create absorption fridges at their most fundamental level have also been developed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Machine Learned Hückel Theory: Interfacing Physics and Deep Neural Networks

The Hückel Hamiltonian is an incredibly simple tight-binding model known for its ability to capture qualitative physics phenomena arising from electron interactions in molecules and materials. Part of its simplicity arises from using only two types of empirically fit physics-motivated parameters: the first describes the orbital energies on each atom and the second describes electronic interactions and bonding between atoms. By replacing these empirical parameters with machine-learned dynamic values, we vastly increase the accuracy of the extended Hückel model. The dynamic values are generated with a deep neural network, which is trained to reproduce orbital energies and densities derived from density functional theory. The resulting model retains interpretability, while the deep neural network parameterization is smooth and accurate and reproduces insightful features of the original empirical parameterization. Altogether, this work shows the promise of utilizing machine learning to formulate simple, accurate, and dynamically parameterized physics models.

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

U Metal Oxidation with Steam (U)

The Savannah Rive National Laboratory 9SRNL) is evaluating a method to produce uranium oxide (UO 2 ) from uranium metal using steam. The typical methods for making UO 2 yield powders with particle-size ranges of 1 to 50 micron. The goal of the steam oxidation is to produce a depleted uranium (DU) product with a larger particle size that can be used as a surrogate for testing physical and mechanical process techniques to refine uranium and/or plutonium bearing residues. Several studies in the literature indicated that reactions between uranium metal and steam will yield UO 2 as fine powders in the 160 to 500 °C range, but granular materials form above 500°C. Above 880 °C, a hard compact scale forms on the U metal. Furthermore, sintering of UO 2 occurs at about 1300 °C, but it has been stated that the presence of steam enhances the sintering characteristics of UO 2 . Initial studies into the oxidation of DU metal in an atmosphere of argon and steam have been completed. Oxidation of depleted uranium metal in Ar/steam at 600, 710, 700, and 830 °C yielded at UO 2 product with little or no impurity at each temperature. Although not enough sample was available to perform sieve analysis for particle size distribution, SEM images reveal that a significant volume fraction of the products at 710, 770, and 830 °C exceeds 100 microns. The material produced at 600 °C was too fine for programmatic needs. Questions associated with the UO 2 include morphology of product, mechanical strength of the particles, and performance in milling operations. The size of the samples produced in these tests did not allow for quantitative assessment of particle size of mechanical strength. A qualitative assessment of the particles formed at 710, 770, and 830 °C was that the samples produced at 710 and 770 °C were friable, and the sample produced 830 °C had more mechanical strength. Sintering of the samples in Ar/steam at 970 °C appeared to increase the strength of the particles without changing their general particle size characteristics. Due to the small sample sizes, the increase in strength could not be quantified. If the sintering at 970 °C is inadequate, higher sintering temperatures in Ar/steam can be evaluated. Future work will increase scale of the steam oxidation conditions which yield the most-favorable particle-size distributions with sufficient mechanical strength for programming usage. Further consideration will be given to oxidation temperature, steam addition to the reaction vessel, U metal pretreatment, and sintering conditions after oxidation. The UO 2 production could involve a two-step process where particles are formed in steam at lower temperatures (700-850 °C) and then sintered in Ar/steam at elevated temperatures (950-1200 °C). Larger samples (50-100 g each) will provide sufficient UO 2 product for sieve analysis and assessment of mechanical strength relatives to subsequent milling operations.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Many-Body Effects in FeN4 Center Embedded in Graphene

We introduce a computational approach to study porphyrin-like transition metal complexes, bridging density functional theory and exact many-body techniques, such as the density matrix renormalization group (DMRG). We first derive a multi-orbital Anderson impurity Hamiltonian starting from first principles considerations that qualitatively reproduce generalized gradient approximation (GGA)+U results when ignoring inter-orbital Coulomb repulsion U ′ and Hund exchange J. An exact canonical transformation is used to reduce the dimensionality of the problem and make it amenable to DMRG calculations, including all many-body terms (both intra- and inter-orbital), which are treated in a numerically exact way. We apply this technique to FeN 4 centers in graphene and show that the inclusion of these terms has dramatic effects: as the iron orbitals become single occupied due to the Coulomb repulsion, the inter-orbital interaction further reduces the occupation, yielding a non-monotonic behavior of the magnetic moment as a function of the interactions, with maximum polarization only in a small window at intermediate values of the parameters. Furthermore, U ′ changes the relative position of the peaks in the density of states, particularly on the iron d z 2 orbital, which is expected to affect the binding of ligands greatly.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Essential barrier height and a probabilistic approach in characterizing potential landscape

In this work we propose a probabilistic approach to investigate the shape of landscapes of multi-dimensional potential functions. Under a suitable coupling scheme, two copies of the overdamped Langevin dynamics associated with the potential function are coupled, and the coupling times are collected. Assuming a set of intuitive yet technically challenging conditions on the coupling scheme, it is shown that the tail distributions of the coupling times exhibit qualitatively different dependencies on the noise magnitude for single-well versus multi-well potential functions. More specifically, for convex single-well potentials, the negative tail exponent of the coupling time distribution is uniformly bounded away from zero by the convexity parameter and is independent of the noise magnitude. In contrast, for multi-well potentials, the negative tail exponent decreases exponentially as the noise vanishes, with the decay rate governed by the essential barrier height, a quantity introduced in this paper to characterize the non-convex nature of the potential function. Numerical investigations are conducted for a variety of examples, including the Rosenbrock function, interacting particle systems, and loss functions arising in artificial neural networks. These examples not only illustrate the theoretical results in various contexts but also provide crucial numerical validation of the conjectured assumptions, which are essential to the theoretical analysis yet lie beyond the reach of standard technical tools.

97 MATHEMATICS AND COMPUTING↗