Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “model compounds”

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
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Augmenting machine learning of energy landscapes with local structural information

We present a machine learning approach for accurately predicting formation energies of binary compounds in the context of crystal structure predictions. The success of any machine learning model depends significantly on the choice of representation used to encode the relevant physical information into machine-learnable data. We test different representation schemes based on partial radial and angular distribution functions (RDF+ADF) on Al–Ni and Cd–Te structures generated using our genetic algorithm for structure prediction. We observe a remarkable improvement in predictive accuracy upon transitioning from global to atom-centered representations, resulting in a threefold decrease in prediction errors. We show that a support vector regression model using a combination of atomic radial and angular distribution functions performs best at the formation energy prediction task, providing small root mean squared errors of 3.9 meV/atom and 10.9 meV/atom for Al–Ni and Cd–Te, respectively. We test the performance of our models against common traditional descriptors and find that RDF- and ADF-based representations significantly outperform many of those in the prediction of formation energies. The high accuracy of predictions makes our machine learning models great candidates for the exploration of energy landscapes.

Honrao, Shreyas J. (ORCID:0000000166155257)↗
DOE OSTI2021

Four-decade follow-up of a plutonium-contaminated puncture wound treated with Ca-DTPA

Contaminated wounds are a common route of internal deposition of radionuclides for nuclear and radiation workers. They may result in significant doses to radiosensitive organs and tissues in an exposed individual's body. The United States Transuranium and Uranium Registries' whole-body donor (Case 0303) accidentally punctured his finger on equipment contaminated with plutonium nitrate. The wound was surgically excised and medically treated with intravenous injections of Ca-DTPA. A total of 16 g Ca-DTPA was administered in 18 treatments during the 2 months following the accident. Ninety-three urine samples were collected and analysed over 14 years following the accident. An estimated 239 Pu activity of 73.7 Bq was excreted during Ca-DTPA treatment. Post-mortem radiochemical analysis of autopsy tissues indicated that 40 years post-accident 21.6 ± 0.2 Bq of 239 Pu was retained in the skeleton, 12.2 ± 0.3 Bq in the liver, and 3.7 ± 0.1 Bq in other soft tissues; 1.35 ± 0.02 Bq of 239 Pu was measured in tissue samples from the wound site. To estimate the plutonium intake, late urine measurements, which were unaffected by chelation, and post-mortem radiochemical analysis results were evaluated using the IMBA Professional Plus software. The application of the National Council on Radiation Protection and Measurements wound model with an assumption of intake material as a predominantly strongly retained soluble plutonium compound with a small insoluble fraction adequately described the data (p = 0.46). The effective intake was estimated to be 50.2 Bq of plutonium nitrate and 1.5 Bq of the fragment. Here, the prompt medical intervention with contaminated tissue excision and subsequent Ca-DTPA decorporation therapy reduced 239 Pu activity available for uptake and long-term retention in this individual's systemic organs by a factor of 38.

239Pu↗
DOE OSTI2024

Inhibition of atomic layer deposition of TiO 2 by functionalizing silicon surface with 4-fluorophenylboronic acid

As the size of the components in electronic devices decreases, new approaches and chemical modification schemes are needed to produce nanometer-size features with bottom-up manufacturing. Organic monolayers can be used as effective resists to block the growth of materials on non-growth substrates in area-selective deposition methods. However, choosing the appropriate surface modification requires knowledge of the corresponding chemistry and also a detailed investigation of the behavior of the functionalized surface in realistic deposition schemes. This study aims to investigate the chemistry of boronic acids that can be used to prepare such non-growth areas on elemental semiconductors. 4-Fluorophenylboronic acid is used as a model to investigate the possibility to utilize the Si(100) surface functionalized with this compound as a non-growth substrate in a titanium dioxide (TiO 2 ) deposition scheme based on sequential doses of tetrakis(dimethylamido)titanium and water. Here, a combination of X-ray photoelectron spectroscopy and time-of-flight secondary ion mass spectrometry allows for a better understanding of the process. The resulting surface is shown to be an effective non-growth area to TiO 2 deposition when compared to currently used H-terminated silicon surfaces but to exhibit much higher stability in ambient conditions.

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

When does global attention help: a unified empirical study on atomistic graph learning

Graph neural networks (GNNs) are widely used as surrogates for costly experiments and first-principles simulations to study the behavior of compounds at atomistic scale, and their architectural complexity is constantly increasing to enable the modeling of complex physics. While most recent GNNs combine more traditional message passing neural networks (MPNNs) layers to model short-range interactions with more advanced graph transformers (GTs) with global attention mechanisms to model long-range interactions, it is still unclear when global attention mechanisms provide real benefits over well-tuned MPNN layers due to inconsistent implementations, features, or hyperparameter tuning. We introduce the first unified, reproducible benchmarking framework–built on HydraGNN–that enables seamless switching among four controlled model classes: MPNN, MPNN with chemistry/topology encoders, GPS-style hybrids of MPNN with global attention, and fully fused localglobal models with encoders. Using seven diverse open-source datasets for benchmarking across regression and classification tasks, we systematically isolate the contributions of message passing, global attention, and encoder-based feature augmentation. Our study shows that encoder-augmented MPNNs form a robust baseline, while fused localglobal models yield the clearest benefits for properties governed by long-range interaction effects. We further quantify the accuracycompute trade-offs of attention, reporting its overhead in memory. Together, these results establish the first controlled evaluation of global attention in atomistic graph learning and provide a reproducible testbed for future model development.

Equivariant graph neural networks↗
DOE OSTI2023

Global seamless tidal simulation using a 3D unstructured-grid model (SCHISM v5.10.0)

We present a new 3D unstructured-grid global ocean model to study both tidal and nontidal processes, with a focus on the total water elevation. Unlike existing global ocean models, the new model resolves estuaries and rivers down to ~8 m without the need for grid nesting. The model is validated with both satellite and in situ observations for elevation, temperature, and salinity. Tidal elevation solutions have a mean complex root-mean-square error (RMSE) of 4.2 cm for M2 and 5.4 cm for all five major constituents in the deep ocean. The RMSEs for the other four constituents, S2, N2, K1, and O1, are, respectively, 2.05, 0.93, 2.08, and 1.34 cm). The nontidal residual assessed by a tide gauge dataset (GESLA) has a mean RMSE of 7 cm. For the first time ever, we demonstrate the potential for seamless simulation on a single mesh from the global ocean into several estuaries along the US West Coast. The model is able to accurately capture the total elevation, even at some upstream stations. The model can therefore potentially serve as the backbone of a global tide surge and compound flooding forecasting framework.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Accurate determination of solvation free energies of neutral organic compounds from first principles

The main goal of molecular simulation is to accurately predict experimental observables of molecular systems. Another long-standing goal is to devise models for arbitrary neutral organic molecules with little or no reliance on experimental data. While separately these goals have been met to various degrees, for an arbitrary system of molecules they have not been achieved simultaneously. For biophysical ensembles that exist at room temperature and pressure, and where the entropic contributions are on par with interaction strengths, it is the free energies that are both most important and most difficult to predict. We compute the free energies of solvation for a diverse set of neutral organic compounds using a polarizable force field fitted entirely to ab initio calculations. The mean absolute errors (MAE) of hydration, cyclohexane solvation, and corresponding partition coefficients are 0.2 kcal/mol, 0.3 kcal/mol and 0.22 log units, i.e. within chemical accuracy. The model (ARROW FF) is multipolar, polarizable, and its accompanying simulation stack includes nuclear quantum effects (NQE). The simulation tools’ computational efficiency is on a par with current state-of-the-art packages. The construction of a wide-coverage molecular modelling toolset from first principles, together with its excellent predictive ability in the liquid phase is a major advance in biomolecular simulation.

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

Estimating source-sink distributions and fluxes of reactive nitrogen and sulfur within a mixed forest canopy

The vertical source-sink distribution of air pollutants within and above forested canopies is necessary for describing the biological, physical, and chemical processes influencing the soil-vegetation-atmosphere exchange. Here, this study implemented inverse modeling methods to estimate the source-sink and flux profiles of reactive nitrogen (N) and sulfur (S) compounds from measurements of the mean concentration profiles of ammonia (NH 3 ), nitric acid (HNO 3 ), sulfur dioxide (SO 2 ), and particulate ammonium (NH 4 + ), nitrate (NO 3 − ), and sulfate (SO 4 2− ) at a forest site in the southern Appalachian Mountains. Three inverse approaches utilizing different approximations to scalar transport within the canopy were developed and evaluated against sensible heat flux measurements. The Eulerian model (EUL), which incorporates vertical velocity skewness, performed well in reproducing the turbulent heat fluxes and was subsequently used to calculate the chemical source-sink and flux profiles. Above-canopy fluxes of NH 3 were downward, indicating that the forest was a net sink of NH 3 . The soil/litter layer was both a source and a sink for NH 3 but the exchange rate at the forest floor was small. Fluxes of HNO 3 , SO 2 , NO 3 − , NH 4 + , and SO 4 2- were uni-directional (deposition only) between the air and the canopy/ground and increased monotonically from the forest floor to the canopy top. Crown foliage dominated the uptake of reactive N and S during the growing season, accounting for 80–90% of the total canopy-scale flux. Fluxes and canopy-ground partitioning estimated using the resistance-based Surface Tiled Aerosol and Gas Exchange (STAGE) model were generally comparable to EUL. The comparison highlights the need for improved parameterizations of litter exchange and NH 3 compensation points in resistance models for forest ecosystems. The findings here benefit the application of critical loads in forest ecosystems and guide further development of resistance-based exchange models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

IMPECCABLE: Integrated Modeling Pipeline for COVID Cure by Assessing Better Leads

ABSTRACT The drug discovery process currently employed in the pharmaceutical industry typically requires about 10 years and $2-3 billion to deliver one new drug. This is both too expensive and too slow, especially in emergencies like the COVID-19 pandemic. In silico methodologies need to be improved both to select better lead compounds, so as to improve the efficiency of later stages in the drug discovery protocol, and to identify those lead compounds more quickly. No known methodological approach can deliver this combination of higher quality and speed. Here, we describe an Integrated Modeling PipEline for COVID Cure by Assessing Better LEads (IMPECCABLE) that employs multiple methodological innovations to overcome this fundamental limitation. We also describe the computational framework that we have developed to support these innovations at scale, and characterize the performance of this framework in terms of throughput, peak performance, and scientific results. We show that individual workflow components deliver 100× to 1000× improvement over traditional methods, and that the integration of methods, supported by scalable infrastructure, speeds up drug discovery by orders of magnitudes. IMPECCABLE has screened ~10 11 ligands and has been used to discover a promising drug candidate. These capabilities have been used by the US DOE National Virtual Biotechnology Laboratory and the EU Centre of Excellence in Computational Biomedicine.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Exploring the Links between Structural Distortions, Orbital Ordering, and Multipolar Magnetic Ordering in Double Perovskites Containing Re(VI) and Os(VII)

A combination of high-resolution powder diffraction techniques and solid-state NMR has been employed to explore the links between crystal structure, orbital ordering, and magnetism in three isostructural double perovskites containing transition metal ions with a 5d 1 configuration. In Ba 2 ZnReO 6 , both neutron and synchrotron X-ray powder diffraction data reveal a cubic-to-tetragonal transition at 23 K that breaks the degeneracy of the t 2g orbitals and leads to a pattern of orbital ordering that stabilizes magnetic ordering when the sample is cooled below 16 K. Similar behavior is observed in Ba 2 MgReO 6 , with an orbital ordering temperature of 33 K and a magnetic ordering temperature of 18 K. Prior theoretical works suggest that the pattern of orbital order seen in the P42/mnm space group is needed to stabilize the heavily canted antiferromagnetism of these compounds. Unfortunately, powder diffraction data is not sensitive enough to differentiate between the I4/mmm and P42/mnm structural models, as the distortions are too subtle to be unambiguously identified from either neutron or synchrotron X-ray powder diffraction methods. In contrast, both diffraction and 7 Li NMR data indicate that Ba 2 LiOsO 6 retains the cubic structure down to 1.7 K. The antiferromagnetic ground state and lack of any sign of orbital ordering in Ba 2 LiOsO 6 provide compelling evidence that the electronically driven tetragonal distortion seen in Ba2ZnReO6, and Ba2MgReO6 is intimately linked to the magnetic ordering seen in those compounds. The absence of magnetic reflections in high intensity neutron powder diffraction data collected on Ba 2 MgReO 6 strongly suggests ordering of multipolar moments on Re(VI), likely ferro-octupolar ordering.

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

Monte Carlo Hauser-Feshbach computer code system to model nuclear reactions: YAHFC

A computer program framework, YAHFC, to model low-energy nuclear reactions is presented. The framework allows for reactions with incident particles ranging from protons/neutrons to alphas and is designed to address reactions that ultimately lead to the formation of compound nuclear systems that then decay statistically as outlined in concepts of Hauser and Feshbach. Additionally, instead of a reaction, it is also possible to model the decay of a nuclear system with an initial excitation and population. The code models nuclear decays with a Monte Carlo process that tracks the decay of each state. This allows for an exact representation of the spectra for all emitted particles in each of the final exit channels and the possibility of generating reaction data for simulation purposes. The program is interfaced with the optical model code system FRESCOX to calculate transmission coefficients as well as the effects of coupled channels and other direct excitations via the distorted wave Born approximation (DWBA). Modules are included to account for nuclear processes such as width corrections, pre-equilibrium emission, and fission. The program is controlled by a series of input commands and while a set of input parameters exists for each projectile and target, the input commands allow for complete control over each input parameter. Extensive data files are produced and a program is provided that converts YAHFC data files into nuclear data library entries in the generalized nuclear data structure (GNDS).

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Machine Learning Models for Predicting Molecular UV–Vis Spectra with Quantum Mechanical Properties

Accurate understanding of Ultraviolet–visible (UV–Vis) spectra is critical for highthroughput design of compounds for drug discovery. Experimentally determining UV–Vis spectra can become expensive when dealing with a large quantity of novel molecules. This provides us an opportunity to drive computational advances in molecular property predictions using quantum mechanics and machine learning. In this work, we use both Quantum Mechanically (QM) predicted and measured UV–Vis spectra as input to modify four different machine learning architectures: UVvis-SchNet, UVvis- DTNN, UVvis-Transformer, and UVvis-MPNN. Here we find that the UVvis-MPNN model outperforms the other models when using optimized 3D coordinates and QM predicted spectra as input features. This model has the highest performance for predicting UVVisible spectra with a training RMSE of 0.06 and validation RMSE of 0.08.

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

Initial Uncertainty Analysis of Carbon Tetrachloride Contamination and Remediation in the Ringold A and Lower Mud Units at the Central Plateau

The long-term effectiveness of groundwater cleanup at the Hanford Site Central Plateau depends on predictive models that can capture key uncertainties in contaminant fate and transport. Carbon tetrachloride (CCl 4 ), a persistent and toxic compound, presents particular challenges due to variability in degradation rates, uncertainty in initial plume distribution, and subsurface heterogeneity. These uncertainties directly influence plume persistence, migration pathways, and remedy performance, and thus must be systematically evaluated to support long-term remediation planning. To address these gaps, a large-scale Monte Carlo analysis was conducted using the Plateau to River (P2R) model framework. The modeling approach parameterized three primary uncertainty factors: (1) degradation rate, (2) initial plume distribution, and (3) hydraulic conductivity. Degradation was represented as a first-order process, with half-lives ranging from 70 to 700 years. Initial plume distributions were created using a geostatistical simulation method (sgsim), which generates many equally plausible versions of how contaminants might be distributed underground. From this, 100 different scenarios were mapped onto the P2R grid. Variability in hydraulic conductivity was represented in a similar way, with 100 scenarios each for the Ringold Lower Mud and Ringold A units (layers 6 and 7), based on fitted exponential variograms and conditioned to well data. In total, more than 1000 realizations were simulated to assess plume behavior under uncertainty. Results demonstrate that degradation kinetics exert the strongest control over plume persistence: Shorter half-lives produced rapid mass reduction, while longer half-lives yielded persistent plumes with limited attenuation. A nonlinear response was observed, with steep mass reductions at half-lives greater than 200 years and near-linear declines beyond this threshold, reflecting interactions between degradation and pumping. The initial plume distribution strongly influenced early transport patterns, with broader sources generating larger plume footprints, although pump-and-treat operations constrained plume migration to managed areas. By comparison, hydraulic conductivity variability in the Ringold units had only a secondary influence, modifying spreading behavior without altering the dominant migration pathways governed by source configuration and hydraulic controls. Overall, the analysis highlights that uncertainty in degradation rate and initial plume configuration are the primary drivers of variability in plume predictions, while conductivity heterogeneity plays a limited role. These findings underscore the need for improved site-specific data on degradation processes and source characterization to enhance the reliability of long-term performance assessments and to better inform remedial decision-making at the Central Plateau.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2025

Extreme flood return levels in a U.S. mid-Atlantic estuary using 40-year fluvial-coastal model simulations

Using an integrated watershed-coastal modeling framework, we conducted long-term historical simulations (1980-2019) of fluvial and coastal flooding in the Delaware Bay and River, a vulnerable estuarine system in the U.S., at high spatial resolutions. By applying process-based models, we generated physically consistent and spatially detailed estimates of estuarine, riverine, and surge-driven extreme water level compared to previous studies that used field data only. We then evaluated changes in the magnitude of flood events using the 40-year simulations and detrended Floodwater Depth values with stationary extreme value analysis. Our detailed assessment of spatial-varying extreme values revealed how different flood-generation mechanisms can dominate various zones in the estuary. The datasets produced through this work will be valuable for long-term flood hazard mitigation planning in coastal communities in the Delaware Bay and River region. Additionally, this work will serve as a benchmark for other coastal flood hazard modeling communities worldwide, aiding them in systematically modeling long-term and continuous extreme flood events.

Coastal and Riverine Flooding↗
DOE OSTI2024

Preparation of Neptunyl and Plutonyl Acetates To Access Nonaqueous Transuranium Coordination Chemistry

Uranyl diacetate dihydrate is a useful reagent for the preparation of uranyl (UO 2 2+ ) coordination complexes, as it is a well-defined stoichiometric compound featuring moderately basic acetates that can facilitate protonolysis reactivity, unlike other anions commonly used in synthetic actinide chemistry such as halides or nitrate. Despite these attractive features, analogous neptunium (Np) and plutonium (Pu) compounds are unknown to date. Here, in this study, a modular synthetic route is reported for accessing stoichiometric neptunyl(VI) and plutonyl(VI) diacetate compounds that can serve as starting materials for transuranic coordination chemistry. The new NpO 2 2+ and PuO 2 2+ complexes, as well as a corresponding molecular UO 2 2+ complex, are isomorphous in the solid state, and in solution show similar solubility properties that facilitate their use in synthesis. In both solid and solution state, the +VI oxidation state (O.S.) is maintained, as demonstrated by vibrational and optical spectroscopy, confirming that acetate anions stabilize the oxidizing, high-valent +VI states of Np and Pu as they do for the more stable U(VI). All three acetate salts readily react with a model diprotic ligand, affording incorporation of U(VI), Np(VI), and Pu(VI) cores into molecular coordination compounds that occurs concomitantly with elimination of acetic acid; the new complexes are high-valent, yet overall charge neutral, facilitating entry into nonaqueous chemistry by rational synthesis. Computational studies reveal that the dianionic ligand framework assists in stabilizing the +VI O.S. via donation to the 5f shells of the actinides, highlighting the potential usefulness of protonolysis reactivity toward preparation of stabilized high-valent transuranic species.

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

Integrated Quantum-Classical Protocol for the Realistic Description of Solvated Multinuclear Mixed-Valence Transition-Metal Complexes and Their Solvatochromic Properties

Linear cyanide-bridged polymetallic complexes, which undergo photoinduced metal-to-metal charge transfer, represent prototypical systems for studying long-range electron-transfer reactions and understanding the role played by specific solute–solvent interactions in modulating the excited-state dynamics. Here, to tackle this problem, while achieving a statistically meaningful description of the solvent and of its relaxation, one needs a computational approach capable of handling large polynuclear transition-metal complexes, both in their ground and excited states, as well as the ability to follow their dynamics in several environments up to nanosecond time scales. Here, we present a mixed quantum classical approach, which combines large-scale molecular dynamics (MD) simulations based on an accurate quantum mechanically derived force field (QMD-FF) and self-consistent QMD polarized point charges, with IR and UV–vis spectral calculations to model the solvation dynamics and optical properties of a cyano-bridged trinuclear mixed-valence compound (trans-[(NC) 5 Fe III (μ-CN)Ru II (pyridine) 4 (μ-NC)Fe III (CN) 5 ] 4– ). We demonstrate the reliability of the QMD-FF/MD approach in sampling the solute conformational space and capturing the local solute–solvent interactions by comparing the results with higher-level quantum mechanics/molecular mechanics (QM/MM) MD reference data. The IR spectra calculated along the classical MD trajectories in different solvents correctly predict the red shift of the CN stretching band in the aprotic medium (acetonitrile) and the subtle differences measured in water and methanol, respectively. By explicitly including the solvent molecules around the cyanide ligands and calculating the thermal averaged absorption spectra using time-dependent density functional theory calculations within the Tamm–Dancoff approximation, the experimental solvatochromic shift is quantitatively reproduced going from water to methanol, while it is overestimated for acetonitrile. This discrepancy can likely be traced back to the lack of important dispersion interactions between the solvent cyano groups and the pyridine substituents in our micro solvation model. The proposed protocol is applied to the ground state in water, methanol, and acetonitrile and can be flexibly generalized to study excited-state nonequilibrium solvation dynamics.

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

Evaluation of GlassNet for physics-informed machine learning of glass stability and glass-forming ability

Glassy materials form the basis of many modern applications, including nuclear waste immobilization, touch-screen displays, and optical fibers, and also hold great potential for future medical and environmental applications. However, their structural complexity and large composition space make design and optimization challenging for certain applications. Of particular importance for glass processing and design is an estimate of a given composition's glass-forming ability (GFA). However, there remain many open questions regarding the underlying physical mechanisms of glass formation, especially in oxide glasses. It is apparent that a proxy for GFA would be highly useful in glass processing and design, but identifying such a surrogate property has proven itself to be difficult. While glass stability (GS) parameters have historically been used as a GFA surrogate, recent research has demonstrated that most of these parameters are not accurate predictors of the GFA of oxide glasses. Here, in this work, we explore the application of an open-source pre-trained neural network model, GlassNet, that can predict the characteristic temperatures necessary to compute GS with reasonable performance and assess the feasibility of using these physics-informed machine learning (PIML)-predicted GS parameters to estimate GFA. In doing so, we track the uncertainties at each step of the computation—from the original ML prediction errors to the compounding of errors during GS estimation, and finally to the final estimation of GFA. While GlassNet exhibits reasonable accuracy on all individual properties, we observe a large compounding of error in the combination of these individual predictions for the PIML prediction of GS, finding that random forest models offer similar accuracy to GlassNet. We also break down the performance of GlassNet on different glass families and find that the error in GS prediction is correlated with the error in crystallization peak temperature prediction. Lastly, we utilize this finding to assess the relationship between top-performing GS parameters and GFA for two ternary glass systems: sodium borosilicate and sodium iron phosphate glasses. We conclude that to obtain true ML predictive capability of GFA, significantly more data needs to be collected.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Novel Microbial Routes to Synthesize Industrially Significant Precursor Compounds

Ethylene is the most widely employed organic precursor compound in industry. The potential to impact ethylene formation via recently discovered microbial processes is tenable using plentiful CO2 feedstocks. The overall long-term objective of this project was to develop an industrially compatible microbial process to synthesize ethylene in high yields. The key objective of this project was to fully define and initially characterized a recently discovered and genetically regulated anaerobic pathway to produce high levels of ethylene called the Dihydroxyacetone Phosphate - Ethylene Pathway in phototrophic bacteria. This was addressed through the following specific aims: 1. Fully probe the catalytic potential of all enzymes of the DHAP ethylene pathway and determine the regulatory mechanism of DHAP-ethylene pathway gene expression. 2. Discover effective and active ethylene enzymes encoded in cultured and uncultured organisms from anoxic environments. 3.Model the thermodynamics and kinetics of ethylene synthetic pathways to guide engineering efforts in integrating best performing DHAP-ethylene pathway enzymes into model bacteria chassis for enhance ethylene yields. Through this project we discovered the initially missing genetic and enzyme component of the DHAP-ethylene pathway that directly synthesized ethylene and other important industrial compounds like methane and ethane from specific substrates. We uncovered and partially characterized a nitrogenase-like reductase that functions in DHAP-ethylene pathway specifically and in methionine synthesis in general. This nitrogenase-like system is called the Methylthio-Alkane Reductase (MAR) for its ability to cleave volatile organic sulfur compounds into methanethiol (CH3-SH) for methionine synthesis and a hydrocarbon byproduct. Key to the DHAP-ethylene pathway, MAR is the essential enzyme that cleaves 2-methylthioethanol (CH3-S-CH2-CH2-OH) into ethylene. Coordinately, we uncovered that the MAR genes and genes associated with conversion of methanethiol (CH3-SH) to methionine are under genetic control of a LysR Type Transcriptional Regulator called SalR, whose activity is dependent upon the amount of sulfate available to the cell. When sulfate as the preferred sulfur source for cell growth drops below 200 micromolar, SalR become active for expressing the MAR and methionine biosynthesis genes to enable the cell to grow from volatile organic sulfur compounds and make ethylene. Metabolic thermos-kinetic modeling revealed that these MAR reactions for ethylene and other hydrocarbon production are highly thermodynamically favorable and are one of the largest driving forces for ethylene production by the DHAP-ethylene pathway for high ethylene yields. Modeling also indicated that a key aldolase and to a lesser extent an isomerase of the DHAP-ethylene pathway for production of the ethylene precursor, 2-methylthioethanol, also would increase ethylene yields. Through metagenomic mining and gene synthesis by the JGI DNA synthesis program, over 500 aldolase and isomerase homologs were synthesized and screened. From this, variants were uncovered with substantially higher activity that increased ethylene yields 5-fold via the aldolase reaction and 1.5-fold via the isomerase reaction. Each of these elements that increase ethylene production were integrated together via plasmid under appropriate gene promoter elements in the phototrophic bacterium, Rhodospirillum rubrum, resulting in at least 3 orders of magnitude increase in ethylene yield from carbon dioxide feedstock.

10 SYNTHETIC FUELS↗
DOE OSTI2020

Substrate-explicit Modeling Tutorial

This tutorial Narrative demonstrates how to apply thermodynamic theory, also known as lambda theory, to convert molecular formulas of compounds (derived from FTICR-MS peaks using Formulatiry and R codes) into stoichiometric and kinetic forms of biogeochemical reactions and how to use the resulting kinetic equations to simulate dynamic conversion of compounds in batch and continuous stirred tank reactors. It will address under what conditions respiration rates are driven by thermodynamics, and how respiration rates respond to the variations in parameters and input variables and how to interpret the results.

54 ENVIRONMENTAL SCIENCES↗