Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Direct coupling analysis”

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 595 records · Page 33

General framework for quantifying dissipation pathways in open quantum systems. III. Off-diagonal subsystem–bath couplings

This paper extends the previously reported theory of dissipation pathways [C. W. Kim and I. Franco, J. Chem. Phys. 160, 214111 (2024)] to incorporate off-diagonal subsystem–bath coupling, which is often required to model molecular systems where the environment directly influences transitions and couplings between subsystem states. We systematically derive master equations for both population transfer and dissipation into individual bath components, for which we also rigorously prove energy conservation and detailed balance. The approach is based on second-order perturbation theory with respect to the subsystem–bath couplings, whose form is not limited to any specific model. The accuracy of the developed method is tested by applying it to diverse model Hamiltonians involving linearly coupled harmonic oscillator baths and comparing the outcomes against the hierarchical equations of motion (HEOM) method. Overall, our method accurately quantifies the contributions of specific bath components to the overall dissipation while significantly reducing the computational cost compared to numerically exact methods such as HEOM, thus offering a path to examine how vibronic interactions steer non-adiabatic processes in realistic chemical systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Effective field theory of Stückelberg vector bosons

We explore the effective field theory of a vector field $X^μ$ that has a Stückelberg mass. The absence of a gauge symmetry for $X^μ$ implies Lorentz-invariant operators are constructed directly from Xμ. Beyond the kinetic and mass terms, allowed interactions at the renormalizable level include $X_μX^μH^†H, (X_μX^μ)^2$, and $X_μj^μ$, where $j^μ$ is a global current of the SM or of a hidden sector. We show that all of these interactions lead to scattering amplitudes that grow with powers of $\sqrt{s}/m_X$, except for the case of $X_μj^μ$ where $j^μ$ is a nonanomalous global current. The latter is well-known when $\textit{X}$ is identified as a dark photon coupled to the electromagnetic current, often written equivalently as kinetic mixing between $\textit{X}$ and the photon. The power counting for the energy growth of the scattering amplitudes is facilitated by isolating the longitudinal enhancement. We examine in detail the interaction with an anomalous global vector current $X_μj^μ_{\text{anom}}$, carefully isolating the finite contribution to the fermion triangle diagram. We calculate the longitudinally-enhanced observables $Z → X_γ$ (when $m_X < m_Z), f\bar{f} → X_γ$, and $Z_γ → Z_γ$ when $\textit{X}$ couples to the baryon number current. Introducing a "fake" gauge-invariance by writing $X^μ = A^μ – ∂^μπ/m_X$, the would-be gauge anomaly associated with $A_μj^μ_{\text{anom}}$ is canceled by $j^μ_{\text{anom}}∂_μπ/m_X$; this is the four-dimensional Green-Schwarz anomaly-cancellation mechanism at work. Our analysis demonstrates there is a much larger set of possible interactions that an EFT with a Stückelberg vector field can have, revealing scattering amplitudes that grow with energy. The growth of these amplitudes can be tamed by a dark Higgs sector, but this requires dark Higgs boson interactions (and reintroduces fine-tuning in the dark Higgs sector) that can be separated from X interactions only in the limit $g \ll 1$ .

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Rapid dissolution without elemental fractionation by laser driven hydrothermal processing

Traditional dissolution of geologic samples often requires a significant time investment. Here, we present an alternative method for the dissolution of geologic materials using laser-driven hydrothermal processing (LDHP). LDHP uses laser energy directed onto a submerged sample, which increases the temperature and pressure at the liquid–sample interface and drives the hydrothermal dissolution coupled with photomechanical spallation, an ablative process. This uses focused 527 nm laser energy at 40 W average power, 1 kHz pulse repetition rate, and 115 ns pulse duration. Importantly, when LDHP is performed on basalt geostandards (BCR-2 and BHVO-2) using the conditions outlined, we show that LDHP does not produce significant elemental fractionation and, thus, can be considered an alternative processing method to traditional mechanical crushing and acid digestion. Additionally, it is possible using LDHP to utilize the spatially confined beam to target and selectively isolate individual phases in a rock, potentially alleviating the need for mechanical separation of inclusions that are difficult to physically isolate. Furthermore, using this outlined method of LDHP, we demonstrate full dissolution of 120 mg of obsidian in 85 minu, meaning that LDHP is a potentially very useful method when sample processing is time sensitive.

58 GEOSCIENCES↗

A Realistic Full-Scale 3D Modeling of Turning Using Coupled Smoothed Particle Hydrodynamics and Finite Element Method for Predicting Cutting Forces

Computational modelling is an effective technique for understanding the complex physics of machining. Large deformations, material separation, and high computational requirements are the key challenges faced while simulating machining. This work introduces a full-scale three-dimensional model of turning operations using a combined approach based on the Smoothed Particle Hydrodynamics (SPH) and Finite Element (FE) methods. By exploiting the advantages of each method, this approach leads to high-fidelity coupled SPH-FE machining models. Cutting forces and chip morphology are the primary results of interest. The machining models are validated with the results of turning experiments. Two-dimensional machining model underpredicts the cutting force and feed force by approximately 49% and 70%, respectively. Moreover, passive force cannot be predicted using the two-dimensional model. On the other hand, with the three-dimensional models developed in this manuscript, the difference between the total simulated force and experimentally measured force is ∼17%. The chip morphologies correlate with experiments in terms of the direction of the chip movement and the “long” continuous chips observed while turning Al 6061. This work expands the realm of machining simulations from two-dimensional orthogonal machining or sectional three-dimensional model to a full-scale realistic simulation. The encouraging simulation results show the potential to study more complex phenomena, such as machining stability and tool path modulation.

42 ENGINEERING↗

Charting C–C coupling pathways in electrochemical CO 2 reduction on Cu(111) using embedded correlated wavefunction theory

The electrochemical CO 2 reduction reaction (CO 2 RR) powered by excess zero-carbon-emission electricity to produce especially multicarbon (C 2+ ) products could contribute to a carbon-neutral to carbon-negative economy. Foundational to the rational design of efficient, selective CO 2 RR electrocatalysts is mechanistic analysis of the best metal catalyst thus far identified, namely, copper (Cu), via quantum mechanical computations to complement experiments. Here, we apply embedded correlated wavefunction (ECW) theory, which regionally corrects the electron exchange-correlation error in density functional theory (DFT) approximations, to examine multiple C–C coupling steps involving adsorbed CO (*CO) and its hydrogenated derivatives on the most ubiquitous facet, Cu(111). We predict that two adsorbed hydrogenated CO species, either *COH or *CHO, are necessary precursors for C–C bond formation. The three kinetically feasible pathways involving these species yield all three possible products: *COH–CHO, *COH–*COH, and *OCH–*OCH. The most kinetically favorable path forms *COH–CHO. In contrast, standard DFT approximations arrive at qualitatively different conclusions, namely, that only *CO and *COH will prevail on the surface and their C–C coupling paths produce only *COH–*COH and *CO–*CO, with a preference for the first product. This work demonstrates the importance of applying qualitatively and quantitatively accurate quantum mechanical method to simulate electrochemistry in order ultimately to shed light on ways to enhance selectivity toward C 2+ product formation via CO 2 RR electrocatalysts.

30 DIRECT ENERGY CONVERSION↗

Evidence of cryptic methane cycling and non-methanogenic methylamine consumption in the sulfate-reducing zone of sediment in the Santa Barbara Basin, California

The recently discovered cryptic methane cycle in the sulfate-reducing zone of marine and wetland sediment couples methylotrophic methanogenesis to anaerobic oxidation of methane (AOM). Here we present evidence of cryptic methane cycling activity within the upper regions of the sulfate-reducing zone, along a depth transect within the Santa Barbara Basin, off the coast of California, USA. The top 0–20 cm of sediment from each station was subjected to geochemical analyses and radiotracer incubations using 35 S–SO$^{2-}_{4}$, 14 C–mono-methylamine, and 14 C–CH 4 to find evidence of cryptic methane cycling. Methane concentrations were consistently low (3 to 16 µM) across the depth transect, despite AOM rates increasing with decreasing water depth (from max 0.05 nmol cm -3 d -1 at the deepest station to max 1.8 nmol cm -3 d -1 at the shallowest station). Porewater sulfate concentrations remained high (23 to 29 mM), despite the detection of sulfate reduction activity from 35 S–SO$^{2-}_{4}$ incubations with rates up to 134 nmol cm -3 d -1 . Metabolomic analysis showed that substrates for methanogenesis (i.e., acetate, methanol and methylamines) were mostly below the detection limit in the porewater, but some samples from the 1–2 cm depth section showed non-quantifiable evidence of these substrates, indicating their rapid turnover. Estimated methanogenesis from mono-methylamine ranged from 0.2 to 0.5 nmol cm -3 d -1 . Discrepancies between the rate constants (k) of methanogenesis (from 14 C–mono-methylamine) and AOM (from either 14 C–mono-methylamine-derived 14 C–CH 4 or from directly injected 14 C–CH 4 ) suggest the activity of a separate, concurrent metabolic process directly metabolizing mono-methylamine to inorganic carbon. We conclude that the results presented in this work show strong evidence of cryptic methane cycling occurring within the top 20 cm of sediment in the Santa Barbara Basin. The rapid cycling of carbon between methanogenesis and methanotropy likely prevents major build-up of methane in the sulfate-reducing zone. Furthermore, our data suggest that methylamine is utilized by both methanogenic archaea capable of methylotrophic methanogenesis and non-methanogenic microbial groups. We hypothesize that sulfate reduction is responsible for the additional methylamine turnover, but further investigation is needed to elucidate this metabolic activity.

54 ENVIRONMENTAL SCIENCES↗

Review—In Situ X-ray and Infrared Spectroscopic Studies of Electrochemical Systems

Despite of intense research and a wealth of data, the phenomena occurring during electrocatalysis are still a major obstacle in many chemical processes. Molecular analysis of the electrode/electrolyte interface is needed to correctly describe the reaction through identifying the species involved, their interaction with the environment and kinetics in situ, i.e. while the reaction is taking place. That can be done by coupling the electrochemical system with complementary non-electrochemical techniques. Particularly revealing are in situ X-ray spectroscopic techniques to analyze the electrode itself, providing the information on the changes in the catalyst during the reaction. The synergy of the traditional electrochemical techniques with the complementary spectroscopic methodologies offer understanding of the electrode/electrolyte interface above and beyond traditional experimental mainframe. Here we demonstrate how in situ X-ray absorption spectroscopy (XAS), in situ infrared reflection/absorption spectroscopy (IRRAS), and traditional voltammetric studies can increase our understanding of electrochemical processes during oxidation of ethanol. The results show the pronounced role of electrode surface in determining reaction kineticks and revealed the selectivity of the catalyst to complete oxidation pathway. They further provide understanding of the parameters that enhance its oxidation for future designing catalysts for alcohol oxidation fuel cells.

30 DIRECT ENERGY CONVERSION↗

Design Optimization of A generic Fissile Solution for Mo99 production using Electron Beam-based Neutron Generator using MCNP+CFD

Molybdenum-99 ( 99 Mo) is a medical isotope used in 80% of all medical imaging procedures today. However, 99 Mo is not directly administered to patients for imaging; its decay product, technetium-99m ( 99 mTc), is a metastable isotope with a half-life of ~6 hours. Such a short half-life does not allow for production, separation, and application to a patient before decaying further. Therefore, the solution is to produce the precursor to 99 mTc, 99 Mo. 99 Mo has a half-life of 66 hours, long enough to produce the material and send it to the desired location before it fully decays. Under National Nuclear Security Administration’s (NNSA) highly enriched uranium (HEU) minimization mission, the Material Management and Minimization (M3) office leads the Molybdenum-99 (Mo-99) program. A major goal of the Mo-99 program is to develop methods of producing 99 Mo using low enriched uranium (LEU) opposed to the traditionally used HEU domestically in the U.S. Internationally, LEU is the new standard for use in research reactors and isotope production facility. In this study, we introduce a generic fissile solution system for 99 Mo production facility meeting the design requirement of 1) LEU usage as fissile material and 2) subcritical fissile system coupled with a horizontal electron beam (E-beam) accelerator. The focus of this research is to optimize a fissile solution configuration for maximum 99 Mo production yield using MCNP model, and ultimately implement the optimized design with computational fluid dynamics (CFD) model to understand thermal hydraulic behavior and solution convection characteristic. The report consists of two parts: 1) isotope yield calculation and criticality analysis with a wide range of parametric test matrix, and 2) thermal hydraulic analysis for solution mixing and cooling assessment.

60 APPLIED LIFE SCIENCES↗

Investigating instabilities in magnetized low-pressure capacitively coupled RF plasma using particle-in-cell (PIC) simulations

The effect of a uniform magnetic field on particle transport in low-pressure radio frequency (RF) capacitively coupled plasma (CCP) has been studied using a particle-in-cell model. Three distinct regimes of plasma behavior can be identified as a function of the magnetic field. In the first regime at low magnetic fields, asymmetric plasma profiles are observed within the CCP chamber due to the effect of $\overrightarrow{E}$ × $\overrightarrow{B}$ drift. As the magnetic field increases, instabilities develop and form self-organized spoke-shaped structures that are distinctly seen within the bulk plasma closer to the sheath. In this second regime, the spoke-shaped coherent structures rotate inside the plasma chamber in the -$\overrightarrow{E}$ × $\overrightarrow{B}$ direction, where $\overrightarrow{E}$ × $\overrightarrow{B}$ are the DC electric and magnetic field vectors, respectively, and the DC electric field exists in the sheath and pre-sheath regions. The spoke rotation frequency is in the megahertz range. As the magnetic field strength increases further, the rotating coherent spokes continue to exist near the sheath. The coherent structures are, however, accompanied by new small-scale incoherent structures originating and moving within the bulk plasma region away from the sheath. This is the third regime of plasma behavior. The threshold values of the magnetic field between these regimes were found not to vary with changing plasma reactor geometry (e.g., area ratio between ground and powered electrodes) or the use of an external capacitor between the RF-powered electrode and the RF source. The threshold values of the magnetic field between these regimes shift toward higher values with increasing gas pressure. Analysis of the results indicates that the rotating structures are due to the lower hybrid instability driven by density gradients and electron-neutral collisions. This paper provides guidance on the upper limit of the magnetic field for instability-free operation in low-pressure CCP-based semiconductor deposition and etch systems that use the external magnetic field for plasma uniformity control.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Exploring new physics with O(keV) electron recoils in direct detection experiments

Motivated by the recent XENON 1T results, we explore various new physics models that can be discovered through searches for electron recoils in O(keV)-threshold direct-detection experiments. First, we consider the absorption of axion-like particles, dark photons, and scalars, either as dark matter relics or being produced directly in the Sun. In the latter case, we find that keV mass bosons produced in the Sun provide an adequate fit to the data but are excluded by stellar cooling constraints. We address this tension by introducing a novel Chameleon-like axion model, which can explain the excess while evading the stellar bounds. We find that absorption of bosonic dark matter provides a viable explanation for the excess only if the dark matter is a dark photon or an axion. In the latter case, photophobic axion couplings are necessary to avoid X-ray constraints. Second, we analyze models of dark matter-electron scattering to determine which models might explain the excess. Standard scattering of dark matter with electrons is generically in conflict with data from lower-threshold experiments. Momentum-dependent interactions with a heavy mediator can fit the data with dark matter mass heavier than a GeV but are generically in tension with collider constraints. Next, we consider dark matter consisting of two (or more) states that have a small mass splitting. The exothermic (down)scattering of the heavier state to the lighter state can fit the data for keV mass splittings. Finally, we consider a subcomponent of dark matter that is accelerated by scattering off cosmic rays, finding that dark matter interacting though an O(100 keV)-mass mediator can fit the data. The cross sections required in this scenario are, however, typically challenged by complementary probes of the light mediator. Throughout our study, we implement an unbinned Monte Carlo analysis and use an improved energy reconstruction of the XENON 1T events.

79 ASTRONOMY AND ASTROPHYSICS↗

End-Use Savings Shapes Measure Documentation: Console Water-to-Air Geothermal Heat Pump

Executive Summary Building on the successfully completed effort to calibrate and validate the U.S. Department of Energy’s ResStock™ and ComStock™ models over the past several years, the objective of this work is to produce national datasets that empower analysts working for federal, state, utility, city, and manufacturer stakeholders to answer a broad range of analysis questions. The goal of this work is to develop energy efficiency, electrification, and demand flexibility end-use load shapes (electricity, gas, propane, or fuel oil) that cover most of the high-impact, market-ready (or nearly market-ready) measures. An end-use savings shape is the difference in energy consumption between a baseline building and a building with an energy efficiency, electrification, or demand flexibility measure applied. It results in a time series profile that is broken down by end use and fuel (electricity or on-site gas, propane, or fuel oil use) at each time step. ComStock is a highly granular, bottom-up model that uses multiple data sources, statistical sampling methods, and advanced building energy simulations to estimate the annual sub-hourly energy consumption of the commercial building stock across the United States. The baseline model intends to represent the U.S. commercial building stock as it existed in 2018. This measure models the conversion of an existing heating, ventilating, and air-conditioning (HVAC) system to a series of “console” water-to-air heat pumps served by a ground heat exchanger. Console water-to-air geothermal heat pumps (GHP) are all-in-one packages that have no or minimal ductwork and serve individual spaces. Properly designed ground heat exchanger-coupled systems can offer benefits in energy efficiency relative to “conventional” HVAC systems, as well as facilitating beneficial electrification. Console GHPs can be coupled to a ground loop on the source side and can directly replace electric baseboard heaters or air-source packaged terminal heat pumps. Console GHPs can also bring in outdoor air for ventilation. This measure will be referred to throughout the document as the “Console GHP” upgrade.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Mesoscopic structure of mixed type domain walls in multiaxial ferroelectrics

The structure of a 180° uncharged rotational domain wall in a multiaxial ferroelectric film is studied in the framework of an analytical Landau-Ginzburg-Devonshire (LGD) approach. Finite element modeling (FEM) is used to solve numerically the system of the coupled nonlinear Euler-Lagrange (EL) second-order differential equations for two components of polarization. We show that the structure of the domain wall and corresponding metastable or stable phase of the film are controlled by a single parameter—the dimensionless ferroelectric anisotropy factor μ. We fitted the static profile of a solitary domain wall, calculated by FEM, with kinklike functions for polarization components, and extracted the five μ-dependent parameters from the fitting to FEM curves. The surprisingly high accuracy of the fitting results for two polarization components in the entire μ range allows us to conclude that the analytical functions, which are trial functions in the direct variational method, can be treated as a high-accuracy variational solution of the static EL equations. We further derive exact two-component analytical solutions of the static EL equations for a polydomain 180° domain structure in a multiaxial ferroelectric film. Using these, we derive analytical expressions for the system free energy and analyze its dependence on the film thickness and boundary conditions at the film surfaces. The single-domain state is ground for zero polarization derivative at the surfaces, while the polydomain states minimize the system energy for zero polarization at the surfaces. Counterintuitively, the energy of the polydomain states split into two levels, “0” and “1,” for zero polarization at the surfaces, and each of the levels contains a large number of close-energy sublevels, whose structure is characterized by a different number and type of domain walls. The analytical solutions can become a useful tool for Bayesian analysis of high-resolution scanning transmission electron microscopy images in ferroelectric films.

36 MATERIALS SCIENCE↗

The scalar chemical potential in cosmological collider physics

Non-analyticity in co-moving momenta within the non-Gaussian bispectrum is a distinctive sign of on-shell particle production during inflation, presenting a unique opportunity for the “direct detection” of particles with masses as large as the inflationary Hubble scale ( H ). However, the strength of such non-analyticity ordinarily drops exponentially by a Boltzmann-like factor as masses exceed H . In this paper, we study an exception provided by a dimension-5 derivative coupling of the inflaton to heavy-particle currents, applying it specifically to the case of two real scalars. The operator has a “chemical potential” form, which harnesses the large kinetic energy scale of the inflaton, \( {\overset{\cdot }{\phi}}_0^{1/2}\approx 60H \) , to act as an efficient source of scalar particle production. Derivative couplings of inflaton ensure radiative stability of the slow-roll potential, which in turn maintains (approximate) scale-invariance of the inflationary correlations. We show that a signal not suffering Boltzmann suppression can be obtained in the bispectrum with strength f NL ~ \( \mathcal{O} \) (0 . 01–10) for an extended range of scalar masses \( \lesssim {\overset{\cdot }{\phi}}_0^{1/2} \) , potentially as high as 10 15 GeV, within the sensitivity of upcoming LSS and more futuristic 21-cm experiments. The mechanism does not invoke any particular fine-tuning of parameters or breakdown of perturbation-theoretic control. The leading contribution appears at tree-level , which makes the calculation analytically tractable and removes the loop-suppression as compared to earlier chemical potential studies of non-zero spins. The steady particle production allows us to infer the effective mass of the heavy particles and the chemical potential from the variation in bispectrum oscillations as a function of co-moving momenta. Our analysis sets the stage for generalization to heavy bosons with non-zero spin.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Spatially Accelerated Winding Numbers for Curved Geometry

The generalized winding number (GWN) is a scalar field that supports robust containment queries on curved geometry, including non-watertight, overlapping, and nested boundary representations. While queries can be easily parallelized over samples, direct evaluation on parametric curves and surfaces remains costly for large and complex models. Fast, state-of-the-art GWN approaches leverage a spatial index to approximate the GWN, typically coupled with a Taylor expansion which approximates the GWN contribution for far clusters of geometric primitives. However, such methods operate only on discrete inputs such as triangle meshes and point clouds, and would introduce containment errors near boundaries if applied to curved input. We extend support for fast GWN evaluation over arbitrary collections of NURBS curves in 2D and trimmed NURBS patches in 3D via a Bounding Volume Hierarchy that stores efficiently precomputed moment data in the hierarchy nodes. When querying the hierarchy, approximations for far clusters are used alongside direct evaluation for nearby NURBS primitives, achieving sub-linear complexity while preserving the geometric features in the vicinity of the query point. Central to our performance improvements is an adaptive subdivision strategy for NURBS primitives during a preprocessing phase, creating better spatial partitions while retaining the same accuracy for containment decisions as a direct evaluation. We demonstrate the performance and accuracy of our approach across a large collection of 2D and 3D datasets.

Computer science↗

Measurement of Volatile Compounds for Real-Time Analysis of Soil Microbial Metabolic Response to Simulated Snowmelt

Snowmelt dynamics are a significant determinant of microbial metabolism in soil and regulate global biogeochemical cycles of carbon and nutrients by creating seasonal variations in soil redox and nutrient pools. With an increasing concern that climate change accelerates both snowmelt timing and rate, obtaining an accurate characterization of microbial response to snowmelt is important for understanding biogeochemical cycles intertwined with soil. However, observing microbial metabolism and its dynamics non-destructively remains a major challenge for systems such as soil. Microbial volatile compounds (mVCs) emitted from soil represent information-dense signatures and when assayed non-destructively using state-of-the-art instrumentation such as Proton Transfer Reaction-Time of Flight-Mass Spectrometry (PTR-TOF-MS) provide time resolved insights into the metabolism of active microbiomes. In this study, we used PTR-TOF-MS to investigate the metabolic trajectory of microbiomes from a subalpine forest soil, and their response to a simulated wet-up event akin to snowmelt. Using an information theory approach based on the partitioning of mutual information, we identified mVC metabolite pairs with robust interactions, including those that were non-linear and with time lags. The biological context for these mVC interactions was evaluated by projecting the connections onto the Kyoto Encyclopedia of Genes and Genomes (KEGG) network of known metabolic pathways. Simulated snowmelt resulted in a rapid increase in the production of trimethylamine (TMA) suggesting that anaerobic degradation of quaternary amine osmo/cryoprotectants, such as glycine betaine, may be important contributors to this resource pulse. Unique and synergistic connections between intermediates of methylotrophic pathways such as dimethylamine, formaldehyde and methanol were observed upon wet-up and indicate that the initial pulse of TMA was likely transformed into these intermediates by methylotrophs. Increases in ammonia oxidation signatures (transformation of hydroxylamine to nitrite) were observed in parallel, and while the relative role of nitrifiers or methylotrophs cannot be confirmed, the inferred connection to TMA oxidation suggests either a direct or indirect coupling between these processes. Overall, it appears that such mVC time-series from PTR-TOF-MS combined with causal inference represents an attractive approach to non-destructively observe soil microbial metabolism and its response to environmental perturbation.

59 BASIC BIOLOGICAL SCIENCES↗

Advanced Reactor Designs Security Analysis, Risk, and Recommendations: Risks, Consequences, and Possible by-Design Mitigation Approaches Associated with Select Advanced Reactors

Next-generation advanced reactors (ARs) incorporate enhanced safety systems, have smaller source terms, and feature compact modular designs, which should lessen their collective risk profiles. However, to fully evaluate risk, security needs to be a part of the equation. Without taking security into consideration, safety systems and components in the new ARs may be vulnerable to sabotage. These base attributes, coupled with enhanced security features specific to AR design through sound engineering and security-by-design (SeBD), should provide developers and operators with lower inherent security risk profiles. Building security early into the AR design may remove or passively secure potential critical targets from an adversary’s reach , thereby increasing overall safety and security. An integrated approach and diverse design team that includes engineering, operations, and security experts are fundamental to building security into the design without sacrificing fundamental operational efficiencies and principles. The objective of this project was to evaluate the security and safety interfaces for five classes of reactors, identify potential security vulnerabilities of structures, systems, and components (SSC), and underscore the need to consider security alongside safety in the design o f these concepts. The five reactor classes evaluated in this project and presented in this report are molten-salt reactors (MSR), high temperature gas reactors (HTGR), sodium-fast reactors (SFR), advanced light-water reactors (ALWR), and microreactors. These designs were selected because they reflect the concepts that are closest to market deployment and have received significant resource investments from the public and private sector. This project assesses the inherent security risks posed by common classes of ARs, provides a methodology and framework to assess security along with safety, and offers an analysis of potential mitigation strategies that could be incorporated. For each AR technology, the SSCs that relate to radionuclide source safety functions are discussed to understand the SSC contribution to safety and relative importance in the protective strategy for the design. The assumptions that went into evaluating each reactor concept originated from generic publicly available nonproprietary information and should not directly be used to qualify an absolute risk profile nor to rank specific AR designs. Instead, the purpose of the analysis is to understand and compare the generic inherent security risks of different AR technologies.

98 - NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL↗

DOE new players Carbon cycle (2016-2021)

The overarching scientific goals of our multidisciplinary grant was to further expand understanding about the key microorganisms (players), metabolic strategies (processes), and interspecies relationships (interactions) involved in the formation and oxidation of methane in the environment. This research applied novel environmental metagenomics, transcriptomics, and proteomic techniques, state-of-the-art analytical imaging, stable isotope geochemistry, and reaction-transport modeling to address these goals and develop an ‘ecosystems level’ understanding of the factors which regulate microbial methane cycling in anoxic sedimentary ecosystems. For decades, it was believed that the obligate step in methanogenesis catalyzed by methyl coenzyme M reductase (Mcr) was limited to a specific branch of the archaeal Domain, formerly known as the Euryarchaeota. Less than a decade ago co-I Tyson’s team published a surprising metagenomic-based discovery of divergent Mcr genes in a novel uncultured phylum (Bathyarchaeota), catalyzing a major shift in thinking about the diversity of microorganisms that encode the potential for methane (or higher alkane) metabolism in anoxic environments (Evans et al., 2015). In our work here, we further expand on the groups of archaea harboring the genomic potential for methane or hydrocarbon metabolism using environmental metagenomics and new gene targeted bioinformatics techniques. We additionally advanced understanding about the terminal electron acceptors and metabolic potential supporting the anaerobic oxidation of methane (AOM) in terrestrial ecosystems, specifically focused on new lineages of ANME archaea capable of respiring manganese oxides with methane presumably using large extracellular multi-heme cytochrome complexes. New details about specific syntrophic mechanisms underlying the exchange of electrons during sulfate-coupled methane oxidation between ANME-2 archaea and their sulfate-reducing bacterial partners were also elucidated as part of this funded project. Through a series of experimental ‘omics and single cell stable isotope probing studies with incubated environmental sediment samples and a cultured model electrogenic microorganism combined with model-based predictions. Combined, this work provides strong support for the hypothesis of direct interspecies electron transfer (DIET) is the dominant syntrophic mechanism controlling the anaerobic oxidation of methane with sulfate over other proposed mechanisms and additionally illustrates important spatial constraints and the underlying physico-chemical factors influencing AOM syntrophic consortia structure for DIET and extracellular metal respiration. This collaborative multi-institutional project successfully advanced several of our milestone goals including the identification of new microbial players containing methyl coenzyme M reductases hypothesized to be central to methane or hydrocarbon cycling in anoxic environments and enhancing fundamental knowledge about the role extracellular electron transfer plays in the ecophysiology of methanotrophic archaea respiring metal oxides and in the physical and metabolic structuring of syntrophic interactions in methane-rich sedimentary ecosystems.

03 NATURAL GAS↗

Revisiting Understanding of Electrochemical CO 2 Reduction on Cu(111): Competing Proton-Coupled Electron Transfer Reaction Mechanisms Revealed by Embedded Correlated Wavefunction Theory

Copper (Cu) electrodes, as the most efficacious of CO 2 reduction reaction (CO 2 RR) electrocatalysts, serve as prototypes for determining and validating reaction mechanisms associated with electrochemical CO 2 reduction to hydrocarbons. As in situ electrochemical mechanism determination by experiments is still out of reach, such mechanistic analysis typically is conducted using density functional theory (DFT). The semilocal exchange-correlation (XC) approximations most often used to model such catalysis unfortunately engender a basic error: predicting the wrong adsorption site for CO (a key CO 2 RR intermediate) on the most ubiquitous facet of Cu, namely, Cu(111). This longstanding inconsistency casts lingering doubt on previous DFT predictions of the attendant CO 2 RR kinetics. Here, we apply embedded correlated wavefunction (ECW) theory, which corrects XC functional error, to study the CO 2 RR on Cu(111) via both surface hydride (*H) transfer and proton-coupled electron transfer (PCET). We predict that adsorbed CO (*CO) reduces almost equally to two intermediates, namely, hydroxymethylidyne (*COH) and formyl (*CHO) at –0.9 V vs the RHE. In contrast, semilocal DFT approximations predict a strong preference for *COH. With increasing applied potential, the dominance of *COH (formed via potential-independent surface *H transfer) diminishes, switching to the competitive formation of both *CHO and *COH (both formed via potential-dependent PCET). Furthermore, our results also demonstrate the importance of including explicitly modeled solvent molecules in predicting electron-transfer barriers and reveal the pitfalls of overreliance on simple surface *H transfer models of reduction reactions.

30 DIRECT ENERGY CONVERSION↗