Search NASA⌕ Search

SEARCH · Search NASA

Results for “DFT simulations”

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 595 records · Page 33
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Gas-phase negative ion photoelectron spectroscopy and reactivity of phenylacetylide

Arylacetylides such as phenylacetylide (PhCC–, 1) are important nucleophiles used in synthetic chemistry yet rarely have they been studied as bare carbanions. In this work, the phenylacetylide anion was formed via electrospray ionization (ESI) or multistage mass spectrometry (MSn) experiments and subsequently examined in the gas phase by negative ion photoelectron spectroscopy (NIPES), ion-molecule reactions (IMR), alongside theoretical calculations to probe its fundamental structure and reactivity. Photoelectron spectra of PhCC– (1) revealed vertical (VDE) and adiabatic detachment energies (ADE), both of 3.220 eV. The latter value is also the electron affinity (EA) of the phenylethynyl radical (PhCC•) from which a bond dissociation energy (BDE) of phenylacetylene (PhCCH) was derived to be ca. 131.3 ± 1.7 kcal mol-1 using a gas-phase thermochemical cycle. Detachment of an electron from PhCC– at 266 nm likely results in resonant autodetachment, supported by Franck-Condon Factor (FCF) simulations. The phenylacetylide anion reacts with methyl iodide (CH3I) and allyl iodide (C3H5I) via SN2 nucleophilic displacement with measured rate coefficients of 4.29 and 3.66 × 10-10 cm3 molecule-1 s-1, respectively. The mechanisms associated with these displacement reactions are explained and understood in terms of the reaction kinetics, natural bond orbital (NBO) theory, and Density Functional Theory (DFT) calculations.

Ma, Howard Z.↗
DOE OSTI2023

Atomistic insight into the effects of electrostatic fields on hydrocarbon reaction kinetics

Reactive Molecular Dynamics (MD) and Density Functional Theory (DFT) computations are performed to provide insight into the effects of external electrostatic fields on hydrocarbon reaction kinetics. By comparing the results from MD and DFT, the suitability of the MD method in modeling electrodynamics is first assessed. Results show that the electric field-induced polarization predicted by the MD charge equilibration method is in good agreement with various DFT charge partitioning schemes. Then, the effects of oriented external electric fields on the transition pathways of non-redox reactions are investigated. Results on the minimum energy path suggest that electric fields can cause catalysis or inhibition of oxidation reactions, whereas pyrolysis reactions are not affected due to the weaker electronegativity of the hydrogen and carbon atoms. MD simulations of isolated reactions show that the reaction kinetics is also affected by applied external Lorentz forces and interatomic Coulomb forces since they can increase or decrease the energy of collision depending on the molecular conformation. In addition, electric fields can affect the kinetics of polar species and force them to align in the direction of field lines. These effects are attributed to energy transfer via intermolecular collisions and stabilization under the external Lorentz force. The kinetics of apolar species is not significantly affected mainly due to the weak induced dipole moment even under strong electric fields. The dynamics and reaction rates of species are studied by means of large-scale combustion simulations of n-dodecane and oxygen mixtures. Results show that under strong electric fields, the fuel, oxidizer, and most product molecules experience translational and rotational acceleration mainly due to close charge transfer along with a reduction in their vibrational energy due to stabilization. This study will serve as a basis to improve the current methods used in MD and to develop novel methodologies for the modeling of macroscale reacting flows under external electrostatic fields.

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

DFT-based insight into finite-temperature properties of ferroelectric perovskites with lone-pair: the case of CsGeX 3 (X = Cl, Br, I)

Ferroelectrics remain in the focus of scientific attention for decades owing to their fundamental and practical appeal. Recently, ferroelectricity has been demonstrated in semiconducting halide perovskites (Zhang et al 2022 Sci. Adv. 8 eabj5881), offering both a rare combination of ferroelectricity and semiconductivity in the same material and a possible alternative to the prevailing perovskite oxide ferroelectrics. We propose a route to simulating such materials at finite temperatures capable of reproducing key experimental and first-principle data, such as Curie temperature, phase transition sequence, spontaneous polarization, and soft mode frequencies. The key methodological finding is the superior performance of hybrid exchange correlation functionals in parametrization of effective Hamiltonians for ferroelectrics with lone pair. The parametrization for effective Hamiltonians for CsGeX 3 (X = Cl, Br, I) is reported. The application of methodology to study polarization reversal in CsGeX 3 allows for the development of a ‘minimalistic’ model for polarization reversal in ferroelectrics that provides an insight into the mechanisms of polarization reversal and its key features, such as the relationship between the coercive field, temperature, and AC field frequency. Importantly, the model reveals the origin of the well-known and ever-puzzling overestimation of coercive fields in computations. Furthermore, we report a variety of finite-temperature properties of CsGeX 3 ferroelectrics, such as dielectric susceptibility, pyroelectric coefficients, and energy storage density, which reveal that these halide perovskites possess properties comparable to their oxide counterparts. Here, we believe that our work provides significant methodological advancements, deepens fundamental understanding of ferroelectrics, and reveals the potential of halide perovskite ferroelectrics.

effective Hamiltonian↗
DOE OSTI2023

Using Metadynamics to Reveal Extractant Conformational Free Energy Landscapes

Understanding the impact of extractant functionalization on metal-binding energetics in liquid-liquid extraction is essential to guide the development of better separation processes. Traditionally, computational extractant design uses electronic structure calculations on metal-ligand clusters to determine the metal-binding energy of the lowest energy state. Although highly accurate, this approach does not account for all of the relevant physics encountered under experimental conditions. Such methodologies often neglect entropic contributions such as temperature effects and ligand flexibility, in addition to approximating solvent-extractant interactions with implicit solvent models. In this study, we use classical molecular dynamics simulations with an advanced sampling method, metadynamics, to map out extractant molecule conformational free energies in the condensed phase. Here we generate the complete conformational landscape in solution for a family of bidentate malonamide-based extractants with different functionalizations of the headgroup and the side chains. In particular, we show how such alkyl functionalization reshapes the free energy landscape, affecting the free energy penalty of organizing the extractant into the cis-like metal-binding conformation from the trans-like conformation of the free extractant in solution. Specifically, functionalizing alkyl tails to the center of the headgroup has a greater influence on increasing molecular rigidity and disfavoring the binding conformation than functionalizing side chains. These findings are consistent with trends in metal-binding energetics based on experimentally reported distribution ratios. We also consider a different bidentate extractant molecule, carbamoylmethylphosphine oxide, and show how the choice of solvent can further reshape the conformational energetic landscape. This study demonstrates the feasibility of using molecular dynamics simulations with advanced sampling techniques to investigate extractant conformational energetics in solution, which, more broadly, will enable extractant design that accounts for entropic effects and explicit solvation.

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

Irradiation enhanced diffusion and diffusional creep in U₃Si₂

U₃Si₂ is a candidate as an advanced fuel for light water reactors (LWRs). Compared to UO₂, there are significant data gaps for the thermophysical and thermomechanical properties of U₃Si₂. Atomic scale simulations can help fill in those gaps, while also providing a more mechanistic understanding of the material properties. Creep is an important property for fuel performance as it influences the pellet-cladding gap, which in turn affects the temperature of the fuel (a key parameter in many other aspects of fuel behavior). In this work, we employed a DFT-informed cluster dynamics framework to predict point defect concentrations in U₃Si₂. The concentrations, diffusivities, and volumes of these defects were used to develop a diffusional creep model based on bulk (Nabarro-Herring) and grain boundary (Coble) diffusion processes. It was found that irradiation enhanced Nabarro-Herring creep is dominant at low temperatures, while Coble creep dominates for high temperatures. The model compares well against available experimental data and has been implemented in BISON. A demonstration case using simple power profiles has been carried out, showing that negligible creep occurs due to the low temperatures experienced by U₃Si₂, a consequence of its high thermal conductivity.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Mechanism of Solid Electrolyte Interphase Formation on Lithium Anodes Studied by Molecular Dynamics Simulations

We investigated the initial formation of the SEI on lithium-metal anodes in presence of several electrolyte compositions using ReaxFF, in absence of a bias potential. We observed the formation of well-defined regions with different density and oxidation states of Li. Dissolution of Li metal atoms results in an uneven distribution of Li in the liquid phase. The electrolyte decomposition mechanism was studied for various electrolyte mixtures. Overall, the reaction mechanisms observed with ReaxFF are in good agreement with previous reports using DFT. However, we found some discrepancies in specific predictions by the ReaxFF that should be addressed in future work.

25 ENERGY STORAGE↗
DOE OSTI2025

Fundamental microscopic properties as predictors of large-scale quantities of interest: Validation through grain boundary energy trends

Correlations between fundamental microscopic properties computable from first principles, which we term canonical properties, and complex large-scale quantities of interest (QoIs) provide an avenue to predictive materials discovery. Here, we propose that such correlations can be efficiently discovered through simulations utilizing approximate interatomic potentials (IPs), which serve as an ensemble of “synthetic materials”. As a proof of principle we build a regression model relating canonical properties to the symmetric tilt grain boundary (GB) energy curves in face-centered cubic crystals, characterized by the scaling factor in the universal lattice matching model of Runnels et al. (2016), which we take to be our QoI. Our analysis recovers known correlations of GB energy to other properties and discovers new ones. We also demonstrate, using available density functional theory (DFT) GB energy data, that the regression model constructed from IP data is consistent with DFT results, confirming the assumption that the IPs and DFT belong to same statistical pool and thereby validating the approach. Regression models constructed in this fashion can be used to predict large-scale QoIs based on first-principles data and provide a general method for training IPs for QoIs beyond the scope of first-principles calculations.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Perspective on Many-Body Methods for Molecular Polaritonic Systems

Recent advances in strong light–matter interactions have revealed a wealth of new physical phenomena in molecules embedded in optical cavities, including modified chemical reactivity, altered excitation spectra, and novel quantum correlations. To describe these effects from first-principles, the field of ab initio quantum electrodynamics (QED) has emerged as a compelling extension of quantum chemistry that treats electronic and photonic degrees of freedom on equal footing. In this Perspective, we review the growing landscape of many-body QED methods, including Hartree–Fock, density functional theory (QEDFT), time-dependent DFT (QED-TDDFT), configuration interaction (QED-CI), complete active space (QED-CASSCF), coupled cluster (QED-CC), quantum Monte Carlo (QED-QMC), and density matrix renormalization group (QED-DMRG), highlighting recent developments and implementations. We further explore real-time methods, gradient and Hessian formalisms, and the integration of nonadiabatic nuclear dynamics. Applications range from benchmark simulations of polaritonic chemistry to quantum simulations on emerging quantum hardware. We conclude by outlining future directions for theory development and interdisciplinary efforts at the interface of quantum chemistry, condensed matter, and quantum optics.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Parameterizing empirical interatomic potentials for predicting thermophysical properties via an irreducible derivative approach: the case of ThO 2 and UO 2

The accuracy of classical physical property predictions using molecular dynamics simulations is determined by the quality of the interatomic potentials. Here we introduce a training approach for empirical interatomic potentials (EIPs) which is well suited for capturing phonons and phonon-related properties. Our approach is based on direct comparisons of the second- and third-order irreducible derivatives (IDs) between an EIP and the Born–Oppenheimer potential within density functional theory (DFT) calculations. IDs fully exploit space group symmetry and allow for training without redundant information. We demonstrate the fidelity of our approach in the context of ThO 2 and UO 2 , where we optimize parameters of an embedded-atom method potential in addition to core–shell interactions. Our EIPs provide thermophysical properties in good agreement with DFT and outperform widely utilized EIPs for phonon dispersion and thermal conductivity predictions. Reasonable estimates of thermal expansion and formation energies of Frenkel pairs are also obtained.

empirical interatomic potential↗
NASA NTRS2002

Implementation of a 622 Mbps Digital Modem

This paper presents the implementation and initial test results of an Orthogonal Frequency Division Multiplexing (OFDM) digital modem (modulator and demodulator) with an aggregate information throughput of 622 megabits per second (Mbps). The OFDM waveform is constructed by dividing an incoming data stream into four channels, each channel using either a 16-ary Quadrature Amplitude Modulation (16QAM) scheme or an 8-Phase Shift Keying (8PSK) scheme. The generation and detection of the composite waveform are performed using Discrete Fourier Transform (DFT) and polyphase filtering, to digitally stack and band-limit the individual carriers respectively. The four-channel OFDM approach enables the implementation of a modem that can be both power and bandwidth efficient, with sufficient parallelism to meet higher data rate goals. As a result, the OFDM modem requires only a 240 MHz bandwidth to transmit 622 Mbps. Hardware and simulation results in the form of spectrum diagrams and bit-error-rate (BER) curves are also presented in this paper.

Kifle, Muli↗
DOE OSTI2024

A Molecular View of Methane Activation on Ni(111) through Enhanced Sampling and Machine Learning

A combination of machine learned interatomic potentials (MLIPs) and enhanced sampling simulations is used to investigate the activation of methane on a Ni(111) surface. The work entails the development and iterative refinement of MLIPs, initially trained on a dataset constructed via ab initio molecular dynamics (AIMD) simulations, supplemented by adaptive biasing forces, to enrich the sampling of catalytically relevant configurations. Our results reveal that by incorporating collective variables that capture the behavior of the reactant molecule, as well as additional frames that describe the dynamic response of the catalytic surface, it is possible to enhance considerably the accuracy of predicted energies and forces. By employing enhanced sampling schemes in the refinement of the MLIP, we systematically explore the potential energy surface, leading to a refined MLIP capable of predicting DFT-level energies and forces and replicating key geometric characteristics of the catalytic system. The resulting free energy landscapes at several temperatures provide a detailed view of the thermodynamics and dynamics of methane activation. Specifically, as methane approaches and dissociates on the catalytic surface, the process involves the dynamic interplay of CH 4 and the Ni catalyst that includes both enthalpic and entropic contributions. The progression towards the transition state involves an CH 4 moiety that is increasingly restrained in its ability to rotate or translate, while the stage following the transition state is characterized by a notable rise of the Ni atom that interacts with the cleaved C–H bond. Furthermore, this leads to an increase in the mobility of the adsorbed species, a feature that becomes more pronounced at higher temperatures.

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

Thermodynamic Stability and Site‐Specific Distribution of Graphitic and Pyridinic Nitrogen in Graphene Moiré on Ru(0001)

Abstract Graphene‐like materials are of interest for large‐scale hydrogen storage applications due to their lightweight, durable, and scalable properties. Nitrogen‐doping minimizes kinetic limitations in diffusion and recombination on surfaces, however, the role of graphitic nitrogen (GN) and pyridinic nitrogen (PN) is not well understood. Nitrogen‐doped graphene is synthesized on Ru(0001) using chemical vapor deposition (CVD) of pyridine and ion irradiation. Scanning tunneling microscopy (STM), x‐ray photoelectron spectroscopy (XPS), and density functional theory (DFT) are used to identify the structure, location, and thermodynamic stability of nitrogen species within the graphene moiré. CVD of pyridine results in a low nitrogen concentration (<0.1at%), while the post‐growth nitrogen ion irradiation allows us to increase the concentration further. The concentration of GN and PN is controlled by varying the ion dose and annealing temperature. Comparison of measured and simulated STM images of GN and PN yield an excellent agreement, allowing us to confidently establish that GN is preferentially located near the center of the Atop region, while PN is located in the valley region of the graphene moiré. This report explicitly confirms the site assignments and provides a foundation for the site synthesis and analysis of structural and electronic properties that drive the reactivity of N‐doped graphene.

Gedara, Buddhika S. A. [Physical and Computational↗
DOE OSTI2023

LAVA 1.0: A general-purpose python toolkit for calculation of material properties with LAMMPS and VASP

Here, we introduce LAVA, a general-purpose python toolkit to provide user-friendly and high throughput calculations for material properties using both LAMMPS and VASP. It contains a set of classes as well as pre-processing and post-processing functions to prepare, execute and extract information from LAMMPS/VASP simulations. An overview of the program structure is provided. The current version contains modules to calculate elastic and mechanical properties such as lattice constant, cohesive energy, cold curve, elastic constants, bulk and shear modulus, volume conserving/non-conserving deformation path, vacancy/interstitial formation energy, surface energy, stacking fault energy, melting point, radial distribution function, and thermal expansion. These functionalities are demonstrated for different interatomic potentials and DFT calculations in Al.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Defect modeling in semiconductors: the role of first principles simulations and machine learning

Abstract Point defects in semiconductors dictate their electronic and optical properties. Vacancies, interstitials, substitutional defects, and defect complexes can form in the semiconductor lattice and significantly impact its performance in applications such as solar absorption, light emission, electronics, and catalysis. Understanding the nature and energetics of point defects is essential for the design and optimization of next-generation semiconductor technologies. Here, we provide a comprehensive overview of the current state of research on point defects in semiconductors, focusing on the application of density functional theory (DFT) and machine learning (ML) in accelerating the prediction and understanding of defect properties. DFT has been instrumental in accurately calculating defect formation energies, charge transition levels, and other defect-related properties such as carrier recombination rates and lifetimes, and ion migration barriers. ML techniques, particularly neural networks, have emerged as powerful tools for enabling rapid prediction of defect properties at DFT-accuracy in order to overcome the expense of using large supercells and advanced functionals. We begin this article with a discussion of different types of point defects and complexes, their impact on semiconductor properties, and the experimental and DFT approaches typically used for their characterization. Through multiple case studies, we explore how DFT has been successfully applied to understand defect behavior across a variety of semiconductors, and how ML approaches integrated with DFT can efficiently predict defect properties and facilitate the discovery of new materials with tailored defect behavior. Overall, the advent of ‘DFT+ML’ promises to drive advancements in semiconductor technology, catalysis, and renewable energy applications, paving the way for the development of high-performance semiconductors which are defect-tolerant or have desirable dopability.

Rahman, Md Habibur (ORCID:000000027705984X)↗
DOE OSTI2024

Methyl α‐D‐galactopyranosyl‐(1→3)‐β‐D‐galactopyranoside and methyl β‐D‐galactopyranosyl‐(1→3)‐β‐D‐galactopyranoside: Glycosidic linkage conformation determined from MA'AT analysis

Abstract MA'ATanalysis has been applied to two biologically‐importantO‐glycosidic linkages in two disaccharides, α‐D‐Galp‐(1→3)‐β‐D‐GalpOMe (3) and β‐D‐Galp‐(1→3)‐β‐D‐GalpOMe (4). Using density functional theory (DFT) to obtain parameterized equations relating a group of trans‐O‐glycosidic NMR spin‐couplings to eitherphi(ϕ') orpsi(ψ'), and experimental 3 J COCH , 2 J COC , and 3 J COCC spin‐couplings measured in aqueous solution in 13 C‐labeled isotopomers, probability distributions ofϕ'andψ'in each linkage were determined and compared to those determined by aqueous 1‐μs molecular dynamics (MD) simulation. Good agreement was found between theMA'ATand single‐state MD conformational models of these linkages for the most part, with modest (approximately <15°) differences in the mean values ofϕ'andψ', although the envelope of allowed angles (encoded in circular standard deviations or CSDs) is consistently larger forϕ'determined fromMA'ATanalysis than from MD for both linkages. TheMA'ATmodel of the α‐Galp‐(1→3)‐β‐Galplinkage agrees well with those determined previously using conventional NMR methods ( 3 J COCH values and/or 1 H‐ 1 H NOEs), but some discrepancy was observed for the β‐Galp‐(1→3)‐β‐Galplinkage, which may arise from errors in the conventions used to describe the linkage torsion angles. Statistical analyses of X‐ray crystal structures show ranges ofϕ'andψ'for both linkages that include the mean angles determined fromMA'ATanalyses, although both angles adopt a wide range of values in the crystalline state, withϕ'in β‐Galp‐(1→3)‐β‐Galplinkages showing greater‐than‐expected conformational variability.

Chemistry↗
DOE OSTI2026

Effects of Composition and Oxidation States on the Structures of Chromium-Containing Sodium Silicate Glasses: Molecular Dynamics Simulations using Machine Learning Interatomic Potentials

Chromium represents a significant challenge for the vitrification of high-level nuclear waste into silicate and borosilicate glasses due to its low solubility and variable oxidation states, which can limit the waste loading due to promotion of crystallization or phase separation during processing. In this study, we modeled chromium containing silicate glasses using molecular dynamics simulations with three machine learning interatomic potentials (MLIPs), MACE, CHGNet, and PFP were employed, to gain insights on glass composition and oxidation states on the structures of these glasses. One of the goals is to evaluate their ability of these MLIPs to accurately represent the general structure of silicate glasses and chromium local environments as a function of chromium oxidation states. Density Functional Theory (DFT) based calculations and experimental data such as neutron structure factors were used to validate the structural models. It was found that the foundation models of all three MLIPs are able to reproduce general structural features of the sodium silicate glass structure consistent with experimental and DFT data, but only CHGNet and PFP can accurately capture the oxidation states and local environment of chromium: tetrahedral for Cr6+ and octahedral for Cr3+. Furthermore, we studied the effect of varying Cr3+/ Cr6+ (Cr3+/Crtotal) ratio and total chromium content using PFP. Our results show that Cr6+ enhances network polymerization by reducing non-bridging oxygens through Na? charge compensation required due to the formation of chromate (CrO42-) species, while Cr³? acts as a network modifier that disrupts connectivity. System size effects on the structural characteristics and chromium environments were also tested using the PFP potential. This work highlights the importance of careful validation on the precision, transferability, and potential of MLIPs for modeling glasses containing transition metal elements that can exist in multiple oxidation states. It is also encouraging to see the foundational models are all three MLFFs are able to reproduce the basic sodium silicate glass structures, while suggesting additional training or refining is needed to improve the description of more complex systems containing transition metals.

Puga, Christina L.↗
DOE OSTI2022

Density functional thermodynamic description of spin, phonon and displacement degrees of freedom in antiferromagnetic-to-paramagnetic phase transition in YNiO 3

This work herein demonstrates a direct density functional description of the finite-temperature thermodynamic properties of solids exhibiting phase transitions through positional and spin symmetry breaking degrees of freedom. A classic example addressed here is the rare-earth (R) nickelates RNiO 3 where the ground state is characterized by crystallographic and magnetic (e.g., antiferromagnetic) long-range order (LRO), whereas the higher temperature paramagnetic phase manifests a range of local spin and positional symmetry breaking motifs with short-range order (SRO). Unlike time-dependent simulations of spin and positional degrees of freedom, in the present work, phases are described via a superposition of static configurations constructed by populating a periodic base lattice supercell allowing for the formation of energy lowing distribution of positional and spin local motifs. The thermal populations of the configurations in such a superposition phase are obtained from the energy-minimized Density Functional Theory (DFT)-calculated partition functions at different temperatures. This approach offers flexible inclusion of different physical contributions to the free energy, such as elastic, electronic and phonon free energies, all obtained from the same underlying DFT total energy calculations of periodic structures. The thermodynamic and magnetic properties of both LRO and SRO crystallographic and spin phases, including antiferromagnetic (AFM) to paramagnetic (PM) Néel phase transition in YNiO3 are studied. Including spin and phonon contributions, we find a DFT-calculated Néel temperature to be 144 K in satisfactory agreement with the experimental value of 145 K; whereas omitting the phonon contribution, one obtains a Néel temperature of 81 K. We present phonon contributions to the DFT-calculated temperature-dependent SRO, heat capacities, and the polymorphous distribution of nonzero local magnetic moments in the PM phase. This approach thus extends to finite temperatures the symmetry-broken DFT description of both the AFM and PM phases, demonstrating that a thermodynamic superposition approach based on symmetry broken configurations evaluated by a mean-field like DFT is sufficient to obtain a consistent description of the thermal physics of the AFM, PM phases and their interconversion in 3d oxides illustrated by YNiO 3 .

36 MATERIALS SCIENCE↗
DOE OSTI2026

Maximizing machine learning interatomic potential transferability for the discovery of the novel stellated octadecagon Bi18-Pt24 cage structure

Achieving true transferability remains the central challenge for Machine Learning Interatomic Potentials (ML-IAPs) in modeling complex bimetallic nanoclusters across their vast potential energy surfaces. We systematically investigate data selection strategies to optimize the Chebyshev Interaction Model for Efficient Simulation (ChIMES) potential for the Bi-Pt nanoclusters by comparing three innovative sampling methods: Principal Component Analysis (PCA)/k-means (structural diversity), t-distributedStochasticNeighborEmbedding (t-SNE)/k-means (force-space diversity), and hierarchical clustering. Quantitatively, the PCA/k-means strategy proved most effective for global accuracy, yielding the lowest force errors and achieving energy root mean square errors (RMSE) values competitive with Density Functional Theory (DFT), demonstrating excellent accuracy (19.16meV/atom). Structural validation on 34 unique DFT-optimized isomers further confirmed the potential’s high fidelity, with the best model PCA/k-means reproducing structures with an average root mean square deviation (RMSD) of 0.10 Å. However, the t-SNE methods, by maximizing diversity in the force space, demonstrated superior extrapolative power, leading to the more precise prediction of a novel stellated octadecagon Bi18⁢Pt24 cage structure, demonstrating the potential for exploring previously unseen morphologies. Our results establish a clear methodology for strategic data sampling that successfully maximizes ML-IAP transferability, providing an accurate and computationally efficient tool that accelerates the theoretical discovery of complex bimetallic architectures.

Vangheluwe, Raphaël [Université Paris-Saclay, CNRS↗