Search NASA⌕ Search

SEARCH · Search NASA

Results for “default”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
NASA NTRS2022

Assessment of Retrieved GMI Emissivity Over Land, Snow, and Sea Ice in the GEOS System

Directly assimilating microwave radiances over land, snow and sea ice remains a significant challenge for data assimilation systems. These data assimilation systems are critical to the success of global numerical weather prediction systems including the Global Earth Observing System-Atmospheric Data Assimilation System (GEOS-ADAS). Extending more surface sensitive microwave channels over land, snow and ice could provide a needed source of data for Numerical Weather Prediction particularly in the Planetary Boundary Layer (PBL). Unfortunately, the accuracy of emissivity models currently available within the GEOS-ADAS along with other data assimilation systems are insufficient to simulate and assimilate radiances. Recently, Munchak et al. (2020) published a 5-year climatological database for retrieved microwave emissivity from the GPM Microwave Imager (GMI) aboard the Global Precipitation Mission (GPM). In this work the database is utilized by modifying the GEOS-ADAS to use this emissivity database in place of the default emissivity value available in the Community Radiative Transfer Model (CRTM), which is the fast radiative transfer model used by the GEOS-ADAS. As a first step, the GEOS-ADAS is run in a so-called “stand-alone” mode to simulate radiances from GMI using the default CRTM emissivity, and replacing the default CRTM emissivity models with values from Munchak et al, 2020. The simulated observations using Munchak et al., 2020 agree more closely with observations from GMI. These results are presented along with a discussion of the implication for GMI observations within the GEOS-ADAS.

GPM↗
DOE OSTI2022

Speed vs Accuracy: Effect on Ligand Pose Accuracy of Varying Box Size and Exhaustiveness in AutoDock Vina

Abstract Structure‐based virtual high‐throughput screening involves docking chemical libraries to targets of interest. A parameter pertinent to the accuracy of the resulting pose is the root mean square deviation (RMSD) from a known crystallographic structure, i. e., the ‘docking power’. Here, using a popular algorithm, Autodock Vina, as a model program, we evaluate the effects of varying two common docking parameters: the box size (the size of docking search space) and the exhaustiveness of the global search (the number of independent runs starting from random ligand conformations) on the RMSD from the PDBbind v2017 refined dataset of experimental protein‐ligand complexes. Although it is clear that exhaustiveness is an important parameter, there is wide variation in the values used, with variation between 1 and >100. We, therefore, evaluated a combination of cubic boxes of different sizes and five exhaustiveness values (1, 8, 25, 50, 75, 100) within the range of those commonly adopted. The results show that the default exhaustiveness value of 8 performs well overall for most box sizes. In contrast, for all box sizes, but particularly for large boxes, an exhaustiveness value of 1 led to significantly higher median RMSD (mRMSD) values. The docking power was slightly improved with an exhaustiveness of 25, but the mRMSD changes little with values higher than 25. Therefore, although low exhaustiveness is computationally faster, the results are more likely to be far from reality, and, conversely, values >25 led to little improvement at the expense of computational resources. Overall, we recommend users to use at least the default exhaustiveness value of 8 for virtual screening calculations.

Agarwal, Rupesh↗
DOE OSTI2022

Design and implementation of dynamic I/O control scheme for large scale distributed file systems

In this paper, we have analyzed the input/output (I/O) activities of Cori, which is a high-performance computing system at the National Energy Research Scientific Computing Center at Lawrence Berkeley National Laboratory. Our analysis results indicate that most users do not adjust storage configurations but rather use the default settings. In addition, owing to the interference from many applications running simultaneously, the performance varies based on the system status. To configure file systems autonomously in complex environments, we developed DCA-IO, a dynamic distributed file system configuration adjustment algorithm that utilizes the system log information to adjust storage configurations automatically. Our scheme aims to improve the application performance and avoid interference from other applications without user intervention. Moreover, DCA-IO uses the existing system logs and does not require code modifications, an additional library, or user intervention. To demonstrate the effectiveness of DCA-IO, we performed experiments using I/O kernels of real applications in both an isolated small-sized Lustre environment and Cori. Our experimental results shows that our scheme can improve the performance of HPC applications by up to 263% with the default Lustre configuration.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Linear and Nonlinear Solvers for Simulating Multiphase Flow within Large-Scale Engineered Subsurface Systems

Simulation of multiphase flow in the subsurface is well-known to be computationally challenging. While there have been many studies that have explored approaches to overcoming these challenges, they often utilize relatively simple case studies. In this paper, we focus on the unique numerical challenges posed by modeling large-scale engineered subsurface systems, characterized by discrete features embedded in a heterogeneous natural subsurface setting. The man-made features such as shafts, tunnels, and barriers often cause multiple challenges in modeling the domain for multiphase porous media flow. This flow scenario can have a wide range of applications such as nuclear waste repositories, enhanced recovery of a petroleum reservoir, geothermal engineering, and carbon sequestration. An example of these severe numerical challenges is the case of performance assessment (PA) for Waste Isolation Pilot Plant (WIPP), the only operating deep geological repository in the US, which simulates extreme material properties of bedded salt rock formation and extreme contrast due to open excavation next to the formation. The models have extremes not only of permeability and porosity but also of the constitutive models needed for multiphase flow; additionally, they have process models like salt creep closure reducing porosity over time, fracturing in clay and anhydrite interbeds of the bedded salt, gas generation from the waste materials, and unintentional human borehole intrusions in some scenarios. Numerical simulations require the solution of coupled systems of nonlinear PDEs; in our work, we use the open-source simulator PFLOTRAN which is based on Finite Volume discretization. The solution of the nonlinear equations requires use of the Newton-Raphson iteration at each time step, which entails the solution of the linearized Jacobian system at each iteration. The effects of all the processes (i.e., large number of unknowns, highly nonlinear constitutive relations, large contrasts in material properties in short distances) lead to an ill-conditioned Jacobian matrix that severely challenges traditional linear solver, i.e., stabilized biconjugate gradient with block Jacobi incomplete LU preconditioner (BCGS-ILU) leading to non-convergence for traditional Newton-Raphson nonlinear solver causing unacceptably long computation time for each model. This paper presents linear solvers such as constrained pressure residual (CPR) two-stage preconditioner with alternate-block-factorization (ABF) and quasi- implicit pressure and explicit saturation (QIMPES) decouplers and flexible generalized residual solver (FGMRES). The new general-purpose nonlinear solver, Newton trust-region dogleg Cauchy (NTRDC), is also introduced to resolve extreme nonlinearities in the models. We demonstrate the effectiveness of each method relative to the default BCGS-Newton solver. The two best cases had nearly 50 times speed-up and achieved completion of a simulation in 14 hours that never completed due to non-convergence with the default solver. We also investigate the strong scalability of each method and discuss some of the deficiencies found for Block Jacobi preconditioner using parallel domain decomposition, and node packing effects of modern processor architecture.

Preconditioner, Nonlinear, Porous media, Multiphas↗
DOE OSTI2025

Data-Driven Recommendation of Optimal Tuning Scheme for Range-Separated Hybrid Functionals in Solution-Phase UV/Vis Absorption Energy Prediction

Time-dependent density functional theory (TDDFT) combined with range-separated hybrid (RSH) functionals and a tuned range-separation parameter γ offers a computationally economical approach for high-throughput excited- state property predictions. The γ-tuning procedure in the gas phase is well established. However, no agreement on the best γ- tuning procedure has been made when considering the solvent effect with implicit solvent models like the polarizable continuum model (PCM). To answer that question, this study created a diverse dataset with 937 molecules with experimental solutionphase UV/vis absorption spectra. Three γ-tuning methods, the gasphase γ-tuning (GPγT), the partial vertical γ-tuning (PVγT), and the strict vertical γ-tuning (SVγT), were evaluated for the ωPBEh functional over the entire dataset. Additional benchmarks are done for the optimally tuned screened range-separated hybrid combined with the PCM approach (SRSH-PCM) and the solvation-mediated tuning procedure (sol-med-OT). Our findings revealed that the optimal γ-values obtained by the PVγT and the SVγT are significantly smaller than the GPγT. This trend holds consistently across all molecules in our dataset, and we explained the origin of this phenomenon. TDDFT calculations with PVγTand SVγT-tuned γ-values and default global Fock exchange fraction achieve superior performance compared to those using GPγTtuned or default γ and slightly outperform SRSH-PCM and sol-med-OT with similar or lesser computational cost. Furthermore, we found that the smaller γ-values from SVγT captured the expected 1/(εR) asymptotic behavior in the solution phase, resulting in accurate prediction of solution-phase CT excitations, consistent with the screened asymptote behavior encoded in SRSH-PCM. These results show that SVγT is the best scheme for high-throughput UV/vis absorption spectrum calculations using the ωPBEh functional from a data-driven perspective.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Improved Convective Ice Microphysics Parameterization in the NCAR CAM Model

Partitioning deep convective cloud condensates into components that sediment and detrain, known to be a challenge for global climate models, is important for cloud vertical distribution and anvil cloud formation. In this study, we address this issue by improving the convective microphysics scheme in the National Center for Atmospheric Research Community Atmosphere Model version 5.3 (CAM5.3). The improvements include: (1) considering sedimentation for cloud ice crystals that do not fall in the original scheme, (2) applying a new terminal velocity parameterization that depends on the environmental conditions for convective snow, (3) adding a new hydrometeor category, “rimed ice,” to the original four-class (cloud liquid, cloud ice, rain, and snow) scheme, and (4) allowing convective clouds to detrain snow particles into stratiform clouds. Results from the default and modified CAM5.3 models were evaluated against observations from the U.S. Department of Energy Tropical Warm Pool-International Cloud Experiment (TWP-ICE) field campaign. The default model overestimates ice amount, which is largely attributed to the underestimation of convective ice particle sedimentation. By considering cloud ice sedimentation and rimed ice particles and applying a new convective snow terminal velocity parameterization, the vertical distribution of ice amount is much improved in the midtroposphere and upper troposphere when compared to observations. The vertical distribution of ice condensate also agrees well with observational best estimates upon considering snow detrainment. Further, comparison with observed convective updrafts reveals that current bulk model fails to reproduce the observed updraft magnitude and occurrence frequency, suggesting spectral distributions be required to simulate the subgrid updraft heterogeneity.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Impacts of Sub‐Grid Topographic Representations on Surface Energy Balance and Boundary Conditions in the E3SM Land Model: A Case Study in Sierra Nevada

Abstract Sub‐grid topographic heterogeneity has large impacts on surface energy balance and land‐atmosphere interactions. However, the impacts of representing sub‐grid topographic effects in land surface models (LSMs) on surface energy balance and boundary conditions remain unclear. This study analyzed and evaluated the impacts of sub‐grid topographic representations on surface energy balance, turbulent heat flux, and scalar (co‐)variances in the Energy Exascale Earth System Model (E3SM) land model (ELM). Three sub‐grid topographic representations in ELM were compared: (a) the default sub‐grid structure (D), (b) the recently developed sub‐grid topographic structure (T), and (c) high spatial resolution (1KM). Additionally, two different solar radiation schemes in ELM were compared: (a) the default plane‐parallel radiative transfer scheme (PP) and (b) the parameterization scheme (TOP) that accounts for sub‐grid topographic effects on solar radiation. A series of offline simulations with the three grid discretization structures (D, T, and 1KM) and two schemes of solar radiation (TOP and PP) were carried out using the Sierra Nevada, California. 1KM simulations with TOP well capture the spatial heterogeneity of surface fluxes compared to Moderate Resolution Imaging Spectroradiometer remote sensing data. There are significant differences between TOP and PP in the 1‐km simulated surface energy balance, but the differences in mean values and standard deviations become small when aggregated to the grid scale (i.e., 0.5°). The T configuration better mimics the 1KM simulations with TOP than the D configuration and better captures the sub‐grid topographic effects on surface energy balance and boundary conditions. These results underline the importance of representing sub‐grid topographic heterogeneities in LSMs and motivate future research to understand the sub‐grid topographic effects on land‐atmosphere interactions over mountainous areas.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2023

Orographic Cirrus and Its Radiative Forcing in NCAR CAM6

Cirrus clouds play an important role in the Earth's radiative energy budget, thereby affecting the climate state and climate change. Orographic gravity wave (OGW)-induced sub-grid scale vertical velocity (i.e., cooling rate) is not resolved by large-scale models and its impact on ice formation in cirrus clouds is not well quantified. In this study, one sub-grid scale OGW scheme (e.g., McFarlane) is used in the Community Atmosphere Model version 6 (CAM6) to generate vertical velocity variance (σw) for cirrus formation. Results from the default model and simulations with the OGW-induced σw are evaluated against the DOE ARM Small Particles in Cirrus (SPARTICUS) campaign observations. The OGW based on the McFarlane scheme increases the sub-grid scale σw over mountains compared to the default model and improves the model agreement with the SPARTICUS observations. Larger σ w due to OGWs can trigger more frequent homogeneous nucleation in orographic cirrus and generates a higher number concentration of ice crystals observed during the SPARTICUS campaign. Moreover, our evaluation of the model simulations against satellite observations indicates that the McFarlane scheme generates high in-cloud ice number concentrations (>200 L -1 ) in the upper troposphere over mountains and high plateaus at mid- and high-latitudes of the winter hemisphere as shown in the observations. More ice crystals with smaller sizes absorb more infrared radiation (+0.523 ± 0.125 W m -2 ). The net radiative cloud forcing change at the top of the atmosphere is +0.330 W m -2 due to the orographic cirrus.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

The Role of In-Cloud Wet Removal in Simulating Aerosol Vertical Profiles and Cloud Radiative Forcing

Among the physical processes controlling aerosol vertical profiles, in-cloud wet removal is of utmost importance while its representation in global climate models (GCMs) is crude. In this study, we implement into the Community Atmosphere Model version 6 (CAM6) a physically based aerosol wet removal parameterization scheme that explicitly treats aerosol activation, removal and resuspension. Evaluation against in-situ observations shows that the default scheme substantially overestimates the upper tropospheric black carbon (BC) and sea salt mass. Our physically based scheme reduces BC and sea salt mass by a factor of 10 and 1,000, respectively, in better agreement with observations. Also, the new scheme slightly increases number of aerosol particles between 12 nm and 4.8 μm in diameter, thereby mitigating the aerosol number underestimation in the default scheme. Our new scheme reduces the overestimation of coarse-mode aerosol (0.5–4.8 μm) number. Overall, the aerosol property changes (mass decrease and number increase) reduce the cloud condensation nuclei (CCN) concentration at low supersaturation (i.e., 0.02% and 0.1%), and increase CCN at high supersaturations (i.e., 0.5% and 1%). Consequently, the global annual mean cloud liquid water path increases by 1.89 g m -2 and the ice water path increases by 0.51 g m -2 . The global annual mean shortwave, longwave, and net cloud radiative forcing change by -1.06, 0.57, and -0.48 W m -2 , respectively. Further improvement is needed to reflect the real physics that the removal efficiencies for aerosol mass and number are disproportionate and to advect cloud-borne (activated) aerosols for a complete aerosol lifecycle.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Improving BC Mixing State and CCN Activity Representation With Machine Learning in the Community Atmosphere Model Version 6 (CAM6)

Abstract Representing mixing state of black carbon (BC) is challenging for global climate models (GCMs). The Community Atmosphere Model version 6 (CAM6) with the four‐mode version of the Modal Aerosol Module (MAM4) represents aerosols as fully internal mixtures with uniform composition within each aerosol mode, resulting in high degree of internal mixing of BC with non‐BC species and large mass ratio of coating to BC ( R BC , the mass ratio of non‐BC species to BC in BC‐containing particles). To improve BC mixing state representation, we coupled a machine learning (ML) model of BC mixing state index trained on particle‐resolved simulations to the CAM6 with MAM4 (MAM4‐ML). In MAM4‐ML, we use R BC to partition accumulation mode particles into two new modes, BC‐free particles and BC‐containing particles. We adjust R BC to make the modeled BC mixing state index ( χ mode ) match the one predicted by the ML model ( χ ML ). On a global average, the mass fraction of BC‐containing particles in accumulation mode decreases from 100% (MAM4‐default) to 48% (MAM4‐ML). The globally averaged χ mode decreases from 78% (MAM4‐default) to 63% (MAM4‐ML, 19% reduction) and agrees well with χ ML (66%). The R BC decreases by 52% for accumulation mode and better agrees with observations. The hygroscopicity drops by 9% for BC‐containing particles in accumulation mode, leading to a 20% reduction in the BC activation fraction. The surface BC concentration increases most (6.9%) in the Arctic, and the BC burden increases by 4%, globally. Our study highlights the application of the ML model for improving key aerosol processes in GCMs.

Meteorology & Atmospheric Sciences↗
DOE OSTI2025

Integrating Characteristic Arctic Vegetation in a Land Surface Model Improves Representation of Carbon Dynamics Across a Tundra Landscape

Arctic warming is altering vegetation and carbon dynamics with global implications, yet Earth System Model (ESM) predictions in the Arctic remain highly uncertain, in part due to historically limited data for model parameterization and validation. As such, ESMs typically represent Arctic ecosystems in an oversimplified manner. Recently, nine plant functional types (PFTs) designed to realistically represent tundra vegetation were integrated into the Energy Exascale Earth System Model (E3SM) Land Model (ELM) and parameterized using plot-scale observations from a single site. Additional evaluation was needed to determine their transferability across the Arctic. Here, in this study, we evaluated whether refined representation of tundra vegetation improved model accuracy by conducting spatially explicit 100 × 100 m resolution ELM simulations on Alaska's Seward Peninsula. Simulations with the default two-PFT configuration and with the nine Arctic-specific PFTs were benchmarked against observations of net ecosystem exchange, gross primary production, and aboveground biomass from multiple data streams including an eddy covariance flux tower, flux chambers, and aircraft and unoccupied aerial system hyperspectral remote sensing. Evaluation revealed that Arctic-specific PFT simulations produced more realistic landscape-level carbon exchanges, and better captured observed heterogeneity in biomass and productivity, explaining 60%–70% of spatial variance (R 2 = 0.6–0.7) compared to just 12%–18% (R 2 = 0.12–0.18) with the default configuration. However, the refined model failed to reproduce observed aboveground biomass for highly productive alder-willow communities, requiring further evaluation of carbon allocation parameterizations for tall shrubs that are increasingly expanding across tundra landscapes. Our results demonstrate that enhanced representation of vegetation heterogeneity boosts predictive understanding of tundra carbon dynamics, facilitating regional to pan-Arctic model and remote-sensing scaling.

Murphy, Bailey A. [Oak Ridge National Laboratory (↗
DOE OSTI2023

Physics informed deep neural network embedded in a chemical transport model for the Amazon rainforest

Secondary organic aerosols (SOA) are fine particles in the atmosphere, which interact with clouds, radiation and affect the Earth’s energy budget. SOA formation involves chemistry in gas phase, aqueous aerosols, and clouds. Simulating these chemical processes involve solving a stiff set of differential equations, which are computationally expensive steps for three-dimensional chemical transport models. Deep neural networks (DNNs) are universal function approximators that could be used to represent the complex nonlinear changes in aerosol physical and chemical processes; however, key challenges such as generalizability to extended time periods, preservation of mass balance, simulating sparse model outputs, and maintaining physical constraints have limited their use in atmospheric chemistry. Here, we develop an approach of using a physics-informed DNN that overcomes previous such challenges and demonstrates its applicability for the chemical formation processes of isoprene epoxydiol SOA (IEPOX-SOA) over the Amazon rainforest. The DNN is trained with data generated by simulating IEPOX-SOA over the entire atmospheric column, using the Weather Research and Forecasting Model coupled with Chemistry (WRF-Chem). The trained DNN is then embedded within WRF-Chem to replace the computationally expensive default solver of IEPOX-SOA formation. The trained DNN predictions generalizes well with the default model simulation of the IEPOX-SOA mass concentrations and its size distribution (20 size bins) over several days of simulations in both dry and wet seasons. The embedded DNN reduces the computational expense of WRF-Chem by a factor of 2. Our approach shows promise in terms of application to other computationally expensive chemistry solvers in climate models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Performance Evaluation of Adaptive Routing on Dragonfly-based Production Systems

Performance of applications in production environments can he sensitive to network congestion. Cray Aries supports adaptively routing each network packet independently based on the load or congestion encountered as a packet traverses the network. Software can dictate different routing policies, adjusting between minimal and non-minimal bias, for each posted message. We have extensively evaluated the sensitivity of the routing bias selection on application performance as well as whole system performance in both production and controlled conditions. We show that the default routing bias used in Aries-based systems is often sub-optimal and that using a higher bias towards minimal routes will not only reduce the congestion effects on the application but also will decrease the overall congestion on the network. This routing scheme results in not only improved mean performance (by up to 12%) of most production applications hut also reduced run-to-run variability. Our study prompted the two supercomputing facilities (ALCF and NERSC) to change the default routing mode on their Aries-based systems. We present the substantial improvement measured in the overall congestion management and interconnect performance in production after making this change.

Chunduri, Sudheer↗
DOE OSTI2022

Multi-modelling predictions show high uncertainty of required carbon input changes to reach a 4‰ target

Soils store vast amounts of carbon (C) on land, and increasing soil organic carbon (SOC) stocks in already managed soils such as croplands may be one way to remove C from the atmosphere, thereby limiting subsequent warming. The main objective of this study was to estimate the amount of additional C input needed to annually increase SOC stocks by 4‰ at 16 long-term agricultural experiments in Europe, including exogenous organic matter (EOM) additions. We used an ensemble of six SOC models and ran them under two configurations: (1) with default parametrization and (2) with parameters calibrated site-by-site to fit the evolution of SOC stocks in the control treatments (without EOM). We compared model simulations and analysed the factors generating variability across models. The calibrated ensemble was able to reproduce the SOC stock evolution in the unfertilised control treatments. We found that, on average, the experimental sites needed an additional 1.5 ± 1.2 Mg C ha -1 year -1 to increase SOC stocks by 4‰ per year over 30 years, compared to the C input in the control treatments (multi-model median ± median standard deviation across sites). That is, a 119% increase compared to the control. While mean annual temperature, initial SOC stocks and initial C input had a significant effect on the variability of the predicted C input in the default configuration (i.e., the relative standard deviation of the predicted C input from the mean), only water-related variables (i.e., mean annual precipitation and potential evapotranspiration) explained the divergence between models when calibrated. Here, our work highlights the challenge of increasing SOC stocks in agriculture and accentuates the need to increasingly lean on multi-model ensembles when predicting SOC stock trends and related processes. To increase the reliability of SOC models under future climate change, we suggest model developers to better constrain the effect of water-related variables on SOC decomposition.

4 per 1000 initiative↗
DOE OSTI2024

CLDERA-E3SM

SAND2025-03846O CLDERA-E3SM enhances the E3SMv2 model by adding new options such as complete stratospheric volcanic aerosol representation and idealized physics model configuration for the atmosphere component. CLDERA-E3SM also allows for aerosol source emission specification, which is a powerful method for tracking aerosol and its climate impacts in E3SMv2. This new code better represents stratospheric sulfate aerosol and includes modifications such as changes to the default aerosol size distributions and to stratospheric aerosol microphysical treatment, which allows stratospheric sulfate aerosol to grow larger than in the default E3SMv2. New aerosol source tagging can be customized to separate the global emission of select aerosol tracers—such as sulfate, black carbon, and organic carbon aerosol—into a new, user specified list of sources. It also includes idealized simulations of the atmosphere and volcanic aerosol tracers. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.

Hillman, Benjamin↗
DOE OSTI2025

Commercial Fleet Level Emissions and Energy Tracker (COFLEET) v1.0

This tool generates the latest US commercial vehicle fleet composition and operational characteristics using 2021 US VIUS data and assess the fleetwide energy and emission outcomes. The emission rates for commercial trucks and default fleet composition are collected from the U.S. EPA's MOtor Vehicle Emission Simulator (MOVES4). The 2021 US VIUS data is applied to generate fleet characteristics such as the long-haul fraction and the vehicle mileage accumulation rate. The tool also provides fleet turnover and emission forecasts under various forecasted fleet electrification scenarios. This tool helps policymakers and practitioners advance the commercial fleet generation for emission models. This study also provides some sample datasets for state and local transportation/air quality agencies to test, which can reduce the estimation bias associated with using MOVES default fleets.

Xu, Xiaodan [Lawrence Berkeley National Laboratory↗
DOE OSTI2026

matsim-agents v1.0

matsim-agents is a multi-agent AI framework for atomistic materials simulation and discovery. It orchestrates large language models (LLMs), machine-learned interatomic potentials (MLIPs), and DFT codes into a single agentic loop running on laptops and DOE leadership-class supercomputers. MULTI-AGENT ORCHESTRATION A LangGraph state machine with three nodes: a Planner that converts a natural-language research objective into structured tasks; an Executor that dispatches atomistic tools and loops until the queue is empty; and an Analyst that summarizes results into a human-readable report. State is checkpointed after every step and human-in-the-loop gates can be inserted at any edge. HYPOTHESIS-DRIVEN DISCOVERY CHAT An interactive REPL (matsim-agents chat) that couples LLM dialogue with atomistic simulation. Chemical formulas are automatically detected in conversation turns and trigger a full crystal-phase exploration: structure generation → relaxation → stability scoring → result injection back into the conversation, creating a closed hypothesis-refinement loop. CRYSTAL PHASE ENUMERATION Given a composition, the phase explorer enumerates prototypes by stoichiometry: elemental (fcc/bcc/hcp/sc/diamond), binary 1:1 (rocksalt/CsCl/zincblende/ wurtzite/fluorite/rutile), ternary 1:1:3 (cubic perovskite), ternary 1:2:4 (perovskite + spinel), quaternary 1:1:2:6 (Fm-3m double perovskite). 2-D prototypes (graphene, h-BN, MoS2 2H/1T) and multilayer stacking are also supported via --include-2d and --num-layers. SUPERCELL GENERATION AND SITE DECORATION Auto-tiling to a minimum atom count (--min-atoms), explicit NxNxN tiling (--supercell), symmetry-distinct site decorations (--n-orderings), and isotropic lattice-scale sweeps (--lattice-scales) for volume bracketing. MLFF RELAXATION AND STABILITY SCORING HydraGNN (multi-headed GNN) drives structure relaxation via ASE with FIRE, BFGS, or BFGSLineSearch. Stability output: delta-E/atom ranking across phases and a max-residual-force dynamical-stability proxy. Other MLIPs (MACE, NequIP, Orb) can be plugged in through the same interface. DFT BACKENDS Quantum ESPRESSO pw.x and VASP 6.6 are first-class labellers. Both have validated GPU builds and SLURM/PBS launchers for three DOE platforms: Frontier (AMD MI250X, ROCm), Aurora (Intel PVC, oneAPI), Perlmutter (NVIDIA A100, CUDA). QE produces ~100 binaries (pw.x, ph.x, epw.x, ...). VASP supports scf, relax, vc-relax, and vc-relax-shape run types. ACTIVE-LEARNING LOOP matsim-agents al run CONFIG.yaml drives an iterative HydraGNN-DFT loop: MD generates candidates → ensemble/MC-dropout uncertainty selects the most informative → DFT labels them in parallel inside one allocation → dataset grows → HydraGNN retrains → repeat. DFT backend is a single YAML toggle (dft.backend: vasp | qe). LLM-generated seed structures are supported (no curated POSCAR library needed). Config uses ${VAR}, ${VAR:-default}, ${VAR:?msg} shell-style substitution for cross-user/cross-site portability. LLM BACKENDS Ollama (local, default), vLLM (HPC multi-GPU serving), OpenAI, Anthropic, HuggingFace Transformers+Accelerate. Selected at runtime via flag or env var with no code changes. HPC PORTABILITY Same Python entry points run on Frontier (ROCm 7.2), Aurora (oneAPI), and Perlmutter (CUDA 12). DFT and ML stacks are never co-loaded in the same shell; they couple through the scheduler and filesystem. Advanced multi-node launchers (serve, discovery-chat, single-relaxation, active-learning, QE warm-start) are provided for all three platforms. CODABENCH COMPETITION BUNDLE A self-contained benchmark: 159 atomistic test structures across 11 material classes, 5 tasks (formation energy, forces, ML relaxation, AI-DFT relaxation, phase stability ranking), public/private leaderboard split (30/70), and four ready-to-run baselines: MACE-MP-0, HydraGNN, UMA, AllScAIP.

Lupo Pasini, Massimiliano [Oak Ridge National Labo↗
DOE OSTI2022

Conservation of Dry Air, Water, and Energy in CAM and Its Potential Impact on Tropical Rainfall

Abstract For the Community Atmosphere Model version 6 (CAM6), an adjustment is needed to conserve dry air mass. This adjustment exposes an inconsistency in how CAM6’s energy budget incorporates water—in CAM6 water in the vapor phase has energy, but condensed phases of water do not. When water vapor condenses, only its latent energy is retained in the model, while its remaining internal, potential, and kinetic energy are lost. A global fixer is used in the default CAM6 model to maintain global energy conservation, but locally the energy tendency associated with water changing phase violates the divergence theorem. This error in energy tendency is intrinsically tied to the water vapor tendency, and reaches its highest values in regions of heavy rainfall, where the error can be as high as 40 W m −2 annually averaged. Several possible changes are outlined within this manuscript that would allow CAM6 to satisfy the divergence theorem locally. These fall into one of two categories: 1) modifying the surface flux to balance the local atmospheric energy tendency and 2) modifying the local atmospheric tendency to balance the surface plus top-of-atmosphere energy fluxes. To gauge which aspects of the simulated climate are most sensitive to this error, the simplest possible change—where condensed water still does not carry energy and a local energy fixer is used in place of the global one—is implemented within CAM6. Comparing this experiment with the default configuration of CAM6 reveals precipitation, particularly its variability, to be highly sensitive to the energy budget formulation. Significance Statement This study examines and explains spurious regional sources and sinks of energy in a widely used climate model. These energy errors result from not tracking energy associated with water after it transitions from the vapor phase to either liquid or ice. Instead, the model used a global fixer to offset the energy tendency related to the energy sources and sinks associated with condensed water species. We replace this global fixer with a local one to examine the model sensitivity to the regional energy error and find a large sensitivity in the simulated hydrologic cycle. This work suggests that the underlying thermodynamic assumptions in the model should be revisited to build confidence in the model-simulated regional-scale water and energy cycles.

Harrop, Bryce E.↗