Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “average atom”

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

Thermodynamic origin of the pressure-induced Invar effect: General criterion and experimental study of Fe 68 ⁢Pd 32

Synchrotron X-ray diffraction measurements were performed on 57 Fe 68 Pd 32 at multiple pressures and two temperatures in a diamond-anvil cell. Between 4 and 11GPa, the thermal expansion was zero or slightly negative. This pressure-induced Invar effect was studied further with 57 Fe nuclear forward scattering and nuclear resonant inelastic X-ray scattering to obtain information on the pressure-induced changes in both the magnetization and the phonon density of states. Magnetic entropy and phonon entropy were obtained from these results, the latter with additional measurements from inelastic neutron scattering to account for the contributions from Pd atoms. The dependencies of these entropies on pressure gave the magnetic and phonon contributions to thermal expansion. These canceled in the region of pressure-induced Invar behavior, even though they individually increased by more than a factor of 2 below the Curie pressure. The behavior of phonons gives evidence for spin-phonon interactions in Fe 68 Pd 32 . A general explanation of the pressure-induced Invar effect is presented, showing that Invar behavior is typically expected at a pressure P* below a magnetic transition at pressure PC. The difference in pressure P C − P* scales with the fractional reduction in magnetic exchange interaction divided by the average Gr¨uneisen parameter.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Electronic structure and exchange interactions in altermagnetic MnGeP 2 in the quasiparticle self-consistent G W approach

The quasiparticle self-consistent 𝐺⁡𝑊 method is used to study the electronic band structure, optical dielectric function, and exchange interactions in chalcopyrite, 𝐼⁢ ‾ 4 ⁢2⁢𝑑, structure MnGeP 2 . The material is found to be an antiferromagnetic semiconductor with a lowest direct gap of 2.44 eV at the Γ point and a lower indirect gap of 1.87 eV from Γ to 𝑀. The material is an altermagnet because the two magnetic atoms of opposite spin are related by a twofold rotation operation perpendicular to the main fourfold rotation inversion axis. The spin splittings along a low symmetry line like 𝑃⁢𝑁 is sizable while at k points on the diagonal mirror planes or on the twofold symmetry axes the spin splitting is zero. The exchange interactions are calculated using a linear response approach. The antiferromagnetic exchange interaction between nearest neighbors in the primitive unit cell is dominating and found to be slightly decreasing upon carrier doping but not sufficiently to change the interaction to become ferromagnetic. The bare (noninteracting) and interacting transverse spin susceptibilities, which provide interatomic site exchange interactions after averaging over the muffin-tin spheres, are calculated from the 𝐺⁡𝑊 band structure and wave functions. From these exchange interactions, the spin wave spectra are obtained along the high symmetry lines and the Néel temperature is calculated using the mean-field and Tyablikov estimations. The dielectric function and the optical absorption spectra are calculated including excitonic effects using the Bethe Salpeter equation. The exchange interactions around Mn Ge defect sites is also studied. While we find it can generate ferromagnetic interactions with neighboring spins, we did not find direct evidence of producing an overall ferromagnetic phase. First, if Mn antisites are introduced by exchanging Mn with a nearby Ge, the interactions stay largely antiferromagnetic. Second, when we add additional Mn, in other words in Mn-rich stoichiometry, the Mn Ge antisites produce a strong ferromagnetic interaction primarily with the Mn in the same basal plane but weaker ferromagnetic interaction with adjacent plane Mn. The interactions between regular lattice Mn stay antiferromagnetic as before and thus favor keeping the antiferromagnetic order along the [001] direction. Adding Mn antisites, however, does lead to a metallic band structure.

36 MATERIALS SCIENCE↗
NASA NTRS1987

Oscillator strengths for S I, S II, and S III

A series of calculations for atomic data of various sulfur and oxygen ions is examined. Recent observations of the Io plasma torus obtained with the Voyager UV Spectrometer, the IUE satellite short wavelength spectrograph, and the rocket-borne faint object telescope are discussed. The calculation of oscillator strengths for S II, the P I sequence, S I, and S III in terms of configuration interaction effects is described. The derivation of orbital wave functions is considered. The use of the close coupling method to estimate collision strengths is studied. The accuracy of these calculations depends on: (1) the number of states used in the close coupling expansion; (2) resonance contributions to the thermally averaged collision strength; and (3) the quality of the target state wave functions. Tables of the derived oscillator strengths are presented.

Ho, Y. K.↗
NASA NTRS1994

On the apparent velocity of integrated sunlight. 2: 1983-1992 and comparisons with magnetograms

We report additional results in our program to monitor the wavelength stability of lines in the 2.3 micrometer spectrum of integrated sunlight. We use the McMath Fourier transform spectrometer (FTS) of the National Solar Observatory to monitor 16 delta V = 2 lines of (12)C(16)O, as well as five atomic lines. Wavenumber calibration is achieved using a low-pressure N2O absorption cell and checked against terrestrial atmospheric lines. Imperfect optical integration of the solar disk remains the principal source of error, but this error has been reduced by improved FTS/telescope collimation and observing procedures. The present results include data from an additional 13 quarterly observing runs since 1985. We continue to find that the apparent velocity of integrated sunlight is variable, in the sense of having a greater reshift at solar maximum. This is supported by the temporal dependence of the integrated light velocity, and by the presence of a correlation between velocity and the disk-averaged magnetic flux derived from Kitt Peak magnetograms. The indicated peak-to-peak apparent velocity amplitude over a solar cycle is approximately the same as the velocity amplitude of the Sun's motion about the solar system barycenter. This represents about half the amplitude which we inferred in Paper I (Deming et al. 1987), but the present result has a much greater statistical significance. Our results have implications for those investigations which search for the Doppler signatures of planetary-mass companions to solar-type stars. We contrast our results to the recent finding by McMillan et al. 1993 that solar absorption lines in the violet spectral region are wavelength-stable over the solar cycle.

Deming, Drake↗
NASA NTRS2001

In Situ Measurements of Meteoric Ions

Extraterrestrial material is the source of metal ions in the Earth's atmosphere, Each year approx. 10(exp 8) kg of material is intercepted by the Earth. The origin of this material is predominantly solar orbiting interplanetary debris from comets or asteroids that crosses the Earth's orbit. It contains a very small amount of interstellar material. On occasion the Earth passes through enhanced amounts of debris associated with the orbit of a decaying comet. This leads to enhanced meteor shower displays for up to several days. The number flux of shower material is typically several times the average sporadic background influx of material. Meteoric material is some of the earliest material formed in the solar system. By studying the relative elemental abundances of atmospheric metal ions, information can be gained on the chemical composition of cometary debris and the chemical makeup of the early solar system. Using in situ sampling with rocket-borne ion mass spectrometers; there have been approximately 50 flights that made measurements of the metal ion abundances at attitudes between 80 and 130 km. It is this altitude range where incoming meteoric particles am ablated, the larger ones giving rise to visible meteor. displays. In several rocket measurements isotopic ratios of different atomic ion mass components and metal molecular ion concentrations have been determined and used to identify unambiguously the measured species and to investigate the processes controlling the metal ion distributions The composition of the Earth's ionosphere was first sampled by an ion mass spectrometer flown an a rocket in 1956. In 1958 a rocket-borne ion spectrometer identified, fbr the first time, a layer of metal ions near 95 km. These data were interpreted as evidence of an extraterrestrial rather than a terrestrial source. Istomin predicted: "It seems probable that with some improvement in the method that analysis of the ion composition in the E-region may be used for determining the chemical composition of those meteors which do not reach the ground. Particularly, we hope to get information about the composition difference between particles of different meteor showers and also sporadic and shower meteoroids". These visions categorized the aims of many subsequent rocket-borne ion mass spectrometer experiments in the lower ionosphere, Although the use such measurements to deduce the composition of different classes of meteoroids has not been successful, the past four decades of rocket observations have provided po%erful sets of data for advancing our understanding of meteor ablation, meteoric composition, metal neutral and ion chemistry as well as ionospheric dynamics.

Grebowsky, Joseph M.↗
DOE OSTI2026

QCalEval: Benchmarking Vision-Language Models for Quantum Calibration Plot Understanding

Quantum computing calibration depends on interpreting experimental data, and calibration plots provide the most universal human-readable representation for this task, yet no systematic evaluation exists of how well vision-language models (VLMs) interpret them. We introduce QCalEval, the first VLM benchmark for quantum calibration plots: 243 samples across 87 scenario types from 22 experiment families, spanning superconducting qubits and neutral atoms, evaluated on six question types in both zero-shot and in-context learning settings. The best general-purpose zero-shot model reaches a mean score of 72.3, and many open-weight models degrade under multi-image in-context learning, whereas frontier closed models improve substantially. A supervised fine-tuning ablation at the 9-billion-parameter scale shows that SFT improves zero-shot performance but cannot close the multimodal in-context learning gap. As a reference case study, we release NVIDIA Ising Calibration 1, an open-weight model based on Qwen3.5-35B-A3B that reaches 74.7 zero-shot average score.

Cao, Shuxiang↗
DOE OSTI2020

Experimental investigation of laser ablation of stone polycrystalline targets

We report the results of an experimental investigation of ablation of stone polycrystalline targets of complex multicomponent composition, which imitate the substance of asteroids. The targets were irradiated by nanosecond pulses of a neodymium laser at an energy density Φ{sub L} of up to 5 × 10{sup 4} J cm{sup −2}. The experiments demonstrated the existence of two ablation regimes, with the boundary between them lying at Φ{sub L} ≈ 4000 J cm{sup −2}. The regime change is characterised by a change in the form of the dependence of the surface mass density of removed target material on the laser energy density and by the appearance of a minimum in the dependence of specific energy of destruction on Φ{sub L}. This is supposedly related to the passage from a one-dimensional plasma plume expansion to the three-dimensional one and the corresponding decrease in the efficiency of energy transfer from the laser beam to the target due to a lowering of laser-produced plasma density. Our experiments also showed the existence of a maximum in impulse coupling coefficient C{sub m} as a function of laser energy density (C{sub m} ≈ 6.3 × 10{sup −5} N W{sup −1} for Φ{sub L} = Φ{sub opt} ≈ 25 J cm{sup −2}). Maxima were also recorded in the dependences of the ablation efficiency and average ablation flow velocity on Φ{sub L}. For Φ{sub L} > Φ{sub opt}, the decrease in the function C{sub m}(Φ{sub L}) turns out to be much steeper than for metals and polymer materials. The difference is supposedly due to the lower strength and lower plasticity of the polycrystalline stone targets. (paper)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

The lattice parameter – composition relationship of the body centered cubic uranium-niobium alloys

In this paper, in-situ time-of-flight neutron diffraction experiments were performed on uranium-niobium alloy with 6 wt% Nb to study the lattice parameter-composition relationship for the body centered cubic (bcc) alloys. Based on lattice-parameter measurements of Nb-supersaturated γs phase over the range of 450 – 790 °C, an improved Vegard's-type relationship was established for determination of Nb concentrations at elevated temperatures. Neutron diffraction data were also collected as a function of aging time over the 450 - 600 °C range when γs phase decomposed isothermally into the mixture of orthorhombic α-U and Nb-rich bcc phase, as well as at room-temperature on the alloys aged ex-situ at 500 °C up to five years. From these measurements, the average Nb concentrations in the bcc phase were determined based on the Rietveld refinements of weight fraction and mass conservation relations (lever rule). Over the 15at% – 78at% range of Nb concentrations that correspond to different experimental aging times, the lattice parameters at constant temperatures exhibit a nonlinear S-shaped variation with Nb concentration, and the associated excess volumes of mixing can be described by a subregular solution model of the Redlich–Kister type of polynomial. Over the full range of composition the S-shaped deviation from Vegard's law can be modeled using a combination of an elastic continuum model and a perturbation to the radii of the solute atoms in the solvent, suggesting that electronic interactions between solute and solvent atoms could play an important role in the compositional dependence of lattice parameter for the γ-phase U-Nb alloys. While Vegard's law is a straightforward and reasonably good approximation for the bcc solid solutions in the U-Nb system, the Nb concentrations determined from the weight fraction refinements of diffraction data provide internally consistent, mass-conserving estimates of solute redistribution for the monotectoid reaction.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Single-size and cluster dynamics modeling of intra-granular fission gas bubbles in UO 2

For this work, we perform simulations of intra-granular fission gas bubble evolution in UO 2 using both a relatively simple, computationally inexpensive single-size model and a detailed cluster dynamics model. Simulations encompass 36 experimental cases from 4 different databases, covering various temperature and burnup levels. We systematically compare results from the two models to each other and to post-irradiation experimental data of bubble average size and number density. Overall, the model-to-model comparisons reveal an excellent agreement across the set of simulations. This outcome indicates that, in spite of the underlying assumptions, the single-size model provides a good approximation of the complex physical behavior that is more rigorously described by the cluster dynamics model. Qualitatively, both models reproduce the trends of the experimental data with temperature and burnup correctly. Quantitatively, calculated results are either in good agreement with the data or within errors that appear consistent with the inherent uncertainties. Moreover, for the single-size model, we demonstrate and assess a multiscale approach whereby values for the fission gas atom diffusion coefficient from separate atomistic calculations are used. Systematic comparisons to experimental data point out a credible accuracy of the multiscale model.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Earth-abundant Li-ion cathode materials with nanoengineered microstructures

Manganese-based materials have tremendous potential to become the next-generation lithium-ion cathode as they are Earth abundant, low cost and stable. Here we show how the mobility of manganese cations can be used to obtain a unique nanosized microstructure in large-particle-sized cathode materials with enhanced electrochemical properties. By combining atomic-resolution scanning transmission electron microscopy, four-dimensional scanning electron nanodiffraction and in situ X-ray diffraction, we show that when a partially delithiated, high-manganese-content, disordered rocksalt cathode is slightly heated, it forms a nanomosaic of partially ordered spinel domains of 3–7 nm in size, which impinge on each other at antiphase boundaries. The short coherence length of these domains removes the detrimental two-phase lithiation reaction present near 3 V in a regular spinel and turns it into a solid solution. This nanodomain structure enables good rate performance and delivers 200 mAh g –1 discharge capacity in a (partially) disordered material with an average primary particle size of ~5 µm. The work not only expands the synthesis strategies available for developing high-performance Earth-abundant manganese-based cathodes but also offers structural insights into the ability to nanoengineer spinel-like phases.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Shafranov-shift destabilization of ballooning-type micro-instabilities

Contrary to common wisdom, we argue that the Shafranov shift is destabilizing for the ambient ballooning-type instabilities, which account for most of the transport flux in tokamak core plasmas. Higher Shafranov shift indeed reduces the magnetic drift frequency ω d , i.e. improves the bad curvature mildly, around the outboard midplane. However, this improvement is limited to a very narrow spatial region while ω d is increased over the remaining poloidal space. The eigenfunction averaged $\langle$ω d $\rangle$ is effectively enhanced due to the finite mode width of the eigenfunction and is thus destabilizing to the ballooning-type mode, as demonstrated by gyrokinetic simulations with the CGYRO code using local Miller equilibrium geometry. The predicted nonlinear flux also increases with Shafranov shift, consistent with linear simulations. The reduced transport model TGLF can capture the physics reasonably well.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Molecular-level design of alternative media for energy-saving pilot-scale fibrillation of nanocellulose

The outstanding mechanical properties, light weight, and biodegradability of cellulose nanofibrils (CNFs) make them promising components of renewable and sustainable next-generation reinforced composite biomaterials and bioplastics. Manufacturing CNFs at a pilot scale requires disc-refining fibrillation of dilute cellulose fibers in aqueous pulp suspensions to shear the fibers apart into their nanodimensional forms, which is, however, an energy-intensive process. Here, in this study, we used atomistic molecular dynamics (MD) simulation to examine media that might facilitate the reduction of interactions between cellulose fibers, thereby reducing energy consumption in fibrillation. The most suitable medium found by the simulations was an aqueous solution with 0.007:0.012 wt.% NaOH:urea, and indeed this was found in pilot-scale experiments to reduce the fibrillation energy by ~21% on average relative to water alone. The NaOH:urea-mediated CNFs have similar crystallinity, morphology, and mechanical strength to those formed in water. The NaOH and urea act synergistically on CNFs to aid fibrillation but at different length scales. NaOH deprotonates hydroxyl groups leading to mesoscale electrostatic repulsion between fibrils, whereas urea forms hydrogen bonds with protonated hydroxyl groups thus disrupting interfibril hydrogen bonds. This suggests a general mechanism in which an aqueous medium that contains a strong base and a small organic molecule acting as a hydrogen-bond acceptor and/or donor may be effectively employed in materials processes where dispersion of deprotonable polymers is required. The study demonstrates how atomic-detail computer simulation can be integrated with pilot-scale experiments in the rational design of materials processes for the circular bioeconomy.

36 MATERIALS SCIENCE↗
DOE OSTI2021

A Valentine’s tribute to Lab love: Watch clips from 1943 wedding of two staff

She was a secretary for Director J. Robert Oppenheimer. He was a physicist who helped create the world’s first atomic bombs. And for 65 years they were Mr. and Mrs. Bradner. Hugh Bradner was the first person that Marjorie Hall met upon arrival to join the Lab’s clandestine war-time efforts. He happened to be in Dorothy McKibbin’s Santa Fe office when Marjorie walked in. It was love at first sight, Marjorie recalled years later in an interview. The Bradners were married in September 1943 at McKibbin’s home. It was one of many weddings there during the Manhattan Project era. Los Alamos was fenced and guarded during this time, as staff raced to secretly develop nuclear weapons to help end World War II. Most were young – the average age was late 20s – and many, like the Bradners, met their lifelong partners then. The Bradners left Los Alamos after the war ended in 1945. Hugh is credited with inventing the neoprene wetsuit in 1951. The couple had one daughter. Marjorie, 89, and Hugh, 92, died just a few weeks apart in 2008 in California. Many of their home movies from their time in Los Alamos are preserved at the Lab. They are part of the collections in the National Security Research Center, which is the Lab’s classified library and also houses unclassified artifacts.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2024

Assessment of the Impact of Realistic Sensor Physics and the Integration of Ex-Core Sensors on Reactor Power Synthesis

In the work documented in this report, a weighting function–based core power synthesis method was applied to multiple Monte Carlo N-Particle (MCNP) reactor models, which are informed based on simulated self-powered neutron detector (SPND) responses. The weighting function method used has been coined the point-based iterative (PBI) method. The goal of this application is to assess the impact of considering realistic sensor physics in the generation of the simulated SPND outputs as well as to consider how the synthesis is impacted based on the inclusion of ex-core detectors in the model. The NuScale small modular reactor (SMR) and Westinghouse AP1000 pressurized water reactor (PWR) are the models that served as the testbeds for the assessment of realistic sensor physics; this was achieved by using Geant4 SPND models in comparison with analytical models, such that the effect of electron transport in realistic SPND geometries in the Geant4 model can be understood in terms of synthesis error and convergence time. The comparison was considered for fuel burnup–induced perturbations, for a range of sensor string densities and synthesized power distribution axial fidelities. The Texas A&M Testing, Research, Isotopes, General Atomics Reactor (TAMU TRIGA) reactor MCNP model was used to assess the impact of ex-core sensors; this was done by performing synthesis with and without the ex-core detectors and by quantifying the synthesis error and number of iterations associated with Gaussian-type perturbations in many locations in the core. The TAMU TRIGA model was particularly pertinent for this study because of the interest in future experimental tests with SPNDs in this reactor, as well as the ease of modifying the MCNP model to include ex-core detectors with heterogeneously described response functions. Results from the comparison between the Geant4 and analytical SPND models indicate that similar average and maximum synthesis errors were obtained for burnup-induced perturbations in both the NuScale SMR and the AP1000. This was true for a range of sensor string densities and axial fidelities. However, there were marked differences between both the Geant4 and analytically informed models in terms of the iterations required to converge on the synthesized power distribution. Namely, the Geant4-informed models tended to lead to fewer iterations, except for a few sensor–core configurations that had particularly numerous iterations. Results from the ex-core sensor assessment with the TAMU TRIGA model indicate that the inclusion of ex-core sensors drastically reduces the synthesis error of Gaussian-type perturbations close to the edge of the core, and it slightly reduces synthesis errors for perturbations closer to the center of the core. This was achieved with a minimal increase in computational cost—that is, the number of iterations required for convergence. The errors were identified to be in the same location as the perturbation in the core, indicating that the methodology remains robust for unperturbed regions of the core. A secondary result from this study with the TAMU TRIGA was yielded by analysis of the neutron flux levels in the in-core and ex-core sensor locations of the core; these flux levels indicate that SPNDs could be used as both in-core and ex-core sensors, so long as the emitter material is sensitive to thermal neutrons. The results from these studies provide a quantitative understanding of the importance of considering realistic sensor physics and including ex-core sensors to perform accurate and timely power distribution synthesis of a reactor core.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS1995

The neutral atomic phases of the interstellar medium

We calculate the thermal equilibrium gas temperature of the diffuse interstellar medium. Our method incorporates a new photoelectric heating rate from small grains and polycyclic aromatic hydrocarbons (PAHs) that accounts for a size distribution of particles extending from 100 to 3 A radius. We also include a detailed treatment of the ionization rates and heating due to the soft X-ray background and due to cosmic rays. Phase diagrams (thermal pressure P versus hydrogen density n) are presented for gas that is illuminated by local interstellar far-ultraviolet (FUV) and X-ray radiation fields. A stable two-phase medium is produced with thermal pressure in the range P/k approximately = to 10(exp 3-4) K/cc. We demonstrate that photoelectric heating from PAHs dominates in the warm neutral phase (WNM) and cold neutral phase (CNM). If the C II (158 micrometers cooling per hydrogen nucleus in the solar neighborhood represents an average value for the Galaxy, we predict L(sub CII) approximately = to 7 x 10(exp 7) solar luminosities from the CNM in the Galaxy, comparable to that observed by the Cosmic Background Explorer (COBE). We discuss the dependence of the results on absorbing column density, gas phase abundances, dust abundances and metallicity, FUV field, and the X-ray radiation field. These results will be useful in modeling the multiphase structure of high-velocity clouds in the halo, the interstellar matter (ISM) at other galactocentric radii, and the ISM in external galaxies and galactic nuclei.

Wolfire, M. G.↗
DOE OSTI2022

Thermodynamics of the Metal Carbonates and Bicarbonates of Mn, Co, Ni, Cu, and Zn Relevant to Mineral Energetics

The gas phase heats of formation of ground-state MCO 3 , M(HCO 3 ) 2 , and M(HCO 3 )(OH), where M = Mn, Co, Ni, Cu, and Zn, have been predicted using the correlated molecular orbital theory at the CCSD(T) level extrapolated to the complete basis set limit using the Feller–Peterson–Dixon (FPD) approach. Cohesive energies of the carbonates were predicted based on the calculated gas phase and experimental solid heats of formation. Coulombic dissociation energies (CDEs) between metal cations and anions show a near-linear correlation with Shannon metal cation atomic radii, yet no correlation is found with the hardness of these cations. The total reaction dissociation energies (TRDEs) of transition metals are higher than their CDEs for the di-bicarbonates, in contrast to those for Mg and Ca based on our prior work. In addition to differences in the energies needed to prepare the transition metal dications, electron donation from the ligands to the 3d orbitals of open-shell transition metal dications from lone pairs of adjacent O atoms also plays a role. No electron donation from the ligands to the fully occupied 3d orbitals of Zn and Cd was found. Decomposition energies for generating MO, CO 2 , and/or H 2 O were calculated. Gas phase metal exchange energies only partially correlate with the electrochemical series for M(s) → M 2+( aq). The FPD heats of formation were used to benchmark a range of density functional theory exchange–correlation functionals, including those commonly used in solid-state mineral calculations. None of the functionals provided chemical accuracy agreement (±1 kcal/mol) with the FPD results. So, the best agreement with the FPD results is predicted for the $\tau$-HCTH functional with an average unsigned error of 8.3 kcal/mol.

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

Propagation of priors for more accurate and efficient spectroscopic functional fits and their application to ferroelectric hysteresis

Multi-dimensional spectral-imaging is a mainstay of the scanning probe and electron microscopies, micro-Raman, and various forms of chemical imaging. In many cases, individual spectra can be fit to a specific functional form, with the model parameter maps, providing direct insight into material properties. Since spectra are often acquired across a spatial grid of points, spatially adjacent spectra are likely to be similar to one another; yet, this fact is almost never used when considering parameter estimation for functional fits. On datasets tried here, we show that by utilizing proximal information, whether it be in the spatial or spectral domains, it is possible to improve the reliability and increase the speed of such functional fits by ~2-3x, as compared to random priors. We explore and compare three distinct new methods: (1) spatially averaging neighborhood spectra, and propagating priors based on functional fits to the averaged case, (2) hierarchical clustering-based methods where spectra are grouped hierarchically based on response, with the priors propagated progressively down the hierarchy, and (3) regular clustering without hierarchical methods with priors propagated from fits to cluster means. Our results highlight that utilizing spatial and spectral neighborhood information is often critical for accurate parameter estimation in noisy environments, which we show for ferroelectric hysteresis loops acquired on a prototypical PbTiO3 thin film with piezoresponse spectroscopy. This method is general and applicable to any spatially measured spectra where functional forms are available. Examples include exploring the superconducting gap with tunneling spectroscopy, using the Dynes formula, or current-voltage curve fits in conductive atomic force microscopy mapping. Here we explore the problem for ferroelectric hysteresis, which, given its large parameter space, constitutes a more difficult task than, for example, fitting current-voltage curves with a Schottky emission formula.

42 ENGINEERING↗
DOE OSTI2023

Towards Realistic and High Fidelity Models for Nuclear Reactor Power Synthesis Simulation with Self-Powered Neutron Detectors

As presented in this report, a weighting function–based inferencing method is being applied to synthesize the power distribution in next-generation and university research reactors based on simulated self power neutron detector (SPND) responses. The overall goal is to assess the impacts of sensor uncertainty and true power distribution perturbations on the error in the synthesized power distribution. Regarding sensor uncertainty, the NuScale Small Modular Reactor (SMR) and the Westinghouse AP1000 serve as testbeds for analyzing the impact of varying the sensor uncertainty, as well as varying the number of sensors per sensor string in the reactor core. The reactor models are informed by Monte Carlo N-Particle (MCNP) neutron flux tallies. For the NuScale SMR and Westinghouse AP1000, the SPND response functions (i.e., the response of the SPNDs to individual segments of fuel) were determined homogeneously. Regarding an analysis of power distribution perturbation detection, the Texas A&M Testing, Research, Isotopes, General Atomics Reactor (TAMU TRIGA) reactor was used as a demonstration case with one particular arrangement of SPNDs; the response functions for this reactor model were determined heterogeneously, making this a uniquely high-fidelity demonstration of perturbation detection. Finally, SPND models generated in the Geometry and Tracking 4 (Geant4) code have been generated and tested for comparison with traditionally implemented analytical SPND models, with the intent for Geant4 integration with the full methodological framework. SPND current outputs as a function of distance from some fuel assembly segment in the NuScale SMR are compared with the analytically determined currents. Results from the sensor uncertainty simulations for the NuScale SMR and AP1000 indicate that the average error in the inferred power distribution on the fuel assembly segment level is reasonably low, being slightly less than the random uncertainty applied to all respective SPNDs in both cores. For example, if all SPNDs in the core have a random uncertainty of 5%, then the corresponding fuel assembly segment level error (i.e. difference between the true and inferred local power) is ~2–3%. However, the maximum error in the inferred power distribution on the fuel assembly segment level can be considerably high (>15%) when SPND random uncertainties start to exceed ~3%. In general, the average and maximum errors in the inferred power distribution were slightly higher in the AP1000 as opposed to the NuScale SMR for the sensor string configurations considered herein. Another result determined from analysis of the sensor uncertainty simulations was that increasing the number of SPNDs per string does not clearly reduce inferred power distribution error and can in fact make the error large in some cases; however, this assessment may skewed due to imposed iteration limits. Results from the perturbation detection demonstration using the high-fidelity TAMU TRIGA model indicate that, given the arrangement of 17 SPND strings and 4 SPNDs per string considered herein, there is a clear, provable ability to infer a localized Gaussian-type peak perturbation in the 3D power distribution. Such a perturbation was detected with an average fuel assembly segment level error of 0.19%, and the general visualization of the detected perturbation clearly indicates that the magnitude and shape were appropriately resolved. Finally, the electrical current output generated by the Geant4 modeled SPND indicates significant magnitude differences than the analytically modeled SPND, demonstrating the need for accurate SPND models which account for finite sensor geometry effects to inform the power synthesis work described herein.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗