Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “vector fields”

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 883 records · Page 49
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Timelike meson form factors beyond the elastic region from lattice QCD

We present a calculation of the vector-isovector timelike form factors of the pion and the kaon using lattice quantum chromodynamics. We calculate two-point correlation functions with m π ∼ 280 MeV , extracting both the finite-volume spectrum and matrix elements for these states created from the vacuum by a vector current. After determining the coupled-channel π π , K K ¯ scattering amplitudes, we perform the necessary correction for the significant finite-volume effects present in the current matrix elements, leading to the timelike form factors. We find these to be dominated by the presence of the ρ resonance, and we extract its decay constant by an analytic continuation of the amplitudes to the resonance pole. In addition, the spacelike pion form factor is determined on the same lattice configurations, and a dispersive parametrization is used to simultaneously describe the spacelike and elastic timelike regions. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Structure formation in turbulence as an instability of effective quantum plasma

Structure formation in turbulence can be understood as an instability of “plasma” formed by fluctuations serving as effective particles. These “particles” are quantumlike in the sense that their wavelengths are non-negligible compared to the sizes of background coherent structures. The corresponding “kinetic equation” describes the Wigner matrix of the turbulent field, and the coherent structures serve as collective fields. This formalism is usually applied to manifestly quantumlike or scalar waves. Here, we show how to systematically extend it to more complex systems using compressible Navier–Stokes turbulence as an example. In this case, the fluctuation Hamiltonian is a five-dimensional matrix operator and diverse modulational modes are present. As an illustration, we calculate these modes for a sinusoidal shear flow and find two modulational instabilities. One of them is specific to supersonic flows, and the other one is a Kelvin–Helmholtz-type instability that is a generalization of the known zonostrophic instability. Our calculations are readily extendable to other types of turbulence, for example, magnetohydrodynamic turbulence in plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Efficient ptychographic phase retrieval via a matrix-free Levenberg-Marquardt algorithm

The phase retrieval problem, where one aims to recover a complex-valued image from far-field intensity measurements, is a classic problem encountered in a range of imaging applications. Modern phase retrieval approaches usually rely on gradient descent methods in a nonlinear minimization framework. Calculating closed-form gradients for use in these methods is tedious work, and formulating second order derivatives is even more laborious. Additionally, second order techniques often require the storage and inversion of large matrices of partial derivatives, with memory requirements that can be prohibitive for data-rich imaging modalities. We use a reverse-mode automatic differentiation (AD) framework to implement an efficient matrix-free version of the Levenberg-Marquardt (LM) algorithm, a longstanding method that finds popular use in nonlinear least-square minimization problems but which has seen little use in phase retrieval. Furthermore, we extend the basic LM algorithm so that it can be applied for more general constrained optimization problems (including phase retrieval problems) beyond just the least-square applications. Since we use AD, we only need to specify the physics-based forward model for a specific imaging application; the first and second-order derivative terms are calculated automatically through matrix-vector products, without explicitly forming the large Jacobian or Gauss-Newton matrices typically required for the LM method. We demonstrate that this algorithm can be used to solve both the unconstrained ptychographic object retrieval problem and the constrained “blind” ptychographic object and probe retrieval problems, under the popular Gaussian noise model as well as the Poisson noise model. We compare this algorithm to state-of-the-art first order ptychographic reconstruction methods to demonstrate empirically that this method outperforms best-in-class first-order methods: it provides excellent convergence guarantees with (in many cases) a superlinear rate of convergence, all with a computational cost comparable to, or lower than, the tested first-order algorithms.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Modern Muon Physics: Selected Issues

The muon was discovered almost hundred years ago, but still remains largely a mystery. For a long time, the only difference between a muon and an electron was the approximately 200 times greater mass of the muon. However, with the later discovery of the charm quantum number, the muon played a decisive role in the formulation of the Standard Model of elementary particles, determining together with muon neutrino and the strange s and the charm c quarks the second generation of particles of this model. The equality of the reaction constant of muon decay with the constant of the vector interaction of beta-decay played a decisive role in the formulation of the conservation of the vector current (CVC) of the weak interaction of elementary particles. Processes involving muons testified in favor of the partial conservation of the axial current of weak interaction (PCAC), that is, one pion exchange in the interaction of the lepton and hadron currents. It was the muon that gave a rise to hopes of creating an alternative energy source within the framework of Muon Catalyzed Fusion (μCF). The possible realization of this process in the future as the intensity of muon beams increases can revolutionired energy supply. Nuclear fission by muons and the study of heavy-nucleus mesoatoms introduced muon physics into several separate branches of physical research.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

A compute-bound formulation of Galerkin model reduction for linear time-invariant dynamical systems

This work aims to advance computational methods for projection-based reduced-order models (ROMs) of linear time-invariant (LTI) dynamical systems. For such systems, current practice relies on ROM formulations expressing the state as a rank-1 tensor (i.e., a vector), leading to computational kernels that are memory bandwidth bound and, therefore, ill-suited for scalable performance on modern architectures. This weakness can be particularly limiting when tackling many-query studies, where one needs to run a large number of simulations. This work introduces a reformulation, called rank-2 Galerkin, of the Galerkin ROM for LTI dynamical systems which converts the nature of the ROM problem from memory bandwidth to compute bound. We present the details of the formulation and its implementation, and demonstrate its utility through numerical experiments using, as a test case, the simulation of elastic seismic shear waves in an axisymmetric domain. We quantify and analyze performance and scaling results for varying numbers of threads and problem sizes. In conclusion, we present an end-to-end demonstration of using the rank-2 Galerkin ROM for a Monte Carlo sampling study. We show that the rank-2 Galerkin ROM is one order of magnitude more efficient than the rank-1 Galerkin ROM (the current practice) and about 970 times more efficient than the full-order model, while maintaining accuracy in both the mean and statistics of the field.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Jet-veto resummation at N 3 LL p + NNLO in boson production processes

Vetoing energetic jet activity is a crucial tool for suppressing backgrounds and enabling new physics searches at the LHC, but the introduction of a veto scale can introduce large logarithms that may need to be resummed. We present an implementation of jet-veto resummation for color-singlet processes at the level of N 3 LL p matched to fixed-order NNLO predictions. Our public code MCFM allows for predictions of a single boson, such as Z/γ*, W ± or H, or with a pair of vector bosons, such as W + W – , W ± Z or ZZ. The implementation relies on recent calculations of the soft and beam functions in the presence of a jet veto over all rapidities, with jets defined using a sequential recombination algorithm with jet radius R. However one of the ingredients that is required to reach full N 3 LL accuracy is only known approximately, hence N 3 LL p . We describe in detail our formalism and compare with previous public codes that operate at the level of NNLL. Our higher-order predictions improve significantly upon NNLL calculations by reducing theoretical uncertainties. We demonstrate this by comparing our predictions with ATLAS and CMS results.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Large Kohn Anomaly and Phonon Collapse Induced by Charge Density Wave in UPt 2⁢ Si 2

Using high-energy-resolution inelastic x-ray scattering, we observe anomalous softening and damping of the transverse acoustic phonon in UPt 2 ⁢Si 2 as the system is cooled toward the charge density wave (CDW) transition temperature 𝑇 CDW . The phonon exhibits a marked Kohn-type anomaly around the CDW wave vector 𝑄 CDW and becomes overdamped within a finite momentum range already well above 𝑇 CDW . The dispersion anomaly is consistent with potential Fermi surface nesting, which together with the extended phonon collapse indicates strong electron-phonon coupling. The transition temperature estimated from the phonon softening is markedly lower than 𝑇 CDW , consistent with a primarily electronic instability rather than a phonon-driven transition. In conclusion, our results establish UPt 2⁢ Si 2 as a prime example of a strongly correlated electron CDW system with exceptionally strong electron-phonon coupling driving phonon softening and collapse.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

A Performance and Energy Study of GPU-Resident Preconditioners for Conjugate Gradient Solvers: In the Context of Existing and Novel Approaches

Optimizing a particular subprogram out of the set of Basic (sparse) Linear Algebra Subprograms (BLAS) for a given architecture is a common topic of research. In applications, however, these BLAS functions rarely appear in isolation; usually, many of them are used together, in various combinations and with varying inputs. As the need to solve a large, sparse linear system is ubiquitous throughout HPC applications, linear solvers constitute a realistic, sufficiently complex and well-defined representative use case for composite BLAS routines. To this end, based on a representative set of matrices drawn from a diverse set of fields, we present a framework to study, from the performance and energy perspective, the efficacy of GPU- resident parallel Conjugate Gradient (CG) linear solver with different preconditioner options, including Gauss-Seidel, Jacobi, and incomplete Cholesky. We also propose a novel GPU-based preconditioner, in which the triangular solves are approximated by an iterative process. The development of this preconditioner was motivated by solving large graph Laplacian linear systems, for which the existing preconditioners either perform slow on GPU-based platforms or are not applicable. We compare the performance of these preconditioners on different hardware accelerator architectures, i.e., AMD MI250X, MI100, Nvidia A100, V100, and Jetson. Our experiments reveal performance trade-offs and provide information on how to select the best strategy for the given linear system, dictated by its properties, and the platform of interest. We demonstrate the application of our novel preconditioner for solving CG and graph Laplacian systems. Overall, the framework can be utilized as a benchmark to guide informed decisions in choosing a specific preconditioner, i.e., whether it is better to rely on the performance of a triangular solver or on the performance of sparse matrix-vector product. Finally, by considering power consumption to solve the linear systems, we report the energy footprint for the solvers.

Preconditioned Conjugate Gradient, GPUs, iterative↗
DOE OSTI2026

Nucleon axial-vector form factor and radius from radiatively corrected antineutrino scattering data

The nucleon axial-vector form factor, 𝐺 𝐴 , is critical to determine the electroweak interactions of leptons with nucleons. Important examples of processes influenced by 𝐺 𝐴 are elastic (anti)neutrino-nucleon scattering and muon capture by the proton. Sparse experimental data results in a large uncertainty on the momentum dependence of 𝐺 𝐴 and has motivated the consideration of new experimental probes and first-principles lattice quantum chromodynamics (QCD) evaluations. The comparison of new and precise theoretical predictions for 𝐺 𝐴 with future experimental data necessitates the application of radiative corrections to experimentally observable processes. We apply these corrections in the extraction of 𝐺 𝐴 and the associated axial-vector radius from the recent MINERvA antineutrino-hydrogen data, compare the effects from radiative corrections to other uncertainties in neutrino scattering experiments, and discuss the comparison of lattice QCD evaluations to experimental measurements.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Direct time-domain determination of electron-phonon coupling strengths in chromium

In this work, we report the results of an ultrafast, direct structural measurement of optically pumped phonons in a Cr thin film using ultrashort x-ray pulses from a free-electron laser. In addition to measuring and confirming the known long-wavelength dispersion relation of Cr along a particular acoustic branch, we are able to determine the relative phase of the phonons as they are generated. The Cr sample exhibits two generation mechanisms for the phonons: the releasing of a preexisting charge density wave at higher frequencies, and the creation of an acoustic strain pulse via laser heating that dominates at lower frequencies. Furthermore, for the latter mechanism, we are able to measure the frequency dependence of the time required to generate the phonons. To explain the observed magnitude and slope of the delays, we perform first-principles simulations in the framework of density functional perturbation theory and ab initio molecular dynamics to fit anharmonic phonon models. In conclusion, these results show that the wave-vector dependence of the electron-phonon coupling is the driving mechanism behind the delay times: Phase-space limitation leads to higher times near the zone center. The absolute magnitudes of the delay times measured are found to be much shorter than the equilibrium electron-phonon coupling times we compute, indicating that the coupling strength is greatly enhanced when the electronic system is out of equilibrium with the lattice, as has been seen in bismuth and other systems.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Leading-color two-loop amplitudes for four partons and a $\textit{W}$ boson in QCD

We present the leading-color two-loop QCD corrections for the scattering of four partons and a W boson, including its leptonic decay. The amplitudes are assembled from the planar two-loop helicity amplitudes for four partons and a vector boson decaying to a lepton pair, which are also used to determine the planar two-loop amplitudes for four partons and a $Z/γ^\ast$ boson with a leptonic decay. The analytic expressions are obtained by setting up a dedicated Ansatz and constraining the free parameters from numerical samples obtained within the framework of numerical unitarity. The large linear systems that must be solved to determine the analytic expressions are constructed to be in Vandermonde form. Such systems can be very efficiently solved, bypassing the bottleneck of Gaussian elimination. Our results are expressed in a basis of one-mass pentagon functions, which opens the possibility of their efficient numerical evaluation.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2019

Anatomy of deconfinement

In the weak coupling limit of ${\rm SU}(N)$ Yang-Mills theory and the ${\rm O}(N)$ vector model, explicit state counting allows us to demonstrate the existence of a partially deconfined phase: $M$ of $N$ colors deconfine, and $\frac{M}{N}$ gradually grows from zero (confinement) to one (complete deconfinement). We point out that the mechanism admits a simple interpretation in the form of spontaneous breaking of gauge symmetry. In terms of the dual gravity theory, such breaking occurs during the formation of a black hole. We speculate whether the breaking and restoration of gauge symmetry can serve as an alternative definition of the deconfinement transition in theories without center symmetry, such as QCD. We also discuss the role of the color degrees of freedom in the emergence of the bulk geometry in holographic duality.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Conservative Spin-Magnitude Change in Orbital Evolution in General Relativity

We show that physical scattering observables for compact spinning objects in general relativity can depend on additional degrees of freedom in the spin tensor beyond those described by the spin vector alone. The impulse, spin kick, and leading-order waveforms exhibit such a nontrivial dependence. A signal of this additional structure is the change in the magnitude of the spin vector under conservative Hamiltonian evolution, similar to our previous studies in electrodynamics. These additional degrees of freedom describe dynamical mass multipoles of compact objects and decouple for black holes. We also show that the conservative impulse, spin kick, and change of the additional degrees of freedom are encoded in the eikonal phase.

Classical black holes↗
DOE OSTI2021

$J/\psi $-p scattering length from the total and differential photoproduction cross sections

The J/ψ-p scattering length, αJ/ψp, can be extracted from the J/ψ photoproduction cross section near threshold using the Vector Meson Dominance (VMD) model to relate the reaction γp→J/ψp to J/ψp→J/ψp. Such estimates based on experimental data result in values for |αJ/ψp|, which are much lower than most of the theoretical predictions. In this work, we study the relations between the different results, depending on the use of the total or the differential cross sections, and the method of extrapolating the data to threshold in the case of a low-statistics data sample, such as the near threshold J/ψ photoproduction dataset. We estimate a range for |αJ/ψp| of 0.003–0.025 fm as extracted from experimental data within the VMD model and discuss possible reasons for such lower values compared to the theoretical results.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Light dark matter search with nitrogen-vacancy centers in diamonds

We propose an approach to directly search for light dark matter, such as the axion or the dark photon, by using magnetometry with nitrogen-vacancy centers in diamonds. If the dark matter couples to the electron spin, it affects the evolution of the Bloch vectors consisting of the spin triplet states, which may be detected through several magnetometry techniques. We give several concrete examples with the use of dc and ac magnetometry and estimate the sensitivity on dark matter couplings.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Light new physics in XENON1T

We examine the recently-reported low-energy electron recoil spectrum observed at the XENON1T underground dark matter direct detection experiment, in the context of new interactions with solar neutrinos. In particular we show that scalar and vector mediators with masses $\lesssim 50$ keV coupled to leptons could already leave a visible signature in the XENON1T experiment, similar to the observed peak below 7 keV. This signals that dark matter detectors are already competing with neutrino scattering experiments such as GEMMA, CHARM-II and Borexino. If these results from XENON1T are interpreted as a new signal of such physics, the parameters which fit the excess face challenges from astrophysics which seem very difficult to overcome. If they are rather viewed as a constraint on new couplings, they herald the start of an era of novel precise probes of physics beyond the standard model with dark matter detectors.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Loops and trees in generic EFTs

We consider aspects of tree and one-loop behavior in a generic 4d EFT of massless scalars, fermions, and vectors, with a particular eye to the high-energy limit of the Standard Model EFT at operator dimensions 6 and 8. First, we classify the possible Lorentz structures of operators and the subset of these that can arise at tree-level in a weakly coupled UV completion, extending the tree/loop classification through dimension 8 using functional methods. Second, we investigate how operators contribute to tree and one-loop helicity amplitudes, exploring the impact of non-renormalization theorems through dimension 8. We further observe that many dimension 6 contributions to helicity amplitudes, including rational parts, vanish exactly at one-loop level. This suggests the impact of helicity selection rules extends beyond one loop in non-supersymmetric EFTs.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Constraining light thermal inelastic dark matter with NA64

A vector portal between the Standard Model and the dark sector is a predictive and compelling framework for thermal dark matter. Through co-annihilations, models of inelastic dark matter (iDM) and inelastic Dirac dark matter (i2DM) can reproduce the observed relic density in the MeV to GeV mass range without violating cosmological limits. In these scenarios, the vector mediator behaves like a semi-visible particle, evading traditional bounds on visible or invisible resonances, and uncovering new parameter space to explain the muon (g–2) anomaly. By means of a more inclusive signal definition at the NA64 experiment, we place new constraints on iDM and i2DM using a missing energy technique. With a recast-based analysis, we contextualize the NA64 exclusion limits in parameter space and estimate the reach of the newly collected and expected future NA64 data. Our results motivate the development of an optimized search program for semi-visible particles, in which fixed-target experiments like NA64 provide a powerful probe in the sub-GeV mass range.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗