Engineering topics
Publications and source records attributed to St. John, Peter C..
Computational predictions of vaporization properties aid the de novo design of green chemicals, including clean alternative fuels, working fluids for efficient thermal energy recovery, and polymers that are easily degradable and recyclable.
Bond dissociation energetics underpin the thermodynamics of chemical transformations where bonds are broken or formed and can also be used to predict reaction rates and selectivities.
Not Available
Designing alternative fuels for advanced compression ignition engines necessitates a predictive model for cetane number (CN). In this study, the physics-informed graph neural networks are introduced for a reliable CN prediction by considering molecular features pertinent to the physical properties of molecules that affect CN. The reliability of measured data is another key factor to consider for improving the predictive model. Various experimental instruments for measuring CN exist, including standard and non-standard methods. In this regard, a systematic data quality analysis was carried out for the total 630 CNs collected from literature and new measurements in this study using Advanced Fuel Ignition Delay Analyzer (AFIDA). The results from this data curation process were reflected in the model by imposing lower sample weights on the data coming from less reliable measurement techniques. This approach effectively maximized the prediction accuracy while incorporating data from all available sources. Using the sample weights decreased the mean absolute error (MAE) up to 0.8 CN units. The accuracy was also improved by introducing the CN-related physical properties (the number of hydrogen bond donors and acceptors); the test set MAE is 5.74 and 7.01 for the model with and without such properties, respectively. Investigating molecular structural effects on CN was also carried out to gain chemical insights into factors used to design new fuel candidates. The dimensionality reduction analysis of feature vectors showed a clear clustering in terms of functional groups and CN and the structural effect derived from the model was consistent with the physicochemical insights. Finally, this physics-informed model and data curation would be helpful for accurate CN prediction and inform rational fuel design.
This study developed the design principle of bio-derived ether fuels with low soot emission and high reactivity based on an experiment-theory combined approach.
Muconic acid is a bioprivileged molecule that can be converted into direct replacement chemicals for incumbent petrochemicals and performance-advantaged bioproducts. In this study, Pseudomonas putida KT2440 is engineered to convert glucose and xylose, the primary carbohydrates in lignocellulosic hydrolysates, to muconic acid using a model-guided strategy to maximize the theoretical yield. Using adaptive laboratory evolution (ALE) and metabolic engineering in a strain engineered to express the D-xylose isomerase pathway, we demonstrate that mutations in the heterologous D-xylose:H + symporter (XylE), increased expression of a major facilitator superfamily transporter (PP_2569), and overexpression of aroB encoding the native 3-dehydroquinate synthase, enable efficient muconic acid production from glucose and xylose simultaneously. Using the rationally engineered strain, we produce 33.7 g L -1 muconate at 0.18 g L -1 h -1 and a 46% molar yield (92% of the maximum theoretical yield). This engineering strategy is promising for the production of other shikimate pathway-derived compounds from lignocellulosic sugars.
Molecular descriptors encoding kinetic and thermodynamic stabilization capture the difference between transient and persistent organic radicals.
Optimizing the metabolism of microbial cell factories for yields and titers is a critical step for economically viable production of bioproducts and biofuels. In this process, tuning the expression of individual enzymes to obtain the desired pathway flux is a challenging step, in which data from separate multiomics techniques must be integrated with existing biological knowledge to determine where changes should be made. Following a design-build-test-learn strategy, building on recent advances in Bayesian metabolic control analysis, we identify in this work key enzymes in the oleaginous yeast Yarrowia lipolytica that correlate with the production of itaconate by integrating a metabolic model with multiomics measurements. To this extent, we quantify the uncertainty for a variety of key parameters, known as flux control coefficients (FCCs), needed to improve the bioproduction of target metabolites and statistically obtain key correlations between the measured enzymes and boundary flux. Based on the top five significant FCCs and five correlated enzymes, our results show phosphoglycerate mutase, acetyl-CoA synthetase (ACSm), carbonic anhydrase (HCO3E), pyrophosphatase (PPAm), and homoserine dehydrogenase (HSDxi) enzymes in rate-limiting reactions that can lead to increased itaconic acid production.
Particulate matter (PM) indices - those linking PM emissions from gasoline engines to the composition and properties of the fuel - have been a topic of significant study over the last decade. It has long been known that fuel composition has a significant impact on particulate emissions from gasoline engines. Since gasoline direct injection (GDI) engines have become the market-leading technology, this has become more significant because the evaporative behavior of fuel increases in importance. Several PM indices have been developed to provide metrics describing this behavior and correlating PM emissions. In this article, 16 different PM indices are identified and collected - to the authors' knowledge, all of the indices are available at the time of writing. The indices are reviewed and discussed in the context of the information required to calculate them, as well as their utility. Additionally, the authors believe that there is a need for indices that provide both a detailed and robust correlation, as well as those that are less sophisticated yet sufficient for specific use cases. Future research is suggested to guide the technical community toward improvements in the indices' methods and equations for both high and low fidelity and high and low time investment.
From quantum chemical and experimental NMR data, a 3D graph neural network, CASCADE, has been developed to predict carbon and proton chemical shifts. Stereoisomers and conformers of organic molecules can be correctly distinguished.
Prior engineering of the ethanologen Zymomonas mobilis has enabled it to metabolize xylose and to produce 2,3-butanediol (2,3-BDO) as a dominant fermentation product. When co-fermenting with xylose, glucose is preferentially utilized, even though xylose metabolism generates ATP more efficiently during 2,3-BDO production on a BDO-mol basis. To gain a deeper understanding of Z. mobilis metabolism, we first estimated the kinetic parameters of the glucose facilitator protein of Z. mobilis by fitting a kinetic uptake model, which shows that the maximum transport capacity of glucose is seven times higher than that of xylose, and glucose is six times more affinitive to the transporter than xylose. With these estimated kinetic parameters, we further compared the thermodynamic driving force and enzyme protein cost of glucose and xylose metabolism. It is found that, although 20% more ATP can be yielded stoichiometrically during xylose utilization, glucose metabolism is thermodynamically more favorable with 6% greater cumulative Gibbs free energy change, more economical with 37% less enzyme cost required at the initial stage and sustains the advantage of the thermodynamic driving force and protein cost through the fermentation process until glucose is exhausted. Glucose-6-phosphate dehydrogenase (g6pdh), glyceraldehyde-3-phosphate dehydrogenase (gapdh) and phosphoglycerate mutase (pgm) are identified as thermodynamic bottlenecks in glucose utilization pathway, as well as two more enzymes of xylose isomerase and ribulose-5-phosphate epimerase in xylose metabolism. Acetolactate synthase is found as potential engineering target for optimized protein cost supporting unit metabolic flux. Pathway analysis was then extended to the core stoichiometric matrix of Z. mobilis metabolism. Growth was simulated by dynamic flux balance analysis and the model was validated showing good agreement with experimental data. Dynamic FBA simulations suggest that a high agitation is preferable to increase 2,3-BDO productivity while a moderate agitation will benefit the 2,3-BDO titer. Taken together, this work provides thermodynamic and kinetic insights of Z. mobilis metabolism on dual substrates, and guidance of bioengineering efforts to increase hydrocarbon fuel production.
Machine-readable chemical structure representations are foundational in all attempts to harness machine learning for the prediction of reactivities, selectivities, and chemical properties directly from molecular structure. The featurization of discrete chemical structures into a continuous vector space is a critical phase undertaken before model selection, and the development of new ways to quantitatively encode molecules is an active area of research. Here, we highlight the application and suitability of different representations, from expert-guided “engineered” descriptors to automatically “learned” features, in different prediction tasks relevant to organic and organometallic chemistry, where differing amounts of training data are available. These tasks include statistical models of stereo- and enantioselectivity, thermochemistry, and kinetics developed using experimental and quantum chemical data. The use of expert-guided molecular descriptors provides an opportunity to incorporate chemical knowledge, domain expertise, and physical constraints into statistical modeling. In applications to stereoselective organic and organometallic catalysis, where data sets may be relatively small and 3D-geometries and conformations play an important role, mechanistically informed features can be used successfully to obtain predictive statistical models that are also chemically interpretable. We provide an overview of several recent applications of this approach to obtain quantitative models for reactivity and selectivity, where topological descriptors, quantum mechanical calculations of electronic and steric properties, along with conformational ensembles, all feature as essential ingredients of the molecular representations used. Alternatively, more flexible, general-purpose molecular representations such as attributed molecular graphs can be used with machine learning approaches to learn the complex relationship between a structure and prediction target. This approach has the potential to out-perform more traditional representation methods such as “hand-crafted” molecular descriptors, particularly as data set sizes grow. One area where this is particularly relevant is in the use of large sets of quantum mechanical data to train quantitative structure–property relationships. A general approach toward curating useful data sets and training highly accurate graph neural network models is discussed in the context of organic bond dissociation enthalpies, where this strategy outperforms regression using precomputed descriptors. Finally, we describe how graph neural network predictions can be incorporated into mechanistically informed statistical models of chemical reactivity and selectivity. Once trained, this approach avoids the expensive computational overhead associated with quantum mechanical calculations, while maintaining chemical interpretability. We illustrate examples for which fast predictions of bond dissociation enthalpy and of the identities of radicals formed through cleavage of a molecule’s weakest bond are used in simple physical models of site-selectivity and reactivity.
Not Available
Hydroxymethylfurfural (HMF) is one of the important renewable platform compounds that can be obtained from biomass feedstocks through glucose conversion catalyzed by Brønsted and Lewis acids. However, it is challenging to enhance the HMF yield due to side reactions. In this study, a systematic approach combining theory and experiment was performed to investigate the influence of Lewis acids and organic solvents on the HMF yield. For the Lewis acid effect, a relationship between chemical hardness and experimental HMF yields was found in the rate-limiting step of glucose-to-fructose isomerization for six metal chlorides; HMF production was promoted when the metal chloride and a substrate had a similar chemical hardness. To study the organic solvent effect, a multivariate model was developed based on the insights gained from the mechanistic study of fructose dehydration, to predict HMF yields in a given water-organic cosolvent system. It showed a reliable accuracy in evaluating HMF yields with a mean absolute error (MAE) of 3.0% with respect to experimental HMF yields for 13 solvents, and also predicted HMF yields with a MAE of 10.7% for four new solvents. Chemical interpretation of the model revealed that it is desirable to use a solvent capable of stabilizing the carbocation intermediates with low proton transfer activity and high hydrogen bond basicity, to maximize the HMF yield. This multivariate model informs experimentalists about rational selection of solvents with very low computational costs needed to calculate only six variables for each solvent. It can be expanded to other catalytic systems such as heterogeneous Brønsted–Lewis bifunctional catalysts and enables optimization of reaction conditions to obtain other useful platform molecules through biomass conversion.
$\phi$-sensitivity is the change in ignition delay time (IDT) with respect to the fuel-to-air equivalence ratio ($\phi$). High $\phi$-sensitivity is a desirable fuel property for applications in advanced compression ignition and multi-mode engine designs. Understanding how $\phi$-sensitivity depends on chemical structure is essential for selecting promising biofuels from the ever-growing list of proposed candidates. Here, we investigate the effect of chemical structure on $\phi$-sensitivity with experiment, simulation, and theory. Experimental Advanced Fuel Ignition Delay Analyzer (AFIDA) measurements for 2,4-dimethylpentane and diisopropyl ether provide evidence that branching and functional groups strongly impact $\phi$-sensitivity. Further insights into this dependence are obtained with 0-D kinetic simulations with existing mechanisms for n-pentane, diethyl ether, 3-pentanone, n-heptane, 2-methylhexane, 2,4-dimethylpentane, and 2,2,3-trimethylbutane. Quantum mechanical (QM) G4 calculations of low-temperature reactions help explain the observed experimental and simulation trends. Specifically, these QM calculations provide theoretical estimates of the ketohydroperoxide (KHP) dissociation rates, the HO 2 formation rates from peroxy radical (ROO), and the “cross-over” temperatures, i.e., the temperature at which ROO dissociation is favored compared to hydroperoxyl radical (QOOH) formation. Each of these reaction rates is compared to the n-alkane reference point to determine the impact of branching and different functional groups. Although kinetic mechanisms typically assume that KHP dissociation rates are invariant of chemical environment, our QM results suggest that this rate can span a range of roughly two orders of magnitude. We also discuss the importance of including the peroxy-hydroperoxy (OO-OOH) hydrogen transfer reaction for branched ethers. Finally, the insights gained assist in proposing a highly $\phi$-sensitive compound, namely, isopropyl propyl ether.
Sooting propensity, a measurement of how much particulate matter is produced when a fuel is burned, is a property of significant interest among researchers who are striving to discover the next generation of cleaner, more efficient fuels and fuel additives. Many compounds are not viable as fuels and/or fuel additives, and as a result, designing cleaner-burning biofuels using only experimental techniques is inefficient. Predictive models have been instrumental in reducing this inherent difficulty, providing researchers with a tool to preemptively screen compounds before production and testing. The present work compares the accuracies and interpretabilities of existing models used to predict a particular measure of sooting propensity, Yield Sooting Index (YSI). These models include artificial neural networks, graph neural networks, and multivariate equations. A novel equation for predicting YSI based on atom path count and bond order is proposed, which can highlight key structural components that contribute to YSI. It was found that artificial neural networks slightly outperform graph neural networks and greatly outperform multivariate equations in blind (test set) prediction accuracy; however, graph neural networks and multivariate equations provide significantly more interpretability as to how compound structure relates to YSI. Predictions of YSI are compared to experimental measurements for previously un-tested compounds with cetane numbers comparable to diesel fuel (50-60) (butyl decanoate, ethyl decanoate, 1,4-bis(ethenoxymethyl)cyclohexane, and 5-heptyloxolan-2-one), and it was found that these compounds produce significantly less soot compared to diesel fuel.
Yield Sooting Index (YSI) measurements have shown that oxygenated aromatic compounds (OACs) tend to have lower YSI than aromatic hydrocarbon (AHC) compounds. For example, typical AHCs such as toluene and ethyl benzene have a YSI of 170 and 216, respectively, in contrast, OACs such as phenol and anisole have a YSI of 81 and 111, respectively. However, this trend is not always true as was observed for the structural isomers 1-phenylethanol (1PE, YSI=142) and 2-phenylethanol (2PE, YSI=207), where 2PE contains a YSI more representative of AHCs than OACs. We applied flow reactor experiments and density functional theory (DFT) calculations to examine how oxygen functionality present in 1PE and 2PE alters the reaction pathways leading to the observed difference in soot formation. It was determined that the proximity of the oxygen functional group to the aromatic ring determines whether the oxygen remains attached to the primary reacting species (for 1PE) or is eliminated early in the combustion sequence (for 2PE). For these alcohols, preservation of the oxygen in the molecule leads to further OACs, while loss of the oxygen leads to AHCs and benzyl radical. The direct pathways to AHCs and benzyl radical result in the higher YSI observed for 2PE.