SEARCH · Engineering Papers

Results for “perturbation analysis”

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 793 records · Page 44

All sources NASA NTRS DOE OSTI

DOE OSTI2024

An extension of Thwaites’ method for turbulent boundary layers

Thwaites developed an approximate method for determining the evolution of laminar boundary layers. The approximation follows from an assumption that the growth of a laminar boundary layer in the presence of pressure gradients could be parameterized solely as a function of the Holstein–Bohlen flow parameter, thus reducing the von Kármán momentum integral to a first-order ordinary differential equation. This method is useful for the analysis of laminar flows, and in computational potential flow solvers to account for the viscous effects. In this work, an approximate method for determining the momentum thickness of a two-dimensional, turbulent boundary layer is proposed following Thwaites’ work. It is shown that the method provides good estimates of the momentum thickness for multiple boundary layers, including both favourable and adverse pressure gradient effects, up to the point of separation. In the limit of high Reynolds numbers, it is possible to derive a criterion for the onset of separation from the proposed model, which is shown to be in agreement with prior empirical observations. The sensitivity of the separation location with respect to upstream perturbations is also analysed through this model for the NASA/Boeing speed bump and the transonic Bachalo–Johnson bump.

17 WIND ENERGY↗

DOE OSTI2023

DESI mock challenge: Halo and galaxy catalogues with the bias assignment method

We present a novel approach to the construction of mock galaxy catalogues for large-scale structure analysis based on the distribution of dark matter halos obtained with effective bias models at the field level. We aim to produce mock galaxy catalogues capable of generating accurate covariance matrices for a number of cosmological probes that are expected to be measured in current and forthcoming galaxy redshift surveys (e.g. two- and three-point statistics). The construction of the catalogues shown in this paper is part of a mock-comparison project within the Dark Energy Spectroscopic Instrument (DESI) collaboration. We use the bias assignment method ( BAM ) to model the statistics of halo distribution through a learning algorithm using a few detailed N-body simulations, and approximated gravity solvers based on Lagrangian perturbation theory. We introduce cosmic-web-dependent corrections to modelling redshift-space distortions at the N-body level – both in the halo and galaxy distributions –, as well as a multi-scale approach for accurate assignment of halo properties. Using specific models of halo occupation distributions to populate halos, we generate galaxy mocks with the expected number density and central-satellite fraction of emission-line galaxies, which are a key target of the DESI experiment. BAM generates mock catalogues with per cent accuracy in a number of summary statistics, such as the abundance, the two- and three-point statistics of halo distributions, both in real and redshift space. In particular, the mock galaxy catalogues display ~3%-10% accuracy in the multipoles of the power spectrum up to scales of k ~ 0.4 h -1 Mpc. We show that covariance matrices of two- and three-point statistics obtained with BAM display a similar structure to the reference simulation. BAM offers an efficient way to produce mock halo catalogues with accurate two- and three-point statistics and is able to generate a variety of multi-tracer catalogues with precise covariance matrices of several cosmological probes. We discuss future developments of the algorithm towards mock production in DESI and other galaxy-redshift surveys.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2024

Excited state electronic structure of dimethyl disulfide involved in photodissociation at ∼200 nm

Dimethyl disulfide (DMDS), one of the smallest organic molecules with an S–S bond, serves as a model system for understanding photofragmentation in polypeptides and proteins. Prior studies of DMDS photodissociation excited at ∼266 nm and ∼248 nm have elucidated the mechanisms of S–S and C–S bond cleavage, which involve the lowest excited electronic states S 1 and S 2 . Far less is known about the dissociation mechanisms and electronic structure of relevant excited states of DMDS excited at ∼200 nm. Herein we present calculations of the electronic structure and properties of electronic states S 1 –S 6 accessed when DMDS is excited at ∼200 nm. Our analysis includes a comparison of theoretical and experimental UV spectra, as well as theoretically predicted one-dimensional cuts through the singlet and triplet potential energy surfaces along the S–S and C–S bond dissociation coordinates. Finally, we present calculations of spin–orbit coupling constants at the Franck–Condon geometry to assess the likelihood of ultrafast intersystem crossing. We show that choosing an accurate yet computationally efficient electronic structure method for calculating the S 0 –S 6 potential energy surfaces along relevant dissociation coordinates is challenging due to excited states with doubly excited character and/or mixed Rydberg-valence character. Our findings demonstrate that the extended multi-state complete active space second-order perturbation theory (XMS-CASPT2) balances this computational efficiency and accuracy, as it captures both the Rydberg character of states in the Franck–Condon region and multiconfigurational character toward the bond-dissociation limits. In conclusion, we compare the performance of XMS-CASPT2 to a new variant of equation of motion coupled cluster theory with single, double, and perturbative triple corrections, EOM-CCSD(T)(a)*, finding that EOM-CCSD(T)(a)* significantly improves the treatment of doubly excited states compared to EOM-CCSD, but struggles to quantitatively capture asymptotic energies along bond dissociation coordinates for these states.

Rishi, Varun↗

DOE OSTI2026

Understanding polaronic transport in complex oxides by combining precise synthesis and first-principles many-body theory

In complex oxides, charge carriers often couple strongly with lattice vibrations to form polarons–entangled electron–phonon quasiparticles whose transport properties remain difficult to characterize. Experimental access to intrinsic polaronic transport requires ultraclean samples, while theoretical description demands methods beyond low-order perturbation theory. Here, we show a predictive theory–experiment workflow to study polaron transport in complex oxides. Focusing on a prototypical polaronic oxide, anatase TiO 2 , we combine growth of high-quality oxygen-vacancy-doped films using hybrid molecular beam epitaxy with a first-principles electron–phonon diagrammatic Monte-Carlo (FEP-DMC) framework recently developed for accurate polaron predictions. Our films exhibit record-high electron mobility for anatase TiO 2 , in excellent agreement with FEP-DMC calculations conducted prior to experiment, which predict a room-temperature mobility of 45 ± 15 cm −2 V −1 s −1 and a mobility-temperature scaling of μ ∝ T −1.9 ± 0.077 . Microscopic analysis using scanning transmission electron microscopy and x-ray photoelectron spectroscopy reveals the role of oxygen vacancies in modulating transport at lower temperatures. FEP-DMC further provides quantitative insight into polaron formation energy, phonon cloud distribution, lattice distortion around the polaron, and the polaronic contribution to mobility. Together, these results provide a deeper microscopic understanding of large-polaron transport in a complex oxide and provide the blueprint to characterize other polaronic materials.

anatase TiO2↗

DOE OSTI2026

Coupling of high-resolution mass spectrometer and photosynthesis system for comprehensive leaf volatile metabolite profiling

Background Leaf-level biogenic volatile organic compounds (BVOCs) emissions represent a major source of organic gases in the atmosphere, influencing both climate and air quality. These emissions are strongly driven by environmental perturbations, which affect individual plant- to ecosystem-level processes. Uncovering all the BVOCs and understanding how their emissions respond to altered environmental conditions provide critical insights into vegetation-driven changes in atmospheric chemistry. We developed a tandem instrumentation setup that integrates a proton transfer reaction time-of-flight mass spectrometer (PTR-ToF-MS) with parts-per-trillion detection limits and a photosynthetic infrared gas exchange system for the untargeted survey of all the BVOCs. This novel system enables simultaneous, real-time monitoring of BVOC emissions and photosynthetic parameters at the leaf level, offering new opportunities to disentangle the physiological and environmental drivers of VOC release. Furthermore, we established the VOC Analysis and Processing Optimization Resource (VAPOR), an open-access software tool designed for rapid data post-processing and the analysis of the variability of hundreds of BVOCs. We assessed the performance of the tandem system under varying background conditions, using standard gas mixtures and a range of environmental factors. Results Blank emissions were substantially lower for major BVOCs (e.g., isoprene) compared to those observed in plant emissions. Despite this, the observation of background-level VOCs highlights the importance of routinely acquiring and accounting for blank measurements in analyses using the coupled instrumentation. Introduction of known VOC concentrations to the system demonstrated a linear response across different compounds with varying molecular compositions, indicating minimal gas loss regardless of chemical moieties within the coupled instrumentation. We applied the optimized system to investigate the physiological mechanisms driving BVOC emissions across different genotypes of poplar and pennycress. The high mass resolution capabilities of the PTR-ToF-MS, coupled with comprehensive VAPOR-driven data analysis, enabled the identification of several important BVOCs, including methanol and methanethiol; these BVOCs displayed substantial variation across pennycress genotypes and showed concentrations ~ 100–350% higher than the blank. Moreover, isoprene emissions varied significantly among poplar genotypes grown in different potting media. Conclusions Tandem instrumentation offers a powerful tool for profiling volatile molecular markers and elucidating their genetic and environmental underpinnings. This approach enhances our ability to predict BVOC emissions in response to genotype by environmental interactions and contributes to a deeper understanding of vegetation responses to environmental changes.

Biogenic volatile organic compounds↗

DOE OSTI2020

Analysis into Asymptotic Convergence to Full Nonlinear Solutions and Exploration of the Implication of Numerical Operator Mutation of Differential Systems

A robust, sufficiently accurate and practical hydrodynamic simulation toolset is required as a key component of the modeling and simulation of air-gap electrostatic discharge events. This work was performed to complement these ongoing efforts. In particular, hydrodynamic simulations must be vetted to ensure they are robust and sufficiently accurate over relevant characteristic scales. Verification models were generated in order to cultivate the technical knowledge and expertise needed to properly create, implement and execute numerical simulations. Furthermore, this effort was utilized extensively to educate students on the mathematical and numerical principles underlying hydrodynamic simulations. This education opportunity, provided in a holistic and rigorous manner, has greatly benefited developing scientists and engineers with the necessary understandings and toolsets required to excel at accomplishing the task at hand, and, more generally, it has enabled them to generate key programmatic deliverables. This report articulates several subtilties; specifically, how perturbations, nonlinear behavior, and dissipative mechanisms influence numerical stability, how to properly structure mathematical and numerical solutions, and how to properly generate error estimation/assignment. A more rigorous discussion of the consequences of such topics can be found in the body of this report in Chapters 2 and 3 with qualitative findings discussed in Chapter 4.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2022

Tidal frequencies and quasiperiodic subsurface water level variations dominate redox dynamics in a salt marsh system

Salt marshes are hotspots of nutrient processing en route to sensitive coastal environments. Whilst our understanding of these systems has improved over the years, we still have limited knowledge of the spatiotemporal variability of critical biogeochemical drivers within salt marshes. Sea-level rise will continue to force change on salt marsh functioning, highlighting the urgency of filling this knowledge gap. Our study was conducted in a central California estuary experiencing extensive marsh drowning and relative sea-level rise, making it a model system for such an investigation. Here we instrumented three marsh positions subjected to different degrees of tidal inundation (6.7%, 8.9%, and 11.2% of the time for the upper, middle, and lower marsh positions, respectively), providing locations with varied biogeochemical characteristics and hydrological interactions at the site. We continuously monitored redox potential (Eh) at depths of 0.1, 0.3, and 0.5 m, subsurface water levels (WL), and temperature at 0.7 m depth at each marsh position. To understand how drivers of subsurface biogeochemical processes fluctuate across tidal cycles, we used wavelet analyses to explain the interactions between Eh and WL. We found that tidal forcing significantly affects key drivers of biogeochemical processes by imparting controls on Eh variability, likely driving subsurface hydro-biogeochemistry of the salt marsh. Wavelet coherence showed that the Eh-WL relationship is nonlinear, and their lead–lag relationship is variable. We found that precipitation events perturb Eh at depth over timescales of hours, even though WL shows relatively minimal change during events. Furthermore, this work highlights the importance of high frequency in situ measurements, such as Eh, to help explain factors that govern subsurface biogeochemistry and hydrological processes in salt marshes.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

Efficient 3-D velocity model building using joint inline and crossline plane-wave wave-equation migration velocity analyses

SUMMARY Wave-equation migration velocity analysis (WEMVA) is an image-domain inversion method for velocity model building. Automatic plane-wave WEMVA (PWEMVA) calculates the moveouts of plane-wave common-image gathers (CIGs) by searching a best-fitting parabola with semblance analysis and backprojects residual CIG moveouts into wavefield wave paths with a reflection tomographic kernel. However, 3-D PWEMVA is very computationally expensive because 3-D reflection tomographic inversion requires at least five 3-D reverse-time migrations per iteration and stores two types of source wavefields at model boundaries. We develop a joint inline and crossline PWEMVA method for efficient 3-D velocity model building. We alternatively implement the inline and crossline PWEMVAs with a constraint for each other, in which we iteratively construct the 3-D velocity model update through 1-D spline interpolation of 2-D gradients. The inline and crossline joint inversion is practical since PWEMVA only inverts for low-wavenumber velocity perturbations along wave paths, and the method can take less than 1 per cent of the computational cost of full 3-D PWEMVA. To construct unaliased plane waves for our joint inline and crossline PWEMVA, we develop a 3-D data interpolation method in the frequency–wavenumber (FK) domain to recover regularly and randomly missing traces. The method minimizes the misfit on sufficiently localized data subsets with iterative optimal step lengths and a gradient preconditioner that iteratively selects dominant dips along different azimuths. In numerical experiments, we use a 3-D synthetic seismic data set and a land 3-D field seismic data set acquired at the Farnsworth CO2-EOR (enhanced oil recovery) field to demonstrate the efficacy of our velocity model building and data interpolation methods.

Liu, Xuejian↗

DOE OSTI2026

Overview of the Nephele Perturbed Parameter Ensemble for Aerosol‐Cloud Interactions in E3SMv3

Aerosol-cloud interactions (aci) are the leading source of uncertainty in inferring climate sensitivity from the historical record. Earth system models (ESMs) struggle to represent aci because the processes responsible for these phenomena occur at much finer time and space scales than can be resolved by any ESM. Observational constraints provide key benchmarks to test ESMs, but cannot be used alone to fully understand aci processes except in very specific cases where causality is controlled; some degree of modeling is required to infer aci and estimate radiative forcing. Here, we generate and characterize a perturbed parameter ensemble (PPE) in version 3 of the Energy Exascale ESM (E3SMv3). We perturb 25 parameters that govern aci processes over 250 members and integrate the model over present-day and preindustrial aerosol emissions. We find that the process representation in E3SMv3 is flexible and can generate global-mean effective radiative forcings due to aci (ERFaci) ranging from −3.0 to +0.9 W m −2 . The positive ERFaci values simulated by a portion of the PPE are implausible and result from parameter combinations that produce unrealistic top-of-atmosphere energy fluxes. While global-mean cloud droplet number concentration always increases in response to anthropogenic aerosol, cloud liquid water path can both increase and decrease, suggesting that precipitation suppression is not the only aerosol-cloud adjustment represented by E3SMv3. Analysis of which processes control liquid cloud adjustment in the PPE points toward stratiform precipitation processes and aerosol activation, which is consistent with many previous ESMs, as well as the new two-moment convective cloud microphysics in E3SMv3.

Nugent, Jacqueline M. [Univ. of Wyoming, Laramie, ↗

DOE OSTI2020

Crossing symmetry, transcendentality and the Regge behaviour of 1d CFTs

We develop the technology for Polyakov-Mellin (PM) bootstrap in one- dimensional conformal field theories (CFT 1 ). By adding appropriate contact terms, we bootstrap various effective field theories in AdS 2 and analytically compute the CFT data to one loop. The computation can be extended to higher orders in perturbation theory, if we ignore mixing, for any external dimension. We develop PM bootstrap for O(N) theories and derive the necessary contact terms for such theories (which also involves a new higher gradient contact term absent for N = 1). We perform cross-checks which include considering the diagonal limit of the 2d Ising model in terms of the 1d PM blocks. As an independent check of the validity of the results obtained with PM bootstrap, we propose a suitable basis of transcendental functions, which allows to fix the four-point correlators of identical scalar primaries completely, up to a finite number of ambiguities related to the number of contact terms in the PM basis. We perform this analysis both at tree level (with and without exchanges) and at one loop. We also derive expressions for the corresponding CFT data in terms of harmonic sums. Finally, we consider the Regge limit of one-dimensional correlators and derive a precise connection between the latter and the large-twist limit of CFT data. Exploiting this result, we study the crossing equation in the three OPE limits and derive some universal constraints for the large-twist limit of CFT data in Regge-bounded theories with a finite number of exchanges.

1/N expansion↗

DOE OSTI2022

Supervised learning and the finite-temperature string method for computing committor functions and reaction rates

A central object in the computational studies of rare events is the committor function. Though costly to compute, the committor function encodes complete mechanistic information of the processes involving rare events, including reaction rates and transition-state ensembles. Under the framework of transition path theory, Rotskoff et al. [Proceedings of the 2nd Mathematical and Scientific Machine Learning Conference, Proceedings of Machine Learning Research (PLMR, 2022), Vol. 145, pp. 757–780] proposes an algorithm where a feedback loop couples a neural network that models the committor function with importance sampling, mainly umbrella sampling, which collects data needed for adaptive training. Here, in this work, we show additional modifications are needed to improve the accuracy of the algorithm. The first modification adds elements of supervised learning, which allows the neural network to improve its prediction by fitting to sample-mean estimates of committor values obtained from short molecular dynamics trajectories. The second modification replaces the committor-based umbrella sampling with the finite-temperature string (FTS) method, which enables homogeneous sampling in regions where transition pathways are located. We test our modifications on low-dimensional systems with non-convex potential energy where reference solutions can be found via analytical or finite element methods, and show how combining supervised learning and the FTS method yields accurate computation of committor functions and reaction rates. We also provide an error analysis for algorithms that use the FTS method, using which reaction rates can be accurately estimated during training with a small number of samples. The methods are then applied to a molecular system in which no reference solution is known, where accurate computations of committor functions and reaction rates can still be obtained.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Altered Cardiac Energetics and Mitochondrial Dysfunction in Hypertrophic Cardiomyopathy

Background: Hypertrophic cardiomyopathy (HCM) is a complex disease partly explained by the effects of individual gene variants on sarcomeric protein biomechanics. At the cellular level, HCM mutations most commonly enhance force production, leading to higher energy demands. Despite significant advances in elucidating sarcomeric structure–function relationships, there is still much to be learned about the mechanisms that link altered cardiac energetics to HCM phenotypes. Here, we test the hypothesis that changes in cardiac energetics represent a common pathophysiologic pathway in HCM. Methods: We performed a comprehensive multiomics profile of the molecular (transcripts, metabolites, and complex lipids), ultrastructural, and functional components of HCM energetics using myocardial samples from 27 HCM patients and 13 normal controls (donor hearts). Results: Integrated omics analysis revealed alterations in a wide array of biochemical pathways with major dysregulation in fatty acid metabolism, reduction of acylcarnitines, and accumulation of free fatty acids. HCM hearts showed evidence of global energetic decompensation manifested by a decrease in high energy phosphate metabolites (ATP, ADP, and phosphocreatine) and a reduction in mitochondrial genes involved in creatine kinase and ATP synthesis. Accompanying these metabolic derangements, electron microscopy showed an increased fraction of severely damaged mitochondria with reduced cristae density, coinciding with reduced citrate synthase activity and mitochondrial oxidative respiration. These mitochondrial abnormalities were associated with elevated reactive oxygen species and reduced antioxidant defenses. However, despite significant mitochondrial injury, HCM hearts failed to upregulate mitophagic clearance. Conclusions: Overall, our findings suggest that perturbed metabolic signaling and mitochondrial dysfunction are common pathogenic mechanisms in patients with HCM. These results highlight potential new drug targets for attenuation of the clinical disease through improving metabolic function and reducing mitochondrial injury.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2022

Doppler-Free Saturation Spectroscopy (DFSS) for Magnetic and Electric Field Measurements in an FRC plasma

The Cooperative Research and Development Agreement (CRADA) between Oak Ridge National Laboratory (ORNL) and TAE Technologies (TAE) focused on a feasibility study of Doppler-free saturation spectroscopy (DFSS) as a diagnostic strategy to measure the internal magnetic and electric fields of TAE’s C-2W field-reversed configuration (FRC) fusion device. Non-perturbative measurements of magnetic and electric fields in the core of FRCs can provide critical information needed to constrain equilibrium models and for the continued advancement of the fusion research program at TAE. DFSS is a well-established, laser-based absorption spectroscopy technique used to suppress inhomogeneous (Doppler) broadening linewidths and reveal Zeeman/Stark split hyperfine structures of target emission lines. The collaboration involved the use of ORNL’s suite of experimentally validated quantum mechanical modelling codes to simulate theoretical atomic spectra for a variety of machine parameters and laser beam configurations. TAE’s expertise was heavily used to provide measurements/models of the required machine parameters needed as inputs to the simulations. Results of an all-encompassing error analysis demonstrate the favourable capabilities of the diagnostic and provided invaluable insights as to how DFSS can be best used to accurately measure internal B/E-fields in C-2W.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2022

The information content of jet quenching and machine learning assisted observable design

Jets produced in high-energy heavy-ion collisions are modified compared to those in proton-proton collisions due to their interaction with the deconfined, strongly-coupled quark-gluon plasma (QGP). In this work, we employ machine learning techniques to identify important features that distinguish jets produced in heavy-ion collisions from jets produced in proton-proton collisions. We formulate the problem using binary classification and focus on leveraging machine learning in ways that inform theoretical calculations of jet modification: (i) we quantify the information content in terms of Infrared Collinear (IRC)-safety and in terms of hard vs. soft emissions, (ii) we identify optimally discriminating observables that are in principle calculable in perturbative QCD, and (iii) we assess the information loss due to the heavy-ion underlying event and background subtraction algorithms. We illustrate our methodology using Monte Carlo event generators, where we find that important information about jet quenching is contained not only in hard splittings but also in soft emissions and IRC-unsafe physics inside the jet. This information appears to be significantly reduced by the presence of the underlying event. We discuss the implications of this for the prospect of using jet quenching to extract properties of the QGP. Since the training labels are exactly known, this methodology can be used directly on experimental data without reliance on modeling. We outline a proposal for how such an experimental analysis can be carried out, and how it can guide future measurements.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2020

Error field correction strategies in preparation to MAST-U operation

In magnetic fusion devices, unwanted non-axisymmetric magnetic field perturbations, known as error fields (EF), can have detrimental effects on plasma stability and confinement. Such EFs may originate from several sources, such as axi-symmetric coil misalignments, coil feeds, 3D structures in the wall surrounding the plasma, presence of ferromagnetic materials near the plasma surface, blanket materials in future devices. To minimize their impact on plasma performance and on the available operational space, it is important to identify the EF sources and develop EF control strategies. MAST Upgrade (MAST-U) is a spherical tokamak which will operate in the near future after a series of enhancements from the previous MAST experiment. To deliver a machine with EF amplitude low enough to allow a high quality experimental programme, systematic analysis of the intrinsic EF sources has been carried out for poloidal eld (P) and divertor (D) coils. To deliver a machine with EF amplitude low enough to allow a high quality experimental programme, systematic analysis of the intrinsic EF sources has been carried out for poloidal field (P) and divertor (D) coils, whose magnetic eld measurements were available when writing this paper. Such measurements reveal that P and D coils are 3D deformed and thus are responsible for intrinsic EFs, with mainly n=1 and n=2 toroidal mode numbers. In preparation to MAST-U operation, both passive and active EF control strategies have been adopted for n=1 and n=2 EFs compensation. Passive EF control consisted of installing finely each D and P coil within MAST-U device so as to minimize the intrinsic n=1 EF amplitude. The optimal coil alignment has been determined based on magnetic eld measurements and the corresponding 3D electro-magnetic modelling, and envisaged coil shift and tilt of 3.2 mm and 0.7 mrad, respectively, in the case of P coil named P4. Conversely, active control will be use during MAST-U operation to reduce the n=2 EF which is associated mainly with P4 and P5, as well. Since these coils have been re-used from the MAST device, studies attempting n=2 EF control, based on MAST plasmas, have been modelled utilizing the MARS-F code to interpret experimental results in MAST and to give hints for future n=2 EF control studies in MAST-U. Here, a model-based control set for n=2 EF control has been identified which would allow for the minimization of rotation braking, of the resonant magnetic field at the q=2 and of the plasma displacement in MAST-U, simultaneously.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2026

Evaluating the limitations of Bayesian metabolic control analysis

Bayesian Metabolic Control Analysis (BMCA) is a promising framework for inferring metabolic control coefficients in data-limited scenarios, combining Bayesian inference with linear-logarithmic (lin-log) rate laws. These metabolic control coefficients quantify how changes in enzyme activities affect steady-state fluxes and metabolite concentrations across a metabolic network. However, its predictive accuracy and limitations remain underexplored. This study systematically evaluates BMCA’s ability to infer elasticity values, flux control coefficients (FCC), and concentration control coefficients (CCC) under varying data availability conditions using three synthetic metabolic network models. We demonstrate that BMCA predictions are highly dependent on the inclusion of flux and enzyme concentration data, with the omission of these datasets leading to severe inaccuracies. In our synthetic, enzyme-perturbation datasets, external metabolite concentrations had minimal impact and, in some cases, their exclusion improved predictions; when external-nutrient perturbations were introduced and those concentrations were observed, gains were at most modest. Additionally, we find that posterior estimation with both ADVI and HMC can underestimate large-magnitude elasticities in our synthetic settings, with ADVI showing somewhat higher variance under strong up-regulation; thus, recovering |elasticity| ≳ 1.5 remains challenging regardless of the inference engine. ADVI also fails to accurately infer allosteric interactions, even when regulatory effects are strong. While BMCA maintains reasonable accuracy in partially recovering the rankings of the highest FCC values, its estimates of absolute values remain constrained by prior assumptions and data limitations. Our findings reveal the BMCA algorithm’s strengths and weaknesses, providing guidance on its application in metabolic engineering, and highlighting the need for methodological refinements to enhance its predictive capabilities.

59 BASIC BIOLOGICAL SCIENCES↗

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

Strongly coupled $\mathrm{QFT}$ dynamics via $\mathrm{TQFT}$ coupling

We consider a class of quantum field theories and quantum mechanics, which we couple to $\mathbb{Z}$ N topological QFTs, in order to classify non-perturbative effects in the original theory. The $\mathbb{Z}$ N TQFT structure arises naturally from turning on a classical background field for a $\mathbb{Z}$ N 0- or 1-form global symmetry. In SU(N) Yang-Mills theory coupled to $\mathbb{Z}$ N TQFT, the non-perturbative expansion parameter is exp[-S I /N] = exp[-8π 2 /g 2 N] both in the semi-classical weak coupling domain and strong coupling domain, corresponding to a fractional topological charge configurations. To classify the non-perturbative effects in original SU(N) theory, we must use PSU(N) bundle and lift configurations (critical points at infinity) for which there is no obstruction back to SU(N). These provide a refinement of instanton sums: integer topological charge, but crucially fractional action configurations contribute, providing a TQFT protected generalization of resurgent semi-classical expansion to strong coupling. Monopole-instantons (or fractional instantons) on T 3 x $S^1_L$ can be interpreted as tunneling events in the ’t Hooft flux background in the PSU(N) bundle. The construction provides a new perspective to the strong coupling regime of QFTs and resolves a number of old standing issues, especially, fixes the conflicts between the large-N and instanton analysis. We derive the mass gap at θ = 0 and gaplessness at θ = π in $\mathbb{CP}$ 1 model, and mass gap for arbitrary θ in $\mathbb{CP}$ N-1 , N ≥ 3 on $\mathbb{R}$ 2 .

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗