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 973 records · Page 54

Maximizing efficiency of dataset compression for machine learning potentials with information theory

Machine learning interatomic potentials (MLIPs) balance high accuracy and lower costs compared to density functional theory calculations, but their performance often depends on the size and diversity of training datasets. Large datasets improve model accuracy and generalization but are computationally expensive to produce and train on, while smaller datasets risk discarding rare but important atomic environments and compromising MLIP accuracy/reliability. Here, we develop an information-theoretical framework to quantify the efficiency of dataset compression methods and propose an algorithm that maximizes this efficiency. By framing atomistic dataset compression as an instance of the minimum set cover (MSC) problem over atom-centered environments, our method identifies the smallest subset of structures that contains as much information as possible from the original dataset while pruning redundant information. The approach is extensively demonstrated on the GAP-20 and TM23 datasets and validated on 64 varied datasets from the ColabFit repository. Across all cases, MSC consistently retains outliers, preserves dataset diversity, and reproduces the long-tail distributions of forces even at high compression rates, outperforming other subsampling methods. Furthermore, MLIPs trained on MSC-compressed datasets exhibit reduced error for out-of-distribution data even in low-data regimes. We explain these results using an outlier analysis and show that such quantitative conclusions could not be achieved with conventional dimensionality reduction methods. The algorithm is implemented in the open-source QUESTS package and can be used for several tasks in atomistic modeling, from data subsampling, outlier detection, and training improved MLIPs at a lower cost.

36 MATERIALS SCIENCE↗

Landscape fragmentation overturns classical metapopulation thinking

Habitat loss and isolation caused by landscape fragmentation represent a growing threat to global biodiversity. Existing theory suggests that the process will lead to a decline in metapopulation viability. However, since most metapopulation models are restricted to simple networks of discrete habitat patches, the effects of real landscape fragmentation, particularly in stochastic environments, are not well understood. To close this major gap in ecological theory, we developed a spatially explicit, individual-based model applicable to realistic landscape structures, bridging metapopulation ecology and landscape ecology. This model reproduced classical metapopulation dynamics under conventional model assumptions, but on fragmented landscapes, it uncovered general dynamics that are in stark contradiction to the prevailing views in the ecological and conservation literature. Notably, fragmentation can give rise to a series of dualities: a) positive and negative responses to environmental noise, b) relative slowdown and acceleration in density decline, and c) synchronization and desynchronization of local population dynamics. Furthermore, counter to common intuition, species that interact locally (“residents”) were often more resilient to fragmentation than long-ranging “migrants.” This set of findings signals a need to fundamentally reconsider our approach to ecosystem management in a noisy and fragmented world.

54 ENVIRONMENTAL SCIENCES↗

Spectroscopy-guided discovery of three-dimensional structures of disordered materials with diffusion models

Spectroscopy techniques such as x-ray absorption near edge structure (XANES) provide valuable insights into the atomic structures of materials, yet the inverse prediction of precise structures from spectroscopic data remains a formidable challenge. In this study, we introduce a framework that combines generative artificial intelligence models with XANES spectroscopy to predict three-dimensional atomic structures of disordered systems, using amorphous carbon (a-C) as a model system. In this work, we introduce a new framework based on the diffusion model, a recent generative machine learning method, to predict 3D structures of disordered materials from a target property. For demonstration, we apply the model to identify the atomic structures of a-C as a representative material system from the target XANES spectra. We show that conditional generation guided by XANES spectra reproduces key features of the target structures. Furthermore, we show that our model can steer the generative process to tailor atomic arrangements for a specific XANES spectrum. Finally, our generative model exhibits a remarkable scale-agnostic property, thereby enabling generation of realistic, large-scale structures through learning from a small-scale dataset (i.e. with small unit cells). Our work represents a significant stride in bridging the gap between materials characterization and atomic structure determination; in addition, it can be leveraged for materials discovery in exploring various material properties as targeted.

36 MATERIALS SCIENCE↗

Importance of bulk excitations and coherent electron-photon-phonon scattering in photoemission from PbTe(111): Ab initio theory with experimental comparisons

This paper presents a fully ab initio many-body photoemission framework that includes coherent three-body electron-photon-phonon scattering to predict the transverse momentum distributions and the mean transverse energies (MTEs) of bulk photoelectrons from single-crystal photocathodes. The need to develop such a theory stems from the lack of studies that provide complete understanding of the underlying fundamental processes governing the transverse momentum distribution of photoelectrons emitted from single crystals. For example, initial predictions based on density-functional theory calculations of effective electron masses suggested that the (111) surface of PbTe would produce very small MTEs (≤ 15 meV), whereas our experiments yielded MTEs ten to twenty times larger than these predictions and also exhibited a lower photoemission threshold than predicted. The ab initio framework presented in this paper correctly reproduces both the magnitude of the MTEs from our measurements in PbTe(111) and the observed photoemission below the predicted threshold. Our results show that both photoexcitations into states that propagate in the bulk of the material and coherent many-body electron-photon-phonon scattering processes, which initial predictions ignored, play surprisingly important roles in photoemission from PbTe(111). Lastly, from the lessons learned, we recommend a procedure for rapid computational screening of potential single-crystal photocathodes for applications in next-generation ultrafast electron diffraction and X-ray free-electron lasers, which will enable significant advances in condensed matter research.

36 MATERIALS SCIENCE↗

Beirut Explosion Yield and Mushroom Cloud Height - Effects of the Source Environment

I use crater dimensions to estimate the yield of the August 4th, 2020 Beirut explosion to be equivalent to approximately 1.4 kilotons of TNT with a lower bound of about 0.7 kilotons. Based on the amount of ammonium nitrate reported to have been stored at the Beirut harbor, I assume an upper bound for the yield of 2.75 kilotons. However, it is highly likely that the yield was less than 2.75 kilotons, since reported values for TNT equivalence of ammonium nitrate are typically much less than one hundred percent. The crater-size based yield estimates are based on crater radius estimates from satellite imagery and empirical curves and data for scaled crater radius from past chemical and nuclear explosions. I present evidence that suggests that the relatively large crater radius is due to a high degree of coupling of shock wave energy to the surrounding medium and a reduction of the effective stress because of a high level of saturation of the geologic media beneath the explosion. I provide yield estimates based on seismic body-wave magnitude and crater depth as corroborating evidence. I compare preliminary estimates for the maximum debris cloud height, based on cell phone videos/images, with predicted maximum heights for this yield range from empirical formulas and numerical cloud-rise models. Based on a preliminary analysis of cell phone footage, the observed maximum cloud height appears to be approximately 1600 m. This is much lower than that predicted using standard empirical formulas and buoyant cloud rise models. I present results from a modified buoyant cloud rise model that more accurately predicts the maximum cloud height by allowing for the inclusion of a fixed amount of air and/or water into the fireball at the start of cloud rise. The amount of mass that needs to be added at the start, to reproduce the observed maximum cloud height, is relatively small compared to the total mass entrained during cloud rise. A much greater amount of dry air or debris is required, relative to water, for an equivalent reduction in maximum cloud height. The ammonium nitrate is one possible source for water in the fireball since it was being stored in a very humid environment and ammonium nitrate is known to be hygroscopic. The ground beneath the explosion, especially if it were saturated, and the nearby harbor could also have been sources for water or debris in the fireball.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

Exploring the process-structure-property relationship of Aerosol Deposition to phosphor coatings for non-contact thermometry

Full-field, multi-measurand diagnostics provide rich validation data necessary to improve the product life cycle time of nuclear safety components. Thermophosphor digital image correlation (TP+DIC) is a method of simultaneously measuring strain and temperature fields using patterned phosphor coatings deposited with aerosol deposition (AD). While TP+DIC produces a functional diagnostic, the coating’s reproducibility and the effect of the patterned features on the inferred temperature remains uncharacterized. This NSR&D project provided the opportunity to study two areas: 1) the tunability and repeatability of aerosol deposition and 2) the robustness of aerosol deposition phosphor on deforming substrates. The first area explores the process-property relationship of parameters elucidating the significance of each on the coating. The second area explores the relationship between the features’ characteristics (namely thickness) and the phosphor emission and inferred temperature. Together, the results will lead to the improved accuracy and functionality of TP+DIC for qualification testing of nuclear safety components.

36 MATERIALS SCIENCE↗

Earth System Model Aerosol–Cloud Diagnostics (ESMAC Diags) package, version 2: assessing aerosols, clouds, and aerosol–cloud interactions via field campaign and long-term observations

Abstract. Poor representations of aerosols, clouds, and aerosol–cloud interactions (ACIs) in Earth system models (ESMs) have long been the largest uncertainties in predicting global climate change. Huge efforts have been made to improve the representation of these processes in ESMs, and the key to these efforts is the evaluation of ESM simulations with observations. Most well-established ESM diagnostics packages focus on the climatological features; however, they lack process-level understanding and representations of aerosols, clouds, and ACIs. In this study, we developed the Earth System Model Aerosol–Cloud Diagnostics (ESMAC Diags) package to facilitate the routine evaluation of aerosols, clouds, and ACIs simulated the Energy Exascale Earth System Model (E3SM) from the US Department of Energy (DOE). This paper documents its version 2 functionality (ESMAC Diags v2), which has substantial updates compared with version 1 (Tang et al., 2022a). The simulated aerosol and cloud properties have been extensively compared with in situ and remote-sensing measurements from aircraft, ship, surface, and satellite platforms in ESMAC Diags v2. It currently includes six field campaigns and two permanent sites covering four geographical regions: the eastern North Atlantic, the central US, the northeastern Pacific, and the Southern Ocean. These regions produce frequent liquid- or mixed-phase clouds, with extensive measurements available from the DOE Atmospheric Radiation Measurement user facility and other agencies. ESMAC Diags v2 generates various types of single-variable and multivariable diagnostics, including percentiles, histograms, joint histograms, and heatmaps, to evaluate the model representation of aerosols, clouds, and ACIs. Select examples highlighting the capabilities of ESMAC Diags are shown using E3SM version 2 (E3SMv2). In general, E3SMv2 can reasonably reproduce many observed aerosol and cloud properties, with biases in some variables such as aerosol particle and cloud droplet sizes and number concentrations. The coupling of aerosol and cloud number concentrations may be too strong in E3SMv2, possibly indicating a bias in processes that control aerosol activation. Furthermore, the liquid water path response to a perturbed cloud droplet number concentration behaves differently in E3SMv2 and observations, which warrants further study to improve the cloud microphysics parameterizations in E3SMv2.

54 ENVIRONMENTAL SCIENCES↗

Causality guided machine learning model on wetland CH 4 emissions across global wetlands

Wetland CH 4 emissions are among the most uncertain components of the global CH 4 budget. The complex nature of wetland CH 4 processes makes it challenging to identify causal relationships for improving our understanding and predictability of CH 4 emissions. In this study, we used the flux measurements of CH 4 from eddy covariance towers (30 sites from 4 wetlands types: bog, fen, marsh, and wet tundra) to construct a causality-constrained machine learning (ML) framework to explain the regulative factors and to capture CH 4 emissions at sub-seasonal scale. We found that soil temperature is the dominant factor for CH 4 emissions in all studied wetland types. Ecosystem respiration (CO 2 ) and gross primary productivity exert controls at bog, fen, and marsh sites with lagged responses of days to weeks. Integrating these asynchronous environmental and biological causal relationships in predictive models significantly improved model performance. More importantly, modeled CH 4 emissions differed by up to a factor of 4 under a +1°C warming scenario when causality constraints were considered. These results highlight the significant role of causality in modeling wetland CH 4 emissions especially under future warming conditions, while traditional data-driven ML models may reproduce observations for the wrong reasons. Our proposed causality-guided model could benefit predictive modeling, large-scale upscaling, data gap-filling, and surrogate modeling of wetland CH 4 emissions within earth system land models.

54 ENVIRONMENTAL SCIENCES↗

Nonaxisymmetric Shapes of Biological Membranes from Locally Induced Curvature

In various biological processes such as endocytosis and caveolae formation, the cell membrane is locally deformed into curved morphologies. Previous models to study membrane morphologies resulting from locally induced curvature often only consider the possibility of axisymmetric shapes-an indeed unphysical constraint. Past studies predict that the cell membrane buds at low resting tensions and stalls at a flat pit at high resting tensions. In this work, we lift the restriction to axisymmetry to study all possible membrane morphologies. Only if the resting tension of the membrane is low, we reproduce axisymmetric membrane morphologies. When the resting tension is moderate to high, we show that 1) axisymmetric membrane pits are unstable and 2) nonaxisymmetric ridge-shaped structures are energetically favorable. Furthermore, we find the interplay between intramembrane viscous flow and the rate of induced curvature affects the membrane's ability to transition into nonaxisymmetric ridges and axisymmetric buds. In particular, we show that axisymmetric buds are favored when the induced curvature is rapidly increased, whereas nonaxisymmetric ridges are favored when the curvature is slowly increased. Overall, our results hold relevant implications for biological processes such as endocytosis and physical phenomena like phase separation in lipid bilayers.

59 BASIC BIOLOGICAL SCIENCES↗

Modeling the Effects of Artificial Drainage on Agriculture-Dominated Watersheds Using a Fully Distributed Integrated Hydrology Model

In agriculture-dominated watersheds where natural drainage is poor, agricultural ditches (narrow engineered channels) and tile drains (perforated pipes) are widely employed to enhance surface and subsurface drainage, respectively. Despite their relatively small scale, these features exert substantial control over the hydro-biogeochemical function of watersheds and their effects need to be represented in the models. We introduce a novel strategy to incorporate the effects of artificial agricultural drainage into a fully distributed basin-scale integrated surface-subsurface hydrology models. In our approach, narrow agriculture ditches for surface drainage are resolved efficiently using ditch-aligned computational meshes that are hydrologically conditioned to ensure connectivity in the stream/ditch network. For tile drainage in the subsurface, we use the physically based Hooghoudt's drainage equation as a subgrid model and route the water drained through tiles to the nearest ditch. Without site-specific calibration, this model reproduced observed streamflow in the Portage River Watershed (>1,000 km 2 ) as recorded by a USGS gauge with good accuracy (normalized KGE = 0.81) and outperformed a calibrated SWAT model (normalized KGE = 0.68). Numerical experiments confirm that artificial drainage reduces surface inundations and effectively controls the water table. At the watershed scale, artificial drainage increases baseflow but has little effect on watershed discharges above the 90th percentile. The strong physical underpinnings and reduced need for calibration allow us to study the impacts of artificial drainage on distributed hydrological response in terms of fluxes and states and provide a platform for investigating watershed-scale nutrient transport.

54 ENVIRONMENTAL SCIENCES↗

MISPR : an open-source package for high-throughput multiscale molecular simulations

Computational tools provide a unique opportunity to study and design optimal materials by enhancing our ability to comprehend the connections between their atomistic structure and functional properties. However, designing materials with tailored functionalities is complicated due to the necessity to integrate various computational-chemistry software (not necessarily compatible with one another), the heterogeneous nature of the generated data, and the need to explore vast chemical and parameter spaces. The latter is especially important to avoid bias in scattered data points-based models and derive statistical trends only accessible by systematic datasets. Here, we introduce a robust high-throughput multi-scale computational infrastructure coined MISPR (Materials Informatics for Structure–Property Relationships) that seamlessly integrates classical molecular dynamics (MD) simulations with density functional theory (DFT). By enabling high-performance data analytics and coupling between different methods and scales, MISPR addresses critical challenges arising from the needs of automated workflow management and data provenance recording. The major features of MISPR include automated DFT and MD simulations, error handling, derivation of molecular and ensemble properties, and creation of output databases that organize results from individual calculations to enable reproducibility and transparency. In this work, we describe fully automated DFT workflows implemented in MISPR to compute various properties such as nuclear magnetic resonance chemical shift, binding energy, bond dissociation energy, and redox potential with support for multiple methods such as electron transfer and proton-coupled electron transfer reactions. The infrastructure also enables the characterization of large-scale ensemble properties by providing MD workflows that calculate a wide range of structural and dynamical properties in liquid solutions. MISPR employs the methodologies of materials informatics to facilitate understanding and prediction of phenomenological structure–property relationships, which are crucial to designing novel optimal materials for numerous scientific applications and engineering technologies.

36 MATERIALS SCIENCE↗

Changes in Future Precipitation Mean and Variability across Scales

Changes in precipitation variability can have large societal consequences, whether at the short time scales of flash floods or the longer time scales of multiyear droughts. Recent studies have suggested that in future climate projections, precipitation variability rises more steeply than does its mean, leading to concerns about societal impacts. Here we evaluate changes in mean precipitation over a broad range of spatial and temporal scales using a range of models from high-resolution regional simulations to millennial-scale global simulations. Results show that changes depend on the scale of aggregation and involve strong regional differences. On local scales that resolve individual rainfall events (hours and tens of kilometers), changes in precipitation distributions are complex and variances rise substantially more than means, as is required given the well-known disproportionate rise in precipitation intensity. On scales that aggregate across many events, distributional changes become simpler and variability changes smaller. At regional scale, future precipitation distributions can be largely reproduced by a simple transformation of present-day precipitation involving a multiplicative shift and a small additive term. The “extra” broadening is negatively correlated with changes in mean precipitation: in strongly “wetting” areas, distributions broaden less than expected from a simple multiplicative mean change; in “drying” areas, distributions narrow less. Precipitation variability changes are therefore of especial concern in the subtropics, which tend to dry under climate change. Outside the tropics, variability changes are similar on time scales from days to decades (i.e., show little frequency dependence). This behavior is highly robust across models, suggesting it may stem from some fundamental constraint.

54 ENVIRONMENTAL SCIENCES↗

Additive manufacturing and processing of architected materials

Architected materials are a unique and emerging class of materials where performance is fundamentally controlled by geometry at multiple length scales, from the nano- to the macroscale, rather than chemical composition alone. As a result, the realization of these remarkable materials is contingent upon the ability to faithfully reproduce the designed architecture. Furthermore, this presents fundamental challenges in fabrication due to the required three-dimensional complexity, multiple length scales, range of material constituents, possibility of multiple materials in a single architecture, and overall manufacturing throughput. Additive manufacturing (AM) processes can provide solutions to some of these challenges and are discussed in this article. Specifically, light-based and extrusion-based processes and associated materials are presented with an emphasis on recent developments, including volumetric additive manufacturing, and on-the-fly mixing of materials in extrusion-based printing systems. While remarkable advancements have been made in AM for architected materials, bringing these materials and processes to industrial realization remains a significant challenge.

36 MATERIALS SCIENCE↗

Alteration of phycobilisome excitation energy transfer properties in response to attenuations in peripheral electron flow

In Synechocystis sp. PCC 6803 ( S . 6803), two types of phycobilisome (PBS) complexes, CpcG-PBS and CpcL-PBS, function to harvest light energy for photosynthetic reaction centers (RCs), photosystem I (PSI) and photosystem II (PSII). The compositional differences between these two forms of PBS and their specificity for RCs have led to suggestions that they may differ in function. To address this question, we examined how PBS-RC interactions, and the transfer of excitation energy from PBS to RCs, might be adjusted under conditions where electron demand and photon availability are modulated. The CpcG-PBS, CpcL-PBS, and RC complexes were isolated from a S . 6803 strain defective in expression of flavodiiron 1 (oxygen reduction reaction 1, ORR1) grown under varied light regimes. The energy transfer preference from CpcL-PBS to either PSI or PSII was investigated by in vitro crosslinking and 77 K fluorescence emission spectroscopy to assess energy transfer efficiency under photoexcitation. While the results demonstrate that the transfer of excitation energy from CpcL-PBS favors PSI over PSII in WT strains as previously shown, the preference of CpcL-PBS switches from PSI to PSII in ORR1 strains. Surprisingly, this change in preference was reproduced when ORR1 CpcL-PBS was crosslinked with WT RCs, or when WT CpcL-PBS was cross-crosslinked with ORR1 RCs, indicating there are physical modifications to both PBS and RCs that mediate the preference switch. In contrast, the analysis with ORR1 CpcG-PBS shows similar preferences to WT. Additionally, PBS populations in ORR1 shifted to a greater proportion of CpcL-PBS relative to CpcG-PBS. These results demonstrate that under conditions where electron utilization changes, there is a tuning of the excitation energy allocation from CpcL-PBS to RCs to manage the energy distribution for photosynthesis under dynamic flux conditions.

59 BASIC BIOLOGICAL SCIENCES↗

A Lake Biogeochemistry Model for Global Methane Emissions: Model Development, Site‐Level Validation, and Global Applicability

Abstract Lakes are important sentinels of climate change and may contribute over 30% of natural methane (CH 4 ) emissions; however, no earth system model (ESM) has represented lake CH 4 dynamics. To fill this gap, we refined a process‐based lake biogeochemical model to simulate global lake CH 4 emissions, including representation of lake bathymetry, oxic methane production (OMP), the effect of water level on ebullition, new non‐linear CH 4 oxidation kinetics, and the coupling of sediment carbon pools with in‐lake primary production and terrigenous carbon loadings. We compiled a lake CH 4 data set for model validation. The model shows promising performance in capturing the seasonal and inter‐annual variabilities of CH 4 emissions at 10 representative lakes for different lake types and the variations in mean annual CH 4 emissions among 106 lakes across the globe. The model reproduces the variations of the observed surface CH 4 diffusion and ebullition along the gradients of lake latitude, depth, and surface area. The results suggest that OMP could play an important role in surface CH 4 diffusion, and its relative importance is higher in less productive and/or deeper lakes. The model performance is improved for capturing CH 4 outgassing events in non‐floodplain lakes and the seasonal variability of CH 4 ebullition in floodplain lakes by representing the effect of water level on ebullition. The model can be integrated into ESMs to constrain global lake CH 4 emissions and climate‐CH 4 feedback.

54 ENVIRONMENTAL SCIENCES↗

Structural and genetic convergence of HIV-1 neutralizing antibodies in vaccinated non-human primates

A primary goal of HIV-1 vaccine development is the consistent elicitation of protective, neutralizing antibodies. While highly similar neutralizing antibodies (nAbs) have been isolated from multiple HIV-infected individuals, it is unclear whether vaccination can consistently elicit highly similar nAbs in genetically diverse primates. Here, we show in three outbred rhesus macaques that immunization with Env elicits a genotypically and phenotypically conserved nAb response. From these vaccinated macaques, we isolated four antibody lineages that had commonalities in immunoglobulin variable, diversity, and joining gene segment usage. Atomic-level structures of the antigen binding fragments of the two most similar antibodies showed nearly identical paratopes. The Env binding modes of each of the four vaccine-induced nAbs were distinct from previously known monoclonal HIV-1 neutralizing antibodies, but were nearly identical to each other. The similarities of these antibodies show that the immune system in outbred primates can respond to HIV-1 Env vaccination with a similar structural and genotypic solution for recognizing a particular neutralizing epitope. These results support rational vaccine design for HIV-1 that aims to reproducibly elicit, in genetically diverse primates, nAbs with specific paratope structures capable of binding conserved epitopes.

59 BASIC BIOLOGICAL SCIENCES↗

Earth System Model Overestimation of Cropland Temperatures Scales With Agricultural Intensity

Intensive crop growth can modify regional climate by partitioning energy to latent heating through transpiration, cooling growing season temperatures. Recent work shows that cooling associated with agriculture can dampen anthropogenic warming over breadbasket regions. However, it is unknown whether climate models reproduce crop influences on regional climate, and thus the future risk of extreme climate events over global breadbasket regions. Here, we show that models overestimate growing season temperatures and underestimate evapotranspiration (ET) over global croplands, and that these differences increase with cropped area. We trace this warm and dry difference through each model's representation of the surface energy budget, showing that model differences in transpiration, leaf area index, and the ratio of transpiration to total ET drive the overall effect. While the implications of these model deficiencies for future projections are uncertain, they point to the importance of improving representations of crop-climate processes to better assess breadbasket vulnerability to climate change.

54 ENVIRONMENTAL SCIENCES↗