SEARCH · Search NASA

Results for “Branch and bound”

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 73 records · Page 4

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Uranium and Nickel Partitioning in a Contaminated Riparian Wetland

Uranium (U) and nickel (Ni) released 50 years ago have been immobilized in the Tims Branch wetlands located on the Savannah River Site in the United States. Sediments were collected from seven locations to identify the factors responsible for this attenuation. Ni and U contents in the solids were significantly correlated, suggesting that depositional as opposed to chemical processes contributed to their spatial distribution. Based on sequential extractions, 63 ± 16% of the U was partitioned into the organic fraction, whereas Ni was distributed between several sediment fractions. An inverse pH-organic matter (OM) correlation and positive correlations of OM with total U and organic-bound U/Ni suggest that increased OM preservation and binding to the mineral surfaces were likely responsible for Ni- and especially U-sediment retention (Tims Branch pH = 4.84 ± 0.68). EXAFS analysis indicated the predominance of U(VI) coordinated with clay minerals (~65%), together with ~35% coordinated to either OM (in areas with elevated OM levels) or iron oxides. The desorption-Kd coefficients of U (3972 ± 1370 L/kg) and Ni (30 ± 8 L/kg) indicate that dissolved Ni poses a greater long-term risk than dissolved U for migrating downstream. This study suggests that a delicate balance of geochemical properties controls whether wetlands behave as sinks or sources of contaminants.

63 RADIATION, THERMAL, AND OTHER ENVIRON. POLLUTAN↗

DOE OSTI2025

Search for the lepton-flavor-violating τ − → e ∓ ℓ ± ℓ − decays at Belle II

We present the result of a search for the charged-lepton-flavor violating decays τ − → e ∓ ℓ ± ℓ − , where ℓ is a muon or an electron, using a data sample with an integrated luminosity of 428 fb −1 recorded by the Belle II experiment at the SuperKEKB e + e − collider. The selection of e + e − → τ + τ − events containing a signal candidate is based on an inclusive-tagging reconstruction and on a boosted decision tree to suppress background. Upper limits on the branching fractions between 1.3 and 2.5 × 10 −8 are set at the 90% confidence level. These results are the most stringent bounds to date for four of the modes.

flavour physics↗

DOE OSTI2021

Burning velocities of R-32/O 2 /N 2 mixtures: Experimental measurements and development of a validated detailed chemical kinetic model

This work entails characterizing the flammability of the refrigerant R-32 (CH 2 F 2 ) by both experimental measurements and modeling. Burning velocities S u were measured using a constant-volume spherical-flame method for R-32/O2/N 2 mixtures with O 2 /N 2 ratios ranging from 21% (synthetic air) to 40%, pressures of (1 to 3) bar, and equivalence ratios $\phi$ of (0.8 to 1.3). Based on a critical assessment of available data, and extended by our own calculations, a detailed chemical kinetic model was developed and key reactions determined using reaction path and sensitivity analyses. Initiation and combustion were identified as distinct kinetic regimes and burning velocities were found to be controlled by two primary reactions: unimolecular decomposition of CH 2 F 2 → CHF + HF and the subsequent reaction, CHF + O 2 → CHFO + O, the latter reaction initiating the radical chain propagating and branching by producing O atoms. Sensitive rate constants in the kinetic model were critically adjusted within their uncertainties and current knowledge bounds to best fit the experimental burning velocities. We found that rate constants in the model could be adjusted to match a given experimental S u for specific conditions (O 2 loading, P, T, $\phi$). This, however, then fixes predicted burning velocities for other all conditions within (3 to 4)% if physically realistic rate parameters are maintained. Thus, the entire set of experimental data is fit, not just to particular conditions. Relative random uncertainties in the experimental Su measurements were (4 to 6)%, but assumptions made for thermal radiation lost by the burned gas in the spherical-flame experiments add an additional systematic uncertainty. Systematic differences between the limiting cases of adiabatic (no thermal radiation lost) and optically-thin (all thermal radiation lost) varied significantly with conditions and ranged from (4 to 30)% at high to low velocities, respectively, translating into uncertainties of (2 to 15)% considering the average of two limiting cases. Comparison of experimental and kinetically modeled Su values suggests that the burned gas tends towards the optically-thin limit at the lowest pressures and fuel loadings and toward the adiabatic limit at the highest pressures and loadings. We tested and found support for this conclusion with a detailed analysis as a function of all the conditions (T, P, % O 2 , $\phi$). This behavior appears to transition from optically-thin to adiabatic as the density of the initial fuel increases, which results in increased CO 2 in the burned gas and thus increased absorption of the thermal radiation (consistent with the Beer-Lambert Law). The validated detailed model based on evaluated kinetics is shown to accurately predict burning velocities for R-32 O 2 /N 2 mixtures over a wide range of conditions and provides a reliable basis for extrapolation to other conditions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Improved analysis of arabinoxylan-bound hydroxycinnamate conjugates in grass cell walls

Background: Arabinoxylan in grass cell walls is acylated to varying extents by ferulate and p-coumarate at the 5-hydroxy position of arabinosyl residues branching off the xylan backbone. Some of these hydroxycinnamate units may then become involved in cell wall radical coupling reactions, resulting in ether and other linkages amongst themselves or to monolignols or oligolignols, thereby crosslinking arabinoxylan chains with each other and/or with lignin polymers. This crosslinking is assumed to increase the strength of the cell wall, and impedes the utilization of grass biomass in natural and industrial processes. A method for quantifying the degree of acylation in various grass tissues is, therefore, essential. We sought to reduce the incidence of hydroxycinnamate ester hydrolysis in our recently introduced method by utilizing more anhydrous conditions. Results: The improved methanolysis method minimizes the undesirable ester-cleavage of arabinose from ferulate and p-coumarate esters, and from diferulate dehydrodimers, and produces more methanolysis vs. hydrolysis of xylan-arabinosides, improving the yields of the desired feruloylated and p-coumaroylated methyl arabinosides and their diferulate analogs. Free ferulate and p-coumarate produced by ester-cleavage were reduced by 78% and 68%, respectively, and 21% and 39% more feruloyl and p-coumaroyl methyl arabinosides were detected in the more anhydrous method. The new protocol resulted in an estimated 56% less combined diferulate isomers in which only one acylated arabinosyl unit remained, and 170% more combined diferulate isomers conjugated to two arabinosyl units. Conclusions: Overall, the new protocol for mild acidolysis of grass cell walls is both recovering more ferulate- and p-coumarate-arabinose conjugates from the arabinoxylan and cleaving less of them down to free ferulic acid, p-coumaric acid, and dehydrodiferulates with just one arabinosyl ester. This cleaner method, especially when coupled with the orthogonal method for measuring monolignol hydroxycinnamate conjugates that have been incorporated into lignin, provides an enhanced tool to measure the extent of crosslinking in grass arabinoxylan chains, assisting in identification of useful grasses for biomass applications.

09 BIOMASS FUELS↗

DOE OSTI2023

Barium stars as tracers of s -process nucleosynthesis in AGB stars: II. Using machine learning techniques on 169 stars

Barium (Ba) stars are characterised by an abundance of heavy elements made by the slow neutron capture process (s-process). This peculiar observed signature is due to the mass transfer from a stellar companion, bound in a binary stellar system, to the Ba star observed today. The signature is created when the stellar companion is an asymptotic giant branch (AGB) star. We aim to analyse the abundance pattern of 169 Ba stars using machine learning techniques and the AGB final surface abundances predicted by the FRUITY and Monash stellar models. We developed machine learning algorithms that use the abundance pattern of Ba stars as input to classify the initial mass and metallicity of each Ba star’s companion star using stellar model predictions. We used two algorithms. The first exploits neural networks to recognise patterns, and the second is a nearest-neighbour algorithm that focuses on finding the AGB model that predicts the final surface abundances closest to the observed Ba star values. In the second algorithm, we included the error bars and observational uncertainties in order to find the best-fit model. The classification process was based on the abundances of Fe, Rb, Sr, Zr, Ru, Nd, Ce, Sm, and Eu. We selected these elements by systematically removing s-process elements from our AGB model abundance distributions and identifying the elements whose removal had the biggest positive effect on the classification. We excluded Nb, Y, Mo, and La. Our final classification combined the output of both algorithms to identify an initial mass and metallicity range for each Ba star companion. With our analysis tools, we identified the main properties for 166 of the 169 Ba stars in the stellar sample. The classifications based on both stellar sets of AGB final abundances show similar distributions, with an average initial mass of M = 2.23 M ⊙ and 2.34 M ⊙ and an average [Fe/H] = –0.21 and –0.11, respectively. We investigated why the removal of Nb, Y, Mo, and La improves our classification and identified 43 stars for which the exclusion had the biggest effect. We found that these stars have statistically significant and different abundances for these elements compared to the other Ba stars in our sample. We discuss the possible reasons for these differences in the abundance patterns.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

A Multi-Branch Decoder Network Approach to Adaptive Temporal Data Selection and Reconstruction for Big Scientific Simulation Data

A key challenge in scientific simulation is that the simulation outputs often require intensive I/O and storage space to store the results for effective post hoc analysis. This article focuses on a quality-aware adaptive temporal data selection and reconstruction problem where the goal is to adaptively select simulation data samples at certain key timesteps in situ and reconstruct the discarded samples with quality assurance during post hoc analysis. This problem is motivated by the limitation of current solutions that a significant amount of simulation data samples are either discarded or aggregated during the sampling process, leading to inaccurate modeling of the simulated phenomena. Two unique challenges exist: 1) the sampling decisions have to be made in situ and adapted to the dynamics of the complex scientific simulation data; 2) the reconstruction error must be strictly bounded to meet the application requirement. To address the above challenges, we develop DeepSample , an error-controlled convolutional neural network framework, that jointly integrates a set of coherent multi-branch deep decoders to effectively reconstruct the simulation data with rigorous quality assurance. The results on two real-world scientific simulation applications show that DeepSample significantly outperforms other state-of-the-art methods on both sampling efficiency and reconstructed simulation data quality.

Zhang, Yang↗

DOE OSTI2022

Joint routing of conventional and range-extended electric vehicles in a large metropolitan network

Range-extended electric vehicles combine the higher efficiency and environmental benefits of battery-powered electric motors with the longer mileage and autonomy of conventional internal combustion engines. This combination is particularly advantageous for time-constrained delivery routing in dense urban areas, where battery recharging along routes can be too time-consuming to economically justify the use of all-electric vehicles. However, switching from electric to conventional fossil fuel modes also results in higher costs and emissions and lower efficiency. This paper analyzes this heterogeneous vehicle routing problem and describes two solution methods: an exact branch-price-and-cut algorithm and an iterated tabu search metaheuristic. From a methodological perspective, we find that the exact algorithm consistently obtains tight lower bounds that also serve to certify the metaheuristic solutions as near-optimal. From a policy standpoint, we examine a large-scale real-world case study concerning parcel deliveries in the Chicago metropolitan area and quantify various operational metrics including energy costs and vehicle miles traveled. We find that by deploying roughly 20% of range -extended vehicles with a modest all-electric range of 33 miles, parcel distributors can save energy costs by up to 17% while incurring less than 0.5% increase in vehicle miles traveled. Increasing the range to 60 miles further reduces costs by only 4%, which can alternatively be achieved by decreasing the average service time by 1 minute or increasing driver working time by 1 hour. In conclusion, our study reveals several key areas of improvement on which vehicle manufacturers, distributors, and policy makers can focus their attention.

33 ADVANCED PROPULSION SYSTEMS↗

DOE OSTI2020

Wave function analysis with a maximum flow algorithm

An efficient algorithm for computing the maximum-flow path in a network is applied to the identification of the dominant configuration state functions (CSFs) in a graphically contracted function (GCF), configuration interaction, wave function. The flow network is a space of spin-adapted CSFs represented by a Shavitt graph, wherein the nodes correspond to orbital occupations and spin quantum numbers. The graph nodes are connected by arcs, and an arc density is defined as sums of the associated squared CSF coefficients. A max-min approach determines an upper bound to the maximum possible incoming flow for each graph node. A backtracking step generates a candidate walk and is followed by a limited search of alternative branching paths for the dominant CSF. The arc density contributions are removed from the graph, and the algorithm is reapplied to the updated graph. This list of generated walks can be partitioned in order to guarantee that the dominant CSFs have been identified. All of the steps in this algorithm are computationally efficient and do not depend on the potentially large dimension of the underlying linear CSF expansion space. An analysis of low-lying valence states of C-2 illustrates the method.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2024

Search for the decay B s 0 → J / ψ π 0 at Belle experiment

We have analyzed 121.4 fb − 1 of data collected at the ϒ ( 5 S ) resonance by the Belle experiment using the KEKB asymmetric-energy e + e − collider to search for the decay B s 0 → J / ψ π 0 . We observe no signal and report an upper limit on the branching fraction B ( B s 0 → J / ψ π 0 ) of 1.21 × 10 − 5 at 90% confidence level. This result is the most stringent, improving the previous bound by 2 orders of magnitude. Published by the American Physical Society 2024

Astronomy & Astrophysics↗

DOE OSTI2021

HMGA2 contributes to vascular development and sprouting angiogenesis by promoting IGFBP2 production

Highlights: • Our findings demonstrates that HMGA2 promotes central arteries formation in zebrafish. • HMGA2 induces vascular development through modulation of IGFBP2 whatever in vivo or in vitro. • The mRNA expressions of HMGA2 and IGFBP2 increased significantly in peripheral blood of Hemangioma patients. • HMGA2-IGFBP2 are promising candidates towards therapeutical vascularization in vascular diseases like Hemangioma. Angiogenesis is the process by which new blood vessels form from preexisting vessels and regulates the processes of embryonic development, wound healing and tumorigenesis. HMGA2 is involved in the occurrence of several cancers, but its biological role and the exact downstream genes involved in vascular development and sprouting angiogenesis remain largely unknown. Here, we first found that HMGA2 knockdown in zebrafish embryos resulted in defects of central artery formation. RNA sequencing revealed that IGFBP2 was significantly downregulated by interference with HMGA2, and IGFBP2 overexpression reversed the inhibition of brain vascular development caused by HMGA2 deficiency. In vitro, we further found that HMGA2 knockdown blocked the migration, tube formation and branching of HUVECs. Similarly, IGFBP2 protein overexpression attenuated the impairments induced by HMGA2 deficiency. Moreover, the promotion of angiogenesis by HMGA2 overexpression was verified in a Matrigel plug assay. We next found that HMGA2 bound directly to a region in the IGFBP2 promoter and positively regulated IGFBP2 expression. Interestingly, the mRNA expression levels of HMGA2 and IGFBP2 were increased significantly in the peripheral blood of hemangioma patients, indicating that overexpression of HMGA2 and IGFBP2 results in vessel formation, consistent with the results of the in vivo and in vitro experiments. In summary, our findings demonstrate that HMGA2 promotes central artery formation by modulating angiogenesis via IGFBP2 induction.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2020

Turbulence/flame/wall interactions in non-premixed inclined slot-jet flames impinging at a wall using direct numerical simulation

In this paper, three-dimensional turbulent non-premixed oblique slot-jet flames impinging at a wall were investigated using direct numerical simulation (DNS). Two cases are considered with the Damköhler number (Da) of case A being twice that of case B. A 17 species and 73-step mechanism for methane combustion was employed in the simulations. It was found that flame extinction in case B is more prominent compared to case A. Reignition in the lower branch of combustion for case A occurs when the scalar dissipation rate relaxes, while no reignition occurs in the lower branch for case B due to excessive scalar dissipation rate. A method was proposed to identify the flame quenching edges of turbulent non-premixed flames in wall-bounded flows based on the intersections of mixture fraction and OH mass fraction iso-surfaces. The flame/wall interactions were examined in terms of the quenching distance and the wall heat flux along the quenching edges. There is essentially no flame/wall interaction in case B due to the extinction caused by excessive turbulent mixing. In contrast, significant interactions between flames and the wall are observed in case A. The quenching distance is found to be negatively correlated with wall heat flux as previously reported in turbulent premixed flames. The influence of chemical reactions and wall on flow topologies was identified. The FS/U and FC/U topologies are found near flame edges, and the NNN/U topology appears when reignition occurs. The vortex-dominant topologies, FC/U and FS/S, play an increasingly important role as the jet turbulence develops.

42 ENGINEERING↗

DOE OSTI2022

Late day measurement of excised branches results in uncertainty in the estimation of two stomatal parameters derived from response curves in Populus deltoides Bartr. x Populus nigra L.

Many terrestrial biosphere models depend on an understanding of the relationship between stomatal conductance and photosynthesis. However, unlike the measurement of photosynthetic parameters, such as the maximum carboxylation capacity, where standard methods (e.g., CO 2 response or ACi curves) are widely accepted, a consensus method for empirically measuring parameters representing stomatal response has not yet emerged. Most models of stomatal response to environment represent stomatal conductance as being bounded by a lower intercept parameter (g 0 ), and linearly scaled based on a multivariate term described by the stomatal slope parameter (g 1 ). Here we employ the widely used Unified Stomatal Optimization model, to test whether g 1 and g 0 parameters are impacted by the choice of measurement method, either on an intact branch, or a cut branch segment stored in water. We measured paired stomatal response curves on intact and excised branches of a hybrid poplar clone (Populus deltoides Bartr. x Populus nigra L. OP367), measured twice over a diurnal period. We found that predawn branch excision did not significantly affect measured g 0 and g 1 when measured within 4 h of excision. Measurement in the afternoon resulted in significantly higher values of g 1 and lower values of g0, with values changing by 55% and 56% respectively. Excision combined with afternoon measurement resulted in a marked effect on parameter estimates, with g1 increasing 89% from morning to afternoon and a 25% lower g 1 for cut branches than those measured in situ. We also show that in hybrid poplar the differences in parameter estimates obtained from plants measured under different conditions can directly impact models of canopy function, reducing modeled transpiration by 18% over a simulated 12.5-hour period. While these findings are only for a single isohydric woody species, our findings suggest that stomatal optimality parameters may not remain constant throughout the day.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Binders for silicon electrodes in lithium-ion batteries

An electrode for a lithium-ion electrochemical cell comprises silicon particles and carbon particles coated on a conductive current collector. The silicon and carbon particles being bound to each other and to the current collector by a cross-linked binder formed from a combination of a poly(carboxylic acid) such as poly(acrylic acid) and a branched polyethyleneimine. A method of preparing the anode also is described.

Shi, Zhangxing↗

DOE OSTI2026

Extending Symmetry-Preserving Attention Networks (SPANet) for Jet Assignment in Fully Hadronic $t\bar{t}$ Events

Fully hadronic $t\bar{t}$ reconstruction requires assigning reconstructed jets to the two top-quark decay branches, $t\bar{t}\to(Wb)(Wb)\to(qqb)(qqb)$. This is a combinatorially large, symmetry-rich, and frequently underconstrained set-assignment problem due to detector effects and partial reconstructability. SPANet provides a strong symmetry-preserving baseline for this task, but its standard inference eagerly commits to a single hypothesis and does not explicitly treat branch-level reconstructability decisions. We present a two-stage extension of the SPANet pipeline for fully hadronic $t\bar{t}$ events: (i) SPANet is modified into a proposal model that generates a shortlist of candidate jet assignments; (ii) a custom dual-head set transformer is trained to re-rank these candidates and identify reconstructible branches. Our extension gives modest improvements in candidate-selection efficiency, with SPANet achieving $69.9\%$ and our extension achieving $71.1\%$. Finally, an oracle study of the generated shortlist shows a substantial upper bound of $88.5\%$ on the evaluated fully reconstructable subset, suggesting that the shortlist contains significant residual information that is not fully exploited.

Lisondodi Tada, Mateo [Puerto Rico U., Mayaguez] (↗

DOE OSTI2024

Search for Baryon Number Violation in Top Quark Production and Decay Using Proton-Proton Collisions at $\sqrt{𝑠}$ = 13 TeV

A search is presented for baryon number violating interactions in top quark production and decay. The analysis uses data from proton-proton collisions at a center-of-mass energy of 13 TeV, collected with the CMS detector at the LHC with an integrated luminosity of 138 fb −1 . Candidate events are selected by requiring two oppositely charged leptons (electrons or muons) and exactly one jet identified as originating from a bottom quark. Multivariate discriminants are used to separate the signal from the background. No significant deviation from the standard model prediction is observed. Upper limits are placed on the strength of baryon number violating couplings. For the first time the production of single top quarks via baryon number violating interactions is studied. This allows the search to set the most stringent constraints to date on the branching fraction of the top quark decay to a lepton, an up-type quark (𝑢 or 𝑐), and a down-type quark (𝑑, 𝑠, or 𝑏). The results improve the previous bounds by 3 to 6 orders of magnitude based on the fermion flavor combination of the baryon number violating interactions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2021

3d N=4 Bootstrap and Mirror Symmetry

We investigate the non-BPS realm of 3d {N} = 4 N = 4 superconformal field theory by uniting the non-perturbative methods of the conformal bootstrap and supersymmetric localization, and utilizing special features of 3d {N} = 4 N = 4 theories such as mirror symmetry and a protected sector described by topological quantum mechanics (TQM). Supersymmetric localization allows for the exact determination of the conformal and flavor central charges, and the latter can be fed into the mini-bootstrap of the TQM to solve for a subset of the OPE data. We examine the implications of the Z_2 Z 2 mirror action for the SCFT single- and mixed-branch crossing equations for the moment map operators, and apply numerical bootstrap to obtain universal constraints on OPE data for given flavor symmetry groups. A key ingredient in applying the bootstrap analysis is the determination of the mixed-branch superconformal blocks. Among other results, we show that the simplest known self-mirror theory with SU(2) \times SU(2) S U ( 2 ) × S U ( 2 ) flavor symmetry saturates our bootstrap bounds, which allows us to extract the non-BPS data and examine the self-mirror Z_2 Z 2 symmetry thereof.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2021

25 Siβ + -decay spectroscopy

Background: β-decay spectroscopy provides valuable information on exotic nuclei and a stringent test for nuclear theories beyond the stability line. Purpose: To search for new β-delayed protons and γ rays of 25 Si to investigate the properties of 25 Al excited states. Method: 25 Si β decays were measured by using the Gaseous Detector with Germanium Tagging system at the National Superconducting Cyclotron Laboratory. The protons and γ rays emitted in the decay were detected simultaneously. Here, a Monte Carlo method was used to model the Doppler broadening of 24 Mg γ-ray lines caused by nuclear recoil from proton emission. Shell-model calculations using two newly developed universal sd-shell Hamiltonians were performed. Results: The most precise 25Si half-life to date has been determined. A new proton branch at 724(4) keV and new proton-γ-ray coincidences have been identified. Three 24 Mg γ-ray lines and eight 25 Al γ-ray lines are observed for the first time in 25 Si decay. The first measurement of the 25 Si β-delayed γ-ray intensities through the 25 Al unbound states is reported. All the bound states of 25 Al are observed to be populated in the β decay of 25 Si. Several inconsistencies between the previous measurements have been resolved, and new information on the 25 Al level scheme is provided. An enhanced decay scheme has been constructed and compared to the mirror decay of 25 Na and the shell-model calculations. Conclusions: The measured excitation energies, γ-ray and proton branchings, log $ft$ values, and Gamow-Teller transition strengths for the states of 25 Al populated in the β decay of 25 Si are in good agreement with the shell model calculations, offering gratifyingly consistent insights into the fine nuclear structure of 25 Al.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2022

HȮ 2 + HȮ 2 : High level theory and the role of singlet channels

At low temperatures and high pressures, the HȮ 2 + HȮ 2 reaction is an important reaction in combustion. There are significant unresolved discrepancies between prior theoretical and experimental studies of this reaction. It has generally been presumed to occur as an abstraction on the triplet surface to produce H 2 O 2 + 3 O 2 . Here, we employ a combination of high-level electronic structure theory (the ANL0 composite method or multi-reference methods as appropriate), sophisticated transition state theory (vibrationally adiabatic torsions, variational, and variable reaction coordinate), and master equation analyses to predict the thermal kinetics on the H 2 O 4 surface. Notably, this analysis suggests a significant branching to HȮ 3 + ȮH near 1000 K via reaction on the singlet surface. This channel does not appear to have been considered in prior combustion models. HȮ 3 itself is a metastable complex (bound by only 3 kcal mol -1 ) that rapidly dissociates to ȮH + 3 O 2 . Thus, the net reaction for this channel, HȮ 2 + HȮ 2 → ȮH + ȮH + O 2 , converts two low reactivity HȮ 2 radicals into two highly reactive ȮH radicals. The ramifications of these newly derived rate expressions are highlighted through kinetic modeling studies of H 2 , CH 3 OH, $\textit{n}$-heptane, and isooctane oxidation; all at the high pressures of relevance to combustion devices.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗