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 1,189 records · Page 66
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

A Process-based Model with Temperature, Water, and Lab-derived Data Improves Predictions of Daily Culex pipiens/restuans Mosquito Density

While the number of human cases of mosquito-borne diseases has increased in North America in the last decade, accurate modeling of mosquito population density has remained a challenge. Longitudinal mosquito trap data over the many years needed for model calibration, and validation is relatively rare. In particular, capturing the relative changes in mosquito abundance across seasons is necessary for predicting the risk of disease spread as it varies from year to year. We developed a discrete, semi-stochastic, mechanistic process-based mosquito population model that captures life-cycle egg, larva, pupa, adult stages, and diapause for Culex pipiens (Diptera, Culicidae) and Culex restuans (Diptera, Culicidae) mosquito populations. This model combines known models for development and survival into a fully connected age-structured model that can reproduce mosquito population dynamics. Mosquito development through these stages is a function of time, temperature, daylight hours, and aquatic habitat availability. The time-dependent parameters are informed by both laboratory studies and mosquito trap data from the Greater Toronto Area. The model incorporates city-wide water-body gauge and precipitation data as a proxy for aquatic habitat. This approach accounts for the nonlinear interaction of temperature and aquatic habitat variability on the mosquito life stages. We demonstrate that the full model predicts the yearly variations in mosquito populations better than a statistical model using the same data sources. This improvement in modeling mosquito abundance can help guide interventions for reducing mosquito abundance in mitigating mosquito-borne diseases like West Nile virus.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Characterization and rank assignment criteria for the anaerobic fungi (Neocallimastigomycota)

Establishing a solid taxonomic framework is crucial for enabling discovery and documentation efforts. This ensures effective communication between scientists as well as reproducibility of results between laboratories, and facilitates the exchange and preservation of biological material. Such framework can only be achieved by establishing clear criteria for taxa characterization and rank assignment. Within the anaerobic fungi (phylum Neocallimastigomycota), the need for such criteria is especially vital. Difficulties associated with their isolation, maintenance and long-term storage often result in limited availability and loss of previously described taxa. To this end, we provide here a list of morphological, microscopic, phylogenetic and phenotypic criteria for assessment and documentation when characterizing newly obtained Neocallimastigomycota isolates. We also recommend a polyphasic rank-assignment scheme for novel genus-, species- and strain-level designations for newly obtained Neocallimastigomycota isolates.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Tropospheric aerosols over the western North Atlantic Ocean during the winter and summer deployments of ACTIVATE 2020: life cycle, transport, and distribution

The Aerosol Cloud meTeorology Interactions oVer the western ATlantic Experiment (ACTIVATE) is a NASA mission to characterize aerosol–cloud interactions over the western North Atlantic Ocean (WNAO). Such characterization requires understanding of life cycle, composition, transport pathways, and distribution of aerosols over the WNAO. This study uses the GEOS-Chem model to simulate aerosol distributions and properties that are evaluated against aircraft, ground-based, and satellite observations during the winter and summer field deployments in 2020 of ACTIVATE. Transport in the boundary layer (BL) behind cold fronts was a major mechanism for the North American continental outflow of pollution to the WNAO in winter. Turbulent mixing was the main driver for the upward transport of sea salt within and ventilation out of BL in winter. The BL aerosol composition was dominated by sea salt, which increased in the summer, followed by organics and sulfate. Aircraft in situ aerosol measurements provided useful constraints on wet scavenging in GEOS-Chem. The model generally captured observed features such as continental outflow, land–ocean gradient, and mixing of anthropogenic aerosols with sea salt. Model sensitivity experiments with elevated smoke injection heights to the mid-troposphere (versus within BL) better reproduced observations of smoke aerosols from the western US wildfires over the WNAO in the summer. Model analysis suggests strong hygroscopic growth of sea salt particles and their seeding of marine BL clouds over the WNAO (< 35° N). Future modeling efforts should focus on improving parameterizations for aerosol wet scavenging, implementing realistic smoke injection heights, and applying high-resolution models that better resolve vertical transport.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Diurnal cycle of precipitation over global monsoon systems in CMIP6 simulations

Accurate simulation of the diurnal cycle of precipitation over different climate regimes, including global monsoon systems, is an ongoing challenge for current state-of-the-art climate models. Here, in this study, the diurnal cycles of summertime precipitation from 21 CMIP6 historical simulations are evaluated with the NASA IMERG data via the diurnal harmonic amplitude and phase from the Fourier analysis in seven regional monsoon systems: Northern and Southern Africa monsoons, South Asia monsoon, East Asia monsoon, Australia monsoon, and North and South America monsoons. Particularly, we focus on four subregions including ocean, coastal ocean, land, and coastal land considering the different underlying physical processes. It is found that the CMIP6 historical simulations perform fairly well on reproducing the diurnal cycle of precipitation over ocean and coastal ocean, while the simulated diurnal harmonic phase and amplitude over land and coastal land show larger biases and model spread compared to biases over ocean and coastal ocean. The simulation of precipitation diurnal cycle is substantially improved in the multi-model mean of both CMIP6 historical and CMIP6 AMIP simulations compared with that of CMIP5 historical simulations. But the model bias of precipitation diurnal cycle is in general similar between CMIP6 AMIP and CMIP6 historical simulations, which suggests that the impact of the interactive ocean on the simulated diurnal precipitation averaged over a large domain is generally minor. The influence of model resolution on the diurnal cycle is shown in mixed directions across various subregions, and significant improvement from the low- to high-resolution models is only notable over coastal ocean.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

A FAIR-compliant parts catalogue for genome engineering and expression control in Saccharomyces cerevisiae

The synthetic biology toolkit for baker's yeast, Saccharomyces cerevisiae, includes extensive genome engineering toolkits and parts repositories. However, with the increasing complexity of engineering tasks and versatile applications of this model eukaryote, there is a continued interest to expand and diversify the rational engineering capabilities in this chassis by FAIR (findable, accessible, interoperable, and reproducible) compliance. In this study, we designed and characterised 41 synthetic guide RNA sequences to expand the CRISPR-based genome engineering capabilities for easy and efficient replacement of genomically encoded elements. Moreover, we characterize in high temporal resolution 20 native promoters and 18 terminators using fluorescein and LUDOX CL-X as references for GFP expression and OD600 measurements, respectively. Additionally, all data and reported analysis is provided in a publicly accessible jupyter notebook providing a tool for researchers with low-coding skills to further explore the generated data as well as a template for researchers to write their own scripts. We expect the data, parts, and databases associated with this study to support a FAIR-compliant resource for further advancing the engineering of yeasts.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Accelerating Finite-temperature Kohn-Sham Density Functional Theory\ with Deep Neural Networks

We present a numerical modeling workflow based on machine learning (ML) which reproduces the the total energies produced by Kohn-Sham density functional theory (DFT) at finite electronic temperature to within chemical accuracy at negligible computational cost. Based on deep neural networks, our workflow yields the local density of states (LDOS) for a given atomic configuration. From the LDOS, spatially-resolved, energy-resolved, and integrated quantities can be calculated, including the DFT total free energy, which serves as the Born-Oppenheimer potential energy surface for the atoms. We demonstrate the efficacy of this approach for both solid and liquid metals and compare results between independent and unified machine-learning models for solid and liquid aluminum. Our machine-learning density functional theory framework opens up the path towards multiscale materials modeling for matter under ambient and extreme conditions at a computational scale and cost that is unattainable with current algorithms.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Understanding flavin electronic structure and spectra

Abstract Flavins have emerged as central to electron bifurcation, signaling, and countless enzymatic reactions. In bifurcation, two electrons acquired as a pair are separated in coupled transfers wherein the energy of both is concentrated on one of the two. This enables organisms to drive demanding reactions based on abundant low‐grade chemical fuel. To enable incorporation of this and other flavin capabilities into designed materials and devices, it is essential to understand fundamental principles of flavin electronic structure that make flavins so reactive and tunable by interactions with protein. Emerging computational tools can now replicate spectra of flavins and are gaining capacity to explain reactivity at atomistic resolution, based on electronic structures. Such fundamental understanding can moreover be transferrable to other chemical systems. A variety of computational innovations have been critical in reproducing experimental properties of flavins including their electronic spectra, vibrational signatures, and nuclear magnetic resonance (NMR) chemical shifts. A computational toolbox for understanding flavin reactivity moreover must be able to treat all five oxidation and protonation states, in addition to excited states that participate in flavoprotein's light‐driven reactions. Therefore, we compare emerging hybrid strategies and their successes in replicating effects of hydrogen bonding, the surrounding dielectric, and local electrostatics. These contribute to the protein's ability to modulate flavin reactivity, so we conclude with a survey of methods for incorporating the effects of the protein residues explicitly, as well as local dynamics. Computation is poised to elucidate the factors that affect a bound flavin's ability to mediate stunningly diverse reactions, and make life possible. This article is categorized under: Structure and Mechanism > Computational Biochemistry and Biophysics Electronic Structure Theory > Combined QM/MM Methods Theoretical and Physical Chemistry > Spectroscopy

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Sprain energy consequences for damage localization and fracture mechanics

The 2023 smooth Lagrangian Crack-Band Model (slCBM), inspired by the 2020 invention of the gap test, prevented spurious damage localization during fracture growth by introducing the second gradient of the displacement field vector, named the “sprain,” as the localization limiter. The key idea was that, in the finite element implementation, the displacement vector and its gradient should be treated as independent fields with the lowest ( C 0 ) continuity, constrained by a second-order Lagrange multiplier tensor. Coupled with a realistic constitutive law for triaxial softening damage, such as microplane model M7, the known limitations of the classical Crack Band Model were eliminated. Here, we show that the slCBM closely reproduces the size effect revealed by the gap test at various crack-parallel stresses. To describe it, we present an approximate corrective formula, although a strong loading-path dependence limits its applicability. Except for the rare case of zero crack-parallel stresses, the fracture predictions of the line crack models (linear elastic fracture mechanics, phase-field, extended finite element method (XFEM), cohesive crack models) can be as much as 100% in error. We argue that the localization limiter concept must be extended by including the resistance to material rotation gradients. We also show that, without this resistance, the existing strain-gradient damage theories may predict a wrong fracture pattern and have, for Mode II and III fractures, a load capacity error as much as 55%. Finally, we argue that the crack-parallel stress effect must occur in all materials, ranging from concrete to atomistically sharp cracks in crystals.

Science & Technology - Other Topics↗
DOE OSTI2021

Qualification of 3-D Printed Mortar With Electrical Conductivity Measurements

Additive manufacturing (AM) or 3D printing of concrete allows for construction of arbitrary shape structures without a mold. Since reproducibility of 3D printed concrete is lower than that of conventional fabrication, each 3D printed structure should be monitored for proper curing. Conventional qualification of concrete is based on several tests, including destructive compressive strength measurements. Because a structure is 3-D printed layer-by-layer, the surfaces of AM concrete structures have significant surface roughness. This limits the applicability of conventional nondestructive testing methods. We investigated qualification of 3D printed mortar by monitoring curing with nondestructive electrical conductivity measurements. Bulk resistance of concrete was extracted from electrochemical impedance spectroscopy (EIS measurements made using custom flexible self-adhesive electrodes, which contour to rough surfaces. We show that bulk resistivity of concrete increases linearly with time. This allows for developing a calibration curve for compressive strength lookup from nondestructive electrical conductivity measurements. Conductivity measurements also allow for estimation of formation factor, which is an indicator of mortar permeability.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Measurement of isolated prompt photon production in pp and p–Pb collisions at the LHC

This paper presents the measurement of the isolated prompt photon inclusive production cross section in pp and p–Pb collisions by the ALICE Collaboration at the LHC. The measurement is performed in p–Pb collisions at centre-of-mass energies per nucleon pair of $\sqrt{s_{NN}}$ = 5.02 Te and 8.16 TeV, as well as in pp collisions at $\sqrt{s}$ = 5.02 TeV and 8 TeV. The cross section is obtained at midrapidity (|γ| < 0.7) using a charged-track based isolation momentum p$^{iso,~ch}_{T}$ < 1.5 GeV/c in a cone with radius R = 0.4. The data for both collision systems are well reproduced by perturbative QCD (pQCD) calculations at next-to-leading order (NLO) using recent parton distribution functions for free (PDF) and bound (nPDF) nucleons. Furthermore, the nuclear modification factor R pA for both collision energies is consistent with unity for p T > 20 GeV/c. However, deviations from unity (R pA < 1) of up to 20% are observed for p T < 20 GeV/c with limited significance, indicating the possible presence of nuclear effects in the initial state of the collision.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

FtsZ-mediated fission of a cuboid bacterial symbiont

Less than a handful of cuboid and squared cells have been described in nature, which makes them a rarity. Here, we show how Candidatus Thiosymbion cuboideus, a cube-like gammaproteobacterium, reproduces on the surface of marine free-living nematodes. Immunostaining of symbiont cells with an anti-fimbriae antibody revealed that they are host-polarized, as these appendages exclusively localized at the host-proximal (animal-attached) pole. Moreover, by applying a fluorescently labeled metabolic probe to track new cell wall insertion in vivo, we observed that the host-attached pole started septation before the distal one. Similarly, Ca. T. cuboideus cells immunostained with an anti-FtsZ antibody revealed a proximal-to-distal localization pattern of this tubulin homolog. Although FtsZ has been shown to arrange into squares in synthetically remodeled cuboid cells, here we show that FtsZ may also mediate the division of naturally occurring ones. This implies that, even in natural settings, membrane roundness is not required for FtsZ function.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Hexagonal Lattices of HIV Capsid Proteins Explored by Simulations Based on a Thermodynamically Consistent Model

HIV capsid proteins (CAs) may self-assemble into a variety of shapes under in vivo and in vitro conditions. Here, we employed simulations based on a residue-level coarse-grained (CG) model with full conformational flexibility to investigate hexagonal lattices, which are the underlying structural pattern for CA aggregations. Facilitated by enhanced sampling simulations to rigorously calculate CA dimerization and polymerization affinities, we calibrated our model to reproduce the experimentally measured affinities. Here using the calibrated model, we performed unbiased simulations on several large systems consisting of 1512 CA subunits, allowing reversible binding and unbinding of the CAs in a thermodynamically consistent manner. In one simulation, a preassembled hexagonal CA sheet developed spontaneous curvatures reminiscent of those observed in experiments, and the edges of the sheet exhibited local curvatures larger than those of the interior. In other simulations starting with randomly distributed CAs at different concentrations, existing CA assemblies grew by binding free capsomeres to the edges and by merging with other assemblies. At high CA concentrations, rapid establishment of predominant aggregates was followed by much slower adjustments toward more regular hexagonal lattices, with increasing numbers of intact CA hexamers and pentamers being formed. Our approach of adapting a general CG model to specific systems by using experimental binding data represents a practical and effective strategy for simulating and elucidating intricate protein aggregations.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Reproducibility of Radiokrypton in Deep Desert Aquifers: Insights from a Decade of Research

Great technical advances have been achieved since the first atom-trap trace analysis (ATTA) -based radiokrypton application in Egypt, where 1 Myr old groundwater was discovered. Beyond advances in ATTA measurement capabilities, including reduction in sample size, analysis duration, and analytical uncertainty, major progress has been achieved over the past two decades in the sample collection and preparation techniques. These advances paved the expansion of ATTA-based noble gas applications to many other aquifers worldwide, illuminating the nature and flow pattern of deep groundwater systems. While the potential of this new analytical technique for old groundwater dating is well recognized, another important aspect yet to be examined is the reproducibility of radiokrypton in aquifers over time, i.e., how representative is a discrete groundwater sample, collected at a specific time and location, for the natural groundwater system? The likelihood of a negative answer is increased by flow-field disturbance in aquifers following massive groundwater abstraction. Here, in this work, we present repeated 81 Kr sampling and measurements in twenty-one sites over Israel, mostly of deep (up to 1 km) wells tapping confined aquifers in the arid to hyperarid Negev desert. The results demonstrate that radiokrypton measurements are indeed reproducible, even in cases where samples were collected as long as nine years apart and from highly productive (∼1 Mm 3 /yr order) pumping wells. Furthermore, many of the repeated measurements in this study (17 out of the 21 sites) were conducted with different ATTA Instruments in two different laboratories using slightly different sampling, preparation, and analysis techniques, yet with an overall good agreement. The consistency in the ATTA-based 81 Kr-dating results over time highlights the robustness of this state-of-the-art technique as a tool to unravel groundwater flow patterns and encourages further applications to many other yet-to-be-explored deep aquifers.

atom-trap trace analysis↗
DOE OSTI2022

Non‐Linear Dimensionality Reduction With a Variational Encoder Decoder to Understand Convective Processes in Climate Models

Abstract Deep learning can accurately represent sub‐grid‐scale convective processes in climate models, learning from high resolution simulations. However, deep learning methods usually lack interpretability due to large internal dimensionality, resulting in reduced trustworthiness in these methods. Here, we use Variational Encoder Decoder structures (VED), a non‐linear dimensionality reduction technique, to learn and understand convective processes in an aquaplanet superparameterized climate model simulation, where deep convective processes are simulated explicitly. We show that similar to previous deep learning studies based on feed‐forward neural nets, the VED is capable of learning and accurately reproducing convective processes. In contrast to past work, we show this can be achieved by compressing the original information into only five latent nodes. As a result, the VED can be used to understand convective processes and delineate modes of convection through the exploration of its latent dimensions. A close investigation of the latent space enables the identification of different convective regimes: (a) stable conditions are clearly distinguished from deep convection with low outgoing longwave radiation and strong precipitation; (b) high optically thin cirrus‐like clouds are separated from low optically thick cumulus clouds; and (c) shallow convective processes are associated with large‐scale moisture content and surface diabatic heating. Our results demonstrate that VEDs can accurately represent convective processes in climate models, while enabling interpretability and better understanding of sub‐grid‐scale physical processes, paving the way to increasingly interpretable machine learning parameterizations with promising generative properties.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Hybrid molecular beam epitaxy growth of BaTiO 3 films

The ability to reproducibly synthesize thin films with precise composition and controlled structure is essential for fundamental study and mass production. Here, we demonstrate the hybrid molecular beam epitaxy (MBE) growth of epitaxial, single crystalline BaTiO3 films with different thicknesses on Nb-doped SrTiO3 substrates with atomically smooth surfaces. By combining scanning transmission electron microscopy, temperature-dependent high-resolution x-ray diffraction, reflection high-energy electron diffraction, and atomic force microscopy, we study the effect of growth conditions and the interplay between stoichiometry and epitaxial strain on the resulting structure. Furthermore, we demonstrate a close to bulk-like ferroelectric phase transition in thicker films and highlight the effect of strain on the phase transition temperature. This work establishes the hybrid MBE approach for the growth of heteroepitaxial BaTiO3 films on conducting substrates with scalable thickness and controlled stoichiometry.

Materials Science↗
DOE OSTI2026

Hierarchical Reinforcement Learning of a Short-Range Bond-Order Potential for Silica: Analytic Embedding of Coordination with Classical Efficiency

Reinforcement learning (RL) has recently emerged as a data-efficient strategy to parametrize short-range interatomic potentials. Building on our past RL optimization of pairwise silica models, we extend the framework to a bond-order (Tersoff-type) potential that provides an analytic embedding of local coordination through a three-body term. A hierarchical RL workflow combining continuous-action Monte Carlo Tree Search and property-based rewards efficiently explores the 26-dimensional parameter space, sequentially optimizing lattice parameters, densities, angles, and cohesive energies of 21 silica polymorphs. The resulting models, Q-Tersoff and ML-Tersoff, reproduce the energetic ordering of low-energy phases and capture the angular correlations and amorphous structure factors of silica with improved fidelity over pairwise force fields, while remaining orders of magnitude faster than high-dimensional machine-learned potentials. Both models underperform for elastic constants and high-energy frameworks, delineating the limits of the current analytic form. The approach establishes a general and interpretable route to angle-aware, short-range potentials that bridge physics-based and machine-learned descriptions of silicate materials.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Electric dipole polarizability of 58 Ni

The electric dipole strength distribution in 58 Ni between 6 and 20 MeV has been determined from proton inelastic scattering experiments at very forward angles at RCNP, Osaka. The experimental data are rather well reproduced by quasiparticle random-phase approximation calculations including vibration coupling, despite a mild dependence on the adopted Skyrme interaction. They allow an estimate of the experimentally inaccessible high-energy contribution above 20 MeV, leading to an electric dipole polarizability α D ⁡ ( 58 Ni) = 3.48 (31)⁢ fm 3 . This serves as a test case for recent extensions of coupled-cluster calculations with chiral effective field theory interactions to nuclei with two nucleons on top of a closed-shell system.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Contrasting Stratospheric Smoke Mass and Lifetime From 2017 Canadian and 2019/2020 Australian Megafires: Global Simulations and Satellite Observations

Stratospheric injections of carbonaceous aerosols and combustion gases by extreme wildfires have become increasingly common. Recent “megafires,” particularly large and intense fires, delivered particulate burdens to the lower stratosphere comparable to those of moderate volcanic eruptions. The 2017 Canadian megafire generated four large Pyrocumulonimbi (pyroCbs), injecting up to ≈0.3 Tg of smoke in the lower stratosphere. Even more extreme, the 2019/2020 Australian event produced a pyroCb activity resulting in stratospheric smoke intrusions of ≈1 Tg. To understand their contrasting behavior, we present global climate simulations of the atmospheric response to these events, applying smoke burdens informed by remote observations. Model outcomes, compared to satellite data of smoke transport, reproduce reasonably well the initial plume rise, at 0.2–0.3 km/day, attaining heights of ≈20 km in Canada and above 30 km in Australia. Global dispersal of the plume occurs within about 3 weeks in both cases, consistent with observations. Smoke removal timescales, ≈5 months for the Canadian megafire, agree with remote measurements. During the Australian megafire, observations indicate stratospheric injections three times as large, and models predict comparatively longer smoke lifetimes, ≈16 months. After the latter event, atmospheric optical depths and radiative cooling achieved values close to those measured following the Pinatubo eruption. Sensitivity tests of model assumptions indicate, in accord with prior studies, that smoke burden, injection heights, and black carbon content can determine plume evolution and possible climate impacts. An empirical relation between peak heights of stratospheric plumes and lifetimes is derived that can help assess megafire impacts on the stratosphere, climate and the Earth system

54 ENVIRONMENTAL SCIENCES↗