Search NASA⌕ Search

SEARCH · Search NASA

Results for “potential”

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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

The Realistic Potential of Soil Carbon Sequestration in U.S. Croplands for Climate Mitigation

Existing estimates of the climate mitigation potential from cropland carbon sequestration (C-sequestration) are limited because they tend to assume constant rates of soil organic carbon change over all available cropland area, use relatively coarse land delineations, and often fail to adequately consider the agronomic and socioeconomic dimensions of agricultural land use. This results in an inflated estimate of the C-sequestration potential. We address this gap by defining a more appropriate land base for cover cropping in the United States for C-sequestration purposes: stable croplands in annual production systems that can integrate cover cropping without irrigation. Our baseline estimate of this suitable stable cropland area is 32% of current U.S. cropland extent. Even an alternative, less restrictive definition of stability results in a large reduction in area (44% of current U.S. croplands). Focusing cover crop implementation to this constrained land base would increase durability of associated C-sequestration and limit soil carbon loss from land conversion to qualify for carbon-specific incentives. Applying spatially-variable C-sequestration rates from the literature to our baseline area yields a technical potential of 19.4 Tg CO 2 e yr –1 annually, about one-fifth of previous estimates. We also find the cost of realizing about half (10 Tg CO 2 e yr –1 ) of this potential could exceed 100 USD Mg CO 2 e –1 , an order of magnitude higher than previously thought. While our economic analyses suggest that financial incentives are necessary for large-scale adoption of cover cropping in the U.S., they also imply any C-sequestration realized under such incentives is likely to be additional.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Atomic cluster expansion potential for large scale simulations of hydrocarbons under shock compression

We present an Atomic Cluster Expansion (ACE) machine learned potential developed for high-fidelity atomistic simulations of hydrocarbons, targeting pressures and temperatures near and above supercritical fluid regimes for molecular fluids. A diverse set of stoichiometries were covered in training, including 1:0 (pure carbon), 1:4 (methane), and 1:1 (benzene), and rich bonding environments sampled at supercritical temperatures, hydrogen rich, reactive mixtures where metastable stoichiometries arise, including 1:2 (ethylene) and 1:3 (ethane). A high-fidelity training database was constructed by performing large-scale quantum molecular dynamic simulations [density functional theory (DFT) MD] of diamond, graphite, methane, and benzene. A novel approach to selecting structures from DFT MD is also presented, which allows for the rapid selection of unique DFT MD frames from complex trajectories. Comparisons to DFT and experimental data demonstrate that the presented ACE potential accurately reproduces isotherms, carbon melting curves, radial distribution functions, and shock Hugoniots for carbon and hydrocarbon systems for pressures up to 100 GPa and temperatures up to 6000 K for hydrocarbon systems and up to 9000 K for pure carbon systems. This work delivers a potential that can be used for accurate, large-scale simulations of shocked hydrocarbons and demonstrates a methodology for fitting and validating machine learning interatomic potentials to complex molecular environments, which can be applied to energetic materials in future works.

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

Benchmarking machine learning interatomic potentials via phonon anharmonicity

Abstract Machine learning approaches have recently emerged as powerful tools to probe structure-property relationships in crystals and molecules. Specifically, machine learning interatomic potentials (MLIPs) can accurately reproduce first-principles data at a cost similar to that of conventional interatomic potential approaches. While MLIPs have been extensively tested across various classes of materials and molecules, a clear characterization of the anharmonic terms encoded in the MLIPs is lacking. Here, we benchmark popular MLIPs using the anharmonic vibrational Hamiltonian of ThO 2 in the fluorite crystal structure, which was constructed from density functional theory (DFT) using our highly accurate and efficient irreducible derivative methods. The anharmonic Hamiltonian was used to generate molecular dynamics (MD) trajectories, which were used to train three classes of MLIPs: Gaussian approximation potentials, artificial neural networks (ANN), and graph neural networks (GNN). The results were assessed by directly comparing phonons and their interactions, as well as phonon linewidths, phonon lineshifts, and thermal conductivity. The models were also trained on a DFT MD dataset, demonstrating good agreement up to fifth-order for the ANN and GNN. Our analysis demonstrates that MLIPs have great potential for accurately characterizing anharmonicity in materials systems at a fraction of the cost of conventional first principles-based approaches.

interatomic potentials↗
DOE OSTI2025

Deep potential molecular dynamics simulations of ion-enhanced etching of silicon by atomic chlorine

The continued development of plasma-assisted processing techniques requires a fundamental understanding of plasma-surface interactions. Molecular dynamics (MD) simulations have been employed to complement experimental studies and better understand the properties of such systems. Recently, machine learning (ML) methods have enabled the development of ab initio-based interatomic potentials, which can be generalized to complex combinations of multiple atom types. In this work, we use ML potentials developed using the Deep Potential Molecular Dynamics (DeepMD) framework to provide a model of ion-enhanced etching of Si by Cl atoms. We demonstrate the importance of proper selection of the training data set to the accuracy of the DeepMD model and compare our results to MD results using empirical potentials, as well as to experimental measurements. Exposure of undoped Si at 300 K to thermal Cl atoms yields a steady-state Cl coverage of 1.25 monolayers, which is slightly lower than the value obtained in previous experimental studies. Predictions of Si etch yields by simultaneous Cl atom and Ar + ion impacts as a function of ion energy, neutral to ion flux ratio, and angle of incidence of the ions are in reasonably good agreement with classical MD results and experimental measurements. Finally, etch yields and SiCl x mixed layer thicknesses during simultaneous bombardment of the Si(100) surface by Cl atoms and Cl + ions are in good agreement with experimental data. In conclusion, the present work is a necessary condition for the extension of the DeepMD procedure to more complex systems of interest in plasma-surface interactions.

Artificial neural networks↗
DOE OSTI2026

Surface water nitrogen and sediment potential nitrate reduction rates, nutrient stocks, and stable isotopes from nine wetlands at the Tanglewood Biological Station, Alabama

This dataset supports a broader study investigating wetland hydrologic and biogeochemical responses to inundation disturbances. Bimonthly surface water and sediment sampling events were conducted at nine wetland sites situated within the Tanglewood Biological Station in Alabama from April 2023 to February 2024. The contents included in the data package include surface water nitrogen (nitrogen oxides and ammonium) and sediment potential nitrate reduction rates (measured as potential denitrification and dissimilatory nitrate reduction to ammonium processing), nutrient stocks (total carbon, total nitrogen, and organic matter), and stable isotopes (carbon and nitrogen). Water level data related to each wetland location can be found at https://data.ess-dive.lbl.gov/view/doi:10.15485/2530253 (Kirker et al., 2024) and related water geochemistry data can be found at https://data.ess-dive.lbl.gov/datasets/doi:10.15485/3001967 (Forbes et al., 2025). In addition to a readme, this data package also includes a file-level metadata (FLMD) file that describes each file and a data dictionary (DD) that describes all column/row headers and variable definitions. This dataset is comprised of (1) file-level metadata; (2) data dictionary; (3) field metadata and international generic sample numbers (IGSNs); (4) readme; (5) the field protocol; and (6) a subfolder with sample data. The sample data subfolder contains (1) sediment potential denitrification rate, (2) sediment potential dissimilatory nitrate reduction to ammonium (DNRA) rate, (3) sediment total carbon and nitrogen content, (4) sediment stable isotopes (delta nitrogen-15 and delta carbon-13), (5) sediment percent organic matter, (6) surface water nitrous oxides, (7) surface water ammonium, and (8) methods codes. All files are .csv or .pdf.

Ammonium↗
DOE OSTI2024

Evaluation of Rooftop Solar Potential in Chernihiv and Lviv, Ukraine, and Efficacy of High-Resolution 3D Data Digital Twins [Slides]

NREL evaluated rooftop solar photovoltaic (PV) siting opportunities in the cities of Chernihiv and Lviv, Ukraine, leveraging very-high-resolution 3D elevation data to calculate technical potential. The study assessed total rooftop solar capacity and annual energy production. In Chernihiv and Lviv, 116,503 buildings were analyzed for their rooftop solar potential. The buildings in the study areas include a mixture of residential, commercial, and industrial buildings, characterized by diverse roof shapes and sizes. The total estimated rooftop solar capacity is 332 MWDC in Chernihiv and 873 MWDC in Lviv with annual energy production up to 376.2 GWhDC in Chernihiv and 995.5 GWhDC in Lviv. This accounting provides a clear estimate of the potential rooftop solar installations that could be realized under optimal conditions, considering both technical constraints and the geographic distribution of available rooftop space. Related to Russia's invasion of Ukraine, NREL estimates loss from buildings damaged or destroyed of 2,754 buildings, 20.15 MWDC of capacity, and 22,869 MWhDC of annual energy production in Chernihiv as well as a loss of 1,316 buildings, 34.49 MWDC of capacity, and 39,369 MWhDC of annual energy production in Lviv. The study also assessed the feasibility of adapting this methodology on a national scale using either simulated digital surface models (DSMs) or a digital twin approach. While the very-high-resolution DSM provided more precise results, the simulated DSM demonstrated reasonable accuracy for broader applications in modeling aggregated distributed solar supply. The feasibility of using a digital twin for Ukraine's national rooftop solar potential is considered promising, with certain limitations in areas with highly variable building stock and heavy war damage.

14 SOLAR ENERGY↗
NASA NTRS2003

Ventricular filling slows epicardial conduction and increases action potential duration in an optical mapping study of the isolated rabbit heart

INTRODUCTION: Mechanical stimulation can induce electrophysiologic changes in cardiac myocytes, but how mechanoelectric feedback in the intact heart affects action potential propagation remains unclear. METHODS AND RESULTS: Changes in action potential propagation and repolarization with increased left ventricular end-diastolic pressure from 0 to 30 mmHg were investigated using optical mapping in isolated perfused rabbit hearts. With respect to 0 mmHg, epicardial strain at 30 mmHg in the anterior left ventricle averaged 0.040 +/- 0.004 in the muscle fiber direction and 0.032 +/- 0.006 in the cross-fiber direction. An increase in ventricular loading increased average epicardial activation time by 25%+/- 3% (P < 0.0001) and correspondingly decreased average apparent surface conduction velocity by 16%+/- 7% (P = 0.007). Ventricular loading did not significantly alter action potential duration at 20% repolarization (APD20) but did at 80% repolarization (APD80), from 179 +/- 7 msec to 207 +/- 5 msec (P < 0.0001). The dispersion of APD20 was decreased with loading from 19 +/- 2 msec to 13 +/- 2 msec (P = 0.024), whereas the dispersion of APD80 was not significantly changed. These electrophysiologic changes with ventricular loading were not affected by the nonspecific stretch-activated channel blocker streptomycin (200 microM) and were not attributable to changes in myocardial perfusion or the presence of an electromechanical decoupling agent (butanedione monoxime) during optical mapping. CONCLUSION: Acute loading of the left ventricle of the isolated rabbit heart decreased apparent epicardial conduction velocity and increased action potential duration by a load-dependent mechanism that may not involve stretch-activated channels.

NASA Program Biomedical Research and Countermeasur↗
NASA NTRS2001

Dipole potentials indicate restructuring of the membrane interface induced by gadolinium and beryllium ions

The dipole component of the membrane boundary potential, phi(d), is an integral parameter that may report on the conformational state of the lipid headgroups and their hydration. In this work, we describe an experimental approach to measurements of the dipole potential changes, Deltaphi(d), and apply it in studies of Be(2+) and Gd(3+) interactions with membranes composed of phosphatidylserine (PS), phosphatidylcholine (PC), and their mixtures. Deltaphi(d) is determined as the difference between the changes of the total boundary potential, phi(b), measured by the IFC method in planar lipid membranes and the surface potential, phi(s), determined from the electrophoretic mobility of liposomes. The Gouy-Chapman-Stern formalism, combined with the condition of mass balance, well describes the ion equilibria for these high-affinity cations. For the adsorption of Be(2+) and Gd(3+) to PC membranes, and of Mg(2+) to PS membranes, the values of Deltaphi(b) and Deltaphi(s) are the same, indicative of no change of phi(d). Binding of Gd(3+) to PS-containing membranes induces changes of phi(d) of opposite signs depending on the density of ionized PS headgroups in the bilayer. At low density, the induced Deltaphi(d) is negative (-30 mV), consistent with the effect of dehydration of the surface. At maximal density (pure PS, neutral pH), adsorption of Be(2+) or Gd(3+) induces an increase of phi(d) of 35 or 140 mV, respectively. The onset of the strong positive dipole effect on PS membranes with Gd(3+) is observed near the zero charge point and correlates with a six-fold increase of membrane tension. The observed phenomena may reflect concerted reorientation of dipole moments of PS headgroups as a result of ion adsorption and lipid condensation. Their possible implications to in-vivo effects of these high-affinity ions are discussed.

Non-NASA Center↗
DOE OSTI2025

Floating Photovoltaic Technical Potential: A Novel Geospatial Approach on Federally Controlled Reservoirs in the United States

Floating photovoltaic generation is a rapidly expanding sector of the solar energy industry, and understanding the quantity that can feasibly be installed is a crucial step to understand its role in future energy systems. This paper presents a novel spatially explicit methology of FPV potential for federally owned and managed reservoirs in the United States that uses site-specific attributes of reservoirs to estimate available area and potential generation capacity. The analysis finds that the proportion reservoir area that is found to be available for FPV development is similar to assumed values used in previous research on average, however there is a wide variability in this proportion on a site by site basis. Potential FPV generation capacity on these reservoirs is estimated to be in the range of 861 to 1,042 GWdc depending on input assumptions, likely representing a significant portion of future US solar generation needs. This work represents an advancement in methods used to estimate FPV potential that presents many natural extensions for further research.

floating solar↗
DOE OSTI2024

Multi-level Monte Carlo methods in chemical applications with Lennard-Jones potentials and other landscapes with isolated singularities

We describe and compare outcomes of various Multi-Level Monte Carlo (MLMC) method variants, motivated by the potential of improved computational efficiency over rejection based Monte Carlo, which scales poorly with problem dimension. With an eye toward its application to computational chemical physics, we test MLMC's ability to sample trajectories on two problems — a familiar double-well potential, with known stationary distributions, and a Lennard-Jones solid potential (a Galton Board). By sampling Brownian motion trajectories, we are able to compute expectations of observable averages. These multi-basin potential energy problems capture the essence of the challenges with using MLMC, namely, maintaining correspondence of sample paths as time-resolution is varied. Addressing this challenge properly can lead to MLMC significantly outperforming standard Monte Carlo path sampling. We describe the essence of this problem and suggest strategies that circumvent diverging multilevel sample paths for an important class of problems. In the tests we also compare the computational cost of several, “adaptive,” variants of MLMC. Our results demonstrate that MLMC overcomes the collision, time scale limitation of the more familiar Brownian path MC samplers, and our implementation provides tunable error thresholds, making MLMC a promising candidate for application to larger and more complex molecular systems.

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

Hydropower potential derived from streamflow extremes for Alaska, USA

Alaska is an expansive region known for its abundant natural resources, including thousands of miles of streams and rivers. These rivers represent potential opportunities for future hydropower development that could provide reliable energy supply for local communities. There is limited long-term high temporal resolution streamflow data available for the region, making data-driven estimates of potential hydropower and its variability across the state challenging. This study provides a novel data-driven approach for hydropower capacity estimation across Alaska. We use supervised machine learning to develop a relationship between the daily and peak flow duration curves in order to augment the size of our dataset from 44 sites to 67 sites. We perform a stochastic hydropower estimation across the 67 sites and identify approximately 1000 MW of total potential hydropower capacity distributed across these sites. Our study provides the first step towards more comprehensive hydropower estimation for this critical region, highlighting the need for future work integrating high-resolution spatial data, community needs, and economic constraints in estimates of potential hydropower development in Alaska.

Hydropower↗
DOE OSTI2024

Machine learning for the redox potential prediction of molecules in organic redox flow battery

Here, organic redox flow batteries (ORFB) are recognized as an innovative technology for the large-scale storage of renewable energy. The redox potential of organic redox-active molecules plays a vital role in their performance. Advanced screening techniques like high-throughput experiment and machine learning (ML) have significantly enhanced organic material performance and transformed the field of ORFB. However, the scarcity of experimental data poses a considerable challenge for ML model development in this domain. In our study, we developed lightweight graph-based Gaussian process regression (GPR) models with GPU-accelerated marginalized graph kernel and hybrid kernel to predict the redox potentials of organic redox-active molecules for ORFBs, specifically focusing on small datasets. To evaluate model accuracy, we created a new experimental database of organic redox-active molecules by the data from hundreds of published papers and assembled previous computational datasets. We also considered some key parameters, such as pH conditions and solvent type, to assess their impact on redox potential prediction. Our GPR model predicted redox potentials with high accuracy across all datasets using minimal training data. The study provides powerful tools for molecule screening and design and delivers valuable guidance on designing training datasets for costly experiments.

25 ENERGY STORAGE↗
DOE OSTI2026

A Transferable Force Field for Simulating Adsorption in Metal–Organic Frameworks with Open Metal Sites Based on the 12–6–4 Lennard-Jones Potential

Metal−organic frameworks (MOFs) that contain coordinatively unsaturated open metal sites (OMSs) provide strong host− guest interactions, making them promising sorbents for low-concentration gas adsorption applications such as direct air capture and atmospheric water harvesting. However, accurately modeling host−guest interactions involving OMSs remains challenging for classical force fields (FFs) based on the 12−6 Lennard−Jones (LJ) potential, as the polarization effect of the guest molecule induced by the positively charged OMS is not considered. Here, we introduce an FF based on the 12−6−4 LJ potential, which incorporates charge−induced dipole interactions and is parametrized against a diverse set of host−guest potential energy surfaces (PESs) obtained from density functional theory (DFT). The resulting FF, trained on a generic trimetallic cluster, performs well in both host−guest binding energetics and gas adsorption isotherms across different OMS-containing MOFs, including MOF-74 series and Cu-BTC. These results highlight the excellent transferability of our approach and its potential to enhance the accuracy and robustness of high-throughput MOF discovery workflows, particularly for gas adsorption and separation in large and diverse MOF databases.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Hierarchical Reinforcement Learning of a Short-Range Bond-Order Potential for Silica: Analytic Embedding of Coordination with Classical Efficiency

Reinforcement learning (RL) has recently emerged as a data-efficient strategy to parametrize short-range interatomic potentials. Building on our past RL optimization of pairwise silica models, we extend the framework to a bond-order (Tersoff-type) potential that provides an analytic embedding of local coordination through a three-body term. A hierarchical RL workflow combining continuous-action Monte Carlo Tree Search and property-based rewards efficiently explores the 26-dimensional parameter space, sequentially optimizing lattice parameters, densities, angles, and cohesive energies of 21 silica polymorphs. The resulting models, Q-Tersoff and ML-Tersoff, reproduce the energetic ordering of low-energy phases and capture the angular correlations and amorphous structure factors of silica with improved fidelity over pairwise force fields, while remaining orders of magnitude faster than high-dimensional machine-learned potentials. Both models underperform for elastic constants and high-energy frameworks, delineating the limits of the current analytic form. The approach establishes a general and interpretable route to angle-aware, short-range potentials that bridge physics-based and machine-learned descriptions of silicate materials.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Liquid–Vapor Phase Equilibrium in Molten Aluminum Chloride (AlCl 3 ) Enabled by Machine Learning Interatomic Potentials

Molten salts are promising candidates in numerous clean energy applications, where knowledge of thermophysical properties and vapor pressure across their operating temperature ranges is critical for safe operations. Due to challenges in evaluating these properties using experimental methods, fast and scalable molecular simulations are essential to complement the experimental data. In this study, we developed machine learning interatomic potentials (MLIP) to study the AlCl 3 molten salt across varied thermodynamic conditions (T = 473–613 K and P = 2.7–23.4 bar), which allowed us to predict temperature-surface tension correlations and liquid–vapor phase diagram from direct simulations of two-phase coexistence in this molten salt. Two MLIP architectures, a Kernel-based potential and neural network interatomic potential (NNIP), were considered to benchmark their performance for AlCl 3 molten salt using experimental structure and density values. The NNIP potential employed in two-phase equilibrium simulations yields the critical temperature and critical density of AlCl 3 that are within 10 K (∼3%) and 0.03 g/cm 3 (∼7%) of the reported experimental values. An accurate correlation between temperature and viscosities is obtained as well. In doing so, we report that the inclusion of low-density configurations in their training is critical to more accurately represent the AlCl 3 system across a wide phase-space. The MLIP trained using PBE-D3 functional in the ab initio molecular dynamics (AIMD) simulations (120 atoms) also showed close agreement with experimentally determined molten salt structure comprising Al 2 Cl 6 dimers, as validated using Raman spectra and neutron structure factor. Furthermore, the PBE-D3 as well as its trained MLIP showed better liquid density and temperature correlation for AlCl 3 system when compared to several other density functionals explored in this work. Overall, the demonstrated approach to predict temperature correlations for liquid and vapor densities in this study can be employed to screen nuclear reactors-relevant compositions, helping to mitigate safety concerns.

Ab initio molecular dynamics↗
DOE OSTI2025

Crystal Nucleation Kinetics and Mechanism: Influence of Interaction Potential

Modulating liquid-to-solid transitions and the resulting crystalline structure for tailored properties is much desired. Colloidal systems are exemplary to this end, but the fundamental knowledge gaps in relating the influence of intermolecular interactions to crystallization behavior continue to hinder progress. In this study, we address this knowledge gap by studying nucleation and growth in systems with modified Lennard-Jones potential. Specifically, we study the commonly used 12-6 potential and a softer 7-6 potential. The thermodynamic state point for the study is chosen such that both systems are investigated at the same level of supercooling and pressure. Under these conditions, we find that the nucleation rate for both systems is comparable. Interestingly, the nucleation pathways and resulting crystal structures are different. In the 12-6 system, nucleation and growth occur predominantly through the FCC structure. Softening the potential alters the critical nucleus composition and introduces two distinct nucleation pathways. One pathway predominantly leads to the nucleus with a body-centered cubic (BCC) structure, while the other favors the face-centered cubic (FCC) arrangement. Our study illustrates that polymorph selection can be achieved through modifications to intermolecular interactions without impacting nucleation kinetics. The results have significant implications in designing approaches for polymorph selection and modulating self-assembly mechanisms.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Tracking Spatiotemporal Electric Potential in Batteries Using High-Resolution Operando X‑ray Transmission Imaging

The formation of significant concentration gradients across electrolytes in batteries affects the rate at which electrochemical reactions occur. In this work, we use high-resolution operando X-ray transmission imaging to capture spatiotemporal salt concentration profiles c(x,t) in a symmetric cell comprising a polymer electrolyte sandwiched between two lithium–indium alloy electrodes during a constant-current experiment followed by open-circuit relaxation. The decay of open-circuit potential is related to the concentration dependence of the potential across concentration cells, U. We show how operando c(x,t) data can be used to calculate the spatiotemporal electric potential “inside” the polarized electrolyte. We track the spatial- and time-dependent cell potential during the constant-current step and distinguish its two contributions: a concentration overpotential governed by U. and an ohmic contribution governed by ionic conductivity. Over most of the time window, the concentration overpotential dominates. At steady state, it is a factor of 7 larger than the ohmic contribution. Such findings indicate that efforts to design new polymer electrolytes should focus on minimizing concentration gradients.

Electrical conductivity↗
DOE OSTI2026

Machine Learning an Ab-Initio Based Bond-Order Potential for Bismuthene

Bismuthene is a heavy 2D material whose strong spin–orbit coupling and recently observed single-element ferroelectricity have intensified interest in its structural, vibrational, and transport properties. Accurate modeling of these behaviors requires a short-range interatomic potential that can reproduce the underlying bonding physics at a fraction of the computational cost of first-principles methods. However, such a potential is currently unavailable. Here, in this work, we construct a Tersoff bond-order potential for β-bismuthene using a reinforcement-learning framework that integrates a continuous Monte Carlo Tree Search with a simplex-based local optimizer. The optimized parameter sets reproduce first-principles lattice constants, cohesive energy, the equation of state, elastic constants, and phonon dispersion. We validate the models by performing thermal-conductivity calculations and uniaxial fracture simulations our findings confirm the reliability of the resulting models across multiple thermomechanical regimes. Comparison of the three best solutions reveals how differences in pairwise interactions, angular terms, and bond-order behavior govern phonon features and mechanical responses. We demonstrate an interpretable and computationally efficient potential for bismuthene and demonstrate a general reinforcement-learning strategy for developing bond-order models in emerging 2D materials.

deformation↗