Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “predictive synthesis”

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,315 records · Page 73

Effective medium approximation for the refractive index of stratified metal oxide composites synthesized by atomic layer deposition

Atomic layer deposition (ALD) is a unique method for synthesizing conformal layers with precise composition. It is especially useful for the synthesis of mixed metal oxides for functional materials. One application of interest is the use of ALD for tailoring the refractive index of coatings. In homogeneously distributed composites of metal oxides, the refractive properties can be approximated as the average of the indices of the components. This is known as an effective medium approximation (EMA) and can be used to design the macroscopic properties of composites. ALD produces layered, anisotropic films, so the validity of an EMA in describing these films is not clear. Here, we use optical simulations and experimental characterization of stratified composites of TiO 2 and Al 2 O 3 to study the application of an EMA to the transmission and reflection behavior of ALD-prepared mixed metal oxide thin films. We found that when the characteristic layer thickness is smaller than roughly 10 nm, the optical spectra of theoretical and experimental ALD films match the equivalent theoretical spectra for a material with a refractive index calculated from a simple, compositionally weighted EMA. Hence, this EMA could describe the transmission and reflectance spectra in ALD-derived TiO2–Al 2 O 3 nanolaminates when the films were sufficiently stratified even without thorough characterization of the real layers in the films. As a result, we demonstrated that ALD can be used to prepare effectively homogenous mixed metal oxide films with a predictable, tailorable refractive index of any value between those of the TiO 2 and Al 2 O 3 component materials.

42 ENGINEERING↗

Synthesis, crystal and electronic structure of the Zintl phase Ba 16 Sb 11 . A case study uncovering greater structural complexity via monoclinic distortion of the tetragonal Ca 16 Sb 11 structure type

The binary Zintl phase Ba 16 Sb 11 has been synthesized and structurally characterized. Detailed studies via single-crystal X-ray diffraction methods indicate that although Ba 16 Sb 11 appears to crystallize in the tetragonal Ca 16 Sb 11 structure type (space group $P\bar{4}2_1$m with a=13.5647(9) Å, c=12.4124(12)Å, Z=2, R 1 = 3.14%; wR 2 = 4.77%), there exists an extensive structural disorder. Some Ba 16 Sb 11 crystals were found to be monoclinic and the structure was solved and refined in space group P2 1 (a=18.3929(12) Å, b=13.5233(8) Å, c=18.3978(12) Å, β=94.6600(10)°; Z=4, R 1 =5.84 %; wR 2 =9.58 %). The latter corresponds to a 2-fold superstructure of the tetragonal one, which provides a disorder-free structural model. In both descriptions, the disordered tetragonal and the ordered monoclinic superstructure, the basic building units that make up the structure of this Ba-rich compound are pairs of face-shared square antiprisms of Ba atoms, which are centered by Sb atoms. The dimerized antiprisms are linked into parallel chains via square prisms of Ba atoms, which are also centered by Sb atoms. The Zintl concept can be applied in a straightforward manner and as result, the structure of Ba 32 Sb 22 (=2×Ba 16 Sb 11 ) can be rationalized as (Ba 2+ ) 32 (Sb 3– ) 20 [Sb 2 ] 4– . Notably, the partitioning of the valence electrons is done taking into an account the homoatomic Sb–Sb contacts (d=3.01 Å), which can be clearly distinguished in the lower symmetry space group. Electronic structure calculations of Ba 16 Sb 11 are in good accordance with the Zintl rationalization and predict a semiconductor with a band gap of 0.77 eV.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Collaborative Research: Unravelling the Physics Associated with the Production of Extremely Dense Plasma States of Microscale (Final Report) Nanosecond-pulsed Discharges

The aim of this project is to study and establish the physical mechanisms that contribute to the formation of anomalously dense plasmas in high-pressure nanosecond-pulsed discharges. These discharges have a broad range of applications such as plasma-assisted combustion, plasma flow actuators, biomedical sterilization and exotic materials synthesis. The structure and formation of these discharges, producing high plasma densities of ~10 14 -10 15 cm -3 , are well-studied and understood. Fast-pulsed microscale high-pressure discharges can be driven to even higher densities of > 10 19 cm -3 , approaching warm dense matter conditions. The mechanisms that generate these plasmas have not been understood. Analysis of the warm dense matter state under laboratory conditions is an expensive and non-trivial endeavor. For instance, dense plasmas can be generated by electrical explosions of metal foils and wires. Plasmas generated after the explosion have a short lifetime and often present difficult conditions for diagnostics. Generation of dense plasmas was also achieved during high-voltage nanosecond pulsed discharges when the so-called explosive electron emission is obtained. Unfortunately, this process is very difficult to control for the studies of warm dense matter. In our recent study, we have shown that additional heating of plasma by lasers can further increase the density of plasma and even lead to the fully ionized state. This method, potentially, allows better control of the plasma parameters. In this work, we studied a second stage laser-heated micro-discharge using a self-consistent one-dimensional particle-in-cell Monte Carlo-collision (1D PIC-MCC) model coupled with Maxwell’s equations. We predicted the generation of a fully ionized plasma on the picosecond time scale. However, this model considered the plasma as an ideal gas despite the high pressure and the nearly fully ionized state. The ideal plasma model assumes that the dilute gas approximation is valid, where the inter-particle interactions are negligible. For charged particles this assumption holds as long as the shielded Coulomb potential assumption is valid. For very high plasma densities, this concept breaks down since the Debye sphere surrounding each charged particle no longer contains enough electrons to statistically provide the shielding of the single particle Coulomb interaction potential. At such densities, the plasma can no longer be described as ideal and non-ideal coupling effects need to be considered. In this report, we elucidate our recent work of developing a PIC-MCC model with improvements for non-ideal plasma conditions due to Coulomb coupling at high densities. In particular, we study the interaction of green light radiation and a dense microplasma, and explore the non-ideal plasma effects in this interaction. In this computational model, we implement the two most important non-ideal effects: ionization potential depression (IPD) and enhanced collision cross sections. Our primary goal is to study the physics associated with electromagnetic (EM) wave heating, also called the second-stage wave-heating, and establish the role of plasma non-ideality in this phenomenon. Our secondary goal is to improve the chemistry mechanism of the 1D PIC-MCC model by including a more detailed excited species collision treatment. At high pressures, stepwise ionization from excited species might play an important role in the ionization process. Previously, this ionization mechanism was neglected due to the excitation collision cross section of xenon being smaller than that of ionization. However, a preliminary study showed that the excited species density in the initial microplasma was an order of magnitude higher than the electron density. Therefore, my aim is to determine the significance of this additional ionization pathway to the plasma generation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A previously uncharacterized divisome-associated lipoprotein, DalA, is needed for normal cell division in Rhodobacterales

ABSTRACT The bacterial cell envelope is a key subcellular compartment with important roles in antibiotic resistance, nutrient acquisition, and cell morphology. We seek to gain a better understanding of proteins that contribute to the function of the cell envelope in Alphaproteobacteria . Using Rhodobacter sphaeroides , we show that a previously uncharacterized protein, RSP_1200, is an outer membrane (OM) lipoprotein that non-covalently binds peptidoglycan (PG). Using a fluorescently tagged version of this protein, we find that RSP_1200 undergoes a dynamic repositioning during the cell cycle and is enriched at the septum during cell division. We show that the position of RSP_1200 mirrors the location of FtsZ rings, leading us to propose that RSP_1200 is a newly identified component of the R. sphaeroides ’ divisome. Additional support for this hypothesis includes the co-precipitation of RSP_1200 with FtsZ, the Pal protein, and several predicted PG L,D-transpeptidases. We also find that a ∆ RSP_1200 mutation leads to defects in cell division, sensitivity to PG-active antibiotics, and results in the formation of OM protrusions at the septum during cell division. Based on these results, we propose to name RSP_1200 DalA (for division-associated lipoprotein A) and postulate that DalA serves as a scaffold to position or modulate the activity of PG transpeptidases that are needed to form envelope invaginations during cell division. We find that DalA homologs are present in members of the Rhodobacterales order within Alphaproteobacteria . Therefore, we propose that further analysis of this and related proteins will increase our understanding of the macromolecular machinery and proteins that participate in cell division in Gram-negative bacteria. IMPORTANCE Multi-protein complexes of the bacterial cell envelope orchestrate key processes like growth, division, biofilm formation, antimicrobial resistance, and production of valuable compounds. The subunits of these protein complexes are well studied in some bacteria, and differences in their composition and function are linked to variations in cell envelope composition, shape, and proliferation. However, some envelope protein complex subunits have no known homologs across the bacterial phylogeny. We find that Rhodobacter sphaeroides RSP_1200 is a newly identified lipoprotein (DalA) and that loss of this protein causes defects in cell division and changes the sensitivity to compounds, affecting cell envelope synthesis and function. We find that DalA forms a complex with proteins needed for cell division, binds the cell envelope polymer peptidoglycan, and colocalizes with enzymes involved in the assembly of this macromolecule. The analysis of DalA provides new information on the cell division machinery in this and possibly other Alphaproteobacteria.

59 BASIC BIOLOGICAL SCIENCES↗

Modeling MTS pyrolysis and SiC deposition kinetics using principal component analysis and neural networks

Accurate chemical kinetics modeling is crucial for improving the efficiency of chemical processing and synthesis of ceramic matrix composites. Detailed kinetic models are computationally expensive due to the large number of transported chemical species, while the simplified physics-based models, such as single-step global mechanisms, are efficient but often overlook key chemical intermediates and pathways. Recent deep learning approaches promise accurate and cost-effective models. Yet, they require additional closures for the transported nonlinear latent variables, complicating integration with existing solvers. In this work, we develop a hybrid linear—nonlinear reduced model for silicon carbide deposition from methyltrichlorosilane precursor by combining principal component analysis (PCA) and autoencoder (AE) neural network (NN) approaches. PCA is used to identify a smaller set of linear transport variables, enabling direct reuse of conventional transport solvers. NNs then reconstruct the full chemical state from these reduced variables. We demonstrate the method on a chemical vapor deposition reactor—comprising a gas-phase pyrolysis plug flow reactor and a heterogeneous surface reactor—over a wide range of temperatures, pressures, and residence times. Our PCA–AE model achieves high accuracy with only five transported scalars, achieving an eightfold cost reduction compared to detailed mechanisms, in both a priori (using data from the test set only) and a posteriori (coupled with a differential equation solver). In conclusion, notable errors arise primarily near training domain boundaries and for long residence times, indicating the need for domain shift indicators and better long-horizon predictions in future reduced chemistry model development.

autoencoder neural networks↗

Adaptive Discovery and Mixed-Variable Optimization of Next Generation Synthesizable Microelectronic Materials

Design of new microelectronic materials is characterized by several challenges such as high-dimensionality of the atomic structure-composition variable space, formidable cost of directly using high-fidelity simulations for design optimization, dispersity in literature-reported similar materials and synthesis methods, complex physical mechanisms, and mixed qualitative and quantitative design variables that lead to a disjointed design space. Even though machine learning (ML) techniques have been employed to expedite materials innovation, existing methods treat ML and design optimization as two separate processes, failing to resolve the fundamental challenges associated with high dimensionality and mixed-variable complexity. We have developed a ML enhanced mixed-variable material design optimization framework to efficiently extract useful information from existing data in literature and physics-based simulations to guide the autonomous search for optimal materials. Our proposed framework is composed of four computational modules: (1) a natural language processing (NLP) based virtual screening module, (2) classification based concept exploration module, (3) a density functional theory (DFT)-based high-fidelity evaluation model, and (4) a novel latent-variable Gaussian process (LVGP) ML model for mixed-variable problems with uncertainty quantification, which seamlessly integrates with Bayesian Optimization (BO) and achieves superb efficiency through embedded physics-based dimension reduction. Our approach is demonstrated and validated using the testbed of functional materials exhibiting metal-insulation transitions (MITs), with the targeted reversible resistivity changes (∼10^5) near room temperature. At the end of the 30-month project, we have developed a series of new ML techniques using NLP, conditional variational autoencoders, active learning, latent-variable Gaussian processes, integrated with Bayesian optimization. Our project has resulted in new predicted MITs compounds and improved understanding of MITs microscopic mechanisms, which in turn will revolutionize microelectronics science to provide energy-saving solutions. Our research has improved both creativity and efficiency in transforming rare-event discoveries of new functional materials to persistent innovations. In addition to open-sourcing the online MIT database and the classification model, the LVGP open source code has been downloaded more than 15,000 times within two years. More than 40 MIT compounds have been identified and many have been pursued experimentally via collaborators. The research results are published in close to 20 collaborative papers in high-impact journals, such as Chem. Mater., Appl. Phys. Rev., Sci. Rep., among others of design space.

36 MATERIALS SCIENCE↗

Cell wall pH and auxin transport velocity

According to the chemiosmotic polar diffusion hypothesis, auxin pulse velocity and basal secretion should increase with decreasing cell wall pH. Experiments were designed to test this prediction. Avena coleoptile sections were preincubated in either fusicoccin (FC), cycloheximide, pH 4.0, or pH 8.0 buffer and subsequently their polar transport capacities were determined. Relative to controls, FC enhanced auxin (IAA) uptake while CHI and pH 8.0 buffer reduced IAA uptake. Nevertheless, FC reduced IAA pulse velocity while cycloheximide increased velocity. Additional experiments showed that delivery of auxin to receivers is enhanced by increased receiver pH. This phenomenon was overcome by a pretreatment of the tissue with IAA. Our data suggest that while acidic wall pH values facilitate cellular IAA uptake, they do not enhance pulse velocity or basal secretion. These findings are inconsistent with the chemiosmotic hypothesis for auxin transport.

Plant Growth Regulators/metabolism↗

Quasicubic model for metal halide perovskite nanocrystals

In this work, we present an analysis of quantum confinement of carriers and excitons, and exciton fine structure, in metal halide perovskite (MHP) nanocrystals (NCs). Starting with coupled-band k · P theory, we derive a nonparabolic effective mass model for the exciton energies in MHP NCs valid for the full size range from the strong to the weak confinement limits. We illustrate the application of the model to CsPbBr 3 NCs and compare the theory against published absorption data, finding excellent agreement. We then apply the theory of electron-hole exchange, including both short- and long-range exchange interactions, to develop a model for the exciton fine structure. We develop an analytical quasicubic model for the effect of tetragonal and orthorhombic lattice distortions on the exchange-related exciton fine structure in CsPbBr 3 , as well as some hybrid organic MHPs of recent interest, including formamidinium lead bromide (FAPbBr 3 ) and methylammonium lead iodide (MAPbI3). Testing the predictions of the quasicubic model using hybrid density functional theory (DFT) calculations, we find qualitative agreement in tetragonal MHPs but significant disagreement in the orthorhombic modifications. Moreover, the quasicubic model fails to correctly describe the exciton oscillator strength and with it the long-range exchange corrections in these systems. Introducing the effect of NC shape anisotropy and possible Rashba terms into the model, we illustrate the calculation of the exciton fine structure in CsPbBr 3 NCs based on the results of the DFT calculations and examine the effect of Rashba terms and shape anisotropy on the calculated fine structure.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Unexpected features in the optical vibrational spectra of δ-UO 3

Uranium trioxide displays a complex chemical phase space, with at least six structurally distinct polymorphs accessible via different synthetic routes. Remarkably, despite its technological importance, full structural and electronic characterization of these polymorphs remains an open area of study. δ-UO 3 in particular has attracted significant theoretical attention due to its high point group and space group symmetries, having U (VI) in octahedral coordination with polyhedra interconnected through corner-sharing to build a 3-D cubic lattice with space group symmetry Pm-3m and Z = 1. Critical experimental information, such as its optical vibrational spectra, are not known. Here, we study the Raman and infrared (IR) spectra of δ-UO 3 together with the support of density functional theory (DFT) calculations for spectral interpretation. A symmetry analysis of the DFT-predicted phonon eigenmodes indicates that δ-UO 3 should have two IR active modes and no Raman active modes. Experimental results, however, indicate significant Raman scattering from δ-UO 3 . We therefore propose four potential explanations for this apparent contradiction: a possible tetragonal distortion to the cubic cell, the existence of a surface impurity layer, vacancy scattering, and structural activation of Raman signal. We use powder X-ray diffraction and confocal Raman spectroscopy with depth profiling to investigate these possibilities and suggest future experiments to explore this phenomenon in more detail. Understanding the lattice dynamics of δ-UO 3 is important for identification of technogenic U phases via Raman and infrared spectroscopy and our results indicate that the simple understanding of δ-UO 3 as a high-symmetry cubic structure should be reconsidered.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗

Retrospective Analysis of Backwater Habitat Availability Using Remote Sensing

Low-velocity channel-margin habitats serve as important nursery habitats for the endangered Colorado pikeminnow (Ptychocheilus lucius) in the middle Green River between Jensen and Ouray, Utah. These habitats, known as backwaters, are associated with emergent sandbars, and are shaped and reformed annually by peak flows. Our recent knowledge about backwater characteristics and dynamics that is summarized in the synthesis report (Grippo et al. 2015) was based on detailed annual survey data from a relatively small sample of backwaters that were collected from 2003 to 2014 and reach-wide evaluations of backwater surface area based on manual interpretation of aerial and satellite imagery. Methods that bridge the gap between the detailed surveys from a small number of backwaters and the reach-wide assessment of their surface area would enable an assessment of the availability of backwater habitats that meet the minimum depth requirements for suitable habitat for Colorado pikeminnow. In 2015 Argonne National Laboratory (Argonne) tested three regression models—linear, multiple, and partial least square (PLS) regression models—for estimating backwater depth using National Agriculture Imagery Program (NAIP) imagery collected in July 2006 that covered the Jensen-Ouray reach of the Green River (Hamada and LaGory 2016). The results suggested that a PLS regression model showed high correlation with the reference depth (R 2 = 0.69) and had the most unbiased and consistent estimate of backwater depth. The results also indicated that the PLS model would provide reasonable estimates of the amount of habitat providing a minimum suitable depth of 30 cm for young-of-the-year Colorado pikeminnow, even though absolute depth estimates may be uncertain for backwater areas deeper than approximately 40 cm. The study also provided insights regarding the amount and the selection of calibration and validation (cal-val) data needed for improving the accuracy of depth prediction.

54 ENVIRONMENTAL SCIENCES↗

Experimental Calibration of Flow Rule Coefficients for LS-Dyna Mat_213

Composite materials have become a popular choice in aerospace engineering applications due to their low weight, high specific stiffness, high specific energy absorption, and economic synthesis. As use of composite materials has proliferated, the need for accurate simulation under extreme loading conditions is of increasing importance. *MAT_213 is a tabulated plasticity and damage material model in LS-DYNA designed for such use cases. To account for the nonlinear material response, a plasticity algorithm is used, which employs a non-associative flow rule with coefficients requiring experimental calibration. This work outlines an experimental methodology for calibrating these coefficients using the macroscopic stress-strain response and microscopic in-plane strain measurements from coupon scale experiments. The model material used for the development and demonstration of this methodology is a plain weave woven fabric with T300 carbon fiber tows in the warp direction and Kevlar© 49 fiber tows in the weft direction, and an Epon© 828/Epikure© 3223 thermoset epoxy matrix. Coupon-scale tensile testing performed at multiple fiber angles (carbon fibers oriented 0º, 22.5º, 45º, 67.5º, and 90º relative to the loading axis) provides data for multiple calibration processes. The measured plastic strain increments in the principal material coordinate are used to set limits for ratios between flow rule coefficients. The identified ranges are used to bound the search for a set of coefficients, leading to convergence of different coupon responses in the effective stress – effective plastic strain space. Flow rule coefficients found by this methodology are then used as (i) input for LS-DYNA simulations of uniaxial tension coupon tests at varying fiber angles, and (ii) comparison of predicted values for the plastic Poisson’s ratio in global coordinates with the experimentally measured one for all off-axis orientations. These results are quantitatively and qualitatively compared to macroscopic stress-strain curves and experimental strain components. The conformity of the experimentally measured quantities with those derived from the simulations provide a means for assessing the quality of calibrated coefficients.

Polymer Matrix Composites↗

The γ -process nucleosynthesis in core-collapse supernovae

The γ-process nucleosynthesis in core-collapse supernovae is generally accepted as a feasible process for the synthesis of neutron-deficient isotopes beyond iron. However, crucial discrepancies between theory and observations still exist: the average yields of γ-process nucleosynthesis from massive stars are still insufficient to reproduce the solar distribution in galactic chemical evolution calculations, and the yields of the Mo and Ru isotopes are a factor of ten lower than the yields of the other γ-process nuclei. We investigate the γ-process in five sets of core-collapse supernova models published in the literature with initial masses of 15, 20, and 25 M ⊙ at solar metallicity. We compared the γ-process overproduction factors from the different models. To highlight the possible effect of nuclear physics input, we also considered 23 ratios of two isotopes close to each other in mass relative to their solar values. Further, we investigated the contribution of C–O shell mergers in the supernova progenitors as an additional site of the γ-process. Our analysis shows that a large scatter among the different models exists for both the γ-process integrated yields and the isotopic ratios. We find only ten ratios that agree with their solar values, all the others differ by at least a factor of three from the solar values in all the considered sets of models. The γ-process within C–O shell mergers mostly influences the isotopic ratios that involve intermediate and heavy proton-rich isotopes with A > 100. We conclude that there are large discrepancies both among the different data sets and between the model predictions and the solar abundance distribution. More calculations are needed; particularly updating the nuclear network, because the majority of the models considered in this work do not use the latest reaction rates for the γ-process nucleosynthesis. Moreover, the role of C–O shell mergers requires further investigation.

79 ASTRONOMY AND ASTROPHYSICS↗

Experimental study of the partitioning of some platinum group elements (Pd and Ir) between orthopyroxene and silicate melt

Past experiments and observations on natural samples have largely focused on the roles of olivine and chromite in controlling the behaviour of the platinum-group elements (PGE) during melting and solidification, whereas other phases, such as pyroxene, have gone largely uncharacterized. Here, to address this, experiments have been done to measure the partitioning of Pd (with a subset of results for Ir), between orthopyroxene and silicate melt at 1340 °C, 0.1 MPa and log fO 2 of FMQ - 1 to FMQ + 6 (FMQ = Fayalite-Magnetite-Quartz). The X-ray Absorption Near-Edge Structure (XANES) was measured in a subset of experiment glasses. Glass concentrations of Pd (corrected to unit Pd activity) increase from ~6 to ~650 ug/g over the fO 2 range of experiments. The slope of the solubility-fO 2 relation is consistent with Pd 1+ as the dominant oxidation state, with evidence for Pd 0 and Pd 2+ at the lowest and highest experiment fO 2 , respectively. Consistent with this result, the XANES reveal spectral features similar to Pd 0 and Pd 2+ spectral reference materials (specRM) at the most reduced and oxidized synthesis conditions, respectively. Other lines of evidence require the presence of a third melt species, here interpreted to be Pd 1+ . Values of orthopyroxene/melt partition coefficients for Pd (D Pd Opx/melt ) are 0.0051 (+/-0.006) at log fO 2 < ΔFMQ + 3, increasing with fO 2 to a maximum of 0.013 at ~FMQ + 6. Sodium partition coefficients, expected to be similar to Pd, range from 0.0061 (+/-0.00061) at FMQ + 3, increase to 0.007–0.009 at higher fO 2 , but with no clear systematic trend. A value for D Ir Opx-melt of ~0.6 was measured at ~FMQ + 4, indicating significantly more compatible behaviour for Ir relative to Pd. Partitioning results are interpreted in the context of the Blundy-Wood elastic strain model in which the variation in partitioning is related to ionic radius mismatch to an optimal crystallographic site size. Based on the trend in ionic radius with oxidation state, the estimated ionic radius of Pd 1+ in octahedral coordination is similar to Na 1+ , and comparison to previous orthopyroxene-melt partitioning experiments suggests D Pd1+ opx/melt and DNaopx/melt should be nearly identical, consistent with the results of this study. The ionic radius of VI-fold Pd 2+ is close to the optimal M2 site size, so an increased proportion of this species with fO 2 accounts for the larger values of D Pd opx/melt at the highest fO 2 investigated. The much larger partition coefficient for Ir is consistent with the presence of Ir 2+ , whose estimated ionic radius is close to Fe 2+ and Mg 2+ , as well as predictions for the optimal M1 site size. With the assumption that D Pd opx/melt = D Na opx/melt , combined with a revised value for the Pd content of the primitive mantle, a melting model is presented that better reproduces the Pd concentration of high degree melts from sulfide-free mantle sources.

58 GEOSCIENCES↗

Deep VLA imaging of twelve extended 3CR quasars

A representative sample of 12 extended quasars from the 3CR catalog has been imaged at 4.9 GHz using the Very Large Array (VLA). These full synthesis observations typically achieve an rms noise of 20 micro-Jy per beam. Jets are detected on at least one side of every source. The jets are well collimated compared with those in less powerfull sources, but spreading is detected in most of them. The opening angles of several jets are not constant, but show recollimation after an initial regime of rapid spreading. Many of the jets contain quasiperiodic strings of knots, of which the knot closest to the central feature is usually the brightest (until the jet nears its hot spot). The degrees of linear polarization at the jet knots range from less than 5% to approx. 50%, but show no common trend with distance along the jets. In knots that are elongated in directions close to that of the jet, the E vectors tend to be orthogonal to the jet axis. The prominence of the inner, straighter jet segments relative to the extended lobes correlates significantly with the prominence of the milliarcsecond-scale central features, but the prominence of the more bent jet segments does not. Candidates for counterjet emission are detected in seven sources, but there is no unambiguous, continuous counterjet in any of them. Estimates of the flux density ratios between the straighter jet segments and the counterjets based on these tentative detections range from 1.2:1 to greater than 175:1. There is no evidence in this sample that counterjet detectability correlates with such putative inclination indicators as central feature prominence or projected linear size. There is also no evidence that the prominence of the counterjets anticorrelates with that of the jets as predicted by simple relativistic-beaming models for the jet/counterjet asymmetry. There is, however, strong evidence that large bends in the main jet favor counterjet detection, and there are no counterjet candidates opposite long, uninterrupted straight segments of the main jets. The detectability of the counterjets in these quasars may therefore be strongly influenced by interactions between the underlying beams and inhomogeneities in the surrounding material. We offer a new empirical definition of the term 'hot spot' that is intended to improve the distinction between such features and 'jet knots.' Both the compactness of hot spots and their position in the lobe are affected by whether they are fed by a detectable jet. When the hot spots differ significantly in compactness, the more compact one is always on the jetted side. Jetted hot spots are also more likely to be recessed deeply from the outer edge of their lobes than are their counterjetted counterparts. The jetted hot spot is less prominent relative to other extended emission if the jet bends through a large angle, particularly if a large bend occurs abruptly. The counterjetted hot spot is also less well defined if the jet is more bent. The lobes of several sources show considerable inhomogeneity, including filamentation.

Bridle, Alan H.↗

Coupled Microbial-Conversion and Computational-Fluid-Dynamics (CFD) Models for Butanediol Production in Micro-Aerated Reactors

Microbial conversion of substrates to macromolecules has been widely used in the synthesis of value-added products in pharmaceutical and biotechnology industries. These bioreactions are also being investigated in the production of low-value commodities such as biofuels [1]. Gas-liquid mass-transfer and transport-reaction coupling are important challenges when designing and scaling up these reactor systems. Experiments in wellmixed small-scale reactors have enabled characterization of microbial reactivity, while their coupling with macroscale transport remains relatively unexplored. In this work, we use a coupled metabolic-CFD model to study the action of a genetically engineered microbe Zymomonas mobilis [2] on sugars to produce 2,3-Butanediol (BDO). BDO is an important hydrocarbon intermediate that can be catalytically upgraded to several fuels and chemicals [3]. An important aspect to this particular microbial conversion is the need for micro-aerated environments as opposed to traditional aerobic fermentation. Slight variations in oxygen concentration can result in competing reaction pathways that disable BDO production. Hence, gas-liquid mass transfer and transport need to be optimized in large-scale reactors to maximize BDO production, for which CFD is a valuable tool. The aerobic-fermentation CFD model previously developed by the authors [5] for simulating bubble-column and airlift reactors at scale was used in this study. The Reynolds-averaged mass, momentum and energy transport equations for interpenetrating gas and liquid phase are solved in this model along with the transport and interphase mass transfer of oxygen. Our previous work used a phenomenological model for microbial oxygen uptake that neglected microbial growth and other reaction pathways. In this work, a detailed metabolic model enabled prediction of product formation and inhibition pathways. In order to manage computational cost, we used a subcycling technique [6] that takes advantage of the clear separation in transport (~ 200 sec) and reaction (~ 2-3 hours) timescales. The CFD model is first solved to steady state, after which the metabolic model is advanced at every cell in the computational domain using the local oxygen concentration. The CFD model is then run to achieve a new steady state that provides a new oxygen distribution for the metabolic model. This process, where reaction and fluid updates are interleaved together, is iterated until reactants are completely exhausted. This work will examine the performance of different reactor designs such as bubble column and airlift reactors at scale (250-500 m3). Oxygen mass-transfer coefficient and distribution are critically analyzed among reactors, and optimization studies pertaining to aeration is presented. Furthermore, it has been observed in experiments that high BDO production may be achieved by manipulating the aerobic environment over the course of reaction, such that oxygen concentration is high during the growth phase, and very low as sugar is depleted. This characteristic will be addressed by our simulations for which a timedependent scheduling strategy for aeration is presented that maximizes BDO production. [1] Humbird, D., Davis, R., and McMillan, J., Aeration costs in stirred-tank and bubble column bioreactors, Biochemical Engineering Journal, 127, 161—166, 2017 [2] Yang, S., Mohagheghi, A., Franden, M. A., Chou, Y.-C., Chen, X., Dowe, N., Himmel, M. E., and Zhang, M., Metabolic engineering of zymomonas mobilis for 2, 3-butanediol production from lignocellulosic biomass sugars. Biotechnology for biofuels, 9(1):189, 2016 [3] Kim, S. J., Sim, H. J., Kim, J. W., Lee, Y. G., Park, Y. C., and Seo, J. H., Enhanced production of 2,3-butanediol from xylose by combinatorial engineering of xylose metabolic pathway and cofactor regeneration in pyruvate decarboxylase-deficient Saccharomyces cerevisiae. Bioresource Technology, 245:1551–1557, 2017 [4] Weller, H., Tabor, G., Jasak, H. and Fureby, C., A tensorial approach to computational continuum mechanics using object-oriented techniques, Computers in physics, 12, 6, 620--631, 1998

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

An Operational Overview of the EXport Processes in the Ocean from RemoTe Sensing (EXPORTS) Northeast Pacific Field Deployment

The goal of the EXport Processes in the Ocean from RemoTe Sensing (EXPORTS) field campaign is to develop a predictive understanding of the export, fate, and carbon cycle impacts of global ocean net primary production. To accomplish this goal, observations of export flux pathways, plankton community composition, food web processes, and optical, physical, and biogeochemical (BGC) properties are needed over a range of ecosystem states. Here we introduce the first EXPORTS field deployment to Ocean Station Papa in the Northeast Pacific Ocean during summer of 2018, providing context for other papers in this special collection. The experiment was conducted with two ships: a Process Ship, focused on ecological rates, BGC fluxes, temporal changes in food web, and BGC and optical properties, that followed an instrumented Lagrangian float; and a Survey Ship that sampled BGC and optical properties in spatial patterns around the Process Ship. An array of autonomous underwater assets provided measurements over a range of spatial and temporal scales, and partnering programs and remote sensing observations provided additional observational context. The oceanographic setting was typical of late-summer conditions at Ocean Station Papa: a shallow mixed layer, strong vertical and weak horizontal gradients in hydrographic properties, sluggish sub-inertial currents, elevated macronutrient concentrations and low phytoplankton abundances. Although nutrient concentrations were consistent with previous observations, mixed layer chlorophyll was lower than typically observed, resulting in a deeper euphotic zone. Analyses of surface layer temperature and salinity found three distinct surface water types, allowing for diagnosis of whether observed changes were spatial or temporal.The 2018 EXPORTS field deployment is among the most comprehensive biological pump studies ever conducted. A second deployment to the North Atlantic Ocean occurred in spring 2021, which will be followed by focused work on data synthesis and modeling using the entire EXPORTS data set.

Biological pump↗

Methane fluxes in tidal marshes of the conterminous United States

Abstract Methane (CH 4 ) is a potent greenhouse gas (GHG) with atmospheric concentrations that have nearly tripled since pre‐industrial times. Wetlands account for a large share of global CH 4 emissions, yet the magnitude and factors controlling CH 4 fluxes in tidal wetlands remain uncertain. We synthesized CH 4 flux data from 100 chamber and 9 eddy covariance (EC) sites across tidal marshes in the conterminous United States to assess controlling factors and improve predictions of CH 4 emissions. This effort included creating an open‐source database of chamber‐based GHG fluxes ( https://doi.org/10.25573/serc.14227085 ). Annual fluxes across chamber and EC sites averaged 26 ± 53 g CH 4 m −2 year −1 , with a median of 3.9 g CH 4 m −2 year −1 , and only 25% of sites exceeding 18 g CH 4 m −2 year −1 . The highest fluxes were observed at fresh‐oligohaline sites with daily maximum temperature normals (MATmax) above 25.6°C. These were followed by frequently inundated low and mid‐fresh‐oligohaline marshes with MATmax ≤25.6°C, and mesohaline sites with MATmax >19°C. Quantile regressions of paired chamber CH 4 flux and porewater biogeochemistry revealed that the 90th percentile of fluxes fell below 5 ± 3 nmol m −2 s −1 at sulfate concentrations >4.7 ± 0.6 mM, porewater salinity >21 ± 2 psu, or surface water salinity >15 ± 3 psu. Across sites, salinity was the dominant predictor of annual CH 4 fluxes, while within sites, temperature, gross primary productivity (GPP), and tidal height controlled variability at diel and seasonal scales. At the diel scale, GPP preceded temperature in importance for predicting CH 4 flux changes, while the opposite was observed at the seasonal scale. Water levels influenced the timing and pathway of diel CH 4 fluxes, with pulsed releases of stored CH 4 at low to rising tide. This study provides data and methods to improve tidal marsh CH 4 emission estimates, support blue carbon assessments, and refine national and global GHG inventories.

54 ENVIRONMENTAL SCIENCES↗

PPI DataHub Project Data Package: S. elongatus PCC 7942 Circadian Control Bioproduction Transcriptomics (PB-DP3)

The purpose of this experiment was to evaluate how circadian clock regulation impacts carbon partitioning between storage, growth, and product synthesis in Synechococcus elongatus PCC 7942 in providing insights to strategies for enhanced bioproduction. Sample data was acquired using a Illumina HiSeq sequencer system and processed for RNA sequencing (RNA-Seq) expression analysis. Transcriptomic differential expression analysis revealed coordinated circadian clock-driven adjustment of the cell cycle and rewiring of energy and carbon metabolism. Processed RNA-Seq datasets are openly accessible from the PNNL DataHub project dataset download page and contain secondary processed RNA-seq results files and supporting metadata materials linked to relevant source code information supporting data transparency and reuse.

59 BASIC BIOLOGICAL SCIENCES↗