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

A Bilevel Approach for Identifying the Worst Contingencies for Nonconvex Alternating Current Power Systems

We address the bilevel optimization problem of identifying the most critical attacks to an alternating current (AC) power flow network. The upper-level binary maximization problem consists of choosing an attack that is treated as a parameter in the lower-level defender minimization problem. Instances of the lower-level global minimization problem by themselves are NP-hard due to the nonconvex AC power flow constraints, and bilevel solution approaches commonly apply a convex relaxation or approximation to allow for tractable bilevel reformulations at the cost of underestimating some power system vulnerabilities. Our main contribution is to provide an alternative branch-and-bound algorithm whose upper bounding mechanism (in a maximization context) is based on a reformulation that avoids relaxation of the AC power flow constraints in the lower-level defender problem. Lower bounding is provided with semidefinite programming (SDP) relaxed solutions to the lower-level problem. We establish finite termination with guarantees of either a globally optimal solution to the original bilevel problem, or a globally optimal solution to the SDP-relaxed bilevel problem which is included in a vetted list of upper-level attack solutions, at least one of which is a globally optimal solution to the bilevel problem. We demonstrate through computational experiments applied to IEEE case instances both the relevance of our contribution, and the effectiveness of our contributed algorithm for identifying power system vulnerabilities without resorting to convex relaxations of the lower-level problem. We conclude with a discussion of future extensions and improvements.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Parallel Solver Framework for Mixed-Integer PDE-Constrained Optimization

ROL-PEBBL is a C++, MPI-based parallel code for mixed-integer PDE-constrained optimization (MIPDECO). In these problems we wish to optimize (control, design, etc.) physical systems, which must obey the laws of physics, when some of the decision variables must take integer values. ROL-PEBBL combines a code to efficiently search over integer choices (PEBBL = Parallel Enumeration Branch-and-Bound Library) and a code for efficient nonlinear optimization, including PDE-constrained optimization (ROL = Rapid Optimization Library). In this report, we summarize the design of ROL-PEBBL and initial applications/results. For an artificial source-inversion problem, finding sources of pollution on a grid from sparse samples, ROL-PEBBLs solution for the nest grid gave the best optimization guarantee for any general solver that gives both a solution and a quality guarantee.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Extreme-scale stochastic optimization and simulation via learning-enhanced decomposition and parallelization (Final Technical Report)

Stochastic optimization and simulation models ubiquitously arise in designing and operating complex service/engineering systems. They can be extreme in scale due to high-dimensional data and decisions, and can also involve decisions made sequentially in response to newly revealed data, both causing significant computational challenge. The objective of this research is to explore a unified framework that integrates machine learning with discrete optimization and risk-averse modeling, to improve the efficiency of decomposition paradigms for stochastic optimization and simulations at extreme scale. The models we consider represent a broad class of complex decision-making problems, where 0-1 or continuous decisions are made before and/or after knowing multiple sources of uncertainties that could be correlated. We will employ machine learning methods to dynamically decide and prioritize computational procedures, including cut generation, branching, and bounding of the optimal objective. Furthermore, the research will shed new lights on the traditional decomposition algorithms for extreme-scale computing. Deliverables of the research include new modeling and computational methods for advancing the state-of-the-art research in optimization and simulation, bringing many relevant risk-averse, data-driven optimization problems in practice within the range of tractability. Examples include distributed computing server scheduling and sensor deployment for monitoring critical infrastructures. Success in this effort will enable progress in solving multiple extreme-scale problems in the complex system design and operations arising from DoE missions in energy, environment, and national security.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2024

Search for exotic decays of the Higgs boson to a pair of pseudoscalars in the $μμ bb$ and $ττ bb$ final states

A search for exotic decays of the Higgs boson (H) with a mass of 125 GeV to a pair of light pseudoscalars a 1 is performed in final states where one pseudoscalar decays to two quarks and the other to a pair of muons or leptons. A data sample of proton–proton collisions at $\sqrt{s}$ = 13 TeV corresponding to an integrated luminosity of 138 fb -1 recorded with the CMS detector is analyzed. No statistically significant excess is observed over the standard model backgrounds. Upper limits are set at 95% confidence level (CL) on the Higgs boson branching fraction to μ μ b b and to τ τ b b via a pair of a 1 s. The limits depend on the pseudoscalar mass m a 1 and are observed to be in the range (0.17–3.3) x 10 -4 and (1.7–7.7) x 10 -2 in the μ μ b b and τ τ b b final states, respectively. In the framework of models with two Higgs doublets and a complex scalar singlet (2HDM+S), the results of the two final states are combined to determine upper limits on the branching fraction ${\mathcal {B}}(\text {H} \rightarrow \text {a}_{1} \text {a}_{1} \rightarrow \ell \ell \text{ b } \text{ b})$ at 95% CL, with ℓ being a muon or a τ lepton. For different types of 2HDM+S, upper bounds on the branching fraction ${\mathcal {B}}(\text {H} \rightarrow \text {a}_{1} \text {a}_{1} )$ are extracted from the combination of the two channels. In most of the Type II 2HDM+S parameter space ${\mathcal {B}}(\text {H} \rightarrow \text {a}_{1} \text {a}_{1} )$, values above 0.23 are excluded at 95% CL for m a 1 values between 15 and 60 GeV.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Search for flavour-changing neutral-current couplings between the top quark and the photon with the ATLAS detector at s = 13 TeV

This letter documents a search for flavour-changing neutral currents (FCNCs), which are strongly suppressed in the Standard Model, in events with a photon and a top quark with the ATLAS detector. The analysis uses data collected in pp collisions at TeV during Run 2 of the LHC, corresponding to an integrated luminosity of 139 fb –1 . Both FCNC top-quark production and decay are considered. The final state consists of a charged lepton, missing transverse momentum, a b-tagged jet, one high-momentum photon and possibly additional jets. A multiclass deep neural network is used to classify events either as signal in one of the two categories, FCNC production or decay, or as background. No significant excess of events over the background prediction is observed and 95% CL upper limits are placed on the strength of left- and right-handed FCNC interactions. The 95% CL bounds on the branching fractions for the FCNC top-quark decays, estimated (expected) from both top-quark production and decay, are B(t → uγ) < 0.85 (0.88$^{+0.37}_{-0.25}$) × 10 –5 and B(t → cγ) < 4.2 (3.40$^{+1.35}_{-0.95}$) × 10 –5 for a left-handed tqγ coupling, and B(t → uγ) < 1.2 (1.20$^{+0.50}_{-1.33}$) ×10 –5 and B(t → cγ) < 4.5 (3.70$^{+1.47}_{-1.03}$) × 10 –5 for a right-handed coupling.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

RF antenna helicity dependent particle heating in a helicon source

Experiments have demonstrated that ion phenomena, such as the lower hybrid resonance, play an important role in helicon source operation. Damping of the slow branch of the bounded whistler wave at the edge of a helicon source (i.e. the Trivelpiece-Gould mode) has been correlated with the creation of energetic electrons, heating of ions at the plasma edge, and anisotropic ion heating. Here we present ion velocity distribution function measurements, electron density and temperature measurements, and magnetic fluctuation measurements on both sides of an m = |1| helical antenna in a helicon source as a function of the driving frequency, magnetic field strength, and magnetic field orientation relative to the antenna helicity. Significant electron and ion heating (up to two times larger) occurs on the side of the antenna consistent with the launch of the m = +1 mode. The electron and ion heating occurs within one electron skin depth of the plasma edge, where slow wave damping is expected. The source parameters for enhanced particle heating are also consistent with lower hybrid resonance effects, which can only occur for Trivelpiece-Gould wave excitation.

wave modes↗
DOE OSTI2025

SNoGloDe: A Structured Nonlinear Global Decomposition Solver

Large-scale optimization problems often require decomposition strategies and customized algorithms to achieve optimal solutions within a reasonable time. Building on the work of Cao and Zavala (2019) for solving nonlinear two-stage stochastic programs to global optimality, we implement and extend their approach. We generalize to optimization problems reformulated with a block-angular constraint structure (e.g., temporal decomposition). Our framework, written in Python using Pyomo, is highly customizable and enables parallel execution of the decomposition. SNoGloDe allows tailored branching strategies, lower bounding problems, and candidate generators to leverage problem-specific knowledge. To demonstrate effectiveness, we compare SNoGloDe’s performance with Gurobi on a temporally decomposed produced water case study.

algorithms↗
DOE OSTI2022

Scalable branching on dual decomposition of stochastic mixed-integer programming problems

In this work, we present a scalable branching method for the dual decomposition of stochastic mixed-integer programming. Our new branching method is based on the branching method proposed by Caroe and Schultz that creates branching disjunctions on first-stage variables only. We propose improvements to the process for creating branching disjunctions, including (1) branching on the optimal solutions of the Dantzig-Wolfe reformulation of the restricted master problem and (2) using a more comprehensive (yet simple) measure for the dispersions associated with subproblem solution infeasibility. We prove that the proposed branching process leads to an algorithm that terminates finitely, and we provide conditions under which globally optimal solutions can be identified after termination. We have implemented our new branching method, as well as the Caroe-Schultz method and a branch-and-price method, in the open-source software package DSP. Using SIPLIB test instances, we present extensive numerical results to demonstrate that the proposed branching method significantly reduces the number of node subproblems and solution times.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

$Ξ$ 𝑏 → $Ξ$ form factors from lattice QCD and standard-model predictions for $Ξ$ 𝑏 → $Ξ$⁢𝜇 + ⁢𝜇 − and $Ξ$ 𝑏 → $Ξ$⁢𝛾 decays

We present the first lattice QCD determination of the $Ξ$ 𝑏 → $Ξ$ vector, axial-vector, and tensor form factors, which are relevant for the theory of rare decays including $Ξ$ 𝑏 → $Ξ$⁢ℓ + ⁢ℓ − and $Ξ$ 𝑏 → $Ξ$⁢𝛾. The calculation is performed with 2+1 flavors of domain-wall fermions at three different lattice spacings and pion masses in the range from approximately 430 to 230 MeV. The bottom quark is implemented using an anisotropic clover action. Three-point functions with a wide range of source-sink separations and model averaging are used to extract the ground-state contributions. We fit the dependence of the form factors on the momentum transfer, the pion mass, and the lattice spacing using modified 𝑧 expansions that account for subthreshold branch cuts, and apply dispersive bounds and asymptotic behavior constraints to achieve controlled uncertainties in the full semileptonic kinematic region. Using our form factor results, we present standard model predictions for the $Ξ$$^{−}_{𝑏}$ → $Ξ$ − ⁢𝛾 and $Ξ$$^{−}_{𝑏}$ → $Ξ$ − ⁢𝜇 + ⁢𝜇 − branching fractions and two angular observables.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Production of muonic kaon atoms at high-energy colliders

We develop a framework for the formation of exotic muonic kaon atoms ($K\mu$) in semileptonic $D^{0}$ decays, using the effective weak Hamiltonian, a helicity-based treatment of the leptonic current, and a nonrelativistic bound-state projection. The resulting branching ratio, $\mathrm{BR}(D^{0}\!\to(K\mu )\nu_{\mu})=2.29\times10^{-10}$, is implemented in a ROOT-based code to estimate yields at RHIC, LHC, and STCF. We show quantitatively that $K\mu$ atoms—also produced through coalescence in the quark--gluon plasma (QGP)—provide a sensitive probe of low-momentum primordial muons and early-time electromagnetic radiation, offering complementary constraints in an otherwise unexplored phase space for thermal dilepton and photon emission. Newly estimated dissociation cross sections in detector material indicate that secondary‑vertex reconstruction should be experimentally feasible, allowing clean experimental identification of the atoms. Projected yields from QGP coalescence in LHC and RHIC heavy-ion collisions, and from $D^{0}$ decays in LHC high-luminosity $p+p$ collisions indicate that the first observation of $K\mu$ atoms is within reach. Content from this work may be used under the terms of the Creative Commons Attribution 3.0 licence. Any further distribution of this work must maintain attribution to the author(s) and the title of the work, journal citation and DOI. Article funded by SCOAP3 and published under licence by Chinese Physical Society and the Institute of High Energy Physics of the Chinese Academy of Science and the Institute of Modern Physics of the Chinese Academy of Sciences and IOP Publishing Ltd.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Analytic error analysis of cross section interpolation methods in nodal diffusion codes - II: Numerical results

This paper is the second part of a two-part paper that documents the numerical results for the partial derivatives model presented in part I. In this paper, we derive the error bounds for the analytical point-wise error expression and verify our bounds with numerical experiments. The point-wise error expressions make available, and bound, the sources that contribute to the total error of the interpolated cross section in terms of the Lagrange interpolation errors and the model form error. MPACT is used to generate two-group homogenized cross sections for Westinghouse's AP1000 Region 4 lattice to evaluate the accuracy of the bounds. Error bounds calculated over a grid are compared to numerical data for uni-variate and multi-variate interpolation. The point-wise error bounds of a typical case matrix - two branches in each state variable - are displayed for bi-variate interpolation in the state variables: moderator density, fuel temperature, and boron concentration. The error bounds are shown to be highly accurate compared to numerical results, and in accordance with the underlying physics. We then discuss and show how the sources of error contribute to the total error, and consider the improvement of each error source. Finally, we mention future work such as propagating our cross section error bounds through a reactivity calculation. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

Beta Decay of Molecular Tritium

The beta decay of tritium in the form of molecular T 2 is the basis of sensitive experiments to measure neutrino mass. The final-state electronic, vibrational, and rotational excitations modify the beta spectrum significantly and are obtained from theory. We further determine measurements of the branching ratios to specific ionization states for the isotopolog HT. Two earlier, concordant measurements gave branching ratios of HT to the bound HHe + ion of 89.5% and 93.2%, in sharp disagreement with the theoretical prediction of 55%–57%, raising concerns about the theory’s reliability in neutrino mass experiments. Our finding, 56.5(6)%, is compatible with the theoretical expectation and disagrees strongly with the previous measurements.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Kaon physics without new physics in $$ \varepsilon _K$$

Abstract Despite the observation of significant suppressions of $$b\rightarrow s\mu ^+\mu ^-$$ b → s μ + μ - branching ratios no clear sign of New Physics (NP) has been identified in $$\Delta F=2$$ Δ F = 2 observables $$\Delta M_{d,s}$$ Δ M d , s , $$\varepsilon _K$$ ε K and the mixing induced CP asymmetries $$S_{\psi K_S}$$ S ψ K S and $$S_{\psi \phi }$$ S ψ ϕ . Assuming negligible NP contributions to these observables allows to determine CKM parameters without being involved in the tensions between inclusive and exclusive determinations of $$|V_{cb}|$$ | V cb | and $$|V_{ub}|$$ | V ub | . Furthermore this method avoids the impact of NP on the determination of these parameters present likely in global fits. Simultaneously it provides SM predictions for numerous rare K and B branching ratios that are most accurate to date. Analyzing this scenario within $$Z^\prime $$ Z ′ models we point out, following the 2009 observations of Monika Blanke and ours of 2020, that despite the absence of NP contributions to $$\varepsilon _K$$ ε K , significant NP contributions to $$K^+\rightarrow \pi ^+\nu {\bar{\nu }}$$ K + → π + ν ν ¯ , $$K_{L}\rightarrow \pi ^0\nu {\bar{\nu }}$$ K L → π 0 ν ν ¯ , $$K_S\rightarrow \mu ^+\mu ^-$$ K S → μ + μ - , $$K_L\rightarrow \pi ^0\ell ^+\ell ^-$$ K L → π 0 ℓ + ℓ - , $$\varepsilon '/\varepsilon $$ ε ′ / ε and $$\Delta M_K$$ Δ M K can be present. In the simplest scenario, this is guaranteed, as far as flavour changes are concerned, by a single non-vanishing imaginary left-handed $$Z^\prime $$ Z ′ coupling $$g^L_{sd}$$ g sd L . This scenario implies very stringent correlations between the Kaon observables considered by us. In particular, the identification of NP in any of these observables implies automatically NP contributions to the remaining ones under the assumption of non-vanishing flavour conserving $$Z^\prime $$ Z ′ couplings to $$q{\bar{q}}$$ q q ¯ , $$\nu {\bar{\nu }}$$ ν ν ¯ , and $$\mu ^+\mu ^-$$ μ + μ - . A characteristic feature of this scenario is a strict correlation between $$K^+\rightarrow \pi ^+\nu {\bar{\nu }}$$ K + → π + ν ν ¯ and $$K_{L}\rightarrow \pi ^0\nu {\bar{\nu }}$$ K L → π 0 ν ν ¯ branching ratios on a branch parallel to the Grossman-Nir bound. Moreover, $$\Delta M_K$$ Δ M K is automatically suppressed as seems to be required by the results of the RBC-UKQCD lattice QCD collaboration. Furthermore, there is no NP contribution to $$K_L\rightarrow \mu ^+\mu ^-$$ K L → μ + μ - which otherwise would bound NP effects in $$K^+\rightarrow \pi ^+\nu {\bar{\nu }}$$ K + → π + ν ν ¯ . Of particular interest are the correlations of $$K^+\rightarrow \pi ^+\nu {\bar{\nu }}$$ K + → π + ν ν ¯ and $$K_{L}\rightarrow \pi ^0\nu {\bar{\nu }}$$ K L → π 0 ν ν ¯ branching ratios and of $$\Delta M_K$$ Δ M K with the ratio $$\varepsilon '/\varepsilon $$ ε ′ / ε . We investigate the impact of renormalization group effects in the context of the SMEFT on this simple scenario.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Three exceptions to the Grossman-Nir bound

A bstract We show that the Grossman-Nir (GN) bound, Br( K L → $$ {\pi}^0\nu \overline{\nu} $$ π 0 ν ν ¯ ) ≤ 4 . 3 Br( K + → $$ {\pi}^{+}\nu \overline{\nu} $$ π + ν ν ¯ ), can be violated in the presence of light new physics with flavor violating couplings. We construct three sample models in which the GN bound can be violated by orders of magnitude, while satisfying all other experimental bounds. In the three models the enhanced branching ratio Br( K L → π 0 + inv) is due to K L → π 0 Φ 1, K L → π 0 Φ 1 Φ 1 , K L → $$ {\pi}^0{\psi}_1{\overline{\psi}}_1 $$ π 0 ψ 1 ψ ¯ 1 transitions, respectively, where Φ 1 ( ψ 1 ) is a light scalar (fermion) that escapes the detector. In the three models Br( K + → π + + inv) remains very close to the SM value, while Br( K L → π 0 + inv) can saturate the present KOTO bound. Besides invisible particles in the final state (which may account for dark matter) the models require additional light mediators around the GeV-scale.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Communication-Constrained Expansion Planning for Resilient Distribution Systems

Distributed generation and remotely controlled switches have emerged as important technologies to improve the resiliency of distribution grids against extreme weather-related disturbances. Therefore it becomes important to study how best to place them on the grid in order to meet a resiliency criteria, while minimizing costs and capturing their dependencies on the associated communication systems that sustain their distributed operations. This paper introduces the Optimal Resilient Design Problem for Distribution and Communication Systems (ORDPDC) to address this need. The ORDPDC is formulated as a two-stage stochastic mixed-integer program that captures the physical laws of distribution systems, the communication connectivity of the smart grid components, and a set of scenarios that specifies which components are affected by potential disasters. The paper proposes an exact branch-and-price algorithm for the ORDPDC that features a strong lower bound and a variety of acceleration schemes to address degeneracy. The ORDPDC model and branch-and-price algorithm were evaluated on a variety of test cases with varying disaster intensities and network topologies. The results demonstrate the significant impact of the network topologies on the expansion plans and costs, as well as the computational benefits of the proposed approach.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Charm meson and charm-meson molecule in an expanding hadron gas

We present our studies on the evolution of charm mesons after the kinetic freeze-out of the expanding hadron gas produced in a heavy-ion collision and the thermal correction to a loosely bound charm-meson molecule in a pion gas. The πD* → πD* reaction rates have t-channel singularities that give contributions inversely proportional to the thermal width of the D. The ratio of the D 0 and D + production rate can differ significantly from those predicted using the measured D* branching fractions. The thermal correction to a loosely bound charm-meson molecule in a pion gas comes primarily from the complex thermal energy shift of the charm-meson constituents.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗