Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Reliability”

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 1,297 records · Page 72
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Standardized nomenclature for photovoltaic connectors

Photovoltaic (PV) systems rely on discrete connectors for the efficient and safe flow of power from module to module and from strings to combiner boxes and inverters. Despite their functional importance, no common nomenclature for PV connectors currently exists, resulting in confusion and miscommunication. Misunderstood terms like "MC4 compatible", "cross-mating", "intermating", "female", and "male" can lead to installation and maintenance errors and compromise system reliability. We believe a standardized terminology will reduce confusion, help support installation best practices, aid in maintenance and repair, inform next-generation designs, and provide a technical basis for improved codes and standards. To that end, we are proposing a standardized glossary for 4 mm PV connectors (the most common type of connector used in PV applications) based on, and validated by, a Sandia National Laboratories' investigation that included the following sources: 1) a comprehensive review of official documents from 20 connector manufacturers, including schematics, datasheets, installation manuals, and catalogs, as well as relevant patents; 2) two rounds of surveys distributed to stakeholders, including connector manufacturers, engineers, asset owners, test labs, and researchers; and 3) visual examination of 25 different models of 4 mm single-pole DC PV connectors to document variations in design and functionality. This work provides a foundation for establishing a clear and consistent terminology for PV connectors that will in turn enable progress toward greater reliability and collaboration across the industry.

14 SOLAR ENERGY↗
DOE OSTI2022

Overcoming significant challenges in extracting off-stoichiometric thermodynamics using the compound energy formalism through complementary use of experimental and first principles data: A case study of Ba 1-x Sr x FeO 3-δ

The compound energy formalism (CEF) is a powerful framework to describe the thermodynamics of metal oxides as a function of off-stoichiometry, temperature, and composition. The thermodynamic properties are crucial materials design attributes in metal oxide-based oxygen-exchange chemical processes. Despite the richness of information an accurate CEF model provides, a method to determine a unique and accurate fit for oxygen exchange materials remains elusive. This contribution details a method for fitting the CEF model that overcomes the current fitting challenges through three innovations: 1) the combination of density functional theory calculations with experimental data decorrelates excess terms and delineates the enthalpic/entropic contributions to the Gibbs free energy; 2) a systematic determination of the important CEF model terms, removing thermodynamic predetermining human intervention; 3) a self-consistent solution of the starting oxygen offstoichiometry (δ 0 ) of thermogravimetric measurements. Thus, our method enables the reliable extraction of off-stoichiometric metal oxide thermodynamic properties and facilitates rapid materials compositional screening, and reliable process design of systems dependent on off-stoichiometric redox-active metal oxides. We apply this method to a Ba x Sr 1-x FeO 3-δ test case. We find by systematically examining the performance of the CEF model fit with and without each innovation that all three innovations are necessary for an accurate fit. We determined that reduction enthalpy is higher and more sensitive to off-stoichiometry when the Sr fraction is large (139.5 and 185.3 kJ/mol O 2 for SrFeO 3 at δ = 0 and δ = 0.5, respectively vs. nearly constant 83 kJ/mol O 2 for BaFeO 3 ). However, the reduction entropy is mostly insensitive to Sr fraction, but highly dependent on δ suggesting larger contributions in the non-configurational entropy parameters.

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

Mechanisms of Degradation of Toxic Nerve Agents: Quantum-chemical Insight into Interactions of Sarin and Soman with Molybdenum Dioxide

With the urgent need to provide an efficient and reliable protection for people from fatal chemical and biological weapons, our fundamental understanding of how toxins interact with filters is far from complete. The situation is further complicated by natural difficulties of performing experimental measurements with lethal toxins. Unlike experiments, computational modeling offers an attractive and safe yet reliable way of studying behavior of toxic agents on a variety of substrates at a great level of detail. Here, we report DFT-based quantum-chemical calculations of adsorption and decomposition of DMMP, sarin and soman on MoO 2 (011) surface. Our calculations show that MoO 2 strongly adsorbs toxic nerve agents and quickly decomposes them. Decomposition of DMMP on the MoO 2 (011) surface proceeds via the PO-CH 3 bond breaking and a formation of a surface methoxy group. The calculated activation barrier for this reaction is 131.5 kJ mol -1 . Unlike DMMP, decomposition of sarin and soman proceeds via the dealkylation reaction yielding propene and 3,3-dimethyl-1-butene, respectively. Decomposition of sarin requires a remarkably low energy (53.7 kJ mol- 1 ), whereas the similar reaction in soman requires 50 kJ mol -1 more energy. We also make specific predictions to guide Ambient-Pressure X-ray Photoelectron Spectroscopy (APXPS) experiments on sarin interaction with MoO 2 samples. We conclude that MoO 2 serves as an efficient substrate able of degrading toxins.

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

Utility Regulation Through Legislation: A Cautionary Tale for Legislators, Regulators, Stakeholders, and Utilities

State legislatures in the U.S. are becoming increasingly involved in utility regulation. While state utility legislation can expand regulator authority, for example in the case of renewable energy or energy efficiency standard administration, it can also reduce it. Given the challenging nature of for-profit monopoly regulation, state utility legislation often reduces regulator authority unintentionally. This editorial begins with an examination of the outcomes of state legislation related to baseload generation in the U.S. and the resulting impacts to customers, shareholders, and state economies. The record indicates that utility regulation through state legislation has caused significant economic harm, generally by encouraging for-profit utilities to make riskier investments. A description of more recent state legislation involving distribution grids is presented, along with evidence that increases in grid investment encouraged by state legislation has not improved the reliability of U.S. for-profit utilities. The authors claim that features common to state utility legislation, particularly the requirement that utilities provide investment plans to regulators in advance, practically eliminates cost disallowance risk. This, in turn, reduces regulator authority despite legislative provisions intended to protect consumers. The authors hypothesize that the practical elimination of cost disallowance risk, combined with information and expertise asymmetry among regulators and stakeholders, encourages for-profit monopoly utilities to make grid investments that are not cost effective. The editorial concludes that legislators should avoid utility regulation through legislation, thereby preserving regulator authority, whenever possible. The authors also present their fifth annual Customer Value Ranking of U.S. investor-owned electric utilities. The ranking represents a comparison of the benefits utilities deliver to customers (measured by reliability performance and customer satisfaction) to the costs customers pay (measured by O&M and capital spending per customer).

ENERGY PLANNING, POLICY, AND ECONOMY,POWER TRANSMI↗
DOE OSTI2022

Machine learning for ultrasonic nondestructive examination of welding defects: A systematic review

Recent years have seen a substantial increase in the application of machine learning (ML) for automated analysis of nondestructive examination (NDE) data. One of the applications of interest is the use of ML for the analysis of data from in-service inspection of welds in nuclear power and other industries. These types of inspections are performed in accordance with criteria described in the ASME Boiler and Pressure Vessel Code and require the use of reliable NDE techniques. The rapid growth in ML methods and the diversity of possible approaches indicate a need to assess the current capabilities of ML and automated data analysis for NDE and identify any gaps or shortcomings in current ML technologies as applied to the automated analysis of NDE data. In particular, there is a need to determine the impact of ML on the NDE reliability. This paper discusses the findings from a literature survey on the current state of ML for the automated analysis of data from ultrasonic NDE of weld flaws. It discusses an overview of ultrasonic NDE as used for weld inspections in nuclear power and other industries. Herein, data sets and ML models used in the literature are summarized, along with a generally applicable workflow for ML. Findings on the capabilities, limitations and potential gaps in feature selection, data selection, and ML model optimization are discussed. The paper identified several needs for quantifying and validating the performance of ML methods for ultrasonic NDE, including the need for common data sets.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Development and testing of coarse-grained models for ultrasonic simulations of cast austenitic stainless steel

Ultrasonic inspection of cast austenitic stainless steel (CASS) in the nuclear industry is particularly challenging because of sound field scatter and attenuation caused by the coarse-grained microstructure. Modeling and simulation are important tools in ultrasonic testing, as they can be used to help address key aspects of inspections, such as developing new probe designs, predicting inspection reliability, and testing phased-array focal laws. However, developing a useful and reliable CASS model is challenging due to the many grain interfaces and crystalline orientations that must be captured. We demonstrate a method of creating a realistic CASS model that is usable in CIVA, a commercially available modeling and simulation software platform. Using polished and chemically etched sections, we generate models of a coarse-grained equiaxed specimen and a columnar specimen. We also test an alternative method of generating a coarse-grained model using Voronoi regions. We qualitatively compare sound field scatter and quantitatively compare sound field attenuation and beam partitioning in simulated sound fields to those of laboratory-measured sound fields. Results show that the Voronoi models perform as well as or better than the models based on actual grain morphology. Here we also show that model-to-model randomness in Voronoi grain structure can impact the magnitude of a simulated echo response by a factor of two or more. Although CASS models are potentially a good depiction of reality for a given scenario, they should not be considered representative since CASS morphology can change significantly from specimen to specimen or within the same specimen.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Techno-economic assessment of distributed wellhead RO water treatment for nitrate removal and salinity reduction: A field study in small disadvantaged communities

Techno-economic analysis of distributed wellhead water treatment and desalination (DWTD) systems was carried out based on a three-year field study in three small, disadvantaged communities (DACs) to evaluate the reliability and affordability of upgrading their impaired well water. The local water supplies of the three study DACs, located in Salinas Valley, California, were contaminated with nitrate at levels (~ 12–87 mg/L NO$^{-}_{3}$ - N) ) above the California maximum contaminant level (MCL) of 10 mg/L NO$^{-}_{3}$ - N , and had elevated water salinity (~600–1,600 mg/L total dissolved solids(TDS)) above its secondary MCL (SMCL) of 500 mg/L TDS. Well water nitrate removal and salinity reduction were accomplished via reverse osmosis (RO) based DWTD systems that operated autonomously, supported by remote monitoring and supervisory cyberinfrastructure. Reliable DWTD operation provided treated water quality, with respect to nitrate and salinity, in the range of 0.5–6.3 mg/L NO$^{-}_{3}$ - N and 57–161 mg/L TDS, respectively, which were well below the respective MCL and SMCL. The levelized cost of water treatment was in the range of ~$$2/m 3 - $$2.9/m 3 which aligns with typical residential water costs in California and in the study region, and monthly residential water costs (39 dollars-74 dollars/residential unit/month) were also within the range in California. The study showcased the DWTD approach as a viable and potentially scalable solution for upgrading impaired local potable water supply of communities lacking centralized water delivery infrastructure. However, streamlined permitting processes and standardized regulatory frameworks are critical to promoting wider adoption and maximizing the socio-economic benefits of the DWT approach. Moreover, DACs are likely to require government subsidies in order to cover the CapEx of DWTD systems in addition to upgrade of site infrastructure.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2023

Using low-coverage whole genome sequencing (genome skimming) to delineate three introgressed species of buffalofish ( Ictiobus )

Consumption of buffalofish has been sporadically associated with Haff disease-like illnesses involving sudden onset muscle pain and weakness due to skeletal muscle rhabdomyolysis, but determination of precisely which species are associated with these illnesses has been impeded by a lack of species-specific DNA-based markers. Here, three closely related species of buffalofish native to the Mississippi River Basin (Ictiobus bubalus, Ictiobus cyprinellus and Ictiobus niger) that have previously proven genetically indistinguishable using both mitochondrial and nuclear single-locus sequencing were reliably discriminated using low-coverage whole genome sequencing (‘genome skimming’). Using 44 specimens representing the three species collected from the mid/upper (Missouri) and lower (Louisiana) regions of the species’ native ranges, the SISRS (Site Identification from Short Read Sequences) bioinformatics pipeline was adapted to (1) identify over 620Mbp of putatively homologous nuclear sequence data and (2) isolate over 140,000 single-nucleotide polymorphisms (SNPs) that supported accurate species delimitation, all without the use of a reference genome or annotation data. These sites were used to classify Ictiobus spp. samples with genome-skim data, along with a larger set (n = 67) where ultraconserved elements (UCEs) were sequenced. Analyses of whole mitochondrial data revealed more limited signal. Nearly all samples matched their purported species based on morphologic identification, but two Missouri samples morphologically identified as I. niger grouped with samples of I. bubalus, albeit with significant enrichment of I. niger SNPs. To our knowledge this is the first report of a DNA-based tool to reliably discriminate these three morphologically distinct species.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Neural network ensembles and uncertainty estimation for predictions of inelastic mechanical deformation using a finite element method-neural network approach

The finite element method (FEM) is widely used to simulate a variety of physics phenomena. Approaches that integrate FEM with neural networks (NNs) are typically leveraged as an alternative to conducting expensive FEM simulations in order to reduce the computational cost without significantly sacrificing accuracy. However, these methods can produce biased predictions that deviate from those obtained with FEM, since these hybrid FEM-NN approaches rely on approximations trained using physically relevant quantities. In this work, an uncertainty estimation framework is introduced that leverages ensembles of Bayesian neural networks to produce diverse sets of predictions using a hybrid FEM-NN approach that approximates internal forces on a deforming solid body. The uncertainty estimator developed herein reliably infers upper bounds of bias/variance in the predictions for a wide range of interpolation and extrapolation cases using a three-element FEM-NN model of a bar undergoing plastic deformation. This proposed framework offers a powerful tool for assessing the reliability of physics-based surrogate models by establishing uncertainty estimates for predictions spanning a wide range of possible load cases.

42 ENGINEERING↗
DOE OSTI2022

Focused Ion Beam Preparation of Low Melting Point Metals: Lessons Learned From Indium

Indium (In) and other low melting point metals are used as interconnects in a variety of hybridized circuits and a full understanding of the metallurgy of these interconnects is important to the reliability and performance of the devices. Here, this paper shows that room temperature focused ion beam (FIB) preparation of cross-sections, using Ga + or Xe + can result in artifacts that obscure the true In microbump structure. The use of modified milling strategies to minimize the increased local sample temperature are shown to produce cross-sections that are representative of the In bump microstructure in some sample configurations. Furthermore, cooling of the sample to cryogenic temperatures is shown to reliably eliminate artifacts in FIB prepared cross-sections of In bumps allowing the true bump microstructure to be observed.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Trust Not Verify? The Critical Need for Data Curation Standards in Materials Informatics

The importance of data curation has been recognized in multiple areas of research; however, the discussion of this important issue is only beginning to emerge in materials science. In this Perspective, we highlight the benefits of using the standardized data curation protocols in materials science and discuss current gaps in accurate and reproducible data reporting using case studies drawn from high-impact materials science papers and well-known databases such as the Crystallography Open Database (COD) and the Cambridge Structural Database (CSD). We argue that both experimental and computational materials scientists need to embrace a culture of rigorous data curation as part of modern research data management. We propose a sample data curation pipeline for materials chemistry and illustrate its use by creating two new materials chemistry databases. Here, we hope that this perspective will serve to catalyze further discussion and promote the continuous development of rigorous data curation practices within the materials science research community. We posit that adherence to best practices of data curation will promote and enhance the reliability, reproducibility, and integrity of materials research and enable the development of reliable AI and machine learning models that critically depend on the use of quality data.

Chemical structure↗
DOE OSTI2025

Thin-Film Fracture Behavior for Diketopyrrolopyrrole Semiconducting Polymeric Films

Fracture energy, which quantifies a material’s resistance to the propagation of a pre-existing crack, is a key parameter for ensuring the mechanical reliability of stretchable organic electronic devices. However, most existing methods, such as a four-point bending fracture energy, utilized for measuring the fracture energy of semiconducting polymeric thin films are complicated by substrate effects, making it challenging to isolate the intrinsic behavior of the film from interfacial influences. In this study, we employed a pseudo free-standing pure shear method to systematically investigate the cohesive fracture energy of poly(diketopyrrolopyrrole-terthiophene) P(DPP-T)-based thin films to examine the effects of nanoconfinement, side chain length, degree of crystallinity, and strain rates. This method effectively eliminates substrate interference, enabling a direct assessment of the cohesive fracture energy of P(DPP-T) thin films. We found that thinner films and those with lower molecular weights exhibited significantly reduced fracture energies due to diminished chain entanglements. Additionally, films with shorter side chains displayed notably higher fracture energies, which were attributed to an increase in the degree of crystallinity. Finally, slower strain rates led to higher fracture energies, consistent with an enhanced stress relaxation. These insights offer practical guidelines for designing mechanically robust semiconducting polymers, contributing to the advancement of reliable, durable, flexible, and wearable electronic devices.

Zhang, Song [Univ. of Southern Mississippi, Hattie↗
DOE OSTI2026

Dara: Automated Multiple-Hypothesis Phase Identification and Refinement from Powder X-ray Diffraction

Powder X-ray diffraction (XRD) is a foundational technique for characterizing crystalline materials. However, the reliable interpretation of XRD patterns, particularly in multiphase systems, remains a manual and expertise-demanding task. As a characterization method that only provides structural information, multiple reference phases can often be fit to a single pattern, leading to potential misinterpretation when alternative solutions are overlooked. To ease humans’ efforts and address the challenge, we introduce Dara (data-driven automated Rietveld analysis), a framework designed to automate the robust identification and refinement of multiple phases from powder XRD data. Dara performs an exhaustive tree search over all plausible phase combinations within a given chemical space and validates each hypothesis using the BGMN Rietveld refinement routine. Key features include structural database filtering, automatic clustering of isostructural phases during tree expansion, and peak-matching-based scoring to identify promising phases for refinement. When ambiguity exists, Dara generates multiple hypothesis which can then be decided between by human experts or with further characterization tools. By enhancing the reliability and accuracy of phase identification, Dara enables scalable analysis of realistic complex XRD patterns and provides a foundation for integration into multimodal characterization workflows, moving toward fully self-driving materials discovery.

Biological databases↗
DOE OSTI2025

Improving Bond Dissociations of Reactive Machine Learning Potentials through Physics-Constrained Data Augmentation

In the field of computational chemistry, predicting bond dissociation energies (BDEs) presents well-known challenges, particularly due to the multireference character of reactive systems. Many chemical reactions involve configurations where single-reference methods fall short, as the electronic structure can significantly change during bond breaking. As generating training data for partially broken bonds is a challenging task, even state-of-the-art reactive machine learning interatomic potentials (MLIPs) often fail to predict reliable BDEs and smooth dissociation curves. By contrast, simple and inexpensive physics-based models, such as the well-established Morse potential, do not suffer from any such limitations. This work leverages the Morse potential to improve reactive MLIPs by augmenting the training data set with inexpensive Morse data along the dissociation pathways. Further, this physics-constrained data augmentation (PCDA) approach results in MLIPs with smooth bond dissociation curves as well as near coupled-cluster level BDEs, all without requiring any expensive multireference quantum mechanical calculations. A case study for methane combustion demonstrates how the PCDA approach can improve an existing reactive MLIP, namely, ANI-1xnr. In conclusion, not only are the BDEs and bond dissociation curves for all radicals and molecules significantly improved compared to ANI-1xnr but the PCDA-trained MLIP retains the reliability of ANI-1xnr when performing reactive molecular dynamics simulations.

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

MS25: Materials Science-Focused Benchmark Data Set for Machine Learning Interatomic Potentials

Here, we present MS25, a benchmark data set for evaluating machine learning interatomic potentials (MLIPs) across diverse materials-relevant systems including MgO surfaces, liquid water, zeolites, a catalytic Pt surface reaction, high-entropy alloys (HEAs), and disordered Zr-oxides. Five MLIP architectures (MACE, NequIP, Allegro, MTP, and Torch-ANI) are trained and tested, focusing not only on traditional metrics (energies, forces, and stresses) but also explicitly validating derived physical observables such as lattice constants, volumes, and reaction barriers. We find that most models reach comparable accuracy on standard error metrics across the simple systems, although equivariant MLIPs offer 1.5–2× improvements over nonequivariant MLIPs in energy and force error for structurally complex or compositionally disordered environments such as HEAs and Zr–O systems. Our analysis highlights that low errors in energy and force predictions do not guarantee reliable observables, emphasizing the necessity of explicit validation. We demonstrate limitations in cross-framework transferability, as models trained on one zeolite framework (CHA) fail to reliably generalize to predictions of structurally distinct frameworks (e.g., MFI). Size-extensive tests show some dependence on system size for MgO, resulting from forced periodicity. The HEA and Zr–O data sets are identified as challenging tests for future benchmarks and MLIP model architecture developments as they show significant differentiation in error between MLIP architectures and are still relatively difficult at 1000 training images. Moving forward, we recommend that benchmarking efforts shift their focus from marginal accuracy improvements in energy and force errors toward identifying and understanding model failure modes, rigorously assessing transferability, and evaluating how their errors affect observable predictions. For researchers looking to choose an MLIP architecture, we suggest selecting equivariant MLIP architectures if the complexity of the system is a challenge. For simple materials problems, auxiliary features such as integration with molecular dynamics engines, trade-offs between computational data set generation cost vs MLIP inference speed, and framework integration may play a more important decision factor than small differences in error metrics that are unlikely to matter for production-level research.

chemical structure↗
DOE OSTI2020

Utilizing Essential Symmetry Breaking in Auxiliary-Field Quantum Monte Carlo: Application to the Spin Gaps of the C 36 Fullerene and an Iron Porphyrin Model Complex

We present three distinct examples where phaseless auxiliary-field quantum Monte Carlo (ph-AFQMC) can be reliably performed with a single-determinant trial wave function with essential symmetry breaking. Essential symmetry breaking was first introduced by Lee and Head-Gordon [ Phys. Chem. Chem. Phys. 2019, 21, 4763-4778, 10.1039/C8CP07613H]. We utilized essential complex and time-reversal symmetry breaking with ph-AFQMC to compute the triplet–singlet energy gap in the TS12 set. We found statistically better performance of ph-AFQMC with complex-restricted orbitals than with spin-unrestricted orbitals. We then showed the utilization of essential spin symmetry breaking when computing the singlet–triplet gap of a known biradicaloid, C 36 . ph-AFQMC with spin-unrestricted Hartree–Fock (ph-AFQMC+UHF) fails catastrophically even with spin-projection and predicts no biradicaloid character. With approximate Brueckner orbitals obtained from regularized orbital-optimized second-order Møller–Plesset perturbation theory (κ-OOMP2), ph-AFQMC quantitatively captures strong biradicaloid character of C 36 . Lastly, we applied ph-AFQMC to the computation of the quintet–triplet gap in a model iron porphyrin complex where brute-force methods with a small active space fail to capture the triplet ground state. We show unambiguously that neither triplet nor quintet is strongly correlated using UHF, κ-OOMP2, and coupled-cluster with singles and doubles (CCSD) performed on UHF and κ-OOMP2 orbitals. There is no essential symmetry breaking in this problem. By virtue of this, we were able to perform UHF+ph-AFQMC reliably with a cc-pVTZ basis set and predicted a triplet ground state for this model geometry. The largest ph-AFQMC in this work correlated 186 electrons in 956 orbitals. Our work highlights the utility, scalability, and accuracy of ph-AFQMC with a single-determinant trial wave function with essential symmetry breaking for systems mainly dominated by dynamical correlation with little static correlation.

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

CHARMM-GUI Nanomaterial Modeler for Modeling and Simulation of Nanomaterial Systems

Molecular modeling and simulation are invaluable tools for nanoscience that predict mechanical, physicochemical, and thermodynamic properties of nanomaterials and provide molecular-level insight into underlying mechanisms. However, building nanomaterial-containing systems remains challenging due to the lack of reliable and integrated cyberinfrastructures. Here we present Nanomaterial Modeler in CHARMM-GUI, a web-based cyberinfrastructure that provides an automated process to generate various nanomaterial models, associated topologies, and configuration files to perform state-of- the-art molecular dynamics simulations using most simulation packages. The nanomaterial models are based on the interface force field, one of the most reliable force fields (FFs). The transferability of nanomaterial models among the simulation programs was assessed by single-point energy calculations, which yielded 0.01% relative absolute energy differences for various surface models and equilibrium nanoparticle shapes. Three widely used Lennard-Jones (LJ) cutoff methods are employed to evaluate the compatibility of nanomaterial models with respect to conventional biomolecular FFs: simple truncation at r = 12 Å (12 cutoff), force-based switching over 10 to 12 Å (10-12 fsw), and LJ particle mesh Ewald with no cutoff (LJPME). The FF parameters with these LJ cutoff methods are extensively validated by reproducing structural, interfacial, and mechanical properties. As a result, we find that the computed density and surface energies are in good agreement with reported experimental results, although the simulation results increase in the following order: 10-12 fsw <12 cutoff < LJPME. Nanomaterials in which LJ interactions are a major component show relatively higher deviations (up to 4% in density and 8% in surface energy differences) compared with the experiment. Nanomaterial Modeler's capability is also demonstrated by generating complex systems of nanomaterial-biomolecule and nanomaterial-polymer interfaces with a combination of existing CHARMM-GUI modules. We hope that Nanomaterial Modeler can be used to carry out innovative nanomaterial modeling and simulations to acquire insight into the structure, dynamics, and underlying mechanisms of complex nanomaterial-containing systems.

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

Asynchronous Reciprocal Coupling of Martini 2.2 Coarse-Grained and CHARMM36 All-Atom Simulations in an Automated Multiscale Framework

The appeal of multiscale modeling approaches is predicated on the promise of combinatorial synergy. However, this promise can only be realized when distinct scales are combined with reciprocal consistency. Here, we consider multiscale molecular dynamics (MD) simulations that combine the accuracy and macromolecular flexibility accessible to fixed-charge all-atom (AA) representations with the sampling speed accessible to reductive, coarse-grained (CG) representations. AA-to-CG conversions are relatively straightforward because deterministic routines with unique outcomes are achievable. Conversely, CG-to-AA conversions have many solutions due to a surge in the number of degrees of freedom. While automated tools for biomolecular CG-to-AA transformation exist, we find that one popular option, called Backward, is prone to stochastic failure and the AA models that it does generate frequently have compromised protein structure and incorrect stereochemistry. Although these shortcomings can likely be circumvented by human intervention in isolated instances, automated multiscale coupling requires reliable and robust scale conversion. In this work, we detail an extension to Multiscale Machine-learned Modeling Infrastructure (MuMMI), including an improved CG-to-AA conversion tool called sinceCG. This tool is reliable (~98% weakly correlated repeat success rate), automatable (no unrecoverable hangs), and yields AA models that generally preserve protein secondary structure and maintain correct stereochemistry. We describe how the MuMMI framework identifies CG system configurations of interest, converts them to AA representations, and simulates them at the AA scale while on-the-fly analyses provide feedback to update CG parameters. Application to systems containing the peripheral membrane protein RAS and proximal components of RAF kinase on complex eight-component lipid bilayers with ~1.5 million atoms is discussed in the context of MuMMI.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗