Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Model assignment”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

A Vehicle-to-Grid planning framework incorporating electric vehicle user equilibrium and distribution network flexibility enhancement

The rapid surge in electric vehicle (EV) adoption, coupled with advancements in charging technologies, emphasizes the critical necessity for expanding EV recharging infrastructure. Simultaneously, the Distribution Network (DN) encounters escalating challenges in meeting charging demand during peak traffic periods. Consequently, there is a mounting demand for the deployment of innovative Vehicle-to-Grid (V2G) technologies to augment the DN’s flexibility in power dispatch and alleviate travel costs for EV users. Hence, this paper proposes an EV-user-equilibrium-(UE)-constrained V2G planning framework that enhances flexibility in the DN. The framework aims to ascertain the optimal placement and capacity of EV charging stations (EVCSs) and V2G charging piles within the Transportation Network (TN). It takes into account the equilibrium condition stemming from competitive EV charging and routing behaviors alongside the optimal expansion of DN energy resources to accommodate the electricity supplied by the V2G piles. This study commences by analyzing EV drivers’ travel decisions, considering the influence of charging and V2G pile locations and sizes. Subsequently, we tackle the Traffic Assignment Problem with User Equilibrium (TAP-UE) model to characterize the steady-state traffic flow distribution of EVs. Following this, we formulate the optimization model for the Coordinated Power and Transportation Network (CPTN), which encompasses the optimal expansion of DN facilities and traffic flow regulation under UE conditions. To mitigate the computational complexity associated with the V2G planning model, we introduce a series of linearization methods to obtain a manageable Mixed-Integer Linear Programming (MILP) solution. Finally, to validate the efficacy of our proposed planning framework, we apply it to two test systems, including a real-world case study. Through these case studies, we explore the necessity and potential benefits of V2G technologies.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2025

The three cellulose synthase isoforms for secondary cell wall make specific contributions to microfibril synthesis

Cellulose is synthesized at the plasma membrane by the cellulose synthase complex, a structure that contains three distinct isoforms of the catalytic subunit, cellulose synthase A (CESA). The division into three subunits appears early in land plant evolution and is highly conserved, particularly for the secondary cell wall. However, what if any unique roles each isoform plays in the complex remain unclear. Here, we assessed the contributions of specific isoforms to microfibril synthesis. First, we expressed CESA isoforms of the primary cell wall or the moss Physcomitrium patens in Arabidopsis thaliana backgrounds missing a secondary cell wall CESA. While the primary cell wall isoforms rescued the cesa knockout phenotype with partial isoform specificity, those from the moss rescued with fewer restrictions. Then, we recreated various CESA missense mutations in all three of the secondary cell wall isoforms; while results are consistent with isoform specificity, they are difficult to interpret further without molecular structures. Finally, we show that catalytically inactive CESA isoforms restore growth and cellulose content in the corresponding knockout in an isoform-specific manner; along with partial rescue of the growth and cellulose content of the inflorescence stem, the replacement lines have fiber cells with partially disorganized microfibrils and secondary cell wall cellulose with narrow crystal width. Generally, effects were more pronounced in lines where CESA8 was inactivated compared with inactivating CESA4 or 7, which tended to have similar phenotypes to each other. Here, we account for these results with a model for cellulose synthase structure with the isoforms assigned specific localization within the cellulose synthase complex.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

CrossLink: Meshing a 3D part from a STEP file [Slides]

This example shows how to use CrossLink to create a mesh for a 3D geometry part read in from a STEP file. A STEP file (Standard for the Exchange of Product data) is a common file format used for 3D modelling that can be written out from a CAD program (e.g., Creo Parametric). Being able to read and mesh STEP file geometries is essential for meshing parts from engineering models. Creating the mesh typically involves: (1) Importing the geometry; (2) Creating geometry groups; (3) Assigning curves and surfaces to the geometry groups; (4) Building a topology; (5) Applying geometric and mesh constraints; and (6) Generating the final mesh. Current issues with this process in CrossLink include: (1) The GUI does not display trimmed surfaces; (2) The user must manually create and assign geometry groups; and (3) The user must be aware of duplicate curves and surfaces from the CAD model.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Grid-Aware Charging and Operational Optimization for Mixed-Fleet Public Transit

The rapid growth of urban populations and the increasing need for sustainable transportation solutions have prompted a shift towards electric buses in public transit systems. However, the effective management of mixed fleets consisting of both electric and diesel buses poses significant operational challenges. One major challenge is coping with dynamic electricity pricing, where charging costs vary throughout the day. Transit agencies must optimize charging assignments in response to such dynamism while accounting for secondary considerations such as seating constraints. This paper presents a comprehensive mixed-integer linear programming (MILP) model to address these challenges by jointly optimizing charging schedules and trip assignments for mixed (electric and diesel bus) fleets while considering factors such as dynamic electricity pricing, vehicle capacity, and route constraints. We address the potential computational intractability of the MILP formulation, which can arise even with relatively small fleets, by employing a hierarchical approach tailored to the fleet composition. By using real-world data from the city of Chattanooga, Tennessee, USA, we show that our approach can result in significant savings in the operating costs of the mixed transit fleets.

Sen, Rishav↗
DOE OSTI2021

MOOSE Framework Meshing Enhancements to Support Reactor Analysis

MOOSE-based physics codes require an input finite element mesh on which the physics solution is calculated, reported, and transferred to other physics codes. The use of difficult-touse, external licensed software is often required to generate high quality meshes for reactor geometries. High-fidelity geometry modeling also requires elaborate tracking of groups of elements for material property assignment and output reporting which can be considerably complex for the user to identify and maintain. Under the U.S. Department of Energy Office of Nuclear Energy Advanced Modeling and Simulation (NEAMS) program, several meshingrelated enhancements have been developed for the MOOSE framework to address user challenges in creating finite element meshes for advanced reactor geometries. MOOSE mesh generators have been developed to mesh hexagonal geometries (pins, ducted assemblies, and cores) commonly found in liquid-metal cooled fast reactor concepts. The mesh generator used for hexagonal pin cells is generic for regular polygons and therefore may also be used for Cartesian pin cells. Hexagonal pin cells can be stitched into ducted assemblies, and assemblies can be stitched together into a core. The user may specify region ids, region names, and other preferences on the mesh. This control is useful for later material mapping in the MOOSE-based physics codes input. A capability was also developed for meshing rotating control drums including determination of material volume fractions in each mesh element as a function of time. Control drum meshes may be stitched to other hexagonal assemblies to create a core configuration. Additional mesh generators were developed that wrap around the hexagonal meshing capabilities and utilize “extra element integer” ID values on each element. In regular Cartesian or hexagonal assemblies or cores, the bookkeeping of element groups for both material assignment and output reporting can now be automated through assignment of pin, assembly, core, axial and depletion id values stored as extra element integers. The extra element tags on the mesh greatly speed the reactor analyst’s efforts to map materials to meshes, track depletion zones, and parse output such as axial pin power distributions. At the highest level, pin, assembly, and core mesh generators (with this reactor terminology) have also been developed to easily generate regular Cartesian and hexagonal cores, including axial extrusion. These reactor geometry builders call upon the previously mentioned capabilities to produce analysis-ready 3D meshes including material assignments. Open source mesh triangulation capabilities were also investigated for integration into the MOOSE framework to address the need for meshing the core periphery region which extends from the irregular outer assembly border to a cylindrical boundary. Options are limited due to licensing constraints, and the recommendation is pursue building a native MOOSE Delaunay triangulator routine with full functionality. Finally, a series of verification problems were performed with NEAMS physics tools. All developed capabilities will be available in the new open-source “Reactor” module of the MOOSE framework, which is accessible to any MOOSE-based NEAMS physics tool.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Estimating the Likelihood of GHG Concentration Scenarios From Probabilistic Integrated Assessment Model Simulations

The climate scenarios that form the basis for current climate risk assessments have no assigned probabilities, and this impedes the analysis of future climate risks. This paper proposes an approach to estimate the probability of carbon dioxide (CO 2 ) concentration scenarios used in key climate change modeling experiments. It computes the CO 2 emissions compatible with the concentrations prescribed by Coupled Model Intercomparison Project Phase 5 (CMIP5) and CMIP6 experiments. The distribution of these compatible cumulative emissions is interpreted as the likelihood of future emissions given a concentration pathway. Using Bayesian analysis, the probability of each pathway can be estimated from a probabilistic sample of future emissions. The approach is demonstrated with five probabilistic CO 2 emission simulation ensembles from four Integrated Assessment Models (IAM), leading to independent estimates of the likelihood of the CO 2 concentration of Representative Concentration Pathways (RCP) and Shared Socioeconomic Pathways (SSP). Results suggest that SSP5-8.5 is unlikely for the second half of the 21st century, but offer no clear consensus on which of the remaining scenarios is most likely. Estimates of likelihoods of CO 2 concentrations associated with RCP and SSP scenarios are affected by sampling errors, differences in emission sources simulated by the IAMs, and a lack of a common experimental framework for IAM simulations. These shortcomings, along with a small IAM ensemble size, limit the applicability of the results presented here. Novel joint IAM and the Earth System Model experiments are needed to deliver actionable probabilistic climate risk assessments.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Emergence of nosocomial associated opportunistic pathogens in the gut microbiome after antibiotic treatment

Introduction: According to the Centers for Disease Control’s 2015 Hospital Acquired Infection Hospital Prevalence Survey, 1 in 31 hospital patients was infected with at least one nosocomial pathogen while being treated for unrelated issues. Many studies associate antibiotic administration with nosocomial infection occurrence. However, to our knowledge, there is little to no direct evidence of antibiotic administration selecting for nosocomial opportunistic pathogens. Aim: This study aims to confirm gut microbiota shifts in an animal model of antibiotic treatment to determine whether antibiotic use favors pathogenic bacteria. Methodology: We utilized next-generation sequencing and in-house metagenomic assembly and taxonomic assignment pipelines on the fecal microbiota of a urinary tract infection mouse model with and without antibiotic treatment. Results: Antibiotic therapy decreased the number of detectable species of bacteria by at least 20-fold. Furthermore, the gut microbiota of antibiotic treated mice had a significant increase of opportunistic pathogens that have been implicated in nosocomial infections, like Acinetobacter calcoaceticus/baumannii complex, Chlamydia abortus, Bacteroides fragilis, and Bacteroides thetaiotaomicron. Moreover, antibiotic treatment selected for antibiotic resistant gene enriched subpopulations for many of these opportunistic pathogens. Conclusions: Oral antibiotic therapy may select for common opportunistic pathogens responsible for nosocomial infections. In this study opportunistic pathogens present after antibiotic therapy harbored more antibiotic resistant genes than populations of opportunistic pathogens before treatment. Our results demonstrate the effects of antibiotic therapy on induced dysbiosis and expansion of opportunistic pathogen populations and antibiotic resistant subpopulations of those pathogens. Follow-up studies with larger samples sizes and potentially controlled clinical investigations should be performed to confirm our findings.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Mascon distribution techniques for asteroids and comets

The mass-concentration model is an approach that has been used to model the gravitational fields of irregularly shaped bodies such as asteroids and comets. By this approach, the body is treated as a collection of point masses. The method is conceptually simple, easy to program, valid down to the surface, and capable of modeling arbitrary density heterogeneities. How the mass concentrations are distributed as well as how mass is assigned to these concentrations is, however, nontrivial. These aspects significantly affect the accuracy and efficiency of the gravitational model. In this paper, we frame the distribution process in terms of numerical integration applied to finite volume meshes. We describe a new method using unstructured, curvilinear, finite volume meshes to significantly improve the accuracy of the mass-concentration model. We then compare the accuracy and efficiency of several variations of our distribution technique to those from literature using Asteroid Eros and Bennu as example bodies. Our results show that the mascon model can be as accurate as the analytic polyhedral model at the surface using an equivalent number of computational elements—i.e., mascon to surface facets. We report the improvement in the model’s performance can be mainly attributed to the volume mesh topology while mesh curving can provide modest case-dependent improvements.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Distinctive signals of frustrated dark matter

We study a renormalizable model of Dirac fermion dark matter (DM) that communicates with the Standard Model (SM) through a pair of mediators — one scalar, one fermion — in the representation (6, 1, $\frac{4}{3}$) of the SM gauge group SU(3) c x SU(2) L x U(1) Y . While such assignments preclude direct coupling of the dark matter to the Standard Model at tree level, we examine the many effective operators generated at one-loop order when the mediators are heavy, and find that they are often phenomenologically relevant. We reinterpret dijet and pair-produced resonance and jets + E$^{miss}_{T}$ searches at the Large Hadron Collider (LHC) in order to constrain the mediator sector, and we examine an array of DM constraints ranging from the observed relic density Ω χ h$^{2}_{Planck}$ to indirect and direct searches for dark matter. Tree-level annihilation, available for DM masses starting at the TeV scale, is required in order to produce Ω χ h$^{2}_{Planck}$ through freeze-out, but loops — led by the dimension-five DM magnetic dipole moment — are nonetheless able to produce signals large enough to be constrained, particularly by the XENON1T experiment. In some benchmarks, we find a fair amount of parameter space left open by experiment and compatible with freeze-out. In other scenarios, however, the open space is quite small, suggesting a need for further model-building and/or non-standard cosmologies.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Machine learning based rate optimization under geologic uncertainty

We propose a novel approach for rate optimization during a waterflood under geologic uncertainty in reservoir properties such as permeability and porosity. The traditional approach typically involves several runs of the forward simulator. This may not scale well when the optimization is to be performed at the full field-level and over multiple geologic realizations. A machine-learning (ML) based approach which is quick and scalable for rate optimization over multiple geologic realizations is proposed instead. The training data for the model is generated by running the forward simulator with randomly assigned well rates using multiple geologic realizations. A reduced order representation of the permeability heterogeneity in each of the realizations is derived using a grid connectivity transformation (GCT). This step involves finding basis functions corresponding to the different modal frequencies of the grid connectivity represented by the grid Laplacian. The projection of the heterogeneous property field along these basis functions gives the basis coefficients that form the reduced order representation. Subsequently, for each training datapoint, streamlines are traced and the minimum time of flight (TOF) representing the tracer breakthrough time at each producer is recorded. The basis coefficients and well rates are fed to a machine learning model as input and the minimum TOF at the producers forms the output of the model. This trained model can then be used along with an optimizer for computing the optimal injection rates to maximize the injection sweep efficiency. This corresponds to minimizing the variance in the minimum TOF within each well group. Different architectures of neural network are tested using 5-fold cross validation to decide the best ML model to compute the streamline time of flight. The trained model is used to perform well rate optimization over multiple realizations of geology by using a risk tolerance penalty. The optimal well rates thus obtained are compared with two cases: a) equal well rates assigned to all injectors and producers and b) well rates obtained by optimizing over a single realization without considering the uncertainty in geology. The optimal well rates are seen to offer better oil recovery and sweep efficiency than both cases.

02 PETROLEUM↗
DOE OSTI2022

Why is it still too warm or cold in my house? Examining the relationships between energy efficient capital and household energy insecurity

Here, this paper examines the relationships between energy efficient (EE) capital technology and household energy insecurity in the United States. The theoretical model of these relationships employs household production theory to capture the demand for and production of household energy services, and a stochastic production frontier approach to describe how having access to and the usage intensity of EE capital technology could help alleviate inefficiency in the production of household energy services. A working hypothesis formulated from the theoretical model posits that having EE capital technology in the home will reduce the level of household energy insecurity experienced. The extent of energy insecurity experienced is inferred from an energy insecurity index value assigned to each household, generated via the application of a dichotomous Rasch model to questions contained in the 2015 Residential Energy Consumption Survey. Noting the potential simultaneous relationship that exists between a household having access to and the usage intensity of EE capital technology and the experience of being energy insecure, an instrumental variables approach was employed to estimate a series of ordered logit models. Results suggest access to EE capital technology in the form of Energy Star® appliances, Energy Star® windows, or a SMART thermostat does not reduce the probability of experiencing a greater level of energy insecurity. Nor does the usage intensity of EE capital. Thus, policy instruments designed to alleviate household energy insecurity may need to go beyond simply helping households obtain EE capital technology.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2020

SAM Plug-in Development (Phase I Final Report)

The DOE Office of Nuclear Energy (NE) has created an extensive set of advanced modeling and simulation tools for nuclear engineering analysis. The advanced capabilities of these newer analysis codes require more in-depth training, skills, and knowledge in order to effectively utilize them for the design, analysis, and licensing of advanced nuclear systems and experiments. A high learning curve for inexperienced users may deter organizations from incorporating these tools into their internal processes. This project involved development of a plug-in to the Symbolic Nuclear Analysis Package (SNAP) for the System Analysis Module (SAM) tool. SAM is an advanced system analysis tool for reactor transient analyses being developed at Argonne National Laboratory under the U.S. DOE Office of Nuclear Energy’s Nuclear Energy Advanced Modeling and Simulation (NEAMS) program. SAM utilizes an object-oriented application framework (MOOSE), and its underlying meshing and finite-element library (libMesh) and linear and non-linear solvers (PETSc), to leverage modern advanced software environments and numerical methods. SNAP provides a highly flexible framework for creating, modifying and documenting input for engineering analysis codes such as SAM as well as extensive functionality for submitting, monitoring, and interacting with the codes through an intuitive graphical user interface (GUI). The common user interface provided by SNAP minimizes the learning curve for engineers starting with a new analysis code and provides an intuitive framework for transitioning between different analysis codes. SNAP provides a powerful but intuitive interface to facilitate access to advanced modeling and simulation tools for inexperienced users. Unlike many “form based” GUI’s, SNAP maps each engineering code’s component input to an internal database which manages all component input parameters along with component interconnections. This level of abstraction permits SNAP to support several advanced capabilities such as renodalization, model validation and consistency checks, embedded documentation, model notebook generation, data ownership and reviewer tracking, and variable assignment for inputs to name a few. SNAP includes a built-in Python interpreter and is interfaced to several commercial and open source packages including CPython, MATLAB/OCTAVE, Microsoft Office, Open Office, and SANDIA’s DAKOTA package which provides Uncertainty Quantification analysis through the SNAP plug-ins. Phase I of this project involved development a fully functional basic SAM plug-in to SNAP. This plug-in provides the ability to import existing models, graphically construct, edit and submit models using SNAP’s extensive functionality.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2023

2021 Radionuclide Air Emissions Report for Los Alamos National Laboratory, Revised October 2022

Provided is a revision to the Calendar Year 2021 Radionuclide Air Emissions Report for the Los Alamos National Laboratory (LANL). This report is a regulatory compliance deliverable for LANL under the Radionuclide NESHAP, 40 CFR 61 Subpart H. The revised report is needed to correct an error identified regarding generation of wind data files. These wind data files are used for plume modeling as part of calculating off-sites doses from monitored radionuclide air emissions sources. The error resulted from a software update by LANL’s meteorology team after new towers were added to the tower network. As a result of this error, incorrect wind frequencies were assigned to the different stability classes. This slightly affected the plume model calculations and subsequent off-site dose determinations for analyses using these wind files.

2021 Radionuclide, Air Emissions↗
DOE OSTI2020

Optimal Heliostat Assignment Strategy for Multiple-Receiver Systems

Next-generation concentrating solar power concepts currently being developed under the U.S. Department of Energy's Gen3 program include several novel heliostat field and receiver arrangements. Among these concepts are solar fields utilizing multiple receivers or several distinct receiver flux control zones, and each of which has particular power or peak flux constraints. In this paper, we present a linear program that identifies optimal heliostat-to-receiver assignment strategies for multi-receiver or multi-zone solar fields based on detailed calculations made by an optics modeling package. Two distinct problems are addressed: namely, (1) a design problem that chooses the optimal set of heliostats to include in a final layout given receiver constraints and performance, and (2) an operations problem that assigns existing heliostats to a receiver, maximizing field power output while maintaining receiver power constraints. The optimization methodology is integrated into the concentrating solar power tower modeling package SolarPILOTTM, and results are presented showing that heliostat assignment strategies can impact both initial heliostat field layout characteristics and the cost of thermal energy produced by the field. The methodology succeeds in maintaining power requirements but at the expense of producing different flux profiles on each receiver surface. Consequently, the paper also discusses additional measures that are applied to ensure that desirable flux distribution properties are maintained, reducing local flux mean-absolute-deviation values from a baseline of 17% or greater to less than 5% compared to a desired reference flux profile.

41 EE - Solar Energy Technologies Office (EE-4S)↗
DOE OSTI2025

3D modeling of n = 1 RMP driven heat fluxes on the SPARC tokamak PFCs using HEAT

3D heat flux calculations at the lower outer divertor plate of SPARC using the HEAT code show that 3D fields generated from error field correction coils can lead to enhanced peak heat fluxes up to 15 times larger compared to the axisymmetric case. Previously employed to simulate axisymmetric heat flux on 3D plasma facing components, the HEAT code can now predict 3D heat flux generated by non-axisymmetric plasmas. This is achieved via a new HEAT module which leverages the 3D field line tracing capabilities of MAFOT starting from an M3D-C1 (MHD resistive code) perturbed equilibria. The resulting heat flux is assigned using the magnetic footprint and the heat flux layer model, an extension of the 2D heat flux model also known as the Eich, to 3D non-axisymmetric plasmas. For SPARC, the new capabilities of HEAT are used to calculate the 3D heat loads resulting from n = 1 perturbation fields (with n indicating the toroidal periodicity) applied through a toroidal array of six picture frame coils with different amplitude. The comparison with the unperturbed case shows significant changes in shape and intensity of the heat flux profile. The results show that the application of n = 1 3D field leads to a localized enhancement of the heat flux peak, influenced by the wetted area impacted by the magnetic footprint, and the appearance of a secondary heat flux peak, whose intensity depends on amplitude of the applied 3D field and toroidal location.

3D heat flux↗
DOE OSTI2020

Defining Intermediates of Nitrogenase MoFe Protein during N 2 Reduction under Photochemical Electron Delivery from CdS Quantum Dots

Coupling the nitrogenase MoFe protein to light-harvesting semiconductor nanomaterials replaces the natural electron transfer complex with Fe protein and ATP and provides low potential photoexcited electrons for photocatalytic N 2 reduction. A central question is how direct photochemical electron delivery from nanocrystals to MoFe protein is able to support the multi-electron ammonia production reaction. In this study, low photon flux conditions were used to identify the initial reaction intermediates of CdS quantum dot (QD):MoFe protein nitrogenase complexes under photochemical activation. Illumination of CdS QD:MoFe protein complexes led to redox changes in the MoFe protein active site FeMo-co observed as the gradual decline in E0 resting state intensity that was accompanied by an increase in intensity of a new “g eff = 4.5” EPR signal. The magnetic properties of the g eff = 4.5 signal support assignment as a reduced S = 3/2 state, and reaction modeling was used to define it as a two-electron reduced “E 2 ” intermediate. Use of a MoFe protein variant, ß-188 Cys , which poises the P cluster in the oxidized P + state, demonstrated the P cluster can function as a site of photoexcited electron delivery from CdS to MoFe protein. Overall, the results establish the initial steps for how photoexcited CdS delivers electrons into MoFe protein during reduction of N 2 to ammonia, and the role of electron flux in the photochemical reaction cycle.

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

Identification of mobile genetic elements with geNomad

Identifying and characterizing mobile genetic elements in sequencing data is essential for understanding their diversity, ecology, biotechnological applications and impact on public health. Here we introduce geNomad, a classification and annotation framework that combines information from gene content and a deep neural network to identify sequences of plasmids and viruses. geNomad uses a dataset of more than 200,000 marker protein profiles to provide functional gene annotation and taxonomic assignment of viral genomes. Using a conditional random field model, geNomad also detects proviruses integrated into host genomes with high precision. In benchmarks, geNomad achieved high classification performance for diverse plasmids and viruses (Matthews correlation coefficient of 77.8% and 95.3%, respectively), substantially outperforming other tools. Leveraging geNomad’s speed and scalability, we processed over 2.7 trillion base pairs of sequencing data, leading to the discovery of millions of viruses and plasmids that are available through the IMG/VR and IMG/PR databases. geNomad is available at https://portal.nersc.gov/genomad.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Machine learning of 27Al NMR electric field gradient tensors for crystalline structures from DFT

NMR crystallography has emerged as a promising technique for the determination and refinement of atomic coordinates in crystal structures. The crystal structure of compounds containing quadrupolar nuclei, such as 27Al, can be improved by directly comparing solid-state NMR measurements to DFT computations of the electric field gradient (EFG) tensor. The non-negligible computational cost of these first-principles calculations limits the applicability of this method to all but the most well-defined structures. We developed a fast, low-cost machine learning model to predict EFG parameters based on local structural motifs and elemental parameters. We computed 8081 EFG tensors from 1681 27Al crystalline solids using DFT and benchmarked them against 105 experimentally measured 27Al sites. Surprisingly, simple local geometric features dominate the predictive performance of the resulting random-forest model, yielding an R2 value of 0.98 and an RMSE of 0.61 MHz for CQ, the quadrupolar coupling constant. This model accuracy should enable pre-refining future structural assignments before finally validating with first-principles calculations. Such a catalogue of 27Al NMR tensors can serve as a tool for researchers assigning complex NMR spectra influenced by the nuclear electric quadrupole interaction.

Sun, He↗