Search NASA⌕ Search

SEARCH · Search NASA

Results for “computational challenges”

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 937 records · Page 52
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

MgSiP 2 : An Infrared Nonlinear Optical Crystal with a Large Non‐Resonant Phase‐Matchable Second Harmonic Coefficient and High Laser Damage Threshold

Superior infrared nonlinear optical (NLO) crystals are in urgent demand in the development of lasers and optical technologies for communications and computing. The critical challenge is to find a crystal with large non‐resonant phase‐matchable NLO coefficients and high laser damage threshold (LDTs) simultaneously, which however scale inversely. This work reports such a material, MgSiP 2, that exhibits a large second harmonic generation (SHG) coefficient of d 14 ≈ d 36 = 89 ± 5 pm V −1 at 1550 nm fundamental wavelength, surpassing the commercial NLO crystals AgGaS 2 , AgGaSe 2 , and ZnGeP 2 . First principles theory reveals the polarizability and geometric arrangement of the [SiP 4 ] tetrahedral units as the origin of this large nonlinear response. Remarkably, it also exhibits a high LDT value of 684 GW cm −2 , which is six times larger than ZnGeP 2 and three times larger than CdSiP 2 . It has a wide transparency window of 0.53–10.35 µm, allowing broadband tunability. Further, it is Type I and Type II phase‐matchable with large effective SHG coefficients of d eff,I ≈80.2 pm V −1 and d eff,II ≈73.4 pm V −1 . The outstanding properties of MgSiP 2 make it a highly attractive candidate for optical frequency conversion in the infrared.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Ligand enhanced activity of in situ formed nanoparticles for photocatalytic hydrogen evolution

One consistent challenge of both computational and empirical catalyst screening is ensuring that the variables chosen for the screen are driving the performance analyzed. Furthermore, we compare photocatalytic hydrogen evolution from in situ formed Au and Au/Cu nanoparticles to nanoparticles of the same composition synthesized prior to being used as catalysts, as well as compare them to in situ formed Au and Au/Cu nanoparticles in the presence of exogeneous ligand. For all experiments, we observed that ligand-terminated nanoparticles performed better than un-stabilized in situ formed nanoparticles. We tested the generality of this result by studying Co, Ni, and Pd in the same system and also observed that the introduction of nanoparticle ligand leads to enhanced catalytic activity. Taken together, these results suggest 1) nanoparticle ligands can be beneficial, and even necessary, to produce catalysts with sustained activity 2) For computational prediction of these catalysts, factors relating to particle formation and stability need to be considered when both generating predictions as well as interpreting experimental results based on those predictions

08 HYDROGEN↗
DOE OSTI2016

Charge distribution and hydrogen bonding of a collagen α 2 ‐chain in vacuum, hydrated, neutral, and charged structural models

A challenging task in computational biophysics is to ascertain the solvent effect on the electronic structure and interatomic bonding at the atomistic level. Simulations must be carried out on reasonably large biomolecules for accurate calculations to yield valid results. We report the results of a calculation on collagen model in the form of a peptide under three different environments: vacuum, solvated and with neutral and charged sites. Quantitative results and analysis of the partial charge (PC) distribution on each amino acid are discussed. A significant charge transfer of more than 1 electron from protein to water molecules is found with similar results when the model contains charged sites. The main contributions to the interatomic bonding are from hydrogen bonds (HBs) between water‐water and water‐protein pairs. A connection between PC and HBs can be established since the nonpolar amino acids form no HBs and have the smallest PC and vice versa. The ab initio PC obtained are used in the NAMD simulation showing significant improvement over the default values as reflected in the root mean square deviation of atomic positions in the MD steps and the total free energy in energy minimization. These results could facilitate the interpretation of data on interaction of various ligands in charged proteins in relation to isoelectric points. © 2016 Wiley Periodicals, Inc.

Eifler, Jay↗
DOE OSTI2024

Rethinking Programming Paradigms in the QC-HPC Context

Programming for today’s quantum computers is making significant strides toward modern workflows compatible with high performance computing (HPC), but fundamental challenges still remain in the integration of these vastly different technologies. Quantum computing (QC) programming languages share some common ground, as well as their emerging runtimes and algorithmic modalities. In this short paper, we explore avenues of refinement for the quantum processing unit (QPU) in the context of many-tasks management, asynchronous or otherwise, in order to understand the value it can play in linking QC with HPC. Through examples, we illustrate how its potential for scientific discovery might be realized.

Wong, Elaine↗
DOE OSTI2022

Approximate symmetries in d = 4 CFTs with an Einstein gravity dual

By applying the stress-tensor-scalar operator product expansion (OPE) twice, we search for algebraic structures in d = 4 conformal field theories (CFTs) with a pure Einstein gravity dual. We find that a rescaled mode operator defined by an integral of the stress tensor T ++ on a d = 2 plane satisfies a Virasoro-like algebra when the dimension of the scalar is large. The structure is enhanced to include a Kac-Moody-type algebra if we incorporate the T -- component. In our scheme, the central terms are finite. It remains challenging to directly compute the stress-tensor sector of d = 4 scalar four-point functions at large central charge, which, based on holography and bootstrap methods, were recently shown to have a Virasoro/ $\mathcal{W}$-algebra vacuum block-like structure.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Geophysical Constraints on Decarbonized Systems—Building Spatio-Temporal Uncertainties into Future Electricity Grid Planning

Future electricity grids will be characterized by the high penetration of renewables to support the decarbonization process. Yet, this transition will further expose grids to a broad spectrum of geophysical forces, such as weather and climate, land and space, minerals and materials. Understanding their role is thus important to support the grid decarbonization process. Here, we synthesize the current body of knowledge on the relationship between geophysical constraints and electricity grid planning. Additionally, we show that there have been promising advances in the data, methods, and modelling tools needed to incorporate the effect of geophysical constraints on demand, resource availability, and grid operations. There are, however, multiple research avenues that deserve further attention. More system-specific and finer-scale analyses are necessary to better understand how spatio-temporal variability---and associated uncertainty---in geophysical forces affect grid planning. Moreover, we need a broader focus on the multi-sectoral implications of grid decarbonization efforts, so as to limit the risk of unintended consequences. Importantly, all these efforts are challenged by the computational requirements of existing power system models, which often limit our ability to characterize uncertainty, scale analyses across larger domains, or study the multi-sector implications of grid expansion decisions.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2021

Emulating the impact of additional proton-proton interactions in the ATLAS simulation by pre-sampling sets of inelastic Monte Carlo events

The accurate simulation of additional interactions at the ATLAS experiment for the analysis of proton-proton collisions delivered by the Large Hadron Collider presents a significant challenge to the computing resources. During the LHC Run (2015-2018) there were up to 70 inelastic interactions per bunch crossing, which need to be accounted for in Monte Carlo (MC) production. In this document, a new method to account for these additional interactions in the simulation chain is described. Instead of sampling the inelastic interactions and adding their energy deposits to a hard-scatter interaction one-by-one, the inelastic interactions are presampled, independent of the hard scatter, and stored as combined events. Consequently, for each hard-scatter interaction only one such presampled event needs to be added as part of the simulation chain. For the Run 2 simulation chain, with an average of 35 interactions per bunch crossing, this new method provides a substantial reduction in MC production CPU needs of around 20%, while reproducing the properties of the reconstructed quantities relevant for physics analyses with good accuracy.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

A deep learning and finite element approach for exploration of inverse structure–property designs of lightweight hybrid composites

Hybrid composites have important applications, such as high-performance and lightweight materials in aerospace and automotive industries. Hybrid composites utilize the synergy of diverse fillers to achieve desired material properties, but usually have more complicated microstructures. While topology optimization can optimize a particular property, designing hybrid composites for customized mechanical performances, e.g. full-range stress–strain curve, remains challenging. Here, a computational framework that integrated finite element analysis (FEA) and artificial intelligence (AI) methods of Conditional Generative Adversarial Networks (cGAN) deep learning and transfer learning was developed to establish inverse structure–property relationships and design tailor-made hybrid composites. Based on FEA-generated datasets of hybrid fiber-particle–matrix microstructures and their corresponding full-range stress–strain curves, a cGAN architecture was trained to generate tailored microstructures and establish structure–property relationships. Similarity in microstructural features and well-matched stress–strain curves based on the AI-generated composites were achieved. In conclusion, transfer learning was used to expand the pre-trained model for designing different materials systems.

Hybrid composites↗
DOE OSTI2015

A sub-pore model for multi-scale reaction–diffusion problems in porous media

Applications of reaction–diffusion systems in porous media pose a challenging problem for computational modeling approaches due to their multi-physics and multi-scale nature. The length scales usually span 3–4 orders of magnitude while physical phenomena involved include heat and mass transfer processes, and surface reactions. In this paper, a novel methodology that accounts for all the length scales and physical phenomena involved in a single framework is described. A length scale based dual approach is proposed – the larger pore channels (macro-pores) are resolved using conventional numerical techniques and a novel ‘sub-pore’ model is used to account for the unresolved pore channels (sub-pores) and the important physics therein. The porous network in the sub-pore system is composed of a fractal-like hierarchical system of straight cylindrical pores. Simplified governing equations for mass and energy transport are solved within the sub-pore system along with a reaction kinetics model to account for surface adsorption. An implicit coupling strategy is used to couple the macro-pore and the sub-pore systems so as to ensure conservation. The developed methodology is then applied to a few test cases and it is established that the proposed framework is necessary for problems where the adsorption time scale is much smaller than (diffusion-limited) or comparable to the diffusion time scale. Furthermore, it is also demonstrated that the framework can be potentially used to model the network of porous channels in its entirety thus significantly reducing computational costs.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Assessing Implicit and Explicit Polarizable Solvation Models for Nuclear–Electronic Orbital Systems: Quantum Proton Polarization and Solvation Energetics

Accurate simulations of many chemical processes require the inclusion of both nuclear quantum effects and a solvent environment. The nuclear–electronic orbital (NEO) approach, which treats electrons and select nuclei quantum mechanically on the same level, combined with a polarizable continuum model (PCM) for the solvent environment, addresses this challenge in a computationally practical manner. In this work, the NEO-PCM approach is extended beyond the IEF-PCM (integral equation formalism PCM) and C-PCM (conductor PCM) approaches to the SS(V)PE (surface and simulation of volume polarization for electrostatics) and ddCOSMO (domain decomposed conductor-like screening model) approaches. IEF-PCM, SS(V)PE, C-PCM, and ddCOSMO all exhibit similar solvation energies as well as comparable nuclear polarization within the NEO framework. The calculations show that the nuclear density does not leak out of the molecular cavity because it is much more localized than the electronic density. Finally, the polarization of quantized protons is analyzed in both continuum solvent and explicit solvent environments described by the polarizable MB-pol model, illustrating the impact of specific hydrogen-bonding interactions captured only by explicit solvation. Furthermore, these calculations highlight the relationship among solvation formalism, nuclear polarization, and energetics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Uncovering Highly-Excited State Mixing in Acetone Using Ultrafast VUV Pulses and Coincidence Imaging Techniques

Here, understanding the ultrafast dynamics of highly-excited electronic states of small molecules is critical for a better understanding of atmospheric and astrophysical processes, as well as for designing coherent control strategies for manipulating chemical dynamics. In highly excited states, nonadiabatic coupling, electron-electron interactions, and the high density of states govern dynamics. However, these states are computationally and experimentally challenging to access. Fortunately, new sources of ultrafast vacuum ultraviolet pulses, in combination with electron-ion coincidence spectroscopies, provide new tools to unravel the complex electronic landscape. Here we report time-resolved photoelectron-photoion coincidence experiments using 8 eV pump photons to study the highly excited states of acetone. We uncover for the first time direct evidence that the resulting excited state consists of a mixture of both n y → 3p and π → π* character, which decays with a time constant of 330 fs. In the future, this approach can inform models of VUV photochemistry and aid in designing coherent control strategies for manipulating chemical reactions.

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

Topological Superconductivity Based on Antisymmetric Spin–Orbit Coupling

Topological superconductivity (TSC) has drawn much attention for its fundamental interest and application in quantum computation. An outstanding challenge is the lack of intrinsic TSC materials with a p-wave pairing gap, which has led to the development of an effective p-wave theory of coupling s-wave gap with Rashba spin–orbit coupling (RSOC). However, the RSOC-strict mechanism and materials pose still both fundamental and practical limitations. Here, we generalize this theory to antisymmetric SOC (ASOC). Using k·p perturbation theory, we demonstrate that 2D crystals, with point groups of C 2 , C 4 , C 6 , C 2v , C 4v , C 6v , D 2 , D 4 , D 6 , S 4 , or D 2d , can all facilitate the desired ASOC. Remarkably, this enables us to discover 314 TSC candidates by screening 2D material databases, which are further confirmed by first-principles calculations of Majorana boundary modes and the topological invariant of the superconducting gap. As a result, our work fundamentally enriches TSC theory and greatly expands the classes of TSC materials for experimental exploration.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Unraveling Adsorbate-Induced Structural Evolution of Iron Carbide Nanoparticles

Iron carbide (Fe x C y ) nanoparticles (NPs) are promising candidates for replacing platinum group metals in industrial applications, such as high-temperature Fischer–Tropsch synthesis. However, due to their amorphous nature, characterization of the active sites has been challenging experimentally and computationally. Here, using a combined density functional theory (DFT), neural network interatomic potential-assisted global optimization, and ensemble learning study, we evaluate dynamic surface changes associated with syngas (H and CO) interactions. For this purpose, we have developed a general procedure that we use to model an experimentally relevant 270-atom Fe 182 C 88 NP using the neural network-assisted stochastic surface walk global optimization algorithm (SSW-NN). Once generated, the Fe 182 C 88 NP active sites and particle morphology are thoroughly characterized before the effects of syngas adsorbate interactions are explored by using DFT and molecular dynamics simulations. Lastly, we explore correlations between geometric and electronic features of the active sites and the adsorption of H (H ads ), using a regularized random forest machine learning algorithm. In doing so, we identified the Fe–C coordination number and p orbital occupancy as the most important descriptors affecting H ads . Furthermore, using a combined ML and quantum chemistry approach, our work demonstrates a general and efficient procedure for generating and probing complex surface phenomena on binary nanoparticles.

Adsorption↗
DOE OSTI2024

A Diffusion‐Based Uncertainty Quantification Method to Advance E3SM Land Model Calibration

Abstract Calibrating land surface models and accurately quantifying their uncertainty are crucial for improving the reliability of simulations of complex environmental processes. This, in turn, advances our predictive understanding of ecosystems and supports climate‐resilient decision‐making. Traditional calibration methods, however, face challenges of high computational costs and difficulties in accurately quantifying parameter uncertainties. To address these issues, we develop a diffusion‐based uncertainty quantification (DBUQ) method. Unlike conventional generative diffusion methods, which are computationally expensive and memory‐intensive, DBUQ innovates by formulating a parameterized generative model and approximates this model through supervised learning, which enables quick generation of parameter posterior samples to quantify its uncertainty. DBUQ is effective, efficient, and general‐purpose, making it suitable for site‐specific ecosystem model calibration and broadly applicable for parameter uncertainty quantification across various earth system models. In this study, we applied DBUQ to calibrate the Energy Exascale Earth System Model land model at the Missouri Ozark AmeriFlux forest site. Results indicated that DBUQ produced accurate parameter posterior distributions similar to those from Markov Chain Monte Carlo sampling but with 30 times less computing time. This significant improvement in efficiency suggests that DBUQ can enable rapid, site‐level model calibration at a global scale, enhancing our predictive understanding of climate impacts on terrestrial ecosystems.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2018

Epitaxial growth of ultraflat stanene with topological band inversion

Here, two-dimensional (2D) topological materials, including quantum spin/anomalous Hall insulators, have attracted intense research efforts owing to their promise for applications ranging from low-power electronics and high-performance thermoelectrics to fault-tolerant quantum computation. One key challenge is to fabricate topological materials with a large energy gap for room-temperature use. Stanene—the tin counterpart of graphene—is a promising material candidate distinguished by its tunable topological states and sizeable bandgap. Recent experiments have successfully fabricated stanene, but none of them have yet observed topological states. Here we demonstrate the growth of high-quality stanene on Cu(111) by low-temperature molecular beam epitaxy. Importantly, we discovered an unusually ultraflat stanene showing an in-plane s–p band inversion together with a spin–orbit-coupling-induced topological gap (~0.3 eV) at the Γ point, which represents a foremost group-IV ultraflat graphene-like material displaying topological features in experiment. The finding of ultraflat stanene opens opportunities for exploring two-dimensional topological physics and device applications.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Evidence of hydrogen-helium immiscibility at Jupiter-interior conditions

The phase behaviour of warm dense hydrogen–helium (H–He) mixtures affects our understanding of the evolution of Jupiter and Saturn and their interior structures. For example, precipitation of He from a H–He atmosphere at about 1–10 megabar and a few thousand kelvin has been invoked to explain both the excess luminosity of Saturn, and the depletion of He and neon (Ne) in Jupiter’s atmosphere as observed by the Galileo probe. But despite its importance, H–He phase behaviour under relevant planetary conditions remains poorly constrained because it is challenging to determine computationally and because the extremes of temperature and pressure are difficult to reach experimentally. Here we report that appropriate temperatures and pressures can be reached through laser-driven shock compression of H 2 –He samples that have been pre-compressed in diamond-anvil cells. This allows us to probe the properties of H–He mixtures under Jovian interior conditions, revealing a region of immiscibility along the Hugoniot. A clear discontinuous change in sample reflectivity indicates that this region ends above 150 gigapascals at 10,200 kelvin and that a more subtle reflectivity change occurs above 93 gigapascals at 4,700 kelvin. Considering pressure–temperature profiles for Jupiter, these experimental immiscibility constraints for a near-protosolar mixture suggest that H–He phase separation affects a large fraction—we estimate about 15 per cent of the radius—of Jupiter’s interior. Furthermore, this finding provides microphysical support for Jupiter models that invoke a layered interior to explain Juno and Galileo spacecraft observations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Variational multiscale reinforcement learning for discovering reduced order closure models of nonlinear spatiotemporal transport systems

Abstract A central challenge in the computational modeling and simulation of a multitude of science applications is to achieve robust and accurate closures for their coarse-grained representations due to underlying highly nonlinear multiscale interactions. These closure models are common in many nonlinear spatiotemporal systems to account for losses due to reduced order representations, including many transport phenomena in fluids. Previous data-driven closure modeling efforts have mostly focused on supervised learning approaches using high fidelity simulation data. On the other hand, reinforcement learning (RL) is a powerful yet relatively uncharted method in spatiotemporally extended systems. In this study, we put forth a modular dynamic closure modeling and discovery framework to stabilize the Galerkin projection based reduced order models that may arise in many nonlinear spatiotemporal dynamical systems with quadratic nonlinearity. However, a key element in creating a robust RL agent is to introduce a feasible reward function, which can be constituted of any difference metrics between the RL model and high fidelity simulation data. First, we introduce a multi-modal RL to discover mode-dependant closure policies that utilize the high fidelity data in rewarding our RL agent. We then formulate a variational multiscale RL (VMRL) approach to discover closure models without requiring access to the high fidelity data in designing the reward function. Specifically, our chief innovation is to leverage variational multiscale formalism to quantify the difference between modal interactions in Galerkin systems. Our results in simulating the viscous Burgers equation indicate that the proposed VMRL method leads to robust and accurate closure parameterizations, and it may potentially be used to discover scale-aware closure models for complex dynamical systems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Data-driven causal model discovery and personalized prediction in Alzheimer's disease

Abstract With the explosive growth of biomarker data in Alzheimer’s disease (AD) clinical trials, numerous mathematical models have been developed to characterize disease-relevant biomarker trajectories over time. While some of these models are purely empiric, others are causal, built upon various hypotheses of AD pathophysiology, a complex and incompletely understood area of research. One of the most challenging problems in computational causal modeling is using a purely data-driven approach to derive the model’s parameters and the mathematical model itself, without any prior hypothesis bias. In this paper, we develop an innovative data-driven modeling approach to build and parameterize a causal model to characterize the trajectories of AD biomarkers. This approach integrates causal model learning, population parameterization, parameter sensitivity analysis, and personalized prediction. By applying this integrated approach to a large multicenter database of AD biomarkers, the Alzheimer’s Disease Neuroimaging Initiative, several causal models for different AD stages are revealed. In addition, personalized models for each subject are calibrated and provide accurate predictions of future cognitive status.

Zheng, Haoyang (ORCID:0000000168358242)↗