Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “calculator”

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

Quantum many-body calculations using body-centered cubic lattices

It is often computationally advantageous to model space as a discrete set of points forming a lattice grid. This technique is particularly useful for computationally difficult problems such as quantum many-body systems. For reasons of simplicity and familiarity, nearly all quantum many-body calculations have been performed on simple cubic lattices. Since the removal of lattice artifacts is often an important concern, it would be useful to perform calculations using more than one lattice geometry. In this paper we show how to perform quantum many-body calculations using auxiliary-field Monte Carlo simulations on a three-dimensional body-centered cubic (BCC) lattice. As a benchmark test we compute the ground state energy of 33 spin-up and 33 spin-down neutrons in the unitary limit, which is an idealized limit where the interaction range is zero and scattering length is infinite. As a fraction of the free Fermi gas energy E FG , we find that the ground state energy is E 0 /E FG =0.369(2),0.371(2), using two different definitions of the finite-system energy ratio. This is in excellent agreement with recent results obtained on a cubic lattice [He et al., Phys. Rev. A 101, 063615 (2020)]. We find that the computational effort and performance on a BCC lattice is approximately the same as that for a cubic lattice with the same number of lattice points. We discuss how the lattice simulations with different geometries can be used to constrain the size of lattice artifacts in simulations of continuum quantum many-body systems.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Ab initio calculation of the β decay from 11 Be to a 10 Be + $\textit{p}$ resonance

In this work, the exotic β -delayed proton emission is calculated in 11 Be from first principles using chiral two- and three-nucleon forces. To investigate the unexpectedly large branching ratio measured in Ayyad et al. [Phys. Rev. Lett. 123, 082501 (2019)], we calculate the proposed (1/2 + , 1/2) proton resonance in 11 B using the no-core shell model with continuum. This calculation helps to address whether this enhancement is caused by unknown dark decay modes or an unobserved proton resonance. We report a branching ratio of $b_p$ = (1.3 ± 0.5) × 10 –6 , suggesting that its unexpectedly large value is caused by an unobserved proton resonance in 11 B.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Detailed analysis of excited-state systematics in a lattice QCD calculation of 𝑔 𝐴

Excited state contamination remains one of the most challenging sources of systematic uncertainty to control in lattice QCD calculations of nucleon matrix elements and form factors: early time separations are contaminated by excited states and late times suffer from an exponentially bad signal-to-noise problem. High-statistics calculations at large time separations ≳ 1 fm are commonly used to combat these issues. In this work, focusing on g A , we explore the alternative strategy of utilizing a large number of relatively low-statistics calculations at short to medium time separations (0.2–1 fm), combined with a multistate analysis. On an ensemble with a pion mass of approximately 310 MeV and a lattice spacing of approximately 0.09 fm, we find this provides a more robust and economical method of quantifying and controlling the excited state systematic uncertainty. A quantitative separation of various types of excited states enables the identification of the transition matrix elements as the dominant contamination. The excited state contamination of the Feynman-Hellmann correlation function is found to reduce to the 1% level at approximately 1 fm while, for the more standard three-point functions, this does not occur until after 2 fm. Critical to our findings is the use of a global minimization, rather than fixing the spectrum from the two-point functions and using them as input to the three-point analysis. We find that the ground state parameters determined in such a global analysis are stable against variations in the excited state model, the number of excited states, and the truncation of early-time or late-time numerical data.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Relativistic effects in Green's function Monte Carlo calculations of neutrino-nucleus scattering

Here, microscopic calculations of neutrino-nucleus scattering cross sections are critical for the success of the neutrino-oscillation program. In addition to retaining nuclear correlations in the initial and final state of the reaction, they are based on consistent nuclear interactions and transition current operators, thereby enabling robust uncertainty quantification. In this work, we address a significant limitation of these microscopic methods, which arises from their nonrelativistic nature. By performing the calculations in a reference frame that minimizes nucleon momenta and utilizing the so-called two-fragment model, we extend the applicability of Green's function Monte Carlo calculations of neutrino-nucleus scattering to higher momenta than currently possible. To validate this approach, we compare our theoretical predictions against inclusive data measured by the MiniBooNE, T2K, and MINER ν A experiments.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Fully ab-initio all-electron calculation of dark matter-electron scattering in crystals with evaluation of systematic uncertainties

We calculate target-material responses for dark matter–electron scattering at the all-electron level using atom-centered Gaussian basis sets. The all-electron effects enhance the material response at high momentum transfers from dark matter to electrons, q ≳ O ( 10 α m e ) , compared to calculations using conventional plane wave methods, including those used in ; this enhances the expected event rates at energy transfers E ≳ 10 eV , especially when scattering through heavy mediators. We carefully test a range of systematic uncertainties in the theory calculation, including those arising from the choice of basis set, exchange-correlation functional, number of unit cells in the Bloch sum, k -mesh, and neglect of scatters with very high momentum transfers. We provide state-of-the-art crystal form factors, focusing on silicon and germanium. Our code and results are made publicly available as a new tool, called (“”). Published by the American Physical Society 2024

Dreyer, Cyrus E.↗
DOE OSTI2024

Calculating the two-photon exchange contribution to K L → μ + μ − decay

We present a theoretical framework within which both the real and imaginary parts of the complex, two-photon exchange amplitude contributing to K L → μ + μ − decay can be calculated using lattice quantum chromodynamics. The real part of this two-photon amplitude is of approximately the same size as that coming from a second-order weak strangeness-changing neutral-current process. Thus a test of the standard model prediction for this second-order weak process depends on an accurate result of this two-photon amplitude. A limiting factor of our proposed method comes from low-energy three-particle π π γ states. The contribution from these states will be significantly distorted by the finite volume of our calculation—a distortion for which there is no available correction. However, a simple estimate of the contribution of these three-particle states suggests their contribution to be at most a few percent allowing their neglect in a lattice calculation with a 10% target accuracy. Published by the American Physical Society 2024

Astronomy & Astrophysics↗
DOE OSTI2020

Bypassing the Energy Functional in Density Functional Theory: Direct Calculation of Electronic Energies from Conditional Probability Densities

Density functional calculations can fail for want of an accurate exchange-correlation approximation. The energy can instead be extracted from a sequence of density functional calculations of conditional probabilities (CP DFT). Simple CP approximations yield usefully accurate results for two-electron ions, the hydrogen dimer, and the uniform gas at all temperatures. CP DFT has no self-interaction error for one electron, and correctly dissociates H2, both major challenges. For warm dense matter, classical CP DFT calculations can overcome the convergence problems of Kohn-Sham DFT.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Lattice QCD Calculation of Electroweak Box Contributions to Superallowed Nuclear and Neutron Beta Decays

We present the first lattice QCD calculation of the universal axial γ W -box contribution □ γ W V A to both superallowed nuclear and neutron beta decays. This contribution emerges as a significant component within the theoretical uncertainties surrounding the extraction of | V u d | from superallowed decays. Our calculation is conducted using two domain wall fermion ensembles at the physical pion mass. To construct the nucleon four-point correlation functions, we employ the random sparsening field technique. Furthermore, we incorporate long-distance contributions to the hadronic function using the infinite-volume reconstruction method. Upon performing the continuum extrapolation, we arrive at □ γ W V A = 3.65 ( 7 ) lat ( 1 ) PT × 10 − 3 . Consequently, this yields a slightly higher value of | V u d | = 0.973 86 ( 11 ) exp ( 9 ) RC ( 27 ) NS , reducing the previous 2.1 σ tension with the CKM unitarity to 1.8 σ . Additionally, we calculate the vector γ W -box contribution to the axial charge g A , denoted as □ γ W V V , and explore its potential implications. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Optimally tuned starting point for single-shot GW calculations of solids

The dependence of ab initio many-body perturbation theory within the $\textit{GW}$ approximation on the eigensystem used in calculating quasiparticle corrections limits this method's predictive power. Here, we investigate the accuracy of the recently developed Wannier-localized optimally tuned screened range-separated hybrid (WOT-SRSH) functional as a generalized Kohn-Sham starting point for single-shot $\textit{GW}$ ($G_0W_0$) calculations for a range of semiconductors and insulators. Comparison to calculations based on well-established functionals, namely, PBE, PBE0, and HSE, as well as to self-consistent GW schemes and to experiment, shows that band gaps computed via $G_0W_0$@WOT-SRSH have a level of precision and accuracy that is comparable to that of more advanced methods such as quasiparticle self-consistent $\textit{GW}$ and eigenvalue self-consistent $\textit{GW}$. In this work, we also find that $G_0W_0$@WOT-SRSH improves the description of states deeper in the valence band manifold. Finally, we show that $G_0W_0$@WOT-SRSH significantly reduces the sensitivity of computed band gaps to ambiguities in the underlying WOT-SRSH tuning procedure.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Calculated electron paramagnetic resonance g-tensor and hyperfine parameters for zinc vacancy and N related defects in ZnO

Various defects in ZnO, focused on substitutional N O and N 2 in various sites, O-site, interstitial and Zn-site are studied using first-principles calculations with the goal of understanding the electron paramagnetic resonance (EPR) center reported for N 2 in ZnO and substitutional N on the O-site. The g tensors are calculated using the gauge including projector augmented wave (GIPAW) method and compared with experiments. The g-tensor of the free $N^{+}_{2}$ and $N^{–}_{2}$ radicals and their various contributions within the GIPAW theory are analyzed first to provide a baseline reference for the accuracy of the method and for understanding the N 2 behavior in ZnO. Previous controversies on the site location of N 2 in ZnO for this EPR center and on the shallow or deep nature and donor or acceptor nature of this center are resolved. Here, we find that the N 2 on the Zn site is mostly zinc-vacancy like in its spin density and g-tensor, while for the O-site, a model with the N 2 axis lying in-the basal plane and the singly occupied π g -orbital along the c axis provides good agreement with experiment. For the interstitial location, if the N 2 is not strongly interacting with the surroundings, no levels in the gap are found and hence also no possible EPR center. The calculated g-tensors for N O and V Zn are also found to be in good agreement with experiment. The effects of different functionals affecting the localization of the spin density are shown to affect the g-tensor values.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Large-scale calculations of 𝛽-decay rates and implications for 𝑟-process nucleosynthesis

Nuclear 𝛽 decay is a key element of the astrophysical rapid neutron capture process (𝑟 process). In this work, we present state-of-the-art global 𝛽-decay calculations based on the quantified relativistic nuclear energy density functional theory and the deformed proton-neutron quasiparticle random-phase approximation. Our analysis considers contributions from allowed and first-forbidden transitions. We used two point-coupling functionals with carefully calibrated time-odd terms and isoscalar pairing strength. The new calculations display consistent results for both employed functionals, especially near the neutron drip line, suggesting slower 𝛽 decays past the 𝑁=126 neutron shell closure than in commonly used 𝛽-decay models. The new rates, along with the existing rates based on the recent nonrelativistic global calculations, are found to slow down the synthesis of heavy elements in the 𝑟 process and significantly reduce the contribution of neutron-induced fission.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Glauber-Theory Calculations of High-Energy Nuclear Scattering Observables Using Variational Monte Carlo Wave Functions

Experiments using intermediate- to high-energy radioactive nuclear beams present numerous findings. Extracting important properties of physical observables relies on a firm theoretical analysis. Though Glauber theory is believed to work well, no convincing calculation has so far been done. Here, we perform ab initio Glauber theory calculations of both elastic differential cross sections and total reaction cross sections for p+ 12 C, 12 C+ 12 C, and 6 He+ 12 C systems. The wave functions of both 6 He and 12 C are generated by variational Monte Carlo calculations with spatial and spin-isospin correlations induced by realistic two- and three-nucleon potentials. Glauber’s phase-shift function is computed by Monte Carlo integration up to all orders of nucleon-nucleon multiple scatterings. We show an excellent performance of the Glauber description to the selected data on the above systems. We also find that the cumulant expansion of the phase-shift function converges rapidly up to the second order for the above systems. This finding will open up interesting applications for the analysis of high-energy nuclear experiments.

Horiuchi, W. [Osaka Metropolitan University (Japan↗
DOE OSTI2021

PyCrystalField : software for calculation, analysis and fitting of crystal electric field Hamiltonians

PyCrystalField is a Python software package for calculating single-ion crystal electric field (CEF) Hamiltonians. This software can calculate a CEF Hamiltonian ab initio from a point charge model for any transition or rare earth ion in either the J basis or the LS basis, perform symmetry analysis to identify nonzero CEF parameters, calculate the energy spectrum and observables such as neutron spectrum and magnetization, and fit CEF Hamiltonians to any experimental data. Here, the theory, implementation and examples of its use are discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2020

A ray-tracing algorithm for ab-initio calculation of thermal load in undulator-based synchrotron beamlines

The OASYS suite and its powerful integration features are used to implement a ray-tracing algorithm to accurately calculate the thermal load in any component of an undulator-based synchrotron beamline. This is achieved by sampling and converting the SRW source of a given energy into a Shadow source and using the latter code to ray trace the full beamline. The accuracy of the algorithm is proved by reconstructing the full undulator radiation distribution through an aperture and comparing the result with direct calculaton of the total power using SRW. The algorithm is particularly suited to analyze cases with complex beamline layouts and optical elements, such as crystals, multilayers, and compound refractive lenses. Examples of its use to calculate the power load on elements of two of the feature beamlines at the Advanced Photon Source Upgrade Project and a comparison of the results with analytical calculations are presented.

43 PARTICLE ACCELERATORS↗
DOE OSTI2020

Evaluating the Incident Energy of Arcs in Photovoltaic DC Systems: Comparison Between Calculated and Experimental Data

Solar Photovoltaic (PV) systems have permeated the energy generation world at a very high rate, some of the safety codes and standards are still lagging in accurately assessing the hazards and risks associated with PV array arcing energies. Safety professionals and maintenance workers using NFPA 70E have utilized the Doan, Stokes & Oppenlander or Enrique models, meant to determine arc energies in DC power systems using the maximum power method. These methods may lead to an overestimate of energy available in PV systems during a fault. Since PV modules/arrays are non-linear, current limited DC devices, some of these calculation methods may not accurately predict fault energy. This paper will validate current arc energy models for PV systems by comparing experimental and calculated data. Additionally, this data will help modify the current NFPA 70E models related to smaller solar arrays. Understanding where the real safety threshold for DC arc flash in PV systems exists will help maintenance and safety professionals better prepare for a variety of work related activities. This paper will analyze real arc data taken for PV systems <1000VDC and <60amps and compare this to the calculated incident energy models, to include 70E. Thus, using these comparisons, it may be possible to reduce the safety hazard severity and thus relax the PPE requirements for installation and maintenance crews.

14 SOLAR ENERGY↗
DOE OSTI2023

Shell-model calculations for two-neutron transfer near the $N=20$ island of inversion

We describe connections between shell-model calculation results and predictions for two-nucleon transfer reactions. Measurements of the 30 Mg(t,p) 32 Mg reaction were used to identify a low-lying shape-coexisting 0 + state in 32 Mg. Interpretations of those early, limited data were based on simple empirical models. The cross sections for two-nucleon transfer are, however, extremely sensitive to the details of the nuclear-structure. Motivated by the possibility of new high-resolution two-nucleon transfer measurements around the N = 20 Island of Inversion, we have performed shell-model calculations wave functions for states in the nuclei 30,32 Mg using the SDPF-M and SDPF-MU interactions, and from those calculations obtain the two-nucleon transfer amplitudes that we use to predict cross sections for the 28,30 Mg(t,p) 30,32 Mg reactions. The combined results show how the nuclear structure properties influence experimental observables for the (t,p) reaction. We compare our results with the existing data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Graph Contractions for Calculating Correlation Functions in Lattice QCD

Computing correlation functions for many-particle systems in Lattice QCD is vital to extract nuclear physics observables like the energy spectrum of hadrons such as protons. However, this type of calculation has long been considered to be very challenging and computing-resource intensive because of the complex nature of a hadron composed of quarks with many degrees of freedom. In particular, a correlation function can be calculated through a sum of all possible pairs of quark contractions, each of which is a batched tensor contraction, dictated by Wick's theorem. Because the number of terms of this sum can be very large for any hadronic system of interest, fast evaluation of the sum faces several challenges: an extremely large number of contractions, a huge memory footprint at runtime, and the speed of tensor contractions. In this paper, we present a Lattice QCD analysis software suite, Redstar, which addresses these challenges by utilizing novel algorithmic and software engineering methods targeting modern computing platforms such as many-core CPUs and GPUs. In particular, Redstar represents every term in the sum of a correlation function by a graph, applies efficient graph algorithms to reduce the number of contractions to lower the cost of computations, and minimizes the total memory footprint. Moreover, Redstar carries out the contractions on either CPUs or GPUs utilizing an internal and highly efficient Hadron contraction library. Specifically, we illustrate some important algorithmic optimizations of Redstar, show various key design features of Hadron library, and present the speedup values due to the optimizations along with performance figures for calculating six correlations functions on four computing platforms.

Chen, Jie↗
DOE OSTI2020

Code for Calculating Regional Seismic Travel Time

The RSTT software computes predictions of the travel time of seismic energy traveling from a source to a receiver through 2.5D models of the seismic velocity distribution within the Earth. The two primary applications for the RSTT library are tomographic inversion studies and seismic event location calculations. In tomographic inversions studies, a seismologist begins with number of source-receiver travel time observations and an initial starting model of the velocity distribution within the Earth. A forward travel time calculator, such as the RSTT library, is used to compute predictions of each observed travel time and all of the residuals (observed minus predicted travel time) are calculated. The Earth model is then modified in some systematic way with the goal of minimizing the residuals. The Earth model obtained in this way is assumed to be a better model than the starting model if it has lower residuals. The other major application for the RSTT library is seismic event location. Given an Earth model, an initial estimate of the location of a seismic event, and some number of observations of seismic travel time thought to have originated from that event, location codes systematically modify the estimate of the location of the event with the goal of minimizing the difference between the observed and predicted travel times. The second application, seismic event location, is routinely implemented by the military as part of its effort to monitor the Earth for nuclear tests conducted by foreign countries.

BALLARD, SANFORD↗