SEARCH · Search NASA

Results for “Non-equilibrium thermodynamics”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 91 records · Page 5

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Kinetics of solid–liquid interface motion in molecular dynamics and phase-field models: crystallization of chromium and silicon

The results of molecular dynamics (MD) simulations of the crystallization process in one-component materials and solid solution alloys reveal a complex temperature dependence of the velocity of the crystal–liquid interface featuring an increase up to a maximum at 10–30% undercooling below the equilibrium melting temperature followed by a gradual decrease of the velocity at deeper levels of undercooling. At the qualitative level, such non-monotonous behaviour of the crystallization front velocity is consistent with the diffusion-controlled crystallization process described by the Wilson–Frenkel model, where the almost linear increase of the interface velocity in the vicinity of melting temperature is defined by the growth of the thermodynamic driving force for the phase transformation, while the decrease in atomic mobility with further increase of the undercooling drives the velocity through the maximum and into a gradual decrease at lower temperatures. At the quantitative level, however, the diffusional model fails to describe the results of MD simulations in the whole range of temperatures with a single set of parameters for some of the model materials. The limited ability of the existing theoretical models to adequately describe the MD results is illustrated in the present work for two materials, chromium and silicon. It is also demonstrated that the MD results can be well described by the solution following from the hodograph equation, previously found from the kinetic phase-field model (kinetic PFM) in the sharp interface limit. The ability of the hodograph equation to describe the predictions of MD simulation in the whole range of temperatures is related to the introduction of slow (phase field) and fast (gradient flow) variables into the original kinetic PFM from which the hodograph equation is obtained. The slow phase-field variable is responsible for the description of data at small undercoolings and the fast gradient flow variable accounts for local non-equilibrium effects at high undercoolings. The introduction of these two types of variables makes the solution of the hodograph equation sufficiently flexible for a reliable description of all nonlinearities of the kinetic curves predicted in MD simulations of Cr and Si. This article is part of the theme issue ‘Transport phenomena in complex systems (part 1)’.

Science & Technology - Other Topics↗

DOE OSTI2021

Computational Fermi Level Engineering and Doping-Type Conversion of Ga2O3 Via Three-Step Processing

Ga2O3 is being actively explored for power electronics, deep-ultraviolet optoelectronics, and other applications due to its ultra-wide bandgap and low projected fabrication cost of large-size and high-quality crystals. N-type doping of Ga2O3 can be achieved and tuned, but p-type doping faces fundamental obstacles due to deep character of acceptor levels and polaron transport of resulting holes. However, successful engineering of Ga2O3 based devices requires critical control of doping density, Fermi level position, and free carrier concentration, providing opportunities for predictive process simulation. We use first-principles defect theory and defect equilibrium calculations to simulate a 3-step growth-annealing-quench protocol for hydrogen assisted Mg doping in Ga2O3, taking into account the hydrogen-oxygen-water gas phase equilibrium. We predict type conversion to a net p-type regime following O-rich annealing after growth under reducing conditions in the presence of H2. This process is similar to the Mg acceptor activation by H removal in GaN. We show that there is an optimal temperature that maximizes the net acceptor density during the equilibrium annealing step for a given Mg doping level. Quenching of non-equilibrium annealed samples then results in a Fermi level EF below mid-gap down to about EV +1.5 eV, creating a significant number of uncompensated neutral MgGa0 acceptors. The resulting free hole concentration in quenched samples is very low (109 – 1012 cm-3) due to deep energy level of these Mg acceptors, but this type converted Ga2O3 material can create a significant built-in field in a p-n junction with an adjoining n-type material. Additionally, the electron concentration is greatly suppressed in such acceptor-doped Ga2O3, which could enable the use as a current blocking layer to fabricate normally-off (enhancement-mode) vertical Ga2O3 based metal-oxide-semiconductor field effect transistors (MOSFETs), and as a guard ring for edge termination in MOSFETs and Schottky barrier diodes (SBDs) with increased breakdown voltage.

CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSICS,M↗

DOE OSTI2020

First-principles study and experimental characterization of metal incorporation in germanium telluride

Germanium telluride is a well-known phase change material (PCM) used in non-volatile memory cells and radio frequency switches. Controlling the properties of GeTe for improved PCM device performance has sometimes been achieved by doping and/or alloying with metals, often at concentrations greater than 10 at. % and using non-equilibrium methods. Here, since switching PCMs between the low-resistance crystalline and high-resistance amorphous states requires a heating cycle, the stability of metal-incorporated GeTe (Ge 0.5–x M x Te 0.5 ) films is also critical to practical implementation of these materials in electronic and optoelectronic devices. In this work, we use both density-functional theory and experimental characterization methods to probe the solubility and critical properties of Ge 0.5–x M x Te 0.5 films. Using first-principles calculations, we determine the enthalpy of formation for GeTe with 2.08, 4.17, and 6.25 at. % of Cu, Fe, Mn, Mo, and Ti and show trends between the stability of the Ge 0.5–x M x Te 0.5 systems and the atomic position, composition, and distribution of the metal atoms in the GeTe matrix.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Grain growth kinetics of the gamma phase metallic uranium

We report metallic uranium is a leading fuel form for sodium cooled fast reactors as an enabling technology of future nuclear energy systems. Mechanistic understanding of fuel behaviors and kinetics under thermodynamic equilibrium and highly non-equilibrium conditions are essential for evaluating fuel performance. It is important to understand and predict the grain and pore evolutions of metallic fuels under thermal and irradiation conditions. However, very limited data are available on the grain growth kinetics and mechanisms of pure gamma phase uranium. In this paper, the pure gamma uranium pellets with different grain structures were fabricated by combining high-energy ball milling and spark plasma sintering. Isothermal annealing tests were performed to investigate the grain growth behavior of the pure gamma phase uranium with different initial grain sizes. A parabolic relationship in grain growth with time was identified for the submicron-sized (374 nm) sample. In contrast, for the nano-sized (137 nm) sample, the grain growth shows a linear relationship with time. The activation energies of grain growth were determined as 199.5 KJ/mol and 80.6 KJ/mol for nano-sized and submicron-sized grain structures, respectively. For the nano-sized sample, the rate-control step of grain growth is dominated by the triple-junction migration, in which the grain boundary triple junction drags the grain growth, leading to a higher activation energy than the bulk diffusion. The dominating mechanism for the submicron-sized sample is grain boundary diffusion. The mechanistic understanding and critical data obtained on the kinetics of pure uranium phases will be useful to evaluate fuel behavior under thermodynamic equilibrium conditions and develop a high fidelity model to predict fuel performance.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2026

A modified scheil approach for nucleation-dependent solidification pathways

As-solidified microstructures of near-eutectic alloys often contain multiple primary phases that are not expected from equilibrium phase diagrams. Such microstructures are caused by cooling-rate-dependent solidification pathways, a factor not captured by the Scheil–Gulliver model or variations thereof. Here, we present a model and algorithm that incorporate the critical nucleation undercooling for each solid phase into the Scheil–Gulliver model. We hypothesize that the non-equilibrium microstructure formation is primarily governed by a nucleation-competition mechanism. This mechanism accounts for both stable/metastable phase selection and primary-phase formation within eutectic regions driven by asymmetric nucleation barriers. The model is validated against a hypereutectic Al-Fe alloy, where it successfully reproduces the observed microstructural constituents, revealing the key dependencies of solidification microstructure on nucleation kinetics. Applicability to multicomponent systems is demonstrated through a hypereutectic Al–Fe–Si ternary alloy, where the model successfully predicts divorced eutectic microstructures and the associated oscillatory solidification pathways along univariant lines. As a result, the proposed framework establishes a nucleation-dependent computational approach for interpreting and predicting solidification microstructures.

Alloy design↗

DOE OSTI2024

Flame kinetics at scramjet-engine-relevant conditions: Role of prompt dissociation of weakly-bound radicals

Combustion in high-speed ram-based propulsion engines occurs under distinct thermodynamic conditions of high reactant temperatures (greater than 1000 K) and relatively low pressures (<5 atm). There is a lack of fundamental flame measurements at such conditions that result in adiabatic flame temperatures (T ad ) exceeding 2500 K. In this work, we have measured laminar flame speeds of oxygen-enriched CH 4 /oxidizer mixtures at sub-atmospheric conditions to probe kinetics at high T ad using the isobaric spherically expanding flame approach. Simulations with recent kinetic models revealed increasing differences between data and model predictions with increasing T ad , reaching up to 25 %. Kinetic analyses reveal that at the thermodynamic conditions in these O 2 -enriched flames, i.e., lower pressures and higher T ad , the effects of HCO prompt dissociation are accentuated. In addition to HCO, the prompt dissociations of CH 2 OH and C 2 H 5 are also considered. Here, the prompt dissociations of all three radicals were evaluated and their effects considered in flame speed simulations. Reaction path analysis for the present flames revealed that approximately half of the reaction flux for HCO formation undergoes prompt dissociation to H + CO. Furthermore, these analyses also revealed that the pathways and sensitive reactions are similar between oxygen-enriched fuel/oxidizer mixtures and preheated fuel/air mixtures, if both have similar T ad . Thus, flames of oxygen-enriched mixtures could be a surrogate to probe the flame chemistry of highly preheated mixtures at relatively low pressures that are often encountered in ram-based propulsion engine combustors.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Effect of Mn on eutectic phase equilibria in Al-rich Al-Ce-Ni alloys

Microstructural analysis of additively manufactured (AM) Al-Ce-Ni-Mn alloys has identified phases not predicted from existing ternary liquidus projections in the Al-Ce-Ni system. Because the rapid cooling rate of AM is orders of magnitude above that of traditional casting, it is unclear if these additional phases arose from the non-equilibrium processing conditions of AM, a drastic shift in phase stability in the system due to the addition of 1 wt% Mn, or some combination of these two influences. The phases and microstructure of cast samples of Al-Ce-Ni and Al-Ce-Ni-Mn alloys were characterized for several annealing conditions which revealed the equilibrium phases at different temperatures. Phase analysis confirmed that minute levels of Mn substituted for Ni in the system drastically shifts the liquidus projection in the Al-rich corner of the ternary phase diagram such that the eutectic Al 3 Ni phase is suppressed in favor of the Al 23 Ni 6 (Ce,Mn) 4 phase. Further addition of Mn promotes the formation of Al 20 Mn 2 Ce and Al 10 Mn 2 Ce phases. The phase analysis data was then used to improve the CALPHAD modeling of the liquidus projection and isothermal sections for the Al-rich Al-Ce-Ni-Mn quaternary system. Thermodynamic modeling and experimental analysis on phases in the AM sample of Al-Ce-Ni with Mn confirmed that the phases present are consistent with Mn-containing Al-Ce-Ni cast samples. Here, this investigation demonstrates the potential for using secondary alloying elements to drastically alter phase stability and microstructure in alloy systems.

36 MATERIALS SCIENCE↗

DOE OSTI2024

MaterialsMap: A CALPHAD-based tool to design composition pathways through feasibility map for desired dissimilar materials, demonstrated with resistance spot welding joining of Ag-Al-Cu

Assembly of dissimilar metals can be achieved by different methods, for example, casting, welding, and additive manufacturing (AM). However, undesired phases formed in liquid-phase assembling processes due to solute segregation during solidification diminish mechanical and other properties of the processed parts. In the present work, an open-source software named MaterialsMap, has been developed based on the CALculation of Phase Diagrams (CALPHAD) approach. The primary objective of MaterialsMap is to facilitate the design of an optimal composition pathway for assembling dissimilar alloys with liquid-phases based on the formation of desired and undesired phases along the pathway. In MaterialsMap, equilibrium thermodynamic calculations are used to predict equilibrium phases formed at slow cooling rate, while Scheil-Gulliver simulations are employed to predict non-equilibrium phases formed during rapid cooling. Additionally, by combining these two simulations, MaterialsMap offers a thorough guide for understanding phase formation in various manufacturing processes, assisting users in making informed decisions during material selection and production. As a demonstration of this approach, a compositional pathway was designed from pure Al to pure Cu through Ag using MaterialsMap. The design was experimentally verified using resistance spot welding (RSW).

36 MATERIALS SCIENCE↗

DOE OSTI2025

Design of novel refractory equiatomic multi-principal elemental alloys based on Mo-Nb-Ti system for Gen IV reactor applications

Excellent irradiation damage resistance demonstrated by multi-principal elemental alloys (MPEAs) has sparked significant interest among researchers, prompting exploration into their vast compositional space, to validate their suitability for nuclear applications. A combined approach of thermodynamic and empirical parameters calculations alongside CALPHAD (CALculation of PHAse Diagrams) for phase formation predictions enable high-throughput material selection for sophisticated applications like nuclear, overcoming laborious and time-consuming experiments. Key thermodynamic and empirical parameters for eight novel equiatomic MPEAs, based on seven low thermal neutron cross section refractory elements, for predicting phase formation were calculated, and equilibrium and non-equilibrium simulations in CALPHAD were employed to comprehensively model the systems. Pseudo binary phase diagram simulations showed that Zr, V or equiatomic CrV additions to the base MoNbTi alloy (MoNbTi-Zr, MoNbTi-V and MoNbTi-CrV alloys) favor the formation of isomorphous body-centered cubic (BCC) phase at high temperatures, while Cr, Al, equiatomic ZrV, or equiatomic CrAl additions (MoNbTi-Cr, MoNbTi-Al, MoNbTi-ZrV or MoNbTi-CrAl alloys) limit the solubility of them. Equilibrium CALPHAD simulations at 750 oC were consistent with XRD results on MoNbTi, MoNbTiZr and MoNbTiCr alloys, and partially for others. Notably, elemental segregation observed in the backscattered electron (BSE) scanning electron microscopy (SEM) images of the alloys was accurately simulated through non-equilibrium Scheil solidification calculations in CALPHAD, further verified by experiments. The precipitation of TiCr2 Laves phase in Cr containing MoNbTiCr and MoNbTiCrAl was accurately predicted while discrepancies were noted in MoNbTiCrV. The equilibrium simulations also provided insights into phase compositions at specific temperatures offering a pathway for tailoring the desired microstructure and properties of these systems. Empirical parameters calculations successfully predicted random solid solution in the base MoNbTi alloy, and with an exception in MoNbTiV and MoNbTiAl, predicted intermetallic precipitation in the rest, especially, Laves phase precipitation in Cr containing alloys.

36 - MATERIALS SCIENCE↗

DOE OSTI2025

Pseudo-equilibrium theory for extrinsic doping control of the topological semimetal Cd 3 As 2

The standard approach for predicting defect equilibria from first principles assumes that the solid-state system is initially in a thermodynamic equilibrium with the external atomic reservoirs. This “growth step” is then often followed by a temperature quench in a “pseudo-equilibrium” in which some or all defect concentrations are frozen in until only the Fermi level E F remains to be equilibrated. However, this protocol does not account for the possibility of site exchanges which can create important defect redistributions as long as short-range defect migration is kinetically permissible. To model this redistribution, we developed an approach to solve for the non-equilibrium chemical potentials as a function of temperature while maintaining the overall defect stoichiometry. We then apply this approach to the Dirac semimetal Cd 3 As 2 to model extrinsic doping with group 1/11 and 14 elements. Undoped Cd 3 As 2 exhibits an undesirable mismatch between E F and the Dirac point. This unintentional electron doping originates from intrinsic defects and is difficult to overcome through adjustment of synthesis conditions alone. Employing our pseudo-equilibrium modeling, we identify extrinsic doping strategies for realizing doping-balanced Cd 3 As 2 at the relatively low temperatures accessible in thin-film growth of this material.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Fluctuating hydrodynamics of chiral active fluids

Active materials are characterized by continuous injection of energy at the microscopic level and typically cannot be adequately described by equilibrium thermodynamics. Here we study a class of active fluids in which equilibrium-like properties emerge when fluctuating and activated degrees of freedom are statistically decoupled, such that their mutual information is negligible. We analyse three paradigmatic systems: chiral active fluids composed of spinning frictional particles that are free to translate, oscillating granular gases and active Brownian rollers. In all of these systems, a single effective temperature generated by activity parameterizes both the equation of state and the emergent Boltzmann statistics. The same effective temperature, renormalized by velocity correlations, relates viscosities to steady-state stress fluctuations via a Green–Kubo relation. To rationalize these observations, we develop a theory for the fluctuating hydrodynamics of these non-equilibrium fluids and validate it through large-scale molecular dynamics simulations. Our work sheds light on the microscopic origin of odd viscosities and stress fluctuations characteristic of parity-violating fluids, in which mirror symmetry and detailed balance are broken.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2020

Data-guided Multi-Map variables for ensemble refinement of molecular movies

Driving molecular dynamics simulations with data-guided collective variables offer a promising strategy to recover thermodynamic information from structure-centric experiments. In this study, the three-dimensional electron density of a protein, as it would be determined by cryo-EM or x-ray crystallography, is used to achieve simultaneously free-energy costs of conformational transitions and refined atomic structures. Unlike previous density-driven molecular dynamics methodologies that determine only the best map-model fits, our work employs the recently developed Multi-Map methodology to monitor concerted movements within equilibrium, non-equilibrium, and enhanced sampling simulations. Construction of all-atom ensembles along the chosen values of the Multi-Map variable enables simultaneous estimation of average properties, as well as real-space refinement of the structures contributing to such averages. Using three proteins of increasing size, we demonstrate that biased simulation along the reaction coordinates derived from electron densities can capture conformational transitions between known intermediates. The simulated pathways appear reversible with minimal hysteresis and require only low-resolution density information to guide the transition. The induced transitions also produce estimates for free energy differences that can be directly compared to experimental observables and population distributions. The refined model quality is superior compared to those found in the Protein Data Bank. We find that the best quantitative agreement with experimental free-energy differences is obtained using medium resolution density information coupled to comparatively large structural transitions. Practical considerations for probing the transitions between multiple intermediate density states are also discussed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Multiphase phase-field approach for solid–solid phase transformations via propagating interfacial phase in HMX

In this study, a thermodynamically consistent multiphase phase-field theory has been formulated to describe temperature-induced solid–solid phase transformations via an interfacial phase. This theory satisfies all thermodynamic equilibrium and stability conditions. It allows us to formulate the analytical solutions of the interface profile, energy, width, and velocity for each of the propagating solid–melt and solid–solid interfaces. The formation of the intermediate melt during solid–solid phase transformation has been studied hundreds of degrees below melting temperature in the HMX energetic crystal. The effect of the penalizing term in the phase-field model and two nanoscale material parameters on the appearance and degree of disordering of intermediate melt has been analyzed for the non-equilibrium interface. The developed approach is applicable to phase transitions in geological, pharmaceutical, ferroelectric, colloidal, and superhard materials where phase transformations occur via the intermediate phase.

Roy, Arunabha M. (ORCID:000000034790726X)↗

DOE OSTI2022

Treating random sequential addition via the replica method

While many physical processes are non-equilibrium in nature, the theory and modeling of such phenomena lag behind theoretical treatments of equilibrium systems. The diversity of powerful theoretical tools available to describe equilibrium systems has inspired strategies that map non-equilibrium systems onto equivalent equilibrium analogs so that interrogation with standard statistical mechanical approaches is possible. In this work, we revisit the mapping from the non-equilibrium random sequential addition process onto an equilibrium multi-component mixture via the replica method, allowing for theoretical predictions of non-equilibrium structural quantities. We validate the above approach by comparing the theoretical predictions to numerical simulations of random sequential addition.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Out of time order correlation of the Hubbard model with random local disorder

The out-of-time-order correlator (OTOC) serves as a powerful tool for investigating quantum information spreading and chaos in complex systems. We present a method employing non-equilibrium dynamical mean-field theory and coherent potential approximation combined with diagrammatic perturbation on the Schwinger–Keldysh contour to calculate the OTOC for correlated fermionic systems subjected to both random disorder and electron interaction. Furthermore, our key finding is that random disorder enhances the OTOC decay in the Hubbard model for the metallic phase in the weakly interacting limit. However, the current limitation of our perturbative solver restricts the applicability to weak interaction regimes.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2025

Thermodynamic consistency and fluctuations in mesoscopic stochastic simulations of reactive gas mixtures

It is essential that mesoscopic simulations of reactive systems reproduce the correct statistical distributions at thermodynamic equilibrium. By considering a compressible fluctuating hydrodynamics (FHD) simulation method of ideal gas mixtures undergoing reversible reactions described by the chemical Langevin equations, we show that thermodynamic consistency in reaction rates and the use of instantaneous temperatures for the evaluation of reaction rates is required for fluctuations for the overall system to be correct. We then formulate the required properties of a thermodynamically consistent reaction (TCR) model. As noted in the literature, while reactions are often discussed in terms of forward and reverse rates, these rates should not be modeled independently because they must be compatible with thermodynamic equilibrium for the system. Using a simple TCR model where each chemical species has constant heat capacity, we derive the explicit condition that the forward and reverse reaction rate constants must satisfy in order for the system to be thermodynamically consistent. We perform equilibrium and non-equilibrium simulations of ideal gas mixtures undergoing a reversible dimerization reaction to measure the fluctuational behavior of the system numerically. We confirm that FHD simulations with the TCR model give the correct static structure factor of equilibrium fluctuations. For the statistically steady simulation of a gas mixture between two isothermal walls with different temperatures, we show using the TCR model that the temperature variance agrees with the corresponding thermodynamic-equilibrium temperature variance in the interior of the system, whereas noticeable deviations are present in regions near walls, where chemistry is far from equilibrium.

Polimeno, Matteo [University of California, Merced↗

DOE OSTI2026

Spectral deconvolution without the deconvolution: Extracting temperature from x-ray Thomson scattering spectra without the source-and-instrument function

X-ray Thomson scattering (XRTS) probes the dynamic structure factor of the system, but the measured spectrum is broadened by the combined source-and-instrument function (SIF) of the setup. In order to extract properties such as temperature from an XRTS spectrum, the broadening by the SIF needs to be removed. Recent work [Dornheim et al. Nat. Commun. 13 , 7911 (2022)] has suggested that the SIF may be deconvolved using the two-sided Laplace transform. However, the extracted information can depend strongly on the shape of the input SIF, and the SIF is in practice challenging to measure accurately. Here, we propose an alternative approach: we demonstrate that considering ratios of Laplace-transformed XRTS spectra collected at different scattering angles is equivalent to performing the deconvolution, but without the need for explicit knowledge of the SIF. From these ratios, it is possible to directly extract the temperature from the scattering spectra, when the system is in thermal equilibrium. We find the method to be generally robust to spectral noise and physical differences between the spectrometers, and we explore situations in which the method breaks down. Furthermore, the fact that consistent temperatures can be extracted for systems in thermal equilibrium indicates that non-equilibrium effects could be identified by inconsistent temperatures of a few eV between the ratios of three or more scattering angles.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

Enzyme activities predicted by metabolite concentrations and solvent capacity in the cell

Experimental measurements or computational model predictions of the post-translational regulation of enzymes needed in a metabolic pathway is a difficult problem. Consequently, regulation is mostly known only for well-studied reactions of central metabolism in various model organisms. In this study, we use two approaches to predict enzyme regulation policies and investigate the hypothesis that regulation is driven by the need to maintain the solvent capacity in the cell. The first predictive method uses a statistical thermodynamics and metabolic control theory framework while the second method is performed using a hybrid optimization–reinforcement learning approach. Efficient regulation schemes were learned from experimental data that either agree with theoretical calculations or result in a higher cell fitness using maximum useful work as a metric. As previously hypothesized, regulation is herein shown to control the concentrations of both immediate and downstream product concentrations at physiological levels. Model predictions provide the following two novel general principles: (1) the regulation itself causes the reactions to be much further from equilibrium instead of the common assumption that highly non-equilibrium reactions are the targets for regulation; and (2) the minimal regulation needed to maintain metabolite levels at physiological concentrations maximizes the free energy dissipation rate instead of preserving a specific energy charge. The resulting energy dissipation rate is an emergent property of regulation which may be represented by a high value of the adenylate energy charge. In addition, the predictions demonstrate that the amount of regulation needed can be minimized if it is applied at the beginning or branch point of a pathway, in agreement with common notions. The approach is demonstrated for three pathways in the central metabolism of E. coli (gluconeogenesis, glycolysis-tricarboxylic acid (TCA) and pentose phosphate-TCA) that each require different regulation schemes. It is shown quantitatively that hexokinase, glucose 6-phosphate dehydrogenase and glyceraldehyde phosphate dehydrogenase, all branch points of pathways, play the largest roles in regulating central metabolism.

59 BASIC BIOLOGICAL SCIENCES↗