Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “reproducible science”

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 793 records · Page 44

Measurement of charge and light yields for 127 Xe L -shell electron captures in liquid xenon

Dark matter searches using dual-phase xenon time-projection chambers (LXe-TPCs) rely on their ability to reject background electron recoils (ERs) while searching for signal-like nuclear recoils (NRs). ER response is typically calibrated using β -decay sources, such as tritium, but these calibrations do not characterize events accompanied by an atomic vacancy, as in solar neutrino scatters off inner-shell electrons. Such events lead to emission of x rays and Auger electrons, resulting in higher electron-ion recombination and thus a more NR-like response than inferred from β -decay calibration. We present a cross-calibration of tritium β -decays and Xe 127 electron-capture decays (which produce inner-shell vacancies) in a small-scale LXe-TPC and give the most precise measurements to date of light and charge yields for the Xe 127 L -shell electron-capture in liquid xenon. We observe a 6.9 σ ( 9.2 σ ) discrepancy in the L -shell capture response relative to tritium β decays, measured at a drift field of 363 ± 14 V / cm ( 258 ± 13 V / cm ), when compared to simulations tuned to reproduce the correct β -decay response. In dark matter searches, use of a background model that neglects this effect leads to overcoverage (higher limits) for background-only multi-kiloton-year exposures, but at a level much less than the 1 - σ experiment-to-experiment variation of the 90% C.L. upper limit on the interaction rate of a 50 GeV / c 2 dark matter particle.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Remote heteroepitaxy of GaN microrod heterostructures for deformable light-emitting diodes and wafer recycle

There have been rapidly increasing demands for flexible lighting apparatus, and micrometer-scale light-emitting diodes (LEDs) are regarded as one of the promising lighting sources for deformable device applications. Herein, we demonstrate a method of creating a deformable LED, based on remote heteroepitaxy of GaN microrod (MR) p - n junction arrays on c -Al 2 O 3 wafer across graphene. The use of graphene allows the transfer of MR LED arrays onto a copper plate, and spatially separate MR arrays offer ideal device geometry suitable for deformable LED in various shapes without serious device performance degradation. Moreover, remote heteroepitaxy also allows the wafer to be reused, allowing reproducible production of MR LEDs using a single substrate without noticeable device degradation. The remote heteroepitaxial relation is determined by high-resolution scanning transmission electron microscopy, and the density functional theory simulations clarify how the remote heteroepitaxy is made possible through graphene.

36 MATERIALS SCIENCE↗

System size and energy dependence of the mean transverse momentum fluctuations at the LHC

Event-by-event fluctuations of the event-wise mean transverse momentum, $\langle$p T $\rangle$, of charged particles produced in proton–proton (pp) collisions at $\sqrt{s}$ = 5.02 TeV, Xe–Xe collisions at $\sqrt{s_{NN}}$ = 5.44 TeV, and Pb–Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV are studied using the ALICE detector based on the integral correlator $\langle$$\langle$Δp T Δp T $\rangle$$\rangle$. The correlator strength is found to decrease monotonically with increasing produced charged-particle multiplicity measured at midrapidity in all three systems. In Xe–Xe and Pb–Pb collisions, the multiplicity dependence of the correlator deviates significantly from a simple power-law scaling as well as from the predictions of the HIJING and AMPT models. The observed deviation from power-law scaling is expected from transverse radial flow in semicentral to central Xe–Xe and Pb–Pb collisions. In pp collisions, the correlation strength is also studied by classifying the events based on the transverse spherocity, S 0 , of the particle production at midrapidity, used as a proxy for the presence of a pronounced back-to-back jet topology. Low-spherocity (jetty) events feature a larger correlation strength than those with high spherocity (isotropic). The strength and multiplicity dependence of jetty and isotropic events are well reproduced by calculations with the PYTHIA 8 and EPOS LHC models.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Data-model files associated with the manuscript titled "The Importance of Explicitly Representing the Streambed in Watershed Models" (Shuai et al., 2023 HP)

This data package contains the model inputs and outputs used in the manuscript titled "The Importance of Explicitly Representing the Streambed in Watershed Models" (Shuai et al., 2023 HP). The data.zip file contains the data used to drive the Advanced Terrestrial Simulator (ATS) model simulations. The model.zip file contains the XML input file for ATS. The notebook.zip file contains the Jupyter notebooks for pre- and post- processing model results. The figures.zip file contains the raw figures associated with the manuscript. We generated this data package in support of the manuscript and research reproducibility. The background of this study is that the streambed itself has not been explicitly represented in watershed models, although the streambed characteristics are significantly different from those of its surrounding soil. We aim to answer the following questions: 1) How do streambed properties including hydraulic conductivity, thickness and resolution impact the groundwater-surface water exchange fluxes across the streambed? 2) Is an explicit representation of streambed important in watershed modeling?; 3) Is a high-resolution streambed needed for watershed simulations?

54 ENVIRONMENTAL SCIENCES↗

Optimizing cost-effective and benchmarked industry standards to quantify nutrient bioextraction by seaweed

Interest in the utility of seaweed farms to mitigate coastal eutrophication, or nutrient loading, has grown commensurate with the recent rise of the farmed seaweed industry in the U.S. But economic valuation of this ecosystem service remains elusive in part because of challenges in quantifying this spatiotemporally variable biological process with reproducible and comparable metrics. Regulatory bodies that permit wastewater discharge or lease area for aquaculture farms require water quality testing and reporting of dissolved total nitrogen (N) in nearshore marine environments. These metrics must meet EPA standards for testing and reporting (e.g., Total Kjeldahl Nitrogen - TKN). However, these metrics are inherently highly variable over space and time in dynamic nearshore systems, and expensive to evaluate with sufficient breadth to constrain this variance, creating a critical bottleneck to direct quantification of farmed seaweed net uptake rates in situ.

59 BASIC BIOLOGICAL SCIENCES↗

An ensemble of 48 physically perturbed model estimates of the 1/8° terrestrial water budget over the conterminous United States, 1980–2015

Terrestrial water budget (TWB) data over large domains are of high interest for various hydrological applications. Spatiotemporally continuous and physically consistent estimations of TWB rely on land surface models (LSMs). As an augmentation of the operational North American Land Data Assimilation System Phase 2 (NLDAS-2) four-LSM ensemble, this paper describes a dataset simulated from an ensemble of 48 physics configurations of the Noah LSM with multi-physics options (Noah-MP). The 48 Noah-MP physics configurations are selected to give a representative cross-section of commonly used LSMs for parameterizing runoff, atmospheric surface layer turbulence, soil moisture limitation on photosynthesis, and stomatal conductance. The dataset spans from 1980 to 2015 over the conterminous United States (CONUS) at a monthly temporal resolution and a 1/8° spatial resolution. The dataset variables include total evapotranspiration and its constituents (canopy evaporation, soil evaporation, and transpiration), runoff (the surface and subsurface components), as well as terrestrial water storage (snow water equivalent, four-layer soil water content from the surface down to 2 m, and the groundwater storage anomaly). The dataset is available at https://doi.org/10.5281/zenodo.7109816. Evaluations carried out in this study and previous investigations show that the ensemble per forms well in reproducing the observed terrestrial water storage, snow water equivalent, soil moisture, and runoff. Noah-MP complements the NLDAS models well, and adding Noah-MP consistently improves the NLDAS es timations of the above variables in most areas of CONUS. Besides, the perturbed-physics ensemble facilitates the identification of model deficiencies. The parameterizations of shallow snow, spatially varying groundwater dynamics, and near-surface atmospheric turbulence should be improved in future model versions.

54 ENVIRONMENTAL SCIENCES↗

HIRM v1.0: a hybrid impulse response model for climate modeling and uncertainty analyses

Simple climate models (SCMs) are frequently used in research and decision-making communities because of their flexibility, tractability, and low computational cost. SCMs can be idealized, flexibly representing major climate dynamics as impulse response functions, or process-based, using explicit equations to model possibly nonlinear climate and Earth system dynamics. Each of these approaches has strengths and limitations. Here we present and test a hybrid impulse response modeling framework (HIRM) that combines the strengths of process-based SCMs in an idealized impulse response model, with HIRM's input derived from the output of a process-based model. This structure enables the model to capture some of the major nonlinear dynamics that occur in complex climate models as greenhouse gas emissions transform to atmospheric concentration to radiative forcing to climate change. As a test, the HIRM framework was configured to emulate the total temperature of the simple climate model Hector 2.0 under the four Representative Concentration Pathways and the temperature response of an abrupt 4 times CO 2 concentration step. HIRM was able to reproduce near-term and long-term Hector global temperature with a high degree of fidelity. Additionally, we conducted two case studies to demonstrate potential applications for this hybrid model: examining the effect of aerosol forcing uncertainty on global temperature and incorporating more process-based representations of black carbon into a SCM. The open-source HIRM framework has a range of applications including complex climate model emulation, uncertainty analyses of radiative forcing, attribution studies, and climate model development.

54 ENVIRONMENTAL SCIENCES↗

Diabatization with Electrostatic Embedding for Studying Photophysics in Organic Molecular Crystals

Highly emissive organic molecular crystals find applications in several areas, such as organic electronics, solar cells, and sensors. Understanding the excited-state mechanisms underlying these applications is essential for optimizing and controlling them effectively. Exciton models coupled with nonadiabatic dynamics, particularly quantum dynamics, provide crucial insights into photochemical and photophysical processes in molecular crystals. Nevertheless, there remains a lack of general tools and automated workflows to facilitate such simulations. In this paper, we present a computational strategy to investigate the photoactivated dynamics of organic molecular crystals, bridging methodologies traditionally used for molecular systems and materials science, with a particular focus on the interplay between local excitations and charge transfer (CT) processes. We have implemented an interface between the fromage and Overdia programs, enabling the construction of vibronic Hamiltonians for molecular crystals within an excited-state ONIOM(QM:QM′) framework, incorporating long-range electrostatics through a RESP-based Ewald summation. Fragment-based diabatization provides a route to quantum dynamics simulations in weak-to-intermediate coupling regimes. The method was applied to the photophysics of dibenzo[g,p]chrysene (DBC) crystals using time-dependent DFT. The fromage/ Overdia interface was employed to compute the couplings of local excitations and CT states for 18 unique DBC dimers in the crystal and to quantify the influence of electrostatic embedding, which was found to be modest (10−20%). Simulations on π-stacked dimers reproduced the small red shift observed experimentally from solution to crystal, attributed to electronic interactions among fixed monomers rather than crystal electrostatics. Quantum dynamics simulations revealed ultrafast population transfer from bright local excitations to CT states. This approach establishes a robust framework linking molecular and solid-state excited-state dynamics, with potential applications for studying excitations, defects, and impurities in molecular crystals.

Crystals↗

Multi-objective surrogate-assisted calibration of CPFEM models using macroscopic response and in situ EBSD measurements of grain reorientation trajectories

Crystal plasticity finite element method (CPFEM) models are widely used to simulate the deformation behaviour of polycrystalline materials, but their calibration is often limited by their high computational cost and the non-convexity of the optimisation landscape. Here, this study develops a multi-objective surrogate-assisted calibration workflow that couples a multi-objective genetic algorithm (MOGA) with an adaptively trained deep neural network (DNN) surrogate model to efficiently identify CPFEM parameters from experimental data. The workflow is demonstrated on three crystal plasticity (CP) formulations of increasing complexity — Voce hardening (VH), two-coefficient latent hardening (LH2), and six-coefficient latent hardening (LH6) — using in situ electron backscatter diffraction (EBSD) measurements of Alloy 617 under uniaxial tensile loading. The CPFEM models are calibrated against the experimentally observed stress–strain response and reorientation trajectories of eight grains, then validated against eight additional trajectories and overall texture evolution. Across the CP formulations, the macroscopic response was reproduced reliably, while differences emerged in the robustness and accuracy of the grain-scale predictions. Including grain reorientation trajectories in the multi-objective calibration improved texture evolution predictions and filtered out physically inconsistent parameter sets that can arise from calibrating against only the stress–strain data. The workflow also demonstrates good transferability of calibrated parameters from a low- to a high-fidelity microstructural model. These results provide practical guidance for integrating in situ microstructural data into CPFEM through efficient, repeatable, and physically meaningful multi-objective calibration.

Crystal plasticity finite element method↗

Nanoparticle-enhanced absorptivity of copper during laser powder bed fusion

We report laser powder bed fusion (LPBF) of pure copper for thermal and electrical applications is hampered by its low near-infrared absorptivity and high thermal diffusivity. These material properties make it very difficult to localize the thermal energy needed to produce high density 3D printed parts. Modification of metal powders via nanoparticle additives is a promising approach to increasing absorptivity, but the effect of nanoparticles on absorptivity and melting behavior during LPBF is not well understood. In this study, we developed an in situ calorimetry system to measure effective absorptivity during LPBF on copper substrates. We decorated copper substrates using three nanoparticle systems (CuS, TiB 2 , multilayer graphene flakes) and demonstrated an enhanced absorptivity of the decorated substrates relative to pure copper. Graphene nanoflakes resulted in the highest improved absorption relative to pure copper from 0.09 to 0.48, due to their stability at high laser scanning powers. A thermomechanical model with convective heat transfer provided confidence in the measurements by reproducing the experimental melt pool traces. Full 3D cylindrical prints demonstrated an improvement in relative density of the copper-graphene powder prints (in the range of 0.930–0.992), relative to that of as-purchased copper powder prints (in the range of 0.854–0.972). This work provides a fundamental study of nanoparticle-enabled LPBF of highly reflective metals and demonstrates a viable route for expanding the library of reliably printable metals.

36 MATERIALS SCIENCE↗

First-principles investigation of uranium mononitride (UN): Effect of magnetic ordering, spin-orbit interactions and exchange correlation functional

Uranium mononitride (UN) is a promising nuclear fuel that combines the advantageous properties of readily used UO 2 and uranium alloys. Various properties of UN have been previously studied using different density functional theory (DFT) methodologies; however, there are still inconsistencies when it comes to the dynamical stability and defect properties of UN. We address these inconsistencies by studying the UN phonons and defect properties using DFT calculations employing two generalized gradient approximation (GGA) exchange-correlation functionals: PBE and AM05, with and without an added on-site Coulomb repulsion term (+U). Furthermore, we investigate the importance of spin-orbit coupling (SOC) when calculating the properties of UN. We use the different methodologies to determine the preference of UN to have antiferromagnetic (AFM) ordering, as seen in experiments, or ferromagnetic (FM) ordering of the uranium spins. We compare the crystallographic properties, density of states, the DFT X-ray photoelectron spectra and phonon dispersions calculated using the different methodologies. We demonstrate that GGA + U reproduces the AFM ordering in UN, but the crystal structure is dynamically unstable. We also show that magnetic ordering is important in finding the lowest energy defective structure, and that SOC has a distinct influence on the energy of the different uranium interstitial defects. Lastly, we discuss the point defect formation energies under U-rich and N-rich conditions, and the stoichiometric formation energies calculated with the different methodologies, providing insight into the observed tendency for forming hypostoichiometric UN.

36 MATERIALS SCIENCE↗

Emulation of neutron-irradiated microstructure of austenitic 21Cr32Ni model alloy using dual-ion irradiation

Here, in this study, the capability of heavy-ion irradiation to emulate neutron irradiation was demonstrated on an austenitic 21Cr32Ni type ternary model alloy. The model alloy used in this study is chemically analogous but compositionally simpler than of alloy 800H, which is a candidate austenitic Ni alloy which has been proposed for use in Generation IV reactors. The microstructure of the 21Cr32Ni model alloy irradiated in the BOR-60 fast reactor to 17.1 dpa and 35 dpa at ~380°C was characterized using transmission electron microscopy (TEM). The 17.1 dpa BOR-60 irradiated microstructure was then compared with the microstructure of the same material developed under dual-ion (DI) irradiation using various He/dpa ratios between 0.1 and 16.6 appm He/dpa in the temperature range of 430°C-500°C. The results showed that both neutron and DI irradiation of 21Cr32Ni model alloy produced dislocations in the form of a dislocation network as well as {111}-type faulted dislocation loops, cavities, and radiation-induced Ni enrichment at radiation-induced sinks. When the dose and the He/dpa ratio were kept similar to those in neutron irradiation, DI irradiation of the 21Cr32Ni model alloy at 460°C resulted in over-nucleation of small cavities and in a high density of faulted dislocation loops compared to those observed in the fast-neutron irradiated alloy of the same heat irradiated at ~380°C. The optimal condition for reproducing the neutron-irradiated microstructure was DI irradiation at 460°C and 0.1 appm He/dpa. In that case, the faulted loop and cavity size distributions in the BOR-60 irradiated 21Cr32Ni model alloy samples closely matched with those measured in the DI irradiated 21Cr32Ni model alloy sample. The fact that the He/dpa is an order of magnitude smaller than the helium generation rate for fast neutron irradiation, stops over nucleation and allows for the development of a similar microstructure as for neutron irradiation.

21Cr32Ni model alloy↗

Particle-Phase Diffusion Modulates Partitioning of Semivolatile Organic Compounds to Aged Secondary Organic Aerosol

The diffusivity of semivolatile organic compounds (SVOCs) in the bulk particle phase of a viscous atmospheric secondary organic aerosol (SOA) can have a profound impact on aerosol growth and size distribution dynamics. Here, we investigate the bulk diffusivity of SVOCs formed from photo-oxidation of isoprene as they partition to a bimodal aerosol consisting of an Aitken (potassium sulfate) and accumulation mode (aged α-pinene SOA) particles as a function of relative humidity (RH). The model analysis of the observed size distribution evolution shows that liquid-like diffusion coefficient values of D b > 10 –10 cm 2 s –1 fail to explain the growth of the Aitken mode. Instead, much lower values of Db between 2.5 × 10 –15 cm 2 s –1 at 32% RH and 8 × 10 –15 cm 2 s –1 at 82% RH were needed to successfully reproduce the growth of both modes. The diffusivity within the aged α-pinene SOA remains appreciably slow even at 80% RH, resulting in hindered partitioning of SVOCs to large viscous particles and allowing smaller and relatively less viscous particles to effectively absorb the available SVOCs and grow much faster than would be possible otherwise. Our results have important implications for modeling SOA formation and growth in the ambient atmosphere.

54 ENVIRONMENTAL SCIENCES↗

The cholesterol transporter NPC1 is essential for epigenetic regulation and maturation of oligodendrocyte lineage cells

The intracellular cholesterol transporter NPC1 functions in late endosomes and lysosomes to efflux unesterified cholesterol, and its deficiency causes Niemann–Pick disease Type C, an autosomal recessive lysosomal disorder characterized by progressive neurodegeneration and early death. Here, we use single-nucleus RNA-seq on the forebrain of Npc1 –/– mice at P16 to identify cell types and pathways affected early in pathogenesis. Our analysis uncovers significant transcriptional changes in the oligodendrocyte lineage during developmental myelination, accompanied by diminished maturation of myelinating oligodendrocytes. We identify upregulation of genes associated with neurogenesis and synapse formation in Npc1 –/– oligodendrocyte lineage cells, reflecting diminished gene silencing by H3K27me3. Npc1 –/ – oligodendrocyte progenitor cells reproduce impaired maturation in vitro, and this phenotype is rescued by treatment with GSK-J4, a small molecule inhibitor of H3K27 demethylases. Moreover, mobilizing stored cholesterol in Npc1 –/– mice by a single administration of 2-hydroxypropyl-β-cyclodextrin at P7 rescues myelination, epigenetic marks, and oligodendrocyte gene expression. Our findings highlight an important role for NPC1 in oligodendrocyte lineage maturation and epigenetic regulation, and identify potential targets for therapeutic intervention.

59 BASIC BIOLOGICAL SCIENCES↗

Half-hard copper PTW model parameters for plastic deformation calibrated to strain rates from 10 −3 -10 7 /s

This document demonstrates a calibration of the Preston-Tonks-Wallace (PTW) model of plastic deformation that fairly accurately reproduces the flow stress of half-hard copper across strain rates varying by a factor of more than ten billion. The experiments included are quasistatic compression at low rates, Split Hopkinson Pressure Bar (SHPB) compression at intermediate rates, and arrested Richtmyer-Meshkov instability (RMI) growth at high strain rates. The temperatures vary from liquid nitrogen (77 K) to around 500 K. The largest true strains obtained are approximately 0.5. It would be interesting to compare the results to copper deformation data at higher pressures and higher strains.

36 MATERIALS SCIENCE↗

Calculation of mode lifetimes in weakly anharmonic solids using self-consistent ensemble eigenstates of the Liouvillian

Recent schemes for the efficient calculation of vibrational mode frequencies and lifetimes in weakly anharmonic systems such as solid crystal structures have shown promise, with simple relationships being observed between low order moments of the Liouvillian and the observed lifetime. While only shown empirically, two parameters seemed sufficient to determine the lifetime to high accuracy over a large range of temperatures for a few simple low degree of freedom models and a Lennard–Jones solid. In this paper we analytically reproduce the numerically observed relationship using a few generally applicable assumptions about the density of states. This advance is made by studying the autocorrelation based on an “ensemble eigenstate” of the Liouvillian. We propose the use of this ensemble eigenstate as the basis for autocorrelations to simplify calculations in richer systems. Finally, an analysis is also performed of the relationship between canonical and microcanonical autocorrelation functions and a long standing error in the long time behavior of such functions originally presented by Sen, Sinkovits, and Chakravarti is corrected.

36 MATERIALS SCIENCE↗

Combined analysis of HPK 3.1 LGADs using a proton beam, beta source, and probe station towards establishing high volume quality control

The upgrades of the CMS and ATLAS experiments for the high luminosity phase of the Large Hadron Collider will employ precision timing detectors based on Low Gain Avalanche Detectors (LGADs). We present a suite of results combining measurements from the Fermilab Test Beam Facility, a beta source telescope, and a probe station, allowing full characterization of the HPK type 3.1 production of LGAD prototypes developed for these detectors. We demonstrate that the LGAD response to high energy test beam particles is accurately reproduced with a beta source. Finally, we further establish that probe station measurements of the gain implant accurately predict the particle response and operating parameters of each sensor, and conclude that the uniformity of the gain implant in this production is sufficient to produce full-sized sensors for the ATLAS and CMS timing detectors.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

How Do Climate Model Resolution and Atmospheric Moisture Affect the Simulation of Unprecedented Extreme Events Like the 2021 Western North American Heat Wave?

Abstract Although the 2021 Western North America (WNA) heat wave was predicted by weather forecast models, questions remain about whether such strong events can be simulated by global climate models (GCMs) at different model resolutions. Here, we analyze sets of GCM simulations including historical and future periods to check for the occurrence of similar events. High‐ and low‐resolution simulations both encounter challenges in reproducing events as extreme as the observed one, particularly under the present climate. Relatively stronger amplitudes are observed during the future periods. Furthermore, high‐ and low‐resolution short initialized GCM simulations are both able to reasonably predict such strong events and their associated high‐pressure ridge over the WNA with a 1 week forecast lead time. Moisture sensitivity experiments further indicate a drier atmospheric moisture condition results in substantially higher near‐surface temperatures in the simulated heat events.

54 ENVIRONMENTAL SCIENCES↗