Search NASA⌕ Search

SEARCH · Search NASA

Results for “computational frameworks”

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 451 records · Page 25

Atomic-scale understanding of oxide growth and dissolution kinetics of Ni-Cr alloys

Aqueous corrosion of metals is governed by formation and dissolution of a passivating, multi-component surface oxide. Unfortunately, a detailed atomistic description is challenging due to the compositional complexity and the need to consider multiple kinetic factors simultaneously. To this end, we combine experiments with a first-principles-derived, multiscale computational framework that transcends thermodynamic descriptions to explicitly simulate the kinetic evolution of surface oxides of Ni-Cr alloys as a function of composition, temperature, pH, and applied voltage. In the absence of pitting, we identify three distinct voltage regimes, which are kinetically dominated by oxide growth, dissolution, and competitive dissolution and reprecipitation. Evolving compositional gradients and oxide thickness are revealed, including a transition between a metastable Ni-Cr mixed oxide and a thick, porous Ni-dominated oxide. Beyond elucidating the underlying physics, we highlight the need for competing kinetics in models to properly predict the transition from passivation to corrosion. Our results provide a key step towards co-design of alloy composition alongside environmental conditions for sustainable use across a variety of critical energy and infrastructure applications.

36 MATERIALS SCIENCE↗

Integration of theory and experiment in the modelling of heterogeneous electrocatalysis.

Theoretical modelling is essential to deepen our understanding of heterogeneous electrocatalytic energy conversion processes, such as water splitting. Here, Sharon Hammes-Schiffer and Giulia Galli offer their perspectives on the best strategies for successfully studying such systems. Heterogeneous electrocatalysis is critical to many energy conversion processes. Theoretical and computational approaches are essential to interpret experimental data and provide the mechanistic understanding necessary to design more effective catalysts. However, automated general procedures to build predictive theoretical and computational frameworks are not readily available; specific choices must be made in terms of the atomistic structural model and the level of theory, as well as the experimental data used to inform and validate these choices. Here we outline some best practices for modelling heterogeneous systems and present examples in the context of catalysis at metal electrodes and oxides. The level of theory should be chosen for the specific system and properties of interest, and experimental validation is essential from the beginning to the end of the study. Continuous feedback and ultimate integration between experiment and theory enhances the power of calculations to elucidate mechanisms, identify effective descriptors and clarify design principles.

Hammes-Schiffer, Sharon↗

Adaptive optical correction for in vivo two-photon fluorescence microscopy with neural fields

Adaptive optics restore ideal imaging performance in complex samples by measuring and correcting optical aberrations but often require custom-built microscopes with carefully aligned wavefront sensing/shaping devices and can be susceptible to sample motion. Here we describe NeAT, a computational framework using neural fields for adaptive optics two-photon fluorescence microscopy. NeAT estimates wavefront aberration and recovers sample structure from a 3D image stack without requiring external datasets for training. Incorporating motion correction in learning and correcting conjugation errors commonly found in commercial microscopes, NeAT is designed for deployment in biological laboratories for in vivo imaging. We validate NeAT’s performance using a custom-built microscope with a wavefront sensor under varying signal-to-noise ratios, aberration and motion conditions. With a commercial microscope, we demonstrate real-time aberration correction for in vivo morphological and functional imaging in the living mouse brain, with NeAT improving the signal and accuracy of glutamate and calcium imaging of synapses and neurons.

Kang, Iksung↗

Subnetwork representation learning for discovering network biomarkers in predicting lymph node metastasis in early oral cancer

Cervical lymph node metastasis is the leading cause of poor prognosis in oral tongue squamous cell carcinoma and also occurs in the early stages. The current clinical diagnosis depends on a physical examination that is not enough to determine whether micrometastasis remains. The transcriptome profiling technique has shown great potential for predicting micrometastasis by capturing the dynamic activation state of genes. However, there are several technical challenges in using transcriptome data to model patient conditions: (1) An Insufficient number of samples compared to the number of genes, (2) Complex dependence between genes that govern the cancer phenotype, and (3) Heterogeneity between patients between cohorts that differ geographically and racially. We developed a computational framework to learn the subnetwork representation of the transcriptome to discover network biomarkers and determine the potential of metastasis in early oral tongue squamous cell carcinoma. Our method achieved high accuracy in predicting the potential of metastasis in two geographically and racially different groups of patients. The robustness of the model and the reproducibility of the discovered network biomarkers show great potential as a tool to diagnose lymph node metastasis in early oral cancer.

59 BASIC BIOLOGICAL SCIENCES↗

acCRISPR: an activity-correction method for improving the accuracy of CRISPR screens

Abstract High throughput CRISPR screens are revolutionizing the way scientists unravel the genetic underpinnings of engineered and evolved phenotypes. One of the critical challenges in accurately assessing screening outcomes is accounting for the variability in sgRNA cutting efficiency. Poorly active guides targeting genes essential to screening conditions obscure the growth defects that are expected from disrupting them. Here, we develop acCRISPR, an end-to-end pipeline that identifies essential genes in pooled CRISPR screens using sgRNA read counts obtained from next-generation sequencing. acCRISPR uses experimentally determined cutting efficiencies for each guide in the library to provide an activity correction to the screening outcomes via calculation of an optimization metric, thus determining the fitness effect of disrupted genes. CRISPR-Cas9 and -Cas12a screens were carried out in the non-conventional oleaginous yeast Yarrowia lipolytica and acCRISPR was used to determine a high-confidence set of essential genes for growth under glucose, a common carbon source used for the industrial production of oleochemicals. acCRISPR was also used in screens quantifying relative cellular fitness under high salt conditions to identify genes that were related to salt tolerance. Collectively, this work presents an experimental-computational framework for CRISPR-based functional genomics studies that may be expanded to other non-conventional organisms of interest.

59 BASIC BIOLOGICAL SCIENCES↗

Regularization by denoising diffusion models for solving inverse PDE problems with application to full waveform inversion

Partial differential equation (PDE)-governed inverse problems are fundamental across various scientific and engineering applications; yet they face significant challenges due to nonlinearity, ill-posedness, and sensitivity to noise. Here, we introduce a computational framework, regularization by denoising using diffusion models for partial differential equations (RED-DiffEq), by integrating physics-driven inversion and data-driven learning. RED-DiffEq leverages pretrained diffusion models as a regularization mechanism for PDE-governed inverse problems. We apply RED-DiffEq to solve the full waveform inversion problem in geophysics, a challenging seismic imaging technique that seeks to reconstruct high-resolution subsurface velocity models from seismic measurement data. Our method shows enhanced accuracy and robustness compared to benchmark methods. Additionally, it exhibits strong generalization and domain decomposition capacity, enabling the inversion of more complex velocity models with larger domains than those used in training the diffusion model. Our framework can also be directly applied to diverse PDE-governed inverse problems.

Shan, Siming [Yale University, New Haven, CT (Unit↗

Scientific machine learning benchmarks

Deep learning has transformed the use of machine learning technologies for the analysis of large experimental datasets. In science, such datasets are typically generated by large-scale experimental facilities, and machine learning focuses on the identification of patterns, trends and anomalies to extract meaningful scientific insights from the data. In upcoming experimental facilities, such as the Extreme Photonics Application Centre (EPAC) in the UK or the international Square Kilometre Array (SKA), the rate of data generation and the scale of data volumes will increasingly require the use of more automated data analysis. Furthermore, at present, identifying the most appropriate machine learning algorithm for the analysis of any given scientific dataset is a challenge due to the potential applicability of many different machine learning frameworks, computer architectures and machine learning models. Historically, for modelling and simulation on high-performance computing systems, these issues have been addressed through benchmarking computer applications, algorithms and architectures. Extending such a benchmarking approach and identifying metrics for the application of machine learning methods to open, curated scientific datasets is a new challenge for both scientists and computer scientists. Here, we introduce the concept of machine learning benchmarks for science and review existing approaches. As an example, we describe the SciMLBench suite of scientific machine learning benchmarks.

42 ENGINEERING↗

ServiceX A Distributed, Caching, Columnar Data Delivery Service

We will describe a component of the Intelligent Data Delivery Service being developed in collaboration with IRIS-HEP and the LHC experiments. ServiceX is an experiment-agnostic service to enable on-demand data delivery specifically tailored for nearly-interactive vectorized analysis. This work is motivated by the data engineering challenges posed by HL-LHC data volumes and the increasing popularity of python and Spark-based analysis workflows.ServiceX gives analyzers the ability to query events by dataset metadata. It uses containerized transformations to extract just the data required for the analysis. This operation is colocated with the data to avoid transferring unnecessary branches over the WAN. Simple filtering operations are supported to further reduce the amount of data transferred.Transformed events are cached in a columnar datastore to accelerate delivery of subsequent similar requests. ServiceX will learn commonly related columns and automatically include them in the transformation to increase the potential for cache hits by other users.Selected events are streamed to the analysis system using an efficient wire protocol that can be readily consumed by a variety of computational frameworks. This reduces time-to-insight for physics analysis by delegating to ServiceX the complexity of event selection, slimming, reformatting, and streaming.

Galewsky, B.↗

Informatics guided discovery of surface structure-chemistry relationships in catalytic nanoparticles

A data driven discovery strategy based on statistical learning principles is used to discover new correlations between electronic structure and catalytic activity of metal surfaces. From the quantitative formulations derived from this informatics based model, a high throughput computational framework for predicting binding energy as a function of surface chemistry and adsorption configuration that bypasses the need for repeated electronic structure calculations has been developed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

An informatics based analysis of the impact of isotope substitution on phonon modes in graphene

It is shown by informatics that the high frequency short ranged modes exert a significant influence in impeding thermal transport through isotope substituted graphene nanoribbons. Using eigenvalue decomposition methods, we have extracted features in the phonon density of states spectra that reveal correlations between isotope substitution and phonon modes. This study also provides a data driven computational framework for the linking of materials chemistry and transport properties in 2D systems.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Non-equilibrium reaction and relaxation dynamics in a strongly interacting explicit solvent: F + CD{sub 3}CN treated with a parallel multi-state EVB model

We describe a parallelized linear-scaling computational framework developed to implement arbitrarily large multi-state empirical valence bond (MS-EVB) calculations within CHARMM and TINKER. Forces are obtained using the Hellmann-Feynman relationship, giving continuous gradients, and good energy conservation. Utilizing multi-dimensional Gaussian coupling elements fit to explicitly correlated coupled cluster theory, we built a 64-state MS-EVB model designed to study the F + CD{sub 3}CN → DF + CD{sub 2}CN reaction in CD{sub 3}CN solvent (recently reported in Dunning et al. [Science 347(6221), 530 (2015)]). This approach allows us to build a reactive potential energy surface whose balanced accuracy and efficiency considerably surpass what we could achieve otherwise. We ran molecular dynamics simulations to examine a range of observables which follow in the wake of the reactive event: energy deposition in the nascent reaction products, vibrational relaxation rates of excited DF in CD{sub 3}CN solvent, equilibrium power spectra of DF in CD{sub 3}CN, and time dependent spectral shifts associated with relaxation of the nascent DF. Many of our results are in good agreement with time-resolved experimental observations, providing evidence for the accuracy of our MS-EVB framework in treating both the solute and solute/solvent interactions. The simulations provide additional insight into the dynamics at sub-picosecond time scales that are difficult to resolve experimentally. In particular, the simulations show that (immediately following deuterium abstraction) the nascent DF finds itself in a non-equilibrium regime in two different respects: (1) it is highly vibrationally excited, with ∼23 kcal mol{sup −1} localized in the stretch and (2) its post-reaction solvation environment, in which it is not yet hydrogen-bonded to CD{sub 3}CN solvent molecules, is intermediate between the non-interacting gas-phase limit and the solution-phase equilibrium limit. Vibrational relaxation of the nascent DF results in a spectral blue shift, while relaxation of the post-reaction solvation environment results in a red shift. These two competing effects mean that the post-reaction relaxation profile is distinct from what is observed when Franck-Condon vibrational excitation of DF occurs within a microsolvation environment initially at equilibrium. Our conclusions, along with the theoretical and parallel software framework presented in this paper, should be more broadly applicable to a range of complex reactive systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Contact stiffness and damping of liquid films in dynamic atomic force microscope

The mechanical properties and dissipation behaviors of nanometers confined liquid films have been long-standing interests in surface force measurements. The correlation between the contact stiffness and damping of the nanoconfined film is still not well understood. We establish a novel computational framework through molecular dynamics (MD) simulation for the first time to study small-amplitude dynamic atomic force microscopy (dynamic AFM) in a simple nonpolar liquid. Through introducing a tip driven dynamics to mimic the mechanical oscillations of the dynamic AFM tip-cantilever assembly, we find that the contact stiffness and damping of the confined film exhibit distinct oscillations within 6-7 monolayer distances, and they are generally out-of-phase. For the solid-like film with integer monolayer thickness, further compression of the film before layering transition leads to higher stiffness and lower damping, while much lower stiffness and higher damping occur at non-integer monolayer distances. These two alternating mechanisms dominate the mechanical properties and dissipation behaviors of simple liquid films under cyclic elastic compression and inelastic squeeze-out. Our MD simulations provide a direct picture of correlations between the structural property, mechanical stiffness, and dissipation behavior of the nanoconfined film.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Dynamic measurements and simulations of airborne picolitre-droplet coalescence in holographic optical tweezers

We report studies of the coalescence of pairs of picolitre aerosol droplets manipulated with holographic optical tweezers, probing the shape relaxation dynamics following coalescence by simultaneously monitoring the intensity of elastic backscattered light (EBL) from the trapping laser beam (time resolution on the order of 100 ns) while recording high frame rate camera images (time resolution <10 μs). The goals of this work are to: resolve the dynamics of droplet coalescence in holographic optical traps; assign the origin of key features in the time-dependent EBL intensity; and validate the use of the EBL alone to precisely determine droplet surface tension and viscosity. For low viscosity droplets, two sequential processes are evident: binary coalescence first results from the overlap of the optical traps on the time scale of microseconds followed by the recapture of the composite droplet in an optical trap on the time scale of milliseconds. As droplet viscosity increases, the relaxation in droplet shape eventually occurs on the same time scale as recapture, resulting in a convoluted evolution of the EBL intensity that inhibits quantitative determination of the relaxation time scale. Droplet coalescence was simulated using a computational framework to validate both experimental approaches. The results indicate that time-dependent monitoring of droplet shape from the EBL intensity allows for robust determination of properties such as surface tension and viscosity. Finally, the potential of high frame rate imaging to examine the coalescence of dissimilar viscosity droplets is discussed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A Rayleighian approach for modeling kinetics of ionic transport in polymeric media

Here, we report a theoretical approach for analyzing impedance of ionic liquids (ILs) and charged polymers such as polymerized ionic liquids (PolyILs) within linear response. The approach is based on the Rayleigh dissipation function formalism, which provides a computational framework for a systematic study of various factors, including polymer dynamics, in affecting the impedance. We present an analytical expression for the impedance within linear response by constructing a one-dimensional model for ionic transport in ILs/PolyILs. This expression is used to extract mutual diffusion constants, the length scale of mutual diffusion, and thicknesses of a low-dielectric layer on the electrodes from the broadband dielectric spectroscopy (BDS) measurements done for an IL and three PolyILs. Also, static dielectric permittivities of the IL and the PolyILs are determined. The extracted mutual diffusion constants are compared with the self diffusion constants of ions measured using pulse field gradient (PFG) fluorine nuclear magnetic resonance (NMR). For the first time, excellent agreements between the diffusivities extracted from the Electrode Polarization spectra (EPS) of IL/PolyILs and those measured using the PFG-NMR are found, which allows the use of the EPS and the PFG-NMR techniques in a complimentary manner for a general understanding of the ionic transport.

36 MATERIALS SCIENCE↗

Embedding non-collinear two-component electronic structure in a collinear quantum environment

Spin-containing materials are important for spintronic applications. In this work, we present a computational framework to embed nonrelativistic, two-component calculations in a one-component environment. In this framework, both embedding scalar potential and magnetic field can be included to describe the interaction between quantum subsystems. In this current development, a generalized Kohn-Sham density functional theory electronic structure is embedded in unrestricted Kohn-Sham density functional theory. Two test systems are studied, specifically: a Li3 on a closed-shell He lattice and a Li3 on a He lattice containing a Li atom defect. The noncollinearity of Li3 is unaffected upon embedding in a closed-shell environment through the scalar potential embedding. However, the open-shell nature of the Li atom defect introduces an effective magnetic field that couples to the magnetic components of the generalized Kohn-Sham Hamiltonian. These results show that noncollinear quantum embedding in an open-shell collinear environment may modify the spin structure of the embedded system. The formalism developed herein may serve as a useful tool in the modeling of inhomogeneous magnetic fields in two-component calculations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Predicting metal–metal interactions. II. Accelerating generalized schemes through physical insights

Operando-computational frameworks that integrate descriptors for catalyst stability within catalyst screening paradigms enable predictions of rates and selectivity on chemically faithful representations of nanoparticles under reaction conditions. These catalyst stability descriptors can be efficiently predicted by density functional theory (DFT)-based models. The alloy stability model, for example, predicts the stability of metal atoms in nanoparticles with site-by-site resolution. In this work, we use physical insights to present accelerated approaches of parameterizing this recently introduced alloy-stability model. These accelerated approaches meld quadratic functions for the energy of metal atoms in terms of the coordination number with linear correlations between model parameters and the cohesive energies of bulk metals. By interpolating across both the coordination number and chemical space, these accelerated approaches shrink the training set size for 12 fcc p- and d-block metals from 204 to as few as 24 DFT calculated total energies without sacrificing the accuracy of our model. We validate the accelerated approaches by predicting adsorption energies of metal atoms on extended surfaces and 147 atom cuboctahedral nanoparticles with mean absolute errors of 0.10 eV and 0.24 eV, respectively. This efficiency boost will enable a rapid and exhaustive exploration of the vast material space of transition metal alloys for catalytic applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Less can be more: Insights on the role of electrode microstructure in redox flow batteries from two-dimensional direct numerical simulations

Understanding how to structure a porous electrode to facilitate fluid, mass, and charge transport is key to enhancing the performance of electrochemical devices, such as fuel cells, electrolyzers, and redox flow batteries (RFBs). Here, using a parallel computational framework, direct numerical simulations are carried out on idealized porous electrode microstructures for RFBs. Strategies to improve an electrode design starting from a regular lattice are explored. By introducing vacancies in the ordered arrangement, it is possible to achieve higher voltage efficiency at a given current density, thanks to improved mixing of reactive species, despite reducing the total reactive surface. Careful engineering of the location of vacancies, resulting in a density gradient, outperforms disordered configurations. Our simulation framework is a new tool to explore transport phenomena in RFBs, and our findings suggest new ways to design performant electrodes.

25 ENERGY STORAGE↗

Multi-site reaction dynamics through multi-fragment density matrix embedding

The practical description of disordered chemical reactions, where the reactions involve multiple species at multiple sites, is presently a challenge using correlated electronic structure methods due to their high computational cost and steep scaling. Here, we describe the gradient theory of multi-fragment density matrix embedding theory, which potentially provides a minimal computational framework to model such processes at the correlated electron level. We present the derivation and implementation of the gradient theory, its validation on model systems and chemical reactions using density matrix embedding, and its application to a molecular dynamics simulation of proton transport in a small water cluster, a simple example of multi-site reaction dynamics.

Chemistry↗