Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “system analysis module”

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

Recycling the Energy of Indoor Light: Highly Efficient Organic Photovoltaics via a Ternary Strategy

A ternary approach in organic photovoltaics (OPV) is simple and reliable to effectively tune the optical, electrical and morphological properties of the photoactive layer for high efficiency as well as long-term stability under 1 sun. However, there have been few papers reporting spectra and weak illumination compared to outdoor light. In this study, by using two compatible 16 structures, we demonstrate a simultaneous modulation of light absorption and molecular donor polymers (PM7 and PM7 D1) with slightly-different band gaps and similar chemical packing under ambient conditions, which resulted in efficient indoor OPVs exhibiting power conversion efficiency (PCE) over 20% under 1000 lux of warm white light-emitting diode (2900 K). From morphological analysis, we infer that PM7 serves to seed nucleation of PM7 D1 in the ternary blend, templating its crystallization and alignment along the PM7 backbone. Such templating effect leads to increased domain spacing and relative degree of crystallinity (rDoC) compared to those of each binary system. We further show that higher rDoC helps suppress both bimolecular and trap-assisted recombination of photogenerated charges in the ternary devices. As a result, the complementary absorption and synergistic molecular assembly of the two donor polymers enhance the short-circuit current density as to increase the average PCEs from 9.6 to 10.3% under 1 sun and from 18.7 to 20.0% under 1000 lux. We envision that our strategy of incorporating both a planar and a flexible donor polymer with similar chemical structures can be generally applicable to attain a high performance ternary OPV under both 1 sun and indoor light.

14 SOLAR ENERGY↗

A Unified Design Theory for Multi-Port Polyphase Transformers Enabling Scalable Power-Multiplexed EV Fleet Charging Systems

This paper presents a unified analytical design theory for multi-port polyphase transformers, targeting scalable and isolated high-power Electric Vehicle (EV) fleet charging systems with power multiplexing capability. As fleet electrification accelerates, conventional one-to-one charger architectures face significant challenges in infrastructure cost, peak power demand, and low utilization of installed power electronics. Power-multiplexed charging architectures, which dynamically distribute power from a shared pool of converter modules across multiple vehicles, have emerged as a promising solution. However, such architectures require scalable, isolated multi-port power interfaces capable of routing energy among multiple inputs and outputs, whose design remains complex and dependent on iterative modeling. To address this gap, the proposed theory provides closed-form expressions for self-inductance, leakage inductance, and mutual coupling terms for arbitrary multi-phase, multi-port transformer structures. The formulation enables direct synthesis of isolated multi-input and multi-output resonant converter systems without reliance on geometry-specific finite-element analysis or extensive parameter extraction. This capability is particularly critical for power-multiplexed systems, where modular converter structures must interface with multiple vehicles while maintaining galvanic isolation and flexible power allocation. The effectiveness of the proposed framework is demonstrated through the design of a 360 kW multi-phase system operating over a 700–900 VDC input and 400–1250 VDC output range. PLECS simulation results confirm accurate prediction of system behavior and validate the applicability of the approach to multi-port, power-multiplexed charging scenarios. The proposed method significantly reduces design complexity while enabling scalable, cost-effective, and fully utilized EV fleet charging infrastructure.

Asa, Erdem [ORNL] (ORCID:0000000190884812)↗

FORCE Update 2024

The Framework for Optimization of Resources and Economics (FORCE) tool suite is the U.S. Department of Energy’s Nuclear Integrated Energy Systems (IES) Program flagship tool suite for technoeconomic IES analysis of IES. This tool suite is useful for analysis designed to evaluate and improve the technoeconomics of energy production systems, particularly for systems including nuclear technology. In this report, we document the development activity for the FORCE tool suite to extend its capabilities as performed during fiscal year 2024. In addition to reliability and accessibility, capability is one of the three standards guiding the development of the FORCE tool suite and the software codes that are its constituent parts. Extending the capabilities of the FORCE tool suite allows analysis both within the IES program as well as industry, university, and laboratory partners to perform analysis with more accuracy, insight, and impactful narrative. Four areas of capability development were the focus of activity this year: economic parameter uncertainty quantification, multiresolution analysis, components-to-optimization workflow automation, and statespace construction workflows for real-time optimal control. In economic parameter uncertainty quantification, the ability of HERON to capture risk due to scenarios (weather and energy demand uncertainty) was expanded to also include uncertainties in financial parameters such as capital cost or operation and maintenance costs. By including these sources of uncertainty, which are sometimes very large compared with scenario uncertainty, HERON is better able to capture the risk posed by investment in various IES technology. Because of this, analysts can also consider the reduction in risks that can be realized by choice of some technologies. In multiresolution analysis, development activity extended on work completed previously. In fiscal year 2023, methods for decomposing time series signals, such as demand, solar and wind availability, and price profiles, were analyzed and down-selected to those most effective at splitting signals into different resolutions. These resolutions allow considering the influence of different energy demand and supply behaviors across different time scales. For example, energy demand might be divided into seasonal, weekly, and hourly profiles. In fiscal year 2024, this preliminary work was extended and implemented within the Risk Analysis Virtual Environment (RAVEN) risk and uncertainty analysis platform, which is used throughout the FORCE framework. This development of the “multi-resolution time series analysis” (MR-TSA) module in RAVEN allows training synthetic history generators on complex time series. These synthetic history generators can then be used in HERON for generating scenarios that represent possible market and weather scenarios that can be analyzed on different time scales. We envision completing this work in the future, implementing multiresolution dispatch optimization strategies that can make the most beneficial use of these stratified time histories. In components-to-optimization workflow development, workflows for translating user inputs of components into algorithms for algebraic optimization were selected and implemented. Similar algorithms within the Holistic Energy Resource Optimization Network (HERON) were separated from the main code base of HERON and gathered with the components-to-optimization workflows in the new Dispatch Optimization Variable Engine (DOVE) software library. This modularization allows FORCE users to analyze dispatch optimization and energy system duty cycles independently of HERON, which previously was a burdensome task. Additionally, these dispatch optimization algorithms, set up in an independent library, can now be used across all software applications within FORCE, especially including the real-time optimal control software Optimization of Real-time Capacity Allocation (ORCA). Allowing FORCE software to share dispatch optimization algorithms within a single library allows for improved software maintenance and reliability. In statespace characterization workflow development, alternative workflows for optimizing dispatch with additional technical accuracy was the focus, particularly to improve the real-time optimization decision making in ORCA. Using algorithms and workflows initially developed for the Feasible Actuator Range Modifier (FARM), workflows for determining the statespace representation of IES were identified and demonstrated. The resulting dispatch optimization required a more robust optimization algorithm than that originally used in HERON (and moved to DOVE), which required adding an alternate workflow to DOVE that can more accurately match the behavior of physical systems using a partial differential equation representation. In conclusion, capability developments in the FORCE tool suite in fiscal year 2024 have improved the ability of the FORCE tool suite to perform

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Capture Instead of Release: Defect-Modulated Radionuclide Leaching Kinetics in Metal–Organic Frameworks

Comparison of defect-controlled leaching-kinetics modulation of metal–organic frameworks (MOFs) and porous functionalized silica-based materials was performed on the example of a radionuclide and radionuclide surrogate for the first time, revealing an unprecedented readsorption phenomenon. On a series of zirconium-based MOFs as model systems, we demonstrated the ability to capture and retain >99% of the transuranic 241 Am radionuclide after 1 week of storage. Here, we report the possibility of tailoring radionuclide release kinetics in MOFs through framework defects as a function of postsynthetically installed organic ligands including cation-chelating crown ether-based linkers. Based on comprehensive analysis using spectroscopy (EXAFS, UV–vis, FTIR, and NMR), X-ray crystallography (single crystal and powder), and theoretical calculations (nine kinetics models and structure simulations), we demonstrated the synergy of radionuclide integration methods, topological restrictions, postsynthetic scaffold modification, and defect engineering. This combination is inaccessible in any other material and highlights the advantages of using well-defined frameworks for gaining fundamental knowledge necessary for the advancement of actinide-based material development, providing a pathway for addressing upcoming challenges in the nuclear waste administration sector.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

Structure of phospholipase Cε reveals an integrated RA1 domain and previously unidentified regulatory elements

Phospholipase Cε (PLCε) generates lipid-derived second messengers at the plasma and perinuclear membranes in the cardiovascular system. It is activated in response to a wide variety of signals, such as those conveyed by Rap1A and Ras, through a mechanism that involves its C-terminal Ras association (RA) domains (RA1 and RA2). However, the complexity and size of PLCε has hindered its structural and functional analysis. Herein, we report the 2.7 Å crystal structure of the minimal fragment of PLCε that retains basal activity. This structure includes the RA1 domain, which forms extensive interactions with other core domains. A conserved amphipathic helix in the autoregulatory X–Y linker of PLCε is also revealed, which we show modulates activity in vitro and in cells. The studies provide the structural framework for the core of this critical cardiovascular enzyme that will allow for a better understanding of its regulation and roles in disease.

59 BASIC BIOLOGICAL SCIENCES↗

Rubisco regulation in response to altered carbon status in the cyanobacterium Synechococcus elongatus PCC 7942

Abstract Photosynthetic organisms possess a variety of mechanisms to achieve balance between absorbed light (source) and the capacity to metabolically utilize or dissipate this energy (sink). While regulatory processes that detect changes in metabolic status/balance are relatively well studied in plants, analogous pathways remain poorly characterized in photosynthetic microbes. Here, we explored systemic changes that result from alterations in carbon availability in the model cyanobacterium Synechococcus elongatus PCC 7942 by taking advantage of an engineered strain where influx/efflux of a central carbon metabolite, sucrose, can be regulated experimentally. We observed that induction of a high-flux sucrose export pathway leads to depletion of internal carbon storage pools (glycogen) and concurrent increases in estimates of photosynthetic activity. Further, a proteome-wide analysis and fluorescence reporter-based analysis revealed that upregulated factors following the activation of the metabolic sink are concentrated on ribulose-1,5-bisphosphate carboxylase-oxygenase (Rubisco) and auxiliary modules involved in Rubisco maturation. Carboxysome number and Rubisco activity also increased following engagement of sucrose secretion. Conversely, reversing the flux of sucrose by feeding exogenous sucrose through the heterologous transporter resulted in increased glycogen pools, decreased Rubisco abundance, and carboxysome reorganization. Our data suggest that Rubisco activity and organization are key variables connected to regulatory pathways involved in metabolic balancing in cyanobacteria.

Plant Sciences↗

A tissue‐resolved, network‐based transcriptomic framework for abiotic stress responses in sorghum

Developing climate‐resilient crops requires a detailed understanding of stress‐induced gene expression dynamics, as maladaptive responses can compromise their productivity and survival. Sorghum, a globally important cereal with exceptional tolerance to multiple abiotic stresses, provides a powerful system for investigating these dynamics. However, how stress type, tissue specificity, and temporal progression jointly shape transcriptomic responses in crops remains poorly understood. Here, we present a comparative, time‐resolved transcriptomic atlas of sorghum responses to drought, heat, and salinity stress across shoot and root tissues. Integrative analyses revealed that tissue specificity is the dominant determinant of abiotic stress‐induced gene reprogramming across all three stresses. Building on these global comparisons, we focused on heat stress, as it elicited the most coherent and pronounced transcriptional and regulatory responses, enabling deeper network‐level interrogation. Co‐expression network analysis identified tissue‐specific modules enriched for phytohormone‐responsive genes, while gene regulatory network (GRN) mapping and cistrome analyses uncovered transcription factors (TFs) controlling key hub genes within these modules. Together, this study provides a foundational transcriptomic and network‐based resource for dissecting the regulatory architecture of abiotic stress responses in sorghum and offers prioritized candidates for future functional validation and engineering of climate‐resilient crops.

abiotic stress↗

MOOSE Reactor Module Meshing Enhancements to Support Reactor Physics Analysis

The U.S. Department of Energy Office of Nuclear Energy Advanced Modeling and Simulation (NEAMS) program develops an integrated suite of advanced reactor physics tools built upon the Multiphysics Object-Oriented Simulation Environment (MOOSE) framework. Each code generally requires an input finite element mesh on which the physics solution is calculated, reported, and transferred to other physics codes. The meshing process is often burdensome for the complex geometries present in reactors due to lack of easy-to-use, open-source meshing tools. To address the bottleneck associated with meshing complex geometries found in nuclear reactors, NEAMS initiated the development of the MOOSE Reactor Module starting in FY21. The Reactor Module builds off the existing MOOSE Mesh System to include targeted meshing capabilities such as the ability to generate hexagonal pin cells, assemblies with ducts, rotating control drums, cores, peripheral zones around a core, as well as the automatic labeling (“reporting IDs”) of pin, assembly, and planar zones to simplify post-processing of results. As a Physics Module in MOOSE, the Reactor Module is open-source, available with any MOOSE installation, directly compatible with MOOSE-based tools, and can be invoked from MOOSEbased applications to generate meshes. Functionality from the Reactor Module has been applied to several advanced reactor concepts to demonstrate user workflow improvements and accuracy. The primary objective of the Reactor Module is to improve useability of MOOSEbased tools by streamlining mesh generation and output inspection processes. During FY22, the functionality of the Reactor Module (and accompanying Mesh System) has been expanded based on user needs. First, the Reactor Geometry Mesh Builder capability developed primarily in FY21 has been refactored and merged to the public MOOSE repository. This capability wraps underlying Reactor Module mesh generators into a “Pin – Assembly – Core” workflow appropriate for conventional Cartesian and hexagonal geometries, and notably assigns material IDs during mesh generation stage and generates only the minimal number of blocks needed in order to reduce computational burden. Biasing and boundary layer options have been added to the base mesh generators as required by thermal hydraulics solvers. The reporting ID functionality has been expanded to differentiate ring-wise and azimuthal sectors within a pin for use with depletion algorithms, and VectorPostProcessor and Reporter objects are now available to integrate solution variables across zones based on ID combinations. Functionality to trim hexagonal meshes along the center or periphery has been developed so users may leverage symmetry and reflective boundary conditions to reduce the mesh size. A flexible and powerful tool to fill the space between two sidesets has been introduced to the framework and can be used for transition layers such as stitching two assemblies together with different numbers of pins, or for complex geometries which do not follow conventional Cartesian/hexagonal patterns. Finally, additional verification problems were performed with NEAMS physics tools in complement with existing NEAMS work.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Real-time steerable frequency-stepped Doppler backscattering (DBS) system for local helicon wave electric field measurements on the DIII-D tokamak

A new frequency-stepped Doppler backscattering (DBS) system has been integrated into a real-time steerable electron cyclotron heating launcher system to simultaneously probe local background turbulence (f < 10 MHz) and high-frequency (20–550 MHz) density fluctuations in the DIII-D tokamak. The launcher allows for 2D steering (horizontally and vertically) over wide angular ranges to optimize probe location and wavenumber response. The vertical steering can be optimized during a discharge in real time. The new DBS system employs a programmable frequency synthesizer with adjustable dwell time as a source to launch either O or X-mode polarized millimeter waves. This system can step in real-time over the entire E-band frequency range (60–90 GHz). This combination of capabilities allows for the diagnosis of the complex internal spatial structure of high power (>200 kW) helicon waves (476 MHz) injected from an external antenna during helicon current drive experiments in DIII-D. Broadband density fluctuations around the helicon frequency are observed during real-time scans of measurement location and wavenumber during these experiments. Analysis indicates that these broadband high-frequency fluctuations are a result of backscattering of the DBS millimeter-wave probe beam from plasma turbulence modulated by the helicon wave. It is observed that background turbulence is effectively locally “tagged” with the helicon wave electric field, forming images of the turbulent spectrum in the overall density fluctuation spectrum that appear as high-frequency sidebands of the turbulence. These observations of background turbulence and high-frequency fluctuations open up the possibility of monitoring local helicon wave amplitude by comparing the high-frequency signal amplitude to the simultaneously measured background turbulence. In combination with the real-time measurement location and wavenumber scanning capabilities (offered by real-time frequency-stepping and steering), this allows rapid determination of the spatial distribution of the helicon wave power during steady-state plasma operation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

System-Level Efficiency Study of Modular DC/DC and Fuel Cell Stack Series–Parallel Configurations

This paper presents a system-level efficiency study of modular DC/DC converter and fuel cell stack configurations under series, parallel, and series–parallel connections. The investigation considers selected DC/DC converter topologies, including isolated and non-isolated architectures, boost, non-inverting buck-boost, and resonant converters, integrated with commercially available fuel cell stacks, including the Accelera FCE150, Ballard FCgen-HPS, and Toyota TFCM2. DC/DC converter modules are evaluated in modular configurations rated at 60 kW and 90 kW and two distinct output voltage ranges, specifically 580–730V and 780–930V, examining how different interconnection schemes impact overall system efficiency. The converters are evaluated under full load (100%), partial load (66%), and light load (33%) conditions, providing a comprehensive assessment of efficiency and operational characteristics across varying power demands. Approximately two-thousand efficiency data points are obtained from laboratory prototype–level component measurements and validated design evaluations across multiple converter topologies, modular configurations, voltage ranges, and load conditions, providing a robust dataset for comparative system-level efficiency analysis. The results highlight the effects of modularity and topology selection on system-level efficiency, offering a “playbook” framework for designers to select appropriate DC/DC converter arrangements and fuel cell stack connections for series, parallel, or hybrid configurations based on efficiency considerations.

DC/DC↗

Water foraging with dynamic roots in E3SM; The role of roots in terrestrial ecosystem memory on intermediate timescales (Final Technical Report)

Terrestrial ecosystems can show sustained responses to stress events that may last for years after the initial perturbation. This phenomenon is referred to as legacy or memory and it emerges from a set of ecosystem processes that are poorly captured by Earth System and Land Surface Models. Consequently, these models struggle to predict the impact that extreme climate events have on surface energy, carbon and hydrological exchange once a stressor such as drought has relaxed or in response to repeated exposure to stress. In this project, we set out to test how the addition of dynamic root profiles in land surface models impact the capacity for Earth System Models, such as the Department of Energy’s E3SM, to capture realistic legacy effects. Observations have shown that vegetation shifts their root profiles during stress events to forage for water and these altered root profiles may take years to relax back to the initial state. We hypothesized that the alteration of the belowground root structures may be a source of terrestrial ecosystem legacy that is missing from models. To test this idea we undertook three core activities. (1) We developed a global-scale analysis of the impact of dynamic roots on ecosystem legacy building from a recently developed dynamic root module. (2) We developed new root dynamics in the DOE’s Energy Land Model (ELM) that allowed not only root profiles to shift their profile but also simultaneously alter carbon allocation to fine root pools. We then tested the impact of these new dynamics with intensive sensitivity analysis at four long-term AmeriFlux sites. (3) We undertook detailed isotopic analysis of tree rings from these AmeriFlux sites to assess ecophysiological and ecohydrological legacy to stress events that can be used to benchmark the sensitivity experiments with ELM. From these core activities, we report the following key findings. Firstly, on a global scale, the addition of dynamic roots led to chronically water-stressed ecosystems recovering faster to climate stress while wetter ecosystem showed enhanced legacy. This is because across ecosystems, stress events were almost universally associated with water shortages that led to the development of deeper root profiles. These deeper root profiles proved beneficial for recovery from drought stress. While the root dynamics did not universally improve the modeled representation of legacy it showed complex transient dynamics that emerge from the addition of dynamic roots. Secondly, the sensitivity analysis illustrated long term shifts in rooting depths away from the prescribed default profile suggesting that initializing of root profiles could benefit from spin-up simulations that converge on locally optimized root profiles. In addition, by enabling dynamic allocation some of the sites predicted unrealistically low allocation to roots (and high allocation to leaves). While these changes did not dramatically alter modeled gross primary productivity, they illustrated that without more sophisticated root processes in models, there is little penalty to dramatically disinvest in roots. Thirdly, the isotopic analysis of tree rings showed highly distinct legacy responses across species and sites. For example, T. canadensis showed reduced transpiration the year after stress events illustrated by sustained elevated $\delta^{18}$O. In contrast, A saccharum displayed elevated $\delta^{13}$C associated with reduced stomatal conductance in response to the previous years’ stress event. These geochemical signatures of legacy were present despite tree growth returning to normal the year after the stress. The results show the importance of species-level dynamics in legacy that are absent in modeling that assumes common traits within plant functional types. In summary, the work here established new avenues to explore root dynamics in models while illustrating how additional root processes are needed before dynamic carbon allocation can be implemented. Lastly, this project provided mentorship to a postdoctoral fellow, training for an early career scientist and multiple undergraduate students recruited from a minority serving institution.

54 ENVIRONMENTAL SCIENCES↗

Relating Aerial Infrared Thermography Defects to Photovoltaic Performance: Preprint

In this research, we examine the relationship between aerial IR defect analysis and photovoltaic (PV) performance data for twelve utility- and commercial-scale solar sites in the United States. To do this, we fuse the site diagram geoJSON's, aerial infrared thermography (aIRT) defect analyses, and associated inverter time series, allowing for a direct comparison between site defects and time series data. Defect analyses were provided by Zeitview, under its Solar Insights platform. Following the data fusion process, we look at the relationship between system performance and aIRT defects. We investigate the relationship between degradation and hotspot defects, as well as the relationship between AC power data and offline strings and misaligned modules. In general, system degradation was not affected by long-term or balance-of-system (BoS) defects as they occurred infrequently in the data set. However, for one system, a near statistically significant relationship (p-value=0.057) was found when comparing the degradation of inverter blocks with several multi-hotspot defects to all other inverter blocks without this particular defect. There was strong alignment when comparing short-term recoverable module defects such as stuck trackers and offline strings to time series data. In general, we found that when an inverter block has more than 80% of modules flagged for one of these defects, its AC power time data is flat-lined and the inverter block is not producing.

aerial inspection↗

Geothermal Power Systems Analysis: Outcome of Industry Stakeholders Workshop: Preprint

Geothermal cost and performance evaluation implemented via technoeconomic assessment (TEA) modeling is critical for the Department of Energy (DOE) and other geothermal industry stakeholders in assessing the current state of geothermal technologies and to identify existing hurdles to commercially viable geothermal development. The Geothermal Electricity Technology Evaluation Model (GETEM) is a major TEA tool used in estimating the economic feasibility and levelized cost of energy (LCOE) of conventional hydrothermal systems and enhanced geothermal systems (EGS). Since 2021, GETEM has been transitioning from an intricate spreadsheet model to a user-friendly tool within the System Advisor Model (SAM) developed by the National Renewable Energy Laboratory (NREL). Apart from enabling an expanded visibility of the geothermal model among other renewable resources, having GETEM in SAM has the advantage of simulation automation, better usability, updates tracking, active user inputs/feedback, and extended financial modeling. GETEM is used in developing supply curves for the Annual Technology Baseline (ATB). The ATB data are inputs to the Renewable Energy Potential (reV) and the Regional Energy Deployment System (ReEDS) models. The geothermal module in NREL’s reV model assesses the geothermal energy potential in the conterminous United States by defining the geospatial intersection of geothermal resources with existing grid infrastructure within the constraint of land use characteristics. The ReEDS model is a capacity expansion model used for simulating the long-term build-out and operation of the US generation and transmission system based on current energy costs and policies. To ensure enhanced representation of current industry trends in our model transitions and development, we organized a two-day virtual workshop to elicit geothermal industry stakeholder input and recommendations on our current approaches and assumptions on technoeconomic, resource assessment, and deployment scenarios modeling of geothermal technologies. Participants included developers, operators, investors, regulatory agencies, system modelers, national laboratory researchers, consultants, and other stakeholders. In this workshop, we gained stakeholder insights on current geothermal plant performance (i.e., capacity factors), updated drilling costs and learning curves, and next generation technologies such as closed loop and superhot rock geothermal. Other outcomes from this workshop and its impact on future geothermal development feasibility, resource availability, and capacity expansion studies are compiled and discussed.

Annual Technology Baseline↗

A parallel discrete dislocation dynamics/kinetic Monte Carlo method to study non-conservative plastic processes

Non-conservative processes play a fundamental role in plasticity and are behind important macroscopic phenomena such as creep, dynamic strain aging, loop raft formation, etc. In the most general case, vacancy-induced dislocation climb is the operating unit mechanism. While dislocation/vacancy interactions have been modeled in the literature using a variety of methods, the approaches developed rely on continuum descriptions of both the vacancy population and its fluxes. However, there are numerous situations in physics where point defect populations display heterogeneous concentrations and/or non-smooth kinetics. Here, a kinetic Monte Carlo (kMC) approach for modeling vacancy transport in response to arbitrary stress fields is used. Vacancies are treated as point particles and are coupled to the dislocation substructure representing a deformed material via an advection term defined by the local stress gradients. The stress fields and the dislocation substructure are evolved using a discrete dislocation dynamics (DDD) module. To extend the coupled model to the treatment of large systems, we have implemented it in the massively-parallel DDD code ParaDiS. To avoid numerical incompatibilities associated with merging deterministic (DDD) and stochastic (kMC) integration algorithms, we cast the entire elasto-plastic-diffusive problem within a single stochastic framework, taking advantage of a parallel kMC algorithm to evolve the system as a single event-driven process. The large-scale implementation enables the study of the evolution of a variety of dislocation-defect scenarios governed by non-conservative transport kinetics. After carrying out an exhaustive numerical and computational analysis of our parallel algorithm, we show results that emphasize situations where inhomogeneous vacancy dynamics are of relevance, and compare discrete kinetics to continuum solutions for several cases.

36 MATERIALS SCIENCE↗

Surface vs Homogeneous Organo-Hafnium Catalyst Ion-Pairing and Ligand Effects on Ethylene Homo- and Copolymerizations

Heterogeneous catalysts have long dominated polyethylene and polypropylene production, but understanding their catalysis is challenged by uncertainties in active site structures and percentages. Surface-bound organometallic catalysts are an emerging strategy to combine successful homogeneous catalysts having well-understood structures, relatively high percentages of active sites, and exceptional control of selectivity, with the attractions of surface catalysts, to transition promising homogeneous systems to large-scale heterogeneous ones. Nevertheless, surface-bound olefin polymerization catalysts typically produce ultrahigh M w ’s but with low activity and comonomer selectivity. Here, we report the systematic synthesis and characterization of a series of pyridylamido–Hf complexes and their corresponding surface catalysts chemisorbed on sulfated alumina (AlS) and zirconia (ZrS). Comparative ethylene homo- and 1-octene copolymerizations reveal similar activity and 1-octene selectivity trends in the homogeneous and heterogeneous systems. For the surface pyridylamido–Hf catalyst series, large variations in activity (up to 10×) and 1-octene incorporation (up to 28×) are achieved by ligand and support manipulation. Interestingly, while the homogeneous catalysts exhibit positive comonomer effects in ethylene/1-octene copolymerization, the surface catalysts behave oppositely. Extended X-ray absorption fine structure (EXAFS) reveals significantly elongated Hf···O bond distances vs typical Hf–O covalent bonds (2.06 vs 1.97 Å). Furthermore, density functional theory (DFT) analysis of the heterolytic ion pair separation enthalpies, olefin insertion energetics, and NBO/Bader charges also suggest electrostatic Hf cation–anionic support binding and catalytic patterns, which are modulated by the ion-pairing energetics and ligand architecture.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Dissecting subcycle interference in photoelectron holography

Multipath holographic interference in strong-field quantum tunnel ionization is key to revealing subangstrom attosecond dynamics for molecular movies. This critical subcycle motion is often obscured by longer timescale effects such as ring-shaped patterns that appear in above-threshold ionization (ATI). In the present work, we overcome this problem by combining two techniques in theory and experimental analysis: unit-cell averaging and time-filtering data and simulations. Together these suppress ATI rings and enable an unprecedented highly detailed quantitative match between strong-field ionization experiments in argon and the Coulomb-quantum orbit strong-field approximation (CQSFA) theory. Velocity map images reveal fine modulations on the holographic spiderlike interference fringes that form near the polarization axis. CQSFA theory traces this to the interference of three types of electron pathways. The level of agreement between experiment and theory allows sensitive determination of quantum phase differences and symmetries, providing an important tool for quantitative dynamical imaging in quantum systems.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE JGI Metagenome Workflow

The DOE Joint Genome Institute (JGI) Metagenome Workflow performs metagenome data processing, including assembly; structural, functional, and taxonomic annotation; and binning of metagenomic data sets that are subsequently included into the Integrated Microbial Genomes and Microbiomes (IMG/M) (I.-M. A. Chen, K. Chu, K. Palaniappan, A. Ratner, et al., Nucleic Acids Res, 49:D751–D763, 2021, https://doi.org/10.1093/nar/gkaa939) comparative analysis system and provided for download via the JGI data portal (https://genome.jgi.doe.gov/portal/). This workflow scales to run on thousands of metagenome samples per year, which can vary by the complexity of microbial communities and sequencing depth. Here, we describe the different tools, databases, and parameters used at different steps of the workflow to help with the interpretation of metagenome data available in IMG and to enable researchers to apply this workflow to their own data. We use 20 publicly available sediment metagenomes to illustrate the computing requirements for the different steps and highlight the typical results of data processing. The workflow modules for read filtering and metagenome assembly are available as a workflow description language (WDL) file (https://code.jgi.doe.gov/BFoster/jgi_meta_wdl). The workflow modules for annotation and binning are provided as a service to the user community at https://img.jgi.doe.gov/submit and require filling out the project and associated metadata descriptions in the Genomes OnLine Database (GOLD) (S. Mukherjee, D. Stamatis, J. Bertsch, G. Ovchinnikova, et al., Nucleic Acids Res, 49:D723–D733, 2021, https://doi.org/10.1093/nar/gkaa983).

59 BASIC BIOLOGICAL SCIENCES↗

Disruption of histone acetylation homeostasis reveals multilayered chromatin regulation for transcriptional resiliency

Background Epigenetic modifications, nucleosome occupancy, and three-dimensional chromatin architecture collectively create a multi-layered, highly interactive regulatory system for controlling genomic functionality. Dysregulation of epigenetic processes leads to a plethora of abnormalities including disease states. Therapies focused on epigenetic modulation can alter gene expression to correct dysfunction, though the perpetuation of these states and the relationships among chromatin regulatory layers is not well understood. Results Here, we investigated global and local chromatin structural and functional responses after acute histone deacetylase inhibitor treatment (suberoylanilide hydroxamic acid) in lung cancer cells across time. Treatment substantially increased global histone acetylation resulting in a pervasive but not distinctive signature. The spread of acetylation did not significantly impact global chromatin accessibility, and nucleosome remodeling largely occurred at finer scales in functionally relevant genomic regions. Indeed, both H 3 K 4 trimethylation, a mark of active transcription, and gene expression changes were altered in a controlled locus-specific manner, suggesting aberrant acetylation indirectly leads to balanced and bidirectional gene expression profiles from tighter regulation of other chromatin features. HDACi treatment induced (13%) genomic rearrangement in chromatin compartmentalization and moderate weakening of topologically associating domains. Conclusions Continuous wavelet analysis of these features demonstrates that scale-dependent, locus-specific factors influence the relationship between chromatin architecture and functional output, suggesting that regulation of transcription and nucleosome remodeling is not entirely (nor linearly) dependent upon large scale compartment exchange. Structural and functional responses are most pronounced early after treatment with partial persistence of differential local chromatin features and expression later in time; this highlights the plasticity of chromatin regulation, which may have implications for the efficacy of epigenetic treatments. These results demonstrate the effectiveness of multi-layered regulation of transcription: in resilient systems, disruption of one chromatin feature does not distort the regulation of other features in supporting a transcriptional program that allows for survival.

59 BASIC BIOLOGICAL SCIENCES↗