Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “agent based model”

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

Multiscale modeling of HBV infection integrating intra- and intercellular viral propagation to analyze extracellular viral markers

Chronic infection with hepatitis B virus (HBV) is caused by the persistence of closed circular DNA (cccDNA) in the nucleus of infected hepatocytes. Despite available therapeutic anti-HBV agents, eliminating the cccDNA remains challenging. Thus, quantifying and understanding the dynamics of cccDNA are essential for developing effective treatment strategies and new drugs. However, such study requires repeated liver biopsy to measure the intrahepatic cccDNA, which is basically not accepted because liver biopsy is potentially morbid and not common during hepatitis B treatment. We here aimed to develop a noninvasive method for quantifying cccDNA in the liver using surrogate markers in peripheral blood. We constructed a multiscale mathematical model that explicitly incorporates both intracellular and intercellular HBV infection processes. The model, based on age-structured partial differential equations, integrates experimental data from in vitro and in vivo investigations. By applying this model, we roughly predicted the amount and dynamics of intrahepatic cccDNA within a certain range using specific viral markers in serum samples, including HBV DNA, HBsAg, HBeAg, and HBcrAg. Our study represents a significant step towards advancing the understanding of chronic HBV infection. The noninvasive quantification of cccDNA using our proposed method holds promise for improving clinical analyses and treatment strategies. By comprehensively describing the interactions of all components involved in HBV infection, our multiscale mathematical model provides a valuable framework for further research and the development of targeted interventions.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Trust in News as Moderated by Known Truth Assessments and Source

Disinformation has become a problem in various spaces, from social media to organizations. Given the lack of evidence and transparency surrounding organizational disinformation, we aim to understand issue of information exposure effect and source trust from the lens of traditional information consumption online. Using an online experiment and multi-level mixed-effects modeling, we find that the exposure effect exists even over only two exposures to a headline, as long as participants are not shown fact-checking. Fact-checking can be effective in reducing the impact of the exposure effect. Additionally, we find that participants trusted sources significantly more before seeing the headlines. The possible presence of disinformation reduces trust in sources, even when considered highly neutral and reputable. These findings can inform simulations of agent interactions with misinformation as well as future experimental designs focused on evaluations of belief and trust based on recognizability of a source.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2020

Estimation of Large-Scale Wind Field Characteristics Using Supervisory Control and Data Acquisition Measurements

As the wind energy industry continues to push for increased power production and lower cost of energy, the focus of research has expanded from individual turbines to entire wind farms. Among a host of interesting problems to be solved when considering the wind farm as a whole, we consider the challenge of scalar field estimation, based on information already collected at the individual turbine level. We aim to estimate the large-scale, low-frequency characteristics of the wind field, such as the mean wind direction and the overall decrease in wind speed across the farm, and employ a Kalman filter that models the wind field using a polynomial function. We compare the proposed method’s performance to both a simple averaging technique and filtering of individual turbine measurements. The method presented is not limited to wind turbines and is applicable in other situations where multiple remote agents are used to estimate a scalar field.

estimation↗
DOE OSTI2020

3C-like protease inhibitors block coronavirus replication in vitro and improve survival in MERS-CoV–infected mice

Pathogenic coronaviruses are a major threat to global public health, as exemplified by severe acute respiratory syndrome coronavirus (SARS-CoV), Middle East respiratory syndrome coronavirus (MERS-CoV), and the newly emerged SARS-CoV-2, the causative agent of coronavirus disease 2019 (COVID-19). We describe herein the structure-guided optimization of a series of inhibitors of the coronavirus 3C-like protease (3CLpro), an enzyme essential for viral replication. The optimized compounds were effective against several human coronaviruses including MERS-CoV, SARS-CoV, and SARS-CoV-2 in an enzyme assay and in cell-based assays using Huh-7 and Vero E6 cell lines. Two selected compounds showed antiviral effects against SARS-CoV-2 in cultured primary human airway epithelial cells. In a mouse model of MERS-CoV infection, administration of a lead compound 1 day after virus infection increased survival from 0 to 100% and reduced lung viral titers and lung histopathology. These results suggest that this series of compounds has the potential to be developed further as antiviral drugs against human coronaviruses.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2023

Reinforcement Learning of Structured Stabilizing Control for Linear Systems With Unknown State Matrix

This paper delves into designing feedback control gains for a continuous-time linear quadratic regulator (LQR) problem that is constrained to certain predefined structure with unknown state matrix. We bring forth the ideas from reinforcement learning (RL) in conjunction with sufficient stability and performance guarantees in order to design these structured gains using the trajectory measurements of states and controls. Here we first formulate a model-based framework using dynamic programming (DP) to embed the structural constraint to the LQR gain computation in the continuous-time setting, and then subsequently, formulate a policy iteration RL algorithm that can alleviate the requirement of known state matrix in conjunction with maintaining the feedback gain structure. The design enables a distributed learning control design which is necessary for many large-scale cyber-physical systems. Theoretical guarantees are provided for stability and convergence of the structured reinforcement learning (SRL) algorithm. We validate our theoretical results with numerical simulations on a multi-agent networked linear time-invariant (LTI) dynamic system.

42 ENGINEERING↗
DOE OSTI2024

Investigation of the Potential of Repurposing Medium-Density Fiberboard Waste as an Adsorbent for Heavy Metal Ion Removal

Medium-density fiberboard (MDF) waste generation has increased steadily over the past decades, and therefore, the investigation of novel methods to recycle this waste is very important. The potential of repurposing MDF waste as an adsorbent for the treatment of Cd(II), Cu(II), Pb(II), and Zn(II) ions in water was investigated using MDF offcuts. The highest adsorption potential in single-metal ion solution systems was observed for Pb(II) ions. The experimental data of Pb(II) ions fit well with the Freundlich isotherm and pseudo-second-order kinetic models. Complexation and electrostatic interactions were identified as the adsorption mechanisms. The adsorption behavior of multi-metal ion adsorption systems was investigated by introducing Cd(II) ions as a competitive metal ion. The presence of the Cd(II) ions reduced the adsorption potential of Pb(II) ions, yet the preference for the Pb(II) ions remained. Regeneration studies were performed by using 0.1 M HCl as a regeneration agent for both systems. Even though a significant amount of adsorbed metal ions were recovered, the adsorption potential of the MDF was reduced in the subsequent adsorption cycles. Based on these results, MDF fines have the potential to be used as an economical adsorbent for remediation of wastewater containing heavy metal ions.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Reinforcement Learning for feedback-enabled cyber resilience

The rapid growth in the number of devices and their connectivity has enlarged the attack surface and made cyber systems more vulnerable. As attackers become increasingly sophisticated and resourceful, mere reliance on traditional cyber protection, such as intrusion detection, firewalls, and encryption, is insufficient to secure the cyber systems. Cyber resilience provides a new security paradigm that complements inadequate protection with resilience mechanisms. A Cyber-Resilient Mechanism (CRM) adapts to the known or zero-day threats and uncertainties in real-time and strategically responds to them to maintain the critical functions of the cyber systems in the event of successful attacks. Feedback architectures play a pivotal role in enabling the online sensing, reasoning, and actuation process of the CRM. Reinforcement Learning (RL) is an important gathering of algorithms that epitomize the feedback architectures for cyber resilience. It allows the CRM to provide dynamic and sequential responses to attacks with limited or without prior knowledge of the environment and the attacker. In this work, we review the literature on RL for cyber resilience and discuss the cyber-resilient defenses against three major types of vulnerabilities, i.e., posture-related, information-related, and human-related vulnerabilities. Here we introduce moving target defense, defensive cyber deception, and assistive human security technologies as three application domains of CRMs to elaborate on their designs. The RL algorithms also have vulnerabilities themselves. We explain the major vulnerabilities of RL and present develop several attack models where the attacker target the information exchanged between the environment and the agent: the rewards, the state observations, and the action commands. We show that the attacker can trick the RL agent into learning a nefarious policy with minimum attacking effort. The paper introduces several defense methods to secure the RL-enabled systems from these attacks. However, there is still a lack of works that focuses on the defensive mechanisms for RL-enabled systems. Last but not least, we discuss the future challenges of RL for cyber security and resilience and emerging applications of RL-based CRMs.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Development and Validation of a Simulation Testbed for the Intelligent Building Agents Laboratory (IBAL) Using TRNSYS

This paper documents the development and validation of a dynamic primary cooling and thermal storage system simulation testbed. The system simulation testbed, sIBAL, is based on the Intelligent Building Agents Laboratory (IBAL) at the National Institute of Standards and Technology (NIST), which is a research infrastructure and testbed for the development, evaluation, and demonstration of intelligent control algorithms. The sIBAL testbed developed in this project will serve as a virtual twin of the real facility for future control algorithm development. The details of the methodologies used to develop and validate the simulation testbed, which replicates the dynamic behaviors of the primary cooling and ice storage system in the IBAL facility, are presented. The simulation testbed was developed in TRNSYS using built-in component models and MATLAB functions to replicate the two water-cooled chillers, a thermal storage tank, pumps, valves and other components for four different operation modes. Experiments on IBAL components were designed and executed to generate experimental data for model development and verification of the simulation platform. The validation of the simulation results was carried out in two phases: 1) independent component simulations for the chillers and thermal storage tank, and 2) a combined testbed simulation of the entire hydronic system. Comparison of simulation results to the experimental data obtained from the IBAL facility showed errors within 1 °C for the temperatures outputs of both the chiller and the thermal storage model. The error is within an acceptable range for further intelligent control algorithms development. The findings from the study are summarized and presented along with areas where additional research is needed. In addition, data filtering procedures and model refinement measures utilized to improve the accuracy and accelerate the computation time of the simulation are presented.

cooling systems↗
DOE OSTI2021

Learning Sequential Distribution System Restoration via Graph-Reinforcement Learning

We report a distribution service restoration algorithm as a fundamental resilient paradigm for system operators provides an optimally coordinated, resilient solution to enhance the restoration performance. The restoration problem is formulated to coordinate distribution generators and controllable switches optimally. A model-based control scheme is usually designed to solve this problem, relying on a precise model and resulting in low scalability. To tackle these limitations, this work proposes a graph-reinforcement learning framework for the restoration problem. We link the power system topology with a graph convolutional network, which captures the complex mechanism of network restoration in power networks and understands the mutual interactions among controllable devices. Latent features over graphical power networks produced by graph convolutional layers are exploited to learn the control policy for network restoration using deep reinforcement learning. The solution scalability is guaranteed by modeling distributed generators as agents in a multi-agent environment and a proper pre-training paradigm. Comparative studies on IEEE 123-node and 8500-node test systems demonstrate the performance of the proposed solution.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2022

Diffuse Radiation Forcing Constraints on Gross Primary Productivity and Global Terrestrial Evapotranspiration

Abstract The diffuse radiation fertilization effect—the increase in plant productivity in the presence of higher diffuse radiation ( K ↓,d )—is an important yet understudied aspect of atmosphere‐biosphere interactions and can modify the terrestrial carbon, energy, and water budgets. The K ↓,d fertilization effect links the carbon cycle with clouds and aerosols, all of which are large sources of uncertainties for our current understanding of the Earth system and for future climate projections. Here we establish to what extent observational and modeling uncertainty in sunlight's diffuse fraction ( k d ) affects simulated gross primary productivity (GPP) and terrestrial evapotranspiration ( λE ). We find only 48 eddy covariance sites with simultaneous sufficient measurements of K ↓,d with none in the tropical climate zone, making it difficult to constrain this mechanism globally using observations. Using a land modeling framework based on the latest version of the Community Land Model, we find that global GPP ranges from 114 Pg C year −1 when using k d forcing from the Modern‐Era Retrospective analysis for Research and Applications, version 2 reanalysis to a ∼7% higher value of 122 Pg C year −1 when using the Clouds and the Earth's Radiant Energy System satellite product, with especially strong differences apparent over the tropical region (mean increase ∼9%). The differences in λE , although smaller (−0.4%) due to competing changes in shaded and sunlit leaf transpiration, can be greater than regional impacts of individual forcing agents like aerosols. Our results demonstrate the importance of comprehensively and systematically validating the simulated k d by atmosphere modules as well as the response to differences in k d within land modules across Earth System Models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Innovating the next generation of commercial smart building software

Nearly 30% of commercial building energy use is wasted due to equipment faults and HVAC controls problems. The result is increased emissions, compromised comfort and productivity, and less reliable coordination of building power needs with a clean grid. The energy impact alone represents $17 billion in potential savings. Today’s smart building software provides a robust solution to address these operational deficiencies. Energy management and information systems (EMIS) are saving up to 9% on average, with two-year paybacks. They are being incorporated into energy management processes, commissioning services, and utility programs. As effective as they are, two barriers prevent even deeper benefits; limited personnel to fix problems once they are identified, and the expense and time to manually implement changes in control systems. In partnership with the research community, the EMIS industry is developing new capabilities to overcome these barriers. Moving beyond siloed products for either fault detection and diagnostics, or optimal control, these new capabilities empower users to not only automatically identify faults, but also to push corrective action, and control improvements to their buildings. In this paper, several areas for enhancements are documented: ‘one-time’ correction of faults such as setpoints, schedules, and economizer lockouts; short-term active testing for automated proportional integral derivative (PID) loop tuning and functional testing; and continuous supervisory control for demand flexibility and year-round efficiency. Results are presented from a pair of partner implementations out of a dozen providers integrating these enhancements into their products, including field tests from across the country, and insights into operator acceptance and integration into operations and maintenance practices.

Casillas, Armando↗
DOE OSTI2025

Deep reinforcement learning control for co-optimizing energy consumption, thermal comfort, and indoor air quality in an office building

With the recent demand for decarbonization and energy efficiency, advanced HVAC control using Deep Reinforcement Learning (DRL) becomes a promising solution. Due to its flexible structures, DRL has been successful in energy reduction for many HVAC systems. However, only a few researches applied DRL agents to manage the entire central HVAC system and control multiple components in both the water loop and the air loop, owing to its complex system structures. Moreover, those researches have not extended their applications by incorporating the indoor air quality, especially both CO2 and PM2.5concentrations, on top of energy saving and thermal comfort, as achieving those objectives simultaneously can cause multiple control conflicts. What's more, DRL agents are usually trained on the simulation environment before deployment, so another challenge is to develop an accurate but relatively simple simulator. Therefore, we propose a DRL algorithm for a central HVAC system to co-optimize energy consumption, thermal comfort, indoor CO2 level, and indoor PM2.5 level in an office building. To train the controller, we also developed a hybrid simulator that decoupled the complex system into multiple simulation models, which are calibrated separately using laboratory test data. The hybrid simulator combined the dynamics of the HVAC system, the building envelope, as well as moisture, CO2, and particulate matter transfer. Three control algorithms (rule-based, MPC, and DRL) are developed, and their performances are evaluated on the hybrid simulator environment with a realistic scenario (i.e., with stochastic noises). The test results showed that, the DRL controller can save 21.4 % of energy compared to a rule-based controller, and has improved thermal comfort, reduced indoor CO2 concentration. The MPC controller showed an 18.6 % energy saving compared to the DRL controller, mainly due to savings from comfort and indoor air quality boundary violations caused by unmeasured disturbances, and it also highlights computational challenges in real-time control due to non-linear optimization. Finally, we provide the practical considerations for designing and implementing the DRL and MPC controllers based on their respective pros and cons.

Guo, Fangzhou↗
DOE OSTI2020

Forced Patterns of Sea Level Rise in the Community Earth System Model Large Ensemble From 1920 to 2100

To provide context for observed sea level rise, the forced response (FR) in dynamic sea level (DSL) during the 20th and 21st centuries is examined in the Community Earth System Model Large Ensemble (LE). After accounting for simulation drift, which in the LE is sizable, the DSL FR is found to be complex, both in space and time. Its evolving character is suggested to arise from both the diversity and transient evolution of climate forcing agents and the slow adjustment timescales of the intermediate and deep oceans. Nonetheless, various intervals of spatially coherent change simulated for the recent past and near future are identified, and their characteristics and associated driving mechanisms are identified and discussed. The mid- to late-20th century DSL FR is characterized by an hemispherically asymmetric pattern of change, with depressed rates of rise in the northern oceans. There is also a dipole of change in the Southern Ocean caused by changes in near-surface zonal winds. Through the late 20th and early 21st centuries, a different pattern of rise emerges with elevated rates in the tropics and depressed rates at high latitudes. Zonal and interbasin variations characterize both intervals and involve the pattern and depth of anomalous ocean heat content storage and spatial contrasts in the expansion coefficient tied mainly to base state temperature, with greater rates of rise in warm regions per unit warming. Finally, the relative roles of surface flux and ocean convergence anomalies are examined.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Chronic Cardiotoxicity Assays Using Human Induced Pluripotent Stem Cell-Derived Cardiomyocytes (hiPSC-CMs)

Cardiomyocytes (CMs) differentiated from human induced pluripotent stem cells (hiPSCs) are increasingly used in cardiac safety assessment, disease modeling and regenerative medicine. A vast majority of cardiotoxicity studies in the past have tested acute effects of compounds and drugs; however, these studies lack information on the morphological or physiological responses that may occur after prolonged exposure to a cardiotoxic compound. In this review, we focus on recent advances in chronic cardiotoxicity assays using hiPSC-CMs. We summarize recently published literature on hiPSC-CMs assays applied to chronic cardiotoxicity induced by anticancer agents, as well as non-cancer classes of drugs, including antibiotics, anti-hepatitis C virus (HCV) and antidiabetic drugs. We then review publications on the implementation of hiPSC-CMs-based assays to investigate the effects of non-pharmaceutical cardiotoxicants, such as environmental chemicals or chronic alcohol consumption. We also highlight studies demonstrating the chronic effects of smoking and implementation of hiPSC-CMs to perform genomic screens and metabolomics-based biomarker assay development. The acceptance and wide implementation of hiPSC-CMs-based assays for chronic cardiotoxicity assessment will require multi-site standardization of assay protocols, chronic cardiac maturity marker reproducibility, time points optimization, minimal cellular variation (commercial vs. lab reprogrammed), stringent and matched controls and close clinical setting resemblance. A comprehensive investigation of long-term repeated exposure-induced effects on both the structure and function of cardiomyocytes can provide mechanistic insights and recapitulate drug and environmental cardiotoxicity.

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

Enzalutamide response in a panel of prostate cancer cell lines reveals a role for glucocorticoid receptor in enzalutamide resistant disease

Representative in vitro model systems that accurately model response to therapy and allow the identification of new targets are important for improving our treatment of prostate cancer. Here we describe molecular characterization and drug testing in a panel of 20 prostate cancer cell lines. The cell lines cluster into distinct subsets based on RNA expression, which is largely driven by functional Androgen Receptor (AR) expression. KLK3 , the AR-responsive gene that encodes prostate specific antigen, shows the greatest variability in expression across the cell line panel. Other common prostate cancer associated genes such as TMPRSS2 and ERG show similar expression patterns. Copy number analysis demonstrates that many of the most commonly gained (including regions containing TERC and MYC ) and lost regions (including regions containing TP53 and PTEN) that were identified in patient samples by the TCGA are mirrored in the prostate cancer cell lines. Assessment of response to the anti-androgen enzalutamide shows a distinct separation of responders and non-responders, predominantly related to status of wild-type AR. Surprisingly, several AR-null lines responded to enzalutamide. These AR-null, enzalutamide-responsive cells were characterized by high levels of expression of glucocorticoid receptor (GR) encoded by NR3C1 . Treatment of these cells with the anti-GR agent mifepristone showed that they were more sensitive to this drug than enzalutamide, as were several of the enzalutamide non-responsive lines. This is consistent with several recent reports that suggest that GR expression is an alternative signaling mechanism that can bypass AR blockade. This study reinforces the utility of large cell line panels for the study of cancer and identifies several cell lines that represent ideal models to study AR-null cells that have upregulated GR to sustain growth.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Quantifications of growth rate and particle size of magnesium silicate hydrate cement using small‐angle X‐ray scattering

The focus of recent research on cement binders has shifted to alternatives such as magnesium silicate hydrate (M-S-H) to mitigate the environmental impact of conventional concrete production. However, a detailed understanding of the morphology and growth rate of M-S-H, which informs its applicability as a cementing agent, is lacking. These properties are important as they influence cement hydration and consequently the development of macroscopic properties. Here, in this study, the particle size distribution and growth rate of M-S-H synthesized by mixing aqueous stock solutions were analyzed using synchrotron-based small-angle X-ray scattering. The reaction solutions contained dissolved magnesium and silicon with concentrations ranging from 2 to 500 mM, corresponding to saturation indices with respect to M-S-H endmembers, M 0.75 -S-H and M 1.50 -S-H, between 5.8 and 10.3. The precipitation rates of M-S-H align with affinity-based kinetic models and can be described by the following equations: For M 0.75 -S-H: 𝑟 M 0.75 −S−H =10 −9.54±0.70 ⁢(Ω$^{0.59±0.08}_{M_{0.75}⁢\textrm{SH}}$−1), and for M 1.50 -S-H: 𝑟 M 1.50 −S−H =10 −7.56±0.37 ⁢(Ω$^{0.34±0.04}_{M_{1.50}⁢\textrm{SH}}$−1), where 𝑟 is the rate in mol/L/s, and Ω represents the saturation ratio with respect to M-S-H. Moreover, we observed that changing the [Mg]/[Si] ratio affects particle size, a trend that cannot be fully explained by changes in the saturation index alone. These insights are critical for improving the predictability and control of cementation in concrete that utilizes M-S-H-based binders.

M-S-H↗
DOE OSTI2025

Hydrogen and Electric Charging Infrastructure for Heavy-Duty Trucks: A Nationally Scalable Megaregion Assessment

Decarbonizing regional and long-haul freight is challenging due to the limitations of battery-electric commercial vehicles and infrastructure constraints. Hydrogen fuel cell medium- and heavy-duty vehicles (MHDVs) offer a viable alternative, aligning with the decarbonization goals of the Department of Energy and commercial entities. Historically, alternative fuels like compressed natural gas and liquefied propane gas have faced slow adoption due to barriers like infrastructure availability. To avoid similar issues, effective planning and deploying zero-emission hydrogen fueling infrastructure is crucial. This research develops deployment plans for affordable, accessible, and sustainable hydrogen refueling stations, supporting stakeholders in the decarbonized commercial vehicle freight system. It aims to benefit underserved and rural energy-stressed communities by improving air quality, reducing noise pollution, and enhancing energy resiliency. This research also provides a blueprint for replacing diesel in over-the-road Class 8 freight truck applications with hydrogen fueling solutions. The study focuses on the Texas Triangle Megaregion (I-45, I-35, and I-10), the I-10 corridor between San Antonio, TX, and Los Angeles, CA, and the I-5/CA-99 corridors between Los Angeles, CA, and San Francisco, CA. This area represents a significant portion of U.S. heavy-duty freight movement, carrying ~8.5% of the national freight volume. Using the OR-AGENT (Optimal Regional Architecture Generation for Efficient National Transport) modeling framework, the study conducts an advanced assessment of commercial vehicles, road and freight networks, and energy systems. The framework integrates data on freight mobility, traffic, weather, and energy pathways to deliver a region-specific, optimized vehicles powertrain architectures, infrastructure deployment solutions, operational logistics, and energy pathways. By considering all vehicle origin-destination pairs utilizing these corridors and all feasible fueling station location options, the framework's genetic algorithm identifies the minimum number and optimal locations of hydrogen refueling stations, ensuring no vehicle is stranded. It also determines fuel schedules and quantities at each station. A roadmap for station deployment based on multiple adoption trajectories ensures a strategic rollout of hydrogen refueling infrastructure.

Sujan, Vivek [ORNL] (ORCID:0000000269882342)↗
DOE OSTI2026

Agentic Diagrammatica: Towards Autonomous Symbolic Computation in High Energy Physics

We present Diagrammatica, a symbolic computation extension to the HEPTAPOD agentic framework, which enables LLM agents to plan and execute multi-step theoretical calculations. Symbolic computation poses a distinctive reliability challenge for LLM agents, as correctness is governed by implicit mathematical conventions that are not encoded in a form that can be easily checked in the computational backend. We identify two complementary remedies, tool-constrained computation and targeted knowledge grounding, and pursue the first as the primary architecture. Concretely, we concentrate the agent's action distribution onto tool calls with convention-fixing semantics, in which the agent specifies a compact, human-auditable diagram specification and a trusted backend performs the symbolic or numerical manipulations exactly. The toolkit provides two complementary calculation paths consuming a shared diagram specification: Naive Dimensional Analysis (NDA) for order-of-magnitude rate estimates and Exact Diagrammatic Analysis (EDA) for tree-level symbolic calculations via automatic FeynCalc code generation, both supplemented by automatic Feynman diagram enumeration and a navigable theory knowledge base. The architecture is validated on two benchmarks: (1) an exhaustive catalog of all tree-level, single-vertex $1\to 2$ partial decay widths across scalar, fermion, and vector parents, with complete massless and threshold limits and Standard Model validation; and (2) an NDA sensitivity study of the muon decay multiplicity $μ^+ \to ν_μ\barν_e + n(e^+e^-) + e^-$, determining the maximum observable $n$ at current and planned muon experiments.

Menzo, Tony [Alabama U.; Fermilab] (ORCID:00000002↗