Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “statistical physics”

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,063 records · Page 59

Analysis of differential scanning calorimetry data for aged plutonium

Differential scanning calorimetry data for samples of a 52 year old plutonium alloy with 3.3 at. % Ga that were heated beyond the melting point is analyzed using transition state theory to find activation energies for the δ to ε and ε to liquid phase transitions. A Bayesian statistical method involving a Gaussian process model is used to find mean values and confidence intervals for the activation energies. The activation energy for the δ to ε phase transition increases by 3.3 ± 3.8% per decade, relative to the case when all age related plutonium lattice point defects have been removed through annealing. The corresponding increase in activation energy for the ε to liquid transition is shown to be 7.1 ± 1.8% per decade. It is postulated that the change in activation energy with age for both phase transitions is caused, in part, by the accumulation of the same type of lattice point defects associated with the observed increase in elastic bulk modulus over time.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DESI DR2 Galaxy luminosity functions

We present galaxy luminosity functions (LFs) for the Dark Energy Spectroscopic Instrument (DESI) DR2 Bright Galaxy Survey (BGS) in the g, r, z, and $w1$ bands over $0.002\lt z\lt 0.6$. Our analysis uses updated k-corrections and evolutionary corrections, including new polynomial kcorrection fits derived from BGS Year 1 data that supersede earlier GAMA-based prescriptions. Exploiting the statistical power of DESI, we measure LFs to very faint magnitudes, reaching $^{0.1}M_r-5\log h\sim -10$. Independent measurements from the North and South survey regions show excellent agreement around the LF knee, but the very small statistical uncertainties reveal that simple analytic forms fail to capture the full LF shape. The bright end departs from a pure exponential decline, while the faint end exhibits complex, non-powerlaw behaviour, including a pronounced upturn at $^{0.1}M_r-5\log h\gtrsim -15$, which is stronger for red galaxies than for blue. We show that our LFs are largely complete for galaxies with surface brightness $\mu _{50}\lt 25$, and that an apparent steepening fainter than $-13$ is driven primarily by local overdensity and fragmentation of large galaxies. A systematic North–South offset at the brightest magnitudes is traced to red galaxies and may reflect shallower North photometry underestimating extended earlytype profiles, although this remains inconclusive. We therefore also provide LFs based on model Petrosian magnitudes. Redshift splitting reveals small but significant residuals, indicating limitations of a simple global evolutionary model. Using the redshift limits of J. Loveday et al. (2012), we find excellent agreement with GAMA, with substantially reduced statistical errors. These measurements provide a precise reference for studies of environmental and population-dependent LFs and for testing galaxy formation models.

79 ASTRONOMY AND ASTROPHYSICS↗

Jacobian-free Newton–Krylov method for the simulation of non-thermal plasma discharges with high-order time integration and physics-based preconditioning

A preconditioning framework for the numerical simulation of non-thermal streamer discharges is developed using the Jacobian-free Newton-Krylov (JFNK) method. A reduced plasma fluid model is considered, consisting of electrons, one positive ion, one negative ion, and the electrostatic potential. Here, the plasma kinetics model includes ionization, electron-ion recombination, electron attachment, electron detachment, and ion-ion recombination. The governing equations are made dimensionless, discretized in space with finite differences, and integrated in time with a fully implicit method based on high-order backward differentiation formulas. The preconditioning framework is based on a linearized form of the governing equations and physics-based operator splitting. The efficiency of the preconditioning strategy is assessed through two test cases: streamer propagation between parallel plates and an axisymmetric pin-to-pin discharge. The fully implicit approach overcomes traditional restrictions in the time step size due to processes such as electron drift, electron diffusion, and dielectric relaxation. Excellent performance is observed through relevant statistics of the JFNK solver, although the number of linear iterations increases for the pin-to-pin discharge when nonlinear numerical boundary conditions are imposed at the electrodes. Performance studies show scalability with O(100-1000) processors for O(10M) unknowns with ample room for optimization.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Quantitation of amiodarone and N-desethylamiodarone in single HepG2 cells by single-cell printing-liquid vortex capture-mass spectrometry

Quantitative measure of a drug and its associated metabolite(s) with single-cell resolution is often limited by sampling throughput or other compromises that limit broad use. Here, we demonstrate the use of single-cell printing-liquid vortex capture-mass spectrometry (SCP-LVC-MS) to quantitatively measure the intracellular concentrations of amiodarone (AMIO) and its metabolite, N-desethylamiodarone (NDEA), from thousands of single cells across several AMIO incubation concentrations ranging from 0 to 10 μM. Concentrations obtained by SCP-LVC-MS were validated through comparison with average assays and traditional measurement of cells in bulk. Average of SCP-LVC-MS measurements and aggregate vial collection assay the concentrations differed by < 5%. Both AMIO and NDEA had clear log-normal distributions with similar standard deviation of concentrations in the cell population. The mean of both AMIO and NDEA intracellular concentrations were positively correlated with AMIO incubation concentration, increasing from 0.026 to 0.520 and 0.0055 to 0.048 mM for AMIO and NDEA, respectively. The standard deviation of AMIO and NDEA log-normal distribution fits were relatively similar in value across incubation concentrations, 0.15–0.19 log 10 (mM), and exhibited a linear trend with respect to each other. The single cell-resolved conversion ratio of AMIO to NDEA increased with decreasing incubation concentration, 7 ± 2%, 18 ± 3%, and 20 ± 7% for 10.0, 1.0, and 0.1 μM AMIO incubation concentrations, respectively. Finally, association with simultaneously measured lipids had several ions with statistically significant difference in intensity but no clear correlations with AMIO intracellular content was observed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

O-Mode Instability for Non-Thermal Continuum Radiation Generated Near the Plasmapause

Analysis of the ordinary mode (O-mode) instability is performed to comprehend the nonthermal continuum (NTC) radiation near the plasmapause, taking into account the relativistic wave-electron resonance effect. The energy source is the anisotropy in the velocity of the minority suprathermal electron population. Numerical solutions demonstrate that the O-mode can be unstable with multiple narrow frequency bands located close to harmonics of the electron cyclotron frequency above the local electron plasma frequency. These waves have narrow beaming angle bands of nearly 90° relative to the ambient magnetic field. Our findings indicate that NTC radiation generated by this wave-electron resonance instability near the plasmapause can propagate nearer to the magnetic equator with multiple harmonics, which is in agreement with a recent statistical study using Van Allen Probes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Avalanches in strong imbibition

Slow injection of non-wetting fluids (drainage) and strongly wetting fluids (strong imbibition) into porous media are two contrasting processes in many respects: the former must be forced into the pore space, while the latter imbibe spontaneously; the former occupy pore bodies, while the latter coat crevices and corners. These two processes also produce distinctly different displacement patterns. However, both processes evolve via a series of avalanche-like invasion events punctuated by quiescent periods. Here, we show that, despite their mechanistic differences, avalanches in strong imbibition exhibit all the features of self-organized criticality previously documented for drainage, including the correlation scaling describing the space-time statistics of invasion at the pore scale.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Revisiting trends in the exchange current for hydrogen evolution

Nørskov and collaborators proposed a simple kinetic model to explain the volcano relation for the hydrogen evolution reaction on transition metal surfaces such that j 0 = k 0 f(ΔG H ) where j 0 is the exchange current density, f(ΔG H ) is a function of the hydrogen adsorption free energy ΔG H as computed from density functional theory, and k 0 is a universal rate constant. Herein, focusing on the hydrogen evolution reaction in acidic medium, we revisit the original experimental data and find that the fidelity of this kinetic model can be significantly improved by invoking metal-dependence on k 0 such that the logarithm of k 0 linearly depends on the absolute value of ΔG H . Here, we further confirm this relationship using additional experimental data points obtained from a critical review of the available literature. Our analyses show that the new model decreases the discrepancy between calculated and experimental exchange current density values by up to four orders of magnitude. Furthermore, we show the model can be further improved using machine learning and statistical inference methods that integrate additional material properties.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Variational deep learning of equilibrium transition path ensembles

Here, we present a time-dependent variational method to learn the mechanisms of equilibrium reactive processes and efficiently evaluate their rates within a transition path ensemble. This approach builds off of the variational path sampling methodology by approximating the time-dependent commitment probability within a neural network ansatz. The reaction mechanisms inferred through this approach are elucidated by a novel decomposition of the rate in terms of the components of a stochastic path action conditioned on a transition. This decomposition affords an ability to resolve the typical contribution of each reactive mode and their couplings to the rare event. The associated rate evaluation is variational and systematically improvable through the development of a cumulant expansion. We demonstrate this method in both over- and under-damped stochastic equations of motion, in low-dimensional model systems, and in the isomerization of a solvated alanine dipeptide. In all examples, we find that we can obtain quantitatively accurate estimates of the rates of the reactive events with minimal trajectory statistics and gain unique insights into transitions through the analysis of their commitment probability.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Statistical characteristics of laser interference and its effect on IR optoelectronic observation systems

We consider statistical characteristics of laser interference recorded in the focal region of an IR optoelectronic observation system under field and out-of-field illumination by laser radiation. The applicability of the analytical description of experimental histograms of the laser interference signal distribution at the photodetector array output using the probability density approximations in the form of gamma distribution and the Gaussian function is substantiated. It is shown that the specific form of these functions is determined both by the mean value (mathematical expectation) of the interference signal and by the characteristic parameter M, which depends on the statistical properties of the laser interference. (laser applications and other issues in quantum electronics)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Impact and mitigation of spectroscopic systematics on DESI DR1 clustering measurements

The large scale structure catalogs within DESI Data Release 1 (DR1) use nearly 6 million galaxies and quasars as tracers of the large-scale structure of the universe to measure the expansion history with baryon acoustic oscillations and the growth of structure with redshift-space distortions. In order to take advantage of DESI's unprecedented statistical power, we must ensure that the galaxy clustering measurements are unaffected by non-cosmological density fluctuations. One source of spurious fluctuations comes from variation in galaxy density with spectroscopic observing conditions, lowering the redshift efficiency (and thus galaxy density) in certain areas of the sky. We measure the uniformity of the redshift success rate for DESI luminous red galaxies (LRG), bright galaxies (BGS) and quasars (QSO), complementing the detailed discussion of emission line galaxy (ELG) systematics in a companion paper [1]. We find small but significant fluctuations of up to 3% in redshift success rate with the effective spectroscopic signal-to-noise, and create and describe weights that remove these fluctuations. We also describe the process to identify and remove data from certain poorly performing fibers from DESI DR1, and measure the stability of the redshift success rate with time. Finally, we find small but significant correlations of redshift success rate with position on the focal plane, survey speed, and number of exposures required, and show the impact of weights correcting these trends on the power spectrum multipoles and on cosmological parameters from BAO and RSD fits. These corrections change the best-fit parameters by <15% of their statistical errors, and thus contribute negligibly to the overall DESI error budget.

cosmological parameters from LSS↗

First Exclusive Reconstruction of the 𝐵* + , 𝐵* 0 , and 𝐵$^{*0}_{𝑠}$ Mesons and Precise Measurement of Their Masses

Using proton-proton collision data collected by the CMS experiment at √𝑠 = 13 TeV in 2016–2018, corresponding to an integrated luminosity of 140 fb −1 , the first full reconstruction of the three vector 𝐵 meson states, 𝐵* + , 𝐵* 0 , and 𝐵$^{*0}_{𝑠}$, is performed. The mass differences between the excited mesons and their corresponding ground states are measured to be 𝑚⁡(𝐵* + )−𝑚⁡(𝐵 + ) = 45.277 ± 0.039 ± 0.027 MeV, 𝑚⁡(𝐵 *0 )−𝑚⁡(𝐵 0 ) = 45.471 ± 0.056 ± 0.028 MeV, and 𝑚⁡(𝐵$^{*0}_{𝑠}$)−𝑚⁡(𝐵$^{0}_{𝑠}$) = 49.407 ± 0.132 ± 0.041 MeV, where the first uncertainties are statistical and the second are systematic. These results improve on the precision of previous measurements by an order of magnitude.

Bottom mesons↗

Grid-based diffusion Monte Carlo for fermions without the fixed-node approximation

A diffusion Monte Carlo algorithm is introduced that can determine the correct nodal structure of the wave function of a few-fermion system and its ground-state energy without an uncontrolled bias. This is achieved by confining signed random walkers to the points of a uniform infinite spatial grid, allowing them to meet and annihilate one another to establish the nodal structure without the fixed-node approximation. An imaginary-time propagator is derived rigorously from a discretized Hamiltonian, governing a non-Gaussian, sign-flipping, branching, and mutually annihilating random walk of particles. The accuracy of the resulting stochastic representations of a fermion wave function is limited only by the grid and imaginary-time resolutions and can be improved in a controlled manner. Here, the method is tested for a series of model problems including fermions in a harmonic trap as well as the He atom in its singlet or triplet ground state. For the latter case, the energies approach from above with increasing grid resolution and converge within 0.015 E h of the exact basis-set-limit value for the grid spacing of 0.08 a.u. with a statistical uncertainty of 10 –5 E h without an importance sampling or Jastrow factor.

74 ATOMIC AND MOLECULAR PHYSICS↗

Global Sensitivity Analysis of Bulk Properties of an Atomic Nucleus

We perform a global sensitivity analysis of the binding energy and the charge radius of the nucleus 16 O to identify the most influential low-energy constants in the next-to-next-to-leading order chiral Hamiltonian with two- and three-nucleon forces. For this purpose, we develop a subspace-projected coupled-cluster method using eigenvector continuation [Frame D. et al., Phys. Rev. Lett. 121, 032501 (2018)]. With this method, we compute the binding energy and charge radius of 16 O at more than 106 different values of the 16 low-energy constants in one hour on a standard laptop computer. For relatively small subspace projections, the root-mean-square error is about 1% compared to full-space coupled-cluster results. We find that 58(1)% of the variance in energy can be apportioned to a single contact term in the 3 S 1 wave, whereas the radius depends sensitively on several low-energy constants and their higher-order correlations. The results identify the most important parameters for describing nuclear saturation and help prioritize efforts for uncertainty reduction of theoretical predictions. The achieved acceleration opens up an array of computational statistics analyses of the underlying description of the strong nuclear interaction in nuclei across the Segrè chart.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

NO$\nu$A Far Detector $\nu_e$ Sideband Study Technical Note

Till now the NO$\nu$A experiment has used the Far Detector (FD) electron neutrino (anti-neutrino) events in the energy range 1 < $E_ν$ < 4 GeV in 3-favor neutrino oscillation analysis, to constraint the neutrino oscillation parameters. In this study, we have used the high-energy neutrino (sideband) sample with 4 < $E_ν$ < 12 GeV, dominated heavily by the beam electron neutrino/anti-neutrino background events to constrain the neutrino oscillation parameters. This was done to figure out if an additional power can be obtained in constraining the oscillation parameters by adding this sample to the analysis. This technote discusses the modifications done in the existing framework to add the sideband sample to the standard FD $\nu_e$ predictions and its impact on the neutrino oscillation parameters sensitivity of the NO$\nu$A experiment. Results from the statistics only study and the study including all the existing systematic shifts are reported separately.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Measurements of the 238 U/ 235 U and 239 Pu/ 235 U Fission Cross-Section Ratios Using Monoenergetic Neutron Beams

A quasi-monoenergetic neutron beam was used to measure the neutron-induced fission cross-section ratios for 238 U(n,f)/ 235 U(n,f) (Table 1) and 239 Pu(n,f)/ 235 U(n,f) (Table 2). These results are plotted in comparison with data from the fissionTPC and ENDF/B-VIII.0 in Figs.1-3. For each cross section ratio, the total, systematic, and statistical uncertainties are listed, along with the total beamtime required.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Overview of a Methodology for Calculating the A Priori Scan Minimum Detectable Concentration for Post-Processed Radiological Surveys (Final Report) (Rev.1)

Increased continuous data collection using automated data loggers and autonomous radiological survey devices or vehicles has introduced a need for corresponding guidance and statistical techniques for data that are collected without surveyor vigilance. This report presents a method for calculating the a priori scan minimum detectable concentrations (MDCs) for surveys performed without vigilance similar to methods described in the Multi-Agency Radiation Survey and Site Investigation Manual, NUREG-1507, and NUREG/CR-6364. A priori scan MDCs are calculated during survey planning to ensure that survey parameters (e.g., scanning speed, scanning altitude, detector geometry) will lead to collecting data in which potentially contaminated areas can be detected when data are processed after the survey, within acceptable statistical error probabilities.

54 ENVIRONMENTAL SCIENCES↗

Small-correlated-against-large estimator for the lensing of the cosmic microwave background

Weak gravitational lensing of the cosmic microwave background (CMB) carries imprints of the physics operating at redshifts much lower than that of recombination and serves as an important probe of cosmological structure formation, dark matter physics, and the mass of neutrinos. Reconstruction of the CMB lensing deflection field through use of quadratic estimators has proven successful with existing data but is known to be sub-optimal on small angular scales ($\ell > 3000$) for experiments with low noise levels. Future experiments will provide better observations in this regime, but these techniques will remain statistically limited by their approximations. We show that correlations between fluctuations of the large-scale temperature gradient power of the CMB sourced by $\ell < 2000$, and fluctuations of the local small-scale temperature power reveal a lensing signal which is prominent in even the real-space pixel statistics across a CMB temperature map. We present the development of the Small Correlated Against Large Estimator (SCALE), a novel estimator for the CMB lensing spectrum which offers promising complementary analysis alongside other reconstruction techniques in this regime. The SCALE method computes correlations between both the large/small-scale temperature gradient power in harmonic space, and it is able to quantitatively recover unbiased statistics of the CMB lensing field without the need for map-level reconstruction. SCALE can outperform quadratic estimator signal-to-noise by a factor of up to 1.5 in current and upcoming experiments for CMB lensing power spectra $C_{6000 < L < 8000}^{\phi\phi}$.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Progress Toward Interpretable Machine Learning–Based Disruption Predictors Across Tokamaks

Here in this paper we lay the groundwork for a robust cross-device comparison of data-driven disruption prediction algorithms on DIII-D and JET tokamaks. In order to consistently carry on a comparative analysis, we define physics-based indicators of disruption precursors based on temperature, density, and radiation profiles that are currently not used in many other machine learning predictors for DIII-D data. These profile-based indicators are shown to well-describe impurity accumulation events in both DIII-D and JET discharges that eventually disrupt. The univariate analysis of the features used as input signals in the data-driven algorithms applied on the data of both tokamaks statistically highlights the differences in the dominant disruption precursors. JET with its ITER-like wall is more prone to impurity accumulation events, while DIII-D is more subject to edge-cooling mechanisms that destabilize dangerous magnetohydrodynamic modes. Even though the analyzed data sets are characterized by such intrinsic differences, we show through a few examples that the inclusion of physics-based disruption markers in data-driven algorithms is a promising path toward the realization of a uniform framework to predict and interpret disruptive scenarios across different tokamaks. As long as the destabilizing precursors are diagnosed in a device-independent way, the knowledge that data-driven algorithms learn on one device can be re-used to explain a disruptive behavior on another device.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗