Search NASA⌕ Search

SEARCH · Search NASA

Results for “calibration and fitting 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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Ultra-sensitive radon assay using an electrostatic chamber in a recirculating system

Rare event searches such as neutrinoless double beta decay and Weakly Interacting Massive Particle detection require ultra-low background detectors. Radon contamination is a significant challenge for these experiments, which employ highly sensitive radon assay techniques to identify and select low-emission materials. This work presents the development of ultra-sensitive electrostatic chamber (ESC) instruments designed to measure radon emanation in a recirculating gas loop, for future lower background experiments. Unlike traditional methods that separate emanation and detection steps, this system allows continuous radon transport and detection. This is made possible with a custom-built recirculation pump. A Python-based analysis framework, PyDAn, was developed to process and fit time-dependent radon decay data. Radon emanation rates are given for various materials measured with this instrument. A radon source of known activity provides an absolute calibration, enabling statistically-limited minimal detectable activities of 20 µBq. These devices are powerful tools for screening materials in the development of low-background particle physics experiments.

47 OTHER INSTRUMENTATION↗
DOE OSTI2025

Ultra-sensitive radon assay using an electrostatic chamber in a recirculating system

Rare event searches such as neutrinoless double beta decay and Weakly Interacting Massive Particle detection require ultra-low background detectors. Radon contamination is a significant challenge for these experiments, which employ highly sensitive radon assay techniques to identify and select low-emission materials. This work presents the development of ultra-sensitive electrostatic chamber (ESC) instruments designed to measure radon emanation in a recirculating gas loop, for future lower background experiments. Unlike traditional methods that separate emanation and detection steps, this system allows continuous radon transport and detection. This is made possible with a custom-built recirculation pump. A Python-based analysis framework, PyDAn, was developed to process and fit time-dependent radon decay data. Radon emanation rates are given for various materials measured with this instrument. A radon source of known activity provides an absolute calibration, enabling statistically-limited minimal detectable activities of 20 uBq. These devices are powerful tools for screening materials in the development of low-background particle physics experiments.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2025

Searches for New Physics With Muon Conversion at Fermilab and Triboson Production at the LHC

We report on several efforts to search for physics beyond the standard model of particle physics at broad energy scales. The Mu2e experiment at Fermilab will search for charged lepton flavor violation via the muon to electron conversion process, which is suppressed in the Standard Model. Mu2e will be operated at a low energy, yet can probe New Physics at very high mass scales (O(1e3 - 1e4 ) TeV). At high energies, the CMS experiment at the CERN LHC continues to deliver an impressive suite of Standard Model measurements and limits on a variety of New Physics signatures. Mu2e is under construction and slated to collect its first physics data in the coming years. This thesis describes work done during the construction phase of Mu2e and focuses on two critical areas: magnetic field modeling and statistical analysis. We describe a novel method for field modeling which we validate using a simulated dataset representing the expected magnetic field in the Detector Solenoid. This method blends a standard least-squares fitting technique that utilizes physically motivated analytical model functions with a novel physics informed network that is constructed to obey Maxwell’s equations. We show the technique can model the field with an accuracy of 10−7 despite the presence of injected noise in the pseudo-measurements at the 10−5 level. We then present preliminary results of the calibration of 3D Hall probes at the sub-10−4 level. These probes will be used to directly measure the Mu2e Detector Solenoid magnetic field on a sparse grid; these measurements serve as the input to the field model fitting. Finally, we describe the first implementation of both an unbinned shape analysis and a Bayesian interpretation applied to Mu2e pseudo-data. Up to 20% tighter limits can be set by the shape analysis compared to a standard cut & count analysis. The AlCap experiment collected data at PSI in 2015 to measure several important quantities related to nuclear muon capture on an aluminum target, which is a significant background process for Mu2e. The neutron emission from muon capture can introduce background hits in the Mu2e detectors and can increase radiation damage in various elements of the apparatus. We present measurements of the neutron group fluence and mean neutron multiplicity for muon capture on aluminum nuclei. Finally, we discuss an analysis of triboson production at CMS using an Effective Field Theory framework. Standard Model triboson production, which was first observed at CMS in 2020, has a relatively small cross section and provides direct access to both anomalous triple gauge couplings and quartic gauge couplings. These couplings, interpreted in the Standard Model Effective Field Theory, are studied in the present work. We target the boosted regime where the background rate is low and yields are enhanced when dimension-6 and dimension-8 Wilson coefficients are non-zero. We do not observe an excess in the data and therefore set bounds on the Wilson coefficients. For dimension-6 coefficients the tightest observed (expected) bounds are set on cW /Λ2 where Λ is the mass scale of new physics; the bounds are [−0.13, 0.12] TeV−2 ([−0.12, 0.12] TeV−2 ) at 95% CL. The tightest bounds in dimension-8 are set on fT,0 / Λ4 ; the observed (expected) bounds at 95% CL are [−0.63, 0.69] TeV−4 ([−0.54, 0.62] TeV−4 ). Additional results are presented which include scenarios where multiple Wilson coefficients are non-zero, the application of signal model clipping to address unitarity violation in Effective Field Theories, and a novel template fit developed for easier reinterpretation of our results.

Kampa, Cole Erik [Northwestern U. (main)] (ORCID:0↗
DOE OSTI2025

Current interrupt method for calculating the electrochemical impedance in a solid oxide electrolysis stack

Here, in this work the time domain response of Solid Oxide Electrolysis Cells (SOEC) to a current interruption was transformed into the frequency domain using a carrier function Laplace transform, which is fit to the experimental data using a MATLAB Complex Nonlinear Least Squares (CNLS) solver. The hardware implementation, consisting principally of a high-speed switch and a fast-logging Analog to Digital Converter (ADC), was assembled and tested using a calibration module to assess the accuracy, repeatability, and speed of acquisition of the prototype device as compared against a calibrated commercial impedance spectrometer. Additionally, the current interrupt device and commercial FRA were used to acquire the impedance spectra of a four cell SOEC stack with a large, 300 cm 2 , active cell area.

SOEC↗
DOE OSTI2024

An integrated approach to derive relative permeability from capillary pressure

Surface tension affects all aspects of fluid flow in porous media. Through measurements of surface tension interaction under multiphase conditions, a relative permeability curve can be determined. Relative permeability is a numerical description of the interaction between two or more fluids and the porous media. It is a critical parameter for various tools that characterize subsurface multiphase flow systems, such as numerical simulation for carbon sequestration, oil and gas development, and groundwater contamination remediation. Therefore, it is critical to get a good statistical distribution of relative permeability in the porous media under study. Empirical formula for determining relative permeability from capillary pressure are already well established but do not provide the needed flexibility that is required to match laboratory-derived relative permeability curves. By expanding the existing methods for calculating relative permeability from capillary pressure data, it is possible to create both two and three-phase relative permeability curves. Mercury intrusion capillary pressure (MICP) data from the Morrow 'B' Sandstone coupled with interfacial tension and contact angle measurements were used to create a suite of relative permeability curves. Furthermore, these curves were then calibrated to a small sample of existing laboratory curves to elucidate common fitting parameters for the formation that were then used to create relative permeability curves from MICP data that does not have an associated laboratory-measured relative permeability curve.

58 GEOSCIENCES↗
NASA NTRS2000

Phase Screen Determination for the GeoSAR Interferometric Mapping Instrument

GeoSAR is a dual frequency airborne radar interferometric mapping instrument designed to obtain high resolution and high accuracy digital elevation models. The X-band inteferometric antennas are mounted beneath the fuselage with a separation distance of 2.6 m. Midway between the antennas along the aircraft centerline is a large fairing housing the Laser Baseline metrology System (LBMS) and Inertial Navigtion Units (INUs). The P-band antennas are mounted in wingtip pods approximately 10 m from the aircraft centerline. Multipath signals reflected from the LBMS fairing cause phase distortions corrupting the topograhic height measurements for both the X-band and P-band systems. In addition multipath off the wing affects the P-band phase measurements. In this paper we discuss a technique used to mitigate the amount of X-Band multipath and methods calibrating and correcting the phase for both X-Band and P-Band using a high accuracy DEM. Height differences between the interferometrically derived elevations and a photogrammetrically derived DEM are converted into phase differences. Multipath phase errors are sinusoidal with spatial frequency and amplitude related to the magnitude of reflection of the multipath reflection point and its distance from the interferometric antennas. These phase residuals are fitted using a Chebyshev polynomial sequence and used as a phase screen in the processor to correct the phase prior to height reconstruction. Based on the correction signal it is possible to estimate the location and magnitude of the multipath signal. This information can then be used to extend the domain of the phase screen over mapping angles not covered during the calibration passes. Keywords: SAR,interferometry,phase screen

Hensley, Scott↗
NASA NTRS2019

On the Use of Optimization Techniques for Turbulence Model Calibration

This article discusses the use of numerical optimization procedures to aid in the calibration of turbulence model coefficients. Such methods would increase the rigor and repeatability of the calibration procedure by requiring clearly defined and objective optimization metrics, and could be used to identify unique combinations of coefficient values for specific flow problems. The approach is applied to the re-calibration of an explicit algebraic Reynolds stress model for the incompressible planar mixing layer using the Nelder-Mead simplex algorithm and a micro-genetic algorithm with minimally imposed constraints. Three composite fitness functions, each based upon the error in the mixing layer growth rate and the normal and shear components of the Reynolds stresses, are investigated. The results demonstrate a significant improvement in the target objectives through the adjustment of three pressure-strain coefficients. Adjustments of additional coefficients provide little further benefit. Issues regarding the effectiveness of the fitness functions and the efficiency of the optimization algorithms are also discussed.

Yoder, Dennis A.↗
NASA NTRS2002

Participation in the TOMS Science Team

Because of the nominal funding provided by this grant, some of the relevant research is partially funded by other sources. Research performed for this funding period included the following items: We have investigated errors in TOMS ozone measurements caused by the uncertainty in wavelength calibration, coupled with the ozone cross sections in the Huggins bands and their temperature dependence. Preliminary results show that 0.1 nm uncertainty in TOMS wavelength calibration at the ozone active wavelengths corresponds to approx. 1% systematic error in O3, and thus potential 1% biases among ozone trends from the various TOMS instruments. This conclusion will be revised for absolute O3 Measurements as cross sections are further investigated for inclusion in the HITRAN database at the SAO, but the potential for relative errors remains. In order to aid further comparisons among TOMS and GOME ozone measurements, we have implemented our method of direct fitting of GOME radiances (BOAS) for O3, and now obtain the best fitting precision to date for GOME O3 Columns. This will aid in future comparisons of the actual quantities measured and fitted for the two instrument types. We have made comparisons between GOME ICFA cloud fraction and cloud fraction determined from GOME data using the Ring effect in the Ca II lines. There is a strong correlation, as expected, but there are substantial systematic biases between the determinations. This study will be refined in the near future using the recently-developed GOME Cloud Retrieval Algorithm (GOMECAT). We have improved the SAO Ring effect determination to include better convolution with instrument transfer functions and inclusion of interferences by atmospheric absorbers (e.g., O3). This has been made available to the general community.

Chance, Kelly↗
NASA NTRS2008

Calibrating the Helium Pressurization System for the Space Shuttle Liquid-Hydrogen Tank

Analysis of the results from the STS-114 tanking tests and subsequent launch called into question existing thermal and mass models of helium pressurization of the liquid hydrogen tank. This hydrogen tank, which makes up the bottom two-thirds of the External Tank, is pressurized prior to launch to avoid cavitation in the Shuttle Main Engine pumps. At about 2 minutes prior to launch, the main vent valve is closed, and pressurized helium flows into the tank ullage space to achieve set point pressure. As the helium gas cools, its pressure drops, calling for additional helium. Subsequent helium flows are provided in short, timed pulses. The number of pulses is taken as a rough leak indicator. An analysis of thermal models by Marshall Space Flight Center showed considerable uncertainty in the pressure-versus-time behavior of the helium ullage space and the ability to predict the number of pulses normally expected. Kennedy Space Center proposed to calibrate the dime-sized orifice, which together with valves, controls the helium flow quantity (Figure 1). Pressure and temperature sensors were installed to provide upstream and downstream measurements necessary to compute flow rate based on the orifice discharge coefficient. An assessment of flow testing with helium indicated an extremely costly use of this critical resource. In order to reduce costs, we proposed removing the orifices from each Mobile Launcher Platform (MLP) and asking Colorado Engineering Experiment Station Inc. (CEESI) to calibrate the flow. CEESI has a high-pressure air flow system with traceable flow meters capable of handling the large flow rates. However, literature research indicated that square-edged orifices of small diameters often exhibit significant hysteresis and nonrepeatability in the vicinity of choked or sonic flow. Fortunately, the MLP orifices behaved relatively well in testing (Figure 2). Using curve fitting of the air-flow data, in conjunction with ASME orifice modeling equations, a method of relating the helium mass flow to measured air flow data was obtained. This analysis showed that the highest uncertainty in flow occurred in the vicinity of the choking pressure ratio, as would be expected. In addition, analysis of typical flow pulses showed that most of the helium flow occurred either well below or well above this uncertain area. The final result is the ability to provide postlaunch estimates of helium mass flows that are within 1.5 percent of the actual value.

Source record↗
NASA NTRS1996

Land Surface Temperature Measurements form EOS MODIS Data

We have developed a physics-based land-surface temperature (LST) algorithm for simultaneously retrieving surface band-averaged emissivities and temperatures from day/night pairs of MODIS (Moderate Resolution Imaging Spectroradiometer) data in seven thermal infrared bands. The set of 14 nonlinear equations in the algorithm is solved with the statistical regression method and the least-squares fit method. This new LST algorithm was tested with simulated MODIS data for 80 sets of band-averaged emissivities calculated from published spectral data of terrestrial materials in wide ranges of atmospheric and surface temperature conditions. Comprehensive sensitivity and error analysis has been made to evaluate the performance of the new LST algorithm and its dependence on variations in surface emissivity and temperature, upon atmospheric conditions, as well as the noise-equivalent temperature difference (NE(Delta)T) and calibration accuracy specifications of the MODIS instrument. In cases with a systematic calibration error of 0.5%, the standard deviations of errors in retrieved surface daytime and nighttime temperatures fall between 0.4-0.5 K over a wide range of surface temperatures for mid-latitude summer conditions. The standard deviations of errors in retrieved emissivities in bands 31 and 32 (in the 10-12.5 micrometer IR spectral window region) are 0.009, and the maximum error in retrieved LST values falls between 2-3 K. Several issues related to the day/night LST algorithm (uncertainties in the day/night registration and in surface emissivity changes caused by dew occurrence, and the cloud cover) have been investigated. The LST algorithms have been validated with MODIS Airborne Simulator (MAS) dada and ground-based measurement data in two field campaigns conducted in Railroad Valley playa, NV in 1995 and 1996. The MODIS LST version 1 software has been delivered.

Wan, Zhengming↗
DOE OSTI2025

Results of an interlaboratory study on the working curve in vat photopolymerization II: Towards a standardized method

The working curve measurement in photopolymer additive manufacturing is a ubiquitous measure of the cure depth of a printing resin as a function of radiant exposure of light. The fit parameters from this measurement (the depth of light penetration D p and the critical exposure E c ) are used to evaluate and report a resin’s printability, optimize processing parameters, and inform print and resin quality control. Despite its widespread use, the working curve lacks a standard measurement method. Here, following up on our paper “Results of an Interlaboratory Study on the Working Curve in Vat Photopolymerization” from last year, an interlaboratory study on the working curve was performed using calibrated, reproducible, bandpass filtered light sources. With these light sources, the variability between labs in measured working curves was dramatically reduced from the initial interlaboratory study. Aggregate data from this experiment produced reliable D p and E c measurements at 385 nm of 39.2 ± 3.7 µm and 12.3 ± 3.0 mJ cm −2 , respectively. At 405 nm the values of D p and E c are 69.3 ± 3.8 µm and 17.9 ± 2.3 mJ cm −2 , respectively. The results are agnostic to the thickness measurement tool utilized by participants, ensuring broad applicability across laboratories. We also tested the generalizability of the proposed method of using a filtered light source by filtering a commercial 405 nm light source and obtaining a working curve in agreement with the aggregate data from the interlaboratory study. This interlaboratory study provides a basis for a documentary standard for the working curve, so that the entire photopolymer additive manufacturing industry can share reproducible and interoperable working curve data.

36 MATERIALS SCIENCE↗
NASA NTRS2014

Strain Gage Load Calibration of the Wing Interface Fittings for the Adaptive Compliant Trailing Edge Flap Flight Test

This is the presentation to follow conference paper of the same name. The adaptive compliant trailing edge (ACTE) flap experiment safety of flight requires that the flap to wing interface loads be sensed and monitored in real time to ensure that the wing structural load limits are not exceeded. This paper discusses the strain gage load calibration testing and load equation derivation methodology for the ACTE interface fittings. Both the left and right wing flap interfaces will be monitored and each contains four uniquely designed and instrumented flap interface fittings. The interface hardware design and instrumentation layout are discussed. Twenty one applied test load cases were developed using the predicted in-flight loads for the ACTE experiment.

strain gage calibration↗
NASA NTRS2003

First Year Wilkinson Microwave Anisotropy Probe(WMAP) Observations: Data Processing Methods and Systematic Errors Limits

We describe the calibration and data processing methods used to generate full-sky maps of the cosmic microwave background (CMB) from the first year of Wilkinson Microwave Anisotropy Probe (WMAP) observations. Detailed limits on residual systematic errors are assigned based largely on analyses of the flight data supplemented, where necessary, with results from ground tests. The data are calibrated in flight using the dipole modulation of the CMB due to the observatory's motion around the Sun. This constitutes a full-beam calibration source. An iterative algorithm simultaneously fits the time-ordered data to obtain calibration parameters and pixelized sky map temperatures. The noise properties are determined by analyzing the time-ordered data with this sky signal estimate subtracted. Based on this, we apply a pre-whitening filter to the time-ordered data to remove a low level of l/f noise. We infer and correct for a small (approx. 1 %) transmission imbalance between the two sky inputs to each differential radiometer, and we subtract a small sidelobe correction from the 23 GHz (K band) map prior to further analysis. No other systematic error corrections are applied to the data. Calibration and baseline artifacts, including the response to environmental perturbations, are negligible. Systematic uncertainties are comparable to statistical uncertainties in the characterization of the beam response. Both are accounted for in the covariance matrix of the window function and are propagated to uncertainties in the final power spectrum. We characterize the combined upper limits to residual systematic uncertainties through the pixel covariance matrix.

Hinshaw, G.↗
DOE OSTI2025

Selection function of clusters in Dark Energy Survey year 3 data from cross-matching with South Pole Telescope detections

Context. Galaxy clusters selected based on overdensities of galaxies in photometric surveys provide the largest cluster samples. However, modeling the selection function of such samples is complicated by noncluster members projected along the line of sight (projection effects) and the potential detection of unvirialized objects (contamination). Aims. We empirically constrained the magnitude of these effects by cross-matching galaxy clusters selected in the Dark Energy Survey data with the redMaPPer algorithm with significant detections in three South Pole Telescope surveys (SZ, pol-ECS, pol-500d). Methods. For matched clusters, we augmented the redMaPPer catalog with the SPT detection significance. For unmatched objects we used the SPT detection threshold as an upper limit on the SZe signature. Using a Bayesian population model applied to the collected multiwavelength data, we explored various physically motivated models to describe the relationship between observed richness and halo mass. Results. Our analysis reveals a clear preference for models with an additional skewed scatter component associated with projection effects over a purely log-normal scatter model. We rule out significant contamination by unvirialized objects at the high-richness end of the sample. While dedicated simulations offer a well-fitting calibration of projection effects, our findings suggest the presence of redshift-dependent trends that these simulations may not have captured. Our findings highlight that modeling the selection function of optically detected clusters remains a complicated challenge that requires a combination of simulation and data-driven approaches.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2007

Implementation of Leak Test Methods for the International Space Station (ISS) Elements, Systems and Components

The International Space Station (ISS has Qualification and Acceptance Environmental Test Requirements document, SSP 41172 that includes many environmental tests such as Thermal vacuum & Cycling, Depress/Repress, Sinusoidal, Random, and Acoustic Vibration, Pyro Shock, Acceleration, Humidity, Pressure, Electromatic Interference (EMI)/Electromagnetic Compatibility (EMCO), etc. This document also includes (13) leak test methods for Pressure Integrity Verification of the ISS Elements, Systems, and Components. These leak test methods are well known, however, the test procedure for specific leak test method shall be written and implemented paying attention to the important procedural steps/details that, if omitted or deviated, could impact the quality of the final product and affect the crew safety. Such procedural steps/details for different methods include, but not limited to: - Sequence of testing, f or example, pressurization and submersion steps for Method I (Immersion); - Stabilization of the mass spectrometer leak detector outputs fo r Method II (vacuum Chamber or Bell jar); - Proper data processing an d taking a conservative approach while making predictions for on-orbit leakage rate for Method III(Pressure Change); - Proper Calibration o f the mass spectrometer leak detector for all the tracer gas (mostly Helium) Methods such as Method V (Detector Probe), Method VI (Hood), Method VII (Tracer Probe), Method VIII(Accumulation); - Usage of visibl ility aides for Method I (Immersion), Method IV (Chemical Indicator), Method XII (Foam/Liquid Application), and Method XIII (Hydrostatic/Visual Inspection); While some methods could be used for the total leaka ge (either internal-to-external or external-to-internal) rate requirement verification (Vacuum Chamber, Pressure Decay, Hood, Accumulation), other methods shall be used only as a pass/fail test for individual joints (e.g., welds, fittings, and plugs) or for troubleshooting purposes (Chemical Indicator, Detector Probe, Tracer Probe, Local Vacuum Chamber, Foam/Liquid Application, and Hydrostatic/Visual Inspection). Any isolation of SSP 41172 requirements have led to either retesting of hardware or accepting a risk associated with the potential system or component pressure integrity problem during flight.

Underwood, Steve↗
NASA NTRS2023

The Impact and Estimation of Uncertainty Correlation for Multi-Angle Polarimetric Remote Sensing of Aerosols and Ocean Color

Multi-angle polarimetric (MAP) measurements contain rich information for characterization of aerosol microphysical and optical properties that can be used to improve atmospheric correction in ocean color remote sensing. Advanced retrieval algorithms have been developed to obtain multiple geophysical parameters in the atmosphere-ocean system, although uncertainty correlation among measurements is generally ignored due to lack of knowledge on its strength and characterization. In this work, we provide a practical framework to evaluate the impact of the angular uncertainty correlation from retrieval results and a method to estimate correlation strength from retrieval fitting residuals. The Fast Multi-Angular Polarimetric Ocean coLor (FastMAPOL) retrieval algorithm, based on neural network forward models, is used to conduct the retrievals and uncertainty quantification. In addition, we also discuss a flexible approach to include a correlated uncertainty model in the retrieval algorithm. The impact of angular correlation on retrieval uncertainties is discussed based on synthetic AirHARP and HARP2 measurements using a Monte Carlo uncertainty estimation method. Correlation properties are estimated using auto-correlation functions based on the fitting residuals from both synthetic AirHARP and HARP2 data and real AirHARP measurement, with the resulting angular correlation parameters found to be larger than 0.9 and 0.8 for reflectance and DoLP, respectively, which correspond to correlation angles of 10° and 5°. Although this study focuses on angular correlation from HARP instruments, the methodology to study and quantify uncertainty correlation is also applicable to other instruments with angular, spectral, or spatial correlations, and can help inform laboratory calibration and characterization of the instrument uncertainty structure.

PACE↗
NASA NTRS2024

Toward 2D Dynamo Models Calibrated by Global 3D Relativistic Accretion Disk Simulations

Two-dimensional models assuming axisymmetry are an economical way to explore the long-term evolution of black hole accretion disks, but they are only realistic if the feedback of the nonaxisymmetric turbulence on the mean momentum and magnetic fields is incorporated. Dynamo terms added to the 2D induction equation should be calibrated to 3D MHD simulations. For generality, the dynamo tensors should be calibrated as functions of local variables rather than explicit functions of spatial coordinates in a particular basis. In this paper, we study the feedback of non-axisymmetric features on the 2D mean fields using a global 3D, relativistic, Cartesian simulation from the Illinois GRMHD code. We introduce new methods for estimating overall dynamo alpha and turbulent diffusivity effects as well as measures of the dominance of non-axisymmetric components of energies and fluxes within the disk interior. We attempt closure models of the dynamo EMF using least squares fitting, considering both models where coefficient tensors are functions of space and more global, covariant models. None of these models are judged satisfactory, but we are able to draw conclusions on what sorts of generalizations are and are not promising.

Matthew D Duez↗
DOE OSTI2026

Impact of in situ nuclear networks and atomic opacities on neutron star merger ejecta dynamics, nucleosynthesis, and kilonovae

Context. Binary neutron star merger (BNSM) ejecta are key sites of rapid neutron capture (r-process) nucleosynthesis and they produce kilonovae powered by the radioactive decay of freshly synthesized nuclei. Modeling their evolution requires multi-physics simulations involving hydrodynamics, nuclear reactions, and radiative processes. The impact of nuclear burning and atomic opacity is poorly understood and often treated with simplified prescriptions. Aims. We systematically investigate different treatments of nuclear heating, particle thermalization, and atomic opacities in radiation-hydrodynamics simulations of BNSM ejecta and kilonova light curves. Methods. Ejecta profiles from long-term numerical-relativity simulations of asymmetric neutron star binaries with a massive neutron star remnant were evolved to ∼30 days using a 2D ray-by-ray approach. We compared simplified heating-rate and thermalization prescriptions with in situ Nuclear reaction Network (NN) calculations that track nuclear energy deposition and include a composition-dependent thermalization scheme. We also contrasted various gray opacity models with a frequency-dependent treatment based on atomic calculations. Results. Coupling NN and hydrodynamics significantly affects nucleosynthesis and kilonova emission. Assuming homologous expansion alters abundance evolution and produces a narrower, less populated second r-process peak and a third peak shifted to higher mass numbers. The back-reaction of nuclear heating affects the temperature evolution enough to delay and redden the early (t∼ hours) kilonova peaks. A constant thermalization efficiency underestimates and reddens the early emission while overestimating the late-time luminosity compared to the composition-dependent treatment. Analytical opacity prescriptions yield a more extended, colder photosphere, resulting in dimmer, redder kilonovae at early times (t≲ hour), while the delayed recession of the photosphere prolongs the red emission at t ≳ 5 days. Conclusions. Coupling hydrodynamics to an in situ NN is crucial for reliable nucleosynthesis and kilonova predictions. Resolving the first several hundred milliseconds of the hydrodynamics is essential for robust nucleosynthesis calculations. Composition-dependent thermalization and frequency-dependent, atomic-physics-based opacities are needed to accurately capture the temperature evolution of the ejecta and the brightness and color evolution of the kilonova. Calibrated analytic nuclear-power fits with simplified thermalization and opacity prescriptions can still reproduce the density and temperature evolution of the ejecta.

74 ATOMIC AND MOLECULAR PHYSICS↗