Search NASA⌕ Search

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.

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

At least 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Discovery and Spectroscopic Characterization of a Distant, Compact Milky Way Satellite in Gemini

We present the discovery of a compact Milky Way satellite in the constellation of Gemini. This system was discovered by cross-matching detections from two independent search algorithms applied to Blanco/DECam data from the third data release of the DECam Local Volume Exploration survey (DELVE DR3), and confirmed with deeper imaging from Gemini/GMOS-N. Based on these data, we determine that the system is an ultra-faint ($M_V = -2.1^{+0.4}_{-0.6}$), compact ($r_{1/2} = 8.6^{+1.4}_{-1.2}$ pc) system located at a heliocentric distance of $120^{+7}_{-6}$ kpc. These physical properties place the system in the regime of ambiguous, ultra-faint compact Milky Way halo satellites that cannot be confidently classified as dwarf galaxies or star clusters from morphology alone; we therefore name the system DELVE 8/Gemini I. From medium-resolution Keck/DEIMOS spectroscopy, we securely identify four members including two blue horizontal branch stars, confirming the system as a bound satellite moving at a mean radial velocity of $v_{\rm hel} = -82.7^{+3.7}_{-3.9} {\rm km\,s}^{-1}$. We also use these spectra to place an upper limit of $\rm [Fe/H] \lesssim -2.5$ on the metallicity of DELVE 8/Gemini I's brightest star, supporting the classification of the system as either an ancient star cluster or ultra-faint dwarf galaxy. The discovery of faint, distant systems similar to DELVE 8/Gemini I is expected to become more common with upcoming surveys.

Overdeck, K. [Chicago U., Astron. Astrophys. Ctr.;↗
DOE OSTI2020

Scenario Grouping and Decomposition Algorithms for Chance-Constrained Programs

A lower bound for a finite-scenario-based chance-constrained program is the quantile value corresponding to the sorted optimal objective values of scenario subproblems. This quantile bound can be improved by grouping subsets of scenarios at the expense of solving larger subproblems. The quality of the bound depends on how the scenarios are grouped. In this paper, we formulate a mixed-integer bilevel program that optimally groups scenarios to tighten the quantile bounds. For general chance-constrained programs, we propose a branch-and-cut algorithm to optimize the bilevel program, and for chance-constrained linear programs, a mixed-integer linear-programming reformulation is derived. Here, we also propose several heuristics for grouping similar or dissimilar scenarios. Our computational results demonstrate that optimal grouping bounds are much tighter than heuristic bounds, resulting in smaller root-node gaps and better performance of scenario decomposition for solving chance-constrained 0-1 programs. Also, the optimal grouping bounds can be greatly strengthened using larger group size.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Uranium Immobility in Tims Branch Wetland

The objective of this study was to address an item identified by the Low-Level Waste Disposal Facility Federal Review Group (LFRG) panel during their review of the Savannah River Site’s Composite Analysis. Their concern was that the Tims Branch wetland, which has elevated levels of uranium (U) originating from site operations between 1952 to 1985, may be a secondary source term, posing a human health risk. This report briefly describes three studies dealing with: 1) quantifying the mass of U remaining in the wetland, 2) quantifying the tendency of U to release from Tims Branch sediment, and 3) describing the distribution of U in aqueous and suspended solids in Tims Branch water during episodic storm events. The Tims Branch wetland has been extremely effective for the last 50 years at immobilizing U released to the system. A detailed map including 270,000 gamma spectra of the U distribution in the wetland indicated that 84% of the initial 43,500 kg U introduced in the system remain there after 50 years (Kaplan et al., 2020). The effectiveness of immobilizing U is supported by stream water measurements collected over the last several decades indicating that U concentrations have been consistently below drinking water regulatory limits (i.e., the EPA Maximum Contaminant Level, 30 pg/L (SRNS, 2019)). Furthermore, using contaminated Tims Branch sediments, in situ desorption Kd values (ratio of U sediment/U Water) were extremely high, ranging from 104to 106 L/kg, suggesting the U is strongly bound to these wetland sediments. Finally, a detailed study of Tims Branch stream water properties during episodic rain events indicated that essentially all the U transported in this stream was associated with suspended solids (Batson et al., 1996). Together these findings indicate that the Tims Branch wetland is extremely effective in its present state at immobilizing U and that the small masses of U released downstream will be predominantly in association with suspended solids.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2026

Limiter-aware adaptive excitation control for hydropower generators with legacy voltage regulation systems

Hydropower units often rely on legacy excitation systems with fixed automatic voltage regulator (AVR) settings, but parameter drift, changing operating conditions, and evolving grid-service demands can degrade terminal-voltage regulation over time. This problem should be addressed because poor voltage regulation affects reactive-power support, system stability, and reliable plant operation, and replacing installed AVR and limiter logic is costly and disruptive in practice. To address this issue, this paper proposes a retrofit-friendly enhancement for hydropower generator excitation systems that combines adaptive tuning of the AVR middle lead–lag branch with a small proportional-integral correction for removing steady bias and slow drift. The adaptive laws are normalized, projection-bounded, and equipped with leakage, while adaptation is enabled only in the middle branch. Thus, the intent of the underexcitation and overexcitation limiters, as well as legacy current and rate limits, remains unchanged. Stability is analyzed under incremental performance assumptions, yielding boundedness and asymptotic tracking in the absence of saturation. The method is evaluated using real operational data from a utility-scale hydropower unit at Rocky Reach Dam, Unit C-8. Compared with the fixed AVR baseline, the proposed method reduces terminal-voltage mean squared error by about 70%, removes long periods of steady-state bias, and improves transient tracking while coexisting with limiter actions. These results show that substantial voltage-control improvement can be achieved on existing hydropower units with minimal integration effort and without replacing the installed AVR structure.

ESST5B model↗
DOE OSTI2023

Minimizing Optimal Transport for Functions with Fixed-Size Nodal Sets

Consider the class of zero-mean functions with fixed L ∞ and L 1 norms and exactly N ϵ N nodal points. Which functions f minimize W p (f + ,f – ), the Wasserstein distance between the measures whose densities are the positive and negative parts? We provide a complete solution to this minimization problem on the line and the circle, which provides sharp constants for previously proven “uncertainty principle”-type inequalities, i.e., lower bounds on N • W p (f + ,f – ). We further show that, while such inequalities hold in many metric measure spaces, they are no longer sharp when the non-branching assumption is violated; indeed, for metric star-graphs, the optimal lower bound on W p (f + ,f – ) is not inversely proportional to the size of the nodal set, N. Here, based on similar reductions, we make connections between the analogous problem of minimizing W p (f + ,f – ) for f defined on Ω C R d with an equivalent optimal domain partition problem.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Y-12 Groundwater Protection Program Monitoring Well Inspection and Maintenance Plan

This plan describes the systematic approach for: inspecting the physical condition of monitoring wells at Y-12, determining maintenance needs that extend the life of a well, and identifying those wells that no longer meet acceptable monitoring well design or well construction standards and require plugging and abandonment. The inspection and maintenance of groundwater monitoring wells is one of the primary management strategies of the Y-12 Groundwater Protection Program (GWPP) Management Plan, that is, the “proactive stewardship of the extensive monitoring well network at Y-12" (Consolidated Nuclear Security, L.L.C. [CNS], 2018). Effective stewardship, and a program of routine inspections of the physical condition of each monitoring well, ensures that representative water-quality samples and hydrologic data are obtained from the well network and protects the subsurface environment. In accordance with the Y-12 GWPP Monitoring Optimization Plan (MOP) for Groundwater Monitoring Wells at the Y-12 National Security Complex, Oak Ridge, Tennessee (CNS, 2017), the status designation (active or inactive) for each well determines the scope and extent of well inspections and maintenance activities. This plan, in conjunction with the above document, formalizes the GWPP approach to focus available resources on monitoring wells which provide the most useful data, and for that reason the GWPP inspects and performs maintenance on the wells sampled by the GWPP. This plan applies to groundwater monitoring wells installed at Y-12 and the related waste management facilities located within the three hydrogeologic regimes: (1) the Bear Creek Hydrogeologic Regime (Bear Creek Regime), (2) the Upper East Fork Poplar Creek Hydrogeologic Regime (East Fork Regime), and (3) the Chestnut Ridge Hydrogeologic Regime (Chestnut Ridge Regime). The Bear Creek Regime encompasses the section of Bear Creek Valley (BCV) immediately west of Y-12. The East Fork Regime encompasses most of the Y-12 process, operations, and support facilities in BCV west of Scarboro Road. The Chestnut Ridge Regime is directly south of Y-12 and encompasses a section of Chestnut Ridge that is bounded to the west by a surface drainage feature (Dunaway Branch, located immediately west of Industrial Landfill II) and by Scarboro Road to the east. The GWPP maintains an extensive database of geographic and construction details and related information for the monitoring wells in each hydrogeologic regime in the Updated Subsurface Database for Bear Creek Valley, Chestnut Ridge, and Parts of Bethel Valley on the U.S. DOE Oak Ridge Reservation (CNS, 2019). A detailed description of the hydrogeologic framework at Y-12 can be found in the GWPP Management Plan (CNS, 2018).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Leptonic probes of alternative left-right symmetric models

We explore constraints on the parameter space of the alternative left-right model originating from the leptonic sector. Our analyses focuses on both lepton-flavor-conserving observables, particularly the anomalous magnetic moment of the muon, and lepton-flavor-violating processes like μ → e γ decay and μ − e conversions in nuclei. While contributions to the anomalous magnetic moment fall below the measured values at 2 σ , current and future experimental sensitivities to flavor-violating branching rations of the Standard Model leptons are expected to impose lower bounds on the mass of the peculiar S U ( 2 ) R gauge boson of the model. This provides complementary constraints relative to existing limits, which are indirect and derived from collider bounds on the mass of the associated neutral gauge boson Z ′ . Published by the American Physical Society 2025

Frank, Mariana (ORCID:0000000322684821)↗
DOE OSTI2020

KOTO vs. NA62 dark scalar searches

The two kaon factories, KOTO and NA62, are at the cutting edge of the intensity frontier, with an unprecedented numbers of long lived and charged Kaons,~10 13 , being measured and analyzed. These experiments have currently a unique opportunity to search for dark sectors. In this paper, we demonstrate that searches done at KOTO and NA62 are complementary, both probing uncharted territories. We consider two qualitatively different physics cases. In the first, we analyze models of axion-like-particles (ALP) which couple to gluons or electroweak gauge bosons. In the second, we introduce a model based on an approximate strange flavor symmetry that leads to a strong violation of the Grossman-Nir bound. For the first scenario, we design a new search strategy for the KOTO experiment,K L → π 0 a →4γ. Its expected sensitivity on the branching ratio is at the level of 10 -9 . This demonstrates the great potential of KOTO as a discovery machine. In addition, we revisit other bounds on ALPs from Kaon factories, highlighting the main sources of theoretical uncertainty, and collider experiments, and show new projections. For the second scenario, we show that the model may be compatible with the preliminary analysis of the KOTO-data that shows a hint for New Physics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Occupational disorder as the origin of flattening of the acoustic phonon branches in the clathrate Ba 8 ⁢Ga 16 ⁢Ge 30

In the search for high-performance thermoelectrics, materials such as clathrates have drawn attention due to having both glasslike low phonon thermal conductivity and crystal-like high electrical conductivity. Ba 8 ⁢Ga 16 ⁢Ge 30 (BGG) has a loosely bound guest Ba atom trapped inside rigid Ga-Ge cage structures. Avoided crossings between acoustic phonons and the flat guest atom branches have been proposed to be the source of the low lattice thermal conductivity of BGG. Ga-Ge site disorder with Ga and Ge exchanging places in different unit cells has also been reported. We used time-of-flight neutron scattering to measure the complete phonon spectrum in a large single crystal of BGG and compared these results with predictions of density functional theory to elucidate the effect of the disorder on heat-carrying phonons. Experimental results agreed much better with the calculation assuming the disorder than with the calculation assuming the ordered configuration. Although the atomic masses of Ga and Ge are nearly identical, we found that disorder strongly reduces phonon group velocities, which significantly reduces thermal conductivity. Finally, our work points to a path towards optimizing thermoelectrics.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Search for lepton-flavor-violating τ – → μ – μ + μ – decays at Belle II

We present the result of a search for the charged-lepton-flavor violating decay τ – → μ – μ + μ – using a 424 fb –1 sample of data recorded by the Belle II experiment at the SuperKEKB e + e – collider. The selection of e + e – → τ + τ – events is based on an inclusive reconstruction of the non-signal tau decay, and on a boosted decision tree to suppress background. We observe one signal candidate, which is compatible with the expectation from background processes. We set a 90% confidence level upper limit of 1.9 × 10 –8 on the branching fraction of the τ – → μ – μ + μ – decay, which is the most stringent bound to date.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

New bounds on heavy QCD axions from big bang nucleosynthesis

We study big bang nucleosynthesis (BBN) constraints on heavy QCD axions. BBN offers a powerful probe of new physics that modifies the neutron-to-proton ratio during the process, thanks to the precisely measured primordial Helium-4 abundance. A heavy QCD axion provides an attractive target for this probe, because not only is it a well-motivated hypothetical particle by the strong 𝐶⁢𝑃 problem, but also it dominantly decays to hadrons if kinematically allowed. A range of its lifetime is thus excluded where the hadronic decays would significantly alter the neutron-to-proton ratio. We compute axion-induced modification of the neutron-to-proton ratio, and obtain robust upper bounds on the axion lifetimes, as low as 0.017 s for the axion mass higher than 300 MeV. Remarkably, this is stronger than projected future cosmic microwave background bounds via 𝑁 eff . Our bounds are largely insensitive to uncertainties in hadronic cross sections and the axion’s branching fractions into various hadrons, as well as to the precise value of the initial axion abundance. We also incorporate, for the first time, several key improvements, such as scattering processes by energetic 𝐾 𝐿 and secondary hadrons, that can also be important for studying general hadronic injections during BBN, not limited to those from axion decays.

Axions↗
DOE OSTI2025

Search for 𝑃 𝑐$\overline{𝑐}$𝑠 ⁢(4459) 0 and 𝑃 𝑐$\overline{𝑐}$𝑠 (4338) 0 in ϒ(1𝑆,2𝑆) Inclusive Decays at Belle

Using data samples of 102 million ϒ⁡(1⁢𝑆) events and 158 million ϒ⁡(2⁢𝑆) events collected by the Belle detector at the KEKB asymmetric-energy 𝑒+⁢𝑒− collider, we search for [𝑢⁢𝑑⁢𝑠⁢𝑐$\overline{𝑐}$ ] pentaquark states decaying to 𝐽/𝜓⁢Λ. Using the first observations of ϒ⁡(1⁢𝑆,2⁢𝑆) inclusive decays to 𝐽/𝜓⁢Λ, we find evidence of the 𝑃 𝑐$\overline{𝑐}$𝑠 ⁢(4459) 0 state with a local significance of 3.3 standard deviations, including statistical and systematic uncertainties. We measure the mass and width of the 𝑃 𝑐$\overline{𝑐}$𝑠 ⁢(4459) 0 to be (4471.7±4.8±0.6) MeV/𝑐 2 and (22 ±13 ±3) MeV, respectively. The branching fractions for 𝑃 𝑐$\overline{𝑐}$𝑠 ⁢(4459) 0 production are measured to be ℬ⁡[ϒ⁡(1⁢𝑆)→𝑃 𝑐$\overline{𝑐}$𝑠 ⁢(4459)0/ $\overline{𝑃}$ 𝑐$\overline{𝑐}$𝑠 ⁢(4459) 0 +anything]=(3.5±2.0±0.2)×10 −6 and ℬ⁡[ϒ⁡(2⁢𝑆)→𝑃 𝑐$\overline{𝑐}$𝑠 ⁢(4459) 0 /$\overline{𝑃}$ 𝑐$\overline{𝑐}$𝑠 ⁢(4459) 0 +anything]=(2.9±1.7±0.4)×10 −6 . The inclusive branching fractions of ϒ⁡(1⁢𝑆,2⁢𝑆) →𝐽/𝜓⁢Λ/$\overline{Λ}$ are measured to be ℬ⁡[ϒ⁡(1⁢𝑆) →𝐽/𝜓⁢Λ/$\overline{Λ}$ +anything] =(36.9 ±5.3 ±2.4) ×10 −6 and ℬ⁡[ϒ⁡(2⁢𝑆) →𝐽/𝜓⁢Λ/$\overline{Λ}$ +anything] =(22.3 ±5.7 ±3.1) ×10 −6 . We measure the visible cross section 𝜎⁡(𝑒+⁢𝑒−→𝐽/𝜓⁢Λ/ $\overline{Λ}$ +anything)=(90±14±6) fb for the continuum production at √𝑠 =10.52 GeV. In all cases, the first uncertainties are statistical and the second are systematic.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

The Inducible Accumulation of Cell Wall-Bound p -Hydroxybenzoates Is Involved in the Regulation of Gravitropic Response of Poplar

Lignin in Populus species is acylated with p-hydroxybenzoate. Monolignol p-hydroxybenzoyltransferase 1 (PHBMT1) mediates p-hydroxybenzoylation of sinapyl alcohol, eventually leading to the modification of syringyl lignin subunits. Angiosperm trees upon gravistimulation undergo the re-orientation of their growth along with the production of specialized secondary xylem, i.e., tension wood (TW), that generates tensile force to pull the inclined stem or leaning branch upward. Sporadic evidence suggests that angiosperm TW contains relatively a high percentage of syringyl lignin and lignin-bound p-hydroxybenzoate. However, whether such lignin modification plays a role in gravitropic response remains unclear. By imposing mechanical bending and/or gravitropic stimuli to the hybrid aspens in the wild type (WT), lignin p-hydroxybenzoate deficient, and p-hydroxybenzoate overproduction plants, we examined the responses of plants to gravitropic/mechanical stress and their cell wall composition changes. We revealed that mechanical bending or gravitropic stimulation not only induced the overproduction of crystalline cellulose fibers and increased the relative abundance of syringyl lignin, but also significantly induced the expression of PHBMT1 and the increased accumulation of p-hydroxybenzoates in TW. Furthermore, we found that although disturbing lignin-bound p-hydroxybenzoate accumulation in the PHBMT1 knockout and overexpression (OE) poplars did not affect the major chemical composition shifts of the cell walls in their TW as occurred in the WT plants, depletion of p-hydroxybenzoates intensified the gravitropic curving of the plantlets in response to gravistimulation, evident with the enhanced stem secant bending angle. By contrast, hyperaccumulation of p-hydroxybenzoates mitigated gravitropic response. These data suggest that PHBMT1-mediated lignin modification is involved in the regulation of poplar gravitropic response and, likely by compromising gravitropism and/or enhancing autotropism, negatively coordinates the action of TW cellulose fibers to control the poplar wood deformation and plant growth.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Tetraquarks in the Born-Oppenheimer approximation

The conventional loosely bound molecule interpretation of the X(3872) is not compatible with the recent LHCb experimental measurement of the ratio of branching fractions $\mathcal{R}$ = Br(X → ψ′γ)/Br(X → ψγ). We systematically determine the entire tetraquark spectrum for J = 0, 1, 2 and refine the calculation of $\mathcal{R}$ in an improved Born-Oppenheimer description of the X(3872) compact tetraquark. This refinement yields a significantly better agreement with experimental data on $\mathcal{R}$ and on the spectroscopy of the states themselves. Extending the diquark-antidiquark paradigm to encompass tetraquarks that are linear super-positions of open charm singlets and color octets, we discover that these exotic resonances manifest as compact shallow bound states of quarks in color force potentials.

Effective Field Theories of QCD↗
DOE OSTI2022

The Euler characteristic as a basis for teaching topology concepts to crystallographers

The simple Euler polyhedral formula, expressed as an alternating count of the bounding faces, edges and vertices of any polyhedron, V – E + F = 2, is a fundamental concept in several branches of mathematics. Obviously, it is important in geometry, but it is also well known in topology, where a similar telescoping sum is known as the Euler characteristic χ of any finite space. The value of χ can also be computed for the unit polyhedra (such as the unit cell, the asymmetric unit or Dirichlet domain) which build, in a symmetric fashion, the infinite crystal lattices in all space groups. In this application χ has a modified form (χ m ) and value because the addends have to be weighted according to their symmetry. Although derived in geometry (in fact in crystallography), χ m has an elegant topological interpretation through the concept of orbifolds. Alternatively, χ m can be illustrated using the theorems of Harriot and Descartes, which predate the discovery made by Euler. Those historical theorems, which focus on angular defects of polyhedra, are beautifully expressed in the formula of de Gua de Malves. In a still more general interpretation, the theorem of Gauss–Bonnet links the Euler characteristic with the general curvature of any closed space. This article presents an overview of these interesting aspects of mathematics with Euler's formula as the leitmotif. Finally, a game is designed, allowing readers to absorb the concept of the Euler characteristic in an entertaining way.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Non-Hermitian quantum mechanics approach for extracting and emulating continuum physics based on bound-state-like calculations: Detailed description

Here, this work applies a reduced basis method to study the continuum physics of a finite quantum system—either few or many-body. Specifically, I develop reduced-order models, or emulators, for the underlying inhomogeneous Schrödinger equation and train the emulators against the equation's bound-state-like solutions at complex energies. The emulators rapidly and accurately interpolate and extrapolate the matrix elements of the Hamiltonian resolvent operator (Green's function) across a parameter space that includes both complex energy and other real-valued physical inputs in the Schrödinger equation. The spectra, discretized and compressed as the result of emulation, and the associated resolvent matrix elements (or amplitudes), have the defining characteristics of non-Hermitian quantum mechanics calculations, featuring complex eigenenergies with negative imaginary parts and branch cuts moved below the real axis in the complex energy plane. Therefore, one now has a method that extracts continuum physics from bound-state-like calculations and emulates those extractions in the input parameter space. Building on a prior Letter [Zhang, Phys. Rev. Lett. 135, 242501 (2025)], this article provides the full theoretical details, a comprehensive analysis of the method's performance, and a brief discussion of how it can be coupled with existing continuum approaches to perform emulations in their input parameter spaces.

ab initio calculations↗
DOE OSTI2025

Why the electron chooses one of two symmetry-related paths in the type II bacterial photosynthetic reaction centers

A series of electron and proton transfer reactions are carried out in photosynthetic reaction centers that convert the energy of light into high energy reduced products and add to the proton gradient. All photosynthetic reaction centers (photosystems I, II, and different bacterial reaction centers (bRCs)) have their charge separating cofactors arranged with two, c 2 -symmetric paths from a primary donor to a terminal acceptor. In Type II reaction centers (PSII and bRCs), electrons utilize only one active branch. In bRCs the electron transfer occurs in the A branch from the excited state of the bacteriochlorophyll (BChl) dimer P 860 through BChl A , bacteriopheophytin A (BPh A ), to ubiquinone Q A , and then to the B-branch Q B . B-branch BChl B and BPh B do not participate. A similar path is found in PSII. These proteins challenge our understanding of how proteins can turn on or turn off the activity of bound ligands. The shift of the electrochemical midpoint potential (E m ) of each cofactor by the protein environment is calculated using continuum electrostatics within the Multi-Conformation Continuum Electrostatic (MCCE) program. In bRCs from Rhodobacter sphaeroides, the electrostatic potentials favor reduction of the active branch cofactors. The residues that contribute to a more positive potential on the A branch are identified. These are often distant from the cofactors. The large difference of the potential at the competing BChls may direct the electron transfer to the A branch. Residues TyrM210 and PheL181 are at symmetrical positions on the A and B branches, near P 860 , between BChl A and BChl B . Mutated bRCs, that swap the Tyr and Phe (YF mutant), have been shown to support some electron transfer to the B-branch cofactors. Here, the calculated E m s with the YF swap show more positive B-branch potentials, which help explain why there is now electron transfer along the B branch.

Bacterial Reaction Center↗