Search NASA⌕ Search

SEARCH · Search NASA

Results for “FLEXIBLE BODY”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Toward simulating quantum field theories with controlled phonon-ion dynamics: A hybrid analog-digital approach

Quantum field theories are the cornerstones of modern physics, providing relativistic and quantum mechanical descriptions of physical systems at the most fundamental level. Simulating real-time dynamics within these theories remains elusive in classical computing. This provides a unique opportunity for quantum simulators, which hold the promise of revolutionizing our simulation capabilities. Trapped-ion systems are successful quantum-simulator platforms for quantum many-body physics and can operate in digital, or gate-based, and analog modes. Inspired by the progress in proposing and realizing quantum simulations of a number of relativistic quantum field theories using trapped-ion systems, and by the hybrid analog-digital proposals for simulating interacting boson-fermion models, we propose hybrid analog-digital quantum simulations of selected quantum field theories, taking recent developments to the next level. On one hand, the semi-digital nature of this proposal offers more flexibility in engineering generic model interactions compared with a fully-analog approach. On the other hand, encoding the bosonic fields onto the phonon degrees of freedom of the trapped-ion system allows a more efficient usage of simulator resources, and a more natural implementation of intrinsic quantum operations in such platforms. This opens up new ways for simulating complex dynamics of e.g., Abelian and non-Abelian gauge theories, by combining the benefits of digital and analog schemes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Comparing computational times for simulations when using PBPK model template and stand-alone implementations of PBPK models

Introduction We previously developed a PBPK model template that consists of a single model “superstructure” with equations and logic found in many physiologically based pharmacokinetic (PBPK) models. Using the template, one can implement PBPK models with different combinations of structures and features. Methods To identify factors that influence computational time required for PBPK model simulations, we conducted timing experiments using various implementations of PBPK models for dichloromethane and chloroform, including template and stand-alone implementations, and simulating four different exposure scenarios. For each experiment, we measured the required computational time and evaluated the impacts of including various model features (e.g., number of output variables calculated) and incorporating various design choices (e.g., different methods for estimating blood concentrations). Results We observed that model implementations that treat body weight and dependent quantities as constant (fixed) parameters can result in a 30% time savings compared with options that treat body weight and dependent quantities as time-varying. We also observed that decreasing the number of state variables by 36% in our PBPK model template led to a decrease of 20–35% in computational time. Other factors, such as the number of output variables, the method for implementing conditional statements, and the method for estimating blood concentrations, did not have large impacts on simulation time. In general, simulations with PBPK model template implementations of models required more time than simulations with stand-alone implementations, but the flexibility and (human) time savings in preparing and reviewing a model implemented using the PBPK model template may justify the increases in computational time requirements. Conclusion Our findings concerning how PBPK model design and implementation decisions impact computational speed can benefit anyone seeking to develop, improve, or apply a PBPK model, with or without the PBPK model template.

Bernstein, Amanda S.↗
DOE OSTI2026

Ferrimagnetism of ultracold fermions in a multiband Hubbard system

Strongly correlated materials feature multiple electronic orbitals, which are crucial to accurately understanding their many-body properties. In such multiband models, quantum interference can lead to flat energy bands with large degeneracy that gives rise to itinerant magnetic phases. Here, we report on signatures of a ferrimagnetic state realized in a Lieb lattice with ultracold fermions, characterized by antialigned magnetic moments with antiferromagnetic correlations, and concomitant with a finite spin polarization. The signatures remain robust when increasing repulsive interactions from the weakly interacting to the Heisenberg regime and emerge when continuously tuning the lattice unit cell from a square to a Lieb geometry. Our flexible approach paves the way toward exploring exotic phases, such as quantum spin liquids in kagome lattices and heavy fermion behavior in Kondo models.

Lebrat, Martin [Harvard Univ., Cambridge, MA (Unit↗
DOE OSTI2026

Advancing density functional tight-binding method for large organic molecules through equivariant neural networks

Semi-empirical quantum-mechanical (QM) methods have become valuable tools for studying complex (bio)molecular systems due to their balance between computational efficiency and accuracy. A key aspect of these methods is their parameterization, which not only governs the reliability of the results but also provides an opportunity to enhance their overall performance. In our previous work [J. Phys. Chem. Lett., 2021, 11, 16], we advanced the third-order semi-empirical density functional tight-binding (DFTB3) method for computing multiple properties of small molecules by developing the machine learning (ML) potential NN rep to bridge the gap between DFTB3 electronic components and those of the hybrid DFT-PBE0 functional. To overcome the limitations of NN rep , we introduce the EquiDTB framework, which leverages physics-inspired equivariant neural networks (NN) to parameterize scalable and transferable many-body Δ TB potentials, replacing the standard pairwise DFTB repulsive potential. This advancement extends the applicability of our ML-corrected DFTB approach to larger molecules and non-covalent systems (including only C, N, O, and H atoms), going beyond the chemical space represented in the training QM datasets. The enhanced performance of EquiDTB over the standard TB methods is demonstrated by the accurate computation of the atomic forces of S66x8 molecular dimers, as well as their interaction energies. Moreover, EquiDTB can be effectively employed to explore the potential energy surfaces of large and flexible drug-like molecules—for example, to determine the minimum energy path between isomers, analyze structural transitions during dynamical simulations, compute vibrational modes, and investigate energetic rankings. The performance for single molecules slightly decreases when the DFTB electronic energy is reduced to first-order but remains superior to standard TB methods. Our work thus demonstrates that an optimal integration of an equivariant NN with QM datasets can advance the DFTB method while maintaining high efficiency, paving the way for reliable (bio)molecular simulations.

Medrano Sandonas, Leonardo [Technische Universität↗
DOE OSTI2021

Vapor–liquid equilibrium of water with the MB-pol many-body potential

Among the many existing molecular models of water, the MB-pol many-body potential has emerged as a remarkably accurate model, capable of reproducing thermodynamic, structural, and dynamic properties across water's solid, liquid, and vapor phases. In this work, we assessed the performance of MB-pol with respect to an important set of properties related to vapor-liquid coexistence and interfacial behavior. Through direct coexistence classical molecular dynamics simulations at temperatures of 400 K < T < 600 K, we calculated properties such as equilibrium coexistence densities, vapor-liquid interfacial tension, vapor pressure, and enthalpy of vaporization and compared the MB-pol results to experimental data. We also compared rigid vs fully flexible variants of the MB-pol model and evaluated system size effects for the properties studied. We found that the MB-pol model predictions are in good agreement with experimental data, even for temperatures approaching the vapor-liquid critical point; this agreement was largely insensitive to system sizes or the rigid vs flexible treatment of the intramolecular degrees of freedom. Furthermore, these results attest to the chemical accuracy of MB-pol and its high degree of transferability, thus enabling MB-pol's application across a large swath of water's phase diagram.

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

Changes in polarization dictate necessary approximations for modeling electronic deexcitation intensity: Application to x-ray emission

Accurate simulation of electronic excitations and deexcitations are critical for complementing complex spectroscopic experiments and can provide validation to theoretical approaches. Here, using a generalized framework, we contrast the accuracy and validity of orbital-constrained and linear-response approaches that build upon Kohn-Sham density functional theory (DFT) to simulate emission spectra of electronic origin and propose an efficient approximation, named many-body x-ray emission spectroscopy or MBXES, for simulating such processes. We show analytically as well as with computed examples that for electronic (de)excitation leading to an appreciable change in polarization (i.e., density rearrangement), the adiabatic approximation in a response-based formalism will be inadequate for the calculation of oscillator strength. Thus, such a change (e.g., in the net electrostatic dipole moment of a finite system) can be used as a metric for evaluating the applicability of the adiabatic response-based approach and can be particularly valuable in x-ray emission spectroscopy. On the other hand, MBXES, the flexible method introduced in this paper, can compute oscillator strengths accurately at a much lower computational expense on the basis of two DFT-based self-consistent field calculations. Using illustrative examples of emission spectra, the efficacy of the MBXES method is demonstrated by comparison with its parent theory, orbital-optimized DFT, and with experiments.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

A Lizard-inspired Tube Inspector (LTI) Robot

This final report summarizes the research findings of the current project. This project is a collaboration between New Mexico State University (NMSU) as a lead (Recipient) and Arizona State University (ASU) as a Co-Recipient. Tubular structures are common components of boilers and heat exchangers in power plants. Over time, these components may suffer corrosion, cracks, and stress-corrosion cracks in either the body or the welded connections. A single tube leakage can cause an outage of several weeks. Regular inspection is a key safety factor when ensuring that power plants are maintained in reliable, operational condition. This inspection, however, is challenging, time-consuming, and in many cases impossible due to accessibility issues and safety concerns. Recent developments in robotic-based inspection can offer a great solution. Hard-to-reach places can be inspected without overhauling the unit, saving time and cost. Although several robots have been designed and implemented successfully for inspecting power plant components, in particular tubular structures, their mobility and flexibility are limited. Most of these robots use wheels for mobility which reduces their maneuverability of these robots. Moreover, they usually use magnets to attach to tubes which will not work on tubes with non-ferromagnetic materials. These robots usually carry measurement tools for a particular non-destructive testing (NDT) method such as ultrasound testing (UT) that requires a couplant to perform a point-by-point (scanning) inspection of the tubular structure. In this project, we developed a versatile lizard-inspired tube inspector (LTI) robot with embedded inspection sensing components for tube inspection which removes the need for point-by-point scanning of tube surface for crack and corrosion detection. Inspired by a “lizard”, the novelty of the current project is the integration of couplant-free ultrasound sensing and transmission, advanced data-driven defect detection and imaging, and friction-based mechanical mobility components in a single robot to eliminate a need for smooth surfaces and simple geometry for mobility and scanning. The LTI robot could replace the wheel-based approach with friction-based mobility to significantly increase the flexibility and maneuverability of the robot. The LTI robot can get into a power plant unit, such as a boiler, from a small area allowing it to access a component of interest for inspection (e.g., move on OD, curved and flat surfaces, non-ferromagnetic or ferromagnetic materials, and tubes with rough surfaces and complex geometries). Additionally, an advanced data-driven method using ultrasound data was pioneered to allow the robot to detect defects in the entire area between the robot’s multi-functional mobility system (grippers). Integrating a couplant-free ultrasound sensing in the robot’s grippers as well as using advanced data-driven methods allowed the LTI to detect defects in the entire cross-section of a tube using its grippers when stationary and when mobile.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2021

Creation of a prototype biomimetic fish to better understand impact trauma caused by hydropower turbine blade strikes

Biomimetic model organisms could be useful surrogates for live animals in many applications if the models have sufficient biofidelity. One such application is for use in field and laboratory tests of fish mortality associated with passage through hydropower turbines. Laboratory trials suggest that blade strikes are especially injurious and often causes mortality when fish are struck by thinner blades moving at higher velocities. Dose-response relationships have been created from these data, but the exact relationship between fish mortality and the actual forces enacted on fish during simulated blade strike testing remains unknown. Here, we describe the methods used to create a prototype biomimetic model fish composed of ballistic gelatin and covered with a surrogate skin to better approximate the biomechanical properties of a fish body. Frozen fish were scanned with high-fidelity laser scanners, and a 3D-printed, reusable mold was created from which to cast our gelatin model. Computed tomography scan data, imaged directly or taken from online data repositories, were also successfully used to create CAD models for use in additive manufacturing of molds. One 3-axis accelerometer was embedded into the gelatin to compare accelerometer data to dose-response data from previous laboratory research on live fish. The resulting model ( i.e. , Gelfish) had a statistically indistinguishable tissue durometer to that of real fish tissue and preliminary blade strike impact testing suggested its overall flexibility was similar to that of live fish. Gelfish was designed with biofidelity as its guiding principle and our results suggest initial experimentation was successful. Future research will include replication of initial Gelfish test results, quantitative measurement of model flexibility relative to real fish, and inclusion of surrogate skeletal structures to enhance biofidelity. Use of more sophisticated sensors would also better quantify the physical forces of blade strike impact and help determine how said forces correlate with rates of mortality observed during tests on live fish.

36 MATERIALS SCIENCE↗
DOE OSTI2023

The Synchronic Web

The Synchronic Web is a distributed network for securing data provenance on the World Wide Web. By enabling clients around the world to freely commit digital information into a single shared view of history, it provides a foundational basis of truth on which to build decentralized and scalable trust across the Internet. Its core cryptographical capability allows mutually distrusting parties to create and verify statements of the following form: “I commit to this information—and only this information—at this moment in time.” The backbone of the Synchronic Web infrastructure is a simple, small, and semantic-free blockchain that is accessible to any Internet-enabled entity. The infrastructure is maintained by a permissioned network of well-known servers, called notaries, and accessed by a permissionless group of clients, called journals. Through an evolving stack of flexible and composable semantic specifications, the parties cooperate to generate synchronic commitments over arbitrary data. When integrated with existing infrastructures, adapted to diverse domains, and scaled across the breadth of cyberspace, the Synchronic Web provides a ubiquitous mechanism to lock the world’s data into unique points in discrete time and digital space. This document provides a technical description of the core Synchronic Web system. The distinguishing innovation in our design—and the enabling mechanism behind the model—is the novel use of verifiable maps to place authenticated content into canonically defined locations off-chain. While concrete specifications and software implementations of the Synchronic Web continue to evolve, the information covered in the body of this document should remain stable. We aim to present this information clearly and concisely for technical non-experts to understand the essential functionality and value proposition of the network. In the interest of promoting discourse, we take some liberty in projecting the potential implications of the new model.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Computational studies of impurity migration during induction stirring of molten uranium

Understanding and controlling impurity behavior in actinide metal casting processes are foundational for efficient part production yet remain major challenges for researchers and industry. To help provide insight regarding impurity distribution during actinide metal casting, we have developed computational fluid dynamics (CFD) models for a laboratory-scale system using commercial and open-source code. Multiple simulation frameworks allow for improved confidence in the resulting outputs, while taking advantage of the maturity and convenience offered by commercial platforms and simultaneously maintaining the transparency and flexibility often provided by open-source software. Here, we describe multiple experiment-informed models designed to simulate a specific laboratory system in which uranium melt, containing a known starting concentration of carbon impurity, is electromagnetically stirred in an induction furnace. Here, the goal of the simulations is to predict the motion of uranium carbide microparticles in the velocity field of the melt. Prior to simulating the uranium-carbon system, numerical models were validated using a previously published nonradioactive experimental system. Effects of the size and shape of impurity particles in the models were investigated and agree with experimental findings. Simulation of smaller particles (< 50 µm) shows more homogenous distribution throughout the stirred melt. With increased particle size (100 µm), the body forces, which include the buoyancy force, dominate over the drag force, causing larger particles to move toward the crucible walls and upward in the system.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Exploring finite temperature properties of materials with quantum computers

Abstract Thermal properties of nanomaterials are crucial to not only improving our fundamental understanding of condensed matter systems, but also to developing novel materials for applications spanning research and industry. Since quantum effects arise at the nano-scale, these systems are difficult to simulate on classical computers. Quantum computers can efficiently simulate quantum many-body systems, yet current quantum algorithms for calculating thermal properties of these systems incur significant computational costs in that they either prepare the full thermal state on the quantum computer, or they must sample a number of pure states from a distribution that grows with system size. Canonical thermal pure quantum (TPQ) states provide a promising path to estimating thermal properties of quantum materials as they neither require preparation of the full thermal state nor require a growing number of samples with system size. Here, we present an algorithm for preparing canonical TPQ states on quantum computers. We compare three different circuit implementations for the algorithm and demonstrate their capabilities in estimating thermal properties of quantum materials. Due to its increasing accuracy with system size and flexibility in implementation, we anticipate that this method will enable finite temperature explorations of relevant quantum materials on near-term quantum computers.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Progress from ASDEX Upgrade experiments in preparing the physics basis of ITER operation and DEMO scenario development

An overview of recent results obtained at the tokamak ASDEX Upgrade (AUG) is given. A work flow for predictive profile modelling of AUG discharges was established which is able to reproduce experimental H-mode plasma profiles based on engineering parameters only. In the plasma center, theoretical predictions on plasma current redistribution by a dynamo effect were confirmed experimentally. For core transport, the stabilizing effect of fast ion distributions on turbulent transport is shown to be important to explain the core isotope effect and improves the description of hollow low-Z impurity profiles. The L–H power threshold of hydrogen plasmas is not affected by small helium admixtures and it increases continuously from the deuterium to the hydrogen level when the hydrogen concentration is raised from 0 to 100%. One focus of recent campaigns was the search for a fusion relevant integrated plasma scenario without large edge localised modes (ELMs). Results from six different ELM-free confinement regimes are compared with respect to reactor relevance: ELM suppression by magnetic perturbation coils could be attributed to toroidally asymmetric turbulent fluctuations in the vicinity of the separatrix. Stable improved confinement mode plasma phases with a detached inner divertor were obtained using a feedback control of the plasma β. The enhanced D α H-mode regime was extended to higher heating power by feedback controlled radiative cooling with argon. The quasi-coherent exhaust regime was developed into an integrated scenario at high heating power and energy confinement, with a detached divertor and without large ELMs. Small ELMs close to the separatrix lead to peeling-ballooning stability and quasi continuous power exhaust. Helium beam density fluctuation measurements confirm that transport close to the separatrix is important to achieve the different ELM-free regimes. Based on separatrix plasma parameters and interchange-drift-Alfvén turbulence, an analytic model was derived that reproduces the experimentally found important operational boundaries of the density limit and between L- and H-mode confinement. Feedback control for the X-point radiator (XPR) position was established as an important element for divertor detachment control. Stable and detached ELM-free phases with H-mode confinement quality were obtained when the XPR was moved 10 cm above the X-point. Investigations of the plasma in the future flexible snow-flake divertor of AUG by means of first SOLPS-ITER simulations with drifts activated predict beneficial detachment properties and the activation of an additional strike point by the drifts.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Impact of Control on Availability and Cycling of Residential HVAC in a Real-World Experiment

Demand response is an important emerging part of smart grids and there is a large stream of research from theoretical and modeling perspectives. However, there is relatively little experimental evidence that could help researchers and building operators make informed decisions on best practices for modeling, developing, and deployment of control mechanisms. We contribute to the body of experimental research by providing numerical insights into the role and availability of residential HVAC systems for control. We share the findings on the duration of the cooling cycle, off cycle, and temperature settling time of HVAC systems from data collected from a smart neighborhood located in Atlanta, GA.

cycling↗
DOE OSTI2021

Variational Discrete Action Theory

In this work, we propose the variational discrete action theory (VDAT) to study the ground state properties of quantum many-body Hamiltonians. VDAT is a variational theory based on the sequential product density matrix (SPD) ansatz, characterized by an integer $\mathscr{N}$, which monotonically approaches the exact solution with increasing $\mathscr{N}$. To evaluate the SPD, we introduce a discrete action and a corresponding integer time Green’s function. We use VDAT to exactly evaluate the SPD in two canonical models of interacting electrons: the Anderson impurity model and the d = ∞ Hubbard model. For the latter, we evaluate $\mathscr{N}$ = 2 – 4, where $\mathscr{N}$ = 2 recovers the Gutzwiller approximation (GA), and we show that $\mathscr{N}$ = 3, which exactly evaluates the Gutzwiller-Baeriswyl wave function, provides a truly minimal yet precise description of Mott physics with a cost similar to that of the GA. VDAT is a flexible theory for studying quantum Hamiltonians, competing both with state-of-the-art methods and simple, efficient approaches all within a single framework.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Matilda v1.0: An R package for probabilistic climate projections using a reduced complexity climate model

A primary advantage to using reduced complexity climate models (RCMs) has been their ability to quickly conduct probabilistic climate projections, a key component of uncertainty quantification in many impact studies and multisector systems. Providing frameworks for such analyses has been a target of several RCMs used in studies of the future co-evolution of the human and Earth systems. In this paper, we present Matilda, an open-science R software package that facilitates probabilistic climate projection analysis, implemented here using the Hector simple climate model in a seamless and easily applied framework. The primary goal of Matilda is to provide the user with a turn-key method to build parameter sets from literature-based prior distributions, run Hector iteratively to produce perturbed parameter ensembles (PPEs), weight ensembles for realism against observed historical climate data, and compute probabilistic projections for different climate variables. This workflow gives the user the ability to explore viable parameter space and propagate uncertainty to model ensembles with just a few lines of code. The package provides significant freedom to select different scoring criteria and algorithms to weight ensemble members, as well as the flexibility to implement custom criteria. Additionally, the architecture of the package simplifies the process of building and analyzing PPEs without requiring significant programming expertise, to accommodate diverse use cases. We present a case study that provides illustrative results of a probabilistic analysis of mean global surface temperature as an example of the software application.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Meter-Based Assessment of the Time and Locational Benefits of a Large Utility’s DSM Portfolio

As decarbonization goals drive increasing levels of renewable generation, there is a need to understand the time- and location-based savings benefits of demand-side management (DSM) programs. The challenges of the 'duck curve' are driving the utility industry to consider how programs can be optimized to match demand profiles with low carbon generation resources. From an infrastructure standpoint, time- and location-targeted DSM could serve as a ‘non-wires alternative’ (NWA) to defer equipment upgrades. Additional DSM value streams are motivating innovation in savings evaluation, providing more resolved insights beyond the total annual program impact. Methods grounded in the principles of billing analysis, leveraging hourly metering at the distribution grid, can provide new visibility into the spatial and temporal savings achieved through DSM. A large body of work has investigated related topics including interval meter-based savings analysis, the time- varying nature of efficiency measures, and NWA. A less studied topic concerns the impact of DSM on the grid, based on metered consumption. This paper presents an analysis of interval data across more than 25,000 customers and twelve substations, from the Sacramento Municipal Utility District. The results show for different locations on the grid: achieved savings and the impact on grid consumption; hourly savings shapes for DSM program participants and non-participants, and how those shapes vary with season; and the impact of the programs on peak demand. These findings show the current impact of DSM, with implications for future, more intentional targeting as the utility continues to pursue aggressive electrification, efficiency, load flexibility, and reliable NWA.

Granderson, Jessica↗
DOE OSTI2020

A Thermal Radiation Modulation Platform by Emissivity Engineering with Graded Metal–Insulator Transition

Thermal radiation from a black body increases with the fourth power of absolute temperature (T 4 ), an effect known as the Stefan-Boltzmann law. Typical materials radiate heat at a portion of this limit, where the portion, called integrated emissivity (ε int ), is insensitive to temperature (|dε int /dT| ≈ 10 -4 °C -1 ). The resultant radiance bound by the T 4 law limits the ability to regulate radiative heat. Here, an unusual material platform is shown in which ε int can be engineered to decrease in an arbitrary manner near room temperature (|dε int /dT| ≈ 8 × 10 -3 °C -1 ), enabling unprecedented manipulation of infrared radiation. As an example, ε int is programmed to vary with temperature as the inverse of T 4 , precisely counteracting the T 4 dependence; hence, thermal radiance from the surface becomes temperature-independent, allowing the fabrication of flexible and power-free infrared camouflage with unique advantage in performance stability. In this work, the structure is based on thin films of tungsten-doped vanadium dioxide where the tungsten fraction is judiciously graded across a thickness less than the skin depth of electromagnetic screening.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Permutationally Invariant Polynomial Expansions with Unrestricted Complexity

A general strategy is presented for constructing and validating permutationally invariant polynomial (PIP) expansions for chemical systems of any stoichiometry. Demonstrations are made for three categories of gas-phase dynamics and kinetics: collisional energy-transfer trajectories for predicting pressure-dependent kinetics, three-body collisions for describing transient van der Waals adducts relevant to atmospheric chemistry, and nonthermal reactivity via quasiclassical trajectories. In total, 30 systems are considered with up to 15 atoms and 39 degrees of freedom. Permutational invariance is enforced in PIP expansions with as many as 13 million terms and 13 permutationally distinct atom types by taking advantage of petascale computational resources. The quality of the PIP expansions is demonstrated through the systematic convergence of in-sample and out-of-sample errors with respect to both the number of training data and the order of the expansion, and these errors are shown to predict errors in the dynamics for both reactive and nonreactive applications. Here, the parallelized code distributed as part of this work enables the automation of PIP generation for complex systems with multiple channels and flexible user-defined symmetry constraints and for automatically removing unphysical unconnected terms from the basis set expansions, all of which are required for simulating complex reactive systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗