Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Inference mechanisms”

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 433 records · Page 24

Limited impact of jet-induced feedback in the multi-phase nuclear interstellar medium of 4C12.50

4C12.50 (IRAS 13451+1232) at z = 0.122 is an ultraluminous infared radio galaxy that has often been proposed as a prime candidate for the link between ultraluminous infared galaxies and young radio galaxies. It is also an interesting target to investigate whether and how radio-induced feedback affects the evolution of galaxies in the early phases of radio activity. We study, in detail for the first time, the hot (≥1500 K) molecular gas in 4C12.50. The potential impact of the radio jet on this gas phase, as well as on the star formation activity, are investigated. We also study the ionised (including coronal) gas as traced by the near-infrared lines. Using near-infrared long slit spectroscopy obtained with EMIR on GTC and X-shooter on VLT, we analyse the emission line spectrum of the ionised, coronal, and, specially, the hot molecular gas in the western nucleus hosting the compact radio jet. Based on high spatial resolution ALMA CO(2–1) data, we also revise the location of 4C12.50 in the Kennicutt-Schmidt diagram in order to investigate whether star formation is suppressed. 4C12.50 hosts (2.1 ± 0.4) × 10 4 M ⊙ of hot molecular gas. An unusually high rotational temperature T rot = 3020 ± 160 K is inferred. The molecular gas mass obeys a power-law temperature distribution, dM H2 /dT ∝ T –5 , from T ~300 K and up to ~3000 K. Both results support the idea that shocks (probably induced by the radio jet) contribute to the heating and excitation of the hot molecular gas. A molecular outflow is not detected. The coupling of the outflowing ionised and neutral outflows with the hot molecular gas is poor. Contrary to other studies, we claim that there is no evidence for star formation suppression in this object. If radio-induced feedback can regulate the star formation activity in galaxies, 4C12.50 is a promising candidate to reveal this phenomenon in action. However, we find no solid evidence for a current or past impact of this mechanism on the evolution of this system, neither by clearing out the dusty central cocoon efficiently, nor by suppressing the star formation activity.

79 ASTRONOMY AND ASTROPHYSICS↗

Investigating biological nitrogen fixation via single-cell transcriptomics

The extensive use of nitrogen fertilizers has detrimental environmental consequences, and it is essential for society to explore sustainable alternatives. One promising avenue is engineering root nodule symbiosis, a naturally occurring process in certain plant species within the nitrogen-fixing clade, into non-leguminous crops. Advancements in single-cell transcriptomics provide unprecedented opportunities to dissect the molecular mechanisms underlying root nodule symbiosis at the cellular level. This review summarizes key findings from single-cell studies in Medicago truncatula, Lotus japonicus, and Glycine max. We highlight how these studies address fundamental questions about the development of root nodule symbiosis, including the following findings: (i) single-cell transcriptomics has revealed a conserved transcriptional program in root hair and cortical cells during rhizobial infection, suggesting a common infection pathway across legume species; (ii) characterization of determinate and indeterminate nodules using single-cell technologies supports the compartmentalization of nitrogen fixation, assimilation, and transport into distinct cell populations; (iii) single-cell transcriptomics data have enabled the identification of novel root nodule symbiosis genes and provided new approaches for prioritizing candidate genes for functional characterization; and (iv) trajectory inference and RNA velocity analyses of single-cell transcriptomics data have allowed the reconstruction of cellular lineages and dynamic transcriptional states during root nodule symbiosis.

Lotus japonicus↗

Climate forcing controls on carbon terrestrial fluxes during shale weathering

Climate influences near-surface biogeochemical processes and thereby determines the partitioning of carbon dioxide (CO 2 ) in shale, and yet the controls on carbon (C) weathering fluxes remain poorly constrained. Using a dataset that characterizes biogeochemical responses to climate forcing in shale regolith, we implement a numerical model that describes the effects of water infiltration events, gas exchange, and temperature fluctuations on soil respiration and mineral weathering at a seasonal timescale. Our modeling approach allows us to quantitatively disentangle the controls of transient climate forcing and biogeochemical mechanisms on C partitioning. We find that ~3% of soil CO 2 (1.02 mol C/m 2 /y) is exported to the subsurface during large infiltration events. Here, net atmospheric CO 2 drawdown primarily occurs during spring snowmelt, governs the aqueous C exports (61%), and exceeds the CO 2 flux generated by pyrite and petrogenic organic matter oxidation (~0.2 mol C/m 2 /y). We show that shale CO 2 consumption results from the temporal coupling between soil microbial respiration and carbonate weathering. This coupling is driven by the impacts of hydrologic fluctuations on fresh organic matter availability and CO 2 transport to the weathering front. Diffusion-limited transport of gases under transient hydrological conditions exerts an important control on CO 2 (g) egress patterns and thus must be considered when inferring soil CO 2 drawdown from the gas phase composition. Our findings emphasize the importance of seasonal climate forcing in shaping the net contribution of shale weathering to terrestrial C fluxes and suggest that warmer conditions could reduce the potential for shale weathering to act as a CO 2 sink.

58 GEOSCIENCES↗

Organ-delimited gene regulatory networks provide high accuracy in candidate transcription factor selection across diverse processes

Organ-specific gene expression datasets that include hundreds to thousands of experiments allow the reconstruction of organ-level gene regulatory networks (GRNs). However, creating such datasets is greatly hampered by the requirements of extensive and tedious manual curation. Here, we trained a supervised classification model that can accurately classify the organ-of-origin for a plant transcriptome. This K-Nearest Neighbor-based multiclass classifier was used to create organ-specific gene expression datasets for the leaf, root, shoot, flower, and seed in Arabidopsis thaliana . A GRN inference approach was used to determine the: i. influential transcription factors (TFs) in each organ and, ii. most influential TFs for specific biological processes in that organ. These genome-wide, organ-delimited GRNs (OD-GRNs), recalled many known regulators of organ development and processes operating in those organs. Importantly, many previously unknown TF regulators were uncovered as potential regulators of these processes. As a proof-of-concept, we focused on experimentally validating the predicted TF regulators of lipid biosynthesis in seeds, an important food and biofuel trait. Of the top 20 predicted TFs, eight are known regulators of seed oil content, e.g., WRI1, LEC1, FUS3. Importantly, we validated our prediction of MybS2, TGA4, SPL12, AGL18, and DiV2 as regulators of seed lipid biosynthesis. We elucidated the molecular mechanism of MybS2 and show that it induces purple acid phosphatase family genes and lipid synthesis genes to enhance seed lipid content. This general approach has the potential to be extended to any species with sufficiently large gene expression datasets to find unique regulators of any trait-of-interest.

09 BIOMASS FUELS↗

Selective Amnesia using Contrastive Subnet Erasure for Class Level Unlearning in Vision Models

We study concept-level forgetting in pretrained vision models: removing an entire semantic category so the system no longer recognizes that object in unseen images and contexts, rather than merely forgetting specific training examples. Prior work either applies blunt global projections or fine-tunes parameters, which can introduce collateral damage to unrelated features, add compute, and become unstable as forgetting strength increases. We introduce Contrastive Subnet Erasure (CSE), a training-free, encoder-centric edit that targets a compact set of channels most responsible for the class and attenuates them in a calibrated manner. The modification is algebraically folded into the subsequent layer, yielding no inference-time overhead and leaving task heads unchanged. To evaluate whether forgetting generalizes beyond the data used to specify the class, we introduce a cross dataset protocol in which the class is defined on a source dataset and performance is measured on a disjoint target dataset drawn from a different distribution with no shared images. This setup tests whether the model still fails to recognize the object when it looks different or appears in new scenes, and it helps avoid overfitting to patterns in the source dataset. Across CIFAR 10, CIFAR 100, and ImageNet under this protocol, CSE achieves stronger forgetting of the target class while better preserving non target utility than existing baselines in both single class and multi class settings. Overall, CSE provides a simple, stable, and deployment-ready mechanism for class-level unlearning in vision.

Kotevska, Olivera [ORNL] (ORCID:0000000316772243)↗

Support for Reactor Operators in Case of Cyber-Security Threats (NEUP Final Report)

Increasing use of digital technology in nuclear power plants (NPPs) poses cyber-security as a crucial threat to public safety and to continuous energy production. Cyber-security risks are comprised of complex known and unknown interactions between various entities, system vulnerabilities, network protocols, human users and malicious attacks. There is little understanding or research geared towards plant operators’ response under cyber-security threats and operational procedures to cope with such threats. This is particularly critical when a cyber-security event masquerades as a safety incident or an evolving accident. Instead of leading operators to the remediation of the accident, the masked cyber event may lead them to circumvent the safety mechanisms of the plant. Given these challenges, the objective of this project is to develop a tool to characterize abnormal NPP events as “cyber” or “safety” incidents, and to further develop this tool as a real-time operator aid to assist in response to the incident. This NEUP Final Report summarizes the research efforts towards the above goals and the results in the following three aspects. First, a methodology for classifying abnormal NPP events as cyber-attacks or safety incidents was developed. This provides the basis for responding to potential cyber-attacks. In this methodology, the causal relationships between various variables (i.e., physical variables, network variables, and control variables) for a nuclear system under consideration are modeled using dynamic Bayesian networks. Based on observations of certain variables (i.e., the evidence), inference is performed to identify the state of hidden variables, which helps indicate the occurrence of cyber-attacks or safety incidents. The proposed methodology was applied to the full scope NPP simulator system at The Ohio State University as a case study. The result shows that the proposed methodology is capable of distinguishing between cyber-attacks and safety incidents. Second, a game-theoretic approach for responding to cyber-attacks on nuclear power plants was developed. This is the second major goal of the project, and the method is used to recommend optimal responses to reactor operators in the case of cyber-attacks. We consider two cases in this research. In the first case, the interaction between the defender (i.e., the reactor operator) and the attacker is modeled as a two-player, nonzero-sum, Markov game with an infinite time horizon. In the second case, we consider a semi-Markov game between the defender and the attacker within a finite time horizon. The model in the second case enables us to consider the time sensitive response of the defender. In each case, the method for deriving the Nash equilibrium is developed. The Nash Equilibrium of the game provides the valid prediction of both players’ actions because no single player can benefit from unilaterally deviating from the equilibrium policy if the other player adheres to his/hers, hence it provides the best response of the defender to cyber-attacks. The proposed method in each case is demonstrated through a simplified benchmark digital feedwater control system. Third, a prototype display is developed to integrate the above information. The prototype display depicts the communication information to the operator and incorporates the cyber-security tool developed by the team into the interface. It consists of the regular display in current control rooms, as well as a panel used to indicate the occurrence of cyber-attacks or safety events and a panel used to guide the operator’s response to detected cyber-attacks.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Sterile neutrinos and the global reactor antineutrino dataset

We present results from global fits to the available reactor antineutrino dataset, as of Fall 2019, to determine the global preference for a fourth, sterile neutrino. We have separately considered experiments that measure the integrated inverse-beta decay (IBD) rate from those that measure the energy spectrum of IBD events at one or more locations. The evidence that we infer from rate measurements varies between ≲ 3σ and negligible depending on the reactor antineutrino flux model employed. Moreover, we find that spectral ratios ostensibly imply ≳ 3σ evidence, consistent with previous work, though these measurements are known to be plagued by issues related to statistical interpretation; these results should therefore be viewed cautiously. The software used is the newly developed GLoBESfit tool set which is based on the publicly available GLoBES framework and will be released as open-source software.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Bayesian homodyne and heterodyne tomography

Continuous-variable (CV) photonic states are of increasing interest in quantum information science, bolstered by features such as deterministic resource state generation and error correction via bosonic codes. Data-efficient characterization methods will prove critical in the fine-tuning and maturation of such CV quantum technology. Although Bayesian inference offers appealing properties—including uncertainty quantification and optimality in mean-squared error—Bayesian methods have yet to be demonstrated for the tomography of arbitrary CV states. Here we introduce a complete Bayesian quantum state tomography workflow capable of inferring generic CV states measured by homodyne or heterodyne detection, with no assumption of Gaussianity. As examples, we demonstrate our approach on experimental coherent, thermal, and cat state data, obtaining excellent agreement between our Bayesian estimates and theoretical predictions. Our approach lays the groundwork for Bayesian estimation of highly complex CV quantum states in emerging quantum photonic platforms, such as quantum communications networks and sensors.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A magnon scattering platform

Significance This work describes a general scattering platform that uses magnons to explore the underlying properties of target materials. In this work we show how both phase and amplitude of magnons can be imaged using a nitrogen vacancy center magnetometer and how the scattered pattern of waves can be used to infer geometric and magnetic properties of a target material. To demonstrate this new experimental methodology we use a permalloy disk as our target and show that even with such a simple target unexpected behavior is observed. In addition, we provide a theoretical framework to reconstruct the properties of the target.

36 MATERIALS SCIENCE↗

Enhanced MPM framework with multipatch isogeometric analysis for geotechnical applications

Achieving stable stress solutions at large strains using the Material Point Method (MPM) is challenging due to the accumulation of errors associated with geometry discretization, cell-crossing noise, and volumetric locking. Several simplified attempts exist in the literature to mitigate these errors, including higher-order frameworks. However, the stability of the MPM solution in such frameworks has been limited to simple geometries and the single-phase formulation (i.e., neglecting pore fluid). Although never explored, multipatch isogeometric analysis offers desirable qualities to simulate complex geometries while mitigating errors in the MPM. The degree of required high-order spatial integration has also never been investigated to infer a minimum limit for the stability of the stress solution in MPM. This paper presents a general-purpose numerical framework for simulating stable stresses in porous media, capturing both near incompressibility and multiphase interactions. First, the numerical framework is presented considering Non-Uniform Rational B-splines (NURBS) to perform isogeometric analysis (IGA) in MPM. Additionally, a volumetric strain smoothing algorithm is used to alleviate errors associated with volumetric locking. Second, the manifestation of cell-crossing errors is assessed via a series of problems with orders ranging from linear to cubic interpolation functions. Third, the use of NURBS is investigated and verified for problems with circular geometries. Finally, multipatch analysis is deployed to simulate plane strain and 3D penetration in soils, considering nearly incompressible elastoplastic (total stress) analysis and fully-coupled hydro-mechanical (effective stress) analysis. The stability of the solution is also analyzed for different constitutive models. From the results, it can be concluded that the framework using cubic interpolation functions with strain smoothing is the most convenient, presenting stable stress solutions for a broad range of multiphase geotechnical applications.

58 GEOSCIENCES↗

Advancements in Constitutive Model Calibration: Leveraging the Power of Full‐Field DIC Measurements and In Situ Load Path Selection for Reliable Parameter Inference

Accurate material characterization and model calibration are essential for computationally supported high-consequence engineering decisions. Historically, characterization and calibration methods (1) use simplified test specimen geometries and global data, (2) cannot guarantee that sufficient characterization data are collected for a specific model of interest, (3) use deterministic methods that provide best-fit parameter values with no uncertainty quantification, and (4) are sequential, inflexible, and time-consuming. This work brings together several recent advancements into an improved workflow called interlaced characterization and calibration (ICC) that advances the state-of-the-art in constitutive model calibration. The ICC paradigm (1) employs tools to efficiently use full-field data to calibrate high-fidelity material models, (2) aligns the data needed with the data collected by adopting an optimal experimental design protocol, (3) quantifies parameter uncertainty through Bayesian inference and (4) incorporates these advancements into a quasi real-time feedback loop. The ICC framework is demonstrated here on the calibration of a material model using simulated full-field data for an aluminium cruciform specimen being deformed biaxially. The cruciform is actively driven through the myopically preferred load path using Bayesian optimal experimental design, which selects load steps that yield the maximum expected information gain (EIG). Principal component analysis (PCA) is performed on the model predictions of full-field displacements, and fast surrogate models are built to approximate the input-output relationships of the expensive finite element model. Furthermore, the tools developed and demonstrated here show that high-fidelity constitutive models can be efficiently and reliably calibrated with quantified uncertainty, thus supporting credible decision-making and potentially increasing the agility of solid mechanics modelling by enabling utilization of computational simulations at earlier stages of the design cycle.

Bayesian optimal experimental design↗

SZ–X-Ray Surface Brightness Fluctuations in the SPT-XMM Clusters

The hot plasma in galaxy clusters, the intracluster medium, is expected to be shaped by subsonic turbulent motions, which are key for heating, cooling, and transport mechanisms. The turbulent motions contribute to the nonthermal pressure, which, if not accounted for, consequently imparts a hydrostatic mass bias. Accessing information about turbulent motions is thus of major astrophysical and cosmological interest. Characteristics of turbulent motions can be indirectly accessed through surface brightness fluctuations. This study expands on our pilot investigations of surface brightness fluctuations in the Sunyaev–Zel’dovich and in X-ray data by examining, for the first time, a large sample of 60 clusters using both SPT-SZ and XMM-Newton data and spans the redshift range 0.2 < z < 1.5, thus constraining the respective pressure and density fluctuations within 0.6 R $_{500}$. We deem density fluctuations to be of sufficient quality for 32 clusters, finding mild correlations between the peak of the amplitude spectra of density fluctuations and various dynamical parameters. We infer turbulent velocities from density fluctuations with an average Mach number $\mathcal{M}$ 3D = 0.52 ± 0.14, in agreement with numerical simulations. For clusters with inferred turbulent Mach numbers from fluctuations in both pressure, $\mathcal{M}$ P , and density, $\mathcal{M}$ ρ , we find broad agreement between $\mathcal{M}$ P and $\mathcal{M}$ ρ . Our results suggest either a bimodal or a skewed unimodal Mach number distribution, with the majority of clusters being turbulence-dominated (subsonic) while the remainder are shock-dominated (supersonic).

79 ASTRONOMY AND ASTROPHYSICS↗

Cannabis monoterpene synthases: evaluating structure–function relationships

Terpene synthases catalyze the first committed step in the biosynthesis of terpenes, a structurally diverse class of natural products that also encompasses volatiles derived from precursors in the C10 to C15 range (termed monoterpenes and sesquiterpenes, respectively). In the review section of this article, we are providing information about all functionally characterized monoterpene synthases (MTSs) and sesquiterpene synthases (STSs) of Cannabis sativa L. We are also exploring the locations of MTSs and STSs in the chromosome-level assembly of the reference chemovar CBDRx. A follow-up computational structure–function analysis focuses on MTSs, as there is already a rich literature available on the topic. More specifically, by employing sequence comparisons and homology structural modeling, we infer which amino acid residues are likely to constrain the available space in the active site of cannabis MTSs. The emphasis of these studies was to investigate why some MTSs accept only a C10 diphosphate as substrate, while mixed MTS/STS enzymes also accommodate a C15 diphosphate. Here, by combining a literature review and computational analyses in a hybrid format, we are laying the foundation for future studies to better understand the determinants of substrate and product specificity in these fascinating enzymes.

59 BASIC BIOLOGICAL SCIENCES↗

A numerical study of the meso-structure variability in the compaction process of prepreg platelet molded composites

Compaction deformations in a thermoplastic composite system molded from chopped prepreg platelets were analyzed to obtain statistical distributions of local distortions that define the final consolidated meso-structure. This study provides an insight into and allows quantification of the intrinsic geometric variability of the composite meso-structure that controls the variability of the effective composite mechanical properties. A two-phase viscoelastic finite element model was utilized to treat the behavior of molten prepreg platelets, which were assumed to be incompressible and inextensible in the fiber direction during processing. A two-step simulation procedure with discretely modeled platelets and platelet-to-platelet contacts is proposed. The analysis begins with a drop simulation of the uncontrolled material deposition followed by a 3D finite element analysis of the preform compaction and bulk material flow. The proposed model provides a method for generating a composite meso-structure ab initio starting from a statistical sampling of the initial conditions for platelet deposition. The predicted platelet thickness distribution, intra-platelet fiber waviness and out-of-plane platelet undulations were found to be in good agreement with experimental measurements. An estimate of the volume of resin pockets was inferred from the compaction simulation.

36 MATERIALS SCIENCE↗

Runaway electron seed formation at reactor-relevant temperature

Systematic variation of the pre-disruption core electron temperature (T e ) from 1 to 12 keV using an internal transport barrier scenario reveals a dramatic increase in the production of ‘seed’ runaway electrons (REs), ultimately accessing near-complete conversion of the pre-disruption current into sub-MeV RE current. Injected Ar pellets are observed to ablate more intensely and promptly as T e rises. At high T e , the observed ablation exceeds predictions from published thermal ablation models. Simultaneously, the thermal quench (TQ) is observed to significantly shorten with increasing T e —a surprising result. While the reason for the shorter TQ is not yet understood, candidate mechanisms include: insufficiently accurate thermal ablation models, enhanced ablation driven by the seed RE population, or significant parallel heat transport along stochastic fields. Kinetic modeling that self-consistently treats the plasma cooling via radiation, the induced electric field, and the formation of the seed RE is performed. Including the combined effect of the inherent dependence of hot-tail RE seeding on T e together with the shortened TQ, modeling recovers the progression towards near-complete conversion of the pre-disruption current to RE current as T e rises. Measurement of the HXR spectrum during the early current quench (CQ) reveals a trend of decreasing energy with pre-disruption T e . At the very highest T e (≈ 12 keV), ≈ 100% conversion of the thermal current to runaway current is found. The energy of this peculiar RE beam is inferred to be sub-MeV as it emits vanishingly few MeV hard x-rays (HXRs). These measurements demonstrate novel TQ dynamics as Te is varied and illustrate the limitations of treating the RE seed formation problem without considering the inter-related dependencies of the pellet ablation, radiative energy loss, and resultant variations of the TQ duration. Finally, if the observed shortening of the TQ with increasing T e extends to fusion-grade plasmas, than their propensity to form large quantities of RE seeds at high T e may be far worse than previously thought. Positively, the high T e scenario in DIII-D produces REs so prodigiously that it can serve as a meaningful new platform for demonstrating RE avoidance techniques.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A Magnetized Strongly Turbulent Corona as the Source of Neutrinos from NGC 1068

Abstract The cores of active galactic nuclei are potential accelerators of 10–100 TeV cosmic rays, in turn producing high-energy neutrinos. This picture was confirmed by the compelling evidence of a TeV neutrino signal from the nearby active galaxy NGC 1068, leaving open the question of what is the site and mechanism of cosmic-ray acceleration. One candidate is the magnetized turbulence surrounding the central supermassive black hole. Recent particle-in-cell simulations of magnetized turbulence indicate that stochastic cosmic-ray acceleration is nonresonant, in contrast to the assumptions of previous studies. We show that this has important consequences on a self-consistent theory of neutrino production in the corona, leading to a more rapid cosmic-ray acceleration than previously considered. The turbulent magnetic-field fluctuations needed to explain the neutrino signal are consistent with a magnetically powered corona. We find that strong turbulence, with turbulent magnetic energy density higher than 1% of the rest-mass energy density, naturally explains the normalization of the IceCube neutrino flux, in addition to the neutrino spectral shape. Only a fraction of the protons in the corona, which can be directly inferred from the neutrino signal, are accelerated to high energies. Thus, in this framework, the neutrino signal from NGC 1068 provides a testbed for particle acceleration in magnetized turbulence.

Fiorillo, Damiano F. G. (ORCID:0000000349279850)↗

MOCVD Surface Preparation of V-Groove Si for III-V Growth

V-groove nanopatterning of Si substrates has recently demonstrated promise for achieving high-quality III-V-on-Si epitaxy while providing a lower-cost processing route than chemo-mechanical polishing to produce epi-ready planar wafers. A key factor in determining the crystalline quality of III-V buffer layers is the Si surface structure and its chemical composition. Unlike planar Si surfaces, the surfaces of V-grooves prior to growth have not been studied in detail. Here, we study the surface of V-groove Si prepared for GaP nucleation via X-ray photoelectron spectroscopy and low-energy electron diffraction. We identify several pretreatments, using both 830 and 1000 annealing under an As background pressure, as being suitable for deoxidizing and cleaning the V-groove Si surface. The V-groove Si was found to behave similarly to reference Si(0 0 1) and Si(1 1 1) planar samples, demonstrating that in situ techniques such as reflection anisotropy spectroscopy can be used on reference samples to infer the state of the V-groove surface, and indicating that the extensive research on planar Si surfaces can be directly applied to V-grooves.

III-V epitaxy↗

Bond order predictions using deep neural networks

Machine learning is an extremely powerful tool for the modern theoretical chemist since it provides a method for bypassing costly algorithms for solving the Schrödinger equation. Already, it has proven able to infer molecular and atomic properties such as charges, enthalpies, dipoles, excited state energies, and others. Most of these machine learning algorithms proceed by inferring properties of individual atoms, even breaking down total molecular energy into individual atomic contributions. In this paper, we introduce a modified version of the Hierarchically Interacting Particle Neural Network (HIP-NN) capable of making predictions on the bonds between atoms rather than on the atoms themselves. We train the modified HIP-NN to infer bond orders for a large number of small organic molecules as computed via the Natural Bond Orbital package. We demonstrate that the trained model is extensible to molecules much larger than those in the training set by studying its performance on the COMP6 dataset. This method has applications in cheminformatics and force field parameterization and opens a promising future for machine learning models to predict other quantities that are defined between atoms such as density matrix elements, Hamiltonian parameters, and molecular reactivities.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗