Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “INTERPOLATION”

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 1,675 records · Page 93
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Constraints on the finite volume two-nucleon spectrum at 𝑚𝜋 ≈806 MeV

The low-energy, finite-volume spectrum of the two-nucleon system at a quark mass corresponding to a pion mass of 𝑚𝜋≈806 MeV is studied with lattice quantum chromodynamics (LQCD) using variational methods. The interpolating-operator sets used in [Variational study of two-nucleon systems with lattice QCD, Phys. Rev. D 107, 094508 (2023).] are extended by including a complete basis of local hexaquark operators, as well as plane-wave dibaryon operators built from products of both positive- and negative-parity nucleon operators. Results are presented for the isosinglet and isotriplet two-nucleon channels. In both channels, noticeably weaker variational bounds on the lowest few energy eigenvalues are obtained from operator sets which contain only hexaquark operators or operators constructed from the product of two negative-parity nucleons, while other operator sets produce low-energy variational bounds which are consistent within statistical uncertainties. The consequences of these studies for the LQCD understanding of the two-nucleon spectrum are investigated.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2013

Toward the excited isoscalar meson spectrum from lattice QCD

We report on the extraction of an excited spectrum of isoscalar mesons using lattice QCD. Calculations on several lattice volumes are performed with a range of light quark masses corresponding to pion masses down to about ~400 MeV. The distillation method enables us to evaluate the required disconnected contributions with high statistical precision for a large number of meson interpolating fields. We find relatively little mixing between light and strange in most J PC channels; one notable exception is the pseudoscalar sector where the approximate SU(3) F octet, singlet structure of the η, η' is reproduced. We extract exotic J PC states, identified as hybrid mesons in which an excited gluonic field is coupled to a color-octet qqbar pair, along with non-exotic hybrid mesons embedded in a qq¯-like spectrum.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2015

Phenomenological QCD equation of state for massive neutron stars

Here, we construct an equation of state for massive neutron stars based on quantum chromodynamics phenomenology. Our primary purpose is to delineate the relevant ingredients of equations of state that simultaneously have the required stiffness and satisfy constraints from thermodynamics and causality. These ingredients are (i) a repulsive density-density interaction, universal for all flavors, (ii) the color-magnetic interaction active from low to high densities, (iii) confining effects, which become increasingly important as the baryon density decreases, and (iv) nonperturbative gluons, which are not very sensitive to changes of the quark density. We use the following “3-window” description: At baryon densities below about twice normal nuclear density, 2n 0 , we use the Akmal-Pandharipande-Ravenhall (APR) equation of state, and at high densities, ≥(4–7)n 0 , we use the three-flavor Nambu-Jona-Lasinio (NJL) model supplemented by vector and diquark interactions. In the transition density region, we smoothly interpolate the hadronic and quark equations of state in the chemical potential-pressure plane. Requiring that the equation of state approach APR at low densities, we find that the quark pressure in nonconfining models can be larger than the hadronic pressure, unlike in conventional equations of state. We show that consistent equations of state of stiffness sufficient to allow massive neutron stars are reasonably tightly constrained, suggesting that gluon dynamics remains nonperturbative even at baryon densities ~10n 0 .

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2015

Resonances in coupled πK, ηK scattering from lattice QCD

Coupled-channel πK and ηK scattering amplitudes are determined by studying the finite-volume energy spectra obtained from dynamical lattice QCD calculations. Using a large basis of interpolating operators, including both those resembling a qq-bar construction and those resembling a pair of mesons with relative momentum, a reliable excited-state spectrum can be obtained. Working at mπ = 391 MeV, we find a gradual increase in the JP = 0+ πK phase-shift which may be identified with a broad scalar resonance that couples strongly to πK and weakly to ηK. The low-energy behavior of this amplitude suggests a virtual bound-state that may be related to the κ resonance. A bound state with JP = 1- is found very close to the πK threshold energy, whose coupling to the πK channel is compatible with that of the experimental K*(892). Evidence is found for a narrow resonance in JP = 2+. Isospin–3/2 πK scattering is also studied and non-resonant phase-shifts spanning the whole elastic scattering region are obtained.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2015

Cosmology Constraints from the Weak Lensing Peak Counts and the Power Spectrum in CFHTLenS

Lensing peaks have been proposed as a useful statistic, containing cosmological information from non-Gaussianities that is inaccessible from traditional two-point statistics such as the power spectrum or two-point correlation functions. Here we examine constraints on cosmological parameters from weak lensing peak counts, using the publicly available data from the 154 deg2 CFHTLenS survey. We utilize a new suite of ray-tracing N-body simulations on a grid of 91 cosmological models, covering broad ranges of the three parameters Ω m , σ 8 , and w, and replicating the galaxy sky positions, redshifts, and shape noise in the CFHTLenS observations. We then build an emulator that interpolates the power spectrum and the peak counts to an accuracy of ≤ 5%, and compute the likelihood in the three-dimensional parameter space (Ω m , σ 8 , w) from both observables. We find that constraints from peak counts are comparable to those from the power spectrum, and somewhat tighter when different smoothing scales are combined. Neither observable can constrain w without external data. When the power spectrum and peak counts are combined, the area of the error “banana” in the (Ω m , σ 8 ) plane reduces by a factor of ≈ two, compared to using the power spectrum alone. For a flat Λ cold dark matter model, combining both statistics, we obtain the constraint σ 8 (Ω m /0.27)0.63 = 0.85 +0.03 -0.03 .

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

$B \to \pi \ell \nu$ and $B_s \to K \ell \nu$ form factors and $|V_{ub}|$ from 2+1-flavor lattice QCD with domain-wall light quarks and relativistic heavy quarks

We calculate the form factors for B → πℓν and B s → Kℓν decay in dynamical lattice quantum chromodynamics (QCD) using domain-wall light quarks and relativistic b-quarks. We use the (2+1)-flavor gauge-field ensembles generated by the RBC and UKQCD collaborations with the domain-wall fermion action and Iwasaki gauge action. For the b-quarks we use the anisotropic clover action with a relativistic heavy-quark interpretation. We analyze data at two lattice spacings of a ≈ 0.11, 0.086 fm with unitary pion masses as light as M π ≈ 290 MeV. We simultaneously extrapolate our numerical results to the physical light-quark masses and to the continuum and interpolate in the pion/kaon energy using SU(2) “hard-pion” chiral perturbation theory for heavy-light meson form factors. We provide complete systematic error budgets for the vector and scalar form factors f + (q 2 ) and f 0 (q 2 ) for both B → πℓν and B s → Kℓν at three momenta that span the q 2 range accessible in our numerical simulations. Next we extrapolate these results to q 2 = 0 using a model-independent z-parametrization based on analyticity and unitarity. We present our final results for f + (q 2 ) and f 0 (q 2 )as the coefficients of the series in z and the matrix of correlations between them; this provides a parametrization of the form factors valid over the entire allowed kinematic range. Our results agree with other three-flavor lattice-QCD determinations using staggered light quarks, and have comparable precision, thereby providing important independent cross-checks. Both B → πℓν and B s → Kℓν decays enable determinations of the Cabibbo-Kobayashi-Maskawa matrix element |V ub |. Furthermore, we perform a combined z-fit of our numerical B → πℓν form-factor data with the experimental measurements of the branching fraction from BABAR and Belle leaving the relative normalization as a free parameter; we obtain |V ub | = 3.61(32)×10 -3 , where the error includes statistical and all systematic uncertainties. The same approach can be applied to the decay B s → Kℓν to provide an alternative determination of |V ub | once the process has been measured experimentally. In anticipation of future experimental measurements, we make predictions for B → πℓν and B s → Kℓν differential branching fractions and forward-backward asymmetries in the Standard Model.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2015

Charm and strange quark masses and f D s from overlap fermions

We use overlap fermions as valence quarks to calculate meson masses in a wide quark mass range on the 2 + 1 -flavor domain-wall fermion gauge configurations generated by the RBC and UKQCD Collaborations. The well-defined quark masses in the overlap fermion formalism and the clear valence quark mass dependence of meson masses observed from the calculation facilitate a direct derivation of physical current quark masses through a global fit to the lattice data, which incorporates O($a^2$) and O($m^4a^4$) corrections, chiral extrapolation, and quark mass interpolation. Using the physical masses of D s , D$^*_s$ and J/ψ as inputs, Sommer's scale parameter r 0 and the masses of charm quark and strange quark in the $\overline{MS}$ scheme are determined to be r 0 = 0.465 (4) (9) fm, $m^{\overline{MS}}_c$ (2 GeV ) = 1.118 (6) (24) GeV (or $m^{\overline{MS}}_c(m_c)$ = 1.304 (5) (20) GeV), and $m^{\overline{MS}}_s$ (2 GeV ) = 0.101 (3) (6) GeV, respectively. Furthermore, we observe that the mass difference of the vector meson and the pseudoscalar meson with the same valence quark content is proportional to the reciprocal of the square root of the valence quark masses. The hyperfine splitting of charmonium, M J/ψ -$M_{η_c}$, is determined o be 119(2)(7) MeV, which is in good agreement with the experimental value. Finally, we also predict the decay constant of D s to be $f_{D_S}$ = 254 (2) (4) MeV. The masses of charmonium P-wave states χ c0 , χ c1 and $h_c$ are also in good agreement with experiments.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2016

B*Bπ coupling using relativistic heavy quarks

We report on a calculation of the B*Bπ coupling in lattice QCD. The strong matrix element (Bπ|B*) is directly related to the leading order low-energy constant in heavy meson chiral perturbation theory (HM Χ PT) for B mesons. We carry out our calculation directly at the b-quark mass using a non-perturbatively tuned clover action that controls discretization effects of order |p → a| and (ma) n for all n. Our analysis is performed on RBC/UKQCD gauge configurations using domain-wall fermions and the Iwasaki gauge action at two lattice spacings of a –1 = 1.729(25) GeV, a –1 = 2.281 (28) GeV, and unitary pion masses down to 290 MeV. We achieve good statistical precision and control all systematic uncertainties, giving a final result for the HM Χ PT coupling g b = 0.56(3) stat (7) sys in the continuum and at the physical light-quark masses. Furthermore, this is the first calculation performed directly at the physical b-quark mass and lies in the region one would expect from carrying out an interpolation between previous results at the charm mass and at the static point.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2016

Gradient Flow and Scale Setting on MILC HISQ Ensembles

We report on a scale determination with gradient-flow techniques on the N f = 2 + 1 + 1 HISQ ensembles generated by the MILC collaboration. The ensembles include four lattice spacings, ranging from approximately 0.15 to 0.06 fm, and both physical and unphysical values of the quark masses. The scales p √t 0 /a and w 0 /a and their tree-level improvements,√t 0;imp and w 0;imp , are computed on each ensemble using Symanzik ow and the cloverleaf definition of the energy density E. Using a combination of continuum chiral perturbation theory and a Taylor-series ansatz for the lattice-spacing and strong-coupling dependence, the results are simultaneously extrapolated to the continuum and interpolated to physical quark masses. We also determine the scales p t 0 = 0:1416( +8 -5 ) fm and w 0 = 0:1717( +12 -11 ) fm, where the errors are sums, in quadrature, of statistical and all systematic errors. The precision of w 0 and √t 0 is comparable to or more precise than the best previous estimates, respectively. We also find the continuum mass-dependence of w 0 that will be useful for estimating the scales of other ensembles. Furthermore, we estimate the integrated autocorrelation length of . For long flow times, the autocorrelation length of appears to be comparable to or smaller than that of the topological charge.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2016

Energy Estimation of Cosmic Rays with the Engineering Radio Array of the Pierre Auger Observatory

The Auger Engineering Radio Array (AERA) is part of the Pierre Auger Observatory and is used to detect the radio emission of cosmic-ray air showers. These observations are compared to the data of the surface detector stations of the Observatory, which provide well-calibrated information on the cosmic-ray energies and arrival directions. The response of the radio stations in the 30 to 80MHz regime has been thoroughly calibrated to enable the reconstruction of the incoming electric field. For the latter, the energy density is determined from the radio pulses at each observer position and is interpolated using a two dimensional function that takes into account signal asymmetries due to interference between the geomagnetic and charge excess emission components. We found that the spatial integral over the signal distribution gives a direct measurement of the energy transferred from the primary cosmic ray into radio emission in the AERA frequency range. We measure 15.8MeV of radiation energy for a 1 EeV air shower arriving perpendicularly to the geomagnetic field. This radiation energy – corrected for geometrical effects – is used as a cosmic-ray energy estimator. Performing an absolute energy calibration against the surface-detector information, we observe that this radio-energy estimator scales quadratically with the cosmic-ray energy as expected for coherent emission. Finally we find an energy resolution of the radio reconstruction of 22% for the data set and 17% for a high-quality subset containing only events with at least five radio stations with signal.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

P -wave π π scattering and the ρ resonance from lattice QCD

Here we calculate the parameters describing elastic I = 1, P-wave ππ scattering using lattice QCD with 2 +1 flavors of clover fermions. Our calculation is performed with a pion mass of m π ≈ 320 MeV and a lattice size of L ≈ 3.6 fm. We construct the two-point correlation matrices with both quark-antiquark and two-hadron interpolating fields using a combination of smeared forward, sequential and stochastic propagators. The spectra in all relevant irreducible representations for total momenta |$\vec{P}$| ≤ $\sqrt{3}\frac{2π}{L}$ are extracted with two alternative methods: a variational analysis as well as multiexponential matrix fits. We perform an analysis using Lüscher’s formalism for the energies below the inelastic thresholds, and investigate several phase shift models, including possible nonresonant contributions. We find that our data are well described by the minimal Breit-Wigner form, with no statistically significant nonresonant component. In determining the ρ resonance mass and coupling we compare two different approaches: fitting the individually extracted phase shifts versus fitting the t-matrix model directly to the energy spectrum. We find that both methods give consistent results, and at a pion mass of am π = 0.18295 (36) stat obtain gρππ = 5.69(13) stat (16) sys , am ρ = 0.4609(16) stat (14) sys , and am ρ /am N = 0.7476(38) stat (23) sys , where the first uncertainty is statistical and the second is the systematic uncertainty due to the choice of fit ranges.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Fluctuations in quantum mechanics and field theories from a new version of semiclassical theory. II.

This is the second paper on the semiclassical approach based on the density matrix given by the Euclidean time path integral with fixed coinciding end points. The classical path, interpolating between this point and the classical vacuum (called a “flucton”), as well as systematic one- and two-loop corrections were calculated in the first paper [M. A. Escobar-Ruiz, E. Shuryak, and A. V. Turbiner, Phys. Rev. D 93, 105039 (2016)] for a double-well potential. Here, we extend them for a number of quantum-mechanical problems, such as an anharmonic oscillator and the sine-Gordon potential. The method is based on a systematic expansion in Feynman diagrams and thus can be extended to quantum field theories (QFTs). We show that the loop expansion in quantum mechanics resembles the leading-log approximations in QFT. In this sequel, we present a complete set of results obtained using this method in a unified way. Alternatively, starting from the Schrödinger equation we derive a generalized Bloch equation whose semiclassical-like, iterative solution generates the loop expansion. We rederive the two-loop expansions for all three of the above potentials and extend them to three loops, which has not yet been done via Feynman diagrams. All results for both methods are fully consistent with each other. An asymmetric (tilted) double-well potential (nondegenerate minima) is also studied using the second method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Minkowski conformal blocks and the Regge limit for Sachdev-Ye-Kitaev-like models

We discuss scattering in a conformal field theory (CFT) via the conformal partial-wave analysis and the Regge limit. The focus of this paper is on understanding an operator product expansion (OPE) with Minkowski conformal blocks. Starting with a t-channel OPE, it leads to an expansion for an s-channel scattering amplitude in terms of t-channel exchanges. By contrasting with Euclidean conformal blocks we see a precise relationship between conformal blocks in the two limits without preforming an explicit analytic continuation. We discuss a generic feature for a CFT correlation function having singular F (M) (u,v) ~ u –δ , δ > 0, in the limit u → 0 and v → 1. Here, δ = (ℓ eff – 1)/2, with ℓ eff serving as an effective spin and it can be determined through an OPE. In particular, it is bounded from above, ℓ eff ≤ 2, for all CFTs with a gravity dual, and it can be associated with string modes interpolating the graviton in AdS. This singularity is historically referred to as the Pomeron. This bound is nearly saturated by SYK-like effective d = 1 CFT, and its stringy and thermal corrections have piqued current interests. Our analysis has been facilitated by dealing with Wightman functions. We provide a direct treatment in diagonalizing dynamical equations via harmonic analysis over physical scattering regions. As an example these methods are applied to the SYK model.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2019

Quantum chameleons

We initiate a quantum treatment of chameleonlike particles, deriving classical and quantum forces directly from the path integral. It is found that the quantum force can potentially dominate the classical one by many orders of magnitude.We calculate the quantum chameleon pressure between infinite plates, which is found to interpolate between the Casimir and the integrated Casimir-Polder pressures, respectively in the limits of full screening and no screening. To this end we calculate the chameleon propagator in the presence of an arbitrary number of one-dimensional layers of material. For the Eöt-Wash experiment, the five-layer propagator is used to take into account the intermediate shielding sheet, and it is found that the presence of the sheet enhances the quantum pressure by two orders of magnitude. As an example of implication, we show that in both the standard chameleon and symmetron models, large and previously unconstrained regions of the parameter space are excluded once the quantum pressure is taken into account.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

First-principles equation of state database for warm dense matter computation

We put together a first-principles equation of state (FPEOS) database for matter at extreme conditions by combining results from path integral Monte Carlo and density functional molecular dynamics simulations of the elements H, He, B, C, N, O, Ne, Na, Mg, Al, and Si as well as the compounds LiF , B 4 C , BN , CH 4 , CH 2 , C 2 H 3 , CH , C 2 H , MgO , and MgSiO 3 . For all these materials, we provide the pressure and internal energy over a density-temperature range from ~ 0.5 to 50 g cm - 3 and from ~ 10 4 to 10 9 K, which are based on ~ 5000 different first-principles simulations. We compute isobars, adiabats, and shock Hugoniot curves in the regime of L - and K -shell ionization. Invoking the linear mixing approximation, we study the properties of mixtures at high density and temperature. Furthermore, we derive the Hugoniot curves for water and alumina as well as for carbon-oxygen, helium-neon, and CH-silicon mixtures. We predict the maximal shock compression ratios of H 2 O , H 2 O 2 , Al 2 O 3 , CO , and CO 2 to be 4.61, 4.64, 4.64, 4.89, and 4.83, respectively. Finally we use the FPEOS database to determine the points of maximum shock compression for all available binary mixtures. We identify mixtures that reach higher shock compression ratios than their end members. We discuss trends common to all mixtures in pressure-temperature and particle-shock velocity spaces. In the Supplemental Material, we provide all FPEOS tables as well as computer codes for interpolation, Hugoniot calculations, and plots of various thermodynamic functions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Modeling electrokinetic flows with the discrete ion stochastic continuum overdamped solvent algorithm

In this article we develop an algorithm for the efficient simulation of electrolytes in the presence of physical boundaries. In previous work the discrete ion stochastic continuum overdamped solvent (DISCOS) algorithm was derived for triply periodic domains, and was validated through ion-ion pair correlation functions and Debye-Hückel-Onsager theory for conductivity, including the Wien effect for strong electric fields. In extending this approach to include an accurate treatment of physical boundaries we must address several important issues. First, the modifications to the spreading and interpolation operators necessary to incorporate interactions of the ions with the boundary are described. Next we discuss the modifications to the electrostatic solver to handle the influence of charges near either a fixed potential or dielectric boundary. An additional short-ranged potential is also introduced to represent interaction of the ions with a solid wall. Finally, the dry diffusion term is modified to account for the reduced mobility of ions near a boundary, which introduces an additional stochastic drift correction. Herein, several validation tests are presented confirming the correct equilibrium distribution of ions in a channel. Additionally, the methodology is demonstrated using electro-osmosis and induced-charge electro-osmosis, with comparison made to theory and other numerical methods. Notably, the DISCOS approach achieves greater accuracy than a continuum electrostatic simulation method. We also examine the effect of under-resolving hydrodynamic effects using a “dry diffusion” approach, and find that considerable computational speedup can be achieved with a negligible impact on accuracy.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Atomic dynamics in fluids: Normal mode analysis revisited

Developing microscopic understanding of the thermal properties of liquids is challenging due to their strong dynamic disorder, which prevents characterization of the atomic degrees of freedom. There have been significant research interests in the past few decades to extend the normal mode analysis for solids to instantaneous structures of liquids. However, the nature of normal modes that arise from these unstable structures is still elusive. Here, in this paper, we explore the instantaneous eigenmodes of dynamical matrices of various Lennard-Jones argon liquid and gas systems at high temperatures and show that the normal modes can be interpreted as an interpolation of T →∞ (gas) and T=0 (solid) mode descriptions. We find that normal modes become increasingly collisional and translational, recovering atomistic gaslike behavior rather than vibrational with increase in temperature, suggesting that normal modes in liquids may be described by both solidlike and gaslike modes.

74 ATOMIC AND MOLECULAR PHYSICS↗