Search NASA⌕ Search

SEARCH · Search NASA

Results for “Conditional distribution”

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 415 records · Page 23

Performance Evaluation of Multi-Vendor Grid-Forming Inverters for Grid-Connected Operation Through Hardware Experimentation: Preprint

Existing real-world projects of GFM inverters that operate in parallel to power grids typically are sized between dozens and a few hundreds Megawatt (MW) scale according to a recent NERC GFM inverter white paper. These large systems are often difficult to evaluate prior to deployments because of their large size. The performance of smaller GFM inverters (dozens to a few hundreds MW) that operate parallel with power grids (distribution systems) is even less understood. There is an opportunity to better understand these systems through hardware testing under controlled laboratory conditions. Therefore, this paper presents the functional performance evaluation tests of multiple (three) commercial GFM inverters when they operate in parallel with the grid through hardware experiments. The goal of these tests is to explore the GFM inverters' functionalities and dynamic response when in parallel with power grids to eventually develop universal specifications for GFM inverters. Both steady state (changing the inverter's frequency and voltage droop) and transient (adding step change in grid's frequency/voltage) tests are performed for each GFM inverter with the same testing circuit and testing protocol. The experimental results indicate the bench-marked performance that: 1) the GFM inverters can be dispatched through frequency and voltage droop intercepts to output the target power when paralleled to the grid; 2) the GFM inverters automatically respond to system frequency and voltage events to output the needed power, however, the GFM inverters all show stability issues when absorbing reactive power from the grid.

grid-forming inverters↗

3D TRISO particle-explicit compact meshing

The TRI-structural ISOtropic (TRISO) layered fuel particle is a robust nuclear fuel form offering enhanced safety and performance for advanced reactor concepts, including high-temperature gas-cooled reactors and other Generation IV designs. These poppy-seed-sized particles are embedded in a graphite matrix to form fuel elements that must withstand elevated temperatures and high burn-up levels. The heterogeneous nature of these fuel elements — comprising thousands of randomly distributed TRISO particles — produces complex stress fields and thermal gradients that one- and two-dimensional models cannot accurately capture. While three-dimensional modeling has improved predictions of dimensional changes, internal pressure buildup, and fission product transport under irradiation, current approaches rely on homogenized material properties that are known to have considerable divergence from experimental observations. This work presents a methodology for optimized random packing of TRISO fuel compacts and full three-dimensional mesh generation within the BISON fuel performance code, with each particle coating layer individually discretized. The resulting mesh was demonstrated through heat conduction simulations under representative in-reactor operating conditions, showing strong agreement with expected behavior. This capability enables detailed analysis of particle-to-particle interactions, matrix cracking mechanisms, and the statistical distribution of coating layer failures — all of which directly govern fuel performance and safety margins.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Experimental and Theoretical Evaluation of Feed Flow Collar Design for Shell Fed Hollow Fiber Membrane Modules

An experimental and theoretical study of module collar design is presented here. Hollow fiber membranes are prepared by dip coating a poly(vinylidene) (PVDF) support with a polydimethylsiloxane (PDMS) gutter layer and a Pebax 2533 selective layer. Fiber bundles with a well-defined fiber packing are prepared using a 3D printed module. A parallel fiber bundle consisting of 4-9 uniformly spaced fibers is created with printed tabs that align the fibers and create a tubesheet. The tabs are sealed within a printed case that possesses a series of external ports for gas introduction and removal. Uniquely, both port location and the use of a collar to assist fluid distribution in the shell can be varied for the same fiber bundle. Experimental measurements are compared to computational fluid dynamics (CFD) simulations. The experimental module design allows high-fidelity representation of the fiber bundle and module case in the simulations. Comparisons between experiment and simulation are in good agreement over a broad range of experimental conditions. The detrimental effect of having ports located too close, leading to stagnation regions, is captured as well as the beneficial effects of using a collar for shell-side fluid distribution around the fiber bundle. Such results help validate the use of CFD to develop high-performance module designs.

Tran, Thien↗

Full time-dependent SOLPS-ITER simulation of the SPARC tokamak: actuator design for particle and divertor condition control *

Abstract This paper presents the application of full time-dependent SOLPS-ITER simulations for actuator design in the SPARC tokamak. This study employs both the EIRENE module, a neutral solver, and the B2.5 plasma module in a time-dependent mode. This is in contrast to most SOLPS simulations, which focus on steady-state solutions, where the neutral distribution is evolved without any time limit or for a time step of 1 ⋅ 10 − 3 second, which is several orders of magnitude larger than the fluid plasma time step. The time-dependent EIRENE was tested with a fixed B2.5 background and compared with a simple conductance based model in a simplified pump chamber geometry. This comparison aimed to verify the reliability of the neutral relaxation timescale derived from the time-dependent EIRENE. Subsequently, a full time-dependent simulation was performed in a realistic geometry, with the Monte-Carlo neutral time step synchronized with the plasma fluid time step. The numerical setup of the code, including relative time steps and the size of the census data used to store Monte-Carlo particle information is considered. The full-time dependent simulations are then applied to inform the design of the SPARC louver structure, which affects divertor plasma parameters by regulating the neutral conductance from the divertor to the pump. The response of the plasma and neutral parameters was captured on a timescale that enables the design of the actuator to consider time-dependent control capability. It was found that changing the louver opacity has an equivalent effect as varying the gas throughput via puff actuation. Therefore, equivalent divertor plasma conditions can be obtained from both actuators, while the neutral pressure distribution in the pump and divertor differs for each actuator.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Habitat specialization and edge effects of soil microbial communities in a fragmented landscape

Abstract Soil microorganisms play outsized roles in nutrient cycling, plant health, and climate regulation. Despite their importance, we have a limited understanding of how soil microbes are affected by habitat fragmentation, including their responses to conditions at fragment edges, or “edge effects.” To understand the responses of soil communities to edge effects, we analyzed the distributions of soil bacteria, archaea, and fungi in an experimentally fragmented system of open patches embedded within a forest matrix. In addition, we identified taxa that consistently differed among patch, edge, or matrix habitats (“specialists”) and taxa that showed no habitat preference (“nonspecialists”). We hypothesized that microbial community turnover would be most pronounced at the edge between habitats. We also hypothesized that specialist fungi would be more likely to be mycorrhizal than nonspecialist fungi because mycorrhizae should be affected more by different plant hosts among habitats, whereas specialist prokaryotes would have smaller genomes (indicating reduced metabolic versatility) and be less likely to be able to sporulate than nonspecialist prokaryotes. Across all replicate sites, the matrix and patch soils harbored distinct microbial communities. However, sites where the contrasts in vegetation and pH between the patch and matrix were most pronounced exhibited larger differences between patch and matrix communities and tended to have edge communities that differed from those in the patch and forest. There were similar numbers of patch and matrix specialists, but very few edge specialist taxa. Acidobacteria and ectomycorrhizae were more likely to be forest specialists, while Chloroflexi, Ascomycota, and Glomeromycota (i.e., arbuscular mycorrhizae) were more likely to be patch specialists. Contrary to our hypotheses, nonspecialist bacteria were not more likely than specialist bacteria to have larger genomes or to be spore‐formers. We found partial support for our mycorrhizal hypothesis: arbuscular mycorrhizae, but not ectomycorrhizae, were more likely to be specialists. Overall, our results indicate that soil microbial communities are sensitive to edges, but not all taxa are equally affected, with arbuscular mycorrhizae in particular showing a strong response to habitat edges. In the context of increasing habitat fragmentation worldwide, our results can help inform efforts to maintain the structure and functioning of the soil microbiome.

Winfrey, Claire C. [Department of Ecology and Evol↗

Computing material volume fractions on a superimposed mesh as applied to Monte Carlo particle transport simulations

Here, we present a newly implemented ray tracing algorithm in OpenMC for efficiently computing material volume fractions on superimposed meshes in complex geometries. By firing rays along each coordinate direction through the geometry, the approach accumulates track-length data in each mesh element, thereby determining the fractional composition of each material. Scaling studies on three different models—a random tetrahedra configuration, the Frascati Neutron Generator ITER dose rate benchmark, and a stellarator design—show excellent parallel performance, with nearly linear speedup on modern multi-threaded and distributed-memory systems. An analysis of the residual error relative to high-resolution reference solutions demonstrated that under optimal conditions it decreases as 1/R, where R is the number of rays fired, making it straightforward to achieve user-prescribed accuracy. This new functionality enables practical, mesh-based approaches for detailed nuclear analyses in production Monte Carlo workflows without resorting to expensive, fully conformal or unstructured meshing.

Monte Carlo↗

Assessing the behavioral realism of energy system models in light of the consumer adoption literature

Effective policymaking to achieve net zero greenhouse gas emissions demands an understanding of the complex drivers of, and barriers to, consumer adoption behavior via behaviorally realistic energy system models. Existing models tend to oversimplify by focusing on homogenized financial factors while neglecting consumer heterogeneity and non-monetary influences. This study develops and applies a comprehensive framework for evaluating the behavioral realism of consumer adoption models, informed by the adoption literature. It introduces a typology for factors influencing low-carbon technology adoption decisions: monetary and non-monetary factors relating to household characteristics, psychology, technological attributes, and contextual conditions. Next, reviews of the consumer adoption and decision-making literature identify the most influential adoption factor categories for distributed solar photovoltaics, electric vehicles, and air-source heat pumps. Finally, the extent to which a selection of energy system models accounts for these adoption factors is assessed. Existing models predominantly emphasize the economic aspects of technology, which are generally identified as the most important factors. Where the models fall short — in considering moderately important factor categories — sector-specific and agent-based models can offer more behaviorally realistic insights. This study sheds light on which types of factors are most important for consumer adoption decisions and investigates how well current models rise to the challenge of behavioral realism. The end-to-end analysis presented enables internally consistent comparisons across models and energy technologies. This research advances timely conversations on consumer adoption. It could inform more behaviorally realistic energy system modeling, and thereby more effective decarbonization policymaking.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

CFD Simulation of the Dosing Behavior within the Atomic Layer Deposition Feeding System

The effective operation of atomic layer deposition (ALD) feeding system is the premise of realizing specific ALD processes. In the present work, a detailed computational fluid dynamics (CFD) model of the feeding system has been developed and validated, which accounts for the roles of ALD valves and manifolds. A numerical simulation of the compressible fluid flow and heat/mass transfer within the feeding system was conducted. The dosing amounts and the spatiotemporal distributions of the precursors can be accurately predicted using the CFD model, as validated by experimental results. Different precursors, operating conditions, and structures of the feeding system were simulated and analyzed to examine the operating flexibility of the feeding system. The simulation results can be adopted as the upstream boundary conditions for simulations of the ALD process in the reaction chamber. The substrate-scale simulation indicates that the effect of the feeding system on the film deposition is highly related to the surface kinetics of ALD. The present work can serve as a guide for the development and optimization of different ALD-based processes via proper operation and even the design of the feeding system.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

U redox state tracked in mineralized hydrothermal carbonate with implications for U-Pb geochronology

U-Pb carbonate geochronology can directly constrain the timing and rates of important geological processes. However, the mechanisms and controls on U incorporation, distribution, and retention in carbonate minerals remain unclear, limiting geological interpretations. Here X-ray absorption spectroscopy (µXAS) and in-situ U-Pb carbonate geochronology are combined to temporally track U distribution and redox state in a porphyry-epithermal system. In this setting, multiple generations of carbonate minerals record fluid conditions and processes which control the solubility and deposition of metals, including U. This novel approach provides the first evidence of both oxidized UO 2 2+ and reduced U 4+ species in temporally distinct generations of carbonate within a single sample. Preservation of two different U oxidation states during discrete precipitation events requires U retentivity within older domains, demonstrating that the U-Pb carbonate geochronometer is robust under hydrothermal conditions. Furthermore, crystal zones with abundant fluid/vapour inclusions linked to boiling processes coincide with relatively high levels of U and favourable U/Pb. Targeting carbonate domains with these textures may therefore increase success in U-Pb geochronology. U-Pb carbonate dating combined with µXAS can track the temporal evolution of processes critical for metal deposition in long-lived and multistage hydrothermal-magmatic ore deposit settings.

58 GEOSCIENCES↗

Probabilistic Analysis of Long-Term Degradation of Microwave Cavity Flow Sensor

We are investigating a microwave resonant cavity transducer for flow sensing in the vessel of a high temperature fluid advanced reactor (AR), such as a molten salt cooled reactor (MSCR) or a sodium fast reactor (SFR). This transducer is a hollow metallic cylindrical cavity, with the flat wall of the cylinder flexible enough to undergo microscopic deflection due to dynamic fluid pressure. Membrane deflection leads to a shift in the resonant frequency, which can be detected with a spectrum analyzer. We have performed a proof-of-concept experiment of flow sensing with the transducer in liquid sodium at 340°C in impinging liquid jet geometry. The transducer remained in liquid sodium for 70 days. After removal, no structural damage was observed, and the expected transducer response was verified in a water test. Because long-term (multi-year) experimental tests of transducer resilience to harsh environment are not practical, we have developed a probabilistic model of creep to estimate transducer resilience to the harsh environment. The probabilistic model considers diffusion creep under the condition of high temperature and low stress, where the stress and temperature are allowed to be random variables with Gaussian distributions. Using the probabilistic model, we estimate inelastic membrane deflections due to creep for several temperature ranges. We conclude that for temperatures less than 650°C, creep has negligible long-term effect on the transducer performance. Since a yellowish residue was observed on the transducer surface after 70 days of immersion in liquid sodium, we have investigated possible evidence of corrosion. Chromium depletion is a typical indicator of the corrosion process in stainless steel. Scraping off a residue from the transducer and performing scanning electron microscopy (SEM) with energy dispersive analysis (EDS) did not find any chromium in the residue. Approximately 60% of the residue consisted of copper, which can be attributed to contamination of sodium due to powder residue from machining of copper and brass components of the transducer.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Advancing Grid Resilience through Smart Charge Management: Findings from Maryland’s Pilot

This report presents research findings from a four-year Smart Charge Management (SCM) pilot program conducted by Maryland’s largest electric utilities—Baltimore Gas and Electric (BGE), Potomac Electric Power Company (Pepco), and Delmarva Power & Light (DPL)—to evaluate strategies for optimizing electric vehicle (EV) charging loads and enhancing grid stability. Supported by the U.S. Department of Energy (DOE), Argonne National Laboratory collaborated with all project partners and examined the effectiveness of Time-of-Use (TOU) and Load Balancing (LB) strategies in managing peak demand, deferring costly infrastructure upgrades, and reducing grid constraints at the feeder level. Using charging data from over 4,600 EV drivers, the study analyzed SCM’s impact on the distribution systems of BGE and Pepco, which consists of over 2000 feeders. Unlike prior research that focused on system-wide trends or synthetic feeders, this analysis offers granular, feeder-level insights based on real-world operational data. It highlights how transformer density, load profiles, and infrastructure constraints influence smart charging performance. Results show feeder-level conditions play a crucial role in SCM effectiveness, with most feeders benefiting more from LB, while TOU-based SCM may be sufficient for others. By 2035, LB reduced peak charging loads by 27% on average, compared to 23% under TOU-based SCM, though some feeders saw reductions exceeding 35%, while others experienced minimal impact. Feeders with higher transformer utilization and limited capacity benefited more from LB, which more effectively distributed charging demand during off-peak hours. Beyond reducing grid constraints, SCM offers long-term operational and financial benefits. By shifting EV charging demand strategically, utilities can optimize asset utilization, delay infrastructure investments, and enhance grid performance. In terms of infrastructure upgrade deferrals, at the feeder level, LB consistently reduced peak charging loads and resulting infrastructure upgrade costs, particularly in high EV enrollment areas, decreasing the number of overloaded transformers by up to 35%, while TOU-based SCM achieved 20-30% reductions depending on feeder characteristics. At the system level, LB has the potential to defer total upgrade costs by $\$$186 million for BGE, compared to $\$$159 million under TOU-based SCM. For Pepco, TOU-based SCM performed slightly better, deferring upgrade costs by $\$$30 million, compared to $\$$29 million under LB. Section 4.5 reviews some of the system differences between BGE and Pepco. However, as EV adoption scales, TOU-based SCM will introduce secondary peak charging loads, reinforcing the need for more advanced, adaptive SCM approaches to prevent new grid challenges. As EV adoption continues to grow, feeder-level managed charging strategies will be essential for mitigating grid stress, improving infrastructure efficiency, and maintaining energy affordability for consumers. This report provides critical insights for utilities, Public Utility Commissions (PUCs), and state agencies on the role of feeder-specific smart charging in infrastructure planning, policy development, and grid modernization. The findings underscore the importance of tailored, data-driven SCM solutions that align with local grid conditions, ensuring a resilient, cost-effective transition to increasing EV adoption while safeguarding distribution system performance.

24 POWER TRANSMISSION AND DISTRIBUTION↗

A Physics-Informed Deep Learning Description of Knudsen Layer Reactivity Reduction

A physics-informed neural network (PINN) is used to evaluate the fast ion distribution in the hot spot of an inertial confinement fusion target. The use of tailored input and output layers to the neural network is shown to enable a PINN to learn the parametric solution to the Vlasov–Fokker–Planck equation in the absence of any synthetic or experimental data. As an explicit demonstration of the approach, the specific problem of Knudsen layer fusion yield reduction is treated. Here, the predictions from the Vlasov–Fokker–Planck PINN are used to provide a non-perturbative solution of the fast ion tail in the vicinity of the hot spot, thus allowing the spatial profile of the fusion reactivity to be evaluated for a range of collisionalities and hot spot conditions. Excellent agreement is found between the predictions of the Vlasov–Fokker–Planck PINN and the results from traditional numerical solvers with respect to both the energy and spatial distribution of fast ions and the fusion reactivity profile, demonstrating that the Vlasov–Fokker–Planck PINN provides an accurate and efficient means of determining the impact of Knudsen layer yield reduction across a broad range of plasma conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Dynamic Temporal Graph Sequence Data for Resilience-Oriented Distribution Network Reconfiguration

This dataset comprises temporal dynamic graph sequences generated from power grid simulations focused on grid reconfiguration to enhance resilience. The simulations model failure propagation under varying conditions, with nodes assigned distinct failure probabilities. For each time step, the dataset captures the evolution of node states (functional or failed) and features critical to grid operations, such as pv_output, load_profile, load_dispatch, dg_output, loss, and voltage. Node types include sources, normal loads, and nodes with specific equipment like PVs, micro turbines, or shunt capacitors. The dataset is structured to support the training of dynamic graph neural networks, facilitating research on node feature prediction and edge dynamics under failure scenarios. Three distinct configurations are included, providing a robust foundation for modeling power grid resilience.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Pyrodictium abyssi AbpX reveals a calcium-responsive family of microbial biomatrix proteins that form thermostable hydrogels

Evolutionary pressure on microbial communities propagating under extreme environmental conditions often results in unique structural adaptations to promote cell survival. In this work, we report an investigation of AbpX, a biomatrix protein identified in cultures of the hyperthermophilic archaeon Pyrodictium abyssi. Under ex vivo and in vitro conditions, AbpX assembles into a paracrystalline lattice composed of semiflexible fibrils. CryoEM analysis of recombinant AbpX fibrils reveals that the precursor protein polymerizes through donor strand complementation (DSC), a process previously reported for chaperone-usher fimbriae in Gram-negative bacteria. Unlike the latter DSC protein polymers, AbpX undergoes chaperone-free polymerization in the presence of calcium ions, which are sequestered at the donor strand-acceptor groove interface between protomers in the fibril. Using a combination of cryoEM and crystallographic information, a structural model is proposed for the AbpX lattice that provides insight into its potential role in biofilm formation. These findings suggest that calcium ion coordination may contribute to fibril assembly and preorganize fibrils for incorporation into the protein lattice. Bioinformatic analysis indicates that AbpX exemplifies a distinct and broadly distributed clade of calcium ion responsive biomatrix proteins within the TasA superfamily that can be fabricated into hydrogel biomaterials in vitro under environmentally benign conditions.

59 BASIC BIOLOGICAL SCIENCES↗

A periplasmic zinc capture protein enhances the resistance of Neisseria gonorrhoeae to nutritional immunity

During microbial infection, mammalian hosts reduce the availability of free metals such as zinc in a process known as nutritional immunity. Pathogens counteract nutritional immunity by expressing gene products that enhance growth in metal-limited conditions. One of the most transcriptionally induced genes in zinc-limited Neisseria gonorrhoeae , ngo1049, encodes a DUF4198 family protein we have named Zcp. This family of proteins is widely distributed in Gram-negative bacteria. Here, we provide the first structural, biochemical, and functional characterization of a DUF4198 protein. Zcp is a periplasmic, homodimeric substrate-binding protein (SBP), which binds one zinc ion per subunit with submicromolar affinity. We identified a zinc binding pocket in each subunit, composed of three histidine residues. Zcp enables maximal growth of N. gonorrhoeae in zinc-limited conditions but is dispensable for zinc uptake, in contrast to the cluster A-I SBP ZnuA, which is required for zinc import. The growth defect of zcp mutant N. gonorrhoeae is rescued by zinc supplementation. Zcp associates with proteins with roles in maintaining cell envelope integrity, and N. gonorrhoeae lacking zcp is more sensitive to envelope-targeting antimicrobials. Zcp enables infectivity of human epithelial cells and neutrophils by zinc-limited N. gonorrhoeae . We conclude that N. gonorrhoeae produces Zcp to buffer periplasmic zinc, which enables ZnuA to balance import of different metals and ensures the bioavailability of zinc for extracytoplasmic zinc-requiring proteins, as part of the coordinated response to host-imposed nutritional immunity.

Liyayi, Ian K. [Department of Microbiology, Immuno↗

Integrated Transmission and Distribution Co-Simulation Platform for Demonstration of Bulk Grid Services Using Distributed Energy Resources

In September 2020, the Federal Energy Regulatory Commission (FERC) released Order 2222, which opens wholesale markets to small-capacity distributed energy resources (DERs), recognizing their potential in improving operational efficiency by providing bulk grid services. Therefore, a co-simulation capability that can connect transmission and distribution (T&D) simulations and evaluate the impacts of DER provision of bulk grid services is needed. In this paper, we present a new integrated T&D co-simulation platform that incorporates T&D system simulators, DER aggregator/group strategies, and a co-simulation coordinator. Industry-standard communication protocols are employed to mimic real-world conditions. Secondary frequency regulation is selected as the representative bulk grid service, and we simulate the responses of DERs to the frequency regulation signals. The simulation results for a solar-rich distribution feeder in Colorado, USA, demonstrate how the T&D co-simulation setup is used to evaluate the contributions of DERs to minimize the bulk grid frequency deviation.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Geologic Seawater Air Conditioning (GeoSWAC) System: Resource Assessment and Techno-Economic Evaluation in Puerto Rico

Puerto Rico's hot, humid climate drives a high and persistent demand for cooling, straining an aging and fuel-dependent energy grid while increasing peak electricity loads. Much of this challenge stems from inefficient air-conditioning systems operating in buildings without passive cooling design, making cooling both costly and vulnerable to disruption - especially during hurricanes. To address these issues, researchers evaluated a new alternative called Geologic Seawater Air Conditioning (GeoSWAC), which uses inland wells connected to naturally cold, deep seawater. By avoiding long offshore pipelines and energy-intensive refrigeration cycles, GeoSWAC can cut electricity use by 75-90%, reduce environmental impacts, and operate more reliably during power outages. A case study at the University of Puerto Rico's Rio Piedras campus found that GeoSWAC could deliver significantly lower levelized costs of cooling, reduced operational expenses, and improved water conservation compared to the campus's existing chilled-water plant. The system offers key advantages - including access to a constant cold heat sink, low pumping requirements, and enhanced resilience - for coastal regions such as Puerto Rico, the broader Caribbean, and Florida. However, its effectiveness depends heavily on site-specific geological conditions, such as subsurface connectivity and aquifer characteristics, which require detailed investigation. Future research should refine hydrogeological models, conduct performance and economic analyses, and address regulatory and permitting frameworks to support broader adoption of this promising technology.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Best of both worlds: Enforcing detailed balance in machine learning models of transition rates

The slow microstructural evolution of materials often plays a key role in determining material properties. When the unit steps of the evolution process are slow, direct simulation approaches such as molecular dynamics become prohibitive and Kinetic Monte-Carlo (kMC) algorithms, where the state-to-state evolution of the system is represented in terms of a continuous-time Markov chain, are instead frequently relied upon to efficiently predict long-time evolution. The accuracy of kMC simulations however relies on the complete and accurate knowledge of reaction pathways and corresponding kinetics. This requirement becomes extremely stringent in complex systems such as concentrated alloys where the astronomical number of local atomic configurations makes the a priori tabulation of all possible transitions impractical. Machine learning models of transition kinetics have been used to mitigate this problem by enabling the efficient on-the-fly prediction of kinetic parameters. While conventional KMC methods based on transition state theory naturally yield reversible dynamics that exactly obey the detailed balance criterion, providing strong guarantees on the properties of the stationary distribution, many recently-proposed ML-based approaches to barrier predictions provide no such guarantees. In this study, we derive conditions under which physics-informed ML architectures exactly enforce the detailed balance condition by construction, even when relying on non-extensive descriptions of states in terms of local environments around mobile defects. In conclusion, using the diffusion of a vacancy in a concentrated alloy as an example, we show that such ML architectures also exhibit superior performance in terms of prediction accuracy, demonstrating that the imposition of physical constraints can facilitate the accurate learning of barriers at no increase in computational cost.

36 MATERIALS SCIENCE↗