Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “inversion theory”

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 469 records · Page 26
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Neutron-unbound states in 31 Ne

Background: The Island of Inversion near the N = 20 shell gap is home to nuclei with reordered single-particle energy levels compared to the spherical shell model. Studies of 31 Ne have revealed that its ground state has a halo component, characterized by a valence neutron orbiting a deformed 30 Ne core. This lightly-bound nucleus with a separation energy of only S n = 170 keV is expected to have excited states that are neutron unbound. Purpose: The purpose of this experiment was to investigate the low-lying excited states in 31 Ne that decay by the emission of a single neutron. Methods: An 89 MeV/u 33 Mg beam impinged on a segmented Be reaction target. Neutron unbound states in 31 Ne were populated via a two-proton knockout reaction. The 30 Ne fragment and associated neutron from the decay of 31 Ne were detected by the MoNA-LISA-Sweeper experimental setup at the National Superconducting Cyclotron Laboratory. Invariant mass spectroscopy was used to reconstruct the two-body decay energy ( 30 Ne + n). Results: The two-body decay energy spectrum exhibits two features: a low-lying peak at 0.30 ± 0.17 MeV and a broad enhancement at 1.50 ± 0.33 MeV, each fit with an energy-dependent asymmetric Breit-Wigner line shape representing a resonance in the continuum. Accompanying shell model calculations using the FSU interaction within NuShellX, combined with cross-section cal culations using the eikonal reaction theory, indicate that these peaks in the decay energy spectrum are caused by multiple resonant states in 31 Ne. Conclusions: In this study, excited states in 31 Ne were observed for the first time. Transitions from calculated shell model final states in 31 Ne to bound states in 30 Ne are in good agreement with the measured decay energy spectrum. Cross-section calculations for the two-proton knockout populating 31 Ne states as well as spectroscopic factors pertaining to the decay of 31 Ne into 30 Ne are used to examine the results within the context of the shell model expectations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Inversion of dislocation-impurity interactions in α -Fe under magnetic state changes

Here, in this work, we investigate the dislocation-impurity interaction energies and their profiles for various 3d elements —V, Cr, Mn, Cu, Ni, and Co —in and around 1/2 $\langle111\rangle$ screw dislocations in α -Fe using ab initio methods. We consider the ferromagnetic and paramagnetic states, with the latter being modeled through both the disordered local moment model and a spin-wave approach. Our findings reveal that (1) magnetic effects are large compared to size misfit effects of substitutional impurities, and (2) dislocation-impurity interactions are dependent on the magnetic state of the matrix and thermal lattice expansion. In particular, Cu changes from core-attractive in the ferromagnetic state to repulsive in the paramagnetic state.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Observable proton decay in flipped SU(5)

We explore proton decay in a class of realistic supersymmetric flipped SU(5) models supplemented by a U(1) R symmetry which plays an essential role in implementing hybrid inflation. Two distinct neutrino mass models, based on inverse seesaw and type I seesaw, are identified, with the latter arising from the breaking of U(1) R by nonrenormalizable superpotential terms. Depending on the neutrino mass model an appropriate set of intermediate scale color triplets from the Higgs superfields play a key role in proton decay channels that include p → (e + , μ + ) π 0 , p → (e + , μ + ) K 0 , p → $\bar v$π + , and p → $\bar v$K + . We identify regions of the parameter space that yield proton lifetime estimates which are testable at Hyper-Kamiokande and other next generation experiments. We discuss how gauge coupling unification in the presence of intermediate scale particles is realized, and a Z 4 symmetry is utilized to show how such intermediate scales can arise in flipped SU(5). Finally, we compare our predictions for proton decay with previous work based on SU(5) and flipped SU(5).

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Machine-learning Kohn–Sham potential from dynamics in time-dependent Kohn–Sham systems

Abstract The construction of a better exchange-correlation potential in time-dependent density functional theory (TDDFT) can improve the accuracy of TDDFT calculations and provide more accurate predictions of the properties of many-electron systems. Here, we propose a machine learning method to develop the energy functional and the Kohn–Sham potential of a time-dependent Kohn–Sham (TDKS) system is proposed. The method is based on the dynamics of the Kohn–Sham system and does not require any data on the exact Kohn–Sham potential for training the model. We demonstrate the results of our method with a 1D harmonic oscillator example and a 1D two-electron example. We show that the machine-learned Kohn–Sham potential matches the exact Kohn–Sham potential in the absence of memory effect. Our method can still capture the dynamics of the Kohn–Sham system in the presence of memory effects. The machine learning method developed in this article provides insight into making better approximations of the energy functional and the Kohn–Sham potential in the TDKS system.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Fusion burn-propagation simulations using the collisional and radiative particle-in-cell code TRIFORCE

The ability to accurately model burn propagation in inertial confinement fusion plasmas is crucial for advancing fusion energy research. This work presents enhancements to the triforce hybrid fluid-kinetic multiphysics code, focusing on its kinetic half, which employs the particle-in-cell (PIC) method with Monte Carlo collisions (MCC). We use a moment-preserving collision model that mitigates numerical noise, particularly in spherical geometries where particle weights vary significantly. Additionally, we refine the treatment of inverse bremsstrahlung to account for electron–ion collision frequency reductions in degenerate plasmas and incorporate a blackbody radiation source to enable realistic photon injection. These improvements enable the simulation of 1-dimensional (1D) spherical fusion burn propagation in deuterium–tritium plasmas. Benchmark comparisons with the hydra radiation-hydrodynamics code confirm that triforce accurately captures the dynamics of hot-spot expansion and burn propagation, demonstrating sensitivity to ignition thresholds consistent with theoretical models. Findings show the ignition cliff to be less steep in our work compared to radiation-hydrodynamic modeling. These results highlight the role of kinetic effects in fusion ignition physics and underscore the necessity of hybrid fluid-kinetic models for advancing predictive capabilities in high-energy-density plasma systems.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Possible coexistence of magnetism and paramagnetic singularity in lightly Fe-doped WTe 2

Topological semimetals possess nodal or nodal-line phases where conduction and valence bands touch at points or lines in momentum space, respectively. Such band touching is symmetry protected and gives rise to exotic and interesting electronic properties. Coupling topological order with magnetism provides a platform for exploring time-reversal (TR) symmetry breaking topological physics, such as axion electrodynamics, inverse spin-galvanic effect, and the quantized anomalous Hall effect. The Weyl semimetal (WSM) requires breaking either TR symmetry or lattice inversion symmetry (I). By doping inversion-symmetry-broken WSM with magnetic dopants, one can expect to create a WSM with both symmetries breaking simultaneously. Here, structural, electrical, and magnetic properties of Fe x ⁢W 1–x Te 2 (x = 0 and 0.011) are reported. It is revealed that, with a small Fe doping concentration (x = 0.011), a ferromagnetism is induced at low temperature (<10 K). Scanning tunneling microscopy and spectroscopy measurements in Fe 0.011 ⁢W 0.989 ⁢Te 2 further reveal only substitution and no intercalated dopants being observed. The probabilities of the Fe substitutions at the two nonequivalent W sites are quantified with equal probability. The dl/dV point spectra indicates that the Fe substitution in WTe 2 manifests itself as electron doping regardless of doping sites. The results clearly reveal the possible coexistence of magnetism and Weyl points in the lightly Fe doped WTe2 at low temperature. Furthermore, this provides an ideal system for further study on the interplay between the topological Weyl points and the TR symmetry breaking.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

A closer look at how symmetry constraints and the spin–orbit coupling shape the electronic structure of Bi(111)

Abstract Fully relativistic density-functional-theory calculations of Bi(111) thin films are analyzed to revisit their two metallic surface-states branches. We first contrast these metallic branches with surface states arising at gaps in the valence band opened by the spin–orbit coupling (SOC). We find that the two metallic branches along Γ M ‾ do not overlap with the bulk band at the zone boundary,M. We show that the spin texture observed in such states cannot be traced to the lifting of Kramers’ degeneracy. Instead, we track them to the m j = ± 1 / 2 – m j = ± 3 / 2 SOC splitting, the potential anisotropy for in-plane and out-of-plane states, and the coupling between the opposite surfaces of a slab occurring nearM, which is driven by a spatial redistribution of the four metallic states composing the two metallic branches. Each of these branches appears to be non-degenerate at the tested surface, yet each is degenerate with another state of opposite spin at the other surface. Nevertheless, the four metallic states bear some contribution on both surfaces of the film because of their spatial redistribution nearM. The overlapping among these states nearM, afforded by their spatial redistribution on both surfaces, causes a hybridization that perpetuates the splitting between the two branches, makes the film’s electronic structure thickness dependent nearM, extinguishes the magnetic moment of the metallic states avoiding the magnetic-moment discontinuity atM, and denies the need or expectancy of the metallic branches becoming degenerate atM. We propose that theoppositespin polarization observed for the two metallic branches occurs because the surface atoms retain their covalent bonds and thus cannot afford magnetic polarization. We show that the Rashba-splitting of the metallic states for inversion-asymmetric films does not have a fixed magnitude but can be tuned by changing the perturbation breaking inversion symmetry.

Physics↗
DOE OSTI2026

Importance of finite-size corrections for accurate ab initio modeling of carrier capture at semiconductor defects: A case study of substitutional C N in GaN

In ab initio studies of carrier-capture processes in defective semiconductor materials, the single-effective-mode formalism and the static-coupling approximation have become the predominant theoretical approaches for determining carrier-capture coefficients. The single-mode formalism relies on accurate nonequilibrium defect energies obtained from density-functional theory (DFT), where required inputs are a series of configurationally displaced, defect-containing supercells obtained using an interpolative ansatz, and where the DFT outputs are corresponding total energies that have traditionally been postprocessed using a long-established ground-state formulation of finite-size corrections and defect-formation energies. This formulation remains commonly used even though the defects that form a configuration-coordinate (CC) diagram typically exist as structures that are displaced from the ground state. To remedy this inconsistency, Kumagai has recently proposed novel methods for implementing finite-size corrections specifically intended for DFT calculations of the defect energies used to construct CC diagrams and implement the single-mode formalism [Y. Kumagai, Phys. Rev. B 107, L220101 (2023)]. Kumagai's approach builds on the latest finite-size-correction methods introduced to describe vertical charge-state transitions for charge-localizing point defects in semiconductors and insulators [T. Gake et al., Phys. Rev. B 101, 020102 (2020); S. Falletta et al., Phys. Rev. B 102, 041115 (2020)]. The newly identified finite-size artifact treated in these studies is the polarization charge induced on a configurationally frozen defect and its subsequent interaction with a vertical transition in charge state. In this work, we evaluate Kumagai's proposed methodology by applying it in a high-precision DFT study of carrier capture by substitutional C N in GaN, a well-characterized and technologically relevant defect and material. We have rigorously calculated C N defect energies across various supercell sizes for each defect configuration and charge state on the hole-capture CC diagram of C N (𝑞=−1), enabling a direct comparison of the slopes of the defect energies versus inverse cell size with those predicted by Kumagai. The most consequential prediction of Kumagai's method is that these slopes distinctly vary as the square of the linear-interpolation parameter used to construct the nonequilibrium defect configurations. Our results quantitatively support this prediction. Moreover, with these new finite-size corrections and multiple-cell-size DFT calculations in place, we find that the classical energy barrier for hole capture by C N (𝑞=−1) in GaN decreases to 0.092–0.127 eV. This finding confirms the recent ≈ 0.1 eV prediction of Reshchikov based on the weak temperature dependence for hole capture observed in photoluminescence experiments [M. A. Reshchikov, J. Appl. Phys. 129, 121101 (2021)]. These results stand in stark contrast to previously calculated barriers of 0.486 and 0.73 eV, which also used the single-mode formalism but were obtained by instead using ground-state-based finite-size corrections. Our reduced classical barrier for capture increases the temperature-dependent hole-capture coefficient of a C N (𝑞=−1) defect by more than two to four orders of magnitude for temperatures of 100–600 K, compared to the previous 0.486 eV results. While other defects may not be as dramatically affected as here, we suggest that incorporating proper finite-size corrections for the vertical-transition-like states embedded within CC diagrams is an essential, yet previously unrecognized, component of accurate modeling of carrier-capture when using the single-effective-mode formalism.

dielectric properties↗
DOE OSTI2025

Loosely bound composite dark matter

We investigate loosely bound composite states made of dark matter, where the binding energy for constituent particles is considerably less than the constituent mass. We focus on models of nuclear and molecular dark matter, where constituents are separated by length scales larger than the inverse constituent mass, just like nuclei and atoms in the Standard Model. The cosmology, structure, and interactions at underground experiments are described. We find that loosely bound composites can have a very large cross section for scattering with nuclei that scales with nucleon number like ∼ A 4 . For some couplings, these composites produce extremely soft (≪ keV) individual atomic recoils while depositing a large amount of total recoil energy (≫ keV) in a single passage through a detector, implying an interesting new class of signatures for low threshold direct detection.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Shedding light on shadow generalized parton distributions

The feasibility of extracting generalized parton distributions (GPDs) from deeply virtual Compton scattering (DVCS) data has recently been questioned because of the existence of an infinite set of so-called “shadow GPDs” (SGPDs). These SGPDs depend on the process and manifest as multiple solutions (at a fixed scale 𝑄 2 ) to the inverse problem that needs to be solved to infer GPDs from DVCS data. SGPDs therefore pose a significant challenge for extracting GPDs from DVCS data. With this motivation we study the extent to which QCD evolution can provide constraints on SGPDs. This is possible because the known classes of SGPDs begin to contribute to observables after evolution, and can then be constrained (at the input scale 𝑄$^{2}_{0}$) by data that has a finite 𝑄 2 range. The impact that SGPDs could have on determining the total angular momentum, pressure and sheer force distributions, and tomography is also discussed. We find that 𝑄 2 evolution, given empirical data over a wide range of skewness (𝜉 ) and 𝑄 2 scales, constrains the SGPDs studied at low 𝑥 and low 𝜉. If this finding generalizes to all SGPDs—and we have not been able to find any SGPDs with qualitatively different behavior—then this potentially makes extraction of GPDs from DVCS possible over a limited range of 𝑥 and 𝜉.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Constraints on a massive double-copy and applications to massive gravity

We propose and study a BCJ double-copy of massive particles, showing that it is equivalent to a KLT formula with a kernel given by the inverse of a matrix of massive bi-adjoint scalar amplitudes. For models with a uniform non-zero mass spectrum we demonstrate that the resulting double-copy factors on physical poles and that up to at least 5-particle scattering, color-kinematics duality satisfying numerators always exist. For the scattering of 5 or more particles, the procedure generically introduces spurious singularities that must be cancelled by imposing additional constraints. When massive particles are present, color-kinematics duality is not enough to guarantee a physical double-copy. As an example, we apply the formalism to massive Yang-Mills and show that up to 4-particle scattering the double-copy construction generates physical amplitudes of a model of dRGT massive gravity coupled to a dilaton and a two-form with dilaton parity violating couplings. We show that the spurious singularities in the 5-particle double-copy do not cancel in this example, and the construction fails to generate physically sensible amplitudes. We conjecture sufficient constraints on the mass spectrum, which in addition to massive BCJ relations, guarantee the absence of spurious singularities.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

QuGrav: Bringing gravitational waves to light with qumodes

We propose using qumodes, quantum bosonic modes, for detecting high-frequency gravitational waves via the inverse Gertsenshtein effect, where a gravitational wave resonantly converts into a single photon in a magnetized cavity. For an occupation number n of the photon field in a qumode, the conversion probability is enhanced by a factor of n + 1 due to Bose-Einstein statistics. Unlocking this increased sensitivity entails the ability to continuously prepare the qumode and perform nondemolition measurement on the qumode-qubit system within the qumode coherence time. Our results indicate that, at microwave frequencies and with existing technology, the proposed setup can attain sensitivities within 1.7 orders of magnitude of the cosmological bound. With anticipated near-future improvements, it has the potential to surpass this limit and pave the way for the first exploration of high-frequency cosmological gravitational wave backgrounds. At optical frequencies, it can enhance the sensitivity of current detectors by one order of magnitude. That further enhances their potential in reaching the single-graviton level.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Simulation and Emulation of X-Ray Diffraction from Dynamic Compression Experiments

Many important aspects of the dynamic thermo-mechanical response of materials occur at the mesoscale, i.e. a physical scale of interactions smaller than what can be adequately described by homogenous behaviors, yet larger than the scale of the atomic lattice. Concurrent advancements in computational power, continuum theory, and experimental diagnostics are enabling unprecedented understanding of such interactions. However, we cannot develop a sufficient level of confidence in such mesoscale capability until the constitutive description of the underlying constituents is reliably representative of their actual physical behavior. Therefore, there is a strong need to combine experimental, modeling, and data-science techniques to validate models of the thermomechanical response of individual single crystals. One experimental diagnostic with high potential impact to shock physics and materials science is in-situ x-ray diffraction. This paper is primarily focused on simulation of x-ray diffraction in shock physics, but with an aim toward quantifying parametric uncertainty of simulation models. Here, we develop and demonstrate a data-science and model-driven approach to constrain the parameterization of continuum models of crystal lattice deformation associated with the shock response of crystalline materials. The framework is built around the connection between continuum hydrodynamic simulations of lattice deformation and a new Bragg diffraction simulation code, BarberShop. The dynamic deformation of a crystal lattice is modeled using the DiscoFlux model within an arbitrary Lagrangian-Eulerian hydrodynamic code, FLAG. These detailed continuum simulations of lattice deformation can be computationally slow, thus a statistical model is used to emulate the evolution of lattice deformation fields in time and across the considered model parameter space. Emulated lattice deformation fields can then be generated rapidly for any combination of physics model parameters. In turn, these fields can be fed into BarberShop to realize a rapid prediction of Bragg diffraction patterns associated with particular values of physics model parameters. The framework enables parameterization of the single crystal model to obtain Bragg diffraction patterns that most closely resemble a corresponding measurement. Furthermore, the framework naturally provides sensitivities of the lattice deformation to the physics parameters. We highlight the utility of this framework through the application to a synthetic closed-loop inverse problem leading to the parameterization of a single crystal material model. As a model problem, we consider the dynamic response of the energetic molecular crystal, cyclotrimethylenetrinitramine (or RDX), under dynamic compression induced by simulated flyer plate impact experiments.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Observation of Weyl fermions in a magnetic non-centrosymmetric crystal

The absence of inversion symmetry in non-centrosymmetric materials has a fundamental role in the emergence of a vast number of fascinating phenomena, like ferroelectricity, second harmonic generation, and Weyl fermions. The removal of time-reversal symmetry in such systems further extends the variety of observable magneto-electric and topological effects. Here we report the striking topological properties in the non-centrosymmetric spin-orbit magnet PrAlGe by combining spectroscopy and transport measurements. By photoemission spectroscopy below the Curie temperature, we observe topological Fermi arcs that correspond to projected topological charges of ±1 in the surface Brillouin zone. In the bulk, we observe the linear energy-dispersion of the Weyl fermions. We further observe a large anomalous Hall response in our magneto-transport measurements, which is understood to arise from diverging bulk Berry curvature fields associated with the Weyl band structure. These results establish a novel Weyl semimetal phase in magnetic non-centrosymmetric PrAlGe.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Heavy dark matter in rapidly evolving massive stars

We study the impact of heavy dark matter (DM) captured in massive stars via scattering(s) with the star constituents. We focus on the first stars and use stellar evolution simulations to track down how DM capture evolves over time from the zero-age main sequence to the late metal-rich stages of stellar evolution. During the early hydrogen-helium-dominated phase, the capture process is well described by scattering with two targets. As a star evolves, metal production leads to the formation of a dense core surrounded by a lighter envelope. The core significantly enhances the capture of ultra-heavy DM; in this case, three distinct nuclear species are required to accurately describe multiple-scattering capture. We use the Eddington inversion method to obtain a realistic DM velocity distribution, better suited when the star is near the center of a halo, than the widely used Maxwell-Boltzmann distribution. We find that heavy DM would be able to thermalize and achieve capture-annihilation equilibrium within a massive star's lifetime for regions of the parameter space not excluded by direct detection. For non-annihilating DM, because of the high amount of targets available for capture and despite massive stars being short-lived, it would even be possible for DM to achieve self-gravitation and collapse to a black hole, which eventually could swallow the star from within before the expected end of the star's life, for non-excluded regions of the parameter space. Our results highlight the dependence of DM capture on the stellar evolutionary stage, composition, and halo location, demonstrating that accurate modeling of massive stars is essential for constraining heavy DM with primordial stellar populations.

dark matter theory↗
DOE OSTI2023

Third-order Hall effect in the surface states of a topological insulator

Time reversal and inversion symmetric materials fail to yield linear and nonlinear responses since they possess net zero Berry curvature. However, higher-order Hall response can be generated in these systems upon constraining the crystalline symmetries. Motivated by the recently discovered third-order Hall (TOH) response mediated by Berry connection polarizability, namely, the variation of the Berry connection with respect to an applied electric field, here, we investigate the existence of such a Hall effect in the surface states of hexagonal warped topological insulators (e.g., Bi 2 Te 3 ) under the application of only electric field. Using the semiclassical Boltzmann formalism, we investigate the effect of tilt and hexagonal warping on the Berry connection polarizability tensor and, consequently, the TOH effect, provided the Dirac cone remains gapless. We find that the magnitude of the response increases significantly with increasing tilt strength and warping, and therefore, they can provide the tunability of this effect. In addition, we also explore the effect of chemical doping on the TOH response in this system. Interestingly, we show, based on the symmetry analysis, that the TOH can be the leading-order response in this system, which can directly be verified in experiments.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Orbit-averaging and deposition accuracy for runaway electron beams in hybrid kinetic-MHD simulations of the runaway plateau

We develop a new procedure that combines the kinetic orbit runaway electron code (KORC) and the NIMROD extended-magnetohydrodynamic code to simulate runaway electrons (REs) in the post-disruption plateau. KORC integrates guiding-center orbits, with a barycentric-based binary search strategy providing initial guesses for the Newton–Raphson logical-to-physical coordinate inversion, ensuring reliable particle-to-mesh mapping in NIMROD, whose fields remain static for the present study. Samples are drawn in accord with experimental parallel current profiles of RE beams during the plateau phase. Deposition in NIMROD is verified through comparison with a Python-based finite-element code that ensures periodicity in the poloidal direction and continuity at the magnetic axis. Accurate representation of near-axis fields requires finer mesh resolution to prevent under- and overshoots in current density from orbit inaccuracies. Yet, at a fixed particle count, increasing mesh resolution amplifies statistical noise in the deposited fields. An orbit-averaging method accumulates partial current deposits over multiple kinetic steps and reduces the statistical noise with little added computational cost. By coupling kinetic routines from KORC directly into the NIMROD codebase, these developments lay essential groundwork for future self-consistent KORC–NIMROD coupling.

Algorithms and data structure↗
DOE OSTI2025

Degrees of Rate Control in Interconnected Reaction Networks

Overall reactions in interconnected networks exhibit net, forward, and reverse rates that are governed by both constitutive elementary steps in the pathway of interest and branching elementary steps that lead to alternative products. Accordingly, steps in branching pathways exhibit negative net, forward, and reverse degrees of rate control, as they reduce reaction flux to the desired product. We here contextualize the forward and reverse degrees of rate control in terms of kinetic resistances (inverse of rates) and leverage the additive nature of kinetic resistance to decouple kinetic driving forces contributed by constitutive elementary steps and branching points (nodal species) in interconnected networks. Regardless of the network connectivity, forward and reverse degrees of rate control are shown to converge at equilibrium. Away from equilibrium, we identify two critical features of interconnected networks: stoichiometric regularity─condition where all stoichiometric numbers are unity─and pathway symmetry around nodal species─condition where branching pathways share the same rate constants, stoichiometry, and species concentrations/activities─that result in (i) equal forward, reverse, and consequently net degrees of rate control and (ii) forward and reverse degrees of rate control that exhibit constant offsets, respectively, across all extents of reaction. Furthermore, our discourse further provides a mathematical description for the influence of stoichiometric irregularity and pathway asymmetry on forward and reverse degrees of rate control. Altogether, the presented work details the effects of network (inter)connectivity and stoichiometry on reaction kinetics and, in doing so, establishes general protocols for capturing these effects as additive terms in the formulation of forward and reverse degrees of rate control.

10 SYNTHETIC FUELS↗