Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Accounting”

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 523 records · Page 29

Partitioning Quantum Chemistry Simulations with Clifford Circuits

Current quantum computing hardware is restricted by the availability of only few, noisy qubits which limits the investigation of larger, more complex molecules in quantum chemistry calculations on quantum computers in the near term. Here, in this work, we investigate the limits of their classical and near-classical treatment while staying within the framework of quantum circuits and the variational quantum eigensolver. To this end, we consider naive and physically motivated, classically efficient product ansatz for the parametrized wavefunction adapting the separable-pair ansatz form. We combine it with post-treatment to account for interactions between subsystems originating from this ansatz. The classical treatment is given by another quantum circuit that has support between the enforced subsystems and is folded into the Hamiltonian. To avoid an exponential increase in the number of Hamiltonian terms, the entangling operations are constructed from purely Clifford or near-Clifford circuits. While Clifford circuits can be simulated efficiently classically, they are not universal. In order to account for missing expressibility, near-Clifford circuits with only few, selected non-Clifford gates are employed. The exact circuit structure to achieve this objective is molecule-dependent and is constructed using simulated annealing and genetic algorithms. We demonstrate our approach on a set of molecules of interest and investigate the extent of our methodology’s reach.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The many-body expansion for aqueous systems revisited: 1. Water-water interactions

We revisit the Many-Body Expansion (MBE) for water-water interactions by examining the effects of the basis set, including those resulting from the Basis Set Superposition Error (BSSE) correction, on the various terms for selected sizes of water clusters up to n = 21. The analysis is performed at the second order Møller-Plesset (MP2) perturbation theory with the family of augmented correlation consistent basis sets up to five zeta quality for the (H2O)n, n = 7, 10, 13, 16 and 21 clusters, for which we report either the complete MBE (n = 7, 10) or the one through the 6-body (n = 13) and the 5-body terms (n = 16, 21). Our results suggest that any sizeable contributions to the total cluster binding energy arising from the 5-body and larger terms are solely an artifact of the size of the basis set. Indeed, all terms above the 4-body converge monotonically to practically zero at the Complete Basis Set (CBS) limit. The magnitudes of the 5-body and above terms are accurately reproduced even with the smaller basis set of the series (aug-cc-pVDZ) once the BSSE correction is taken into account. The same level of theory (MP2/aug-cc-pVDZ, BSSEcorrected) also accurately reproduces the MP2/CBS values of the 3- and 4-body terms. The contribution of electron correlation to the 3- and 4-body terms is quite small, so that neglecting the correlation contribution in all terms above the 3-body results in an error of ??(0.1%). The BSSE correction to the largest 2-body term in the MBE was accurately estimated from the function ?? [1 + erf(-?? ??)], which is proportional to the common (overlapping) area between two Gaussian distributions whose centers are separated by R, with the constants ??, ?? fitted to the calculated BSSE corrections for the individual 2-body terms of the clusters with each basis set. Our results unquestionably demonstrate that the MBE for water-water interactions converges to the 4-body term, since any finite terms above the 4-body are artifacts of the size of the basis set and can be safely truncated at the 4-body term when either a very large basis set is used or BSSE corrections are taken into account even with the aug-cc-pVDZ basis set. We expect these findings to have important consequences in the pursuit of accurate many-body molecular dynamics simulations for aqueous systems.

Heindel, Joseph↗

Reactive Molecular Dynamics Simulations and Quantum Chemistry Calculations To Investigate Soot-Relevant Reaction Pathways for Hexylamine Isomers

Sooting tendencies of a series of nitrogen-containing hydrocarbons (NHCs) have been recently characterized experimentally using the yield sooting index (YSI) methodology. This work aims to identify soot-relevant reaction pathways for three selected C6H15N amines, namely, dipropylamine (DPA), diisopropylamine (DIPA), and 3,3-dimethylbutylamine (DMBA) using ReaxFF molecular dynamics (MD) simulations and quantum mechanical (QM) calculations and to interpret the experimentally observed trends. ReaxFF MD simulations are performed to determine the important intermediate species and radicals involved in the fuel decomposition and soot formation processes. QM calculations are employed to extensively search for chemical reactions involving these species and radicals based on the ReaxFF MD results and also to quantitatively characterize the potential energy surfaces. Specifically, ReaxFF simulations are carried out in the NVT ensemble at 1400, 1600, and 1800 K, where soot has been identified to form in the YSI experiment. These simulations account for the interactions among test fuel molecules and pre-existing radicals and intermediate species generated from rich methane combustion, using a recently proposed simulation framework. ReaxFF simulations predict that the reactivity of the amines decrease in the order DIPA > DPA > DMBA, independent of temperature. Both QM calculations and ReaxFF simulations predict that C2H4, C3H6, and C4H8 are the main nonaromatic soot precursors formed during the decomposition of DPA, DIPA, and DMBA, respectively, and the associated reaction pathways are identified for each amine. Both theoretical methods predict that sooting tendency increases in the order DPA, DIPA, and DMBA, consistent with the experimentally measured trend in YSI. This work demonstrates that sooting tendencies and soot-relevant reaction pathways of fuels with unknown chemical kinetics can be identified efficiently through combined ReaxFF and QM simulations. Overall, predictions from ReaxFF simulations and QM calculations are consistent, in terms of fuel reactivity, major intermediates, and major nonaromatic soot precursors.

fuel decomposition↗

Mechanistic Investigation of Dimethylmercury Formation Mediated by a Sulfide Mineral Surface

Mercury (Hg) pollution is a global environmental problem, but many key transformations among 29 various Hg species are not well understood. Although anaerobic microorganisms are the primary 30 producers of monomethylmercury (MMeHg), the abiotic formation of dimethylmercury (DMeHg) 31 from MMeHg have been suggested to account for a large portion of DMeHg in the ocean. It has 32 been shown that abiotic formation of DMeHg from MMeHg can be facilitated by reduced sulfur 33 groups on mineral surfaces. However, a detailed mechanistic understanding of this reaction is 34 lacking. Here we perform density functional theory (DFT) calculations to explore mechanisms of 35 DMeHg formation on the 110 surface of a CdS(s) (hawleyite) nanoparticle. We show that 36 coordination of MMeHg substituents to adjacent reduced sulfur groups protruding from the surface 37 facilitate DMeHg formation and that the reaction proceeds through direct transmethylation from 38 one MMeHg substituent to another. Solvation effects were found to play an important role in the 39 surface reconstruction of the nanoparticle and in decreasing the energetic barrier for DMeHg 40 formation. We also investigate DMeHg formation from MMeHg mediated by organic dithiols and 41 compare the results to those mediated by the mineral surface. In contrast to the surface reaction, 42 the calculations suggest that dithiol-mediated formation of DMeHg instead proceeds through a 43 mechanism in which Hg–Hg interactions activate a methyl group from one MMeHg substituent 44 for intramolecular transfer to the adjacent MMeHg substituent. These findings refine previously 45 proposed mechanisms of abiotic DMeHg formation and fill gaps in our understanding of global 46 Hg cycling in the environment.

Lian, Peng↗

Diffusion-Limited Kinetics in Reactive Systems

A proper representation of chemical kinetics is vital to understanding, modeling, and optimizing many important chemical processes. In liquid and surface phases, where diffusion is slow, the rate at which the reactants diffuse together limits the overall rate of many elementary reactions. Commonly, the textbook Smoluchowski theory is utilized to estimate effective rate coefficients in the liquid phase. On surfaces, modelers commonly resort to much more complex and expensive Kinetic Monte Carlo (KMC) simulations. Here, in this study, we extend the Smoluchowski model to allow the diffusing species to undergo chemical reactions and derive analytical formulas for the diffusion-limited rate coefficients for 3D, 2D, and 2D/3D interface cases. With these equations, we are able to demonstrate that when species react faster than they diffuse they can react orders of magnitude faster than predicted by Smoluchowski theory, through what we term “the reactive transport effect”. We validate the derived steady-state equations against particle Monte Carlo (PMC) simulations, KMC simulations, and non-steady-state solutions. Furthermore, using PMC and KMC simulations, we propose corrections that agree with all limits and the computed data for the 2D and 2D/3D interface steady-state equations, accounting for unique limitations in the associated derived equations. Additionally, we derive equations to handle couplings between diffusion-limited rate coefficients in reaction networks. We believe these equations should make it possible to run much more accurate mean-field simulations of liquids, surfaces, and liquid–surface interfaces accounting for diffusion limitations and the reactive transport effect.

Johnson, Matthew Sean↗

Modeling Pb(II) Adsorption on Mineral Surfaces: Bridging Density Functional Theory and Experiment with Thermodynamic Insights

Despite decades of work on aqueous lead (Pb) adsorption on a-Fe2O3 (hematite) and a-Al2O3 (alumina), gaps between measurements and modeling obscure molecular-level understanding. Achieving well-matched geometries between theory and experimental for mineral-water interfaces is a hurdle, as surface functional group type and distribution must be accounted for in determining mechanisms. Additionally, computational methods that can describe the substrate are often not appropriate to capture aqueous effects. Progress requires focusing on well-studied and relevant systems, such as key facets (001),(012), and (110) of hematite and alumina, and ubiquitous contaminants such as aqueous Pb. In the past, bulk-parametrized bond-valence principles were used to rationalize Pb(II) adsorption trends. These approaches can break down at surfaces, where flexible bonding environments and adsorption-induced surface relaxations play a critical role. Here, we adapt and apply a density functional theory (DFT) and thermodynamics framework, integrating DFT-calculated energies with experimental data and electrochemical principles, to predict Pb(II) adsorption. Our model results capture trends across the full set of surfaces and predict that inner-sphere Pb(II) sorption on (001) alumina varies from unfavorable to weakly favorable across a range of pH conditions. This aligns with experiment insights that Pb(II) interacts at that surface through outer-sphere interactions. Extending to Fe(II) adsorption, we demonstrate a coverage-dependent site preference, potentially explaining disorder in overlayers grown by the oxidative adsorption of Fe(II) on hematite (001).

lead contamination↗

Hydration Free Energies of Polypeptides from Popular Implicit Solvent Models versus All-Atom Simulation Results Based on Molecular Quasichemical Theory

Calculating the hydration free energy of a macromolecule in all-atom simulations has long remained a challenge, necessitating the use of models wherein the effect of the solvent is captured without explicit account of solvent degrees of freedom. This situation has changed with developments in the molecular quasi-chemical theory (QCT)-an approach that enables calculation of the hydration free energy of macromolecules within all-atom simulations at the same resolution as is possible for small molecular solutes. The theory also provides a rigorous and physically transparent framework to conceptualize and model interactions in molecular solutions and thus provides a convenient framework to investigate the assumptions in implicit solvent models. In this study, we compare the results using molecular QCT versus predictions from EEF1, ABSINTH, and GB/SA implicit solvent models for polyglycine and polyalanine solutes covering a range of chain lengths and conformations. The hydration free energies or the differences in hydration free energies between conformers obtained from the implicit solvent models do not agree with explicit solvent results, with the deviations being largest for the group additive EEF1 and ABSINTH models. GB/SA does better in capturing the qualitative trends seen in explicit solvent results. Finally, analysis founded on QCT reveals the critical importance of the cooperativity of hydration that is inherent in the hydrophilic and hydrophobic contributions to hydration-physics that is not well captured in additive models but somewhat better accounted for by means of a dielectric in the GB/SA approach.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Using Metadynamics to Reveal Extractant Conformational Free Energy Landscapes

Understanding the impact of extractant functionalization on metal-binding energetics in liquid-liquid extraction is essential to guide the development of better separation processes. Traditionally, computational extractant design uses electronic structure calculations on metal-ligand clusters to determine the metal-binding energy of the lowest energy state. Although highly accurate, this approach does not account for all of the relevant physics encountered under experimental conditions. Such methodologies often neglect entropic contributions such as temperature effects and ligand flexibility, in addition to approximating solvent-extractant interactions with implicit solvent models. In this study, we use classical molecular dynamics simulations with an advanced sampling method, metadynamics, to map out extractant molecule conformational free energies in the condensed phase. Here we generate the complete conformational landscape in solution for a family of bidentate malonamide-based extractants with different functionalizations of the headgroup and the side chains. In particular, we show how such alkyl functionalization reshapes the free energy landscape, affecting the free energy penalty of organizing the extractant into the cis-like metal-binding conformation from the trans-like conformation of the free extractant in solution. Specifically, functionalizing alkyl tails to the center of the headgroup has a greater influence on increasing molecular rigidity and disfavoring the binding conformation than functionalizing side chains. These findings are consistent with trends in metal-binding energetics based on experimentally reported distribution ratios. We also consider a different bidentate extractant molecule, carbamoylmethylphosphine oxide, and show how the choice of solvent can further reshape the conformational energetic landscape. This study demonstrates the feasibility of using molecular dynamics simulations with advanced sampling techniques to investigate extractant conformational energetics in solution, which, more broadly, will enable extractant design that accounts for entropic effects and explicit solvation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Comment on “Microhydration of Biomolecules: Revealing the Native Structures by Cold Ion IR Spectroscopy”

This paper presents a re-examination of the conclusions of a study reported in The Journal of Physical Chemistry Letters (Saparbaev et al., 2021, 12, 907) that compared the structure of microsolvated ions formed by electrospray ionization to those formed in the gas-phase via a previously published cryogenic ion trap approach. We conducted additional experiments that clearly show that most of the observed differences in the IR spectra can be accounted for by considering the different spectroscopic action schemes used to obtain them. In particular, the presence of the D 2 -tag induces shifts in some of the N-H and O-H peaks which need to be carefully considered before comparing spectra. Furthermore, once these spectral effects are taken into account, we show that both clustering approaches yields similar clusters structures for the small GlyH + (H 2 O)n species. Using unimolecular reaction rate theory, we also show that for the small complexes considered here, only the gas-phase equilibrium distribution of conformers should be expected in both experimental approaches. In addition, the barrier heights necessary to kinetically trap high-energy conformers at 298K is explored using a series of model polyalanine chains.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A Multiscale Simulation Study of Influence of Morphology on Ion Transport in Block Copolymeric Ionic Liquids

In this work, we present results of a multiscale simulation study involving coarse-graining and reverse-mapping steps probing at an atomistic resolution the influence of morphology on ion transport in block copolymer polymeric ionic liquids. We provide a detailed description of the multiscale simulation methodology and then describe the results in the context of four microphase-separated morphologies: two lamella systems, a cylinder, and a gyroid morphology. For the framework adopted in this study (in which the total number of ions was maintained fixed), morphology had little influence on the fraction of ions segregating to the interface and the coordination and hopping characteristics of such interfacially present ions. Such results manifested in anion mobilities being insensitive to the morphology once the dimensionality of the morphology was accounted for. Overall, our results are consistent with the hypothesis that the mobility of ions in such microphase-separated morphologies, after accounting for the dimensionality effects, can be roughly correlated to a linear superposition based on the fraction of bulk anions, which possess mobilities unperturbed from their values in the conducting homopolymers, and the interfacial anions, which exhibit much lower mobilities.

36 MATERIALS SCIENCE↗

Multiscale Modeling and Experimental Insights into High-Temperature Soil Biodegradation Dynamics of Semi-Crystalline Poly(Lactic Acid) Nonwoven Fabrics

This study investigates the biodegradation of semi-crystalline poly(lactic acid) (PLA) nonwovens (NWs) in soil at 58 °C using both experimental and mathematical modeling approaches. The model utilizes a system of parabolic diffusion-reaction partial differential equations (PDEs) to elucidate chemical transformations over time and in space. It accounts for phenomena such as the diffusion of water and lactic acid monomers through the polymer matrix and into the surrounding soil, along with their microbial breakdown. It also accounts for the initial PLA crystallinity and predicts its evolution in time. The model is solved numerically for a single filament, and the results were used to shed light on PLA NW transformations observed in soil over a 180-day incubation period. Various characterization techniques, including scanning electron microscopy (SEM), differential scanning calorimetry (DSC), and Raman spectroscopy, were employed to assess morphological changes, crystallinity, and molecular changes in the NWs throughout the experiment. By comparing the experimental data with the model predictions, the hydrolysis rate coefficient was found to be 3.37 × 10 -7 s -1 , while the rate of microbial degradation of lactic acid monomers was faster, of the order of 9.63 × 10 -7 s -1 . The findings highlight the significant role of crystallinity in the biodegradation process. The PLA degradation ceases when no amorphous material remains, and the crystallinity reaches 0.8, as observed in the experiments by day 120. Furthermore, this research contributes to a deeper understanding of PLA biodegradation dynamics and offers insights for effectively managing biodegradable materials in environmental settings.

Diffusion−reaction modeling↗

Can Li: A Career in Catalysis

Prof. Can Li is distinguished for his seminal achievements in fundamental and applied researches in catalysis, especially in advancing the characterization of catalysts and catalytic reactions related to solar energy conversion and fine chemical synthesis. The scope of his research is wide ranging, and he is leading some of the fields. This account highlights his major scientific achievements in advancing the fundamental sciences of characterizing catalysts and understanding the mechanisms of catalytic reactions by developing in situ spectroscopic techniques (e.g., infrared (IR), ultraviolet (UV), Raman, space- and time-resolved spectroscopies and photoelectric imaging), and applying the obtained knowledge in the rational design and synthesis of catalysts for practical applications in photocatalytic, electrocatalytic and photoelectrocatalytic water splitting and CO 2 reduction, chiral synthesis in nanoreactors and ultradeep desulfurization of fuels, and fine chemical synthesis. In particular, his demonstration of the kiloton scale solar fuel production project will be introduced as a successful model for converting fundamental science into practical and impactful applications to convert CO 2 and water into solar fuels. This account also highlights his services to the catalysis communities and international collaborations for promoting catalysis science worldwide.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Can Isoprene Oxidation Explain High Concentrations of Atmospheric Formic and Acetic Acid over Forests?

Formic and acetic acid concentrations are particularly high over forested areas of the world. However, the gas-phase mechanisms for producing these acids are poorly understood even for isoprene, the globally dominant biogenic hydrocarbon. We quantified formic and acetic acid production from reactions of hydroxyl radical (OH) (between high and low ranges of nitric oxide (NO) levels) with isoprene, methacrolein (MACR), isoprene epoxydiol (IEPOX), isoprene hydroxy hydroperoxide (ISOPOOH), and α-pinene from the focused isoprene experiments at California Institute of Technology (FIXCIT) laboratory chamber study. We find that (i) OH oxidation of MACR, IEPOX, and ISOPOOH are sources of formic acid, (ii) isoprene peroxy radical isomerization and associated photolysis oxidation products are potentially important sources of organic acids, and (iii) high levels of NO generally suppress organic acid formation from OH oxidation of isoprene. We modified existing chemical mechanisms for isoprene oxidation to account for organic acid production pathways observed in the FIXCIT study. We simulated organic acid production during the Southeastern Oxidant and Aerosol Study using the updated chemical mechanisms and represented acetic acid within a factor of 2 but still underpredicted formic acid by a factor of 6. While we cannot explain ambient formic acid with explicit chemical mechanisms, the FIXCIT results suggest that the oxidation of isoprene could account for as much as 70% of the global annual production of formic acid from gas-phase reactions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Modeling the Size Distribution and Chemical Composition of Secondary Organic Aerosols during the Reactive Uptake of Isoprene-Derived Epoxydiols under Low-Humidity Condition

Reactive uptake of isoprene epoxydiols (IEPOX), which are isoprene oxidation products, onto acidic sulfate aerosols is recognized to be an important mechanism for the formation of isoprene-derived secondary organic aerosol (SOA). While a mechanistic understanding of IEPOX-SOA formation exists, several processes affecting their formation remain uncertain. Evaluating mechanistic IEPOX-SOA models with controlled laboratory experiments under longer atmospherically relevant timescales is critical. Here, we implement our latest understanding of IEPOX-SOA formation within a box model to simulate the measured reactive uptake of IEPOX on polydisperse ammonium bisulfate seed aerosols within an environmental Teflon chamber. The model is evaluated with single-particle measurements of size distribution, volume, density, and composition of aerosols due to IEPOX-SOA formation at timescales of hours. Additionally, we find that the model can simulate the growth of particles due to IEPOX multiphase chemistry, as reflected in increases of the mean particle size and volume concentrations, and a shift of number size distribution to larger sizes. The model also predicts the observed evolution of particle number mean diameter and total volume concentrations at the end of the experiment. We show that in addition to self-limiting effects of IEPOX-SOA coatings, accounting for the molar balance between inorganic and organic sulfate and mass accommodation coefficient of IEPOX are important parameters governing modeling of IEPOX-SOA formation. Thus, models which do not account for the molar sulfate balance and/or diffusion limitations within IEPOX-SOA coatings are likely to predict too high IEPOX-SOA formation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Asphericity Can Cause Nonuniform Lithium Intercalation in Battery Active Particles

Uniform intercalation is desired to enable next-generation Li-ion batteries. While we expect nonuniformity in materials undergoing a phase change, single-phase intercalation materials such as nickel manganese cobalt oxide are believed to lithiate uniformly at the particle/electrolyte interface. However, recent imaging reveals nonuniform lithiation. Motivated by this discrepancy, we examine if aspherical particle shape can cause such nonuniformity since the conventional belief is based on spherical particle theory. We obtain real particle geometries using rapid lab-based X-ray computed tomography and subsequently perform physics-based calculations accounting for electrochemical reactions at the particle/electrolyte interface and lithium transport inside the particle bulk. The aspherical geometry breaks the symmetry and causes nonuniform reaction distribution. Such nonuniformity is exacerbated as the particle becomes more aspherical. We report the proposed mechanism represents a fundamental limit on achievable lithiation uniformity in aspherical particles in the absence of other mechanisms causing inhomogeneity, such as grain structure, nonuniform carbon-binder coating, etc.

25 ENERGY STORAGE↗

Beyond Ion Migration in Metal Halide Perovskites: Toward a Broader Photoelectrochemistry Perspective

Ion migration is a broad term used to account for the degradation of halide perovskite materials and devices. However, ion mobility is only one piece of the full picture-mobile ions/defects are first created, then transported, and eventually annihilated or immobilized. In this Perspective, we summarize emerging work that shows how tractable photochemical and Faradaic reactions provide a continuous source of ions to migrate. Furthermore, we discuss strategies to fundamentally manipulate ion migration by targeting specific electrochemical and reduction/oxidation mechanisms. This highlights the important role of defect photoelectrochemistry, as well as the soft nature of the perovskite lattice, in ion migration and self-healing. We conclude that distinguishing more detailed processes involved in “ion migration”, with an emerging focus on the reactions that form mobile ionic defects, is necessary to greatly improve the stability of devices and open up more technological applications.

defects↗

Influence of the Dielectric Constant on the Ionic Current Rectification of Bipolar Nanopores

In this paper, we investigate how the dielectric constant, ϵ, of an electrolyte solvent influences the current rectification characteristics of bipolar nanopores. It is well recognized that bipolar nanopores with two oppositely charged regions rectify current when exposed to an alternating electric potential difference. Here, we consider dilute electrolytes with NaCl only and with a mixture of NaCl and charged nanoparticles. These systems are studied using two levels of description, all-atom explicit water molecular dynamics (MD) simulations and coarse-grained implicit solvent MD simulations. The charge density and electric potential profiles and current-voltage relationship predicted by the implicit solvent simulations with ϵ = 11.3 show good agreement with the predictions from the explicit water simulations. Under nonequilibrium conditions, the predictions of the implicit solvent simulations with a dielectric constant closer to the one of bulk water are significantly different from the predictions obtained with the explicit water model. Further, these findings are closely aligned with experimental data on the dielectric constant of water when confined to nanometric spaces, which suggests that ϵ decreases significantly compared to its value in the bulk. Moreover, the largest electric current rectification is observed in systems containing nanoparticles when ϵ = 78.8. Using enhanced sampling, we have shown that this larger rectification arises from the presence of a significantly deeper minimum in the free energy of the system with a larger ϵ, and when a negative voltage bias is applied. Since implicit solvent models and mean-field continuum theories are often used to design Janus membranes based on bipolar nanopores, this work highlights the importance of properly accounting for the effects of confinement on the dielectric constant of the electrolyte solvent. The results presented here indicate that the dielectric constant in implicit solvent simulations may be used as an adjustable parameter to approximately account for the effects of nanometric confinement on aqueous electrolyte solvents.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Cell-Free Systems Biology: Characterizing Central Metabolism of Clostridium thermocellum with a Three-Enzyme Cascade Reaction

Genetic approaches have been traditionally used to understand microbial metabolism, but this process can be slow in nonmodel organisms due to limited genetic tools. An alternative approach is to study metabolism directly in the cell lysate. This avoids the need for genetic tools and is routinely used to study individual enzymatic reactions but is not generally used to study systems-level properties of metabolism. Here we demonstrate a new approach that we call “cell-free systems biology”, where we use well-characterized enzymes and multienzyme cascades to serve as sources or sinks of intermediate metabolites. This allows us to isolate subnetworks within metabolism and study their systems-level properties. To demonstrate this, we worked with a threeenzyme cascade reaction that converts pyruvate to 2,3-butanediol. Although it has been previously used in cell-free systems, its pH dependence was not well characterized, limiting its utility as a sink for pyruvate. We showed that improved proton accounting allowed better prediction of pH changes and that active pH control allowed 2,3-butanediol titers of up to 2.1 M (189 g/L) from acetoin and 1.6 M (144 g/L) from pyruvate. The improved proton accounting provided a crucial insight that preventing the escape of CO 2 from the system largely eliminated the need for active pH control, dramatically simplifying our experimental setup. We then used this cascade reaction to understand limits to product formation in Clostridium thermocellum, an organism with potential applications for cellulosic biofuel production. We showed that the fate of pyruvate is largely controlled by electron availability and that reactions upstream of pyruvate limit overall product formation.

09 BIOMASS FUELS↗