Engineering PapersSearch

SEARCH · Engineering Papers

Results for “uncertainty estimation”

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.

146 records · Page 8

Heliostat sizing methodology for concentrating solar thermal industrial process heat projects

This study presents a method to obtain a heliostat size that minimizes the levelized cost of heat (LCOH) of a heliostat-based concentrating solar thermal system for applications of solar heating for industrial processes at operating temperatures from 565 to 1550°C. The method extends prior work by embedding a routine for system design that obtains near-optimal subsystem sizes, increasing the fidelity of drive cost functions, and adding an optical performance model to supplement the previously developed cost models, which we update to reflect current pricing trends. An illustrative business case is developed for Daggett, California, targeting specified annual thermal energy outputs of 50 to 400 GWh th . Optical performance is modeled using verified estimates from the literature. A surrogate heliostat cost model, derived from commercial heliostat designs and scaled for production volume, installation, and operations and maintenance costs, is used to develop cost functions. Results show that heliostat size strongly affects the LCOH, producing a characteristic U-shaped trend with a robust near-optimal window of 7-20 m 2 ; the heliostat size producing the lowest project cost in our study grows slightly as the project size increases, and is reduced as the operating temperature increases. The findings in this study are consistent with the general trend of smaller heliostats being deployed at existing projects for high-temperature industrial process heat and reflect the significant reduction in power electronics and other per-heliostat costs. The methodology we propose is general and can be tailored to revised cost curves as the technology continues to evolve.

14 SOLAR ENERGY

Time-dependent scenario modeling for the ST-E1 fusion power plant

ST-E1 is a low aspect ratio fusion power plant being designed by Tokamak Energy targeting 1.5 GW of fusion power. Characterization of the ST-E1 flat-top scenario is described elsewhere McNamara et al (2026 Nucl. Fusion 66 086008); here we focus on addressing the question of how to ramp-up the ST-E1 plasma from an initial state following breakdown and flux-surface formation to the target flat-top state. Being low-aspect ratio, the available solenoid flux of ST-E1 is limited. Therefore, particular consideration is placed on developing ramp-up scenarios that predominantly use inductive flux provided by external vertical field coils. Through time-dependent modeling with METIS, we show that this is possible when the ramp-up is performed at relatively high plasma density: although auxiliary current drive efficiency is reduced, this is significantly outweighed by (1) higher electron-ion collisional equilibration, (2) higher fusion power once ions become sufficiently hot, (3) higher poloidal beta for increased vertical-field flux, and (4) potentially favorable exhaust compatibilities. Ultimately, we show the target ST-E1 flat-top performance can be reached after a ramp-up period lasting 150 s using less than 40 Vs of solenoid flux (with vertical field providing ∼ 90 Vs of flux). The sensitivity to model assumptions are presented, with the general observation that deleterious effects can be mitigated through minor alterations of the auxiliary power temporal waveform and/or total auxiliary power level. The impact of a solenoid and the auxiliary power mix (electron cyclotron heating only versus electron and ion cyclotron heating) on the ST-E1 ramp-up success are also discussed in appendices. On this latter topic, we show that the effect of direct-ion heating during ramp-up is obscured by the uncertainty in the pedestal dynamics, identifying a clear line of future work required to make a definite decision on the ST-E1 auxiliary power mix.

current ramp-up

Dark energy survey year 3 results: likelihood-free, simulation-based w CDM inference with neural compression of weak-lensing map statistics

We present simulation-based cosmological wcold dark matter (wCDM) inference using dark energy survey year 3 weak-lensing maps, via neural data compression of weak-lensing map summary statistics: power spectra, peak counts, and direct map-level compression/inference with convolutional neural networks (CNN). Using simulation-based inference, also known as likelihood-free or implicit inference, we use forward-modelled mock data to estimate posterior probability distributions of unknown parameters. This approach allows all statistical assumptions and uncertainties to be propagated through the forward-modelled mock data; these include sky masks, non-Gaussian shape noise, shape measurement bias, source galaxy clustering, photometric redshift uncertainty, intrinsic galaxy alignments, non-Gaussian density fields, neutrinos, and non-linear summary statistics. We include a series of tests to validate our inference results. This paper also describes the Gower Street simulation suite: 791 full-sky pkdgrav3 dark matter simulations, with cosmological model parameters sampled with a mixed active-learning strategy, from which we construct over 3000 mock dark energy survey lensing data sets. For wCDM inference, for which we allow –1 < w < –$\frac{1}{3}$⁠, our most constraining result uses power spectra combined with map-level (CNN) inference. Using gravitational lensing data only, this map-level combination gives Ω m = 0.283$^{+0.020}_{–0.027}$⁠, S 8 = 0.804$^{+0.025}_{–0.017⁠}$, and w < –0.80 (with a 68 per cent credible interval); compared to the power spectrum inference, this is more than a factor of two improvement in dark energy parameter (Ω⁠ DE , w⁠) precision.

79 ASTRONOMY AND ASTROPHYSICS

Observational evidence for Criegee intermediate oligomerization reactions relevant to aerosol formation in the troposphere

Criegee intermediates are reactive intermediates that are implicated in transforming the composition of Earth’s troposphere and in the formation of secondary organic aerosol, impacting Earth’s radiation balance, air quality and human health. Yet, direct identification of their signatures in the field remains elusive. Here, from particulate and gas-phase mass-spectrometric measurements in the Amazon rainforest, we identify sequences of masses consistent with the expected signatures of oligomerization of the CH 2 OO Criegee intermediate, a process implicated in ozonolysis-driven aerosol formation. In this study, we assess the potential contributions of oligomerization through laboratory ozonolysis experiments, direct kinetic studies of Criegee intermediate reactions, and high-level theoretical calculations. Global atmospheric models built on these kinetics results indicate that Criegee intermediate chemistry may play a larger role in altering the composition of Earth’s troposphere than is captured in current atmospheric models, especially in areas of high humidity. However, the models still capture only a relatively small fraction of the observed signatures, suggesting considerable underestimates of Criegee intermediate concentrations and reactivity and/or the dominance of other, presently uncharacterized, oxidation mechanisms. Resolving the remaining uncertainties in emission inventories and the effects of atmospheric water vapour on key chemical reactions will be required to definitively assess the role of Criegee intermediate oligomerization reactions.

58 GEOSCIENCES

Weibel-mediated filamentary structures observed in the ICF context

Here, in light of novel and past experimental results, we demonstrate how Weibel-mediated filamentary structures can develop in the expanding plasma plume of a laser-irradiated foil. The transverse ballistic cooling that occurs during the quasi-spherical plasma expansion naturally drives an electron pressure anisotropy, resulting in the growth of electron current filaments. This effect competes with electron–ion Coulomb collisions, which tend to isotropize the electron distribution function. Based on theoretical and particle-in-cell modeling, we provide estimates of the dominant wavelength and amplitude of the self-generated magnetic fluctuations, which are found to explain experimental data obtained at the OMEGA and Laser Megajoule facilities.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY

Assessing the Impact of a Forest Canopy on Near-Surface Wind Statistics

Representing the forest canopy in atmospheric numerical models should improve simulated winds within and above the canopy up to a few hundred meters above the ground. Here, in this study, we implement a forest canopy parameterization into the Weather Research and Forecasting (WRF) Model in a large-eddy simulation (LES) mode by applying drag forces across multiple layers within the canopy height. We use unique observations from the Lidar Experiments for Assessing Flow over Forests (LEAFF) field campaign at the Wind River Experimental Forest (WREF) in the U.S. Pacific Northwest to evaluate model performance. In a 2-day case study, the canopy parameterization improved wind predictions both within and above the canopy, particularly during the daytime and at finer grid resolution. Without it, winds were frequently overpredicted above the canopy. Similarly, derived quantities such as the wind shear index also yielded estimates closer to observations with the canopy parameterization implemented. These findings suggest that representing the canopy using drag forces alone can improve simulated mean winds up to 200 m above the surface. Furthermore, second-order statistical moments of wind were more sensitive to canopy density than first-order moments, especially during the daytime. This increased sensitivity and the improved daytime performance in wind speed—evidenced by the lowest bias from observations (3% compared to 20% over diurnal cycle)—imply that winds above the canopy layer are strongly influenced by how well turbulence above the canopy is modeled. The results of this study can serve as a foundation for parameterizing forest canopy effects in coarser weather forecast models.

Energy - Wind

Deep Learning and Photogrammetric Reconstruction for Automated Crack Detection and Dimensional Measurement in Mining Operations

Surface crack detection and dimensional measurement at active mining sites present significant safety and operational challenges. Manual inspection methods are labor-intensive, spatially incomplete, and expose personnel to hazardous environments, while existing automated approaches have been developed primarily for concrete civil infrastructure and have not been validated on the complex, variable surfaces characteristic of mining environments. This dissertation presents an automated pipeline that integrates deep learning semantic segmentation with Structure-from-Motion photogrammetry to detect surface cracks and measure their aperture, length, and vertical displacement from standard RGB imagery acquired during routine Uncrewed Aerial Vehicle (UAV) survey operations, without requiring additional sensor hardware or manual measurement. The pipeline combines a U-Net architecture with an EfficientNet-B0 encoder, pretrained on the SDNET2018 concrete crack dataset and fine-tuned on a mining-specific dataset spanning laboratory concrete specimens, coal refuse impoundment embankments, and post-blast limestone quarry benches. Photogrammetric reconstruction is performed using COLMAP Structure-from-Motion and Multi-View Stereo, with crack segmentation masks projected into the reconstructed point cloud to enable three-dimensional vertical displacement measurement through local plane fitting and bimodal surface detection. The pipeline was validated across 36 controlled laboratory specimens at three imaging distances and four vertical displacement levels, achieving aperture measurement RMSE of 0.047 cm and R² of 0.954, and vertical displacement RMSE of 0.140 cm and R² of 0.966, against independent caliper measurements. Field application at a coal refuse impoundment in southwestern Pennsylvania detected 71 crack components across the embankment crest, with a dominant longitudinal crack exhibiting aperture values reaching 28 cm and a 95th percentile vertical displacement of 35.53 cm, consistent in magnitude and spatial distribution with simultaneously acquired LiDAR-derived estimates. Application across four post-blast limestone quarry bench datasets in California successfully characterized blast-induced fracture networks at ground sampling distances ranging from 0.59 to 1.23 cm/pixel, with detected crack geometries physically consistent with observable surface conditions at each site. The results demonstrate that deep learning-based crack detection and photogrammetric measurement can be integrated into routine UAV inspection workflows at mining sites, providing repeatable, scalable, and quantitative crack characterization across surface types, crack scales, and displacement magnitudes not previously addressed in the literature. The pipeline requires no dedicated surveying equipment beyond the UAV platforms already deployed at mine sites for survey and monitoring purposes, supporting practical adoption within existing operational workflows.

Crack detection, Dimensional Measurement

Linking Community‐Climate Disequilibrium to Ecosystem Function

Turnover in species composition often lags behind the pace of climate change, resulting in mismatches between climate and communities. However, the impact of these community‐climate disequilibria on ecosystem functions is rarely considered, and current methods for measuring disequilibria assume that species ranges were, until recently, in equilibrium with climate. Here, in this work, we develop a simple theoretical model to address both of these problems by linking community‐climate disequilibrium with ecosystem functioning. We show how disequilibrium can impair functioning in the near‐term even when climate change is expected to enhance functioning in the long‐term. Responses are most likely to change over time in communities where turnover is slow, the impact of disequilibrium counteracts the direct effects of climate on ecosystem function, and pre‐existing disequilibrium is large. These findings emphasise the importance of precise and unbiased estimates of community‐climate disequilibria for improving ecological forecasts. By fitting our model to time series of both climate and ecosystem function from a metacommunity simulation, we show the potential for community‐climate disequilibrium to be inferred without direct knowledge about species' distributions or climatic tolerances. We end by outlining a research agenda to apply dynamic disequilibrium concepts and test novel hypotheses across diverse ecosystems.

climate change

Shakeup and shakeoff spectra in the electron capture decay of atomic 7 Be

The most stringent laboratory-based experimental limits on the existence of submegaelectronvolt sterile neutrinos are currently set by decay spectroscopy of radioactive 7 Be embedded into superconducting sensors. The systematic uncertainties are dominated by the modeling of the electron shakeup and shakeoff spectra that are not based on state-of-the-art atomic theory and do not include electron correlations or relativistic effects. We have used the multiconfiguration Dirac-Fock formalism to obtain correlated wave functions ab initio and compute all single and double shake processes in the electron capture decay of atomic 7 Be . The simulations can explain some but not all of the observed spectral features, likely because the wave functions are modified by the Ta sensor material that the 7 Be is embedded into. The new models also show that the 𝐿/𝐾 electron capture ratio of 7 Be in Ta has previously been slightly underestimated revising the previous value of 0.070(7) to a new value of 0.0756(20).

73 NUCLEAR PHYSICS AND RADIATION PHYSICS

Thermal Disorder‐Induced Strain and Carrier Localization Activate Reverse Halide Segregation

The reversal of halide ions is studied under various conditions. However, the underlying mechanism of heat-induced reversal remains unclear. This work finds that dynamic disorder-induced localization of self-trapped polarons and thermal disorder-induced strain (TDIS) can be co-acting drivers of reverse segregation. Localization of polarons results in an order of magnitude decrease in excess carrier density (polaron population), causing a reduced impact of the light-induced strain (LIS – responsible for segregation) on the perovskite framework. Meanwhile, exposing the lattice to TDIS exceeding the LIS can eliminate the photoexcitation-induced strain gradient, as thermal fluctuations of the lattice can mask the LIS strain. Under continuous 0.1 W cm -2 illumination (upon segregation), the strain disorder is estimated to be 0.14%, while at 80 °C under dark conditions, the strain is 0.23%. However, in situ heating of the segregated film to 80 °C under continuous illumination (upon reversal) increases the total strain disorder to 0.25%, where TDIS is likely to have a dominant contribution. Therefore, the contribution of entropy to the system's free energy is likely to dominate, respectively. Various temperature-dependent in situ measurements and simulations further support the results. These findings highlight the importance of strain homogenization for designing stable perovskites under real-world operating conditions.

36 MATERIALS SCIENCE

Hazard and risk analysis framework for nuclear power plant–based integrated energy systems

Employing integrated energy systems (IESs) with nuclear power plants (NPPs) can improve NPP utilization by leveraging dedicated thermal and electric power delivery, but it may also increase operational safety risks. This paper presents a framework to identify and quantify hazards and risks for such IESs. The framework combines accidentology to review past industrial accidents with failure modes and effects analysis (FMEA) to identify potential future incidents. Hydrogen explosion and toxic chemical release hazards are of particular concern. Explosion consequences are quantified using the Bauwens-Dorofeev (Bauwens) and trinitrotoluene equivalent mass (TNT-EM) methods, while chemical release consequences are computed using the Gaussian atmospheric dispersion method. Operational disturbances from direct electrical and thermal integration that may affect NPP safety are modeled using probabilistic risk analysis (PRA). Hazards and risks are then evaluated for regulatory compliance. The framework is applied to IESs comprising pressurized or boiling water reactors supplying three levels of thermal and electrical power to industrial customers. Case studies include high-temperature steam electrolysis hydrogen plants of varying capacities and a synthetic fuel production plant. Sensitivity analysis examines piping component failures in the PRA model as a precursor to cost estimation for thermal extraction line design. Additionally, Fussel-Vessely (FV) and risk increase importance (RII) measures identify risk-informed design improvements for the thermal extraction system. FMEA highlights hazards such as loss of offsite power, prompt loss of electrical load, loss of thermal output, and immediate steam diversion, in addition to hydrogen explosions and toxic chemical releases. Both Bauwens and TNT-EM methods suggest maintaining several hundred meters of separation between the NPP and hydrogen facility to mitigate explosion risks. PRA results show a maximum initiating event frequency increase of 1.15% and an overall risk increase of 0.28%. Importance measure analysis identifies upstream pipe leak isolation components as critical. Evaluating the results against safety regulations, it is concluded that hazards and risks can be managed to comply with regulations through risk-informed thermal and electrical connection designs, component selection, maintenance programs, and safe separation distances between NPPs and integrated industrial facilities.

08 - HYDROGEN

Heterotrophic respiration by soil microbes in a changing climate

Soil microbes strongly influence the soil organic carbon (SOC) pool, which globally stores ~2,000 PgC. Specifically, the balance between microbial heterotrophic respiration (R H ), which degrades SOC, and plant–microbe interactions that stabilize SOC determines whether terrestrial ecosystems are a net source or sink of CO 2 to the atmosphere. Here, in this Review, we evaluate how climate change alters these competing processes. R H is approximately half of total soil respiration, at ~50 PgC yr −1 , with 70% occurring in topsoils. Warming accelerates microbial metabolism, with a 10 °C temperature increase estimated to raise R H by ~50%, an effect that is particularly strong in Arctic soils. Warming also reduces soil moisture, further modulating R H , which responds nonlinearly to soil moisture, being limited by saturation and desiccation and meeting a maximum at intermediate levels. Consequently, R H is highly sensitive to future precipitation changes and drought. However, soil management strategies could enhance SOC stocks and persistence under climate change. Bacterial and fungal inoculants can promote SOC production and stabilization, while deep-rooting plants increase SOC inputs to deeper layers that experience lower R H . Agricultural practices and biochar amendments can also enhance SOC and reduce R H . Expanding field trials across regions, climates and soil types would improve empirical understanding of these responses and support better representation of R H in predictive models, enabling more accurate assessments of climate impacts on SOC storage.

Jansson, Janet K. [Pacific Northwest National Labo

Copacabana: a probabilistic membership assignment method for galaxy clusters

Cosmological analyses using galaxy clusters in optical/near-infrared photometric surveys require robust characterization of their galaxy content. Precisely determining which galaxies belong to a cluster is crucial. In this paper, we present the COlor Probabilistic Assignment of Clusters And BAyesiaN Analysis (Copacabana) algorithm. Copacabana computes membership probabilities for all galaxies within an aperture centred on the cluster using photometric redshifts, colours, and projected radial probability density functions. We use simulations to validate Copacabana and we show that it achieves up to 89 per cent membership accuracy with a mild dependence on photometric redshift uncertainties and choice of aperture size. We find that the precision of the photometric redshifts has the largest impact on the determination of the membership probabilities followed by the choice of the cluster aperture size. We also quantify how much these uncertainties in the membership probabilities affect the stellar mass–cluster mass scaling relation, a relation that directly impacts cosmology. Using the sum of the stellar masses weighted by membership probabilities (⁠μ * ⁠) as the observable, we find that Copacabana can reach an accuracy of 0.06 dex in the measurement of the scaling relation at low redshift for a Legacy Survey of Space and Time type survey. These results indicate the potential of Copacabana and μ * to be used in cosmological analyses of optically selected clusters in the future.

79 ASTRONOMY AND ASTROPHYSICS

First multi-institutional systematic comparison of the neutron ambient dose equivalent produced by proton therapy systems

Objective. Isochronous cyclotrons, synchrocyclotrons, and synchrotrons are used to accelerate protons for proton therapy. An accurate measurement of neutron doses generated by these accelerators and associated delivery systems and its clinical relevance requires systematic protocols and proper neutron dosimetry for a meaningful assessment. We present the first comprehensive comparison of neutron ambient dose equivalent (H*(10)) produced by clinically operational proton therapy systems. Approach. Treatment plans with 10 cm modulation-depth and ranges of 10 cm (R10M10) and 25 cm (R25M10) were created to cover a 10 × 10 × 10 cm 3 water target. The pencil beam scanning proton therapy machines studied were: two gantry-mounted synchrocyclotrons (Hyperscan, Mevion, half-gantry), two isochronous cyclotrons (ProBeam, Varian, full-gantry), one isochronous cyclotron (Proteus, IBA, full-gantry), and two synchrotrons (PROBEAT, Hitachi, full- and half-gantry). Proton beams were delivered to 30 × 30 × 40 cm 3 plastic water phantoms. WENDI-II and LUPIN-BF3-NP neutron rem-meters were positioned at three angles (0°, 45°, 90°) relative to the beam direction to measure the neutron H*(10) at distances between 50–300 cm from the isocenter. Main results. H*(10) showed dependence on beam energy, machine type, and measurement location. The highest reading was for the gantry-mounted synchrocyclotron, whereas other systems produced approximately comparable neutron doses. In all cases, the H*(10) reduced with distance from the isocenter. The H*(10) drop at 2 m distance compared to that at 0.5 m was a factor of ∼5 for the gantry-mounted synchrocyclotron whereas in other systems the decrease was a factor of 10. The WENDI-II device suffered from dead-time-associated under-estimation of the dose by a factor of ∼2–3 under the synchrocyclotron beam due to its high dose-per-pulse. However, WENDI-II and LUPIN-BF3-NP results were within reasonable agreement in isochronous cyclotron and synchrotron beams, indicating that both devices are suitable for those systems. Significance. Neutron H*(10) is dependent on various parameters including beam energy, measurement location, as well as machine design. Caution must be exercised in choosing the appropriate neutron-dose-measurement device to be used for low-duty-factor, particularly in high-instantaneous-rate proton delivery systems. By delivering the same volumetric proton dose across different machines, this work provides a benchmark for inter-system comparisons and serves as a foundation for future studies.

LUPIN

Experimental observation and integrated modelling of proton-beryllium fusion in He and D plasmas at JET

Validated integrated modelling of JET ITER-like wall experiments in which fusion performance is driven by reactions between fast ions and intrinsically present metal wall impurities is presented. A steady-state L-mode plasma with dominant proton-beryllium fusion and neutron yields of up to ≈ 6·10 13 s -1 is developed in He and D, via radiofrequency heating of a H minority. The fusion drive is unambiguously confirmed by the neutral particle analyser, fast ion loss detector, and γ-ray diagnostics. Experiments are analysed via an integrated modelling framework, developed to model the two-stage proton beryllium-fusion chain and produce high-fidelity fusion product source terms. The modelling chain comprises TRANSP and JETTO for plasma core modelling, LOCUST for full orbit product tracking and collisional slowing-down, DRESS to resolve two- and three-body fusion kinematics, and MCNP for neutron transport calculations. Modelling shows that the primary 9 Be(p,n) 9 B reaction is the dominant neutron emitter at naturally present concentrations of beryllium in these experiments. The yield contribution of secondary reactions between fusion products and beryllium, 9 Be(d,n) 10 B and 9 Be(α,n) 12 C, is found to be negligible. The proton-deuteron knock-on effect in D plasmas is modelled, which is calculated to contribute ≈ 25% to the total neutron yield. For both He and D discharges the total computed neutron rates match fission chamber (FC) measurements within the combined experimental and computational uncertainty, with an average discrepancy of ≈ ± 20%. Realistic proton-beryllium neutron sources are propagated through JET’s MCNP neutron transport model which shows that 235 U FCs’ response is sensitive to p–Be source changes, with up to ≈ 10% variation compared to a D–D neutron source. We show that the high-energy tail of the fast proton minority can be studied with multi-foil neutron activation. The framework is also applied to the study of interactions between fast protons and boron impurities, of relevance to ITER. We calculate that in JET conditions a significant alpha source with DT-like energies could be generated through 11 B(p, α)2α fusion, and detected via γ-emission in secondary interactions between fast alphas and boron. The work represents an important step towards validating predictive integrated modelling capabilities for non-standard fusion reactions.

JET

Single-Nucleon Transfer on Unstable 59 Cu Probes the NiCu Cycle in Astrophysical X-Ray Bursts

Recent models of the rapid proton (𝑟⁢𝑝) capture process indicate that a competition between the 59 Cu⁢(𝑝,𝛾) ⁢60 Zn and 59 Cu⁢(𝑝,𝛼) ⁢56 Ni reactions may result in the formation of a nickel-copper (NiCu) cycle that traps the flux of material between 56 Ni and 60 Zn . Here, in this work, we report the identification of 15 proton-unbound levels in 60 Zn , populated via 59 Cu⁢(𝑑,𝑛) transfer, which govern the rate of the 59 Cu⁢(𝑝,𝛾)⁢ 60 Zn reaction in XRBs. Precise excitation energies for levels in 60 Zn were obtained from observed 𝛾 decays, and spectroscopic factors were determined from angle-integrated cross sections. Incorporating these results into stellar-model calculations, we find that with experimentally constrained uncertainties a NiCu cycle in XRBs is indeed possible, though we limit its branching strength to less than 38%. While modest, such a branching has significant impact on the light curve, motivating further studies of the relevant rates. Our calculations also indicate that a significant NiCu cycle leads to an increase in the amount of odd-𝐴 nuclei in the burst ashes, which may affect Urca cooling processes in neutron star crusts.

direct reactions

Study of Key 57 Ni ⁢(𝑝,𝛾)⁢ 58 Cu Resonances to Understand 44 Ti Nucleosynthesis in Supernovae

An important validation of nucleosynthesis models of core-collapse supernovae is the comparison of radioisotope predictions to abundances inferred from observations of 𝛾 rays emitted in remnants. One such isotope, 44 Ti, is especially sensitive to the 57 Ni ⁢(𝑝,𝛾)⁢ 58 Cu reaction rate. Despite this importance, no experimentally constrained rate exists for this reaction. It is thus crucial to determine this rate. GRETINA ORRUBA: Dual Detectors for Experimental Structure Studies (GODDESS) was used along with the first use of the Enge split-pole spectrograph at Notre Dame to measure structure properties of 58 Cu via the 58 Ni ⁢( 3 He,𝑡)⁢ 58 Cu reaction. Here, the combined analysis of these complementary data sets allows for precise determination of 58 Cu level energies, and spins are constrained from analysis of the 𝛾-ray decay. 44 Ti yields were found to change by over 25% compared to previous estimates in model calculations using the new experimentally constrained reaction rate.

energy levels

Metrology for femtosecond pulsed x-ray heating in diamond anvil cell experiments at the European XFEL: Revisiting the iron phase diagram up to 150 GPa

The development of pulsed intense x-ray sources, such as free electron laser, offers new avenues for high pressure experiments. Here, we study the feasibility and metrology of x-ray heating in diamond anvil cells at the European x-ray free electron laser. This method enables one to volumetrically heat the sample while inhibiting chemical migration and probing the crystallographic structure of the sample throughout the heating with a high repetition rate. We focus our study on iron, whose phase diagram is well established up to 100 GPa, to explore the possibilities and limitations of this technique. We volumetrically heat iron samples at starting pressures ranging from 10 to 138 GPa, using the x-ray beam pulsed at 4.5 MHz in a serial pump-and-probe experimental design. Experimental challenges arise from temperature gradients within the sample, changes in temperature at the 100 ns timescale, the difficulty of direct temperature estimates, the effect of thermal pressure, and the presence of metastable crystallites due to rapid cycles of heating and cooling. Hence, we develop a multi-crystal-like data processing method that allows us to account for sample heterogeneity in probed conditions. We then calibrate our measurements using known physical properties of iron under pressure. Thermal pressure in our experiments increases from 4% of the isochoric prediction at 10 GPa to 23% at 138 GPa, and we show that our data are in agreement with most previous observations of iron in this pressure range. The method can now be implemented at higher pressures and temperatures and on materials with unknown phase diagrams.

Materials science