Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “computer systems organization”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,135 records · Page 63

Department of Energy’s Atmospheric System Research (ASR) Program’s Workshop on the Future of Atmospheric Large Eddy Simulation (LES): Workshop Report

Large-eddy simulation (LES) is used as a tool to understand physical processes such as turbulence, aerosols, clouds, precipitation, radiation, the interactions among all these, and their interactions with the underlying surface. Over the next 10 years, LES will drive fundamental progress in open scientific questions in these areas as LES is increasingly used to gain understanding of complex interacting physical processes involving atmospheric turbulence. This growth will be driven both by scientific demand and the expansion of computational resources needed to conduct LES, and the form that the growth takes will largely be determined by how computational resources are leveraged for scientific gain. In particular, we suggest that computational resources are likely to be leveraged in two separate but not necessarily distinct ways. On one hand, growth in computational resources will allow LES to be made more routine, that is, performed more frequently, while on the other hand, the computational expense (measured in total floating point operations) afforded to individual LES will expand dramatically, allowing simulations to increase in both domain size and resolution as well as physical detail. Current U.S. Department of Energy (DOE) projects such as LES ARM Symbiotic Simulation and Observation Activity (LASSO) are leading the way in conducting routine LES, building large, public databases that are accessible for data science, sensitivity studies, and training for machine learning. LES will also become more routine as it becomes more accessible for individual researchers to address their scientific questions of interest. Scientific questions addressed by LES over the next 10 years are likely to include cloud organization and aggregation; aerosol cloud interactions and atmospheric chemistry (including geo-engineering); urban-scale LES; atmospheric extreme events, ranging from small-scale severe weather to wildfires; and ocean-wave-atmosphere interactions. Further LES-related research will likely grow significantly in areas related to societal impact studies of air quality and extreme weather events, applications to renewable energy forecasting and resource assessment, and aid in decision-making processes. The growth in the use of LES in atmospheric science research will drive the need for better physical process representations (e.g., cloud aerosol microphysics, radiation, and atmospheric chemistry) at the scales resolved by LES. To date, many of the process representations used by LES have been taken directly from coarser-resolution models. Promising methods for LES process representations include superdroplet and quadrature methods for microphysics, 3D approaches for radiation, and better representation of chemistry and aerosol processes. At LES resolution, land-atmosphere interactions for complex terrains, land cover/types, biogeochemistry, and plant canopy models are needed as an improvement beyond traditional and widely used Monin-Obuhkov similarity theory.

54 ENVIRONMENTAL SCIENCES↗

Department of Energy’s Atmospheric System Research (ASR) Program’s Workshop on the Future of Atmospheric Large Eddy Simulation (LES) (Workshop Report)

Large-eddy simulation (LES) is used as a tool to understand physical processes such as turbulence, aerosols, clouds, precipitation, radiation, the interactions among all these, and their interactions with the underlying surface. Over the next 10 years, LES will drive fundamental progress in open scientific questions in these areas as LES is increasingly used to gain understanding of complex interacting physical processes involving atmospheric turbulence. This growth will be driven both by scientific demand and the expansion of computational resources needed to conduct LES, and the form that the growth takes will largely be determined by how computational resources are leveraged for scientific gain. In particular, we suggest that computational resources are likely to be leveraged in two separate but not necessarily distinct ways. On one hand, growth in computational resources will allow LES to be made more routine, that is, performed more frequently, while on the other hand, the computational expense (measured in total floating point operations) afforded to individual LES will expand dramatically, allowing simulations to increase in both domain size and resolution as well as physical detail. Current U.S. Department of Energy (DOE) projects such as LES ARM Symbiotic Simulation and Observation Activity (LASSO) are leading the way in conducting routine LES, building large, public databases that are accessible for data science, sensitivity studies, and training for machine learning. LES will also become more routine as it becomes more accessible for individual researchers to address their scientific questions of interest. Scientific questions addressed by LES over the next 10 years are likely to include cloud organization and aggregation; aerosol cloud interactions and atmospheric chemistry (including geo-engineering); urban-scale LES; atmospheric extreme events, ranging from small-scale severe weather to wildfires; and ocean-wave-atmosphere interactions. Further LES-related research will likely grow significantly in areas related to societal impact studies of air quality and extreme weather events, applications to renewable energy forecasting and resource assessment, and aid in decision-making processes.

54 ENVIRONMENTAL SCIENCES↗

Quantum Vibronic Effects on the Electronic Properties of Molecular Crystals

We present a study of molecular crystals, focused on the effect of nuclear quantum motion and anharmonicity on their electronic properties. We consider a system composed of relatively rigid molecules, a diamondoid crystal, and one composed of floppier molecules, NAI-DMAC, a thermally activated delayed fluorescence compound. We compute fundamental electronic gaps at the density functional theory (DFT) level of theory, with the Perdew–Burke–Erzenhof (PBE) and strongly constrained and approximately normed (SCAN) functionals, by coupling first-principles molecular dynamics with a nuclear quantum thermostat. We find a sizable zero-point renormalization (ZPR) of the band gaps, which is much larger in the case of diamondoids (0.6 eV) than for NAI-DMAC (0.22 eV). We show that the frozen phonon (FP) approximation, which neglects intermolecular anharmonic effects, leads to a large error (~50%) in the calculation of the band gap ZPR. Instead, when using a stochastic method, we obtain results in good agreement with those of our quantum simulations for the diamondoid crystal. However, the agreement is worse for NAI-DMAC where intramolecular anharmonicities contribute to the ZPR. Our results highlight the importance of accurately including nuclear and anharmonic quantum effects to predict the electronic properties of molecular crystals.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Creating multifunctional synthetic lichen platforms for sustainable biosynthesis of biofuel precursors

In this project, we were creating a sustainable platform for biofuel production, utilizing carbon-fixing autotrophs to supply oxygen and organic substrates to heterotrophic partners, which in turn produce carbon dioxide to feed the autotrophs. This symbiotic lichen community could lower the input cost, optimize metabolic exchanges and improve the generation of biofuel precursors through multi-omics driven genetic engineering. The cyanobacteria Synechococcus elongatus (S. elongatus) was used as the primary autotroph to provide oxygen and organic substrates, especially sucrose, to a co-culture system. The strain with overexpression of sucrose transporter cscB demonstrated a significant increase in sucrose production under salt stress as what we expected. We also implemented 13C metabolic flux analysis on the sucrose secreting strain S. elongatus cscB-NaCl. Next, transporters proteins like glutamate exporter mscCG from Corynebacterium glutamicum was overexpressed in S. elongatus to improve metabolite exchange. We provided sucrose and glutamate to filamentous fungi to enable their growth and production of biochemicals using substrates from the cyanobacterium. Two fungi ( Aspergillus nidulans and Aspergillus niger ) and two yeast strains ( Rhodotorula toruloides and Lipomyces starkeyi ) served as heterotrophs in this system. They were co-cultured with S.elongatus under different pH condition, and growth on different carbon source to understand the symbiotic system and optimize the parameters for production. The two fungi strained grew much better at pH 7, 8 and 9 with Yeast Nitrogen Base (YNB ) supplementation. Two yeast strains grew well on pH 7 and 8 with YNB supplementation. All four species achieved the highest biomass level when utilizing sucrose as the main carbon sources, which demonstrates a great pairing with S.elongatus . In addition to investigate the metabolite exchange and growth condition of the co-culture system, we incorporated an Aspergillus strain expressing three genes (PAND, BAPAT and HPDH). Expression of these three genes enabled production of 3-hydroxypropionic acid (3HP) with a titer of 3-6 g/L regardless of pH and other nutrient conditions. Moreover, we developed a computational model with different heterotrophic fungi symbiosing with S.elongatus to evaluate the exchange of hundreds of metabolites. The accuracy and sensitivity had been optimized by high-throughput phenotyping assays. This project demonstrated ways to enhance a synthetic lichen platform’s efficiency and scalability to create a robust system for sustainable bioproduct synthesis. The project has enabled us to explore the technological and commercial potential of sustainable lichen co-culture. The project has also trained the next generation of scientists to address sustainability, carbon fixation, and design solutions to contemporary global challenges. The lichen platform represented a novel approach to sustainably producing bioproducts with a reduced carbon footprint and lower costs.

09 BIOMASS FUELS↗

Systems-To-Atoms (S2A): enabling hydrogen for climate security

The project addresses a critical gap in hydrogen infrastructure by integrating system-level energy models with atomic-scale material simulations in a unified Systems-to-Atoms (S2A) framework. The motivation stems from the need to develop efficient, cost-effective, and durable hydrogen transport and utilization technologies to support decarbonization of hard-to-electrify sectors such as heavy-duty transportation. Current system models lack awareness of material performance mechanisms, while material-scale models do not account for system-level usage and variability. To bridge this divide, the team developed a co-simulation capability linking techno-economic analyses, reactor/process-flow modeling, and molecular-scale catalysis simulations. Applied to hydrogen delivery in California, the framework enabled comparative evaluations of compressed, cryogenic, and liquid organic hydrogen carrier (LOHC) pathways, highlighting how catalyst operation and unit process efficiency influence overall performance. The results demonstrate that no single material or transport mode is universally optimal; instead, heterogeneous solutions tuned to specific operational contexts deliver better performance. The project delivers a new capability for cross-scale material co-design, advancing hydrogen infrastructure readiness and informing DOE and LLNL missions in climate and energy resilience.

organic↗

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

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

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Improving the accuracy of the variational quantum eigensolver for molecular systems by the explicitly-correlated perturbative [2] R12 -correction

We provide an integration of the universal, perturbative explicitly correlated [2] R12 -correction in the context of the Variational Quantum Eigensolver (VQE). This approach is able to increase the accuracy of the underlying reference method significantly while requiring no additional quantum resources. The proposed approach only requires knowledge of the one- and two-particle reduced density matrices (RDMs) of the reference wavefunction; these can be measured after having reached convergence in the VQE. This computation comes at a cost that scales as the sixth power of the number of electrons. Here, we explore the performance of the VQE + [2] R12 approach using both conventional Gaussian basis sets and our recently proposed directly determined pair-natural orbitals obtained by multiresolution analysis (MRA-PNOs). Both Gaussian orbital and PNOs are investigated as a potential set of complementary basis functions in the computation of [2] R12 . In particular the combination of MRA-PNOs with [2] R12 has turned out to be very promising – persistently throughout our data, this allowed very accurate simulations at a quantum cost of a minimal basis set. Additionally, we found that the deployment of PNOs as complementary basis can greatly reduce the number of complementary basis functions that enter the computation of the correction at a complexity.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Polarized, V-Shaped, and Conjoined Biscoumarins: From Lack of Dipole Moment Alignment to High Brightness

Eleven conjoined coumarins possessing a chromeno[3,4-c]chromene-6,7-dione skeleton have been synthesized via the reaction of electron-rich phenols with esters of coumarin-3-carboxylic acids, catalyzed by either Lewis acids or 4-dimethylaminopyridine. Furthermore, Michael-type addition to angular benzo[f]coumarins is possible, leading to conjugated helical systems. Arrangement of the electron-donating amino groups at diverse positions on this heterocyclic skeleton makes it possible to obtain π-expanded coumarins with emission either sensitive to, or entirely independent of, solvent polarity with large Stokes shifts. Computational studies have provided a rationale for moderate solvatochromic effects unveiling the lack of collinearity of the dipole moments in the ground and excited states. Depending on the functional groups present, the obtained dyes are highly polarized with dipole moments of ~14 D in the ground state and ~20–25 D in the excited state. Strong emission in nonpolar solvents, in spite of the inclusion of a NO 2 group, is rationalized by the fact that the intramolecular charge transfer introduced into these molecules is strong enough to suppress intersystem crossing yet weak enough to prevent the formation of dark twisted intramolecular charge transfer states. Photochemical transformation of the dye possessing a chromeno[3,4-c]pyridine-4,5-dione scaffold led to the formation of a spirocyclic benzo[g]coumarin.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Nonadiabatic dynamics with classical trajectories: The problem of an initial coherent superposition of electronic states

Advances in coherent light sources and development of pump–probe techniques in recent decades have opened the way to study electronic motion in its natural time scale. When an ultrashort laser pulse interacts with a molecular target, a coherent superposition of electronic states is created and the triggered electron dynamics is coupled to the nuclear motion. A natural and computationally efficient choice to simulate this correlated dynamics is a trajectory-based method where the quantum-mechanical electronic evolution is coupled to a classical-like nuclear dynamics. These methods must approximate the initial correlated electron–nuclear state by associating an initial electronic wavefunction to each classical trajectory in the ensemble. Different possibilities exist that reproduce the initial populations of the exact molecular wavefunction when represented in a basis. We show that different choices yield different dynamics and explore the effect of this choice in Ehrenfest, surface hopping, and exact-factorization-based coupled-trajectory schemes in a one-dimensional two-electronic-state model system that can be solved numerically exactly. Furthermore, this work aims to clarify the problems that standard trajectory-based techniques might have when a coherent superposition of electronic states is created to initialize the dynamics, to discuss what properties and observables are affected by different choices of electronic initial conditions and to point out the importance of quantum-momentum-induced electronic transitions in coupled-trajectory schemes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Controlling Molecular Structure and Spin with Multiconfigurational Quantum Chemistry (Final Technical Report)

For many first-row transition metal complexes, structure-property relationships can be obtained from high-level molecular geometry optimizations and subsequent electronic structure studies. The project funded under this award utilized newly implemented fully internally contracted (FIC) nuclear gradients for extended multi-state (XMS) complete active space second-order multireference perturbation theory (CASPT2) to explore geometry changes in first-row transition metal coordination complexes. At the start of the project, only one full geometry optimization using FIC-CASPT2 analytical gradients had been reported (J. Chem. Theory Comput., 2016, 12 (8), 3781), and many open-questions regarding the performance and achievable accuracy in applying such computations to larger complexes persisted. Key deliverables in the report include the implementation of the numerical Hessian and subsequent vibrational analysis in the BAGEL program package. This includes both full Hessian vibrational analysis (FHVA) and a partial Hessian vibrational analysis (PHVA). The latter of which allows us to reduce the significant cost of computing the vibrational frequencies in a large molecule by focusing on the modes of highest interest. The gradient code also proved useful in the development of an approach to improve the Hubbard-U correction by removing bias in predicting spin-splitting in Fe(II) complexes via plane-wave density functional theory (DFT). Finally, we showed for three families of complexes (spin-crossover (SCO) complexes, metallocorroles complexes, and those with metal-metal bonds) that CASPT2 can result in good molecular geometries and established best practices in undertaking this work. We are now using this approach in collaborative projects with experimental groups, which would not have been possible at the start of this project. A common theme also arose through this work that static and dynamic correlation must be recovered in a balanced way to yield quantitative results. By systematically varying how both effects were included, we were able to make key chemical insights. This work supported the training of two postdoctoral scholars, one graduate student, and one undergraduate student in applying high-level wavefunction based methods to challenging systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Natural Hazard Forecast Alert Grid Risk System

Weather events cause most power outages. Often, we even get notifications on our phones to take cover or be prepared for an imminent event. If electric grid utilities had a similar warning that also included probable scenarios and the equipment involved, they could prepare and minimize the effects. Idaho National Laboratory had a project with the U.S. Department of Energy’s Cybersecurity, Energy Security, and Emergency Response program to develop a grid alert application that receives messages from the existing emergency alert system, filters and determines components possibly affected by the emergency event, calculates probable scenarios using MASTERRI (Modeling And Simulation for Targeted Reliability and Resilience Improvement). For high-risk events, the application can then send alert links to subscribed electric distribution utility operations staff to allow them to see and evaluate the scenarios and the impact in a web based interactive map tool. This proof of concept application used data from utilities and organizations, such as the international regulatory body North American Electric Reliability Corporation, which have complied historical failure data of elements that comprise the U.S. electric grid. Nominal failure rates are obtained from this data. To make this tool possible, estimated failure rates were calculated for different component types given the alert type, severity, and location. Historic weather-related grid element failures were correlated with historic weather events from the Integrated Public Alert & Warning System. These correlated events and failures are used along with Bayesian updates from the historical norms to provide a modified failure rate for grid elements in the alert areas and calculate probable scenarios. Working with an industry collaborator, actual grid models and data were used for demonstration cases. This report outlines the work performed for this project.

24 - POWER TRANSMISSION AND DISTRIBUTION↗

Stakeholder-guided holistic, Adaptive Framework for enhancing community Energy Resilience (SAFER) (Final Technical Report)

The Stakeholder-guided holistic, Adaptive Framework for enhancing community Energy Resilience (SAFER) project advances resilience science and engineering by addressing challenges in rural Kansas communities where aging infrastructure, extreme weather, and socioeconomic disparities heighten vulnerability to energy disruptions. Traditional approaches often focus on technical performance while overlooking community concerns and priorities. SAFER responds by integrating community perspectives with advanced analytical frameworks to create a holistic model for measuring and improving resilience. Project objectives included developing novel resilience metrics, advancing modeling frameworks that capture interdependencies across infrastructures, and embedding community-centric indicators directly into planning processes for distributed energy resources. The key technical innovations included the creation of self-organizing map (SOM)-based indices for objective resilience quantification, hetero-functional graph theory (HFGT) models linking power, water, transportation, and community assets, and graph neural network (GNN) tools for identifying critical nodes in complex systems. Community-centric energy planning was demonstrated through optimal siting and sizing of (photovoltaic) PV and battery storage, ensuring resilience enhancements also addressed energy burden and energy insecurity. SAFER engaged community partners in Dodge City and Ford County through surveys, focus groups, and workshops, generating more than 600 responses that established baseline measures of energy burden, financial insecurity, and willingness-to-pay to avoid outages. This data, organized in terms of a community capitals framework, informed the development of weighted reliability indices that better reflect community costs than traditional utility metrics. SAFER’s GNN-based critical node identification framework identified expert-labelled critical nodes with over 99% accuracy, while also uncovering additional functionalities essential for proactive resilience planning. The project’s models demonstrated that optimal PV and storage deployment could improve resilience indices by over 11 percent, with dispatch strategies further enhancing outcomes, confirming both the technical effectiveness and economic feasibility of these approaches. Through its combined emphasis on rigorous modeling, community-focused planning, and community engagement, SAFER advances the state of resilience research while delivering direct benefits to rural communities. The project provides tools, guidelines, and resilience heatmaps that help utilities, local governments, and residents better anticipate disruptions, prioritize investments, and strengthen the capacity to withstand and recover from energy-related hazards. Furthermore, the developed HFG and GNN frameworks are designed for transferability, allowing them to be adapted for resilience planning in other communities with minimal retraining. This inductive learning capability provides a scalable pathway to extend the SAFER project’s impact. Thus, creating a foundation for a nationally applicable model of infrastructure resilience. Additionally, the HFG can also be extended to include other FEMA community lifelines.

14 SOLAR ENERGY↗

Exploring the frontiers of condensed-phase chemistry with a general reactive machine learning potential

Abstract Atomistic simulation has a broad range of applications from drug design to materials discovery. Machine learning interatomic potentials (MLIPs) have become an efficient alternative to computationally expensive ab initio simulations. For this reason, chemistry and materials science would greatly benefit from a general reactive MLIP, that is, an MLIP that is applicable to a broad range of reactive chemistry without the need for refitting. Here we develop a general reactive MLIP (ANI-1xnr) through automated sampling of condensed-phase reactions. ANI-1xnr is then applied to study five distinct systems: carbon solid-phase nucleation, graphene ring formation from acetylene, biofuel additives, combustion of methane and the spontaneous formation of glycine from early earth small molecules. In all studies, ANI-1xnr closely matches experiment (when available) and/or previous studies using traditional model chemistry methods. As such, ANI-1xnr proves to be a highly general reactive MLIP for C, H, N and O elements in the condensed phase, enabling high-throughput in silico reactive chemistry experimentation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Liquid state properties of SEI components in dimethoxyethane

The solid-electrolyte interphase (SEI) layer is a critical constituent of battery technology, which incorporates the use of lithium metals. Since the formation of the SEI is difficult to avoid, the engineering and harnessing of the SEI are absolutely critical to advancing energy storage. One problem is that much fundamental information about SEI properties is lacking due to the difficulty in probing a chemically complex interfacial system. One such property that is currently unknown is the dissolution of the SEI. This process can have significant effects on the stability of the SEI, which is critical to battery performance but is difficult to probe experimentally. Here, we report the use of ab initio computational chemistry simulations to probe the solution state properties of SEI components LiF, Li 2 O, LiOH, and Li 2 CO 3 in order to study their dissolution and other solution-based characteristics. Ab initio molecular dynamics was used to study the solvation structures of the SEI with a combination of radial distribution functions, discrete solvation structure maps, and vibrational density of states, which allows for the determination of free energies. From the change in free energy of dissolution, we determined that LiOH is the most likely component to dissolve in the electrolyte followed by LiF, Li 2 CO 3 , and Li 2 O although none were favored thermodynamically. Furthermore, this indicates that dissolution is not probable, but Li 2 O would make the most stable SEI with regard to dissolution in the electrolyte.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Library of Advanced Materials for Engineering (LAMÉ) 5.0

Accurate and efficient constitutive modeling remains a cornerstone issue for solid mechanics analysis. Over the years, the LAMÉ advanced material model library has grown to address this challenge by implementing models capable of describing material systems spanning soft polymers to stiff ceramics including both isotropic and anisotropic responses. Inelastic behaviors including (visco)plasticity, damage, and fracture have all incorporated for use in various analyses. This multitude of options and flexibility, however, comes at the cost of many capabilities, features, and responses and the ensuing complexity in the resulting implementation. Therefore, to enhance confidence and enable the utilization of the LAMÉ library in application, this effort seeks to document and verify the various models in the LAMÉ library. Specifically, the broader strategy, organization, and interface of the library itself is first presented. The physical theory, numerical implementation, and user guide for a large set of models is then discussed. Importantly, a number of verification tests are performed with each model to not only have confidence in the model itself but also highlight some important response characteristics and features that may be of interest to end-users. Finally, in looking ahead to the future, approaches to add material models to this library and further expand the capabilities are presented

36 MATERIALS SCIENCE↗

Library of Advanced Materials for Engineering (LAMÉ) 5.4

Accurate and efficient constitutive modeling remains a cornerstone issue for solid mechanics analysis. Over the years, the LAMÉ advanced material model library has grown to address this challenge by implementing models capable of describing material systems spanning soft polymers to stiff ceramics including both isotropic and anisotropic responses. Inelastic behaviors including (visco)plasticity, damage, and fracture have all incorporated for use in various analyses. This multitude of options and flexibility, however, comes at the cost of many capabilities, features, and responses and the ensuing complexity in the resulting implementation. Therefore, to enhance confidence and enable the utilization of the LAMÉ library in application, this effort seeks to document and verify the various models in the LAMÉ library. Specifically, the broader strategy, organization, and interface of the library itself is first presented. The physical theory, numerical implementation, and user guide for a large set of models is then discussed. Importantly, a number of verification tests are performed with each model to not only have confidence in the model itself but also highlight some important response characteristics and features that may be of interest to end-users. Finally, in looking ahead to the future, approaches to add material models to this library and further expand the capabilities are presented.

36 MATERIALS SCIENCE↗

Equation-of-Motion Coupled-Cluster Cumulant Green’s Function for Excited States and X-ray Spectra

Green’s function methods within many-body perturbation theory provide a general framework for treating electronic correlations in excited states and spectra. Conventional methods using the Dyson equation or the cumulant expansion are typically based on the GW self-energy approximation. In order to extend this approximation in molecular systems, a non-perturbative real-time coupled-cluster cumulant Green’s function approach has been introduced, where the cumulant is obtained as the solution to a set of coupled first order, non-linear differential equations. This approach naturally includes non-linear corrections to conventional cumulant Green’s function techniques where the cumulant is linear in the GW self-energy. The method yields the spectral function associated with the core Green’s function, which is directly related to the x-ray photoemission spectra (XPS) of molecular systems. The approach also yields very good results for binding energies and satellite excitations. The x-ray absorption spectrum (XAS) is then obtained as a convolution of the core spectral function and an effective one-body absorption spectrum. Here this approach is extended to include the full coupled-cluster-singles (CCS) core Green’s function by including the complete form of the non-linear contributions to the cumulant as well as all single, double, and triple cluster excitations in the CC amplitude equations. This approach naturally builds in orthogonality and shake-up effects analogous to those in the Mahan-Noizeres-de Dominicis edge singularity corrections that enhance the XAS near the edge. Themethod is illustrated for the XPS and XAS of NH 3 .

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

SAM: A Modern System Code for Advanced Non-LWR Safety Analysis

The System Analysis Module (SAM), developed at Argonne National Laboratory and by collaborators at other organizations, is for advanced non–light water reactor safety analysis. SAM aims to provide fast-running, modest-fidelity, whole-plant transient analysis capabilities that are essential for fast-turnaround design scoping and engineering analyses of advanced reactor concepts. To facilitate code development, SAM utilizes the MOOSE object-oriented application framework, its underlying finite element library, and linear and nonlinear solvers to leverage modern advanced software environments and numerical methods. SAM aims to solve tightly coupled physical phenomena, including fission reaction, heat transfer, fluid dynamics, and thermal-mechanical responses in advanced reactor structures, systems, and components with high accuracy and efficiency. Finally, this paper gives an overview of the SAM code development, including goals and functional requirements, physical models, current capabilities, verification and validation, software quality assurance, and examples of simulations for advanced nuclear reactor applications.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗