Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “computational screening”

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 487 records · Page 27

Small molecule BLVRB redox inhibitor promotes megakaryocytopoiesis and stress thrombopoiesis in vivo

Biliverdin IXβ reductase (BLVRB) is an NADPH-dependent enzyme previously implicated in a redox-regulated mechanism of thrombopoiesis distinct from the thrombopoietin (TPO)/c-MPL axis. Here, we apply computational modeling to inform molecule design, followed by de novo syntheses and screening of unique small molecules retaining the capacity for selective BLVRB inhibition as a novel platelet-enhancing strategy. Two distinct classes of molecules are identified, and NMR spectroscopy and co-crystallization studies confirm binding modes within the BLVRB active site and ring stacking between the nicotinamide moiety of the NADP+ cofactor. A diazabicyclo derivative displaying minimal off-target promiscuity and excellent bioavailability characteristics promotes megakaryocyte speciation in biphenotypic (erythro/megakaryocyte) cellular models and synergizes with TPO-dependent megakaryocyte formation in hematopoietic stem cells. Upon oral delivery into mice, this inhibitor expands platelet recovery in stress thrombopoietic models with no adverse effects. In this work, we identify and validate a cellular redox inhibitor retaining the potential to selectively promote megakaryocytopoiesis and enhance stress-associated platelet formation in vivo distinct from TPO receptor agonists.

36 MATERIALS SCIENCE↗

Energy variation of double K -shell photoionization of Ne

Herein we report on an experimental and theoretical study of double K-shell photoionization of Ne over the 2.3–8.5 keV x-ray energy range. The ratio of double-to-single K-shell photoionization cross sections was determined experimentally by measuring the relative rates of the KK–KL 2,3 L 2,3 ( 2 D) Auger hypersatellite and the K–L 2,3 L 2,3 ( 1 D) diagram Auger transitions. By scaling the hypersatellite-diagram Auger-electron ratios to KK/K cross-section ratios, comparison was made with theoretical cross-section ratios of He-like Ne 8+ determined by the R matrix with pseudostates method. The experimental Ne and theoretical Ne 8+ cross-section ratios show similar variations with energy, but the experimental ratios systematically exceed the calculated ratios and also show a lower threshold energy for double K-shell photoionization onset compared to the computed Ne 8+ threshold. The discrepancy is attributed to effects of L-shell electrons not included in the He-like Ne 8+ calculations. Quantified scaling with nuclear charge Z along the He-like isoelectronic sequence indicates that the measured 10-electron Z = 10 double K-shell photoionization cross section behaves like the computed He-like Z = 8.9 cross section, suggesting an effective nuclear screening parameter of s L =1.1 by the additional eight outer L-shell electrons. Experimental results for the energy variations of Auger electron transitions from other multielectron hole states are also discussed.

74 ATOMIC AND MOLECULAR PHYSICS↗

Automated AI-driven Molecular Design for Therapeutic Discovery

In recent years, artificial intelligence and machine learning (AI/ML) approaches have revolutionized the process of designing new therapeutics, enabling scientists to rapidly respond to emerging threats from various pathogens. A prime example is the SARS-CoV-2 main protease, a key target for the development of antiviral inhibitors. In this study, we employed a novel, integrated approach that combines AI-driven iterative design of inhibitor candidates, screening based on physio-chemical properties and toxicity, physics-based computational modeling of protein-inhibitor interactions, and AI-assisted analysis of Native MS biophysical assay and characterization of designed candidates. Our deep learning 3D-scaffold model, which uses an input scaffold as a starting point, generated tens of thousands of compounds while preserving the key scaffold. To optimize these candidates, we calculated a comprehensive set of 136 descriptors, including both 2D and 3D molecular features, for compounds targeting the SARS-CoV-2 Main protease (Mpro) and a neurodegenerative disease-associated protein, cyclophilin (Cyp). The generated compounds were initially filtered based on their properties and then ranked according to their predicted binding affinity using our automated modeling and ML methods. Experimental validation of the Mpro candidates showing inhibitory activity demonstrates that our workflow can expedite the therapeutic discovery.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The use of computer graphic techniques for the determination of ventricular function.

Description of computer techniques employed to increase the speed, accuracy, reliability, and scope of angiocardiographic analyses determining human heart dimensions. Chamber margins are traced with a Calma 303 digitizer from projections of the angiographic films. The digitized margins of the ventricular images are filed in a computer for subsequent analysis. The margins can be displayed on the television screen of a graphics unit for individual study or they can be viewed in real time (or at any selected speed) to study dynamic changes in the chamber outline. The construction of three dimensional images of the ventricle is described.

Sandler, H.↗

Numerical simulation of vortical flow over an elliptical-body missile at high angles of attack

Numerical solutions to the Reynolds-averaged Navier-Stokes equations are given for the flow about an elliptical body missile (3:1 ellipse) at a Mach number of 2.5 and a unit Reynolds number of 6.56 x 10 to the 6th/m. At high angles of attack, the flow is dominated by large-scale free vortices which occur in the lee-side flow field due to crossflow boundary-layer separation. Emphasis is focused on the accurate prediction of the lee-side vortical flow. Solutions are presented for both symmetric and asymmetric (body rolled 45 deg) configurations at 10 deg and 20 deg angle of attack. The computed results are compared with experimental surface pressure coefficients and vapor-screen photographs. Excellent agreement is obtained in all cases.

Newsome, R. N.↗

Real-time computer-generated hologram by means of liquid-crystal television spatial light modulator

The usefulness of an inexpensive liquid-crystal television) (LCTV) as a spatial light modulator for coherent-optical processing in the writing and reconstruction of a single computer-generated hologram has been demonstrated. The thickness nonuniformities of the LCTV screen were examined in a Mach-Zehnder interferometer, and the phase distortions were successfully removed using a technique in which the LCTV screen was submerged in a liquid gate filled with an index-matching nonconductive mineral oil with refractive index of about 1.45.

Mok, Fai↗

Aero/fluids database system

The AFAS Database System was developed to provide the basic structure of a comprehensive database system for the Marshall Space Flight Center (MSFC) Structures and Dynamics Laboratory Aerophysics Division. The system is intended to handle all of the Aerophysics Division Test Facilities as well as data from other sources. The system was written for the DEC VAX family of computers in FORTRAN-77 and utilizes the VMS indexed file system and screen management routines. Various aspects of the system are covered, including a description of the user interface, lists of all code structure elements, descriptions of the file structures, a description of the security system operation, a detailed description of the data retrieval tasks, a description of the session log, and a description of the archival system.

Reardon, John E.↗

Best Practices for Body Temperature Measurement with Infrared Thermography: External Factors Affecting Accuracy

Infrared thermographs (IRTs) are commonly used during disease pandemics to screen individuals with elevated body temperature (EBT). To address the limited research on external factors affecting IRT accuracy, we conducted benchtop measurements and computer simulations with two IRTs, with or without an external temperature reference source (ETRS) for temperature compensation. The combination of an IRT and an ETRS forms a screening thermograph (ST). We investigated the effects of viewing angle (θ, 0–75°), ETRS set temperature (T ETRS , 30–40 °C), ambient temperature (T atm , 18–32 °C), relative humidity (RH, 15–80%), and working distance (d, 0.4–2.8 m). We discovered that STs exhibited higher accuracy compared to IRTs alone. Across the tested ranges of T atm and RH, both IRTs exhibited absolute measurement errors of less than 0.97 °C, while both STs maintained absolute measurement errors of less than 0.12 °C. The optimal T ETRS for EBT detection was 36–37 °C. When θ was below 30°, the two STs underestimated calibration source (CS) temperature (T CS ) of less than 0.05 °C. The computer simulations showed absolute temperature differences of up to 0.28 °C and 0.04 °C between estimated and theoretical temperatures for IRTs and STs, respectively, considering d of 0.2–3.0 m, T atm of 15–35 °C, and RH of 5–95%. The results highlight the importance of precise calibration and environmental control for reliable temperature readings and suggest proper ranges for these factors, aiming to enhance current standard documents and best practice guidelines. These insights enhance our understanding of IRT performance and their sensitivity to various factors, thereby facilitating the development of best practices for accurate EBT measurement.

47 OTHER INSTRUMENTATION↗

Computational Advances in Ionic Liquid Applications for Green Chemistry: A Critical Review of Lignin Processing and Machine Learning Approaches

The valorization and dissolution of lignin using ionic liquids (ILs) is critical for developing sustainable biorefineries and a circular bioeconomy. This review aims to critically assess the current state of computational and machine learning methods for understanding and optimizing IL-based lignin dissolution and valorization processes reported since 2022. The paper examines various computational approaches, from quantum chemistry to machine learning, highlighting their strengths, limitations, and recent advances in predicting and optimizing lignin-IL interactions. Key themes include the challenges in accurately modeling lignin’s complex structure, the development of efficient screening methodologies for ionic liquids to enhance lignin dissolution and valorization processes, and the integration of machine learning with quantum calculations. These computational advances will drive progress in IL-based lignin valorization by providing deeper molecular-level insights and facilitating the rapid screening of novel IL-lignin systems.

09 BIOMASS FUELS↗

Supercomputer-Based Ensemble Docking Drug Discovery Pipeline with Application to Covid-19

In this work, we present a supercomputer-driven pipeline for in silico drug discovery using enhanced sampling molecular dynamics (MD) and ensemble docking. Ensemble docking makes use of MD results by docking compound databases into representative protein binding-site conformations, thus taking into account the dynamic properties of the binding sites. We also describe preliminary results obtained for 24 systems involving eight proteins of the proteome of SARS-CoV-2. The MD involves temperature replica exchange enhanced sampling, making use of massively parallel supercomputing to quickly sample the configurational space of protein drug targets. Using the Summit supercomputer at the Oak Ridge National Laboratory, more than 1 ms of enhanced sampling MD can be generated per day. We have ensemble docked repurposing databases to 10 configurations of each of the 24 SARS-CoV-2 systems using AutoDock Vina. Comparison to experiment demonstrates remarkably high hit rates for the top scoring tranches of compounds identified by our ensemble approach. We also demonstrate that, using Autodock-GPU on Summit, it is possible to perform exhaustive docking of one billion compounds in under 24 h. Finally, we discuss preliminary results and planned improvements to the pipeline, including the use of quantum mechanical (QM), machine learning, and artificial intelligence (AI) methods to cluster MD trajectories and rescore docking poses.

60 APPLIED LIFE SCIENCES↗

Program Helps Specify Paths For Numerically Controlled Tools

ESDAPT computer program provides graphical programming environment for developing APT (Automatically Programmed Tool) programs for controlling numerically controlled machine tools. Establishes graphical user interface providing user with APT syntax-sensitive text-editing subprogram and windows for displaying geometry and tool paths. APT geometry statements also created by use of menus and screen picks. Written in C language, yacc, lex, and XView for use on Sun4-series computers running SunOS.

Premack, Timothy↗

Prediction of O and OH Adsorption on Transition Metal Oxide Surfaces from Bulk Descriptors

In the search for stable and active catalysts, density functional theory and machine learning (ML) based models can accelerate the screening of materials. While stability is conveniently addressed on the bulk level of computation, the modelling of catalytic activity requires expensive surface simulations. Here, in this work, we develop models for the surface adsorption energy of O and OH intermediates across a consistent and extensive dataset of pure transition metal oxide surfaces. We show that adsorption energies across metal oxidation states of +2 to +6 are well captured from the metal-oxygen bond strength extracted from the bulk level calculation. Specifically, we calculate the integrated crystal orbital Hamiltonian population (ICOHP) of the metal-oxygen bond in the bulk oxide and employ a simple normalization scheme to obtain a strong correlation with adsorption energetics. By combining our ICOHP descriptor with non DFT features in a Gaussian Process regression (GPR) model, we achieve high model accuracy with mean absolute errors of 0.166 and 0.219 eV for OH and O adsorption, respectively. By targeting the O-OH adsorption energy difference with our GPR model, we predict the the oxygen evolution reaction (OER) activity from bulk descriptors only. Furthermore, we utilize the strong correlation between the COHP and metal oxygen bond lengths to rapidly predict adsorption energetics and catalytic activity from the optimized bulk geometry. Our approach can enable an efficient search for active catalysts by eliminating the need for surface calculations in the initial screening phase.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Risk Importance Ranking of Fire Data Parameters to Enhance Fire PRA Model Realism

Fire is historically and analytically a significant contributor to nuclear power plant risk. The level of fire risk and the methods, tools and data for modeling this risk is highly debated by experts. One area of debate is the input data used in fire modeling and how to deal with this data’s high uncertainty. This report outlines initial work performed for determining the key parameters causing this uncertainty and how it propagates into nuclear power plant models. This research paves the way for the development of methods to reduce fire data uncertainty used in modeling. The Nuclear Regulatory Commission has mandated that nuclear power plants perform fire risk modeling. However, there are several issues with the current risk modeling implementation that affect the results. Approved modeling methods can be overly conservative and often do not match plant experience. Also, the data used in the modeling can have high uncertainties and is influenced by expert judgement. To evaluate input data uncertainty, researchers performed an initial review of several fire experiments done at Sandia National Laboratories. Uncertainties for fire data can come from many sources, such as experiment design constraints, environmental conditions, or other plant-specific aspects. There are many different significant and insignificant parameters driving the uncertainty. Additionally, the uncertainty of the different input data used in the fire modeling could have a significant or insignificant effect on the entire plant risk. A four-step methodology was developed to perform Integrated Probabilistic Risk Assessment Importance Ranking. A demonstration case using these steps was set up and three of the four steps were completed in fiscal year (FY) 2019 and the fourth step done FY 2020. These steps are: 1. The qualitative analysis of potential sources was conducted with the following items identified for the demonstration. • Maximum heat release rate • Time to maximum heat release rate • Duration of max heat release rate • Time to decay • Thermal conductivity of concrete • Specific heat of concrete • Density of concrete • Cable jacket thickness 2. A quantitative characterization of dominant sources of uncertainty was performed. A list of distributions and determined values of the dominant sources is shown in Appendix A. 3. A quantitative screening of the potential sources of uncertainty using Morris Elementary Effects Analysis was completed. An experimental model using the physics-based fire modeling tool Fire Dynamics Simulator was developed and coupled with the Risk Analysis Virtual Environment. The Morris analysis identified at least two parameters that can be eliminated as significant contributors (specific heat of concrete and cable jacket thickness). 4. Global importance measure (Global IM) analysis to generate a comprehensive ranking based on their influence on the plant risk. In this research, a moment-independent Global IM is used since it can address (a) uncertainty in the input parameters of the fire model, (b) uncertainty in the risk outputs, and (c) non-linearity and interactions among input parameters in the fire model, more accurately than the correlation-based and variance-based global methods. The observations from the research showed that, depending on the initial and boundary conditions of the fire scenarios, fire-induced damage could have a very small probability and could be dominated by the tail of the uncertainty distribution; hence, the accuracy of the correlation-based and variance-based methods is questionable. The moment-independent Global IM analysis in this research provides a better understanding of how experimental uncertainty data affects industry’s plant models and where improvements in that data will have the largest benefit for improving fire modeling accuracy in causing core damage. Among the five unscreened parameters obtained from the Morris EE analysis, the Global IM analysis results for the case study indicated that max heat release rate and fire location are the most important parameters. The report also outlines benefits of using a unified computational platform that integrates the underlying simulations (e.g., a fire progression model), quantitative screening (using the Morris EE method), and the Global IM analysis. A unified platform can (i) facilitate the ranking of input parameters considering multiple key fire scenarios simultaneously, rather than considering one scenario at a time, (ii) contribute to more explicit and accurate treatment of dependencies at multiple levels of Fire PRA, (iii) facilitate the sampling-based uncertainty quantification for Fire PRA, and (iv) help generating both “industry-wide” and "plant-specific" ranking of uncertainty sources in Fire PRA. Future research should be done to include additional parameters such as detection/suppression or cable fire spread. Adding a PRA software such as SAPHIRE to the RAVEN platform would help with plant model integration and improve treatment of fire-induced dependency. The I-PRA risk importance ranking methodology offered in this report can provide valuable information for efficiently (a) enhancing the realism of Fire PRA for existing plants and (b) supporting the development of Dynamic Fire PRA for advanced reactors and new plants.

97 MATHEMATICS AND COMPUTING↗

Computational Design of Eutectic Molten Salt Mixtures: What Can Thermodynamic Models Do?

In the search for efficient energy storage battery technologies, designing stable electrolytes has been a long-standing challenge. Electrolytes based on molten salt eutectics are known for their stability with minimum parasitic reactions when compared to their widely used organic counterparts. However, the operating temperatures of these molten salt electrolyte-based batteries are dictated by the melting point of the eutectic mixtures. Design and high throughput screening of low melting temperature eutectic molten salt mixtures have been hindered by the lack of computational models. In this work, we develop thermodynamic models to predict the eutectic points of several molten salt mixtures. The framework of the COSMO-SAC model is used for the predictions and is compared with experimental data and other thermodynamic approaches. Rapid thermodynamics-based approaches, as shown in this study, can accelerate the discovery of new materials, complementing experimental techniques.

Ashwin Ravichandran↗

Machine learning dielectric screening for the simulation of excited state properties of molecules and materials

Accurate and efficient calculations of absorption spectra of molecules and materials are essential for the understanding and rational design of broad classes of systems. Solving the Bethe–Salpeter equation (BSE) for electron–hole pairs usually yields accurate predictions of absorption spectra, but it is computationally expensive, especially if thermal averages of spectra computed for multiple configurations are required. We present a method based on machine learning to evaluate a key quantity entering the definition of absorption spectra: the dielectric screening. We show that our approach yields a model for the screening that is transferable between multiple configurations sampled during first principles molecular dynamics simulations; hence it leads to a substantial improvement in the efficiency of calculations of finite temperature spectra. We obtained computational gains of one to two orders of magnitude for systems with 50 to 500 atoms, including liquids, solids, nanostructures, and solid/liquid interfaces. Importantly, the models of dielectric screening derived here may be used not only in the solution of the BSE but also in developing functionals for time-dependent density functional theory (TDDFT) calculations of homogeneous and heterogeneous systems. Overall, our work provides a strategy to combine machine learning with electronic structure calculations to accelerate first principles simulations of excited-state properties.

36 MATERIALS SCIENCE↗

Mechanochemical Solid Form Screening of Zeolitic Imidazolate Frameworks Using Structure-Directing Liquid Additives

We demonstrate a systematic application of the mechanochemical liquid-assisted grinding (LAG) methodology to screen for forms of zinc imidazolate (ZnIm 2 ), of fundamental importance as the simplest member of the zeolitic imidazolate framework materials family. The exploration of 45 different liquid additives, selected based on their molecular structure and physicochemical properties has resulted in eight different ZnIm 2 topological forms, appearing in 13 crystallographically distinct solid forms (including two previously unknown forms of the crb (BCT) topology), amorphous phases, and the interrupted moc-Zn 4 Im 8 HIm. All prepared topological forms were also explored computationally, using dispersion-corrected periodic density functional theory (DFT) calculations, enabling the rationalization of screening outcomes, and setting the stage for future prediction of additive-directed metal–organic framework (MOF) synthesis. This first systematic exploration of LAG in screening for three-dimensional MOFs demonstrates the potential of the liquid additive to not only accelerate materials synthesis, but also to direct it toward topologically different MOFs. The discovery of novel forms of a material that already exhibits at least 21 crystallographically and functionally different forms provides a strong testimony on the power of mechanochemistry in metal–organic materials discovery.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

BEAST DB: Grand-Canonical Database of Electrocatalyst Properties

We present BEAST DB, an open-source database comprised of ab initio electrochemical data computed using grand-canonical density functional theory in implicit solvent at consistent calculation parameters. The database contains over 20,000 surface calculations and covers a broad set of heterogeneous catalyst materials and electrochemical reactions. Calculations were performed at self-consistent fixed potential as well as constant charge to facilitate comparisons to the computational hydrogen electrode. This article presents common use cases of the database to rationalize trends in catalyst activity, screen catalyst material spaces, understand elementary mechanistic steps, analyze the electronic structure, and train machine learning models to predict higher fidelity properties. Users can interact graphically with the database by querying for individual calculations to gain a granular understanding of reaction steps or by querying for an entire reaction pathway on a given material using an interactive reaction pathway tool. BEAST DB will be periodically updated, with planned future updates to include advanced electronic structure data, surface speciation studies, and greater reaction coverage.

database↗

Dielectric screening and vacancy formation for large neutral and charged Si n H m ( n > 1500 ) nanocrystals using real-space pseudopotentials

A commonly used procedure for computing the properties of defects in crystalline materials is to consider a large supercell that includes the defect of interest. This is a straightforward technique as standard energy band codes can be used for such computations. For neutral defects, the only impediment of such an approach is to avoid defect-defect interactions between adjoining cells. However, this procedure can be complex if the defect of interest is charged as the system at large contains Coulombic divergences. Moreover, some have recently argued that the conventional definition of formation energies for charged defects cannot be reconciled with statistical mechanics. Here, we focus on an alternative approach. We consider large nanocrystals wherein a charged defect can be placed. Since the system is confined, a charged defect within the nanocrystal does not result in a Coulombic divergence. The chief impediment is computational, i.e., while no defect-defect or Coulombic divergences are present, the nanocrystal must be sufficiently large to allow the system to properly replicate a bulklike configuration. With the development of new algorithms and hardware advances, computations for systems of sufficient size to address this issue are feasible. In particular, we solve the Kohn-Sham equation in real space using pseudopotential-density-functional theory for large silicon nanocrystals, which contain thousands of atoms. Further, we focus on (i) the screening of a point charge and (ii) the formation of a charged vacancy in hydrogen-terminated silicon nanocrystals. This approach allows us to examine the role of quantum confinement in addition to exploring the bulk limit. Comparisons to other methods confirm the viability of this approach.

36 MATERIALS SCIENCE↗