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 1,081 records · Page 60
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Convergent close-coupling calculations of electrons scattering on electronically excited molecular hydrogen

We use the adiabatic-nuclei molecular convergent close-coupling method to perform calculations of 0.01–1000 eV electrons scattering on the c 3 Π u , a 3 $Σ$$^{+}_{g}$, B 1 $Σ$$^{+}_{u}$, C 1 Π u , and EF 1 $Σ$$^{+}_{g}$ states of H 2 in the v=0 vibrational level. Elastic, superelastic, ionization, and grand-total cross sections are presented, as well as cross sections for excitation of the n=2–3 singlet and triplet states of H 2 (where n is the atomic-limit principle quantum number). Comparison with available theoretical results is made. Good agreement is found with the recent R-matrix results [J. Phys. B 53, 245203 (2020)] for most of the exchange and dipole-forbidden transitions, but not for the dipole-allowed transitions. Here, the sources of disagreement were found to be an unconverged partial-wave expansion and the utilization of the fixed-nuclei approximation (as opposed to adiabatic-nuclei) in the R-matrix calculations.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

Searching for optimal THz generation through calculations of the asymmetry of photoelectron momentum distributions by an improved strong-field-approximation method

From the solution of the time-dependent Schrödinger equation (TDSE), it has been shown that in a two-color laser pulse the relative phase that optimizes terahertz (THz) radiation also optimizes the asymmetry of the angular distribution of photoelectrons. Here, in this work, we show that a second-order strong-field approximation can accurately calculate the asymmetry of photoelectrons and thus can locate the phase delay that optimizes THz generation, doing so four orders of magnitude faster than TDSE calculations. We further trace that this is possible because THz emission originates from a free-free radiative transition between the rescattered electron and the target ion, similar to other electron-ion collisions in the laser field, e.g., high-harmonic generation and nonsequential double ionization. Our results pave the way to locate the optimal phase delay for any typical laser pulses for maximal THz generation in the laboratory.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Ab initio calculations of the rate of carrier trapping and release at dopant sites in NaI: Tl beyond the harmonic approximation

We present ab initio calculations of electron capture and release coefficients as a function of temperature for Tl dopants in the widely used scintillator NaI, which is a soft ionic crystal. The modelled capture and release events occur by transitions mediated by multiple phonon absorption and emission between states with significantly different local geometries around the trapping site. We demonstrate that such transitions are not well-described by the normal harmonic approximation to the nuclear dynamics. We avoid the harmonic approximation by numerically solving the vibrational Schrödinger equation for the motion along/of a single effective phonon coordinate. The localized trapped state is not correctly described by semilocal density functionals, so we employ a hybrid functional tuned to reproduce the band gap and lattice constant of the host material. Our calculations, which combine plane-wave and cluster electronic structure methods in a novel way, are in reasonable agreement with available experimental data for detrapping and predict an unusual temperature dependence for trapping rates.

Prange, Micah P.↗
DOE OSTI2020

Calculations of electronic excitation by protons and alpha particles in silicon

This paper presents an ab initio method of calculating the valence electron contribution to the stopping and interaction cross sections of swift ions in matter using linear response calculations of the dielectric function in crystals. We estimate the contribution to the stopping and cross sections by repeated electron capture and stripping of the projectile ion by the medium. We also predict the electronic stopping of protons and alpha particles in silicon and compare to measured values. The predictions generally agree with the data at ion speeds greater than 1 atomic unit. We find that electron capture and stripping is at most a small contribution to the stopping for both ions studied.

electron excitation, silicon, nuclear interactions↗
DOE OSTI2021

Electron correlations in the cubic paramagnetic perovskite Sr ( V , Mn ) O 3 : Results from fully self-consistent self-energy embedding calculations

In this study, we use the thermodynamically consistent and conserving self-energy embedding theory (SEET) to study the spectra of the prototypical undistorted cubic perovskites SrVO 3 and SrMnO 3 . In the strongly correlated metallic SrVO 3 we find that the usual attribution of the satellite peaks at –1.8 eV to Hund or Hubbard physics in the t 2 g orbitals is inconsistent with our calculations. In the strongly correlated insulator SrMnO 3 we recover insulating behavior due to a feedback effect between the strongly correlated orbitals and the weakly correlated environment. Our calculation shows a systematic convergence of spectral features as the space of strongly correlated orbitals is enlarged, paving the way to a systematic parameter-free study of correlated perovskites.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Effects of substitutional Ga on the phonons of δ -phase Pu based on density functional theory calculations

The delta phase of Pu is stabilized by Ga doping, but the mechanism of this stabilization remains an open question. Density functional theory calculations focused on how Ga doping affects the phonons sheds some light on the phonons' contribution to the stabilization. Here, the calculated phonon modes of Ga-doped delta phase Pu fall into two distinct types: localized, high frequency Ga-dominated phonon modes, and Pu-dominated modes at lower frequencies. Increasing the Ga concentration has an effect on the Pu-dominated phonon modes opposite to that of compression: higher-frequency modes soften, and lower-frequency modes stiffen. The latter provides an indication that the stabilization mechanism is not due to a thermodynamic contribution from the phonons. Furthermore, the stiffened phonon modes include candidate modes that describe possible pathways into low-temperature phases, suggesting that doping with Ga could impede such pathways.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Importance of nonuniform Brillouin zone sampling for ab initio Bethe-Salpeter equation calculations of exciton binding energies in crystalline solids

Excitons are prevalent in semiconductors and insulators, and their binding energies are critical for optoelectronic applications. The state-of-the-art method for first-principles calculations of excitons in extended systems is the ab initio GW-Bethe-Salpeter equation (BSE) approach, which can require a fine sampling of reciprocal space to accurately resolve solid-state exciton properties. Here, in this study, we show, for a range of semiconductors and insulators, that the commonly employed approach of uniformly sampling the Brillouin zone can lead to underconverged exciton binding energies, as impractical grid sizes are required to achieve adequate convergence. We further show that nonuniform sampling of the Brillouin zone, focused on the region of reciprocal space where the exciton wave function resides, enables efficient rapid numerical convergence of exciton binding energies at a given level of theory. We propose a well-defined convergence procedure, which can be carried out at relatively low computational cost and which in some cases leads to a correction of previous best theoretical estimates by almost a factor of 2, qualitatively changing the predicted exciton physics. These results call for the adoption of nonuniform sampling methods for ab initio GW-BSE calculations and for revisiting previously computed values for exciton binding energies of many systems.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Predicting phase transitions in PbTi⁢O 3 using zentropy through quasiharmonic phonon calculations

According to x-ray diffraction (XRD) measurements, PbTi⁢O 3 undergoes a phase transition from a tetragonal ferroelectric (FE) phase to a cubic paraelectric phase at 763 K. However, x-ray absorption fine-structure (XAFS) measurements indicate that PbTi⁢O 3 is locally tetragonal even after the phase transition. The difference in these results is because XAFS measurements can probe local features of a structure, while XRD averages over such local features. For both measurements to be consistent, PbTi⁢O 3 is macroscopically cubic but locally tetragonal after the phase transition. Despite this, most models, such as the Landau-Ginsburg-Devonshire theory and effective Hamiltonians, are still unable to explain this phenomenon. Moreover, these methods involve model parameters fitted to experimental or theoretical data and do not consider other tetragonal configurations, such as domain walls, to predict the phase transition. In our previous study, we used our zentropy approach to predict the phase transition by considering the tetragonal FE ground-state configuration and the tetragonal 90° and 180° domain wall configurations with their total energies at 0 K. Here, in this paper, the Helmholtz energies of the three configurations are obtained from density functional theory calculations through energy-volume curves and phonon calculations. The predicted phase transition temperature using the meta-GGA 𝑟 2⁢ SCAN and revised multiplicities of configurations is 716 K, showing good agreement with the experimental value of 763 K.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Calculating fluctuations and self-correlations numerically for causal charge diffusion in relativistic heavy-ion collisions

Here, we study the propagation and diffusion of electric charge fluctuations in the Bjorken hydrodynamic model with both white and Catteneo noise using purely numerical methods. We show that a lattice of noise fluctuations is required to fully calculate the two-point correlators of charge. We introduce a numerical procedure to solve the stochastic differential equations that arise from the charge conservation equation on the lattice event by event. We explicitly identify the self-correlation term in the case of Catteneo noise and provide a physical interpretation. We provide a numerical recipe to remove this contribution from the full two-point correlators. Finally, we calculate the balance functions for charged hadrons. By limiting the speed of signal propagation, we observe the expected narrowing of the balance functions after removing the self-correlations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Dirac oscillator: An alternative basis for nuclear structure calculations

The isotropic harmonic oscillator supplemented by a strong spin-orbit interaction has been the cornerstone of nuclear structure since its inception more than seven decades ago. In this paper we introduce—or rather re-introduce—the “Dirac oscillator,” a fully relativistic basis that has all the desired attributes of the ordinary harmonic oscillator while naturally incorporating a strong spin-orbit coupling. To assess the power and flexibility of the Dirac oscillator basis in the solution of nuclear structure problems within the framework of covariant density-functional theory. Here, self-consistent calculations of binding energies and ground-state densities for a selected set of doubly magic nuclei are performed using the Dirac oscillator basis and are then compared against results obtained with the often-used Runge-Kutta method. Results obtained using the Dirac oscillator basis reproduce with high accuracy those derived using the Runge-Kutta method and suggest a clear path for a generalization to systems with axial symmetry. Although the harmonic oscillator with spin-orbit corrections has been the staple of the nuclear shell model since the beginning, the Dirac oscillator is practically unknown among the nuclear physics community. In this paper we illustrate the power and flexibility of the Dirac oscillator and suggest extensions to the study of systems without spherical symmetry, as required in constrained calculations of nuclear excitations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Skyrme pseudopotentials at next-to-next-to-leading order: Construction of local densities and first symmetry-breaking calculations

There is an ongoing quest to improve on the spectroscopic quality of nuclear energy density functionals (EDFs) of the Skyrme type through extensions of its traditional form. One direction for such activities is the inclusion of terms of higher order in gradients in the EDF. We report on exploratory symmetry-breaking calculations performed for an extension of the Skyrme EDF that includes central terms with four gradients at next-to-next-to-leading order (N2LO) and for which the high-quality parametrization SN2LO1 has been constructed recently. Up to now, the investigation of such functionals with higher-order terms was limited to infinite matter and spherically symmetric configurations of singly and doubly magic nuclei. We address here nuclei and phenomena that require us to consider axial and nonaxial deformation, both for reflection-symmetric and also reflection-asymmetric shapes, as well as the breaking of time-reversal invariance. Achieving these calculations demanded a number of formal developments. These all resulted from the formulation of the N2LO EDF requiring the introduction of new local densities with additional gradients that are not present in the EDF at NLO. Their choice is not unique, but can differ in the way the gradients are coupled. While designing a numerical implementation of N2LO EDFs in Cartesian three-dimensional coordinate-space representation, we have developed a novel definition and a new unifying notation for normal and pair densities that contain gradients at arbitrary order. Besides having mnemonic advantages, the new notation allows for the easy identification of redundancies and reducibilities in a given set of local densities, and the new definition makes it straightforward to construct densities that automatically adopt the symmetries of the many-body state they are constructed from. The resulting scheme resolves several issues with some of the choices that have been made for local densities in the past, in particular when breaking time-reversal symmetry. Guided by general practical considerations, we propose an alternative form of the N2LO contribution to the Skyrme EDF that is built from a different set of densities. It has exactly the same physics content, but is much more efficient to handle in formal discussions and, compared to the original formulation, leads to a substantial reduction of computational cost and memory requirements in deformed codes. As representative examples for the performance of SN2LO1, we have chosen the ground states of even-even Kr and Nd isotopes, the fission barrier of 240 Pu as well as the superdeformed rotational band of 194 Hg. Overall, for the nuclei and phenomena studied here, the SN2LO1 parametrization does not yet present a systematic improvement over standard NLO parametrizations. This finding calls for improved fit protocols that better discriminate between NLO and N2LO terms and better exploit the unique features of the additional degrees of freedom offered by the latter.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Converged ab initio calculations of heavy nuclei

We propose a novel storage scheme for three-nucleon (3N) interaction matrix elements relevant for the normal-ordered two-body approximation used extensively in ab initio calculations of atomic nuclei. This scheme reduces the required memory by approximately two orders of magnitude, which allows the generation of 3N interaction matrix elements with the standard truncation of E 3max =28, well beyond the previous limit of 18. We demonstrate that this is sufficient to obtain the ground-state energy of 132 Sn converged to within a few MeV with respect to the E 3max truncation. In addition, we study the asymptotic convergence behavior and perform extrapolations to the un-truncated limit. Finally, we investigate the impact of truncations made when evolving free-space 3N interactions with the similarity renormalization group. We find that the contribution of blocks with angular momentum J rel > 9/2 to the ground-state energy is dominated by a basis-truncation artifact, which vanishes in the large-space limit, so these computationally expensive components can be neglected. For the two sets of nuclear interactions employed in this work, the resulting binding energy of 132 Sn agrees with the experimental value within theoretical uncertainties. This work enables converged ab initio calculations of heavy nuclei.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Ab initio coupled-cluster calculations of ground and dipole excited states in He 8

We perform coupled-cluster calculations of ground- and dipole excited-state properties of the 8 He halo nucleus with nucleon-nucleon and three-nucleon interactions from chiral effective field theory, both with and without explicit delta degrees of freedom. By increasing the precision in our coupled-cluster calculations via the inclusion of leading-order three-particle three-hole excitations in the cluster operator, we obtain a ground-state energy and a charge radius that are consistent with experiment, albeit with a slight underbinding. We also investigate the excited states induced by the electric dipole operator and present a discussion on the Thomas-Reiche-Kuhn and cluster sum rules. Lastly, we compute the electric dipole polarizability, providing a theoretical benchmark for future experimental determinations that will study this exotic nucleus.

6 ≤ A ≤ 19↗
DOE OSTI2022

Calculations of the $np→dγ$ reaction in chiral effective field theory

Here, we present a calculation of the radiative capture cross section of the $np→dγ$ reaction in the low-energy range, where the M⁢1 reaction channel dominates. Employing the LENPIC nucleon-nucleon interaction up to the fifth order (N4LO) that is regularized by the semilocal coordinate space regulators, we obtain the initial and final state wave functions, and evaluate the phase shifts of the scattering state and deuteron properties. We derive the transition operator from the chiral effective field theory up to the next-to-next-to leading order (N2LO), where we also regularize the transition operator using regulators consistent with those of the interactions. We compute the capture cross sections, and the results show a converging pattern with the chiral-order expansion of the nucleon-nucleon interaction, where the regulator dependence of the results is weak when higher-order nucleon-nucleon interactions are employed. We quantify the uncertainties of the cross-section results due to the chiral-order truncation. The chirally complete and consistent cross-section results are calculated up to N2LO and they compare well with the experiments and other theoretical predictions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Auxiliary field diffusion Monte Carlo calculations of magnetic moments of light nuclei with chiral effective field theory interactions

We calculate the magnetic moments of light nuclei $( A < 20 )$ using the auxiliary field diffusion Monte Carlo method and local two- and three-nucleon forces with electromagnetic currents from chiral effective field theory. For all nuclei under consideration, we also calculate the ground-state energies and charge radii. We generally find a good agreement with experimental values for all of these observables. For the electromagnetic currents, we explore the impact of employing two different power counting schemes, and study theoretical uncertainties stemming from the truncation of the chiral expansion order by order for select nuclei within these two approaches. In conclusion, we find that it is crucial to employ consistent power counting schemes for interactions and currents to achieve a systematic order-by-order convergence.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Quantum Monte Carlo calculations of electron scattering from 12 C in the short-time approximation

The short-time approximation is a method introduced to evaluate electroweak nuclear response for systems with A ≥ 12, extending the reach of first-principle many-body quantum Monte Carlo calculations. Using realistic two- and three-body nuclear interactions and consistent one- and two-body electromagnetic currents, we calculate longitudinal and transverse response densities and response functions of 12 C. Here, we compare the resulting cross sections with experimental data for electron-nucleus scattering, finding good agreement.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Perturbative treatment of nonlocal chiral interactions in auxiliary-field diffusion Monte Carlo calculations

Nuclear many-body systems, ranging from nuclei to neutron stars, are some of the most interesting physical phenomena in our universe, and quantum Monte Carlo (QMC) approaches are among the most accurate many-body methods currently available to study them. In recent decades, interactions derived from chiral effective field theory (EFT) have been widely adopted in the study of nuclear many-body systems. One drawback of the QMC approach is the requirement that the nuclear interactions need to be local, whereas chiral EFT interactions usually contain nonlocalities. In this work, we leverage the capability of computing second-order perturbative corrections to the ground-state energy in order to develop a self-consistent approach to including nonlocal operators in QMC calculations. In conclusion, we investigate both the deuteron and the neutron-matter equation of state in order to show the robustness of our technique and pave the way for future QMC calculations at higher orders in the EFT, where nonlocal operators cannot be avoided.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Constraining calculations of Pu 239 ( n , γ ) cross sections with multistep cascade γ -ray spectra measured at the DANCE facility

Measurements of multistep cascade gamma-ray spectra from several low-lying 1 + resonances in the Pu 239 ( n , γ ) reaction were made using the DANCE detector. The measured spectra are compared to spectra simulated using several models for the photon strength function and the nuclear level density, representing the first analysis of this nature on a fissile nucleus. In general, models of the photon strength function that included a strong M 1 “scissors mode” contribution coupled with a Gilbert-Cameron nuclear level density provided a reasonable description of the spectra. Simulations using photon strength functions and nuclear level densities specifically for Pu 240 , extracted by the Oslo method, resulted in a very good representation of the measured spectra. Neutron-capture cross section calculations that included the scissors mode provided a good estimate of the previously measured values, without additional normalization. However, the scissors-mode energy that produced the best description of the multistep cascade spectra was different from the value that best reproduced the capture cross section. In addition, the 〈 Γ γ 〉 inferred from the multistep cascade calculations was somewhat less than the tabulated value, which was determined from resonance analysis. Published by the American Physical Society 2025

Ullmann, J. L. (ORCID:0000000207784405)↗