Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “parameter sampling”

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 937 records · Page 52

Energy-resolved neutron imaging and diffraction including grain orientation mapping using event camera technology

Time-of-flight neutron diffraction and energy-resolved imaging each provide unique perspectives into material properties. Neutron diffraction is useful for assessing microstructural parameters such as phase composition, texture, and dislocation densities, though it typically provides averaged data over the sampled volume. Energy-resolved imaging, on the other hand, offers both spatial and spectral information by detecting Bragg edges and neutron absorption resonances, which enables detailed mapping of microstructure and isotopic composition. When combined, these techniques have the potential to enrich our understanding of material behavior across different scales, enhancing our understanding of complex materials. Traditionally, these modalities are conducted on separate instruments, which is time-consuming and poses challenges for data integration. Here, we report the integration of the LumaCam, an event-mode energy-resolved neutron imaging camera with the HIPPO time-of-flight diffractometer at LANSCE. This integration enables simultaneous diffraction and imaging across the full spectrum, with analysis optimized for diffraction and Bragg-edge imaging in the thermal range (0.45–10 Å) and resonance imaging in the epithermal range (0.5–3000 eV), facilitating comprehensive multi-modal analysis. We demonstrate its capabilities through case studies, including spatial mapping of grain orientations in a steel sample and accurate thickness estimations for irregular samples including a depleted uranium cylinder and a natural silver-containing mineral specimen. The combined setup enhances real-time sample alignment and provides comprehensive data for crystal structure, texture, and isotopic composition analysis. This approach opens new possibilities for advanced applications in nuclear engineering, archaeology, and materials science.

36 MATERIALS SCIENCE↗

Imaging dynamic exciton interactions and coupling in transition metal dichalcogenides

Transition metal dichalcogenides (TMDs) are regarded as a possible material platform for quantum information science and related device applications. In TMD monolayers, the dephasing time and inhomogeneity are crucial parameters for any quantum information application. In TMD heterostructures, coupling strength and interlayer exciton lifetimes are also parameters of interest. However, many demonstrations in TMDs can only be realized at specific spots on the sample, presenting a challenge to the scalability of these applications. Here, using multi-dimensional coherent imaging spectroscopy, we shed light on the underlying physics—including dephasing, inhomogeneity, and strain—for a MoSe 2 monolayer and identify both promising and unfavorable areas for quantum information applications. We, furthermore, apply the same technique to a MoSe 2 /WSe 2 heterostructure. Despite the notable presence of strain and dielectric environment changes, coherent and incoherent coupling and interlayer exciton lifetimes are mostly robust across the sample. This uniformity is despite a significantly inhomogeneous interlayer exciton photoluminescence distribution that suggests a bad sample for device applications. This robustness strengthens the case for TMDs as a next-generation material platform in quantum information science and beyond.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Macroscopic trends of linear tearing stability in cylindrical current profiles

Abstract The likelihood of realising tokamak power-plants will be greatly improved by the discovery of high-gain equilibria that resist the formation of small islands and hence avoid the disruptive neoclassical tearing mode. We propose a series of studies to understand how simple tokamak design can leverage aspects of tearing onset physics to maximise passive resistance to island formation. Here we investigate the variation that current profiles can bring about in preventing tearing onset through the cylindrical linear tearing stability parameter Δ ′ . A database of 159148 realistic pilot-plant current profiles was generated with Monte Carlo sampling, and the distribution of Δ ′ values was linked with interpretable profile characteristics. In agreement with prior theoretical and experimental studies, Δ ′ was found to be strongly correlated with the existence and steepness of a local toroidal current well or hill, with the former destabilising and the latter stabilising. In the absence of these two cases, the remaining Δ ′ values were linearly bounded by the toroidal current gradient at the rational surface.

Physics↗

Galaxy blending effects in deep imaging cosmic shear probes of cosmology

ABSTRACT Upcoming deep imaging surveys such as the Vera C. Rubin Observatory Legacy Survey of Space and Time will be confronted with challenges that come with increased depth. One of the leading systematic errors in deep surveys is the blending of objects due to higher surface density in the more crowded images; a considerable fraction of the galaxies which we hope to use for cosmology analyses will overlap each other on the observed sky. In order to investigate these challenges, we emulate blending in a mock catalogue consisting of galaxies at a depth equivalent to 1.3 yr of the full 10-yr Rubin Observatory that includes effects due to weak lensing, ground-based seeing, and the uncertainties due to extraction of catalogues from imaging data. The emulated catalogue indicates that approximately 12 per cent of the observed galaxies are ‘unrecognized’ blends that contain two or more objects but are detected as one. Using the positions and shears of half a billion distant galaxies, we compute shear–shear correlation functions after selecting tomographic samples in terms of both spectroscopic and photometric redshift bins. We examine the sensitivity of the cosmological parameter estimation to unrecognized blending employing both jackknife and analytical Gaussian covariance estimators. An ∼0.025 decrease in the derived structure growth parameter S8 = σ8(Ωm/0.3)0.5 is seen due to unrecognized blending in both tomographies with a slight additional bias for the photo-z-based tomography. This bias is greater than the 2σ statistical error in measuring S8.

79 ASTRONOMY AND ASTROPHYSICS↗

pop-cosmos : redshifts and physical properties of KiDS-1000 galaxies

ABSTRACT Principled Bayesian inference of galaxy properties has not previously been performed for wide-area weak-lensing surveys with millions of sources. We address this gap by applying the pop-cosmos generative model to perform spectral energy distribution (SED) fitting for 4 million KiDS (Kilo-Degree Survey)-1000 galaxies. Calibrated on deep COSMOS2020 photometric data, pop-cosmos specifies a physically motivated prior over the galaxy population up to $z \simeq 6$ in stellar population synthesis (SPS) parameter space. Using the Speculator SPS emulator with GPU (graphics processing unit)-accelerated Markov Chain Monte Carlo sampling, we perform full posterior inference at 8.2 GPU seconds per galaxy, obtaining joint constraints on galaxy redshifts and physical properties. We validate photometric redshifts against $\sim \!185\,\!000$ KiDS galaxies cross-matched to Dark Energy Spectroscopic Instrument Data Release 1 spectroscopic samples, achieving low bias ($2\times 10^{-3}$), scatter ($\sigma _{\mathrm{MAD}}=0.03$), and outlier fraction (3.2 per cent) for the Bright Galaxy Survey, with comparable performance (bias $3\times 10^{-2}$, $\sigma _{\mathrm{MAD}}=0.05$, 1.0 per cent outliers) for luminous red galaxies (LRGs). Within the LRG sample, we identify massive, dusty, star-forming contaminants at $z \simeq 0.4$ satisfying standard colour selections for quenched populations. We infer trends in stellar mass, star formation, metallicity, and dust across five tomographic redshift bins consistent with established scaling relations. Using specific star formation rate constraints, we identify $\sim$7 per cent of KiDS-1000 galaxies as quenched, versus 37 per cent implied by conservative colour cuts. This enables the construction of weak-lensing samples defined by physical properties while mitigating intrinsic alignment systematics and preserving statistical power. Our analysis validates pop-cosmos out of sample, establishing it as a scalable approach for galaxy evolution and cosmological analyses with photometric surveys.

Halder, Anik [Institute of Astronomy and Kavli Ins↗

Sparse expansions of multicomponent oxide configuration energy using coherency and redundancy

We report that compressed sensing has become a widely accepted paradigm to construct high dimensional cluster expansion models used for statistical mechanical studies of atomic configuration in complex multicomponent crystalline materials. However, strict sampling requirements necessary to obtain minimal coherence measurements for compressed sensing to guarantee accurate estimation of model parameters are difficult and in some cases impossible to satisfy due to the inability of physical systems to access certain configurations. Nevertheless, the dependence of energy on atomic configuration can still be adequately learned without these strict requirements by using compressed sensing by way of coherent measurements using redundant function sets known as frames. We develop a particular frame constructed from the union of all occupancy-based cluster expansion basis sets. We illustrate how using this highly redundant frame yields sparse expansions of the configuration energy of complex oxide materials that are competitive and often surpass the prediction accuracy and sparsity of models obtained from standard cluster expansions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Kekulé valence bond order in the honeycomb lattice optical Su-Schrieffer-Heeger model and its relevance to graphene

We perform sign-problem-free determinant quantum Monte Carlo simulations of the optical Su- Schrieffer-Heeger model on a half-filled honeycomb lattice. In particular, we investigate the model’s semi-metal (SM) to Kekulé Valence Bond Solid (KVBS) phase transition at zero and finite temper- atures as a function of phonon energy and interaction strength. Using hybrid Monte Carlo sampling methods we can simulate the model near the adiabatic regime, allowing us to access regions of parameter space relevant to graphene. Our simulations suggest that the SM-KVBS transition is weakly first-order at all temperatures, with graphene situated close to the phase boundary in the SM region of the phase diagram. Furthermore, our results highlight the important role bond-stretching phonon modes play in the formation of KVBS order in strained graphene-derived systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Beyond Fisher forecasting for cosmology

The planning and design of future experiments rely heavily on forecasting to assess the potential scientific value provided by a hypothetical set of measurements. The Fisher information matrix, due to its convenient properties and low computational cost, provides an especially useful forecasting tool. However, the Fisher matrix only provides a reasonable approximation to the true likelihood when data are nearly Gaussian distributed and observables have nearly linear dependence on the parameters of interest. Also, Fisher forecasting techniques alone cannot be used to assess their own validity. Thorough sampling of the exact or mock likelihood can definitively determine whether a Fisher forecast is valid, though such sampling is often prohibitively expensive. Here we propose a simple test, based on the Derivative Approximation for likelihoods (DALI) technique, to determine whether the Fisher matrix provides a good approximation to the exact likelihood. We show that the Fisher matrix becomes a poor approximation to the true likelihood in regions where two-dimensional slices of level surfaces of the DALI approximation to the likelihood differ from two-dimensional slices of level surfaces of the Fisher approximation to the likelihood. We demonstrate that our method accurately predicts situations in which the Fisher approximation deviates from the true likelihood for various cosmological models and several data combinations, with only a modest increase in computational cost compared to standard Fisher forecasts.

79 ASTRONOMY AND ASTROPHYSICS↗

Modest set of simulations of local-type primordial non-Gaussianity in the DESI era

A measurement of a primordial non-Gaussianity (PNG) signal through late- or early-Universe probes has the potential to transform our understanding of the physics of the primordial Universe. While large-scale structure observables in principle contain vital information, interpreting these measurements is challenging due to poorly understood astrophysical effects. Luckily, N-body simulations, such as the abacuspng set presented in this study, consisting of 9 boxes, each of size L box = 2 Gpc/h and particle mass of 1.01 × 10 10 M ⊙ /h, provide a viable path forward. As validation, we find good agreement between the simulations and our expectations from one-loop perturbation theory (PT) and the "separate universe" method for the matter bispectrum, matter power spectrum and the halo bias parameter associated with PNG, b φ . As a science application, here we investigate the link between halo assembly bias and b φ for halo properties known to play a vital role in accurately predicting galaxy clustering: concentration, shear (environment), and accretion rate. We find a strong response for all three parameters, suggesting that the connection between b φ and the assembly history of halos needs to be taken into account by future PNG analyses. We further perform the first study of the b φ parameter from fits to early DESI data of the luminous red galaxy (LRG) and quasistellar object (QSO) samples and comment on the effect on f NL constraints for the allowed galaxy-halo models (note that σ[f NL ]∝ $\frac{σ[b_φ]}{b_φ}$). We find that the error on f NL is 15, 8, 7 for the LRGs at z = 0.5 and z = 0.8 and QSOs at z = 1.4, respectively, suggesting that a thorough understanding of galaxy assembly bias is warranted so as to perform robust high-precision analysis of local-type PNG with future surveys.

79 ASTRONOMY AND ASTROPHYSICS↗

Efficiently Verifiable Quantum Advantage on Near-Term Analog Quantum Simulators

Existing schemes for demonstrating quantum computational advantage are subject to various practical restrictions, including the hardness of verification and challenges in experimental implementation. Meanwhile, analog quantum simulators have been realized in many experiments to study novel physics. In this work, we propose a quantum advantage protocol based on verification of an analog quantum simulation, in which the verifier need only run an O ( λ 2 ) -time classical computation, and the prover need only prepare O ( 1 ) samples of a history state and perform O ( λ 2 ) single-qubit measurements, for a security parameter λ . We also propose a near-term feasible strategy for honest provers and discuss potential experimental realizations. Published by the American Physical Society 2025

Liu, Zhenning (ORCID:000000020794419X)↗

Highly Extensible X-ray Diffraction Toolkit

This code has been previous released as open source under LLNL-CODE-529294. The HEXRD program comprises a library of X-ray diffraction analysis tools and a stand-alone GUI application. While there are many software packages that provide similar capabilities, HEXRD was designed to provide an extremely flexible and extensible basis for describing a broad array of X-ray diffraction instruments; this includes both mono- and poly-chromatic modalities as well as poly- and single-crystal samples. The high-level objects include Material: an abstraction for defining crystalline materials, including unit cell parameters and anisotropic elastic moduli Imageseries: an interface for series of diffraction images containing one or more frames, including provisions for metadata. Instrument: an abstraction of an X-ray diffraction instrument including support for the definition of multiple detector elements as well as powder, Laue, and mono-chromatic single/multi-grain diffraction analysis modalities (including attendant calibration routines). The hexrdgui package is a standalone PySide2/Qt5 GUI built on top of hexrd. The 0.8.1 release includes various bug fixes as well as some new functionality in the GUI, including: o Enhancements to the import tools for LLNL diffraction instruemtns (PXRDIP and TARDIS) o The ability to save 1-d LeBail fits of integrated diffraction spectra.

Barton, NathanR.↗

The Precision of Product Consistency Tests with Environmental Assessment Glass

Several series of triplicate tests conducting following the ASTM C1285 product consistency test method have been performed using samples of Environmental Assessment glass in Teflon™ and steel test vessels. The results of tests with test parameters within the ranges specified in the PCT method A procedure were used to calculate repeatability and reproducibility precision statistics for tests in each type of vessel. These show the test response and test precision attained in tests conducted in Teflon™ and stainless steel test vessels are statistically indistinguishable.

36 MATERIALS SCIENCE↗

Validation of the Public Radiation Exposure Calculation for the Incident at the National Institute of Standards and Technology Center for Neutron Research on February 3, 2021

The Department of Energy/National Nuclear Security Administration Consequence Management Program was contacted by the Health Physics Chief of the National Institute of Standards and Technology Center for Neutron Research (NCNR) to review public radiation exposure calculations for an event that occurred on its Gaithersburg, Maryland, campus on February 3, 2021. Subject matter experts from the Nuclear Emergency Support Team (NEST) assets, specifically the Consequence Management Home Team (CMHT) and the National Atmospheric Release Advisory Center (NARAC), were selected to provide support. CMHT used three separate modeling codes to validate the results the scientists at NCNR calculated using the HotSpot model. The analyses were performed using NARAC’s in-house Lagrangian dispersion codes known as LODI and Aeolus, as well as the Turbo FRMAC software from Sandia National Laboratories. The team used parameters provided by the NCNR scientists regarding the site, applicable observable meteorological data, and environmental survey and sampling data to estimate public exposure. Each model estimated public dose at much less than 0.5 mrem. CMHT concurs with the NCNR public radiation exposure calculations which state that members of the public at the 400-meter boundary would have received a radiological dose of less than 0.5 mrem.

61 RADIATION PROTECTION AND DOSIMETRY↗

3D Printing of Cement-Based Materials Using Seawater for Simulated Marine Environments

Global demand for adaptable and rapidly deployable construction solutions in offshore, coastal, and fluvial environments continues to rise, driven by pressing needs to develop energy platforms, improve coastal resilience, and support emergency response in the face of natural disasters. Increased investment in human-made coastal infrastructure, such as piers, support structures for power lines, offshore wind farms, and seawall protection systems, further underscores this trend. This study investigates the development of printable concrete mixtures for underwater environments using seawater as a replacement for freshwater, using a 3D printing syringe-based extrusion system. The effect of seawater addition and the printing medium (in air vs. underwater) was assessed via rheological and mechanical performance characterization. The results indicate rheological properties are favorable for seawater adoption by producing mixtures with higher yield stress and viscosity with the same levels of admixtures used for freshwater. Seawater-based mixtures demonstrated superior dimensional stability compared to freshwater counterparts, maintaining cross-sectional geometry, while compressive strength results showed no statistical differences between in-air and underwater samples. However, flexural strength was significantly influenced by geometry and printing medium. These findings establish critical rheological parameters for printable underwater mixtures and highlight the need for optimized curing strategies and layer bonding techniques to improve interfacial strength in underwater 3D printing applications.

36 MATERIALS SCIENCE↗

Rejection Sampling with Autodifferentiation -- Case study: Fitting a Hadronization Model

We present an autodifferentiable rejection sampling algorithm termed Rejection Sampling with Autodifferentiation (RSA). In conjunction with reweighting, we show that RSA can be used for efficient parameter estimation and model exploration. Additionally, this approach facilitates the use of unbinned machine-learning-based observables, allowing for more precise, data-driven fits. To showcase these capabilities, we apply an RSA-based parameter fit to a simplified hadronization model.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Temperature dependence and limiting mechanisms of the upper critical field of FeSe thin films

We use magnetoresistance measurements at high magnetic field (μ 0⁢ H ≤ 65 T) and low temperature (T ≥ 500 mK) to gain fresh insights into the behavior of the upper critical field H c⁢2 in superconducting ultrathin FeSe films of varying degrees of disorder, grown by molecular beam epitaxy on SrTiO 3 . Measurements of H c⁢2 across samples with a widely varying superconducting critical temperature (1.2 K ≤ T c ≤ 21 K) generically show similar qualitative temperature dependence. Here, we analyze the temperature dependence of H c⁢2 in the context of Werthamer-Helfand-Hohenberg (WHH) theory. The analysis yields parameters that indicate a strong Pauli paramagnetic pair-breaking mechanism which is also reflected by pseudoisotropic superconductivity in the limit of zero temperature. In the lower T c samples, we observe a spin-orbit scattering-driven enhancement of H c⁢2 above the strongly-coupled Pauli paramagnetic limit. We also observe clear deviations from WHH theory at low temperature, regardless of T c . We attribute this to the multiband superconductivity of FeSe and possibly to the emergence of a low-temperature, high-field superconducting phase.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Risk Importance Ranking of Fire Data Parameters to Enhance Fire PRA Model Realism

Fire is historically and analytically a significant contributor to nuclear power plant risk. The level of fire risk and the methods, tools and data for modeling this risk is highly debated by experts. One area of debate is the input data used in fire modeling and how to deal with this data’s high uncertainty. This report outlines initial work performed for determining the key parameters causing this uncertainty and how it propagates into nuclear power plant models. This research paves the way for the development of methods to reduce fire data uncertainty used in modeling. The Nuclear Regulatory Commission has mandated that nuclear power plants perform fire risk modeling. However, there are several issues with the current risk modeling implementation that affect the results. Approved modeling methods can be overly conservative and often do not match plant experience. Also, the data used in the modeling can have high uncertainties and is influenced by expert judgement. To evaluate input data uncertainty, researchers performed an initial review of several fire experiments done at Sandia National Laboratories. Uncertainties for fire data can come from many sources, such as experiment design constraints, environmental conditions, or other plant-specific aspects. There are many different significant and insignificant parameters driving the uncertainty. Additionally, the uncertainty of the different input data used in the fire modeling could have a significant or insignificant effect on the entire plant risk. A four-step methodology was developed to perform Integrated Probabilistic Risk Assessment Importance Ranking. A demonstration case using these steps was set up and three of the four steps were completed in fiscal year (FY) 2019 and the fourth step done FY 2020. These steps are: 1. The qualitative analysis of potential sources was conducted with the following items identified for the demonstration. • Maximum heat release rate • Time to maximum heat release rate • Duration of max heat release rate • Time to decay • Thermal conductivity of concrete • Specific heat of concrete • Density of concrete • Cable jacket thickness 2. A quantitative characterization of dominant sources of uncertainty was performed. A list of distributions and determined values of the dominant sources is shown in Appendix A. 3. A quantitative screening of the potential sources of uncertainty using Morris Elementary Effects Analysis was completed. An experimental model using the physics-based fire modeling tool Fire Dynamics Simulator was developed and coupled with the Risk Analysis Virtual Environment. The Morris analysis identified at least two parameters that can be eliminated as significant contributors (specific heat of concrete and cable jacket thickness). 4. Global importance measure (Global IM) analysis to generate a comprehensive ranking based on their influence on the plant risk. In this research, a moment-independent Global IM is used since it can address (a) uncertainty in the input parameters of the fire model, (b) uncertainty in the risk outputs, and (c) non-linearity and interactions among input parameters in the fire model, more accurately than the correlation-based and variance-based global methods. The observations from the research showed that, depending on the initial and boundary conditions of the fire scenarios, fire-induced damage could have a very small probability and could be dominated by the tail of the uncertainty distribution; hence, the accuracy of the correlation-based and variance-based methods is questionable. The moment-independent Global IM analysis in this research provides a better understanding of how experimental uncertainty data affects industry’s plant models and where improvements in that data will have the largest benefit for improving fire modeling accuracy in causing core damage. Among the five unscreened parameters obtained from the Morris EE analysis, the Global IM analysis results for the case study indicated that max heat release rate and fire location are the most important parameters. The report also outlines benefits of using a unified computational platform that integrates the underlying simulations (e.g., a fire progression model), quantitative screening (using the Morris EE method), and the Global IM analysis. A unified platform can (i) facilitate the ranking of input parameters considering multiple key fire scenarios simultaneously, rather than considering one scenario at a time, (ii) contribute to more explicit and accurate treatment of dependencies at multiple levels of Fire PRA, (iii) facilitate the sampling-based uncertainty quantification for Fire PRA, and (iv) help generating both “industry-wide” and "plant-specific" ranking of uncertainty sources in Fire PRA. Future research should be done to include additional parameters such as detection/suppression or cable fire spread. Adding a PRA software such as SAPHIRE to the RAVEN platform would help with plant model integration and improve treatment of fire-induced dependency. The I-PRA risk importance ranking methodology offered in this report can provide valuable information for efficiently (a) enhancing the realism of Fire PRA for existing plants and (b) supporting the development of Dynamic Fire PRA for advanced reactors and new plants.

97 MATHEMATICS AND COMPUTING↗

DAmodel: hierarchical Bayesian modelling of DA white dwarfs for spectrophotometric calibration

We use hierarchical Bayesian modelling to calibrate a network of 32 all-sky faint DA white dwarf (DA WD) spectrophotometric standards (⁠16.5 < V , 19.5⁠) alongside three CALSPEC standards, from 912 Å to 32 μm. The framework is the first of its kind to jointly infer photometric zero points and WD parameters (surface gravity log g⁠, effective temperature T eff ⁠, extinction A V ⁠, dust relation parameter R V ) by simultaneously modelling both photometric and spectroscopic data. We model panchromatic Hubble Space Telescope Wide Field Camera 3 (HST/WFC3) UVIS and IR photometry, HST/STIS UV spectroscopy, and ground-based optical spectroscopy to sub-per cent precision. Photometric residuals for the sample are the lowest yet yielding < 0.004 mag RMS on average from the UV to the NIR, achieved by jointly inferring time-dependent changes in system sensitivity and WFC3/IR count-rate nonlinearity. Our GPU-accelerated implementation enables efficient sampling via Hamiltonian Monte Carlo, critical for exploring the high-dimensional posterior space. The hierarchical nature of the model enables population analysis of intrinsic WD and dust parameters. Inferred spectral energy distributions from this model will be essential for calibrating the James Webb Space Telescope as well as next-generation surveys, including Vera Rubin Observatory’s Legacy Survey of Space and Time and the Nancy Grace Roman Space Telescope.

methods: statistical↗