Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “EXPECTATION VALUE”

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 73 records · Page 4

Machine learning approach for vibronically renormalized electronic band structures

Here, we present a machine learning (ML) method for efficient computation of vibrational thermal expectation values of physical properties from first principles. Our approach is based on the nonperturbative frozen phonon formulation in which stochastic Monte Carlo algorithm is employed to sample configurations of nuclei in a supercell at finite temperatures based on a first-principles phonon model. A deep-learning neural network is trained to accurately predict physical properties associated with sampled phonon configurations, thus bypassing the time-consuming ab initio calculations. To incorporate the point-group symmetry of the electronic system into the ML model, group-theoretical methods are used to develop a symmetry-invariant descriptor for phonon configurations in the supercell. We apply our ML approach to compute the temperature dependent electronic energy gap of silicon based on density functional theory (DFT). We show that, with less than a hundred DFT calculations for training the neural network model, an order of magnitude larger number of sampling can be achieved for the computation of the vibrational thermal expectation values. Our work highlights the promising potential of ML techniques for finite temperature first-principles electronic structure methods.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Multifidelity methods for uncertainty quantification of a nonlocal model for phase changes in materials

This study is devoted to the construction of a multifidelity Monte Carlo (MFMC) method for the uncertainty quantification of a nonlocal, non-mass-conserving Cahn-Hilliard model for phase transitions with an obstacle potential. Here, we are interested in estimating the expected value of an output of interest (OoI) that depends on the solution of the nonlocal Cahn-Hilliard model. As opposed to its local counterpart, the nonlocal model captures sharp interfaces without the need for significant mesh refinement. However, the computational cost of the nonlocal Cahn-Hilliard model is higher than that of its local counterpart with similar mesh refinement, inhibiting its use for outer-loop applications such as uncertainty quantification. The MFMC method augments the desired high-fidelity, high-cost OoI with a set of lower-fidelity, lower-cost OoIs to alleviate the computational burden associated with nonlocality. Most of the computational budget is allocated to sampling the cheap surrogate models to achieve speedup, whereas the high-fidelity model is sparsely sampled to maintain accuracy. For the non-mass-conserving nonlocal Cahn-Hilliard model, the use of the MFMC method results in, for a given computational budget, about an order of magnitude reduction in the mean-squared error of the expected value of the OoI relative to that of the Monte Carlo method.

97 MATHEMATICS AND COMPUTING↗

Role of initial conditions in plasma-current coupling of gas-puff Z-pinches

Azimuthal magnetic field measurements obtained during the implosion phase of an oxygen gas-puff Z-pinch on a 500 kA peak current and 180 ns rise time linear transformer driver are presented. While a fraction of the driver current was measured within the imploding plasma, key initial conditions were found to significantly impact the delivery of current to the plasma load. The electrode geometry was modified to assist the initial dielectric breakdown and resulted in improved shot reproducibility. Optimization of the gas injection plenum pressure and timing resulted in an increase in the current coupling parameter, defined as the ratio of the measured value of $B$ $θ$ to the expected value, from 50% to 75%. Furthermore, the degree of radial expansion of the gas puff in the load region, which is suspected to lead to the observed current loss during the implosion, was reduced by shortening the valve opening duration. Additionally, a pre-embedded axial magnetic field of up to 0.2 T was found to have no significant impact on the plasma-current coupling of the oxygen implosions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Polycyclic aromatic hydrocarbons in samples of Ryugu formed in the interstellar medium

Polycyclic aromatic hydrocarbons (PAHs) contain ≲20% of the carbon in the interstellar medium. They are potentially produced in circumstellar environments (at temperatures ≳1000 kelvin), by (~10 kelvin) interstellar clouds, or by processing of carbon-rich dust grains. We report isotopic properties of PAHs extracted from samples of the asteroid Ryugu and the meteorite Murchison. The doubly- 13 C substituted compositions (Δ2× 13 C values) of the PAHs naphthalene, fluoranthene, and pyrene are 9 to 51‰ higher than values expected for a stochastic distribution of isotopes. The Δ2× 13 C values are higher than expected if the PAHs formed in a circumstellar environment, but consistent with formation in the interstellar medium. By contrast, the PAHs phenanthrene and anthracene in Ryugu samples have Δ2× 13 C values consistent with formation by higher-temperature reactions.

Science & Technology - Other Topics↗

Lifting of two-mode states in the D1-D5 CFT

We consider D1-D5-P states in the untwisted sector of the D1-D5 orbifold CFT where one copy of the seed CFT has been excited by a pair of oscillators, each being either bosonic or fermionic. While such states are BPS at the orbifold point, they will in general ‘lift’ as the theory is deformed towards general values of the couplings. We compute the expectation value of this lift at second order in the deformation parameter for the above mentioned states. We write this lift in terms of a fixed number of nested contour integrals on a given integrand; this integrand depends on the mode numbers of the oscillators in the state. We evaluate these integrals to obtain the explicit value of the lift for various subfamilies of states. At large mode numbers one observes a smooth increase of the lift with the dimension of the state h; this increase appears to follow a ~ √h behavior similar to that found analytically in earlier computations for other classes of states.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Weak scale as a trigger

Does the value of the Higgs mass parameter affect the expectation value of local operators in the Standard Model? For essentially all local operators the answer to this question is "no", and this is one of the avatars of the hierarchy problem: nothing is "triggered" when the Higgs mass parameter crosses zero. In this letter, we explore settings in which Higgs mass parameters $can$ act as a "trigger" for some local operators ${\cal O}_T$. In the Standard Model, this happens for ${\cal O}_T = {\rm Tr} (G \tilde G)$. We also introduce a "type-0" two Higgs doublet model, with a $Z_4$ symmetry, for which ${\cal O}_T = H_1 H_2$ is triggered by the Higgs masses, demanding the existence of new Higgs states necessarily comparable to or lighter than the weak scale, with no wiggle room to decouple them whatsoever. Surprisingly, this model is not yet entirely excluded by collider searches, and will be incisively probed by the high-luminosity run of the LHC, as well as future Higgs factories. We also discuss a possibility for using this trigger to explain the origin of the weak scale, invoking a landscape of extremely light, weakly interacting scalars $\phi_i$, with a coupling to ${\cal O}_T$ needed to make it possible to find vacua with small enough cosmological constant. The weak scale trigger links the tuning of the Higgs mass to that of the cosmological constant, while coherent oscillations of the $\phi_i$ can constitute dark matter.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Search for the Higgs boson decay to a Z boson and a photon in pp collisions at $\sqrt{s}$ = 13 TeV and 13.6 TeV with the ATLAS detector

A search for the Higgs boson decay to a Z boson and a photon in the ℓℓγ (ℓ = 𝑒,𝜇) final state is performed using pp collisions at $\sqrt{s}$ = 13.6 TeV recorded with the ATLAS detector at the Large Hadron Collider during 2022–2024, corresponding to an integrated luminosity of 165 fb−1. The signal yield, normalised to the Standard Model prediction, is measured to be 𝜇 = 0.9$^{+0.7}_{−0.6}$, compatible with the expected value of 𝜇 = 1.0 ± 0.7. This corresponds to an observed (expected) signal significance of 1.4 (1.5) standard deviations under the background-only hypothesis. This result is combined with that of a similar search performed with 140 fb −1 of $\sqrt{s}$ = 13 TeV pp collisions to provide the best expected sensitivity to date to this rare decay, namely an observed (expected) signal strength of 𝜇 = 1.3$^{+0.6}_{−0.5}$ (𝜇 = 1.0$^{+0.6}_{−0.5}$), corresponding to an observed (expected) significance of 2.5 (1.9) standard deviations. The measurement is consistent with the Standard Model expectation.

Aad, G. (ORCID:0000000266654934)↗

Time Evolution of Correlation Functions in Quantum Many-Body Systems

We give rigorous analytical results on the temporal behavior of two-point correlation functions—also known as dynamical response functions or Green’s functions—in closed many-body quantum systems. We show that in a large class of translation-invariant models the correlation functions factorize at late times $\langle$A ( t ) B$\rangle$ β → $\langle$ A$\rangle$ β $\langle$ B$\rangle$ β , thus proving that dissipation emerges out of the unitary dynamics of the system. We also show that for systems with a generic spectrum the fluctuations around this late-time value are bounded by the purity of the thermal ensemble, which generally decays exponentially with system size. For autocorrelation functions we provide an upper bound on the timescale at which they reach the factorized late time value. Remarkably, this bound is only a function of local expectation values and does not increase with system size. We give numerical examples that show that this bound is a good estimate in nonintegrable models, and argue that the timescale that appears can be understood in terms of an emergent fluctuation-dissipation theorem. Our study extends to further classes of two point functions such as the symmetrized ones and the Kubo function that appears in linear response theory, for which we give analogous results.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Bayesian Approach to Estimation of Water Table Elevations Using Historical Rasters as Prior Information 2019 - 20430

In cases of complex but only partially known geology and a lack of spatial control in observation well locations, water table elevation estimation is very challenging. In some cases, auxiliary information, such as observations of the movement of tracers, operation of injection and extraction wells, and calibration of groundwater models against the historical elevation data, can be combined with expert judgement to estimate flow directions in areas of sparse data and to aid in the production of more reliable contour maps (and associated rasters) than could be produced by relying on sparse well elevation data alone. Given a historical sequence of these raster maps, the question arises how to automate, to the extent possible, the process of producing new raster maps to reflect data from previous times, the current data and the operation of expert judgement. One solution is to adopt a Bayesian point of view and to regard the historical well elevation data, auxiliary information and historical raster maps as prior information. The well elevations for water table wells, as well as those for injection/extraction wells and the data associated with other relevant variables, can be viewed as predictors for the raster surface. From this prior information, we can, conditional on the values of these predictors for a new time period, compute an expected value map and a standard deviation map for the new raster. These then can be taken to specify a prior predictive distribution for the pixels in the new raster map. Then we condition the pixels, corresponding to water level observation wells within the raster, on the observed values in those wells (which in general will differ from the regression estimate) for the new time period. Given the smoothness of the water table surface, we then smooth the surface of deviations from the mean surface, based on the variograms of the historical rasters, and add this smoothed surface to the regression mean surface. The error structure of the produced raster map is defined by the regression error structure and the error due to smoothing based on the estimated variograms. This methodology has been developed and is being further refined for groundwater monitoring and remediation at LANL. It is a very flexible method that can also be applied with a variety of other predictors applied to model the water level wells in the area of interest over the historical record. The smoothness of the spatial process and its possible evolution over time can then be estimated from the residuals from this regression. This can be augmented by expert hydrogeological opinion based on site topography and hydrogeology. (authors)

54 ENVIRONMENTAL SCIENCES↗

RAPID UNSEPARATED RARE EARTH ELEMENT ANALYSES BY ISOTOPE DILUTION MULTICOLLECTOR INDUCTIVELY COUPLED MASS SPECTROMETRY (ID-MC-ICP-MS)

We present a method for analyzing rare earth elements (REEs) by isotope dilution without separation of individual elements. Utilizing a Nu Plasma 2 multicollector inductively coupled plasma mass spectrometer (MC-ICP-MS) equipped with an Aridus II desolvation nebulizer, sample introduction tuning parameters (torch position, gas flows) were investigated and used to minimize oxide formation (e.g., NdO/Nd <0.05%) and remove the need to analyze the REEs as separated “cuts.” Ratios for spike/sample isotopes in mixed synthetic standards containing varying amounts of Ce, Nd, Sm, Eu, Gd, Dy, Er, and Yb were measured accurately within 1% of IUPAC values. La and Lu in the standards were more than 2% off from expected values due to more demanding interference corrections. A USGS BCR-2 rock standard was also processed with flux fusion, iron coprecipitation and a single column pass bulk separation followed by analysis. The concentrations for Ce, Nd, Sm, Eu, Dy, Er, and Yb fell within 2% of consensus values. A United States Geological Survey (USGS) BIR-1A rock standard was also processed in triplicate by hydrofluoric and nitric acid digestion, yielding rare earth concentrations within 3% of consensus values with 0.3% relative standard deviation (n=3). This demonstrates the ability to produce high quality REE patterns of comparable quality to traditional methods involving labor and time intensive REE separation, with sample preparation and analytical times that are comparable to analyses against isotope dilution by single-collector ICP-MS.

Shen, Steve D.↗

Prototype for Unclassified Algorithm

Continuing previous efforts to investigate and develop the Unclassified Radioisotope Algorithm, the goal of the FY19-FY20 effort was to develop a prototype detector system which uses the algorithm to confirm warhead attributes related to the presence of either weapons grade plutonium (WGPu) or highly enriched uranium (HEU). The final deliverable is a prototype attribute measurement system built with common, commercially available gamma radiation detector components, capable of confirming the presence of specific, complex radioactive sources of interest, without the collection and storage of gamma energy spectra. This is accomplished by processing each gamma pulse as it is received, applying weight values based on the energy and incrementing or decrementing scalar counters which can be compared with expected values to determine if the measured source is consistent with WGPu or HEU. This report documents the design of the prototype system as well as the development of the algorithm and performance testing results. While the previously conceptualized, simple algorithm resulted in a prohibitive amount of false positives, the goal for a simple attribute measurement system capable of verifying Ba-133 and Ra-226 (weapons grade plutonium and highly enriched uranium surrogate testing sources) at over 95% accuracy with sub 5% false positive rate was demonstrated.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

Non-Boolean quantum amplitude amplification and quantum mean estimation

This paper generalizes the quantum amplitude amplification and amplitude estimation algorithms to work with non-Boolean oracles. The action of a non-Boolean oracle $U_\varphi $ on an eigenstate $\mathinner {|{x}\rangle }$ is to apply a state-dependent phase-shift $\varphi (x)$. Unlike Boolean oracles, the eigenvalues $\exp (i\varphi (x))$ of a non-Boolean oracle are not restricted to be $\pm 1$. Two new oracular algorithms based on such non-Boolean oracles are introduced. The first is the non-Boolean amplitude amplification algorithm, which preferentially amplifies the amplitudes of the eigenstates based on the value of $\varphi (x)$. Starting from a given initial superposition state $\mathinner {|{\psi _0}\rangle }$, the basis states with lower values of $\cos (\varphi )$ are amplified at the expense of the basis states with higher values of $\cos (\varphi )$. The second algorithm is the quantum mean estimation algorithm, which uses quantum phase estimation to estimate the expectation $\mathinner {\langle {\psi _0|U_\varphi |\psi _0}\rangle }$, i.e., the expected value of $\exp (i\varphi (x))$ for a random x sampled by making a measurement on $\mathinner {|{\psi _0}\rangle }$. It is shown that the quantum mean estimation algorithm offers a quadratic speedup over the corresponding classical algorithm. Both algorithms are demonstrated using simulations for a toy example. Potential applications of the algorithms are briefly discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Competing antiferromagnetic-ferromagnetic states in a d 7 Kitaev honeycomb magnet

The Kitaev model is a rare example of an analytically solvable and physically instantiable Hamiltonian yielding a topological quantum spin liquid ground state. Here we report signatures of Kitaev spin liquid physics in the honeycomb magnet Li 3 Co 2 SbO 6 , built of high-spin d 7 (Co 2+ ) ions, in contrast to the more typical low-spin d 5 electron configurations in the presence of large spin-orbit coupling. Neutron powder diffraction measurements, heat capacity, and magnetization studies support the development of a long-range antiferromagnetic order space group of C C 2/m, below T N = 11K at μ 0 H = 0T. The magnetic entropy recovered between T = 2 and 50 K is estimated to be 0.6Rln2, in good agreement with the value expected for systems close to a Kitaev quantum spin liquid state. The temperature-dependent magnetic order parameter demonstrates a β value of 0.19(3), consistent with XY anisotropy and in-plane ordering, with Ising-like interactions between layers. Further, we observe a spin-flop-driven crossover to ferromagnetic order with space group of C2/m under an applied magnetic field of μ 0 H ≈ 0.7T at T = 2K. Magnetic structure analysis demonstrates these magnetic states are competing at finite applied magnetic fields even below the spin-flop transition. Both the d 7 compass model, a quantitative comparison of the specific heat of Li 3 Co 2 SbO 6 , and related honeycomb cobaltates to the anisotropic Kitaev model further support proximity to a Kitaev spin liquid state. Furthermore, this material demonstrates the rich playground of high-spin d 7 systems for spin liquid candidates and complements known d 5 Ir- and Ru-based materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Measurement of the $B^0_\mathrm{s}\to J/\psi K^0_\mathrm{S}$ effective lifetime from proton-proton collisions at $\sqrt{s} = 13$ TeV

The effective lifetime of the $ {\textrm{B}}_{\textrm{s}}^0 $ meson in the decay $ {\textrm{B}}_{\textrm{s}}^0\to \textrm{J}/{\uppsi \textrm{K}}_{\textrm{S}}^0 $ is measured using data collected during 2016–2018 with the CMS detector in $ \sqrt{s} $ = 13 TeV proton-proton collisions at the LHC, corresponding to an integrated luminosity of 140 fb$^{−1}$. The effective lifetime is determined by performing a two-dimensional unbinned maximum likelihood fit to the $ {\textrm{B}}_{\textrm{s}}^0 $ meson invariant mass and proper decay time distributions. The resulting value of 1.59 ± 0.07(stat) ± 0.03(syst) ps is the most precise measurement to date and is in good agreement with the expected value.[graphic not available: see fulltext]

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Radiation induced non-linear oscillations in ITER baseline scenario plasmas in DIII-D

Abstract This work shows how the radiation brought about by metals or metal-equivalent radiators such as Kr and Xe produces non-linear dynamics on otherwise stationary β N flattops of DIII-D ITER Baseline Scenario demonstration discharges. The Kr and Xe gases are used to reproduce the radiative loss rates of W in present machines that operate at core temperatures much lower than the expected ITER temperature. Experiments on DIII-D with injection of Kr and Xe, as well as with sources of intrinsic metals reach the range of radiated fraction values expected in the ITER core and experience slow oscillations in temperature and radiated power. In many cases of high radiated fraction, the core temperature decreases enough for the safety factor profile to rise above the 1/1 rational surface, naturally eliminating sawteeth and occasionally producing a persistent helical core. The oscillations can be reproduced by a modified Lotka–Volterra system for temperature and radiated fraction if diffusion and noise are included, which indicates that the interplay between temperature and radiation can be the main cause of the cyclic nature of the system. A new physics based model which includes equations for temperature, density and input power can also reproduce the oscillations observed in the experiments. The present results suggest that the non-linearity of the system can be increased by the inclusion of the inherently non-linear alpha heating term, which is proportional to ∼ n e 2 T i 2 , and obtains oscillations in the model when added to an otherwise more stationary system.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

iBLAST: Incremental BLAST of new sequences via automated e-value correction

Search results from local alignment search tools use statistical scores that are sensitive to the size of the database to report the quality of the result. For example, NCBI BLAST reports the best matches using similarity scores and expect values (i.e., e-values) calculated against the database size. Given the astronomical growth in genomics data throughout a genomic research investigation, sequence databases grow as new sequences are continuously being added to these databases. As a consequence, the results (e.g., best hits) and associated statistics (e.g., e-values) for a specific set of queries may change over the course of a genomic investigation. Thus, to update the results of a previously conducted BLAST search to find the best matches on an updated database, scientists must currently rerun the BLAST search against the entire updated database, which translates into irrecoverable and, in turn, wasted execution time, money, and computational resources. To address this issue, we devise a novel and efficient method to redeem past BLAST searches by introducing iBLAST. iBLAST leverages previous BLAST search results to conduct the same query search but only on the incremental (i.e., newly added) part of the database, recomputes the associated critical statistics such as e-values, and combines these results to produce updated search results. Our experimental results and fidelity analyses show that iBLAST delivers search results that are identical to NCBI BLAST at a substantially reduced computational cost, i.e., iBLAST performs (1 + δ )/ δ times faster than NCBI BLAST, where δ represents the fraction of database growth. We then present three different use cases to demonstrate that iBLAST can enable efficient biological discovery at a much faster speed with a substantially reduced computational cost.

59 BASIC BIOLOGICAL SCIENCES↗

Correlation function distributions for O ( N ) lattice field theories in the disordered phase

Numerical computations in strongly interacting quantum field theories are often performed using Monte Carlo sampling methods. A key task in these calculations is to estimate the value of a given physical quantity from the distribution of stochastic samples that are generated using the Monte Carlo method. Typically, the sample mean and sample variance are used to define the expectation values and uncertainties of computed quantities. However, the Monte Carlo sample distribution contains more information than these basic properties, and it is useful to investigate it more generally. In this work, the exact form of the probability distributions of two-point correlation functions at zero momentum in O ( N ) lattice field theories in the disordered phase and in infinite volume are determined. These distributions allow for a robust investigation of the efficacy of the Monte Carlo sampling procedure and are shown also to allow for improved estimators of the target physical quantity to be constructed. The theoretical expectations are shown to agree with numerical calculations in the O ( 2 ) model. Published by the American Physical Society 2024

Astronomy & Astrophysics↗

Upper bound on the smuon mass from vacuum stability in the light of muon g–2 anomaly

We derive an upper bound on the smuon mass assuming that the muon anomaly is explained by the supersymmetric (SUSY) contribution. In the minimal SUSY standard model, the SUSY contribution to the muon g–2 is enhanced when the Higgsino mass parameter is large. Then, the smuon-smuon-Higgs trilinear coupling is enhanced, which may destabilize the electroweak vacuum. We calculate precisely the decay rate of the electroweak vacuum in such a case. We include one-loop effects which are crucial to determine the overall normalization of the decay rate. Requiring that the theoretical prediction of the muon anomalous magnetic moment is consistent with the observed value at the 1 and 2σ levels (equal to the central value of the observed value), we found that the lightest smuon mass should be smaller than 1.38 and 1.68 TeV (1.20 TeV) for tan β (with tan β being the ratio of the vacuum expectation values of the two Higgs bosons), respectively, and the bound is insensitive to the value of tan β.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗