Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Bounds”

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

Searching for gauge theories with the conformal bootstrap

Infrared fixed points of gauge theories provide intriguing targets for the modern conformal bootstrap program. In this work we provide some preliminary evidence that a family of gauged fermionic CFTs saturate bootstrap bounds and can potentially be solved with the conformal bootstrap. We start by considering the bootstrap for SO( N ) vector 4-point functions in general dimension D . In the large N limit, upper bounds on the scaling dimensions of the lowest SO( N ) singlet and traceless symmetric scalars interpolate between two solutions at Δ = D/ 2 – 1 and Δ = D – 1 via generalized free field theory. In 3D the critical O ( N ) vector models are known to saturate the bootstrap bounds and correspond to the kinks approaching Δ = 1 / 2 at large N . We show that the bootstrap bounds also admit another infinite family of kinks \( {\mathcal{T}}_D \) , which at large N approach solutions containing free fermion bilinears at Δ = D – 1 from below. The kinks \( {\mathcal{T}}_D \) appear in general dimensions with a D -dependent critical N * below which the kink disappears. We also study relations between the bounds obtained from the bootstrap with SO( N ) vectors, SU( N ) fundamentals, and SU( N ) × SU( N ) bi-fundamentals. We provide a proof for the coincidence between bootstrap bounds with different global symmetries. We show evidence that the proper symmetries of the underlying theories of \( {\mathcal{T}}_D \) are subgroups of SO( N ), and we speculate that the kinks \( {\mathcal{T}}_D \) relate to the fixed points of gauge theories coupled to fermions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Causality, unitarity and symmetry in effective field theory

Sum rules in effective field theories, predicated upon causality, place restrictions on scattering amplitudes mediated by effective contact interactions. Through unitarity of the S -matrix, these imply that the size of higher dimensional corrections to transition amplitudes between different states is bounded by the strength of their contributions to elastic forward scattering processes. This places fundamental limits on the extent to which hypothetical symmetries can be broken by effective interactions. All analysis is for dimension 8 operators in the forward limit. Included is a thorough derivation of all positivity bounds for a chiral fermion in SU(2) and SU(3) global symmetry representations resembling those of the Standard Model, general bounds on flavour violation, new bounds for interactions between particles of different spin, inclusion of loops of dimension 6 operators and illustration of the resulting strengthening of positivity bounds over tree-level expectations, a catalogue of supersymmetric effective interactions up to mass dimension 8 and 4 legs and the demonstration that supersymmetry unifies the positivity theorems as well as the new bounds.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Global analysis of $μ → e$ interactions in the SMEFT

We study current experimental bounds on charged lepton flavor violating (CLFV) μ – e interactions in the model-independent framework of the Standard Model Effective Field Theory (SMEFT). Assuming a generic flavor structure in the quark sector, we consider the contributions of CLFV operators to low-energy observables, including μ → eγ and μ → e conversion for quark-flavor conserving operators and CLFV meson decays for quark-flavor violating operators. At high energy, we consider limits on CLFV decays of the Higgs and Z bosons and of the top quark, and obtain bounds on operators with light quarks by recasting searches for production of eμ pairs in pp collisions at the Large Hadron Collider (LHC). We connect observables at low- and high-energy by taking into account renormalization group running and matching between CLFV operators. We also discuss the sensitivity of the future Electron-Ion Collider, where the prospective bounds are derived by imposing simple cuts on final state particles. We find that, in a single operator scenario, bounds on purely leptonic operators are dominated by μ → eγ and μ → e conversion. Semileptonic operators with down-type quarks are also dominantly constrained by low-energy observables, while LHC searches lead the bounds on up-type quark-flavor violating operators. Taking simplified multiple-coupling scenarios, we show that it is easy to evade the strongest low-energy bounds from spin-independent μ → e conversion, and that collider searches are competitive and complementary to constraints from spin-dependent μ → e conversion and other low-energy probes.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

A Celestial route to AdS bulk locality

We prove a precise form of AdS bulk locality by deriving analytical two-sided bounds on bulk Wilson coefficients. Our bounds are on the Wilson coefficients themselves, rather than their ratios, as is typically found in the literature. Inspired by the Celestial amplitudes program in flat space, we perform a Celestial transform of the CFT Mellin amplitude of four identical scalars. Using the crossing symmetric dispersion relation (CSDR), we express the resulting amplitude in terms of crossing symmetric conformal partial waves. The partial waves satisy remarkable positivity properties which along with the unitarity of the CFT prove sufficient to derive the bounds. We then employ our methods in the limit of large AdS radius and recover known bounds on flat space Wilson coefficients and new bounds on their large radius corrections. We check that the planar Mellin amplitude of four stress-tensor multiplets in $\mathcal{N}$ = 4 SYM satisfies our bounds. Finally, using null constraints, we derive a form of low-spin dominance in AdS EFTs. We find that the low-spin dominance is strongest in flat space and weakens as we move away from flat space towards higher AdS curvature.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Multi-material ALE remap with interface sharpening using high-order matrix-free finite element methods

The arbitrary Lagrangian-Eulerian (ALE) technique involves remapping field quantities from a Lagrangian mesh to an optimized mesh in a conservative, accurate and bounds-preserving manner. For methods based on arbitrary order finite elements, as described in a reference, material volume fractions are advected in pseudo-time using flux-corrected transport (FCT) without any form of interface reconstruction. In practice, this can lead to excessive propagation of small volume fractions throughout the domain. In addition, this method requires assembly of a global advection matrix to compute the bounds-preserving low-order FCT solution. In this work, we introduce a new approach for ALE remap using a high-order matrix-free technique which incorporates a flux modification to sharpen material interfaces in a conservative manner. Our approach begins with computing a bounds-preserving low-order solution to the ALE remap equations at the element level. We then compute a sharp interface solution (not guaranteed to be bounds-preserving) which comes from solving an augmented version of the ALE remap equations with a conservative flux modification which acts to sharpen material volume fractions based on their gradients and transport directions. Using the sharp interface solution, we make global corrections to the bounds-preserving solution while maintaining preservation of bounds. By blending with the sharpened solution at the global level we are able to globally conserve mass without hindering the remap pseudo-time step. This new interface-aware ALE remap method is based entirely on partial assembly techniques where globally assembled matrix operators are no longer needed, resulting in a globally matrix-free FCT method for multi-material, multi-field ALE remap with high performance on GPU architectures. We present results of our new remap method on 1D, 2D and 3D benchmarks and describe the algorithmic tailoring for GPU architectures that was developed.

Vargas, Arturo [Lawrence Livermore National Labora↗
DOE OSTI2025

Optimising the processing and storage of visibilities using lossy compression

The next-generation radio astronomy instruments are providing a massive increase in sensitivity and coverage, largely through increasing the number of stations in the array and the frequency span sampled. The two primary problems encountered when processing the resultant avalanche of data are the need for abundant storage and the constraints imposed by I/O, as I/O bandwidths drop significantly on cold storage. An example of this is the data deluge expected from the SKA Telescopes of more than 60 PB per day, all to be stored on the buffer filesystem. While compressing the data is an obvious solution, the impacts on the final data products are hard to predict. In this paper, we chose an error-controlled compressor – MGARD – and applied it to simulated SKA-Mid and real pathfinder visibility data, in noise-free and noise-dominated regimes. As the data have an implicit error level in the system temperature, using an error bound in compression provides a natural metric for compression. MGARD ensures the compression incurred errors adhere to the user-prescribed tolerance. To measure the degradation of images reconstructed using the lossy compressed data, we proposed a list of diagnostic measures, exploring the trade-off between these error bounds and the corresponding compression ratios, as well as the impact on science quality derived from the lossy compressed data products through a series of experiments. We studied the global and local impacts on the output images for continuum and spectral line examples. We found relative error bounds of as much as 10%, which provide compression ratios of about 20, have a limited impact on the continuum imaging as the increased noise is less than the image RMS, whereas a 1% error bound (compression ratio of 8) introduces an increase in noise of about an order of magnitude less than the image RMS. For extremely sensitive observations and for very precious data, we would recommend a 0.1% error bound with compression ratios of about 4. These have noise impacts two orders of magnitude less than the image RMS levels. At these levels, the limits are due to instabilities in the deconvolution methods. We compared the results to the alternative compression tool DYSCO, in both the impacts on the images and in the relative flexibility. MGARD provides better compression for similar error bounds and has a host of potentially powerful additional features.

Techniques: interferometric↗
DOE OSTI2024

Uncovering the binding nature of thiocyanate in contact ion pairs with lithium ions

Ion pair formation is a fundamental molecular process that occurs in a wide variety of systems, including electrolytes, biological systems, and materials. In solution, the thiocyanate (SCN–) anion interacts with cations to form contact ion pairs (CIPs). Due to its ambidentate nature, thiocyanate can bind through either its sulfur or nitrogen atoms, depending on the solvent. This study focuses on the binding nature of thiocyanate with lithium ions as a function of the solvents using FTIR, 2D infrared spectroscopy (2DIR) spectroscopies, and theoretical calculations. The study reveals that the SCN– binding mode (S or N end) in CIPs can be identified through 2DIR spectroscopy but not by linear IR spectroscopy. Linear IR spectroscopy shows that the CN stretch frequencies are too close to one another to separate N- and S-bound CIPs. Moreover, the IR spectrum shows that the S–C stretch presents different frequencies for the salt in different solvents, but it is related to the anion speciation rather than to its binding mode. A similar trend is observed for the anion bend. 2DIR spectra show different dynamics for N-bound and S-bound thiocyanate. In particular, the frequency–frequency correlation function (FFCF) dynamics extracted from the 2DIR spectra have a single picosecond exponential decay for N-bound thiocyanate and a biexponential decay for S-bound thiocyanate, consistent with the binding mode of the anion. Lastly, it is also observed that the binding mode also affects the line shape parameters, probably due to the different molecular mechanisms of the FFCF for N- and S-bound CIPs.

25 ENERGY STORAGE↗
DOE OSTI2020

Signaling and scrambling with strongly long-range interactions

Strongly long-range interacting quantum systems—those with interactions decaying as a power law 1/$r^α$ in the distance $\textit{r}$ on a $\textit{D}$-dimensional lattice for $\textit{α ≤ D}$—have received significant interest in recent years. They are present in leading experimental platforms for quantum computation and simulation, as well as in theoretical models of quantum-information scrambling and fast entanglement creation. Since no notion of locality is expected in such systems, a general understanding of their dynamics is lacking. In a step towards rectifying this problem, we prove two Lieb-Robinson-type bounds that constrain the time for signaling and scrambling in strongly long-range interacting systems, for which no tight bounds were previously known. In this work, our first bound applies to systems mappable to free-particle Hamiltonians with long-range hopping, and is saturable for $\textit{α ≤ D/2}$. Our second bound pertains to generic long-range interacting spin Hamiltonians and gives a tight lower bound for the signaling time to extensive subsets of the system for all $\textit{α < D}$. This many-site signaling time lower bounds the scrambling time in strongly long-range interacting systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Variational quantum algorithm for estimating the quantum Fisher information

The quantum Fisher information (QFI) quantifies the ultimate precision of estimating a parameter from a quantum state and can be regarded as a reliability measure of a quantum system as a quantum sensor. However, estimation of the QFI for a mixed state is in general a computationally demanding task. In this paper we present a variational quantum algorithm called variational quantum Fisher information estimation (VQFIE) to address this task. By estimating lower and upper bounds on the QFI, based on bounding the fidelity, VQFIE outputs a range in which the actual QFI lies. This result can then be used to variationally prepare the state that maximizes the QFI, for the application of quantum sensing. In contrast to previous approaches, VQFIE does not require knowledge of the explicit form of the sensor dynamics. We simulate the algorithm for a magnetometry setup and demonstrate the tightening of our bounds as the state purity increases. For this example, we compare our bounds with literature bounds and show that our bounds are tighter.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Approximate 𝑡-Designs in Generic Circuit Architectures

Unitary 𝑡-designs are distributions on the unitary group whose first 𝑡 moments appear maximally random. Previous work has established several upper bounds on the depths at which certain specific random quantum circuit ensembles approximate 𝑡-designs. Here we show that these bounds can be extended to any fixed architecture of Haar-random two-site gates. This is accomplished by relating the spectral gaps of such architectures to those of one-dimensional brickwork architectures. Our bound depends on the details of the architecture only via the typical number of layers needed for a block of the circuit to form a connected graph over the sites. When this quantity is bounded, the circuit forms an approximate 𝑡-design in at most linear depth. We give numerical evidence for a stronger bound that depends only on the number of connected blocks into which the architecture can be divided. We also give an implicit bound for nondeterministic architectures in terms of properties of the corresponding distribution over fixed architectures.

information scrambling↗
DOE OSTI2025

Mojo: MLIR-based Performance-Portable HPC Science Kernels on GPUs for the Python Ecosystem

We explore the performance and portability of the novel Mojo language for scientific computing workloads on GPUs. As the first language based on the LLVM’s Multi-Level Intermediate Representation (MLIR) compiler infrastructure, Mojo aims to close performance and productivity gaps by combining Python’s interoperability and CUDA-like syntax for compile-time portable GPU programming. We target four scientific workloads: a seven-point stencil (memory-bound), BabelStream (memory-bound), miniBUDE (compute-bound), and Hartree–Fock (compute-bound with atomic operations); and compare their performance against vendor baselines on NVIDIA H100 and AMD MI300A GPUs. We show that Mojo’s performance is competitive with CUDA and HIP for memory-bound kernels, whereas gaps exist on AMD GPUs for atomic operations and for fast-math compute-bound kernels on both AMD and NVIDIA GPUs. Although the learning curve and programming requirements are still fairly low-level, Mojo can close significant gaps in the fragmented Python ecosystem in the convergence of scientific computing and AI.

Godoy, William [ORNL] (ORCID:0000000225905178)↗
DOE OSTI2020

Variational Quantum Fidelity Estimation

Computing quantum state fidelity will be important to verify and characterize states prepared on a quantum computer. In this work, we propose novel lower and upper bounds for the fidelity F ( ρ , σ ) based on the ``truncated fidelity'' F ( ρ m , σ ) , which is evaluated for a state ρ m obtained by projecting ρ onto its m -largest eigenvalues. Our bounds can be refined, i.e., they tighten monotonically with m . To compute our bounds, we introduce a hybrid quantum-classical algorithm, called Variational Quantum Fidelity Estimation, that involves three steps: (1) variationally diagonalize ρ , (2) compute matrix elements of σ in the eigenbasis of ρ , and (3) combine these matrix elements to compute our bounds. Our algorithm is aimed at the case where σ is arbitrary and ρ is low rank, which we call low-rank fidelity estimation, and we prove that no classical algorithm can efficiently solve this problem under reasonable assumptions. Finally, we demonstrate that our bounds can detect quantum phase transitions and are often tighter than previously known computable bounds for realistic situations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Generation of the TSL for Zirconium Hydrides from Ab Initio Methods

Zirconium hydride (ZrHx) is a moderator material used in TRIGA and other reactors that may exist in multiple phases with varying stoichiometry, which include the δ phase and the ϵ phase. Current ENDF/B-VIII.0 ZrHx thermal scattering law (TSL) evaluations do not distinguish between phases. These sub-libraries were generated with the LEAPR module of NJOY using historic phonon spectra derived from a central force model and assume incoherent elastic scattering for both bound hydrogen and zirconium, which neglects the effects of crystal structures important for scattering from zirconium bound in ZrHx. In this work, the TSLs for hydrogen and zirconium bound in δ-ZrHx and ϵ-ZrH2 were generated from phonon spectra derived from modern ab initio lattice dynamics methods and ab initio molecular dynamics. Subsequently, TSLs for hydrogen and zirconium in ZrHx and ZrH2 were generated using the Full Law Analysis Scattering System Hub (FLASSH) code. The built-in generalized coherent elastic routine was used to generate the previously neglected elastic contribution from zirconium for this material. The present TSLs provide both a re-evaluation of the current ZrH sub-libraries and expansion of the set of TSLs available for the examination of neutrons in systems with zirconium hydride, permitting explicit treatment of δ and ϵ phases.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Optical vortex manipulation for topological quantum computation

Topological quantum computation based on Majorana bound states may enable new paths to fault-tolerant quantum computing. Several recent experiments have suggested that the vortex cores of topological superconductors, such as iron-based superconductors, may host Majorana bound states at zero energy. However, quantum computation with these zero-energy vortex bound states requires precise and fast manipulation of individual vortices, which is difficult to do in a scalable manner. To address this issue, in this study we propose a control scheme based on local heating via, for example, scanning optical microscopy to braid vortex-bound Majorana zero modes in a two-dimensional topological superconductor. First, we derive the conditions required for transporting a single vortex between two defects in the superconducting material by trapping it with a hot spot generated by local optical heating. Equipped with critical conditions for the vortex motion, we then establish the ideal material properties for vortex braiding and describe how transition errors resulting from finite speed and/or temperature can be minimized. Our work paves the way toward optical or microscopic control of zero-energy vortex bound states in two-dimensional topological superconductors.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Multilevel Convergence Analysis of Multigrid-Reduction-in-Time

This study presents a multilevel convergence framework for multigrid-reduction-in-time (MGRIT) as a generalization of previous two-grid estimates. The framework provides a priori upper bounds on the convergence of MGRIT V- and F-cycles, with different relaxation schemes, by deriving the respective residual and error propagation operators. The residual and error operators are functions of the time-stepping operator, analyzed directly and bounded in the norm, both numerically and analytically. We present various upper bounds of different computational cost and varying sharpness. These upper bounds are complemented by proposing analytic formulae for the approximate convergence factor of V-cycle algorithms that take the number of fine grid time points, the temporal coarsening factors, and the eigenvalues of the time-stepping operator as parameters. The paper concludes with supporting numerical investigations of parabolic (anisotropic diffusion) and hyperbolic (wave equation) model problems. We assess the sharpness of the bounds and the quality of the approximate convergence factors. Observations from these numerical investigations demonstrate the value of the proposed multilevel convergence framework for estimating MGRIT convergence a priori and for the design of a convergent algorithm. We further highlight that observations in the literature are captured by the theory, including that two-level Parareal and multilevel MGRIT with F-relaxation do not yield scalable algorithms and the benefit of a stronger relaxation scheme. An important observation is that with increasing numbers of levels MGRIT convergence deteriorates for the hyperbolic model problem, while constant convergence factors can be achieved for the diffusion equation. The theory also indicates that L-stable Runge--Kutta schemes are more amendable to multilevel parallel-in-time integration with MGRIT than A-stable Runge--Kutta schemes.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Lepto-axiogenesis and the scale of supersymmetry

If the Peccei-Quinn field containing the QCD axion undergoes rotations in the early universe, the dimension-five operator responsible for neutrino masses can generate a lepton asymmetry that ultimately gives rise to the observed baryon asymmetry of the Universe. This lepto-axiogenesis scenario requires a flat potential for the radial direction of the Peccei-Quinn field, naturally realized in supersymmetric models. We carefully compute the efficiency of this mechanism for the Dine-Fischler-Srednicki-Zhitnitsky (DFSZ) and Kim-Shifman-Vainshtein-Zakharov (KSVZ) axion models and place lower bounds on the masses of scalar superpartners required to reproduce the observed baryon asymmetry. For the KSVZ model, we find an efficiency for generation of the asymmetry six times larger than the previously extant computation after including scattering channels involving superpartners. In this case, the superpartner scale should be above 30 TeV for a domain wall number of one; the lower bound weakens for larger domain wall numbers. We find that the superpartner mass scale may also be as low as ~ 30 TeV for the DFSZ model. In all cases, the lower bound on the superpartner masses is inversely proportional to the sum of the squares of the neutrino masses and so can strengthen as the upper bound on the neutrino mass improves. We identify the parameter space where the axion rotation can simultaneously produce axion dark matter via kinetic misalignment; in this case it is possible to put an upper bound of order PeV on the masses of scalar superpartners.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Aligned yet large dipoles: a SMEFT study

We study a non-universal flavor scenario at the level of the Standard Model Effective Field Theory, according to which the matrix of Wilson coefficients cuW of an up-type electroweak quark dipole operator is aligned with the up-type Yukawa coupling. Such an alignment usually follows from the assumption of Minimal Flavor Violation (MFV), away from which we step by allowing the entries of cuW to be sizable along the first quark generations. A particular example, which we refer to as “inverse hierarchy MFV”, features Wilson coefficients inversely proportional to quark masses, and arises from BSM models respecting MFV and containing heavy fields that replicate the mass hierarchy of SM quarks. We then analyze the phenomenology driven by cuW at colliders and at lower-energy flavor experiments. We show that precision measurements of the process pp → Wh → γγℓν at FCC-hh could set an upper bound on |cuW| ≲ O (10 -2 )(Λ/TeV) 2 , with Λ the cutoff of the effective field theory. This bound is an order of magnitude stronger than the existing LHC bounds. Moreover, we estimate that Wh → b$\overline{b}$ℓν at HL-LHC could also give competitive bounds. In the low-energy regime, we consider bounds arising from rare kaon decays, which turn out to be loose, |c$_{uW}^{11}$ | < O (1)(Λ/TeV) 2 . We finally demonstrate that our flavor and operator assumptions can be derived from a weakly-coupled UV model, which we choose to simultaneously illustrate the UV origin of inverse hierarchy MFV.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Crossing Versus Locking: Bit Threads and Continuum Multiflows

Bit threads are curves in holographic spacetimes that manifest boundary entanglement, and are represented mathematically by continuum analogues of network flows or multiflows. Subject to a density bound, the maximum number of threads connecting a boundary region to its complement computes the Ryu–Takayanagi entropy. When considering several regions at the same time, for example in proving entropy inequalities, there are various inequivalent density bounds that can be imposed. Here, we investigate for which choices of bound a given set of boundary regions can be “locked”, in other words can have their entropies computed by a single thread configuration. We show that under the most stringent bound, which requires the threads to be locally parallel, non-crossing regions can in general be locked, but crossing regions cannot, where two regions are said to cross if they partially overlap and do not cover the entire boundary. We also show that, under a certain less stringent density bound, a crossing pair can be locked, and conjecture that any set of regions not containing a pairwise crossing triple can be locked, analogously to the situation for networks.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗