Search NASA⌕ Search

SEARCH · Search NASA

Results for “Random Interactions”

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 721 records · Page 40
All sourcesNASA NTRSDOE OSTI
NASA NTRS1999

What is the Relationship between the Solar Wind and Storms/Substorms?

The interplanetary magnetic field (IMF) carried past the Earth by the solar wind has long been known to be the principal quantity that controls geomagnetic storms and substorms. Intervals of strong southward IMF with durations of at least a significant fraction of a day produce storms, while more typical, shorter intervals of less-intense southward fields produce substorms. The strong, long-duration southward fields are generally associated with coronal mass ejections and magnetic clouds or else they are produced by interplanetary dynamics initiated by fast solar wind flows that compress preexisting southward fields. Smaller, short-duration southward fields that occur on most days are related to long period waves, turbulence, or random variations in the IMF. Southward IMF enhances dayside reconnection between the IMF and the Earth's dipole with the reconnected field lines supplementing open field lines of the geomagnetic tail and producing an expanded polar cap and increased tail energy. Although the frequent storage of solar wind energy and its release during substorms is the most common mode of solar wind/magnetosphere interaction, under certain circumstances, steady southward IMF seems to produce intervals of relatively steady magnetosphere convection without substorms. During these latter times, the inner magnetosphere remains in a stressed tail-like state while the more distant magnetotail has larger northward field and more dipolar-like field lines. Recent evidence suggests that enhanced magnetosphere particle densities associated with enhanced solar wind densities allow more particles to be accelerated for the ring current, thus creating larger storms.

Fairfield, D. H.↗
NASA NTRS1994

A Study of Parameters Affecting Fibroblast Morphology in Response to an Applied Mechanical Force

A precisely controlled stretch/relaxation regimen (20% elongation at 6.6 cycles/min) was applied to normal human fetal, neonatal and aged dermal fibroblasts cultured on flexible membranes. Culture conditions included poly (NH2) or collagen type I coated substrate membranes; control cultures were grown on the same pliable material in the absence of applied stretch. Direct observation and immunofluorescence analyses revealed a progressive change in cell body orientation limited to the stretched dermal fibroblast cultures. Monolayers gradually (over 4 days) acquired a symmetric, radial distribution equivalent to the biaxial array of the applied force. At high seeding density, alignment was inhibited in the fetal cell cultures. This cell strain required collagen type I coating for optimal attachment to the flexible membrane, preferring growth in three-dimensional cell 'balls' on the poly(NH2) coated substrate. Neonatal cells also required the collagen type I coating, but both neonatal and aged dermal fibroblasts aligned efficiently at all seeding densities examined. The randomly oriented neonatal cells on the unstretched control membranes spontaneously detached at confluence, as a single cell sheet. Their aligned counterparts did not detach until the applied stretch stimulus was removed. Low concentrations of cytochalasin D (62.5 ng/ml) disrupted the stretch-related alignment response. Rhodamine phalloidin staining visualized fewer actin stress fibers in stretched, aligned cells than in controls. Both intercellular interactions and cytoskeletal integrity mediate the response to mechanical strain. Normal rabbit corneal stroma fibroblasts (NRC) were also analyzed, and failed to orient under these conditions. This cell type may require a different regimen, or a longer time period, to demonstrate alignment behavior. Supported by NASA Space Biology RTOP 199-40-22 and the NASA-ARC Director's Discretionary Fund.

Grymes, Rosalind A.↗
NASA NTRS2010

High-Intensity Radiated Field Fault-Injection Experiment for a Fault-Tolerant Distributed Communication System

Safety-critical distributed flight control systems require robustness in the presence of faults. In general, these systems consist of a number of input/output (I/O) and computation nodes interacting through a fault-tolerant data communication system. The communication system transfers sensor data and control commands and can handle most faults under typical operating conditions. However, the performance of the closed-loop system can be adversely affected as a result of operating in harsh environments. In particular, High-Intensity Radiated Field (HIRF) environments have the potential to cause random fault manifestations in individual avionic components and to generate simultaneous system-wide communication faults that overwhelm existing fault management mechanisms. This paper presents the design of an experiment conducted at the NASA Langley Research Center's HIRF Laboratory to statistically characterize the faults that a HIRF environment can trigger on a single node of a distributed flight control system.

Yates, Amy M.↗
NASA NTRS2012

Paramyxovirus Infection Mimics In Vivo Cellular Dynamics in Three-Demensional Human Bronchio-Epithelial Tissue-Like Assemblies

Respiratory syncytial virus and parainfluenza virus cause severe respiratory disease, especially in infants, children and the elderly. An in vitro model that accurately mimics infection of the human respiratory epithelium (HRE) would facilitate vaccine development greatly. Monolayer cultures traditionally used to study these viruses do not accurately and precisely differentiate the replication efficiencies of wild type and attenuated viruses. Therefore, we engineered novel three-dimensional (3D) tissue-like assemblies (TLAs) of human broncho-epithelial (HBE) cells to produce a more physiologically relevant in vitro model of the HRE. TLAs resemble HRE structurally and by expression of differentiated epithelial cell markers. Most significantly, wild type viruses exhibited a clear growth advantage over attenuated strains in TLAs unlike monolayer cultures. In addition, the TLAs responded to virus infection by secreting pro-inflammatory mediators similar to the respiratory epithelia of infected children. These characteristics make the TLA model a valuable platform technology to develop and evaluate live, attenuated respiratory virus vaccine candidates for human use. Respiratory virus diseases, the most frequent and least preventable of all infectious diseases, range in severity from the common cold to severe bronchiolitis and pneumonia . Two paramyxoviruses, respiratory syncytial virus (RSV) and parainfluenza virus type 3 (PIV3), are responsible for a majority of the most severe respiratory diseases of infants and young children. RSV causes 70% of all bronchiolitis cases and is a major cause of morbidity and mortality worldwide, especially in infants. PIV3 causes 10-15% of bronchiolitis and pneumonia during infancy, second only to RSV, and 40% of croup in infants To date, licensed vaccines are not available to prevent these respiratory diseases. At present, traditional monkey kidney (Vero and LLC-MK2) and human (HEp-2) tissue culture cells and small animal models (mouse, cotton rat, guinea pig, ferret, and hamster) fail to accurately imitate viral replication and human disease states (8). Lacking an authentic model has impeded the development and evaluation of live, attenuated vaccine candidates. Development of a physiologically relevant in vitro tissue culture model that reproduces characteristics of the HRE, the primary target of RSV and PIV3, would aid in predicting clinical attenuation and safety of vaccine candidates. Successful tissue engineering of a 3D human intestinal model using novel NASA technology inspired the development of a tri-culture 3D model for the HRE. Sequential layering of primary mesenchymal cells (comprised of normal human fibroblasts and endothelial cells) followed by BEAS-2B epithelial cells derived from human bronchi and tracheae were recapitulated on Cultisphere and/or cytodex3 microcarriers in cylindrical vessels that rotate horizontally creating an organized epithelial structure. Horizontal rotation randomizes the gravity vector modeling aspects of microgravity. Mesenchymal and epithelial cells grown under these conditions reproduce the structural organization, multi-cellular complexity, and differentiation state of the HRE. The opportunity to study respiratory viruses in a nasal epithelium model is invaluable because the most promising respiratory virus vaccine candidates are live attenuated viruses for intranasal administration. Here we characterize the interactions of respiratory viruses and epithelial cells grown under modeled microgravity in comparison to gravity-ladened monolayers. 3D HBE TLAs and traditional monolayers (2D) are infected at 35 C, the upper temperature of the upper HRE, to simulate in vivo infection conditions. Growth kinetics of wild type (wt) RSV and PIV3 viruses were compared in 2D and 3D cells to that of strains attenuated in humans or rhesus macaques. This novel 3D HBE model also offers an opportunity to study whether the epithelial cell function, especially in host defenses recapitulated by mimicking the structural organization of the HRE. In vivo, airway epithelial cells play a significant and dynamic role in host defense by blocking paracellular permeability and modulating airway function through cellular interactions or tight junctions. As regulators of the innate immune response, epithelial cells constitutively express cytokines, chemokines, and colony stimulating factors including RANTES, IL-8, IL-6, GM-CSF, and G-CSF for proactive host defense. In response to viral infection, epithelial cells induce potent immuno-modulatory and pro-inflammatory cytokines that recruit phagocytic and inflammatory cells to clear the virus and enhance protection. Although disease pathogenesis is classically attributed to the cytopathic effects of the pathogen, severe disease states associated with RSV and PIV3 are attributed to the inflammatory response, especially in infants. RSV is a potent inducer of cytokines and pro-inflammatory mediators in epithelial cells in vivo. A differentiated human epithelial model independent of the complete functional immune system will help elucidate the role of epithelial cells in respiratory disease. We reported here, virus and host cell interactions in 3D HBE TLAs are similar to that in vivo. Because the epithelial cell organization of the TLAs impacts not only the expression of airway epithelial characteristics, but also cellular communication, the TLAs represent a more physiologically relevant model of the HRE than BEAS-2B or other non-tumour monolayer models of respiratory disease. As a result, wild type respiratory viruses have a clear growth advantage over attenuated viruses in TLAs unlike traditional monolayers. In addition, the TLAs respond to wild type virus infection by secreting pro-inflammatory mediators characteristic of infected HRE. TLAs expressing microbial defense mechanisms provide an excellent model to study the interactions of respiratory pathogens with their host and to identify the innate immunity mediators. Therefore, 3D HBE TLAs offer advantages for the study of respiratory viruses and the development of viral vaccine candidates.

Deatly, Anne M.↗
DOE OSTI2024

Slow thermalization and subdiffusion in $U$(1) conserving Floquet random circuits

Random quantum circuits are paradigmatic models of minimally structured and analytically tractable chaotic dynamics. We study a family of Floquet unitary circuits with Haar random U(1) charge conserving dynamics; the minimal such model has nearest-neighbor gates acting on spin-1/2 qubits, and a single layer of even/odd gates repeated periodically in time. We find that this minimal model is not robustly thermalizing at numerically accessible system sizes, and displays slow subdiffusive dynamics for long times. We map out the thermalization dynamics in a broader parameter space of charge conserving circuits, and understand the origin of the slow dynamics in terms of proximate localized and integrable regimes in parameter space. In contrast, we find that small extensions to the minimal model are sufficient to achieve robust thermalization; these include (i) increasing the interaction range to three-site gates, (ii) increasing the local Hilbert space dimension by appending an additional unconstrained qubit to the conserved charge on each site, or (iii) using a larger Floquet period comprised of two independent layers of gates. Furthermore, our results should inform future studies of charge conserving circuits, which are relevant for a wide range of topical theoretical questions.

1-dimensional spin chains↗
DOE OSTI2024

The Effect of the Prior and the Experimental Design on the Inference of the Precision Matrix in Gaussian Chain Graph Models

Here, we investigate whether (and how) experimental design could aid in the estimation of the precision matrix in a Gaussian chain graph model, especially the interplay between the design, the effect of the experiment and prior knowledge about the effect. Estimation of the precision matrix is a fundamental task to infer biological graphical structures like microbial networks. We compare the marginal posterior precision of the precision matrix under four priors: flat, conjugate Normal-Wishart, Normal-MGIG and a general independent. Under the flat and conjugate priors, the Laplace-approximated posterior precision is not a function of the design matrix rendering useless any efforts to find an optimal experimental design to infer the precision matrix. In contrast, the Normal-MGIG and general independent priors do allow for the search of optimal experimental designs, yet there is a sharp upper bound on the information that can be extracted from a given experiment. We confirm our theoretical findings via a simulation study comparing (i) the KL divergence between prior and posterior and (ii) the Stein’s loss difference of MAPs between random and no experiment. Our findings provide practical advice for domain scientists conducting experiments to better infer the precision matrix as a representation of a biological network.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Xylan-cellulose core structure of oat water-extractable β-glucan macromolecule: Insight into interactions and organization of the cell wall complex

Water-extractable β-glucan with high molar mass (HM) determines health benefits of oat food. Oat β-glucan was extracted by a standardized in vitro digestion method and co-existing water-extractable polysaccharide (WEP) fraction and its HM-arabinoxylan (HM-AX) subfraction were isolated to identify their highly acid-resistant subunit and investigate molecular interactions between constituent polymers. The WEP and HM-AX samples consisted of arabinoxylans (AXs) (74 and 76 %, respectively), however, cellulose constituted the secondary component (6.6 and 12.8 %, respectively). Multi-detection HPSEC along with specific enzymatic hydrolysis of AXs revealed the presence of the HM-xylan domain (16 and 34 %, respectively) built of numerous single- and multi-component populations with random coil and rod-like conformations, which were embedded in a xylan matrix with spherical conformation and controlled the macromolecular shape. Unlike single-component populations, the multi-component ones were resistant to hydrolytic action of AX-hydrolyzing enzymes and represented the subunits that anchor matrix polysaccharides onto cellulose surface. These results indicate that water-extractable β-glucan macromolecule comprises as integral element a cellulose core with two linking populations, HM-xylan and low molar mass glucomannan, which are surrounded by a feruloylated AX-arabinan-arabinogalactan composite and next laminated by β-glucan matrix. Here, the stiff cellulose-xylan backbone is the basis of HM β-glucan organization, controlled by its cellulose-like segments.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

State densities of heavy nuclei in the static-path plus random-phase approximation

We report that nuclear state densities are important inputs to statistical models of compound-nucleus reactions. State densities are often calculated with self-consistent mean-field approximations that do not include important correlations and must be augmented with empirical collective enhancement factors. Here, we benchmark the static-path plus random-phase approximation (SPA + RPA) to the state density in a chain of samarium isotopes 148–155 Sm against exact results (up to statistical errors) obtained with the shell-model Monte Carlo (SMMC) method. The SPA + RPA method incorporates all static fluctuations beyond the mean field together with small-amplitude quantal fluctuations around each static fluctuation. Using a pairing plus quadrupole interaction, we show that the SPA + RPA state densities agree well with the exact SMMC densities for both the even- and odd-mass isotopes. For the even-mass isotopes, we also compare our results with mean-field state densities calculated with the finite-temperature Hartree-Fock-Bogoliubov (HFB) approximation. We find that the SPA + RPA repairs the deficiencies of the mean-field approximation associated with broken rotational symmetry in deformed nuclei and with the violation of particle-number conservation in the pairing condensate. In particular, in deformed nuclei the SPA + RPA reproduces the rotational enhancement of the state density relative to the mean-field state density.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Damage and Failure in a Statistical Crack Model

In this paper I will take a close look at a statistical crack model (SCM) as is used in engineering computer codes to simulate fracture at high strain rates. My general goal is to understand the macroscopic behavior effected by the microphysical processes incorporated into an SCM. More specifically, I will assess the importance of including local interactions between cracks into the growth laws of an SCM. My strategy will be to construct a numerical laboratory that represents a single computational cell containing a realization of a statistical distribution of cracks. The cracks will evolve by the microphysical models of the SCM, leading to quantifiable damage and failure of the computational cell. I will use the numerical data generated by randomly generated ensembles of the fracture process to establish scaling laws that will modify and simplify the implementation of the SCM into large scale engineering codes.

42 ENGINEERING↗
NASA NTRS1985

A review of critical in-flight events research methodology

Pilot's cognitive responses to critical in-flight events (CIFE's) were investigated, using pilots, who had on the average about 2540 flight hours each, in four experiments: (1) full-mission simulation in a general aviation trainer, (2) paper and pencil CIFE tests, (3) interactive computer-aided scenario testing, and (4) verbal protocols in fault diagnosis tasks. The results of both computer and paper and pencil tests showed only 50 percent efficiency in correct diagnosis of critical events. The efficiency in arriving at a diagnosis was also low: over 20 inquiries were made for 21 percent of the scenarios diagnosed. The information-seeking pattern was random, with frequent retracing over old inquiries. The measures for developing improved cognitive skills for CIFE's are discussed.

Giffin, W. C.↗
DOE OSTI2021

Optimized molecular aggregation via incorporating fluorinated unit in the polymer donor for 17.3% efficiency organic solar cells

Regulating molecular aggregation state via chemical modification of photovoltaic materials to optimize the blend morphology is an effective and challenging strategy to improve the photovoltaic performance of organic solar cells (OSCs). Here, in this work, we demonstrate a random terpolymerization strategy by incorporating 1,3-bis(5-bromo-4-fluorothiophen-2-yl)-5,7-bis(2-ethylhexyl)-4H,8H-benzo[1,2-c:4,5-c']bisthiophene-4,8-dione (BFT, 30 mol%) unit with fluorinated π-bridges into the molecular backbone of PM6 and develop a terpolymer PF1. Benefiting from the strong electronegativity of fluorine atom and noncovalent interaction of F---S, the energy levels, molecular ordering, and aggregation properties of the polymer donors are precisely regulated. Notably, in the active layer film based on PF1: Y6, the π-π stacking distance can reduce to 3.51 Å, which is significantly smaller than that of PM6:Y6 (3.65 Å) and boosts the effective exciton dissociation and charge transport, thereby effectively optimize the photovoltaic properties of the related devices. As a result, the PF1:Y6-based OSCs achieved a power conversion efficiency (PCE) of up to 17.3%, which is much higher than that of the controlled OSCs based on PM6:Y6 (PCE = 16.2%). These results demonstrate high-efficiency OSCs enabled by regulating molecular aggregation via terpolymerization strategy.

14 SOLAR ENERGY↗
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↗
DOE OSTI2024

Transferring predictions of formation energy across lattices of increasing size*

In this study, we show the transferability of graph convolutional neural network (GCNN) predictions of the formation energy of the nickel-platinum solid solution alloy across atomic structures of increasing sizes. The original dataset was generated with the large-scale atomic/molecular massively parallel simulator using the second nearest-neighbor modified embedded-atom method empirical interatomic potential. Geometry optimization was performed on the initially randomly generated face centered cubic crystal structures and the formation energy has been calculated at each step of the geometry optimization, with configurations spanning the whole compositional range. Using data from various steps of the geometry optimization, we first trained our open-source, scalable implementation of GCNN called HydraGNN on a lattice of 256 atoms, which accounts well for the short-range interactions. Using this data, we predicted the formation energy for lattices of 864 atoms and 2048 atoms, which resulted in lower-than-expected accuracy due to the long-range interactions present in these larger lattices. We accounted for the long-range interactions by including a small amount of training data representative for those two larger sizes, whereupon the predictions of HydraGNN scaled linearly with the size of the lattice. Therefore, our strategy ensured scalability while reducing significantly the computational cost of training on larger lattice sizes.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Development of a Random-Forest Cloud-Regime Classification Model Based on Surface Radiation and Cloud Products

Various methods have been developed to characterize cloud type, otherwise referred to as cloud regime. These include manual sky observations, combining radiative and cloud vertical properties observed from satellite, surface-based remote sensing, and digital processing of sky imagers. While each method has inherent advantages and disadvantages, none of these cloud-typing methods actually includes measurements of surface shortwave or longwave radiative fluxes. Here, a method that relies upon detailed, surface-based radiation and cloud measurements and derived data products to train a random-forest machine-learning cloud classification model is introduced. Measurements from five years of data from the ARM Southern Great Plains site were compiled to train and independently evaluate the model classification performance. A cloud-type accuracy of approximately 80% using the random-forest classifier reveals that the model is well suited to predict climatological cloud properties. Furthermore, an analysis of the cloud-type misclassifications is performed. While physical cloud types may be misreported, the shortwave radiative signatures are similar between misclassified cloud types. From this, we assert that the cloud-regime model has the capacity to successfully differentiate clouds with comparable cloud–radiative interactions. Therefore, we conclude that the model can provide useful cloud-property information for fundamental cloud studies, inform renewable energy studies, and be a tool for numerical model evaluation and parameterization improvement, among many other applications.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2010

Hidden Statistics Approach to Quantum Simulations

Recent advances in quantum information theory have inspired an explosion of interest in new quantum algorithms for solving hard computational (quantum and non-quantum) problems. The basic principle of quantum computation is that the quantum properties can be used to represent structure data, and that quantum mechanisms can be devised and built to perform operations with this data. Three basic non-classical properties of quantum mechanics superposition, entanglement, and direct-product decomposability were main reasons for optimism about capabilities of quantum computers that promised simultaneous processing of large massifs of highly correlated data. Unfortunately, these advantages of quantum mechanics came with a high price. One major problem is keeping the components of the computer in a coherent state, as the slightest interaction with the external world would cause the system to decohere. That is why the hardware implementation of a quantum computer is still unsolved. The basic idea of this work is to create a new kind of dynamical system that would preserve the main three properties of quantum physics superposition, entanglement, and direct-product decomposability while allowing one to measure its state variables using classical methods. In other words, such a system would reinforce the advantages and minimize limitations of both quantum and classical aspects. Based upon a concept of hidden statistics, a new kind of dynamical system for simulation of Schroedinger equation is proposed. The system represents a modified Madelung version of Schroedinger equation. It preserves superposition, entanglement, and direct-product decomposability while allowing one to measure its state variables using classical methods. Such an optimal combination of characteristics is a perfect match for simulating quantum systems. The model includes a transitional component of quantum potential (that has been overlooked in previous treatment of the Madelung equation). The role of the transitional potential is to provide a jump from a deterministic state to a random state with prescribed probability density. This jump is triggered by blowup instability due to violation of Lipschitz condition generated by the quantum potential. As a result, the dynamics attains quantum properties on a classical scale. The model can be implemented physically as an analog VLSI-based (very-large-scale integration-based) computer, or numerically on a digital computer. This work opens a way of developing fundamentally new algorithms for quantum simulations of exponentially complex problems that expand NASA capabilities in conducting space activities. It has been illustrated that the complexity of simulations of particle interaction can be reduced from an exponential one to a polynomial one.

Zak, Michail↗
DOE OSTI2022

A new multi-model absolute difference-based sensitivity (MMADS) analysis method to screen non-influential processes under process model and parametric uncertainty

Process-based models have been widely used for hydrologic modeling, and it is a common practice to use sensitivity analysis methods for excluding non-influential hydrologic processes from further investigation and/or model improvement. This study develops a new method called multi-model absolute difference-based sensitivity (MMADS) analysis method to screen non-influential system processes and parameters. MMADS is conceptually similar to the Morris method for addressing parametric uncertainty, but has a unique feature to address both process model uncertainty (i.e., a process may be represented by multiple process models) and process model parameter uncertainty (i.e., parameters associated with a process model are random). MMADS first evaluates absolute differences of a quantity of interest (i.e., a system model output) by varying process models and/or process model parameter values, and then calculates the mean and variance of the differences for investigating process influence. The mean measures overall influence of the process on the quantity of interest, and the variance estimates influence of nonlinear effects of the process and/or its interactions with other processes. MMADS is an extension of the Morris method from a parameter space to a joint parameter-model space for explicitly addressing both process model uncertainty and model parameter uncertainty. The performance of MMADS is evaluated by using two numerical experiments. One experiment is based on Sobol’s G*-function with ten product elements, and has analytical solutions of the MMADS mean and variance of absolute differences. The other experiment is for groundwater flow modeling which considers three processes (i.e., recharge, geology, and snowmelt) that interact with each other. Finally, results indicate that MMADS is computationally efficient and can identify non-influential processes of complex hydrological systems.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Prediction of vacancy defect diffusion paths in high entropy alloys via machine learning on molecular dynamics data

Identifying the diffusion path of point defects is a critical step in understanding their evolution and the mechanisms of related phenomena. Defect diffusion occurs at small length and time scales, with impacts on material properties that may continue to evolve over ns to μs, ms, and the continuum scale (s, min, etc., and cm, m, etc.). The time scale accessible to molecular dynamics (MD) simulations is limited by small step sizes, typically in the fs range. Thus, surrogate models of MD simulations through machine learning (ML)-based algorithms are of great interest, especially for complex systems such as high entropy alloys (HEAs). In this work, dynamics governing vacancy migration in HEA were approximated with graph convolutional network (GCN) models as ansatzes for kinetic Monte Carlo (KMC) rate catalogs. Network design considered that diffusion in crystalline solids generally depends on interactions between defects and their immediate neighbor atoms. Graphs represented the vacancy surroundings, MD-generated trajectories provided training and comparison datasets, and unsupervised GCN models approximated interatomic dynamics governing vacancy migration in HEAs as ansatzes for KMC. A proof-of-concept model trained on MD data for the Fe, Ni, Cr, Co, and Cu HEA environment was used with two different neighbor interactions to assess the feasibility of training a GCN to predict vacancy defect transition rates in the HEA environment. The resulting setup rapidly generated MD-formatted synthetic trajectories based on dynamics learned from the MD training set, with a time acceleration of roughly two orders of magnitude and a similar diffusion coefficient to MD observations. Additionally, Nudged Elastic Band (NEB) calculations were performed on randomly generated FeNiCrCoCu HEA structures to determine vacancy migration barriers across nearest-neighbor sites. Transition probabilities for each jump, categorized by atomic type, were extracted from these calculations. NEB-based and GCN-based approaches led to similar outcomes.

Reimer, C↗
NASA NTRS1985

Fundamentals of microcrack nucleation mechanics

A foundation for ultrasonic evaluation of microcrack nucleation mechanics is identified in order to establish a basis for correlations between plane strain fracture toughness and ultrasonic factors through the interaction of elastic waves with material microstructures. Since microcracking is the origin of (brittle) fracture, it is appropriate to consider the role of stress waves in the dynamics of microcracking. Therefore, the following topics are discussed: (1) microstress distributions with typical microstructural defects located in the stress field; (2) elastic wave scattering from various idealized defects; and (3) dynamic effective-properties of media with randomly distributed inhomogeneities.

Fu, L. S.↗